命令流修改

yangdong28

以下是郝文化主编的《ANSYS土木工程应用实例》里的“回弹”的命令流，运行起来老出错，哪个高手能修改啊？
/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:D
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  !显示节点的等效应力

5dehaoma

俺感觉你那带表情的那行有问题，你看看，是否少打了逗号，
cp，1，uy，，166，170，171，172，173

5dehaoma

不好意思，刚才看错了，上面那没有问题，你的加载时的d,166,uy,-0.00128是不是应该在卡氏坐标下加载阿？你目前是在柱坐标下的，疑问

sysh320

我也有这本书，当时做了没有问题的，我给你对照一下啊

ant008

是不是输入的问题，好好检查一下。