|
|
[求助]
fluent复合材料固化度UDS!!!
求解的方程:
a为固化度,da/dt与a和温度Temp(K为温度的函数,T为固体温度)有关,我想获得a与时间的关系即ai+1=ai+da/dt*dtime
![fluent复合材料固化度UDS!!!]()
捕获.PNG
这个是我的模型,即空气给固体区域加热,固体区域随着温度的升高固化
![fluent复合材料固化度UDS!!!-1]()
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一区文章,我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
现代”学阀”该如何界定
已经有14人回复
师弟论文见刊大半年才想起来申请专利,还能抢救一下吗?
已经有4人回复
|