| 查看: 520 | 回复: 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); } |
回复此楼» 猜你喜欢
有院领导为了换新车，用横向课题经费买了俩车
已经有9人回复
-
酰胺脱乙酰基
已经有13人回复
-
博士延得我，科研能力直往上蹿
已经有8人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有5人回复
-
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优
已经有6人回复
-
面上项目申报
已经有3人回复
-
遇见不省心的家人很难过
已经有22人回复