| 查看: 2685 | 回复: 0 | ||
[求助]
DEFINE_GEOM这个宏怎么用？怎么在fluent中实现曲线拟合，求个例子
|
|
DEFINE_GEOM(nihe,domain,dt,position)这个宏怎么用？怎么在fluent中实现曲线拟合，求个例子 我写了段程序fluent总是出现负体积，程序如下： DEFINE_GEOM(nihe,domain,dt,position) { real x1=0; real x2=0; real x3=0; real x4=0; real y1=0; real xy=0; real x2y=0; real a=0; real b=0; real c=0; real d=0; real m=80; real a0; real a1; real a2; real x01,x02,x03,x11,x12,x13,x21,x22,x23; x1+=position[0]; x2+=position[0]*position[0]; x3+=position[0]*position[0]*position[0]; x4+=position[0]*position[0]*position[0]*position[0]; y1+=position[1]; xy+=position[0]*position[1]; x2y+=position[0]*position[0]*position[1]; a=-1/m*x1*x1+x2; b=-1/m*x1*x2+x3; c=-b*b/a-1/m*x2*x2+x4; d=-b/a/m*x1+1/m*x2; x01=-d/c/a/m*(b*x1-a*x2)+x1*x1/a/m/m+1/m; x02=b*d/a/c-x1/a/m; x03=-d/c; x11=-b/a/a/c/m*(b*x1-a*x2)-x1/a/m; x12=b*b/a/a/c+1/a; x13=-b/a/c; x21=(b*x1/a/m-x2/m)/c; x22=-b/a/c; x23=1/c; a0= x01*y1+x02*xy+x03*x2y; a1= x11*y1+x12*xy+x13*x2y; a2= x21*y1+x22*xy+x23*x2y; Message(" a0=\n", a0 ) ; Message(" a1=\n", a1 ) ; Message(" a2=\n", a2) ; position[1] = a2*position[0]*position[0]+a1*position[0]+a0; } |
回复此楼» 猜你喜欢
申博发邮件
已经有5人回复
-
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后
已经有5人回复
-
国社科系统bug了，是不是要放榜了？
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
点击这里搜索更多相关资源