孔隙率公式22/23？ : a / a 修改意见？ line 24: value: undeclared variable 20
浏览：816
#include "udf.h"
#define X_L 180 // 采空区模型的长度
#define L 15 // 基本顶破碎长度
#define Y_L 120 // 工作面宽度
#define DP 0.15 // 平均粒径
#define T 30 // 倾角
#define TIMESTAMP "06301326" // 时间戳
// 绝对值函数
double my_abs(double x) {
if (x < 0) return -x;
else return x;
}
// 单点的孔隙率函数
double porous(double x, double y, double z) {
{
double valu,va;
valu =((1 + exp(-0.15 * (Y_L / 2 - my_abs(y - Y_L / 2))))) * (1 - 6 / (9.6 - 3.528 * (1 - exp(-x / (L * 2)))));
va =(valu -1);
double value =(va/( 1-0.001464*sin(T)*(Y_L / 2-y))+1);
value = sqrt(value);
return sqrt(value) * 0.8;
}
// 空隙率
DEFINE_PROFILE(porous_profile, thread, position)
{
real r[ND_ND];
real x, y, z, value;
cell_t c;
begin_c_loop(c, thread)
{
C_CENTROID(r, c, thread);
x = r[0];
y = r[1];
z = r[2];
value = 1 * porous(x, y, z);
C_PROFILE(c, thread, position) = value;
}
end_c_loop(c, thread)
}
