实体单元弯矩轴力提取

实体单元弯矩和轴力的提取
对于一个实体梁
长10m，宽0.6m，高1.0m
程序代码
finish$/clear
finish$/prep7
blc4,0,0,10,0.6,1
ET,1,SOLID45$MP,EX,1,200e6$MP,PRXY,1,0.3!定义单元类型、材料特性
ESIZE,0.5$VMESH,ALL
FINISH$/SOLU
asel,s,loc,x,0
nsla,s,1$d,all,all
nsel,s,loc,x,10
*get,Nnode,node,0,count
f,all,fx,100/Nnode
f,all,fy,30/Nnode
f,all,fz,50/Nnode
SOLVE
finish$/post1
! define surface
wpcsys,-1
wpoffs,5,0.3,0.5$wprota,0,0,90
sucr,xsurf1,cplane,3$!supl,xsurf1
!define varible needed
sumap,sx,s,x$!supl,xsurf1,sx
sumap,sy,s,y$!supl,xsurf1,sy
sumap,sz,s,z$!supl,xsurf1,sz
sumap,sxy,s,xy$!supl,xsurf1,sxy
sumap,syz,s,yz$!supl,xsurf1,syz
sumap,sxz,s,xz$!supl,xsurf1,sxz
!--------The following code varies in the normal direction of section definded------
!------------------------This is for the direction of x aixs------------------
! get section force: Ax SFy SFz
sueval,Ax,sx,intg
sueval,SFy,sxy,intg
sueval,SFz,sxz,intg
! get section center
sueval,A,DA,sum! GCX, GCY, GCZ---global Cartesian coordinates at each point on the surface.
sueval,Sx,GCx,intg$dx=Sx/A
sueval,Sy,GCy,intg$dy=Sy/A
sueval,Sz,GCz,intg$dz=Sz/A
! get the section moment: Mz,My,TOR
sucalc,sxdz,sx,mult,GCz$!supl,xsurf1,sxdz
sueval,Mz,sxdz,intg
Mz=Mz-Ax*dz!Mz
sucalc,sxdy,sx,mult,GCy$!supl,xsurf1,sxdy
sueval,My,sxdy,intg
My=My-Ax*dy!My
sucalc,sxydz,sxy,mult,GCz
sucalc,sxzdy,sxz,mult,GCy
sucalc,tordl,sxydz,add,sxzdy
sueval,TOR,tordl,intg
TOR=TOR-SFy*dz-SFz*dy!TOR
!------------------------------------------------------------------------------------
Ax=∫Sx
A
SFy=∫Sxy
A
SFz=∫Sxz
A
My=∫Sxdz
A
Mz=∫Sxdy
A
TOR=∫Sxydz−Sxzdy
A
!-----------------------This is for the direction of x aixs--------------------------
! get section force: Ax SFy SFz
sueval,Ax,sy,intg
sueval,SFy,sxy,intg
sueval,SFz,syz,intg
! get section center
sueval,A,DA,sum! GCX, GCY, GCZ---global Cartesian coordinates at each point on the surface.
sueval,Sx,GCx,intg$dx=Sx/A
sueval,Sy,GCy,intg$dy=Sy/A
sueval,Sz,GCz,intg$dz=Sz/A
! get the section moment: Mz,My,TOR
sucalc,sxdz,sx,mult,GCz$!supl,xsurf1,sxdz
sueval,My,sxdz,intg
My=My-Ax*dz!My
sucalc,sxdy,sx,mult,GCy$!supl,xsurf1,sxdy
sueval,My,sxdy,intg
Mz=Mz-Ax*dy!Mz
sucalc,sxydz,sxy,mult,GCz
sucalc,sxzdy,sxz,mult,GCy
sucalc,tordl,sxydz,add,sxzdy
sueval,TOR,tordl,intg
TOR=TOR-SFy*dz-SFz*dy!TOR
!-----------------------------------------------------------------------------------------。