ansys循环遇到的比较诡异的问题
我循环设置了10步程序计算10后,不停止
继续进行计算
还有,原先赋值的一些常量也变成诸如
Unknown parameter name= *****.A value of 7.888609052E-31 will be
used.
不知道这是什么原因造成的?
我仔细检查过每个循环
*do....
*enddo
都是闭合的 参数名不对 您好
这个情况出现在10次循环结束之后:Unknown parameter name= *****.A value of 7.888609052E-31 will be used.
在10次循环内,各个参数都没有出现异常
所以,我不是很清楚,为什么会继续循环下去?
原帖由 luyulin 于 2008-9-18 10:52 发表 http://forum.simwe.com/images/common/back.gif
参数名不对 以下命令是稳态计算后,进入瞬态计算的命令流
也就是循环10次的命令流
/PREP7
allsel,all
LSCLEAR,fe
LSCLEAR,solid
NEQIT,25,
LNSRCH,1
timestep=0.1 !energy 0.01
QsupAL=0
Usupband=0
current=current0
QsupAL=0
Qsupband=0
Qcoil=0
qmaxz=hdisturbz+cableh/2
qmaxy=hdisturbst+divdegree/2
sumofqc=0
sumofQsupAL=0
sumofQsupband=0
calcno=10
*do,i,1,calcno
/solu
ldread,temp,last,,,,,rth
TIME,timestep*i+0.1
Iteam(i)=current
Rteam(i)=sumofqc/(current**2)
quenchmaxz(i)=qmaxz
quenchmaxy(i)=qmaxy
sumofqc=0
sumofQsupAL=0
sumofQsupband=0
/POST1
set,last
*DO,i,1,vno-tnvno
vsel,s,,,2*(vno-tnvno)+i
NSLV,S
*get,n,node,,num,max
*get,T,node,n,temp
*if,T,gt,(Tdiv+Tcri)/2,then
/solu
*GET,qareaz,NODE,n,LOC,z,
*GET,qareay,NODE,n,LOC,y,
quenchareaz(i)=qareaz
quenchareay(i)=qareay
*else
quenchareaz(i)=0
quenchareay(i)=0
*endif
/post1
!电阻率
*if,T,lt,Tdiv,then
Rc=0
*elseif,T,lt,Tcri
R1=2.86849E-10-1.09501E-13*T-5.92221E-13*T**2
R2=+2.63634E-14*T**3-2.67527E-16*T**4+1.25582E-18*T**5
R3=-2.86424E-21*T**6+2.56929E-24*T**7
R=R1+R2+R3
Rc=R*(T-Tdiv)/(Tcri-Tdiv)
*elseif,T,lt,294
R1=2.86849E-10-1.09501E-13*T-5.92221E-13*T**2
R2=+2.63634E-14*T**3-2.67527E-16*T**4+1.25582E-18*T**5
R3=-2.86424E-21*T**6+2.56929E-24*T**7
Rc=R1+R2+R3
*else
Rc=1.68e-8
*endif
qcteam(i)=current**2*Rc/(Sc**2)
sumofqc=sumofqc+0.92*qcteam(i)*cableh*cablew*2*pi*(inrad+supalh+groinsuh)/tnvno
*enddo
Qcoil=Qcoil+sumofqc*timestep
*DO,i,1,vno-tnvno
vsel,s,,,i
NSLV,S
*get,n,node,,num,max
*get,T,node,n,temp
*if,T,lt,294,then
R1=1.9367E-9-2.9274E-11*T+5.58251E-13*T**2
R2=+2.33302E-14*T**3-2.99186E-16*T**4+1.4937E-18*T**5
R3=-3.45234E-21*T**6+3.07247E-24*T**7
Rsu=R1+R2+R3
*else
Rsu=3.4e-8
*endif
qsuteam_al(i)=(QsupAL/(2*pi*inrad))**2/Rsu
sumofQsupAL=sumofQsupAL+qsuteam_al(i)*Ss_al*2*pi*inrad/tnvno
*enddo
QsupAL=QsupAL+sumofQsupAL*timestep
*DO,i,1,vno-tnvno
vsel,s,,,4*(vno-tnvno)+i
NSLV,S
*get,n,node,,num,max
*get,T,node,n,temp
*if,T,lt,294,then
R1=5.90636E-10-5.08481E-11*T+1.98842E-12*T**2
R2=-1.42089E-14*T**3+4.83747E-17*T**4-7.46509E-20*T**5
R3=+3.7428E-23*T**6
Rsu=R1+R2+R3
*else
Rsu=1.8375e-8
*endif
qsuteam_band(i)=(Usupband/(2*pi*(inrad+supalh+groinsuh+cablew+groinsuh)))**2/Rsu
sumofQsupband=sumofQsupband+qsuteam_band(i)*Ss_band*2*pi*(inrad+supalh+groinsuh+cablew+groinsuh)/tnvno
*enddo
Qsupband=Qsupband+sumofQsupband*timestep
currentn=((Qsum-Qcoil-QsupAL-Qsupband)*2/L)**0.5
QsupAL=M1*(current-currentn)/timestep
Usupband=M2*(current-currentn)/timestep
current=currentn
*if,current,lt,0,then
*exit
*endif
/PREP7
allsel,all
LSCLEAR,fe
LSCLEAR,solid
/solu
ANTYPE,,REST
*DO,i,1,vno-tnvno
vsel,s,,,2*(vno-tnvno)+i
bfv,all,hgen,qcteam(i)
vsel,s,,,i
bfv,all,hgen,qsuteam_al(i)
vsel,s,,,4*(vno-tnvno)+i
bfv,all,hgen,qsuteam_band(i)
*enddo
allsel,all
solve
*enddo
/POST1
SET,LAST
PLNSOL, TEMP,, 0 大家有时间的话
帮忙看看为什么循环10次后
程序不停止的原因
谢谢大家了 多工矿能循环么
页:
[1]