绝对经典－求助一个载荷覆盖的问题！

mashilong

<FONT style="FONT-SIZE: 12px">
<P>毕业设计中的热分析程序：<br>程序介绍：在平板上（面２）加高斯热源（取１／４平板），平板在煤油中加工，所以想加上对流，所以在面２，４，６上施加了对流载荷（表面上生成了surf152单元，经检查确是生成了），结算后表明高斯热源加上了，但是对流载荷总加不上（运算中出了载荷覆盖的警告）。请您帮我看看程序加载上出了什么问题。最好详细点！谢谢<br>finish<br>/clear<br>/title,transient analysis(gauss heat flux)<br>/filnam,transientgauss<br>/TRIAD,RBOT<br>/UNITS,SI <br>/PREP7 <br>et,1,solid70<br>et,2,surf152 <br>keyopt,2,8,2　　　!要在surf152上施加对流<br>X1=0.2　　　　　　　 !输入平板的能量比例<br>U=25　　　　　　　!电压<br>I=3　　　　　　　　　!电流<br>Ton=100      　　　!脉冲宽度<br>pi=3.1415926                     <br>Q=3*X1*U*I                  <br>lisannum=4           ! 控制载荷布的数量   <br>dividenum=20          !控制网格划分的密度<br>block_x=0.0005<br>block_y=0.0005<br>block_z=-0.00005<br>BLOCK,0,block_x,0,block_y,block_z,0, <br>/VIEW,1,1,1,1<br>save<br>mptemp,1,20,250,500,750,1000,1500,1700,2500 　　　　!定义材料属性<br>mpdata,c,1,1,460,480,530,675,670,660,780,820 <br>mpdata,kxx,1,1,50,47,40,27,30,35,140,142 <br>mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090 <br>mpdata,hf,1,1,100,350,520,1000,1500,3000,3100,3500 <br>SAVE<br>FLST,5,2,4,ORDE,2   　　　　!分网部分<br>FITEM,5,4   <br>FITEM,5,-5  <br>CM,_Y,LINE  <br>LSEL, , , ,P51X <br>CM,_Y1,LINE <br>CMSEL,,_Y   <br>!*  <br>LESIZE,_Y1, , ,dividenum,5, , , ,1 <br>!*  <br>FLST,5,2,4,ORDE,2   <br>FITEM,5,1   <br>FITEM,5,8   <br>CM,_Y,LINE  <br>LSEL, , , ,P51X <br>CM,_Y1,LINE <br>CMSEL,,_Y   <br>!*  <br>LESIZE,_Y1, , ,dividenum,5, , , ,1 <br>!*  <br>FLST,5,2,4,ORDE,2   <br>FITEM,5,1   <br>FITEM,5,8   <br>CM,_Y1,LINE <br>LSEL, , , ,P51X <br>*GET,_z1,LINE,,COUNT<br>*SET,_z2,0  <br>*DO,_z5,1,_z1   <br>*SET,_z2,LSNEXT(_z2)<br>*GET,_z3,LINE,_z2,ATTR,NDNX <br>*GET,_z4,LINE,_z2,ATTR,SPNX <br>*get,_z6,line,_z2,attr,kynd <br>*IF,_z3,GT,0,THEN   <br>*IF,_z4,NE,0,THEN   <br>LESIZE,_z2,,,_z3,1/_z4,,,,_z6   <br>*ENDIF  <br>*ENDIF  <br>*ENDDO  <br>CMSEL,S,_Y1 <br>CMDELE,_Y1  <br>!*  <br>FLST,5,3,4,ORDE,3   <br>FITEM,5,9   <br>FITEM,5,-10 <br>FITEM,5,12  <br>CM,_Y,LINE  <br>LSEL, , , ,P51X <br>CM,_Y1,LINE <br>CMSEL,,_Y   <br>!*  <br>LESIZE,_Y1, , ,dividenum/5,5, , , ,1  <br>!*  <br>FLST,5,1,4,ORDE,1   <br>FITEM,5,10  <br>CM,_Y1,LINE <br>LSEL, , , ,P51X <br>*GET,_z1,LINE,,COUNT<br>*SET,_z2,0  <br>*DO,_z5,1,_z1   <br>*SET,_z2,LSNEXT(_z2)<br>*GET,_z3,LINE,_z2,ATTR,NDNX <br>*GET,_z4,LINE,_z2,ATTR,SPNX <br>*get,_z6,line,_z2,attr,kynd <br>*IF,_z3,GT,0,THEN   <br>*IF,_z4,NE,0,THEN   <br>LESIZE,_z2,,,_z3,1/_z4,,,,_z6   <br>*ENDIF  <br>*ENDIF  <br>*ENDDO  <br>CMSEL,S,_Y1 <br>CMDELE,_Y1  <br>!*  <br>CM,_Y,VOLU  <br>VSEL, , , ,       1 <br>CM,_Y1,VOLU <br>CHKMSH,'VOLU'   <br>CMSEL,S,_Y  </P>
<P>MSHKEY,1<br>VATT,1,1,1<br>VMESH,1<br>/PNUM,DEFA<br>EPLOT<br>FINISH　　　　　　　　　　　　　!分网结束</P>
<P>*dim,R,array,lisannum,1,1　　　!高斯热流密度的计算<br>tempr=0.11*I*I+0.73*I+8.52<br>*IF,tempr,GE,Ton,then               <br>    maxR=2.85*I**0.53*Ton**0.38           <br>*ELSE <br>    maxR=2.85*I**0.53*tempr**0.38                   <br>*ENDIF<br>maxenum=dividenum+1<br>*dim,flux2,table,maxEnum,maxEnum,lisannum,x,y,time   <br>*DO,k,1,lisannum,1<br>  t=k*Ton/lisannum       <br>    R(k)=2.85*I**0.53*t**0.38  <br>      *IF,R(k),GE,maxR,then                    <br>       R(k)=maxR<br>      *ENDIF <br>  *do,j,1,maxEnum,1 <br>   flux2(0,j,k)=(j-1)*block_x/dividenum <br>  *enddo   <br>  *do,i,1,maxEnum,1 <br>   flux2(i,0,k)=(i-1)*block_x/dividenum<br>  *enddo <br>  *do,i,1,maxEnum,1 <br>    *do,j,1,maxEnum,1 <br>     xcenter=0 <br>     ycenter=0 <br>     distance=((i-1)*block_x/dividenum-xcenter)**2+((j-1)*block_x/dividenum-ycenter)**2 <br>    *if,distance,le,(R(k)*1E-6)**2,then <br>        flux2(i,j,k)=(Q/(pi*R(k)*R(k)*1E-12)*exp(-3*distance/R(k)*R(k)*1E-12)) !???????? <br>    *else <br>        flux2(i,j,k)=0 <br>    *endif   <br>    *enddo <br>  *enddo <br>       flux2(0,0,k)=(k-1)*Ton/lisannum <br>*enddo <br>/SOLU    <br>ANTYPE,trans<br>TRNOPT,FULL <br>NROPT,AUTO, ,   <br>EQSLV,   　　　　　　<br>SOLCONTROL,ON <br>AUTOTS,on<br>KBC,1                            <br>OUTRES,ALL,ALL,                           <br>TIME,1e-6 <br>DELTIM,1E-6,1E-6,1E-6    <br>TUNIF,25       <br>kbc,1 <br>solve                !先施加25度的开始温度计算，只是第一步起作用<br>type,2 <br>asel,s,,,2,6,2 <br>NSLA,S,1<br>ESLN,S,0<br>ESURF 　　　　　　　!在２，４，６面上生成表面效应单元<br>esel,s,type,,2 <br>sf,all,conv,-1,25 　!施加对流（hf已经在前面材料属性部分定义了）<br>allsel,all <br>nsel,all <br>eplot <br>totaltime=Ton*1E-6<br>inctim=Ton*1E-6/lisannum<br>autots,1<br>*do,t,2e-6,totaltime,inctim <br>time,t  <br>deltim,inctim,inctim,inctim<br>kbc,0 <br>outres,all,all, <br>tsres,erase <br>asel,s,loc,z,0,0<br>sfa,all,1,hflux,%flux2%  !施加高斯热源<br>solve<br>OUTRES,ALL,ALL,<br>/PSF,HFLUX,,2<br>/REPLOT<br>*enddo 　　　　　　<br>finish </P>
<P>毕业设计急用，请大家一定要帮我看看</P><br></FONT>

[此贴子已经被作者于2006-3-14 13:31:02编辑过]