Messages

[AEM 535 HW-7 Taylor Impact Test Simulation]

AEM 535 HW-7 Taylor Impact Test Simulation

http://forums.mbarkey.com/index.php?action=downloads;sa=view;down=17

from Matt F.

During your video I would mention that when editing the Python script which was generated by Abaqus CAE make sure to change ALL lines of code which might save the file to a particular location. While following along with your video I saved my file periodically, but when I changed the Python script to save the file to a new location I only edited the last save command at the very bottom of the file, and left the original file name in some of the previous save commands. This caused me to overwrite some data but it was easily fixed. In hindsight this should have been obvious to me, but perhaps it's worth mentioning in your video.

From Greg L.

Suggestions for an Updated Assignment:
Knowing that the general arrangement of the optimized mesh is accurate with respect to the convergence analysis solution, one possible update to the assignment could be to create one more mesh in Hypermesh but change the radius size of the hole. The same meshing style (number of elements biased towards the hole, transition zone, etc.) should be relatively quick to make. The boundary conditions would be the same and the new loading condition could be provided to ensure the normalized solution. The results of this new analysis could then be compared to the theoretical stress concentration chart provided in class. This effort could be used to show how known meshing solutions can be scaled up for similar applications.

HW2: Use a rectangular plate with 2:1 aspect ratio to get Kt more near 2.5 instead of 2.6.  (Ben W.)

HW1:  pipe in tension: I need to be sure to use U3 instead of Umagnitude for the displacement of the pipe.

[AEM 535 HW-6 part 2 Read and write temperatures to the ODB]

AEM 535 HW-6 part 2 Read and write temperatures to the ODB

http://forums.mbarkey.com/index.php?action=downloads;sa=view;down=16

[AEM 535 HW-6 part 1 Thermal Analysis]

AEM 535 HW-6 part 1 Thermal Analysis

http://forums.mbarkey.com/index.php?action=downloads;sa=view;down=15

[AEM 535 HW-5 Natural Frequencies of a Beam]

AEM 535 HW-5 Natural Frequencies of a Beam

http://forums.mbarkey.com/index.php?action=downloads;sa=view;down=14

These are my notes to update the FEA assignments.

update report format to include element type, software and version.

refer to posts related to each assignment for updates related to that assignment

[HW-4 Plasticity in Plate with a Hole]

HW-4 Plasticity in Plate with a Hole

http://forums.mbarkey.com/index.php?action=downloads;sa=view;down=13

[HW-3 Plate with a hole, optimized mesh]

HW-3 Plate with a hole, optimized mesh

http://forums.mbarkey.com/index.php?action=downloads;sa=view;down=12

[HW-2 Plate with a hole]

HW-2 Plate with a hole

http://forums.mbarkey.com/index.php?action=downloads;sa=view;down=11

[HW-1 Tension on Pipe input file]

HW-1 Tension on Pipe input file

http://forums.mbarkey.com/index.php?action=downloads;sa=view;down=10

[simple XFEM example]

simple XFEM example

http://forums.mbarkey.com/index.php?action=downloads;sa=view;down=9

data files are here:
http://forums.mbarkey.com/index.php?topic=124.0

Code Select Expand



#write selected frame results for nodal temperatures to an ODB file
#or write nodal temperature from a file to an ODB file
#M. E. Barkey
# 14 December 2016

#with reference from python code from Mustaine
#ww3.cad.de/foren/ubb/Forum254/HTML/000798.shmtl
#and
#www.eng-tips.com/viewthread.cfm?qid=306565
#and
#http://www.eng-tips.com/viewthread.cfm?qid=248461




#some of these may not be needed, but I haven't taken the time to sort it out.
import sys
import csv
from odbAccess import *
from abaqusConstants import *
from types import IntType
from numpy import array
from numpy import genfromtxt
import numpy as np
import odbAccess


odbName = 'spt1b'
odb = openOdb(odbName + '.odb', readOnly=False)

#The following is the target ODB
odbWrite=openOdb(path="spt1c.odb", readOnly=False)

print'------------------------------------------------'
print'Define step and frame to be read from existing ODB'
steptoread=odb.steps['Step-2']
frametoread=steptoread.frames[4]
odbSelectResults=frametoread.fieldOutputs['NT11']
t = odbSelectResults

#these can be uncommented for some checks of the data
#print '-------------------------------------------'
#print t.__members__
#print '-------------------------------------------'
#print t.values[1]

print odb.rootAssembly.instances.keys()


#note that this must be the part in your file
instance1 = odbWrite.rootAssembly.instances['PART-1-1']


#the following lines write the entire frame of the first ODB into the target ODB
#the step name, desctription, etc, can all be changed here
newDataSet = odbSelectResults
newResultsStep=odbWrite.Step(name='TestStep3', description = 'user defined', domain=TIME, timePeriod = 0)
newResultsFrame=newResultsStep.Frame(incrementNumber=0, frameValue=0.0)
newResultsField=newResultsFrame.FieldOutput(name='NT11', description='user def', type=SCALAR)
newResultsField.addData(field=newDataSet)
odbWrite.save()



#the following lines read a single column text file of nodes
#then reads a single column text file of corresponding temepatures
#note that the scalar field must be a "sequence of sequences"
#so it must look like this:    (  (temp1,) ,   (temp2,),   (temp3,),   (temp4,) ... )
#and NOT this: (temp1, temp2, temp3, temp4,...)

with open('nodelist.txt') as f:
    nodes = f.read().splitlines()

#makes sure this is a tuple of integer nodes   
n = map(int, nodes)   



#this line reads in the temperature values.  They must be sequentially corressponding to the nodes.
#it puts them in the proper format (sequence of sequence)   
with open('values.txt') as f:
    mylist = [tuple(map(float, i.split(','))) for i in f]   
   

v = mylist

#print v


#the following is the format of the data that will work for 6 nodes
#meblabel = (41312, 41313, 41314, 41315, 41316, 41317)
#mebdata = ( (-100.0,), (-100.0,),(-100.0,),(-100.0,),(-100.0,),(-100.0,))


#un-comment the following line if you want the read-in data to go to a new step
#newResultsStep=odbWrite.Step(name='TestStep4', description = 'user defined', domain=TIME, timePeriod = 0)

newResultsFrame=newResultsStep.Frame(incrementNumber=1, frameValue=1.0)
newResultsField=newResultsFrame.FieldOutput(name='NT11', description='user def', type=SCALAR)

#this is to double check that the correct key has been created--it should be NT11
print '-------------------------------------------'
print newResultsFrame.fieldOutputs.keys('')
print '-------------------------------------------'


newResultsField.addData(position=NODAL, instance=instance1, labels=n, data=v)


odbWrite.save()


print
print '-------------------------------------------'


odb.save()
odb.close()