基于Python的SAP2000计算结果转换为ABAQUS

旧言虐心

本文基于博文”基于Python的SAP2000模型数据和结果数据读取“继续开发了通过Python语言将SAP2000计算结果导入ABAQUS,如下图所示：

                           

代码如下所示：

1. <span style="word-wrap: normal; line-height: 22.5px; font-size: 12px;">#-*-coding:mbcs-*-</span>
   
2. #***********************************************************************************************************
   
3. # <wbr> <wbr> -----------------本程序写入ABAQUS_ODB模型数据和结果数据--------------------------
   
4. # <wbr> <wbr> 1.不同单元类型定义成不同的PART、INSTANCE
   
5. # <wbr> <wbr> 2.不同的PART中写入了相同的全部节点数据
   
6. # <wbr> <wbr> 3.目前只限于矩形截面的截面点，梁、壳单元只限于1个积分点，读入的材料只限于弹性属性
   
7. #***********************************************************************************************************
   
8. <span style="word-wrap: normal; line-height: 22.5px; font-size: 12px;">def FunWriteOdbData(Material,Section,Node,Frame,Area,LoadCase,ResultsPeriod,ResultsJointDispl,ResultsFrameSectionForce</span><wbr>,
   
9. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> ResultsAreaSectionForce,ResultsFrameStress,ResultsAreaStress,FileName):
   
10. <wbr> <wbr> <wbr> #-*-coding:UTF-8-*-
    
11. <wbr> <wbr> <wbr> #导入相应模块
    
12. <wbr> <wbr> <wbr> from odbAccess import *
    
13. <wbr> <wbr> <wbr> from odbMaterial import *
    
14. <wbr> <wbr> <wbr> from odbSection import *
    
15. <wbr> <wbr> <wbr> from abaqusConstants import *
    
16. <wbr> <wbr> <wbr> #----------------------------------------------------------------------
    
17. <wbr> <wbr> <wbr> #1.创建输出数据ODB，同时将创建根装配rootAssembly
    
18. <wbr> <wbr> <wbr> print('4...创建ABAQUS_ODB文件...')
    
19. <wbr> <wbr> <wbr> odb=Odb(name='resultModel',analysisTitle='ODB created with Python ODB API',description='',path=FileName)
    
20. <wbr> <wbr> <wbr>#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    
21. <wbr> <wbr> <wbr>#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    
22. <wbr> <wbr> <wbr> print('5...写入ODB模型数据信息...')
    
23. <wbr> <wbr> <wbr> #2.创建材料
    
24. <wbr> <wbr> <wbr> print('5.1写入材料数据')
    
25. <wbr> <wbr> <wbr> for MatName,MatType,E,u,G in Material:
    
26. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> material_object=odb.Material(name=MatName)
    
27. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> material_object.Elastic(type=ISOTROPIC,table=((E,u),))
    
28. <wbr> <wbr> <wbr>#===========================================================================================================================
    
29. <wbr> <wbr> <wbr> <span style="word-wrap: normal; line-height: 22.5px; font-size: 12px;">#3.创建截面</span>
    
30. <wbr> <wbr> <wbr> print('5.2写入截面数据')
    
31. <wbr> <wbr> <wbr> #a.梁单元截面
    
32. <wbr> <wbr> <wbr> beamSection={} <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #储存为字典类型，便于根据截面名字查找属性
    
33. <wbr> <wbr> <wbr> for ii in range(len(Section['SectionFrame'])):
    
34. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> SectName=Section['SectionFrame'][ii][0]
    
35. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> SectType=Section['SectionFrame'][ii][1]
    
36. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> SectMat =Section['SectionFrame'][ii][2]
    
37. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> if SectType==1: <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #I形截面
    
38. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> h=Section['SectionFrame'][ii][3] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #截面高度
    
39. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> l=h/2 <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #截面局部坐标系原点距底边距离
    
40. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> b1=Section['SectionFrame'][ii][4] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #下翼缘宽度
    
41. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> b2=Section['SectionFrame'][ii][5] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #上翼缘宽度
    
42. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> t1=Section['SectionFrame'][ii][6] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #下翼缘厚度
    
43. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> t2=Section['SectionFrame'][ii][7] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #上翼缘厚度
    
44. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> t3=Section['SectionFrame'][ii][8] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #腹板厚度
    
45. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> odb.IProfile(name=SectName,l=l,h=h,b1=b1,b2=b2,t1=t1,t2=t2,t3=t3)
    
46. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>beamSection[SectName]=odb.BeamSection(name=SectName,integration=BEFORE_ANALYSIS,profile=SectName,material=SectMat)
    
47. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> elif SectType==6: <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #箱形截面
    
48. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> a=Section['SectionFrame'][ii][3] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #截面宽度
    
49. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> b=Section['SectionFrame'][ii][4] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #截面高度
    
50. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> t1=Section['SectionFrame'][ii][5] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #翼缘厚度
    
51. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> t2=Section['SectionFrame'][ii][6] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #腹板厚度
    
52. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> t3=t1
    
53. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> t4=t2
    
54. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>odb.BoxProfile(name=SectName,uniformThickness=OFF,a=a,b=b,t1=t1,t2=t2,t3=t3,t4=t4)
    
55. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>beamSection[SectName]=odb.BeamSection(name=SectName,integration=BEFORE_ANALYSIS,profile=SectName,material=SectMat)
    
56. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> elif SectType==7: <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #圆管形截面
    
57. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> r=Section['SectionFrame'][ii][3] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #截面外半径
    
58. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> t=Section['SectionFrame'][ii][4] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #内壁厚度
    
59. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> odb.PipeProfile (name=SectName,r=r,t=t)
    
60. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>beamSection[SectName]=odb.BeamSection(name=SectName,integration=BEFORE_ANALYSIS,profile=SectName,material=SectMat)
    
61. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> elif SectType==8: <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #矩形截面
    
62. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> a=Section['SectionFrame'][ii][3]
    
63. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> b=Section['SectionFrame'][ii][4]
    
64. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> odb.RectangularProfile(name=SectName,a=a,b=b)
    
65. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>beamSection[SectName]=odb.BeamSection(name=SectName,integration=BEFORE_ANALYSIS,profile=SectName,material=SectMat)
    
66. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> elif SectType==9: <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #圆形截面
    
67. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> r=Section['SectionFrame'][ii][3] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #截面半径
    
68. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> odb.CircularProfile (name=SectName,r=r)
    
69. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>beamSection[SectName]=odb.BeamSection(name=SectName,integration=BEFORE_ANALYSIS,profile=SectName,material=SectMat)
    
70. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> else:
    
71. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> print('警告：缺少截面类型')
    
72. <wbr> <wbr> <wbr> #定义梁截面分类
    
73. <wbr> <wbr> <wbr> sCatRect=odb.SectionCategory(name='Rectangular',description='Space Beam 25 Section Points')
    
74. <wbr> <wbr> <wbr> sp1=sCatRect.SectionPoint(number=1,description='SectionPoint-1')
    
75. <wbr> <wbr> <wbr> sp2=sCatRect.SectionPoint(number=2,description='SectionPoint-5')
    
76. <wbr> <wbr> <wbr> sp3=sCatRect.SectionPoint(number=3,description='SectionPoint-21')
    
77. <wbr> <wbr> <wbr> sp4=sCatRect.SectionPoint(number=4,description='SectionPoint-25')
    
78. <wbr> <wbr> <wbr> #b.壳单元截面
    
79. <wbr> <wbr> <wbr> shellSection={}
    
80. <wbr> <wbr> <wbr> for SectName,SectMat,Thickness in Section['SectionShell']: <wbr>
    
81. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>shellSection[SectName]=odb.HomogeneousShellSection(name=SectName,material=SectMat,thickness=Thickness)
    
82. <wbr> <wbr> <wbr> #定义壳截面分类
    
83. <wbr> <wbr> <wbr> sCatShell=odb.SectionCategory(name='Shell',description='Five-Layered Shell')
    
84. <wbr> <wbr> <wbr> spBot=sCatShell.SectionPoint(number=1,description='Bottom')
    
85. <wbr> <wbr> <wbr> spMid=sCatShell.SectionPoint(number=3,description='Middle')
    
86. <wbr> <wbr> <wbr> spTop=sCatShell.SectionPoint(number=5,description='Top')
    
87. <wbr> <wbr> <wbr>#===========================================================================================================================
    

复制代码

旧言虐心

1. <wbr> <wbr> #4.创建部件
   
2. <wbr> <wbr> <wbr> print('5.3创建部件PART')
   
3. <wbr> <wbr> <wbr> beamPart=odb.Part(name='beamPart', embeddedSpace=THREE_D,type=DEFORMABLE_BODY)
   
4. <wbr> <wbr> <wbr> shellPart=odb.Part(name='shellPart', embeddedSpace=THREE_D,type=DEFORMABLE_BODY)
   
5. <wbr> <wbr> <wbr> #a.创建节点
   
6. <wbr> <wbr> <wbr> print(' <wbr> <wbr> a.写入节点数据')
   
7. <wbr> <wbr> <wbr> Node=sorted(Node,key=lambda x:x[0],reverse=False) <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #排序，节点号从小到大
   
8. <wbr> <wbr> <wbr> beamPart.addNodes(nodeData=tuple(Node),nodeSetName='ALL')
   
9. <wbr> <wbr> <wbr> shellPart.addNodes(nodeData=tuple(Node),nodeSetName='ALL')
   
10. <wbr> <wbr> <wbr> #为便于后面查找节点坐标，把节点数据储存为字典格式
    
11. <wbr> <wbr> <wbr> nodeData={}
    
12. <wbr> <wbr> <wbr> for NodeNum,x,y,z in Node:
    
13. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> nodeData[NodeNum]=(x,y,z)
    
14. <wbr> <wbr> <wbr> #----------------------------------------------------------------------------------
    
15. <wbr> <wbr> <wbr> print(' <wbr> <wbr> b.写入单元数据')
    
16. <wbr> <wbr> <wbr> #b.创建梁单元
    
17. <wbr> <wbr> <wbr> Frame=sorted(Frame,key=lambda x:x[0],reverse=False)
    
18. <wbr> <wbr> <wbr> beamData=[]
    
19. <wbr> <wbr> <wbr> for ii in range(len(Frame)):
    
20. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> beam_el_num=Frame[ii][0]
    
21. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> node1=Frame[ii][1]
    
22. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> node2=Frame[ii][2]
    
23. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> beamData.append((beam_el_num,node1,node2))
    
24. <wbr> <wbr> <wbr> beamData=tuple(beamData)
    
25. <wbr> <wbr> <wbr> beamPart.addElements(elementData=beamData, type='B31',elementSetName='BEAM',sectionCategory=sCatRect) <wbr> #单积分点
    
26. <wbr> <wbr> <wbr> #----------------------------------------------------------------------------------
    
27. <wbr> <wbr> <wbr> #c.创建壳单元
    
28. <wbr> <wbr> <wbr> Area=sorted(Area,key=lambda x:x[0],reverse=False)
    
29. <wbr> <wbr> <wbr> shellData=[]
    
30. <wbr> <wbr> <wbr> for ii in range(len(Area)):
    
31. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> shell=Area[ii][0]
    
32. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> nodes=Area[ii][1]
    
33. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> shellData.append((shell,nodes[0],nodes[1],nodes[2],nodes[3]))
    
34. <wbr> <wbr> <wbr> shellData=tuple(shellData)
    
35. <wbr> <wbr> <wbr> shellPart.addElements(elementData=shellData, type='S4R',elementSetName='SHELL',sectionCategory=sCatShell)#单积分点
    
36. <wbr> <wbr> <wbr>#===========================================================================================================================
    
37. <wbr> <wbr> <wbr> #5.创建装配件
    
38. <wbr> <wbr> <wbr> print('5.4创建装配件Assembly')
    
39. <wbr> <wbr> <wbr> beamInstance=odb.rootAssembly.Instance(name='BeamInstance',object=beamPart)
    
40. <wbr> <wbr> <wbr> shellInstance=odb.rootAssembly.Instance(name='ShellInstance',object=shellPart)
    
41. <wbr> <wbr> <wbr> #指定截面
    
42. <wbr> <wbr> <wbr> #a.梁单元截面
    
43. <wbr> <wbr> <wbr> print(' <wbr> <wbr> a.指定梁截面')
    
44. <wbr> <wbr> <wbr> for ii in range(len(Frame)):
    
45. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> beam_el_num=Frame[ii][0] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #单元编号
    
46. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> beam_el_sect=Frame[ii][3] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #截面名字
    
47. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #计算截面方向余弦
    
48. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> x1=nodeData[Frame[ii][1]][0]
    
49. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> y1=nodeData[Frame[ii][1]][1]
    
50. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> x2=nodeData[Frame[ii][2]][0]
    
51. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> y2=nodeData[Frame[ii][2]][1]
    
52. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> if abs(x2-x1)<1E-5:
    
53. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> cos1=-1
    
54. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> cos2=0
    
55. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> else:
    
56. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> cos1=-(y2-y1)/(x2-x1)
    
57. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> cos2=1
    
58. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> TEMPBEAMLEN=(cos1**2+cos2**2)**0.5
    
59. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> cos1=cos1/TEMPBEAMLEN
    
60. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> cos2=1/TEMPBEAMLEN
    
61. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> elset=beamInstance.ElementSetFromElementLab<wbr>els(name='BEAM_'+str(beam_el_num),elementLabels=(beam_el_num,)) #定义集合
    
62. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> beamInstance.assignSection(region=elset,section=beamSection[beam_el_sect]) <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #赋予截面
    
63. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> beamInstance.assignBeamOrientation(region=elset,method=N1_COSINES,vector=(cos1,cos2,0)) <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #截面方向
    
64. <wbr> <wbr> <wbr> #----------------------------------------------------------------------------------
    
65. <wbr> <wbr> <wbr> #b.壳单元截面
    
66. <wbr> <wbr> <wbr> print(' <wbr> <wbr> b.指定壳截面')
    
67. <wbr> <wbr> <wbr> for ii in range(len(Area)):
    
68. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> shell_el_num =Area[ii][0] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #单元编号
    
69. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> shell_el_sect=Area[ii][2] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #单元截面
    
70. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> elset=shellInstance.ElementSetFromElementLab<wbr>els(name='SHELL_'+str(shell_el_num),elementLabels=(shell_el_num,)) #定义集合
    
71. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> shellInstance.assignSection(region=elset,section=shellSection[shell_el_sect]) <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #赋予截面
    
72. <wbr> <wbr> <wbr>#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    
73. <wbr> <wbr> <wbr>#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
    
74. <wbr> <wbr> <wbr> print('6...写入ODB结果数据信息...')
    
75. <wbr> <wbr> <wbr> print('6.1创建分析步STEP')
    
76. <wbr> <wbr> <wbr> stepsObject={}
    
77. <wbr> <wbr> <wbr> for ii in range(len(LoadCase)):
    
78. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> CaseName=LoadCase[ii][0]
    
79. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> CaseType=LoadCase[ii][1]
    
80. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> if CaseType==1: <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #静力分析
    
81. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>stepsObject[CaseName]=odb.Step(name=CaseName,description='LINEAR_STATIC_ANALYSIS',domain=TIME,timePeriod=1.0)
    
82. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> elif CaseType==3: <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #模态分析
    
83. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>stepsObject[CaseName]=odb.Step(name=CaseName,description='MODAL_EIGEN_ANALYSIS',domain=MODAL)
    
84. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> elif CaseType==7: <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #时程分析
    
85. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> TimePeriod=LoadCase[ii][3]*LoadCase[ii][4]
    
86. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>stepsObject[CaseName]=odb.Step(name=CaseName,description='Time_HISTORY_ANALYSIS',domain=TIME,timePeriod=TimePeriod)
    
87. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> else:
    
88. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> print('警告：缺少分析工况')
    
89. <wbr> <wbr> <wbr>#===========================================================================================================================
    
90. <wbr> <wbr> <wbr> print('6.2写入场数据')
    
91. <wbr> <wbr> <wbr> for kk in range(len(LoadCase)):
    
92. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> CaseName=LoadCase[kk][0]
    
93. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> CaseType=LoadCase[kk][1]
    
94. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> JointDispl=ResultsJointDispl[CaseName]
    
95. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> JointDispl=sorted(JointDispl,key=lambda x:x[0],reverse=False) <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #按单元号从小到大排序 <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>
    
96. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> if CaseType==7: <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #%%%%%时程分析步%%%%%%
    
97. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> FrameSectionForce=ResultsFrameSectionForce<wbr>[CaseName]
    
98. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> AreaSectionForce=ResultsAreaSectionForce[CaseName]
    
99. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> FrameStress=ResultsFrameStress[CaseName]
    
100. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> AreaStress=ResultsAreaStress[CaseName]
     
101. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> FrameSectionForce=sorted(FrameSectionForce,key=lambda x:x[0],reverse=False)
     
102. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> AreaSectionForce=sorted(AreaSectionForce,key=lambda x:x[0],reverse=False)
     
103. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> FrameStress=sorted(FrameStress,key=lambda x:x[0],reverse=False)
     
104. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> AreaStress=sorted(AreaStress,key=lambda x:x[0],reverse=False)
     
105. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> if FrameSectionForce[0][1]>1:
     
106. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> print('警告：梁单元内力输出位置数目大于1')
     
107. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> if AreaStress[0][1]>1:
     
108. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> print('警告：壳单元内力输出位置数目大于1')
     
109. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> IncNum=LoadCase[kk][3] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #增量步总数
     
110. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> DT=LoadCase[kk][4] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #增量步时间间隔
     
111. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> for ii in range(IncNum):
     
112. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #创建帧
     
113. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>frameHistory=stepsObject[CaseName].Frame(incrementNumber=ii+1,frameValue=ii*DT,description='Increment <wbr> '+
     
114. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> str(ii+1)+':'+'Step Time = <wbr> '+str(ii*DT))
     
115. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #---------------------------------------------------------------------------------------------------------
     
116. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #节点位移
     
117. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> nodeLabelData=[] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #节点
     
118. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> dispUData=[] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #位移
     
119. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> dispURData=[] <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> #转角
     
120. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> for jj in range(len(JointDispl)):
     
121. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> nodeLabelData.append(JointDispl[jj][0])
     
122. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> U1=JointDispl[jj][1][ii]
     
123. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> U2=JointDispl[jj][2][ii]
     
124. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> U3=JointDispl[jj][3][ii]
     
125. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> R1=JointDispl[jj][4][ii]
     
126. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> R2=JointDispl[jj][5][ii]
     
127. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> R3=JointDispl[jj][6][ii]
     
128. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> dispUData.append((U1,U2,U3))
     
129. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> dispURData.append((R1,R2,R3))
     
130. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> uField = frameHistory.FieldOutput(name='U',description='U',type=VECTOR,validInvariants=(MAGNITUDE,))#创建位移场
     
131. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>uField.addData(position=NODAL,instance=beamInstance,labels=tuple(nodeLabelData),data=tuple(dispUData)) #写入位移
     
132. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>uField.addData(position=NODAL,instance=shellInstance,labels=tuple(nodeLabelData),data=tuple(dispUData))
     
133. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> urField = frameHistory.FieldOutput(name='UR',description='UR',type=VECTOR,validInvariants=(MAGNITUDE,))
     
134. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>urField.addData(position=NODAL,instance=beamInstance,labels=tuple(nodeLabelData),data=tuple(dispURData))
     
135. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>urField.addData(position=NODAL,instance=shellInstance,labels=tuple(nodeLabelData),data=tuple(dispURData))
     
136. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> stepsObject[CaseName].setDefaultField(uField) <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr>
     
137. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> else:
     
138. <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> <wbr> print('警告：缺少输出工况')
     
139. <wbr> <wbr> <wbr> #----------------------------------------------------------------------------------------------------------------
     
140. <wbr> <wbr> <wbr> #保存数据
     
141. <wbr> <wbr> <wbr> print('6.3...........完成................')
     
142. <wbr> <wbr> <wbr> odb.save()
     
143. <wbr> <wbr> <wbr> odb.close()

复制代码

参考文献：ABAQUS帮助文档