实体单元弯矩轴力提取

合集下载

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

实体单元弯矩和轴力的提取

对于一个实体梁

长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--------------------------