ANSYS中实体单元计算截面弯矩的方法

!!!!悬臂梁中间截面的内力（反力和弯矩）计算。材料属性分别为EX,1,2E11，
!!!!PRXY,1,0.3;单元类型为实体单元，长、宽、高分别为4、0.2、0.3m，一端完全固定，另一端的两个
!!!!上端点处，分别作用有两个大小为1000N的集中荷载，试求跨中截面的反力及弯矩。
!!!!
!!!!采用单元节点力求和法，通过选择节点和单元，然后对单元节点力求和即得到某个
!!!!某个截面的内力。但该法需要所求内力的截面为一列单元的边界，或者说截面
!!!!不穿过单元（节点分布在截面上），这样所求的截面内力是精确的。
!!!!
                          
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
/POST1
NSEL,S,LOC,Z,2-0.05,2    !选择跨中截面及截面右侧的节点和单元
ESLN,,1                  !通过节点选择单元      
NSEL,R,LOC,Z,2           !从中再选择跨中截面的节点  
SPOINT,,2.1,3.15,2       !指定力矩求和中心（跨中截面的中心）
FSUM                     !单元节点力求和，并给出列表结果

      计算结果与手算结果完全一致，说明是正确的，这种方法叫做“单元节点力求和法”，还有另外两种方法，分别叫做“截面分块积分法”“面操作法”，感觉还是单元节点力求和法简单方便。

      单元节点力求和法和板单元中的方法相同，即通过选择节点和单元，然后对单元节点力求和即可得到某个截面的内力。但该法需要所求内力的截面为一列单元的边界，或者说截面不穿过单元，这样所求截面内力是精确的。

      对于复杂结构，由于单元划分控制不可能那么好，就不如面操作准确，除非在划分单元时就决定求解内力的截面，然后将几何实体在此位置切分