| 查看: 2311 | 回复: 7 | |||
[求助]
激光焊接udf，加载温度场不对 已有1人参与
|
|
最近小弟在做激光焊接的流场模拟，属于初学者。编了一个udf，T型接头双侧焊接，通过编译，但是计算之后最高温度却只有320K左右，流场曲线更是不对，现在附上udf程序，希望大家帮助小弟解决一下，不胜感激！ #include "udf.h" #define P 1800 //功率 #define A 0.95 //系数 #define B 0.95 //系数 #define C 0.6 //系数 #define V 0.05 //焊接速度 #define PI 3.14159 DEFINE_SOURCE(xmom_source,c,t,dS,eqn) { int xmom; real X[ND_ND]; real source; real x, y, z, x1, y1, z1, x2, y2, z2; real g1, h1, j1, g2, h2, j2; float z3, z4; // 匙孔深度 real bate; //衰减系数 real Im, Ib, I, PB; real r1, r2, rb, R0; double eta1, eta2; real temp; temp=RP_Get_Real("flow-time"  ; bate=0.01; rb=0.00015+(P-1200)/800*0.000015; R0=10*rb; PB=A*B*C*P; I=8.6E10; Ib=2*bate*PB/(PI*rb*rb); Im=Ib+I; //功率密度 z3=log(Ib/I+1)/bate; z4=log(Ib/I+1)/bate; eta1=135*PI/180; eta2=-135*PI/180; C_CENTROID(X, c, t); x=X[0]; y=X[1]; z=X[2]; x1=x; y1=y-0.0024; z1=z-0.0015; x2=x; y2=y+0.0024; z2=z-0.0015; g1=x1; h1=cos(eta1)*y1+sin(eta1)*z1; j1=cos(eta1)*z1-sin(eta1)*y1; g2=x2; h2=cos(eta2)*y2+sin(eta2)*z2; j2=cos(eta2)*z2-sin(eta2)*y2; r1=sqrt((g1-V*temp)*(g1-V*temp)+h1*h1); r2=sqrt((g2-V*temp)*(g2-V*temp)+h2*h2); //之上为坐标平移及旋转 if(r1>0.00015&&j1==0.0015) { rb=10*rb; source=0.7*2*0.5*P/(PI*rb*rb)*exp(-3*r1*r1/(rb*rb)); dS[eqn]=0; } else if(r1<=0.00015&&j1>-0.0018&&j1<0.0015) { source=Im*exp(-2*r1*r1/(rb*rb))*exp(-bate*z3); dS[eqn]=0; } else if(r2>0.00015&&j1==0.0015) { rb=10*rb; source=0.7*2*0.5*P/(PI*rb*rb)*exp(-3*r2*r2/(rb*rb)); dS[eqn]=0; } else if(r2<=0.00015&&j2>-0.0018&&j2<0.0015) { source=Im*exp(-2*r2*r2/(rb*rb))*exp(-bate*z4); dS[eqn]=0; } else source=dS[eqn]=0; return source; } |
回复此楼» 收录本帖的淘帖专辑推荐
 电弧 |
» 猜你喜欢
301求调剂
已经有18人回复
-
化工学硕 285求调剂
已经有19人回复
-
328求调剂
已经有8人回复
-
求调剂
已经有28人回复
-
304求调剂（085602，过四级，一志愿985）
已经有21人回复
-
288环境专硕,求调材料方向
已经有33人回复
-
285求调剂
已经有16人回复
-
273求调剂
已经有41人回复
-
331求调剂
已经有4人回复
-
298求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 导入udf报致命错误！！！
已经有4人回复
- 曳力函数的的UDF
已经有15人回复
- udf profile update interval 是什么意思？
已经有5人回复
- 英、美、德，哪个学校、研究所焊接（激光焊接/激光复合焊接/金属材料）专业强呢？
已经有7人回复
- UDF应用之后的初始化错误
已经有6人回复
- udf 重力方向变化
已经有17人回复
- 【求助】udf读取混合物密度失败、及平均流速读取？
已经有8人回复
- 【求助】单机多核fluent并行udf
已经有13人回复
- 【求助】Warning: incorrect cg motion UDF value on zone 35
已经有5人回复
- 【交流】激光焊接有哪些前沿性的研究热点呢？
已经有9人回复
- 【讨论】请问大家用UDF的时候，有没有碰到如下代码，很重要的
已经有7人回复
- 【求助】急求激光焊接钛板引弧烧穿的解决方法
已经有7人回复
- 【求助】fluent计算中，udf编译通过，初始化出错
已经有11人回复
- 【求助】物料属性可以用UDF修改吗
已经有8人回复
2楼2012-04-05 15:30:00
3楼2012-04-05 15:48:40
yuchangandy
金虫 (正式写手)
- 应助: 12 (小学生)
- 金币: 3539.1
- 散金: 51
- 红花: 3
- 帖子: 333
- 在线: 615.8小时
- 虫号: 1162965
- 注册: 2010-12-05
- 专业: 金属基复合材料
4楼2012-04-05 15:49:42
5楼2012-04-05 15:50:56
lucklhz
木虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 1447.4
- 散金: 500
- 红花: 6
- 帖子: 753
- 在线: 385.9小时
- 虫号: 1291397
- 注册: 2011-05-10
- 性别: GG
- 专业: 工程热物理相关交叉领域
6楼2012-07-24 14:38:15
7楼2012-07-26 15:26:38
8楼2014-04-09 19:57:06
