求助！！！

finish
/clear
/filname,sandwichpanel
/title,STATIC ANALYSIS OF SANDWICH PANEL
/prep7
b=3000!!板宽
h=3000!!板高
t=60!!板厚
Ex1=2.4e4!!C30混凝土的弹性模量
Ex2=2.06e11!!Q235钢以及拉结件的弹性模量
Ex3=4e6!!eps板弹性模量
Nu1=0.2!!混凝土泊松比
Nu2=0.3!!钢材以及拉结件的泊松比
Nu3=0.05!!保温板泊松比
Fe=4e5!!保温板强度
Fc=25!!混凝土的抗压强度
Ft=1.43e6!!混凝土的抗拉强度
Fy=235e6!!钢材的屈服强度
NN=20!!每条钢筋的分段数
d=8!!钢筋直径
L=2940!!钢筋长度
PI=3.1415926535
Area=PI*d*d/4!!单根钢筋的面积
dl=20!!拉结件直径
s=PI*dl*dl/4!!拉结件面积
rou1=1600!!混凝土密度
rou2=7800!!钢筋密度
rou3=500!!eps板密度
pres=1000!!风压
dL=L/NN!!每段钢筋长度
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************混凝土材料属性************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
et,1,solid65
keyopt,1,1,0!!考虑大变形
keyopt,1,5,1!!给出每个积分点的解
keyopt,1,6,3!!给出积分点的解
keyopt,1,7,1!!考虑应力松弛，有助于收敛
r,1!!混凝土单元实常数
mp,ex,1,Ex1
mp,prxy,1,Nu
mp,dens,1,rou1
tb,concr,1,1,9
tbdata,,0.35,0.75,3.1125,-1!!混凝土材料的相关参数，关闭压碎
tb,miso,1,,15
tbpt,,0.0002,4.8
tbpt,,0.0004,9.375
tbpt,,0.0006,13.51
tbpt,,0.0008,17.02
tbpt,,0.0010,19.83
tbpt,,0.0012,21.95
tbpt,,0.0014,23.43
tbpt,,0.0016,24.365
tbpt,,0.0018,24.856
tbpt,,0.0020,FC
tbpt,,0.0038,FC
tbplot
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************保温材料属性************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
et,4,solid186
mp,ex,4,Ex3
mp,prxy,4,Nu3
mp,dens,4,rou3
tb,bkin,4
tbdata,,Fe
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************钢筋及拉结件材料属性************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
et,2,link180
mp,dens,2,rou2
mp,ex,2,Ex2
mp,prxy,2,Nu2
tb,bkin,2
tbdata,,Fy
r,2,Area
et,3,link180
mp,dens,3,rou2
mp,ex,3,Ex2
mp,prxy,3,Nu2
tb,bkin,3
tbdata,,Fy
r,3,s
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************定义节点************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
*do,i,1,NN+1
*do,j,1,NN+1
n,(i-1)*(NN+1)+j,30+(j-1)*dL,30+(i-1)*dL,30
*enddo
*enddo
*do,i,1,NN+1
*do,j,1,NN+1
n,(i-1)*(NN+1)+j+441,30+(j-1)*dL,30+(i-1)*dL,170
*enddo
*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************定义单元************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!-----------------------钢筋----------------------
mat,2
type,2
real,2
nsel,s,,,1,(NN+1)**2,1
*do,i,1,NN+1
*do,j,1,NN
*if,j,eq,(NN+1)*i,then
e,1,2
*else
e,(i-1)*(NN+1)+j,(i-1)*(NN+1)+j+1
*endif
*enddo
*enddo
*do,i,1,NN
*do,j,1,NN+1
*if,(i-1)*(NN+1)+j,eq,420+i,then
e,1,2
*else
e,(i-1)*(NN+1)+j,(i)*(NN+1)+j
*endif
*enddo
*enddo
*do,i,1,NN+1
*do,j,1,NN
*if,j,eq,(NN+1)*i,then
e,1,2
*else
e,(i-1)*(NN+1)+j+441,(i-1)*(NN+1)+j+1+441
*endif
*enddo
*enddo
*do,i,1,NN
*do,j,1,NN+1
*if,(i-1)*(NN+1)+j,eq,420+i,then
e,1,2
*else
e,(i-1)*(NN+1)+j+441,(i)*(NN+1)+j+441
*endif
*enddo
*enddo
!!!!!!!!!!!!!!!!!!!-----------------------拉结件----------------------
type,3
real,3
mat,3
e,89,(NN+1)**2+89
e,95,(NN+1)**2+95
e,101,(NN+1)**2+101
e,215,(NN+1)**2+215
e,221,(NN+1)**2+221
e,227,(NN+1)**2+227
e,341,(NN+1)**2+341
e,347,(NN+1)**2+347
e,353,(NN+1)**2+353
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************建立几何实体************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
blc4,0,0,b,h,t
wpave,0,0,140
blc4,0,0,b,h,t
wpave,0,0,60
blc4,0,0,b,h,80
wpcsys,-1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************划分几何实体************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
wpave,0,0,30
vsbw,all
wpave,0,0,170
vsbw,all
wpave,0,0,0
wprota,0,90,0
wpave,0,30,0
vsbw,all
wpave,0,2970,0
vsbw,all
wpave,0,30,0
*do,i,1,NN-1
wpoff,0,0,-dL
vsbw,all
*enddo
wpave,0,0,30
wprota,0,0,90
wpave,30,0,0
vsbw,all
wpave,2970,0,0
vsbw,all
wpave,30,0,30
*do,i,1,NN-1
wpoff,0,0,dL
vsbw,all
*enddo
wpcsys,-1!!工作平面复原
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************划分几何边线长度，以便划分生成有限元模型************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
lsel,s,loc,z,30,30
lsel,a,loc,z,170,170
lsel,r,loc,y,30,2970,dL
lesize,all,t*0.5
lsel,s,loc,z,30,30
lsel,a,loc,z,170,170
lsel,r,loc,x,30,2970,dL
lesize,all,t*0.5
lsel,s,loc,z,60,60
lsel,a,loc,z,140,140
lsel,a,loc,z,0,0
lsel,a,loc,z,200,200
lsel,r,loc,x,30,2970,dL
lesize,all,t*0.5
lsel,s,loc,z,60,60
lsel,a,loc,z,140,140
lsel,a,loc,z,0,0
lsel,a,loc,z,200,200
lsel,r,loc,y,30,2970
lesize,all,t*0.5
lsel,all
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************划分生成有限元模型************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
vsel,s,loc,z,0,60
vatt,1,1,1
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,loc,z,140,200
vatt,1,1,1
mshape,0,3d
mshkey,1
vmesh,all
vsel,s,loc,z,60,140
vatt,4,1,1
mshape,0,3d
mshkey,1
vmesh,all
vsel,all
allsel
nummrg,all!!合并
numcmp,all!!压缩
eplot!!重新显示
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!************************求解************************!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
/solu
!!施加约束
lsel,s,loc,y,h
lsel,r,loc,z,0
lsel,r,loc,x,0,b,dL
dl,all,,ux,0
dl,all,,uy,0
dl,all,,uz,0
lsel,all
lsel,s,loc,y,0
lsel,r,loc,z,0
lsel,r,loc,x,0,b,dL
dl,all,,uz,0
lsel,all
dtran
!!施加荷载
!asel,s,loc,z,200,200
!sfa,all,,pres,1.6*120/1000*pres
acel,0,0,-9.8
!sftran
allsel
antype,0
pstres,on
nlgeom,on
nsubst,1e3
outres,all,all
autots,1!!自动时间
lnsrch,1!!线性搜索
cnvtol,f,,0.05,2,0.5!!力收敛准则
arclen,on
allsel
ncnv,2!!求解不收敛仍不跳过
solve
finish