关于弯矩

aforce701

finish
/clear
b0=200
h0=300
l0=3000
ec=3.3e5
p0=0.2
/prep7
et,1,solid95
mp,ex,1,ec
mp,prxy,1,0.167
blc4,,,b0,h0,l0
/view,1,1,1,1
/ang,1
vplot
lsel,s,loc,x,1,b0-1
lsel,r,loc,y,0
dl,all,,uy
ksel,s,loc,x,0
ksel,r,loc,y,0
dk,all,ux
ksel,r,loc,z,0
dk,all,uz
asel,s,loc,y,h0
sfa,all,1,pres,p0
allsel,all
esize,50
vsweep,all
finish
/solu
solve
finish
/post1
pldisp
!------可以在下面添加语句完成(要求使用积分命令)
!A.跨中弯矩
!B.1/4跨剪力
经过摸索，讨论如下：
1.积分命令流如下：
path,kz,2
ppath,1,,b0/2,0,l0/2
ppath,2,,b0/2,h0,l0/2
pdef,sigmz,s,z
plpath,sigmz
!跨中轴力nfor1=-2.55e-5
pcalc,intg,nfor,sigmz,yg,b0
*get,nfor1,path,,last,nfor
!跨中剪力qfor1=-4.039e-9
pdef,sigmyz,s,yz
plpath,sigmyz
pcalc,intg,qfor,sigmyz,yg,1
*get,qfor1,path,,last,qfor
!跨中弯矩mfor1=-45251688误差0.6%
pcalc,mult,md1,sigmz,s !相乘计算
pcalc,intg,mfor,md1,s,b0 !积分计算
*get,mfor1,path,,last,mfor
我看了这么个例子。但是在求剪力的时候，pcalc,intg,qfor,sigmyz,yg,1
最后这个1让我看不懂了。是不是应该是b0呢？