abaqusPython后处理的非典型应用

JingheSu

本来是为上次那个后处理帖子准备的，不过当时没调试出来（主要是不知道怎么写ODB），这周花了些时间弄了一下。

实现的功能：
从一个odb结果文件中copy建立另一个新的ODB

本身可能没有多大意思，不过可以给入门的人作为一个例子用，涉及odb中material的读取，各种set的读取与建立。
另外：这个只针对我所附的INP的odb有效（CAX4R单元，轴对称问题），其他的odb需要对程序作相应修改。

1. 
   
2. from odbAccess import *
   
3. from abaqusConstants import *
   
4. ##
   
5. originODB   =   'testODB.odb'
   
6. odb0    =   openOdb(path=originODB)
   
7. ########
   
8. from odbMaterial import *
   
9. from odbSection import *
   
10. from odbMaterial import *
    
11. ##create odb file name odb1.odb
    
12. analysisTitle0  =   odb0.analysisTitle
    
13. description0    =   odb0.description
    
14. myodb = Odb(name='myData',analysisTitle=analysisTitle0,description=description0,path='odb1.odb')
    
15. root0   =   odb0.rootAssembly
    
16. parts0  =   odb0.parts
    
17. myroot   =   myodb.rootAssembly
    
18. ##copy material from previous odb
    
19. allMaterials = odb0.materials
    
20. for material0 in allMaterials.values():
    
21.     Mname=material0.name
    
22.     Mmaterial=myodb.Material(name=Mname)
    
23.     if hasattr(material0,'elastic'):
    
24.         elastic = material0.elastic
    
25.         MEtype=elastic.type
    
26.         MEtemp=elastic.temperatureDependency
    
27.         MEdep=elastic.dependencies
    
28.         MEdata=elastic.table
    
29.         MEnocomp=elastic.noCompression
    
30.         MEnoten=elastic.noTension
    
31.         MEmodu=elastic.moduli
    
32.         Mmaterial.Elastic(type=MEtype,temperatureDependency=MEtemp, dependencies=MEdep,noCompression=MEnocomp, noTension=MEnoten,moduli=MEmodu, table=MEdata)        
    
33.     if hasattr(material0,'plastic'):
    
34.         plastic = material0.plastic
    
35.         PEtable=plastic.table
    
36.         Mmaterial.Plastic(table=PEtable)
    
37. ##copy sectioncategory from previous odb
    
38. allsectionca = odb0.sectionCategories
    
39. for sectionca0 in allsectionca.values():
    
40.     Sname=sectionca0.name
    
41.     Sdescri=sectionca0.description
    
42.     Spoint=sectionca0.sectionPoints
    
43.     mysectionca=myodb.SectionCategory(name=Sname,description=Sdescri)
    
44. ##copy section from previous odb
    
45. from odbSection import *
    
46. allSections = odb0.sections
    
47. for Section0 in allSections.values():
    
48.     if type(Section0) == HomogeneousSolidSectionType:
    
49.         mysectionN=Section0.name
    
50.         mysectionM=Section0.material
    
51.         mysectiont=Section0.thickness
    
52.         mysection=myodb.HomogeneousSolidSection(name=mysectionN,material=mysectionM,thickness=mysectiont)
    
53. ##create part
    
54. for part0 in parts0.values():
    
55.     mypart=myodb.Part(name=part0.name,embeddedSpace=THREE_D,type=part0.type)
    
56.     instName0=part0.name+'-1'
    
57.     inst0=root0.instances[instName0]
    
58.     mynodes=inst0.nodes
    
59.     myelements=inst0.elements
    
60.     nodelist=[]
    
61.     for i in range(len(mynodes)):
    
62.         idnodes=(mynodes[i].label,mynodes[i].coordinates[0],mynodes[i].coordinates[1],mynodes[i].coordinates[2])
    
63.         nodelist.append(idnodes)
    
64.     elemlist=[]
    
65.     for i in range(len(myelements)):
    
66.         idelem=(myelements[i].label,myelements[i].connectivity[0],myelements[i].connectivity[1],myelements[i].connectivity[2],myelements[i].connectivity[3])
    
67.         elemlist.append(idelem)
    
68.     mypart.addNodes(nodeData=nodelist,nodeSetName=part0.name)
    
69.     mypart.addElements(elementData=elemlist, type='CAX4R', elementSetName=part0.name, sectionCategory=mysectionca)
    
70. ##create instance & nodeset & elementset
    
71. instanceN=[]
    
72. for instance0 in root0.instances.values():
    
73.     instanceN.append(instance0.name)
    
74.     mypart=myodb.parts[instance0.name[:-2]]
    
75.     myinstance=myroot.Instance(name=instance0.name,object=mypart)
    
76. for nodeset0 in root0.nodeSets.values():
    
77.     nodeset0Label=[]
    
78.     for i in range(len(instanceN)):
    
79.         nodeset0Label.append([instanceN[i],[]])
    
80.         for node in nodeset0.nodes[0]:        
    
81.             if node.instanceName==instanceN[i]:
    
82.                 nodeset0Label[i][1].append(node.label)
    
83.     nodeset=myroot.NodeSetFromNodeLabels(name=nodeset0.name,nodeLabels=nodeset0Label)
    
84. for eleset0 in root0.elementSets.values():
    
85.     eleset0Label=[]
    
86.     for i in range(len(instanceN)):
    
87.         eleset0Label.append([instanceN[i],[]])
    
88.         for element in eleset0.elements[0]:        
    
89.             if element.instanceName==instanceN[i]:
    
90.                 eleset0Label[i][1].append(element.label)
    
91.     eleset=myroot.ElementSetFromElementLabels(name=eleset0.name,elementLabels=eleset0Label)
    
92. ## assign section to instance
    
93. eleset0=root0.elementSets[' ALL ELEMENTS']
    
94. eleset0Label=[]
    
95. for i in range(len(instanceN)):
    
96.     eleset0Label=[]
    
97.     for element in eleset0.elements[0]:
    
98.         if element.instanceName==instanceN[i]:
    
99.             eleset0Label.append(element.label)
    
100.     eleset=myroot.instances[instanceN[i]].ElementSetFromElementLabels(name=instanceN[i],elementLabels=eleset0Label)
     
101.     myroot.instances[instanceN[i]].assignSection(region=eleset,section=mysection)
     
102. #for surface0 in root0.surface.values():
     
103. #    print surface0.name
     
104. ##
     
105. for step0 in odb0.steps.values():
     
106.     myname=step0.name
     
107.     mydescription=step0.description
     
108.     mydomain=step0.domain
     
109.     mytimePeriod=step0.timePeriod
     
110.     mypreviousStepName=step0.previousStepName
     
111.     mytotalTime=step0.totalTime
     
112.     mystep=myodb.Step(name=myname,description=mydescription,domain=mydomain,timePeriod=mytimePeriod,previousStepName=mypreviousStepName,totalTime=mytotalTime)
     
113.     firstFrame = step0.frames[0]
     
114.     myframeincr=firstFrame.incrementNumber
     
115.     myframevalue=firstFrame.frameValue
     
116.     myframedescription=firstFrame.description
     
117.     frame0 = mystep.Frame(incrementNumber=myframeincr,frameValue=myframevalue,description=myframedescription)
     
118.     for foutput in firstFrame.fieldOutputs.values():
     
119.         frame0.FieldOutput(foutput)        
     
120.     lastFrame = step0.frames[-1]
     
121.     myframeincr=lastFrame.incrementNumber
     
122.     myframevalue=lastFrame.frameValue
     
123.     myframedescription=lastFrame.description
     
124.     frame1 = mystep.Frame(incrementNumber=myframeincr,frameValue=myframevalue,description=myframedescription)
     
125.     for foutput in lastFrame.fieldOutputs.values():
     
126.         frame1.FieldOutput(foutput)
     
127. ##
     
128. odb0.close()
     
129. myodb.save()
     

复制代码

JingheSu

CreatODB

JingheSu

Surface 2X2 Data plot using abaqus

这两楼先占着，后面续上。
给自己一点压力，不然太懒了。

upbeat928

谢谢楼主的贡献！我相信python会红火起来~，特别有abaqus这个大牛

nature_pma

这个用处不是特别大  因为一般后处理和cae都是在同一个人的  不过楼主另外一个关于python后处理的还是很好的

wenxianzheng

多谢楼主

yuqi9461

楼主这个太厉害了。。我最近在学习Python写脚本，想创建新的odb，看了您上面这个例子，我想问个问题，如果我的模型里面，同时有4节点和8节点的单元，那生成新的odb时候，在create part，读取单元时候应该怎么写呀？（对应您的代码64~69行。）我现在采用的方法是判断节点数来分别读取，但在读取8节点的时候总是报错，该如何处理呀。如果能解决，不胜感激吖~~

yyl2013121245

谢谢楼主

liminggogo

楼主，我想知道怎样用python导出abaqus的。rpt文件，因为我想要那里面的应力，位移之类的数据。