最近在做一个感应加热钎焊的温度场模拟，但是结果出来的温度都是我设置的参考温度，不知道是怎么回事！
模型如下：
红色部分是感应线圈，感应线圈中加载交流电，线圈中通有水；中间绿色的是我要进行加热的试样（从上至下一次上：银触点、钎料层（很薄一层）、铜基体）；这个模型外面会加上一层空气！
我需要的计算的是模型中绿色试样上的温度场分布！
做的时候我是仿照HELP里面那个多物理场耦合的例子做的，在计算的过程中不会出现错误或者警告，但是结果出来不对，加热的试样温度就等于设置的参考温度！
命令流如下：
/prep7
!coil model paramater
r=0.012
r1=0.004
r2=0.0025
ang1=90
ang2=80
!air model paramater
aa=6*r
ba=aa
da=10*r1
!cotact model paramater
ac=0.01
bc=0.01
dc=0.002
!filler model paramater
af=0.01
bf=0.01
df=0.00015
!base model paramater
ab=0.015
bb=0.01
db=0.003
ab1=0.0024
bb1=0.01
db1=0.01
ab2=0.01
bb2=0.01
db2=0.003
!time
ftime=0.2
tinc=.1
time=0
delt=0.1
freq=40000
pi=4*atan(1)
cond=1/(1.75e-8)
muzero=4e-7*pi
mur=1
skind=sqrt(1/(pi*freq*cond*muzero*mur))
emunit,mks
!modeling
!air
k,1,-aa/2,-da/2,-ba/2
k,2,-aa/2,-da/2,ba/2
k,3,aa/2,-da/2,ba/2
k,4,aa/2,-da/2,-ba/2
a,1,2,3,4
voffst,1,da
!water
cyl4,r,0,0,0,r2,360
k,17,0,0.01,0
k,18,0,0,0
vrotat,7,,,,,,17,18,ang1
vrotat,7,,,,,,17,18,-ang1
vrotat,12,,,,,,17,18,ang2
vrotat,17,,,,,,17,18,-ang2
!coil
cyl4,r,0,r2,0,r1,360
vrotat,28,,,,,,17,18,ang1
vrotat,28,,,,,,17,18,-ang1
vrotat,37,,,,,,17,18,ang2
vrotat,46,,,,,,17,18,-ang2
!contact
k,100,ac/2,-r1,-bc/2
k,101,ac/2,-r1,bc/2
k,102,-ac/2,-r1,bc/2
k,103,-ac/2,-r1,-bc/2
a,100,101,102,103
voffst,65,dc
!filler
k,108,af/2,-(r1+dc),-bf/2
k,109,af/2,-(r1+dc),bf/2
k,110,-af/2,-(r1+dc),bf/2
k,111,-af/2,-(r1+dc),-bf/2
a,108,109,110,111
voffst,71,df
!base
k,116,af/2,-(r1+dc+df),-bb/2
k,117,af/2,-(r1+dc+df),bb/2
k,118,af/2-ab,-(r1+dc+df),bb/2
k,119,af/2-ab,-(r1+dc+df),-bb/2
a,116,117,118,119
voffst,77,db
k,124,af/2-ab+ab1,-(r1+dc+db+df),-bb1/2
k,125,af/2-ab+ab1,-(r1+dc+db+df),bb1/2
k,126,af/2-ab,-(r1+dc+db+df),bb1/2
k,127,af/2-ab,-(r1+dc+db+df),-bb1/2
a,124,125,126,127
voffst,83,db1
k,132,af/2-ab,-(r1+dc+df+db+db1),-bb2/2
k,133,af/2-ab,-(r1+dc+df+db+db1),bb2/2
k,134,af/2-ab-ab2,-(r1+dc+df+db+db1),bb2/2
k,135,af/2-ab-ab2,-(r1+dc+df+db+db1),-bb2/2
a,132,133,134,135
voffst,89,-db2
allsel,all
vovlap,all
numcmp,all
! Electromagnetic
et,1,97
et,2,97
et,3,97,1
!air
mp,murx,1,1
!base
mp,murx,2,1
!coil
mp,murx,3,1
!water
mp,murx,4,0.001
!filler
mp,murx,5,1
!contact
mp,murx,6,1
mp,perx,1,1.000585
mp,perx,4,81
mptemp
mptemp,1,0,100,200,300,400,500
mptemp,7,600,700,800,900
mpdata,rsvx,2,1,.174e-7,.229e-7,.298e-7,.366e-7,.435e-7,.504e-7
mpdata,rsvx,2,7,.573e-7,.626e-7,.71e-7,.779e-7
mpdata,rsvx,3,1,.174e-7,.229e-7,.298e-7,.366e-7,.435e-7,.504e-7
mpdata,rsvx,3,7,.573e-7,.626e-7,.71e-7,.779e-7
mpdata,rsvx,5,1,.164e-7,.219e-7,.288e-7,.356e-7,.425e-7,.494e-7
mpdata,rsvx,5,7,.563e-7,.616e-7,.70e-7,.769e-7
mpdata,rsvx,6,1,.1586e-7,.208e-7,.276e-7,.338e-7,.44e-7,.47e-7
mpdata,rsvx,6,7,.536e-7,.585e-7,.664e-7,.735e-7
!base
vsel,s,volu,,7,9,1
vatt,2,,1
vplot,all
esize,0.001,0,
vsweep,all
!filler
vsel,s,volu,,6
vatt,5,,1
vplot,all
esize,0.001,0,
vsweep,all
!contact
vsel,s,volu,,1
vatt,6,,1
vplot,all
lsel,s,,,85,88,1
lesize,all,,,10
lsel,s,,,93,96,1
lesize,all,,,5,5
vsweep,all
!coil
vsel,s,volu,,10,13
vatt,3,,3
vplot,all
esize,0.00075,0,
vsweep,all
!water
vsel,s,volu,,2,5,1
vatt,4,,2
vplot,all
esize,0.00125,0,
vsweep,all
!air
vsel,s,volu,,14
vatt,1,,2
vplot,all
esize,0.01,0
mshape,1,3d
mshkey,0
vmesh,all
allsel,all
asel,s,ext
nsla,s,1
d,all,az,0,,,,ay,ax
nsel,all
allsel,all
asel,s,area,,76
nsla,s,1
d,all,volt,0
asel,s,area,,77
nsla,s,1
cp,1,volt,all
*get,n1,node,,num,min
f,n1,amps,10000
nsel,all
/solu
antype,harm
harfrq,40000
physics,write,magenti
finish
physics,clear
/prep7
et,1,70
et,2,0
et,3,0
mp,dens,2,8900
mp,dens,3,8900
mp,dens,4,1000
mp,dens,6,10800
mptemp
mptemp,1,20,100,200,300,400,600
mptemp,7,800
mpdata,kxx,2,1,396,393,389,384,379,366
mpdata,kxx,2,7,352
mpdata,kxx,3,1,396,393,389,384,379,366
mpdata,kxx,3,7,352
mpdata,kxx,4,1,0.47,0.48,0.49,0.5,0.51,0.52
mpdata,kxx,4,7,0.53
mpdata,kxx,5,1,410,407,404,400,394,381
mpdata,kxx,5,7,367
mpdata,kxx,6,1,429,425,421,416,411,398
mpdata,kxx,6,7,384
mptemp
mptemp,1,20,27,127,327,527,727
mptemp,7,765,927
mpdata,enth,5,1,0,92609056,463285756,1.3748e9,2.3519e9,3.4396e9
mpdata,enth,5,7,3.648547e9,4.4520e9
mptemp
mptemp,1,20,100,200,300,400,500
mptemp,7,600
mpdata,c,2,1,386,391,402,413,426,433
mpdata,c,2,7,452
mpdata,c,6,1,240,245,254,263,273,279
mpdata,c,6,7,295
mp,c,3,390
allsel,
finish
/solu
antype,trans
toffst,273
tunif,30
cnvtol,heat,1
kbc,1
trnopt,full
autos,on
deltim,0.005,0.001,delt,on
outres,basic,all
physics,write,thermal
finish
*do,i,1,ftime/tinc
time=time+tinc
physics,read,magenti
/solu
*if,i,eq,1,then
tunif,30
*else
ldread,temp,last,,,,,rth
*endif
solve
finish
physics,read,thermal
/assign,esav,therm,esav
/assign,emat,therm,emat
/solu
*if,i,gt,1,then
antype,trans,rest
*endif
time,time
esel,s,mat,,2
esel,a,mat,,5,6,1
ldread,hgen,,,,2,,
esel,all
solve
finish
/assign,esav
/assign,emat
*enddo
finish
save,wendu,db
目前还没有考虑热辐射，考虑到计算的花费的时间，我目前把加热时间定为0.2，这样下来计算时间比较短，但我把电流设得比较大，设成10000A了。
麻烦大家帮我看看，问题出在哪里！
谢谢了！
|
IP卡
狗仔卡
发表于 2011-12-7 11:40:37
提升卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
