| 查看: 2673 | 回复: 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; } |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有12人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
点击这里搜索更多相关资源
