finish
/clear
/prep7
l0=1000
xh=400
area=1.84e-2
em=2.06e7
q0=65
qf=10000
h0=9000
err0=1/1000
enum=60
istran=0.9
et,1,link10
r,1,area,istran
mp,ex,1,em
mp,prxy,1,0.3
mp,dens,1,q0/area
!在弦线位置创建模型，施加自重荷载
k,1,l0
k,2,,xh
l,1,2
lesize,all,,,enum
lmesh,all
d,node(l0,0,0),all
d,node(0,xh,0),all
node1=nelem(enum/2,1)
node2=nelem(enum/2,2)


alls

k,3,,,l0

l,2,3
lesize,2,,,enum
lmesh,all


alls

k,4,,,-l0

l,2,4
lesize,3,,,enum
lmesh,all


k,5,-l0


l,2,5
lesize,4,,,enum
lmesh,all


d,node(-l0,0,0),all
d,node(0,0,-l0),all
d,node(0,0,l0),all


d,3,uy,-1.2
d,4,uy,-3.5922
d,5,uy,-7.0613
d,6,uy,-10.407
d,7,uy,-13.631
d,8,uy,-16.731
d,9,uy,-19.708
d,10,uy,-22.563
d,11,uy,-25.294
d,12,uy,-27.903
d,13,uy,-30.390
d,14,uy,-32.754
d,15,uy,-34.995
d,16,uy,-37.114
d,17,uy,-39.111
d,18,uy,-40.985
d,19,uy,-42.737
d,20,uy,-44.367
d,21,uy,-45.875
d,22,uy,-47.261
d,23,uy,-48.525
d,24,uy,-49.667
d,25,uy,-50.687
d,26,uy,-51.586
d,27,uy,-52.363
d,28,uy,-53.018
d,29,uy,-53.552
d,30,uy,-53.964
d,31,uy,-54.256
d,32,uy,-54.474
d,33,uy,-54.402
d,34,uy,-54.208
d,35,uy,-53.894
d,36,uy,-53.459
d,37,uy,-52.903
d,38,uy,-52.227
d,39,uy,-51.429
d,40,uy,-50.512
d,41,uy,-49.474
d,42,uy,-48.315
d,43,uy,-47.037
d,44,uy,-45.638
d,45,uy,-44.119
d,46,uy,-42.480
d,47,uy,-40.721
d,48,uy,-38.842
d,49,uy,-36.843
d,50,uy,-34.725
d,51,uy,-32.487
d,52,uy,-30.130
d,53,uy,-27.653
d,54,uy,-25.057
d,55,uy,-22.341
d,56,uy,-19.507
d,57,uy,-16.553
d,58,uy,-13.480

d,59,uy,-10.288
d,60,uy,-6.9776
d,61,uy,-3.5482


alls

nsel,s,,,3,61

d,all,ux,0
alls





d,121,uy,-1.2
d,120,uy,-3.5922
d,119,uy,-7.0613
d,118,uy,-10.407
d,117,uy,-13.631
d,116,uy,-16.731
d,115,uy,-19.708
d,114,uy,-22.563
d,113,uy,-25.294
d,112,uy,-27.903
d,111,uy,-30.390
d,110,uy,-32.754
d,109,uy,-34.995
d,108,uy,-37.114
d,107,uy,-39.111
d,106,uy,-40.985
d,105,uy,-42.737
d,104,uy,-44.367
d,103,uy,-45.875
d,102,uy,-47.261
d,101,uy,-48.525
d,100,uy,-49.667
d,99,uy,-50.687
d,98,uy,-51.586
d,97,uy,-52.363
d,96,uy,-53.018
d,95,uy,-53.552
d,94,uy,-53.964
d,93,uy,-54.256
d,92,uy,-54.474
d,91,uy,-54.402
d,90,uy,-54.208
d,89,uy,-53.894
d,88,uy,-53.459
d,87,uy,-52.903
d,86,uy,-52.227
d,85,uy,-51.429
d,84,uy,-50.512
d,83,uy,-49.474
d,82,uy,-48.315
d,81,uy,-47.037
d,80,uy,-45.638
d,79,uy,-44.119
d,78,uy,-42.480
d,77,uy,-40.721
d,76,uy,-38.842
d,75,uy,-36.843
d,74,uy,-34.725
d,73,uy,-32.487
d,72,uy,-30.130
d,71,uy,-27.653
d,70,uy,-25.057
d,69,uy,-22.341
d,68,uy,-19.507
d,67,uy,-16.553
d,66,uy,-13.480

d,65,uy,-10.288
d,64,uy,-6.9776
d,63,uy,-3.5482


alls

nsel,s,,,63,121

d,all,uz,0
alls




d,181,uy,-1.2
d,180,uy,-3.5922
d,179,uy,-7.0613
d,178,uy,-10.407
d,177,uy,-13.631
d,176,uy,-16.731
d,175,uy,-19.708
d,174,uy,-22.563
d,173,uy,-25.294
d,172,uy,-27.903
d,171,uy,-30.390
d,170,uy,-32.754
d,169,uy,-34.995
d,168,uy,-37.114
d,167,uy,-39.111
d,166,uy,-40.985
d,165,uy,-42.737
d,164,uy,-44.367
d,163,uy,-45.875
d,162,uy,-47.261
d,161,uy,-48.525
d,160,uy,-49.667
d,159,uy,-50.687
d,158,uy,-51.586
d,157,uy,-52.363
d,156,uy,-53.018
d,155,uy,-53.552
d,154,uy,-53.964
d,153,uy,-54.256
d,152,uy,-54.474
d,151,uy,-54.402
d,150,uy,-54.208
d,149,uy,-53.894
d,148,uy,-53.459
d,147,uy,-52.903
d,146,uy,-52.227
d,145,uy,-51.429
d,144,uy,-50.512
d,143,uy,-49.474
d,142,uy,-48.315
d,141,uy,-47.037
d,140,uy,-45.638
d,139,uy,-44.119
d,138,uy,-42.480
d,137,uy,-40.721
d,136,uy,-38.842
d,135,uy,-36.843
d,134,uy,-34.725
d,133,uy,-32.487
d,132,uy,-30.130
d,131,uy,-27.653
d,130,uy,-25.057
d,129,uy,-22.341
d,128,uy,-19.507
d,127,uy,-16.553
d,126,uy,-13.480

d,125,uy,-10.288
d,124,uy,-6.9776
d,123,uy,-3.5482


alls

nsel,s,,,123,181

d,all,uz,0
alls







d,241,uy,-1.2
d,240,uy,-3.5922
d,239,uy,-7.0613
d,238,uy,-10.407
d,237,uy,-13.631
d,236,uy,-16.731
d,235,uy,-19.708
d,234,uy,-22.563
d,233,uy,-25.294
d,232,uy,-27.903
d,231,uy,-30.390
d,230,uy,-32.754
d,229,uy,-34.995
d,228,uy,-37.114
d,227,uy,-39.111
d,226,uy,-40.985
d,225,uy,-42.737
d,224,uy,-44.367
d,223,uy,-45.875
d,222,uy,-47.261
d,221,uy,-48.525
d,220,uy,-49.667
d,219,uy,-50.687
d,218,uy,-51.586
d,217,uy,-52.363
d,216,uy,-53.018
d,215,uy,-53.552
d,214,uy,-53.964
d,213,uy,-54.256
d,212,uy,-54.474
d,211,uy,-54.402
d,210,uy,-54.208
d,209,uy,-53.894
d,208,uy,-53.459
d,207,uy,-52.903
d,206,uy,-52.227
d,205,uy,-51.429
d,204,uy,-50.512
d,203,uy,-49.474
d,202,uy,-48.315
d,201,uy,-47.037
d,200,uy,-45.638
d,199,uy,-44.119
d,198,uy,-42.480
d,197,uy,-40.721
d,196,uy,-38.842
d,195,uy,-36.843
d,194,uy,-34.725
d,193,uy,-32.487
d,192,uy,-30.130
d,191,uy,-27.653
d,190,uy,-25.057
d,189,uy,-22.341
d,188,uy,-19.507
d,187,uy,-16.553
d,186,uy,-13.480

d,185,uy,-10.288
d,184,uy,-6.9776
d,183,uy,-3.5482




alls

nsel,s,,,183,241

d,all,ux,0
alls




















!求解、更新、判别收敛条件

/solu
antype,0
nlgeom,on
sstif,on
nsubst,20
outres,all,all
solve

/PREP7
!更新几何形状
UPCOORD,1,0


mp,ex,1,1.2e11

r,1,area,1.5e-4


alls

nsel,s,,,3,61
d,all,uy,0
d,all,ux,0
alls


alls

nsel,s,,,63,121
d,all,uy,0
d,all,uz,0
alls

alls

nsel,s,,,123,181
d,all,uy,0
d,all,uz,0
alls

alls

nsel,s,,,183,241
d,all,uy,0
d,all,ux,0
alls
n,242,1030
n,243,-1030
n,244,,,1030
n,245,,,-1030


e,182,243
e,62,244
e,1,242
e,122,245

nsel,s,,,242,245
d,all,all
alls


/solu
antype,0
nlgeom,on
sstif,on
nsubst,20
outres,all,all
solve

/prep7

ddele,all,all
nsel,s,,,242,245
d,all,all
alls












alls

f,2,fx,1000
alls


allsel

FINISH

/SOLU
ANTYPE,STATIC
NLGEOM,on
SSTIF,ON
NSUBST,1
SOLVE
FINISH
!!!********************************


/POST1
PLNSOL,U,SUM,2

不知道什么问题 总是计算结果不收敛