| 查看: 524 | 回复: 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); } |
回复此楼» 猜你喜欢
材料调剂
已经有4人回复
-
面上模板改不了页边距吧？
已经有6人回复
-
307求调剂
已经有6人回复
-
304求调剂
已经有5人回复
-
317一志愿华南理工电气工程求调剂
已经有8人回复
-
272求调剂
已经有3人回复
-
化工专硕348，一志愿985求调剂
已经有6人回复
-
292求调剂
已经有3人回复
-
290求调剂
已经有6人回复
-
295求调剂
已经有5人回复
