孔隙率公式22/23? : a / a 修改意见? line 24: value: undeclared variable 20

浏览:747
#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)
}


邀请回答 我来回答

当前暂无回答

回答可获赠 200金币

没解决?试试专家一对一服务

换一批