| 查看: 545 | 回复: 0 | ||
18428346586铁虫 (初入文坛)
|
[求助]
动网格UDF出错
|
|
FLUENT并行运算导入udf之后，动网格预览就出现如下错误提示：received a fatal signal (Segmentation fault)，请问为什么？ #include"udf.h" #include"dynamesh_tools.h" #define R 90.01 #define r 82 #define e 8 #define w 22.96 static real F(real time,real x,real y) { real alpha,a,b,c; alpha=w*time; a=1+pow(y/x,2); b=-2*e*(cos(alpha)+sin(alpha)*y/x); c=e*e-r*r; if(x>0) { return sqrt(pow((-b+sqrt(b*b-4*a*c))/(2*a),2)+pow((-b+sqrt(b*b-4*a*c))/(2*a)*(y/x),2)); } else if(x=0) { return sqrt(r*r-e*sin(alpha)*e*sin(alpha))-e*sin(alpha); } else { return sqrt(pow((-b-sqrt(b*b-4*a*c))/(2*a),2)+pow((-b-sqrt(b*b-4*a*c))/(2*a)*(y/x),2)); } } DEFINE_GRID_MOTION(motion,domain,dt,time,dtime) { Thread *tf=DT_THREAD((Dynamic_Thread*)dt); face_t f; Node *v; real x,y,z; int n; SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf)); begin_f_loop(f,tf) { f_node_loop(f,tf,n); { v=F_NODE(f,tf,n); x=NODE_X(v); y=NODE_Y(v); z=NODE_Z(v); if(NODE_POS_NEED_UPDATE(v)) { NODE_POS_NEED_UPDATE(v); NODE_X(v)=(R-(R-F(time+dtime,x,y))/(R-F(time,x,y))*(R-sqrt(x*x+y*y)))*x/sqrt(x*x+y*y); NODE_X(v)=(R-(R-F(time+dtime,x,y))/(R-F(time,x,y))*(R-sqrt(x*x+y*y)))*y/sqrt(x*x+y*y); NODE_Z(v)=z; } } } end_f_loop(f,tf); } |
回复此楼» 猜你喜欢
+
已经有3人回复
-
前几天时间戳更新了
已经有11人回复
-
＠_＠
已经有4人回复
-
*^O^*
已经有3人回复
-
重庆理工大学副校长遇刺身亡 传涉案副教授疑因积怨行凶
已经有6人回复
-
面上项目没有好文章就没希望了吗？
已经有20人回复
-
困死了
已经有8人回复
-
还有课题组有博士名额吗
已经有6人回复
-
博士申请
已经有3人回复