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]