|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
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;
}
|
|