ANSYS实体单元求截面内力和弯矩
实例：FINI
/CLEAR
/PREP7
ET,1,SOLID95,
MP,EX,1,2E11
MP,PRXY,1,0.3
BLC4,2,3,0.2,0.3,4
DA,2,ALL !给编号为2的面施加面约束
FK,3,FY,-1000 !给标号为3的关键点施加X方向的集中力为1000
FK,1,FY,-1000 !给标号为4的关键点施加X方向的集中力为1000
ESIZE,0.05,
VMESH,ALL
FINISH
/SOLU
SOLVE
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
感谢楼主分享，回去试试
页:
[1]