上面的附件有点问题，在施加节点荷载之前没有选中节点。我重新贴一下吧。

/CLEAR,START
/prep7

ET,1,SOLID185

MP,EX,1,3.0E10         !材料1弹模
MP,DENS,1,1800       !材料1密度
MP,PRXY,1,0.25         !材料1泊松比

BLOCK,0,2,0,2,0,40

VSEL,ALL
VATT,1,,1,0

MSHAPE,0,3D
MSHKEY,1         !选择映射划分方式

LSEL,ALL
ESIZE,0.1

VSEL,ALL
VMESH,ALL
FINISH

/solu         !进入求解

!底部固定
NSEL,ALL
NSEL,s,loc,z,0
D,all,all,0
!顶部施加水平集中力
NSEL,ALL
F,463,FX,5000
F,443,FX,5000

allsel         !选择所有
solve
finish

/post1
*do,ii,1,21
WPCSYS,-1
WPOFFS,1,1,(ii-1)*2
SUCR,suz%ii%,cplane    !面的建立

SUMAP,mysx,s,x
SUMAP,mysy,s,y
SUMAP,mysz,s,z
SUMAP,mysxy,s,xy
SUMAP,mysyz,s,yz
SUMAP,mysxz,s,xz

SUEVAL,xforce,mysxz,intg    !求截面上FX
SUEVAL,yforce,mysyz,intg    !求截面上FY
SUEVAL,zforce,mysz,intg      !求截面上FZ

SUEVAL,mya,da,sum     !求截面面积并赋给变量mya
SUEVAL,myya,gcy,intg  !求对X轴的面积矩并赋给变量myya
myya=myya/mya         !得到面积重心到X轴的距离=面积矩/面积
SUEVAL,myxa,gcx,intg  !求对Y轴的面积矩并赋给变量myxa
myxa=myxa/mya         !得到面积重心到Y轴的距离=面积矩/面积

SUCALC,szgcy,mysz,mult,gcy    !计算mysz乘以gcy，并赋给面项szgcy
SUEVAL,mx1,szgcy,intg   !对面项szgcy在面上积分，得到X方向弯矩mx1

SUCALC,szgcx,mysz,mult,gcx    !计算mysz乘以gcx，并赋给面项szgcx
SUEVAL,my1,szgcx,intg   !对面项szgcx在面上积分，得到Y方向弯矩my1

*cfopen,result%ii%,txt,,APPEND   !输出计算结果
*VWRITE,xforce,yforce,zforce,mx1,my1
(5E16.8)

SUDEL,suz%ii%

*enddo