/prep7
/title,"ARCH"

*dim,xlx,,36
xlx(1)=15.00,14.50,14.00,13.50,13.00,12.50,12.00,11.50,11.00,10.50
xlx(11)=10.00,9.50,9.00,8.50,8.00,7.50,7.00,6.50,6.00,5.50
xlx(21)=5.00,4.50,4.00,3.50,3.00,2.50,2.00,1.50,1.00,0.50
xlx(31)=0.00

*do,i,1,31,1

angle=atan(0.2597*sinh(1.573*xlx(i)/15.00))
h=0.85
z1=-xlx(i)-(h/2)*sin(angle)
y1=-(3.75/1.514)*(cosh(1.573*(-xlx(i))/15.00)-1)+(h/2)*cos(angle)
z2=-xlx(i)+(h/2)*sin(angle)
y2=-(3.75/1.514)*(cosh(1.573*(-xlx(i))/15.00)-1)-(h/2)*cos(angle)
k,i,0,y1,z1
k,31+i,0,y2,z2
*enddo

*do,i,1,28,3
a,i,i+1,i+2,i+3,i+34,i+33,i+32,i+31
*enddo

arsym,z,1,10,1,100,0,0
allsel
nummrg,all

k,1001,0,0,20
k,1002,4,0,20
l,1001,1002

vdrag,1,2,3,4,5,6,139
vdrag,7,8,9,10,,,139
vdrag,11,12,13,14,15,16,139
vdrag,17,18,19,20,,,139
vdrag,7,14,,,,,139

ldel,139,,,1
allsel
nummrg,all

/view,1,-1,1,-1
allsel
vplo

et,1,plane42
et,2,solid65
mp,ex,2,3.075e10
mp,prxy,2,0.255
mp,dens,2,2780
tb,mkin,2
tbtemp,,strain
tbdata,,10e-4,15e-4,17.5e-4,20e-4,30e-4
tbtemp,,
tbdata,,30.75e6,38.44e6,40.36e6,41e6,34.8e6
r,1
tb,concr,2
tbtemp,,
tbdata,,0.20,0.8,4e6,41e6

asel,s,,,98
aplo
lesize,68,,,4
lesize,272,,,4
lesize,271,,,8
lesize,273,,,8
mshkey,1
amesh,98

allsel
lsel,s,loc,x,0
lsel,a,loc,x,4
lsel,u,,,4,12,4
lsel,u,,,19,75,7
lsel,u,,,79
lsel,u,,,83,132,7
lsel,u,,,150,158,8
lsel,u,,,165,217,13
lsel,u,,,231,246,15
lsel,u,,,259,272,13
lsel,u,,,279
lsel,u,,,288,353,13
lsel,u,,,369,395,13
lsel,u,,,426
lesize,all,,,1

allsel
type,2
mat,2
real,1
vsweep,all
aclear,all

/plopts,date,0
!/triad,ltop
/view,1,-1,1,-1
eplo
/solu
allsel
acel,,9.81
asel,s,,,28
asel,a,,,103
nsla,,1
d,all,all
fini

/solu
allsel
antype,static

nropt,full,,1
autots,on
pred,on

neqit,50
kbc,0
cnvtol,f,100000,0.001,2
cnvtol,u,10,0.001,0

time,1
nsubst,50
allsel
outpr,all,all
outres,all,all
solve

time,2
nsubst,50
!L/4加集中力荷载
a=20

allsel
nsel,s,loc,z,-7.50
F,all,fy,-100000*a
allsel
outpr,all,all
outres,all,all
solve
fini

说明:
建了一片拱肋进行极限承载力计算,模型没有问题,进行静力计算都好的,但是进行非线形分析的时候出问题了,请高手调试一下

qq:13108019

十分感谢,请帮忙修改一下