Abaqus利用region读取结果数据

delay

　　Abaqus计算完成后，为输出某部分模型的结果数据，需要用到region，并且配套使用的是getSubset函数。

　　示例一创建一个变量center，将一个叫做PUNCH的节点set赋值给这个变量，通过整个模型的位移场displacement得到该节点set的输出数据，代码如下：

1. 　　center = odb.rootAssembly.instances['PART-1-1'].nodeSets['PUNCH']
   
2. 
   
3. 　　centerDisplacement = displacement.getSubset(region=center)
   
4. 
   
5. 　　centerValues = centerDisplacement.values
   
6. 
   
7. 　　for v in centerValues:
   
8. 
   
9. 　　print v.nodeLabel, v.data
   

复制代码

　　其中的displacement是已经创建好的位移场数据，centerDisplacement变量得到该节点集的所有数据，centerValues变量得到该节点集的数值数据，通过循环输出节点号和当前数据。


　　示例二如下：

1. 　　#创建变量topCenter，得到单元集CENT
   
2. 
   
3. 　　topCenter =odb.rootAssembly.instances['PART-1-1'].elementSets['CENT']
   
4. 
   
5. 　　#创建变量stressField，得到位移场数据
   
6. 
   
7. 　　stressField = odb.steps['Step-2'].frames[3].fieldOutputs['S']
   
8. 
   
9. 　　#创建变量field，得到该单元集的积分点的数据
   
10. 
    
11. 　　field = stressField.getSubset(region=topCenter,
    
12. 
    
13. 　　position=INTEGRATION_POINT, elementType = 'CAX4')
    
14. 
    
15. 　　#创建变量fieldValues，得到field的数值数据
    
16. 
    
17. 　　fieldValues = field.values
    
18. 
    
19. 　　#循环输出单元号，单元上的积分点号以及积分点的应力分量
    
20. 
    
21. 　　for v in fieldValues:
    
22. 
    
23. 　　print 'Element label = ', v.elementLabel,
    
24. 
    
25. 　　if v.integrationPoint:
    
26. 
    
27. 　　print 'Integration Point = ', v.integrationPoint
    
28. 
    
29. 　　else:
    
30. 
    
31. 　　print
    
32. 
    
33. 　　# For each tensor component.
    
34. 
    
35. 　　for component in v.data:
    
36. 
    
37. 　　# Print using a format. The comma at the end of the
    
38. 
    
39. 　　# print statement suppresses the carriage return.
    
40. 
    
41. 　　print '%-10.5f' % component,
    
42. 
    
43. 　　# After each tuple has printed, print a carriage return.
    
44. 
    
45. 　　Print
    

复制代码

　　最终的输出结果：

1. 　　Element label = 1 Integration Point = 1
   
2. 
   
3. 　　S : 0.01230 -0.05658 0.00892 -0.00015
   
4. 
   
5. 　　Element label = 1 Integration Point = 2
   
6. 
   
7. 　　S : 0.01313 -0.05659 0.00892 -0.00106
   
8. 
   
9. 　　Element label = 1 Integration Point = 3
   
10. 
    
11. 　　S : 0.00619 -0.05642 0.00892 -0.00023
    
12. 
    
13. 　　Element label = 1 Integration Point = 4
    
14. 
    
15. 　　S : 0.00697 -0.05642 0.00892 -0.00108
    
16. 
    
17. 　　Element label = 11 Integration Point = 1
    
18. 
    
19. 　　S : 0.01281 -0.05660 0.00897 -0.00146
    
20. 
    
21. 　　Element label = 11 Integration Point = 2
    
22. 
    
23. 　　S : 0.01183 -0.05651 0.00897 -0.00257
    
24. 
    
25. 　　Element label = 11 Integration Point = 3 ...
    

复制代码

　　另：对于getSubset的position有以下几个可选：
　　INTEGRATION_POINT !积分点
　　NODAL !节点(平均)
　　ELEMENT_NODAL !单元节点(未平均)

　　CENTROID !单元质心点


转自：http://blog.sina.com.cn/s/blog_6465f2ed0102x4eg.html