yydryydr 发表于 2016-5-23 14:59

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

think2015 发表于 2016-5-24 13:25

感谢楼主分享,回去试试
页: [1]
查看完整版本: ANSYS实体单元求截面内力和弯矩