稳态热传递问题求指教
在做一个人二维热仿真问题，E型磁芯中间绕铜导线，简化为二维问题就是中间一些圆，磁芯周边时间30度的边界温度，导体圆面施加功率，想知道怎样验证施加的功率最后导致的温升是正确的。求各位老师指点。理论计算 小七能不能帮我算一下，我没学过传热学理论，可以的话我把模型给你。仿真币我全给都行
华山仙 发表于 2012-10-25 20:00 static/image/common/back.gif
小七能不能帮我算一下，我没学过传热学理论，可以的话我把模型给你。仿真币我全给都行
...
我也只是会模拟 公式好久么有看了 小七工作室 发表于 2012-10-25 20:02 static/image/common/back.gif
我也只是会模拟 公式好久么有看了
那我把我的仿真发上来，帮我看看好吗？这周要交差现在很着急，谢谢了！ 华山仙 发表于 2012-10-25 20:06 static/image/common/back.gif
那我把我的仿真发上来，帮我看看好吗？这周要交差现在很着急，谢谢了！ ...
Finish
/clear
/filname, inductance
/title, ZY_inductance
/PREP7
/number, 0
Pi=acos(-1)
EMUNIT, MKS
core_length=30.25e-3
core_depth=26.3e-3
all_window=19.85e-3
mid_length=10.65e-3
window_length=(all_window-mid_length)/2
window_depth=16.3e-3
left_length=(core_length-all_window)/2
up_depth=5e-3
mid_depth=5e-3
rad_win=0.5e-3
x_dis=1e-3
y_dis=1.5e-3
th_air=1.3e-3 !thickness of air
w_c_dis=1e-3 !distance from winding to core
th_win=5e-3!thickness of winding
ET,1,plane77 !windings
ET,2,plane77 !core
ET,3,plane77 !air
MP,KXX,1,398 !铜导线的热导率
MP,DENS,1,8900 !密度
MP,C,1,385 !比热
MP,KXX,2,0.55 !磁芯的热导率
MP,DENS,2,4700
MP,C,2,0.67
MP,KXX,3,0.02 !空气的热导率
MP,DENS,3,1.2
MP,C,3,1000
rectnag,0,core_length,0,core_depth+th_air
rectnag,left_length,left_length+window_length,up_depth+0.5e-3,up_depth+window_depth+th_air
rectnag,left_length+window_length+mid_length,left_length+2*window_length+mid_length,up_depth+0.5e-3,up_depth+window_depth+th_air
rectnag,0,core_length,up_depth,up_depth+th_air
allsel
aadd,2,3,4
numcmp,all
allsel
asba,1,2,,delete,keep
numcmp,all
numstr,area,50
cyl4,left_length+x_dis+2*rad_win,up_depth+y_dis+2*rad_win,rad_win,0,,360
agen,2,50,,,x_dis+rad_win,,,
agen,9,50,52,,,y_dis
agen,2,50,68,,window_length+mid_length-2*rad_win,,
allsel
aovlap,all
numcmp,all
allsel
asel,s,,,39
aatt,3,,3
smrtsize,3
amesh,all
allsel
asel,s,,,1
aatt,2,,2
smrtsize,3
amesh,all
allsel
asel,s,,,2
aatt,2,,2
smrtsize,3
amesh,all
allsel
asel,s,,,3,38
aatt,1,,1
smrtsize,3
amesh,all
allsel
nsel,s,loc,x,0
nsel,a,loc,x,core_length
nsel,a,loc,y,0
nsel,a,loc,y,core_depth+th_air
d,all,temp,30
allsel
nsel,s,node,,17763
nsel,a,node,,17990
nsel,a,node,,16981
nsel,a,node,,17363
nsel,a,node,,16646
nsel,a,node,,16770
nsel,a,node,,15724
nsel,a,node,,16176
nsel,a,node,,15397
nsel,a,node,,15519
nsel,a,node,,14827
nsel,a,node,,14937
nsel,a,node,,14046
nsel,a,node,,14418
nsel,a,node,,13490
nsel,a,node,,13868
nsel,a,node,,12940
nsel,a,node,,13150
nsel,a,node,,23487
nsel,a,node,,23560
nsel,a,node,,22816
nsel,a,node,,23224
nsel,a,node,,22302
nsel,a,node,,22402
nsel,a,node,,21603
nsel,a,node,,21931
nsel,a,node,,20962
nsel,a,node,,21367
nsel,a,node,,20214
nsel,a,node,,20499
nsel,a,node,,19612
nsel,a,node,,19904
nsel,a,node,,18962
nsel,a,node,,19231
nsel,a,node,,18351
nsel,a,node,,18615
f,all,heat,0.06
/solu
antype,static
nropt,auto
autots,on
nsubst,50
allsel
outres,all,all
solve
/post1
allsel
nsle,s
plnsol,temp
Finish
/clear
/filname, inductance
/title, ZY_inductance
/PREP7
/number, 0
Pi=acos(-1)
EMUNIT, MKS
core_length=30.25e-3
core_depth=26.3e-3
all_window=19.85e-3
mid_length=10.65e-3
window_length=(all_window-mid_length)/2
window_depth=16.3e-3
left_length=(core_length-all_window)/2
up_depth=5e-3
mid_depth=5e-3
rad_win=0.5e-3
x_dis=1e-3
y_dis=1.5e-3
th_air=1.3e-3 !thickness of air
w_c_dis=1e-3 !distance from winding to core
th_win=5e-3!thickness of winding
ET,1,plane77 !windings
ET,2,plane77 !core
ET,3,plane77 !air
MP,KXX,1,398 !铜导线的热导率
MP,DENS,1,8900 !密度
MP,C,1,385 !比热
MP,KXX,2,0.55 !磁芯的热导率
MP,DENS,2,4700
MP,C,2,0.67
MP,KXX,3,0.02 !空气的热导率
MP,DENS,3,1.2
MP,C,3,1000
rectnag,0,core_length,0,core_depth+th_air
rectnag,left_length,left_length+window_length,up_depth+0.5e-3,up_depth+window_depth+th_air
rectnag,left_length+window_length+mid_length,left_length+2*window_length+mid_length,up_depth+0.5e-3,up_depth+window_depth+th_air
rectnag,0,core_length,up_depth,up_depth+th_air
allsel
aadd,2,3,4
numcmp,all
allsel
asba,1,2,,delete,keep
numcmp,all
numstr,area,50
cyl4,left_length+x_dis+2*rad_win,up_depth+y_dis+2*rad_win,rad_win,0,,360
agen,2,50,,,x_dis+rad_win,,,
agen,9,50,52,,,y_dis
agen,2,50,68,,window_length+mid_length-2*rad_win,,
allsel
aovlap,all
numcmp,all
allsel
asel,s,,,39
aatt,3,,3
smrtsize,3
amesh,all
allsel
asel,s,,,1
aatt,2,,2
smrtsize,3
amesh,all
allsel
asel,s,,,2
aatt,2,,2
smrtsize,3
amesh,all
allsel
asel,s,,,3,38
aatt,1,,1
smrtsize,3
amesh,all
allsel
nsel,s,loc,x,0
nsel,a,loc,x,core_length
nsel,a,loc,y,0
nsel,a,loc,y,core_depth+th_air
d,all,temp,30
allsel
nsel,s,node,,17763
nsel,a,node,,17990
nsel,a,node,,16981
nsel,a,node,,17363
nsel,a,node,,16646
nsel,a,node,,16770
nsel,a,node,,15724
nsel,a,node,,16176
nsel,a,node,,15397
nsel,a,node,,15519
nsel,a,node,,14827
nsel,a,node,,14937
nsel,a,node,,14046
nsel,a,node,,14418
nsel,a,node,,13490
nsel,a,node,,13868
nsel,a,node,,12940
nsel,a,node,,13150
nsel,a,node,,23487
nsel,a,node,,23560
nsel,a,node,,22816
nsel,a,node,,23224
nsel,a,node,,22302
nsel,a,node,,22402
nsel,a,node,,21603
nsel,a,node,,21931
nsel,a,node,,20962
nsel,a,node,,21367
nsel,a,node,,20214
nsel,a,node,,20499
nsel,a,node,,19612
nsel,a,node,,19904
nsel,a,node,,18962
nsel,a,node,,19231
nsel,a,node,,18351
nsel,a,node,,18615
f,all,heat,0.06
/solu
antype,static
nropt,auto
autots,on
nsubst,50
allsel
outres,all,all
solve
/post1
allsel
nsle,s
plnsol,temp
就是在一个EI磁芯的中心柱绕铜导线，铜线上消耗的总功率2.3W,在最外边界施加恒温30度，最终目的是想看铜线的发热情况。 华山仙 发表于 2012-10-25 20:15 static/image/common/back.gif
Finish
/clear
/filname, inductance
小七能帮帮忙吗？我在线等！ 华山仙 发表于 2012-10-25 20:19 static/image/common/back.gif
小七能帮帮忙吗？我在线等！
不好意思刚才多粘贴了一回命令。 华山仙 发表于 2012-10-25 20:15 static/image/common/back.gif
Finish
/clear
/filname, inductance
C:\Documents and Settings\Administrator\桌面\新建文件夹 (3)
页:
[1]