| 查看: 2621 | 回复: 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; } |
回复此楼» 猜你喜欢
本9一志愿2 0854低分专硕286求调剂
已经有8人回复
-
070300化学求调剂
已经有15人回复
-
329求调剂
已经有12人回复
-
085701求调剂
已经有6人回复
-
26考研调剂0710 0860
已经有7人回复
-
材料专业383求调剂
已经有8人回复
-
一志愿北京科技大学材料工程085601，求调剂
已经有17人回复
-
一志愿郑州大学材料与化工085600，求调剂
已经有19人回复
-
085600，专业课化工原理，321分求调剂
已经有8人回复
-
求生物学调剂
已经有11人回复
点击这里搜索更多相关资源
