/units,si
/prep7
et,1,plane42
keyopt,1,3,2
mp,ex,1,220e9
mp,nuxy,1,0.3
TB,MISO,1,1,25, 
TBTEMP,0
TBPT,,0.000814,175.88e6,
TBPT,,0.001531,324.98e6,
TBPT,,0.002819,591.26e6,
TBPT,,0.003846,801.457e6,
TBPT,,0.004239,873.22e6,
TBPT,,0.004594,924.399e6,
TBPT,,0.004955,956.88e6,
TBPT,,0.00531,976.99e6,
TBPT,,0.005699,990.2138e6,
TBPT,,0.006261,1000e6,
TBPT,,0.006957,1008.34e6,
TBPT,,0.007804,1015.5e6,
TBPT,,0.008779,1021.185e6,
TBPT,,0.009827,1026.47e6,
TBPT,,0.010929,1031.289e6,
TBPT,,0.012066,1035.444e6,
TBPT,,0.0132,1039.693e6,
TBPT,,0.014354,1044.886e6,
TBPT,,0.015494,1049.796e6,
TBPT,,0.016631,1053.951e6,
TBPT,,0.017785,1058.106e6,
TBPT,,0.019039,1062.827e6,
TBPT,,0.020276,1067.643e6,

csys,1
k,1,0.00512
k,2,0.0064
kgen,36,1,2,1,,10
*do,j,1,69,2
 a,j,j+1,j+3,j+2
*enddo
a,71,72,2,1
save
lesize,all,,,4
amesh,all
nummrg,all
numcmp,all
save
cp,1,uy,166,170,171,172,173
finish

/solu
nsel,s,loc,y,270
d,all,all
nplot
allsel
nplot
nsel,s,loc,y,90
d,all,ux
nplot
allsel

outres,all,all
time,1
NSUBST,100
nlgeom,on 

d,166,uy,-0.00128
solve
*get,def1,node,166,U,Y
nsel,s,loc,y,90
ddele,all,uy
allsel
f,166,fy,0
solve
*get,def2,node,166,U,Y
finish
/post1        !
set,list        !ʾмĽ
set,1,last,1     !ȡһزӲĽ
pldisp,1       !ͬʱʾǰκṹ״
plnsol,s,eqv,1,1  !ʾڵĵЧӦ
set,2,last,1     !ȡڶزӲĽ
pldisp,1       !ͬʱʾǰκṹ״
plnsol,s,eqv,1,1  !ʾڵĵЧӦ





