ANSYS实体单元求截面内力和弯矩

yydryydr

　　实例：

1. 　　FINI
   
2. 
   
3. 　　/CLEAR
   
4. 
   
5. 　　/PREP7
   
6. 
   
7. 　　ET,1,SOLID95,
   
8. 
   
9. 　　MP,EX,1,2E11
   
10. 
    
11. 　　MP,PRXY,1,0.3
    
12. 
    
13. 　　BLC4,2,3,0.2,0.3,4
    
14. 
    
15. 　　DA,2,ALL !给编号为2的面施加面约束
    
16. 
    
17. 　　FK,3,FY,-1000 !给标号为3的关键点施加X方向的集中力为1000
    
18. 
    
19. 　　FK,1,FY,-1000 !给标号为4的关键点施加X方向的集中力为1000
    
20. 
    
21. 　　ESIZE,0.05,
    
22. 
    
23. 　　VMESH,ALL
    
24. 
    
25. 　　FINISH
    
26. 
    
27. 　　/SOLU
    
28. 
    
29. 　　SOLVE
    
30. 
    
31. 　　FINI
    

复制代码

　　需要提取实体矩形梁中间截面的内力和弯矩。
　　内力很好说，主要是弯矩，实体单元因为没有转动自由度，无法直接求出弯矩MX、MY和MZ，需要通过指定截面的中心，根据节点力对该中心进行弯矩求解。

　　指定弯矩求解中心的命令如下：



　　可以通过指定节点号的方式，比如：spoint,100
　　也可以通过给定具体的坐标，该坐标针对总体笛卡尔坐标系而言，比如，spoint,,1.5,2.1,3.5
　　因此，在本次实例中，求解结束后，输入以下命令：
　　ALLS,ALL
　　NSEL,S,LOC,Z,2-0.05,2 !选择跨中截面及截面右侧的节点和单元
　　ESLN,,1 !通过节点选择单元
　　NSEL,R,LOC,Z,2 !从中再选择跨中截面的节点
　　SPOINT,,2.1,3.15,2 !指定力矩求和中心(跨中截面的中心)
　　FSUM !单元节点力求和，并给出列表结果

　　最终得到以下结果



　　此时显式此时的节点和单元，ansys界面分别显示：





　　这里注意一个问题，选择截面节点的时候，不能只是选择截面的节点。否则结果可能会特别小。
　　比如直接输入命令：
　　ALLS,ALL
　　NSEL,S,LOC,Z,2 !从中再选择跨中截面的节点
　　SPOINT,,2.1,3.15,2 !指定力矩求和中心(跨中截面的中心)
　　FSUM

　　得到的结果如下：



　　同样显示此时的单元和节点，ansys界面显示：





　　很明显，结果不一样，单元不一样，这个为什么会这样，还需要去了解具体算法才能知道。


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

think2015

感谢楼主分享，回去试试