找回密码
 注册
查看: 798|回复: 2

运行网上一个流体程序,但为啥对话框说它11行有没被定义的量,有没有大神帮忙看看....

[复制链接]
发表于 2022-9-21 11:27:28 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
这是一个自定义材料的粘度程序如下,也许对你有帮助。

在记事本中编辑的,另存为"visosity1.c"
#include "udf.h"

DEFINE_PROPERTY(cell_viscosity,cell, thread)

{

real mu_lam;

real trial;

rate=CELL_STRAIN_RATE_MAG(cell, thread);

real temp=C_T(cell, thread);

mu_lam=1.e12;

{

if(rate>1.0e-4 && rate<1.e5)

trial=12830000./rate*log(pow((rate*exp(17440.46/temp)/1.535146e8),0.2817)+pow((1.+pow((rate*exp(17440.46/temp)/1.535146e8),0.5634)),0.5));

else if (rate>=1.e5)

trial=128.3*log(pow((exp(17440.46/temp)/1.535146e8),0.2817)+pow((1.+pow((exp(17440.46/temp)/1.535146e8),0.5634)),0.5));

else

trial=1.283e11*log(pow((exp(17440.46/temp)/1.535146e12),0.2817)+pow((1.+pow((exp(17440.46/temp)/1.535146e12),0.5634)),0.5));

}

else if(temp>=855.&&temp<905.)

{

if(rate>1.0e-4 && rate<1.e5)



trial=12830000./rate*log(pow((rate*4.7063),0.2817)+pow((1.+pow((rate*4.7063),0.5634)),0.5))*pow(10.,-0.06*(temp-855.));

else if (rate>=1.e5)

trial=243.654*pow(10.,-0.06*(temp-855.));

else

trial=1.47897e10*pow(10.,-0.06*(temp-855.));

}

else if(temp>=905.)

{

if(rate>1.0e-4 && rate<1.e5)

trial=12830./rate*log(pow((rate*4.7063),0.2817)+pow((1.+pow((rate*4.7063),0.5634)),0.5));

else if (rate>=1.e5)

trial=0.24365;

else

trial=1.47897e7;

}

if(trial<1.e12&&trial>100.)

mu_lam=trial;

else if(trial<=1.)

mu_lam=1.;

else

mu_lam=1.e12;

return mu_lam;

}
 楼主| 发表于 2022-9-21 11:28:19 | 显示全部楼层
上面的应该都是英文不知道咋回事全变中文了
 楼主| 发表于 2022-9-21 11:28:57 | 显示全部楼层
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表