ANSYS中实体单元积分法求截面内力

[code]!!!悬臂梁中间截面的内力（反力和弯矩）计算。材料属性分别为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 !单元节点力求和，并给出列表结果

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