高斯热源分布的udf 编程
[font=宋体]大家好，我是做焊接熔池温度场模拟的，所用软件为[/font][font=Times New Roman]FLUENT[/font][font=宋体]。热源选用高斯分布加到壁面上，可我不会对热源进行编程，希望哪位高手能帮帮忙，指点一下，非常感谢！热源公式为[/font][font=Times New Roman] [/font][font=宋体]高斯函数分布可表示为[/font][font=Times New Roman] q(r)=[/font][font=宋体]（[/font][font=Times New Roman]0.65[/font][font=宋体]×[/font][font=Times New Roman]14[/font][font=宋体]×[/font][font=Times New Roman]150[/font][font=宋体]）÷（π×[/font][font=Times New Roman]0.01[size=10.5pt]2[/size][/font][font=宋体]）×[/font][font=Times New Roman]exp(-(x[size=10.5pt]2[/size][/font][font=宋体]＋[/font][font=Times New Roman][size=10.5pt]z2[/size])[/font][font=宋体]÷[/font][font=Times New Roman]0.01[size=10.5pt]2[/size]) [/font] 我自己仿着帮助文档编了一个定点热源的udf，不知对不对，请各位高手帮忙指点一下，万分感谢！5qU\,pI ^!u K
#include "udf.h"2NW?6LoA!S"H \
DEFINE_PROFILE(immovable_heat,t,i)2FH9c"RB.N8z
{0eb:c%T D
real x[ND_ND];5P/z+r5]Z)a MZ S
real y,z;
face_t f;
begin_f_loop(f,t)bx0Ux$U+h'f+`&[
{^(Qj_en/bv:WN\
F_CENTROID(x,f,t); }#X V0Z2C,Ce7^
z=x[1];4rw-y Uq
y=x[0];
F_PROFILE(f,t,i)=0.65*16*110/(3.1415926*0.003*0.003)*exp(-(y*y+z*z)/(0.003*0.003);O ? Ab?
}W'q.j L:^ w)cH
end_f_loop(f,t)Y K ]]BZGI
}
页:
[1]