下天雄 发表于 2008-9-18 10:47:15

ansys循环遇到的比较诡异的问题

我循环设置了10步
程序计算10后,不停止
继续进行计算
还有,原先赋值的一些常量也变成诸如
Unknown parameter name= *****.A value of 7.888609052E-31 will be
used.


不知道这是什么原因造成的?

我仔细检查过每个循环
*do....
*enddo
都是闭合的

luyulin 发表于 2008-9-18 10:52:52

参数名不对

下天雄 发表于 2008-9-18 10:56:53

您好

这个情况出现在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
参数名不对

下天雄 发表于 2008-9-18 12:27:22

以下命令是稳态计算后,进入瞬态计算的命令流
也就是循环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

下天雄 发表于 2008-9-18 12:28:13

大家有时间的话
帮忙看看为什么循环10次后
程序不停止的原因

谢谢大家了

yuanhang621 发表于 2009-5-7 19:22:11

多工矿能循环么
页: [1]
查看完整版本: ansys循环遇到的比较诡异的问题