| 查看: 714 | 回复: 1 | ||
[求助]
fluent复合材料固化度UDS！！！
|
|
求解的方程： a为固化度，da/dt与a和温度Temp（K为温度的函数，T为固体温度）有关，我想获得a与时间的关系即ai+1=ai+da/dt*dtime  捕获.PNG 这个是我的模型，即空气给固体区域加热，固体区域随着温度的升高固化  slide2.png 这是我编写的UDS但是结算完全是错误的 代码如下 #include "udf.h" #include "math.h" #define A 6.167e20 #define H 77500.0 #define R 8.314 #define E 1.674e5 #define m 0.524 #define n 1.476 DEFINE_SOURCE(dc, c, t, dS, eqn)//固化速率da/dt { real Temp, source,a; Temp = C_T(c, t); a = C_UDSI(c, t, 0);//获取固化度 source = A/exp(E/R/Temp)*pow(a,m)*pow((1-a),n);//本例k2为0 dS[eqn] = 0; return source; } DEFINE_SOURCE(DOC, c, t, dS, eqn)//固化度 { real dc, source; dc=C_UDSI(c,t,1);//获取固话速率 source = C_UDSI(c, t, 0); source = source + dc*CURRENT_TIMESTEP; dS[eqn] = 0; return source; } DEFINE_PROFILE(tt,t,i)//速度入口温度随时间变化 { face_t f; begin_f_loop(f,t) { real time=RP_Get_Real("flow-time"  ;if (time<3400) { F_PROFILE(f,t,i) = 288+0.03*time; } else F_PROFILE(f,t,i) =399; } end_f_loop(f,t) } |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2020-03-30 18:24:36
