求助-热流耦合结果不对
本帖最后由 transient 于 2014-4-28 14:28 编辑各位前辈， 大家好!我在做一个耦合场分析，先是热场（线圈通电产生热），接着通过载荷传递将温度传递到流场（流场流动散去一部分热量，导致降温），接着流场的分析结果再将温度传递到下一个载荷步的热场，如此循环下去。简易的命令流如下： 1建模略2 设置热场分析选项 toffst,273!定义初始温度 !nsel,all!IC,ALL,TEMP,tf !TUNIF,20ESEL,ALL !设定耦合电压 !入口ESEL,S,MAT,,6NSLE,S,ALLNSEL,R,LOC,Y,0cp,1,VOLT,ALLNSEL,ALL !出口ESEL,S,MAT,,6NSLE,S,ALLNSEL,R,LOC,Y,180cp,2,VOLT,ALLD,4881,VOLT,0NSEL,ALL !耦合温度!外层ESEL,S,MAT,,6NSLE,S,ALLnsel,R,LOC,X,R9CP,3,TEMP,all!CPSGEN,3,100,1NSEL,ALL !中间层ESEL,S,MAT,,6NSLE,S,ALLnsel,R,LOC,X,(R9+R3)/2CP,4,TEMP,allNSEL,ALL !里层ESEL,S,MAT,,6NSLE,S,ALLnsel,R,LOC,X,R3CP,5,TEMP,allNSEL,ALL ESEL,ALL !加载 !F,1,AMPS,10*3.415F,5451,AMPS,120*3.415 !对流条件 nsel,s,LOC,X,R8!nsel,a,LOC,X,R7nsel,A,LOC,Z,0nsel,A,LOC,Z,h1 SF,ALL,CONV,10,20 !求解 ANTYPE,transKBC,1trnopt,fullOUTRES,,ALLALLSEL,ALL!AUTOTS,ONDELTIM,0.2,0.2,0.2,on!TIME,1 !一个时间步的时间 physics,write,thermal!physics,clear!savefinish 3 设置流场分析选项 /PREP7 lsclear,allphysics,clear !流场分析et,1,FLUID142 !参量PREP=10 !材料 !定义热传导率MP,KXX,2,1MP,KXX,3,81MP,KXX,4,50.2MP,KXX,5,1MP,KXX,6,400 !COILSMP,KXX,7,237MP,KXX,1,0.54 !非流体材料属性 MP,DENS,1,1e3MP,DENS,2,2e3MP,DENS,3,7.8e3MP,DENS,4,7.8e3MP,DENS,5,1.35e3MP,DENS,6,8.9e3MP,DENS,7,2.7e3 !比热 MP,C,1,4.2e3MP,C,2,600MP,C,3,444MP,C,4,450MP,C,5,1e3MP,C,6,385MP,C,7,880 !流体参数RHO=1e3MU=1.005e-03K=0.54CP=4.2e3 !边界条件!水道两侧nsel,s,LOC,X,R6nsel,a,LOC,X,R7D,ALL,VY,0D,ALL,VX,0NSEL,ALL !上壁面温度!nsel,s,LOC,Y,R1!d,all,temp,ts !NSEL,ALL !入口压强，温度 NSEL,S,LOC,X,R6,R7NSEL,R,LOC,z,h1 D,ALL,VX,0D,ALL,Vy,0!D,all,vz,0.5D,ALL,PRES,10d,all,temp,tf NSEL,ALL !对称面边界!NSEL,S,LOC,y,0!NSEL,r,LOC,y,180!D,ALL,VX,0!D,ALL,Vy,0!D,ALL,Vz,0!NSEL,ALL !出口压强 NSEL,s,LOC,X,R6,R7NSEL,r,LOC,z,0 D,ALL,VX,0D,ALL,Vy,0D,ALL,PRES,0 NSEL,ALL !对流边界条件 nsel,s,LOC,X,R8!nsel,a,LOC,X,R7nsel,A,LOC,Z,0nsel,A,LOC,Z,h1 SF,ALL,CONV,10,20 allselFINISH /SOLU toffst,273!LDREAD,TEMP,,,,2,,RST!LDREAD,HGEN,,,,,,RST !流场求解设置FLDATA,ITER,EXEC,300 !定义迭代次数FLDATA,ITER,CHEC,10 !CHECKPOINT FREQUENCY!FLDATA,TEMP,NOMI,tf !REFERENCE TEMPERATUREFLDATA,PROT,DENS,RHOFLDATA,PROT,VISC,MUFLDATA,PROT,COND,KFLDATA,PROT,SPHT,CP!FLDATA,OUTP,TAUW,T FLDATA4,TIME,STEP,0.2 !每一步(是子步吗?)的时间间隔FLDATA1,SOLU,TRAN,T!FLDATA4,TIME,ISTEP,10FLDATA4,TIME,NUMB,10 !子步数FLDATA4,TIME,glob,50FLDATA4A,STEP,APPE,1 !deltim,10,1,10!ddele,all,temp nsel,all!SOLVE !温度场求解!FLDATA,SOLU,enrg,TFLDATA,ITER,EXEC,50FLDATA,SOLU,TEMP,TFLDATA,SOLU,FLOW,tFLDATA,RELX,TEMP,0.8nsel,all!SOLVEFINISH physics,write,fluid 4 设置耦合分析选项，如下!开始求解 *do,i,1,2time=time+2 physics,read,thermal /solu *if,i,eq,1,then nsel,all IC,ALL,TEMP,tf ! initial temperature*else !ldread,temp,last,,,,,rst !ldread,temp,,,,2,,rst !read thermal analysis temperatures,是否准确？ antype,trans,rest,i,1 !ldread,temp,i,,,2,,rfl !此命令没用 ldread,temp,last,,,2,,rfl !上一步的流场分析温度作为初状态 *endiftime,time solve ! solve harmonicanalysisfinish !读入流场physics,read,fluid/assign,esav,fluidd,esav ! redirect files for use in thermalrestart/assign,emat,fluidd,emat /soluparsav,scalar,parameter,sav !save parameters before multiframe restart*if,i,gt,1,then !FLDATA32,REST,LSTP,0 ! thermal restart !FLDATA32,REST,time,i*2 FLDATA32,REST,LSTP,0*endifparres,new,parameter,sav !FLDATA4,TIME,TEND,time !*if,i,eq,1,then! ldread,temp,,,1,2,,rst !apply coupled joule heating load from !*else! ldread,temp,last,,,2,,rst ! apply coupled joule heating load from ! *endif ldread,temp,i,,,2,,rst ! apply coupled joule heating loadfrom 有问题 esel,allsolvefinish/assign,esav ! reassign files to default/assign,emat*enddo ! end of solutionlooping finishsave 在我试算了2个载荷步后发现1 第一个载荷步是正确的，流场分析前热场分析的温度结果作为初始状态了2 第二个载荷步开始，热场分析的初始温度不是第一步的流场分析结果而是第一步的热场分析结果的温度，这就不对了。请教大家应该怎么设置？ 多谢！！！ 多耦合场还是wb解决比较好
小七工作室 发表于 2014-5-4 16:28
多耦合场还是wb解决比较好
WB能支持命令流吗？
我这个模型必须映射划分网格，所以我是用命令流直接画单元的。
页:
[1]