瞬态分析：请高手指点该怎么做？

ydcow

是一个瞬态分析，不知道怎么做？请高手指点！
time,0.2
nsubst,10
t1=0.2
z01=0.01*SIN(10*t1)
z02=0.1*COS(10*t1)
z03=-SIN(10*t1)
k1=k+a0*m2+a1*c

z1=z4

z11=z1+z01
z12=z2+z02

p1=c*z12+k*z11-m2*g+m2*(a0*z21+a2*z22+a3*z23)+c*(a1*z21+a4*z22+a5*z23)
z51=p1/k1
z53=a0*(z51-z21)-a2*z22-a3*z23
z52=z22+a6*z23+a7*z53
F=m1*(z03+z3)+m1*g+m2*z53+m2*g

fdele,all,all

f,2,fy,-F
solve
z4=uy(2)
z2=uy(2)
z3=uy(2)
z21=z51
z22=z52
z23=z53
*status,z4
parsave,all

!下面我想运行一个循环，知道判断条件，当abs（（z4-z1）/z4）<0.01时循环终止，不知道该怎么做？（循环是重启动分析，不知道这样的循环中，我后一次循环能利用前一次循环的参数不？）
*do,i,1,50
/solu
antype,,rest,1,10
parresume
time,0.2
nsubst,5
t1=0.2
z01=0.01*SIN(10*t1)
z02=0.1*COS(10*t1)
z03=-SIN(10*t1)
k1=k+a0*m2+a1*c

z1=z4
*status,z1
z11=z1+z01
z12=z2+z02

p1=c*z12+k*z11-m2*g+m2*(a0*z21+a2*z22+a3*z23)+c*(a1*z21+a4*z22+a5*z23)
z51=p1/k1
z53=a0*(z51-z21)-a2*z22-a3*z23
z52=z22+a6*z23+a7*z53
F=m1*(z03+z3)+m1*g+m2*z53+m2*g

fdele,all,all

f,2,fy,-F
solve
z4=uy(2)
z2=uy(2)
z3=uy(2)
z21=z51
z22=z52
z23=z53
parsave
fini
*if,abs((z4-z1)/z4),lt,0.01,then
break
*endif
*enddo

[ 本帖最后由 16443 于 2009-9-17 09:20 编辑 ]