snowhite 发表于 2010-8-16 17:15:46

请教：ansys 3D热的流体场仿真问题

模型是一个变压器放在油里冷却，损耗作为激励，最后想仿出温度分布，将油设成流体。

简要程序如下：
/prep7
!单元转换
      et,7,FLUID142   !线圈
      et,2,FLUID142   ! 磁芯
      et,3,FLUID142   ! FR4
      et,4,FLUID142   ! gap
      et,6,FLUID142   ! linker
      et,1,FLUID142   ! oil对油进行流体分析

!定义热参数
       FLData7,PROT,dens,Table
      FLData7,PROT,visc,Table!流体粘性类型
      FLData7,PROT,cond,Table
      FLData7,PROT,spht,Table
      MPTEMP,1,290,300,310,320,330,340,350,360,370,380,390,400,410,420,430
MPDATA,DENS,1,,890,884.1,877.9,871.8,865.8,859.9,853.9,847.8,841.8,836.0,830.6,825.1,818.9,812.1,806.5
MPDATA,VISC,1,,1120E-6,550E-6,288E-6,161E-6,96.6E-6,61.7E-6,41.7E-6,29.7E-6,22.0E-6,16.9E-6,13.3E-6,10.6E-6,8.52E-6,6.94E-6,5.83E-6
MPDATA,KXX,1,,0.145,0.145,0.145,0.143,0.141,0.139,0.138,0.138,0.137,0.136,0.135,0.134,0.133,0.133,0.132
MPDATA,KYY,1,,0.145,0.145,0.145,0.143,0.141,0.139,0.138,0.138,0.137,0.136,0.135,0.134,0.133,0.133,0.132
MPDATA,KZZ,1,,0.145,0.145,0.145,0.143,0.141,0.139,0.138,0.138,0.137,0.136,0.135,0.134,0.133,0.133,0.132
MPDATA,C,1,,1868,1909,1951,1993,2035,2076,2118,2161,2206,2250,2294,2337,2381,2427,2471
      FLData13,vary,dens,1
      FLData13,vary,VISC,1
      FLData13,vary,COND,1
      FLData13,vary,SPHT,1
      *do,i,1,Nc,1
      MP,KXX,6+i,400      !导热系数 W/m.C
MP,KYY,6+i,400      !导热系数 W/m.C
MP,KZZ,6+i,400      !导热系数 W/m.C
      MP,DENS,6+i,8933   !密度kg/m^3
      MP,C,6+i,385      !比热J/g.C
      *enddo
   MP,KXX,6,400      !导热系数 W/m.C
MP,KYY,6,400      !导热系数 W/m.C
MP,KZZ,6,400      !导热系数 W/m.C
      MP,DENS,6,8933   !密度kg/m^3
      MP,C,6,385      !比热J/g.C
      MP,KXX,2,4
      MP,KYY,2,4
      MP,KZZ,2,4
      MP,DENS,2,4600
      MP,C,2,750
      MP,KXX,3,0.294
MP,KYY,3,0.294
MP,KZZ,3,0.294
      MP,DENS,3,1900
      MP,C,3,1150
      MP,KXX,4,0.4
      MP,KYY,4,0.4
      MP,KZZ,4,0.4
      MP,DENS,4,930
      MP,C,4,2300
            
/prep7
！施加激励源
alls
*do,k,1,Ne
esel,s,elem,,k
bfe,all,hgen,,Pv(k)
*enddo

/solu
！设置边界条件   
*SET,T0,273+25
alls
nsel,s,loc,x,0.36   
d,all,vx,0
alls
nsel,s,loc,y,0.27   
d,all,vy,0
alls
nsel,s,loc,z,-0.05
nsel,s,loc,z,0.05   
d,all,vz,0
alls
asel,s,,,461,464
nsla,s,1
d,all,temp,T0
d,all,pres,0
!   Radiation;            
      !线圈
      alls
      vsel,s,type,,7
      nslv,s,1
      SF,All,Rad,0.78,,T0
      !linker
      alls
      vsel,s,type,,6
      nslv,s,1
      SF,All,Rad,0.78,,T0
      !磁芯
      alls
      vsel,s,mat,,2
      nslv,s,1
      SF,All,Rad,0.8,,T0
      !FR4
      alls
      vsel,s,mat,,3
      nslv,s,1
      SF,All,Rad,0.91,,T0
      !gap
      alls
      vsel,s,mat,,4
      nslv,s,1
      SF,All,Rad,0.91,,T0
      !油
      alls
      vsel,s,mat,,1
      nslv,s,1
      SF,All,Rad,0.85,,T0
Acel,0,9.81

/solu
！稳态场   
fldata1,solu,flow,t
fldata1,solu,temp,t
fldata14,temp,nomi,T0   
fldata25,relx,temp,0.8
fldata2,iter,exec,20
fldata18,meth,temp,5
fldata18,meth,press,1   
fldata,solu,tran,0
fldata,conv,pres,1e-10
fldata,conv,temp,1e-3   
alls
solve   
FINISH
      
/POST1
SET,LAST
PLNSOL,TEMP    !绘制温度场分布等值线图
FINISH

/solu
！瞬态场   
fldata1,solu,flow,t
fldata1,solu,temp,t
fldata14,temp,nomi,T0   
fldata18,meth,temp,1
fldata,solu,tran,1
!fldata4,time,step,1
fldata4,time,numb,10
fldata4,time,glob,40
!fldata4,time,appe,1
fldata4,time,tend,180    !end time
alls
solve   
FINISH
/post1
alls
set,last
plnsol,temp

稳态场分析后，所有的网格温度都是初始温度；
瞬态场10s可以得到与实际测量符合的温度分布；
瞬态场180s得到的温度误差就会很大，温度随时间上升不是一个线性的过程；
想请教下，是我的边界条件设置有问题吗？

kuangulugu 发表于 2010-8-16 18:06:44

收藏，好好看看

snowhite 发表于 2010-8-18 17:14:58

困惑，谁能帮个忙吗？仿出来的温度大概只有实测的一半不到。。。

yuguoxingqin 发表于 2010-8-20 15:09:26

问题解决没？我们做的问题比较类似，可以互相讨论！

Albert1984 发表于 2010-9-8 09:59:42

谁有比较好的关于流场分析的书啊
页: [1]
查看完整版本: 请教：ansys 3D热的流体场仿真问题