【分享】高斯分布热源下的两层材料热分析
试着做了高斯分布热源下的两层材料瞬态热分析，上层为金刚石薄膜，下层为SI底。以下是命令流，如果有不合适的地方，欢迎大家提出、讨论！其中高斯分布体热源是通过GUI方式加的 :（Parameters -->Functions-->Define/Edit编辑函数为：
1e13*exp(-({X}^2+{Y}^2))*exp({Z}-2e-4)*exp(-{TIME})
再Parameters -->Functions-->Read From File ，并命名为 myload）
/units,si
/prep7
ET,1,SOLID90
ET,2,MESH200,7
mp,kxx,1,148 !SI参数
mp,c,1,712
mp,dens,1,2330
!mp,kxx,2,2300 !Diamond参数
!mp,c,2,509
!mp,dens,2,3500
mp,kxx,2,170
mp,c,2,870
mp,dens,2,2700
CYLIND,1.25e-3,,2.0e-4,,45
CYLIND,1.25e-3,,2.0e-4,2.1e-4,,45
vglue,all
lesize,6,,,25,2
lesize,1,,,25,2
lesize,9,,,20,0.5
MSHAPE,0,2D
mshkey,1
MAT,1
TYPE,2
AMESH,4
allsel
type,1
mat,1
lesize,4,,,10
VSWEEP,1,4,5
lesize,15,,,25,2
lesize,6,,,25,2
lesize,21,,,2
MSHAPE,0,2D
mshkey,1
MAT,2
TYPE,2
AMESH,13
TYPE,1
MAT,2
LESIZE,13,,,10
VSWEEP,3,13,12
ALLSEL
SAVE
/SOLU
AnTYPE,TRANS
TUNIF,293
kbc,0
asel,s,area,,7
asel,a,area,,11
asel,a,area,,3
nsla,s,1
sf,all,conv,10,293
allsel
nsel,s,loc,z,0,2.0e-4
csys,1
nsel,r,loc,x,0,1.25e-3
csys,0
bf,all,hgen,%myload%
allsel
time,1
deltime,2e-3,2e-3,8e-3
autots,on
outres,all,2
solve
回复： 【分享】高斯分布热源下的两层材料热分析
1/4对称图。回复： 【分享】高斯分布热源下的两层材料热分析
1e13*exp(-({X}^2+{Y}^2))*exp({Z}-2e-4)*exp(-{TIME})请问一下你这个公式是什么意思？这是高斯分布热源的公式吗？我怎么觉得和书上讲的公式不一样呢？请详细说明一下好吗？
而且我运行你的程序时怎么得不出你所显示的结果？
回复： 【分享】高斯分布热源下的两层材料热分析-----------------------------7
我也是从一篇文章中看到这个公式的（脉冲激光下，物体吸收热的高斯分布的体热源公式，只不过前边还有一些常数，还做了时间归一化处理）。而且我运行你的程序时怎么得不出你所显示的结果？
Sorry，原公式应该是这样的：
1e13*exp(-({X}^2+{Y}^2)/9e-10)*exp(1e5*({Z}-2e-4))*({TIME}) *exp(-{TIME})
我在7.0下得到的，另外，不能将体热源在R和Z方向加的太多，否则会计算出散热的结果来，可能是由于R和Z方向太大，数字非常小（10e-50等），Ansys会溢出吧。
再试一试。
回复： 【分享】高斯分布热源下的两层材料热分析
是吗？不过我见到的高斯热源的公式可都是q=qmax*exp(-Kr^2),没见过那么复杂的我还想请教一下
nsel,s,loc,z,0,2.0e-4
csys,1
nsel,r,loc,x,0,1.25e-3
csys,0
这是选的哪一段区域呀？
我想做一个移动的热源，可是在选择实体的时候，我用nsel,s,loc,x这个命令时默认的坐标系为整体坐标系，但整体坐标系是没法移动的，所以如果我的坐标系的原点和要加载的区域的圆心重合的话还可以，如果热源移动，就没法选择了，移动局部坐标系好像没什么效果。到底该如何select呢？请指教
回复： 【分享】高斯分布热源下的两层材料热分析
nsel,s,loc,z,0,2.0e-4csys,1
nsel,r,loc,x,0,1.25e-3
csys,0
选择z方向（0－200 um)且半径方向（R方向0－1.25 mm）区域。选择R方向时，要将坐标系转化为柱坐标下进行（csys，1）。
所以如果我的坐标系的原点和要加载的区域的圆心重合的话还可以，如果热源移动，就没法选择了，移动局部坐标系好像没什么效果。到底该如何select呢？
不明白，讲的详细点吧。
回复： 【分享】高斯分布热源下的两层材料热分析
我想做一个移动热源的例子，可是不知该如何加高斯热源，我想用你这种方法，可是我不知道随着热源的移动，该如何选择实体，不知说明白了没有，如果你知道麻烦你把命令流写出来，谢谢！下面这个链接是一个关于高斯移动热源的例子，我没看懂，我觉得这样做太复杂，你能帮我看看吗？难道真得那么难吗？
http://www.simwe.com/cgi-bin/ut/topic_show.cgi?id=12577&pg=3&h=1&bpg=2&age=0
回复： 【分享】高斯分布热源下的两层材料热分析
To: pzl 好像专门有考虑材料间的界面热阻的单元，在版上搜索一下吧。
To：jnjl
版上的高斯移动热源的例子，用Apd编程，对一般问题都适用。不知道你的移动热源移动规律复不复杂？如果比较简单的话，也许通过GUI方式加载方便一些：
比如热源沿Z方向移动，并且知道移动速度为V，则函数可写为：
Q* exp(-({Z}-V*{TIME}))
这样相当于在每一个时刻，都对热源函数做了更新。由于GUI方式加载时，Ansys在每一个时间步都对函数进行计算，加载，所以也相当于移动了热源。结合你的高斯函数和以上例子，应该不会太难。试一试吧。
这两天在写论文，挺忙的。
回复： 【分享】高斯分布热源下的两层材料热分析
谢谢，非常感谢！回复： 【分享】高斯分布热源下的两层材料热分析
不好意思，再打扰你一下你的例子中加的生热率也是移动的吗？
这是我做的一个例子，关于加高斯热源我还是搞不定，麻烦你帮我看一下，该如何加？
问题描述：
已知板厚3mm,长500mm，宽50mm,现在其长度方向（Z方向）上加一移动高斯热源，速度为0.5m/min，热源的宽度为6mm,热源从z=0处移动到z=20mm处，根据高斯分布公式Q=Qmax*exp(-kr^2),得出Q的函数表达式为Q=58360/exp((10^6)/100*({X}^2+{z}^2)) ，热流密度为q/3.14159*r2,列出命令流如下：
首先利用GUI方式加高斯分布热源
58360/3.14159/({X}^2+{z}^2)/exp((10^6)/100*({X}^2+{z}^2))* exp(-1000*({Z}-0.0083*{TIME}))
再Parameters -->Functions-->Read From File ，并命名为 gaload
/filname,ga1,0
/title,3d moving heat
/prep7
et,1,solid70
keyopt,1,2,0
keyopt,1,4,0
keyopt,1,7,0
keyopt,1,8,0
!*
mptemp,1,0
mpdata,kxx,1,,70
mptemp,1,0
mpdata,c,1,,448
mptemp,1,0
mpdata,dens,1,,7833
block,-0.025,0.025,-0.003,0,0,0.5,
block,-0.005,0.005,-0.003,0,0,0.5,
vsel,s,volu,,1,2
vovlap,all
vglue,all
numcmp,line
numcmp,area
numcmp,volu
/replot
aesize,all,0.002,
lesize,all,0.002, , , ,1, , ,1,
esize,0.002,0,
mshkey,1
mshape,0,3d
vmesh,2,3
aesize,all,0.0005,
lesize,all,0.0005, , , ,1, , ,1,
esize,0.0005,0,
mshkey,1
mshape,0,3d
vmesh,1
/view, 1 ,1,1,1
/ang, 1
/rep,fast
/psf,hflu,,2,0
/pnum,line,1
sfa,1,,conv,50,25,!w/m**2*c
sfa,2,,conv,50,25,
sfa,12,,conv,50,25,
sfa,16,,conv,50,25,
sfa,11,,conv,12.5,25,!bottom
sfa,5,,conv,12.5,25,
sfa,15,,conv,12.5,25,
finish
/solu
antype,trans
trnopt,full
nropt,auto,,
eqslv,
solcontrol,on
autots,on
kbc,1 !step
outres,all,all,
!load step 1, initial conditions 15 c
time,0.001
deltim,0.001,0.001,0.001
timint,off !time integration = off
solve
kbc,0
tunif,15,
!load step 2..... , apply moving heat flux
timint,on !time integration = on
allsel,all
ddele,all,temp
sfdele,all,hflux
!*
!以下四行为加高斯热源处，请帮忙改正一下
nsel,s,loc,x,-0.003,0.003
nsel,r,loc,y,0
nsel,r,loc,z,0,0.02
sf,all,hflu,%gaload%
!*
allsel,all
time,2.4
deltim,0.1,0.1,0.1,
outres,all,2
allsel,all
solve
alls
finish
/post1
/eface,1
set,last
plnsol,temp,
回复： 【分享】高斯分布热源下的两层材料热分析
freund: 怎么计算高斯分布热源下的两层材料热应力分布,能否把命令流传上来,非常感激!!
我信箱: zyp_820@163.com
回复： 【分享】高斯分布热源下的两层材料热分析
谁还会做啊??????????? 呵呵都是不简单！ 有个疑问，初始温度为293k，温升结果怎么那么小？上升了不到3k?
页:
[1]
2