附多荷载步文件命令流，为一型钢混凝土悬挑梁

bjtu2010

我的想法是：采用荷载步文件法

1，荷载步文件1：沿斜向型钢轴向拉力3000KN，同时竖向施加压力0KN，均为斜坡荷载。

2，荷载步文件2：沿斜向型钢轴向拉力3000KN，同时竖向施加压力200KN，均为斜坡荷载。

3，荷载步文件3：沿斜向型钢轴向拉力3000KN，同时竖向施加压力400KN，均为斜坡荷载。

4，荷载步文件4：沿斜向型钢轴向拉力3000KN，同时竖向施加压力600KN，均为斜坡荷载。

5，荷载步文件5：沿斜向型钢轴向拉力3000KN，同时竖向施加压力800KN，均为斜坡荷载。

观察竖向压力对梁受力的影响，并画出各荷载位移曲线在同一图中作比较。以下的命令流正确吗？应该怎么修改？
我算了一次，最后到第五个文件的第6子步时就退出了，怎么只提取出文件1的所有计算结果以用来画文件1的荷载位移图？


!第一个荷载步文件
vsel,s,loc,x,-200,0  !选出垫块
wplane,,0,-(h1/2-wh/2),0,l,-(h2/2-wh/2),0,0,-(h1/2-wh/2),-b/2$wpoff,-50,,  !施加拉力
asel,s,loc,x,0
sfa,all,,pres,-p1$wpstyl,defa$allsel,all
time,1
lnsrch,on
pred,on
nsubst,50,10000,10
outres,all,all
autos,on
neqit,50
cnvtol,u,,0.05  
kbc,0
lswrite,1
!第二个荷载步文件
p2=fc2/(50*w/(l/sqrt(((h2-h1)/2)**2+l**2)))
vsel,s,loc,x,-200,0  !选出垫块
wplane,,0,-(h1/2-wh/2),0,l,-(h2/2-wh/2),0,0,-(h1/2-wh/2),-b/2$wpoff,-50,,  !施加拉力
asel,s,loc,x,0
sfa,all,,pres,-p1$wpstyl,defa$allsel,all

vsel,s,loc,x,-200,0  !选出垫块
asel,s,loc,y,-(h1/2-wh/2) !施加压力
sfa,all,,pres,p2$allsel,all

time,1
lnsrch,on
pred,on
nsubst,50,10000,10
outres,all,all
autos,on
neqit,50
cnvtol,u,,0.05  
kbc,0
lswrite,2
!第三个荷载步文件
p2=fc3/(50*w/(l/sqrt(((h2-h1)/2)**2+l**2)))
vsel,s,loc,x,-200,0  !选出垫块
wplane,,0,-(h1/2-wh/2),0,l,-(h2/2-wh/2),0,0,-(h1/2-wh/2),-b/2$wpoff,-50,,  !施加拉力
asel,s,loc,x,0
sfa,all,,pres,-p1$wpstyl,defa$allsel,all

vsel,s,loc,x,-200,0  !选出垫块
asel,s,loc,y,-(h1/2-wh/2) !施加压力
sfa,all,,pres,p2$allsel,all

time,1
lnsrch,on
pred,on
nsubst,50,10000,10
outres,all,all
autos,on
neqit,50
cnvtol,u,,0.05  
kbc,0
lswrite,3
!第四个荷载步文件
p2=fc4/(50*w/(l/sqrt(((h2-h1)/2)**2+l**2)))
vsel,s,loc,x,-200,0  !选出垫块
wplane,,0,-(h1/2-wh/2),0,l,-(h2/2-wh/2),0,0,-(h1/2-wh/2),-b/2$wpoff,-50,,  !施加拉力
asel,s,loc,x,0
sfa,all,,pres,-p1$wpstyl,defa$allsel,all

vsel,s,loc,x,-200,0  !选出垫块
asel,s,loc,y,-(h1/2-wh/2) !施加压力
sfa,all,,pres,p2$allsel,all

time,1
lnsrch,on
pred,on
nsubst,50,10000,10
outres,all,all
autos,on
neqit,50
cnvtol,u,,0.05  
kbc,0
lswrite,4
!第五个荷载步文件
p2=fc5/(50*w/(l/sqrt(((h2-h1)/2)**2+l**2)))
vsel,s,loc,x,-200,0  !选出垫块
wplane,,0,-(h1/2-wh/2),0,l,-(h2/2-wh/2),0,0,-(h1/2-wh/2),-b/2$wpoff,-50,,  !施加拉力
asel,s,loc,x,0
sfa,all,,pres,-p1$wpstyl,defa$allsel,all

vsel,s,loc,x,-200,0  !选出垫块
asel,s,loc,y,-(h1/2-wh/2) !施加压力
sfa,all,,pres,p2$allsel,all

time,1
lnsrch,on
pred,on
nsubst,50,10000,10
outres,all,all
autos,on
neqit,50
cnvtol,u,,0.05  
kbc,0
lswrite,5

fini
!求解控制设置和求解
/solu
allsel,all                                              !求解之前全选一下
lssolve,1,5,1