| 查看: 9932 | 回复: 19 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
fluent中的体热源加载问题（实际上是热流密度，不是热生成率）
|
|||
|
做焊接激光模拟，采用高斯旋转体热源。随时间移动。融化凝固模型。 查询很多资料和问了若干人后，他们的建议都是加载在能量源项。 /*采用旋转高斯体热源*/ /*能量源项*/ #include "udf.h" DEFINE_SOURCE(heat_flux, cell, thread, dS, eqn) { real x[ND_ND]; real y,I,U,oq,t0,v0,Q,time,cs,HH; real r,rh; real source_heat; cell_t c; C_CENTROID(x,c,thread); y=0.7; /*有效功率因子*/ U=16.; I=100.; v0=2.e-3; /*焊接速度*/ Q=y*I*U; /*有效功率*/ cs=3./pow(oq,2.); /*形状因子*/ HH=1.e-3; /*热源高度*/ oq=2.5e-3; /*热源半径*/ time= RP_Get_Real("flow-time"  ;r=sqrt(pow(x[0]-v0*time,2.)+pow(x[1],2.)); /*坐标与热源中心的距离*/ rh=sqrt(log(HH/x[2])/3.)*oq;/*一定热源高度对应的半径rh*/ if (x[2]==0&&r<=oq) { source_heat=3.*cs*Q/(3.14*HH*(1.-1./exp(3)))*exp(-3.*pow(r,2.)/pow(oq,2.)); dS[eqn] =0; } else if (x[2]==HH&&r==0.) { source_heat=3.*cs*Q/(3.14*HH*(1.-1./exp(3))); dS[eqn] =0; } else if (x[2]>HH) { source_heat=0; dS[eqn] =0; } else if (r<=rh) { source_heat=3.*cs*Q/(3.14*HH*(1.-1./exp(3)))*exp(-3.*cs/log10(HH/x[2])*pow(r,2.)); dS[eqn] =0; } else { source_heat= dS[eqn] =0; } return source_heat; } 条件语句比较复杂，主要怕出问题。呵呵 interpret后没有错误，但加载出不了结果。 希望大家给点意见。 [ 来自科研家族 流体数值模拟 ] [ Last edited by b07988 on 2012-3-4 at 22:31 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 电弧 |
» 猜你喜欢
Cas 72-43-5需要30g，定制合成，能接单的留言
已经有8人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有6人回复
-
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有8人回复
-
磺酰氟产物，毕不了业了！
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有25人回复
-
2026年机械制造与材料应用国际会议 （ICMMMA 2026）
已经有3人回复
-
自荐读博
已经有3人回复
-
不自信的我
已经有5人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助fluent中速度显示问题
已经有22人回复
- fluent中关于拉瓦尔喷管的问题
已经有26人回复
- 求助：fluent中关于壁面传热问题
已经有12人回复
- fluent中的数据统计问题
已经有9人回复
- 关于fluent中残差震荡的问题~诚求各位大侠
已经有10人回复
- fluent中壁面热边界条件的问题
已经有16人回复
- fluent气液双相流中的Patch问题
已经有15人回复
- 【求助】fluent 液体仿真的时候热源怎么加载
已经有8人回复
- 【求助】Fluent中的对称问题
已经有10人回复
- 【求助】fluent中固体与流体区域间的interior界面，出现警告提示该怎么解决
已经有6人回复
- 【求助】FLUENT三维模型中入口速度的方向设置问题【已完成】
已经有17人回复
- 【求助】FLUENT中的内热源怎么设置？在哪里设置？
已经有16人回复
- 【求助】有关fluent中DPM模型的模拟问题
已经有16人回复
10楼2012-11-30 11:23:27