| 查看: 570 | 回复: 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); } |
回复此楼» 猜你喜欢
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
2026年叶企孙基金
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有13人回复
-
基金系统什么内容也没有
已经有9人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有11人回复
-
怎么查啊
已经有6人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有6人回复
-
面上合作单位盖章
已经有5人回复
投票:
30