| 查看: 10169 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 电弧 |
» 猜你喜欢
材料与化工371求调剂
已经有5人回复
-
328分调剂
已经有4人回复
-
295求调剂
已经有7人回复
-
081700，311，求调剂
已经有13人回复
-
283分求调剂
已经有6人回复
-
0854求调剂
已经有3人回复
-
332求调剂
已经有4人回复
-
284求调剂
已经有5人回复
-
304求调剂（085602，过四级，一志愿985）
已经有11人回复
-
0817化学工程与技术求调剂，一志愿中海洋319
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助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人回复
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
b07988(金币+3): ★有帮助 给出了若干建议，谢谢 2012-03-05 12:30:54
xiegangmai(金币+2): 鼓励讨论交流！ 2012-03-05 20:26:19
感谢参与，应助指数 +1
b07988(金币+3): ★有帮助 给出了若干建议，谢谢 2012-03-05 12:30:54
xiegangmai(金币+2): 鼓励讨论交流！ 2012-03-05 20:26:19
|
1, 我记得传热学一个默认的假设是热扰动传递速率无限大，但对与激光加热等加热时间极短的问题，该假设不在适用，这样常见的传热学中微分方程就要做相应的调整，不知道你的激光焊接问题村不存在这个问题？ 2 条件语句中有几个将y坐标和某一值做等值比较，首先C_CENTROID取的是网格中心的坐标值，是一系列离散的值，你等确定你的网格正好有一一个或几个网格的y坐标值是这个值（个人觉得，除非你自己输入坐标生成网格，利用网格生成软件生成的网格是不能保证的）？也就是说你那个几个有等号的条件判断语句有很大的可能是遍历完整个网格也找不到符合条件的，这应该是加载不出结果的原因。所以需要修改判断条件。 3，对于1，请教一下高斯热源定义。 |
2楼2012-03-04 22:58:13
★ ★
xiegangmai(金币+2): 鼓励讨论交流！ 2012-03-05 20:26:43
xiegangmai(金币+2): 鼓励讨论交流！ 2012-03-05 20:26:43
|
1.你所说的加热时间极短，v0=2.e-3; /*焊接速度*/，也就是热源移动的速度，我觉得应该不会有影响。我所写的代码中没有微分方程。除了dS[eqn] =0; 2.这个问题确实没仔细考虑，但仔细想想这样写的等式其实没什么问题。热源是一个倒锥形。所涉及的条件是刚好在那个锥尖处。网格比较密，产生的误差可以忽略。就算关键的地方没有满足坐标要求，但仍旧有部分热源存在的。但实际加载没有效果。 3.如图，因为表达式存在1/log(H/z),其中z不能等于0，H不能等于z,所以条件语句比较复杂。呵呵 残差图就是直线，基本没有波动，说明热源没有加载上。  热源模型表达式  中心那根线的热流密度（处处相等）  整体形状 |
3楼2012-03-05 12:48:52
6楼2012-03-05 18:09:05
【答案】应助回帖
★ ★ ★ ★ ★ ★
b07988(金币+4): ★★★很有帮助 谢谢你的建议，用message做判断是个好方法 2012-03-05 16:09:04
xiegangmai(金币+2): 谢谢参与！ 2012-03-05 20:27:08
b07988(金币+4): ★★★很有帮助 谢谢你的建议，用message做判断是个好方法 2012-03-05 16:09:04
xiegangmai(金币+2): 谢谢参与！ 2012-03-05 20:27:08
|
我没注意到你还有r<=rh这一句。 至于那几个等式比较的语句，并不是你所说的网格较密，误差可忽略。fluent数据是real类型的，小数点位数应该在六位以上，因为是等值比较，只需要最后以为小数不等，条件就不成立。 你可以在你设的有热源的条件语句部分加一个Message语句，把该点的坐标值输出到控制窗口里，根据这个判断，而非根据残差曲线。 |
4楼2012-03-05 15:30:36
5楼2012-03-05 16:19:29
7楼2012-03-06 20:36:41
8楼2012-03-07 12:15:11
9楼2012-11-28 11:20:18
10楼2012-11-30 11:23:27
