回弹.txt

是用apdl语言编的

/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  !显示节点的等效应力