$add,4,4,5$prod,5,3,,,,,,1/1000000
$/axlab,x,strain-eqv$/axlab,y,sigx-eqv(MPa) $PRVAR,4,5
! 剪应变与梁端荷载曲线
$Cmsel,s,gzg$vsel,r,loc,z,0,0.0045$vsel,r,loc,y,-0.125+0.006,0.125-0.006$eslv,s$nsle,s /POST26 NUMVAR,200$ $RFORCE,3,7150,F,Y $esol,4,6351,2661,epel,xy $esol,5,6351,2661,eppl,xy $add,4,4,5$prod,5,3,,,,,,-1/1000 $PRVAR,4,5 /post1
$Cmsel,s,gzg$vsel,r,loc,z,0,0.0045$vsel,r,loc,y,-0.125+0.006,0.125-0.006$$vsel,r,loc,x,-0.175,0.175$eslv,s$nsle,s
SET,LAST$SET,PREVIOUS
*GET,RTIME,ACTIVE,0,SET,TIME !获取最终收敛时间 Pldisp$plnsol,s,x !显示变形图与X方向应力分布图 Plnsol,epel,xy !Xy弹性应变分布图
Plnsol,eppl,xy !xy塑性应变分布图 Plnsol,epto,xy !Xy方向总应变分布图 /post1 Allsel,all
!腹板剪力 $Cmsel,s,gzg $Nslv,s,1
$Nsel,r,loc,x,-0.175,0.175 $Nsel,r,loc,z,-0.0045,0.0045 $Nsel,r,loc,y,0,0.05 $Esln,s,1 $Nsel,r,loc,y,0 $fsum !内混凝土剪力 Allsel,all
$Cmsel,s,zhuzi $Nslv,s,1
$Nsel,r,loc,x,-0.175,0.175 $Nsel,r,loc,z,-0.0625,0.0625 $Nsel,r,loc,y,0,0.05 $Esln,s,1 $Nsel,r,loc,y,0 $fsum /post1 Allsel,all
!外混凝土剪力 $Cmsel,s,zhuzi $Nslv,s,1
$Nsel,r,loc,x,-0.175,0.175
$Nsel,u,loc,z,-0.0625+0.001,0.0625-0.001 $Nsel,r,loc,y,0,0.05 $Esln,s,1 $Nsel,r,loc,y,0 $fsum !混凝土总剪力 Allsel,all
$Cmsel,s,zhuzi $Nslv,s,1
$Nsel,r,loc,x,-0.175,0.175 $Nsel,r,loc,y,0,0.05 $Esln,s,1 $Nsel,r,loc,y,0 $fsum
!总剪力 /post1
Allsel,all
$Nsel,s,loc,x,-0.175,0.175 $Nsel,r,loc,y,0,0.05 $Esln,s,1 $Nsel,r,loc,y,0 $fsum