24小时热门版块排行榜    

查看: 734  |  回复: 1

yumk41

新虫 (初入文坛)

[求助] 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)
}
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

k剑云一木

新虫 (初入文坛)

请问楼主解决了吗?我也遇到同样的问题,知道你是怎样利用uds求解固化度的?
2楼2020-03-30 18:24:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yumk41 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:39 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:20 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:34 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-02 05:34 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:27 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 X5bHxOULxFqb 2026-10-30 22/1100 2026-11-01 16:35 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 12:27 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
信息提示
请填处理意见