| 查看: 1824 | 回复: 3 | |||
[求助]
刚体附近网格随刚体一起运动 已有1人参与
|
|
利用udf计算动网格问题，我想要做的是使得结构体附近的网格随刚体一起运动。 udf介绍：该udf使用newmark算法进行并行运算，使用了define cg motion,并将刚体运动后的坐标输出到txt文件中。 问题：编写的udf没有问题，能够使得结构体运动，但是当把附近网格（此部分网格已经单独划分并命名为fb），将该网格设置为rigid body 使用同一个udf，发现出现运行错误。也就是不设置fb为rigid body时候可以运行，设置成rigid body时候不能运行。 文献中提到，是要编写两个udf来实现。 求助第二个udf该怎么写？ |
回复此楼» 本帖已获得的红花（最新10朵）
小黑牛小白妞» 猜你喜欢
生物学296求调剂
已经有10人回复
-
288资源与环境专硕求调剂，不限专业，有学上就行
已经有19人回复
-
086000生物与医药 初试274求调剂
已经有4人回复
-
315求调剂
已经有7人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有11人回复
-
求调剂
已经有9人回复
-
309求调剂
已经有16人回复
-
322求调剂：一志愿湖南大学 材料与化工（085600），已过六级。
已经有10人回复
-
330分求调剂
已经有6人回复
-
317求调剂
已经有13人回复
2楼2016-07-08 09:58:41
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
木偶小天: 金币+10 2016-07-08 20:53:23
感谢参与，应助指数 +1
木偶小天: 金币+10 2016-07-08 20:53:23
|
#include <stdio.h> #include <math.h> #include "udf.h" #define WALLID 44 #define NUM_CALLS 7 DEFINE_CG_MOTION(cylinder, dt, vel, omega, time, dtime) { FILE *fp=NULL; char filename[]="out.txt"; Domain *domain=Get_Domain(1); Thread *t=Lookup_Thread(domain, WALLID); face_t f; static real lift=0.0,drag=0,lift_coeff=0.0,drag_coeff=0.0; static real a0=0.0,a1=0.0,a2=0.0,a3=0.0,a4=0.0,a5=0.0,a6=0.0,a7=0.0,KK=0.0,DeltT=0.0; static real Y=0.0,YY=0.0,Y1=0.0,YY1=0.0,YYY=0.0,PPt=0.0,Y0=0.0,YY0=0.0,YYY0=0.0; static int last_call=0; static int calls=0; int i; real rho=998.2, U=0.015, Dia=0.01; real m=0.2, k=0.710612, c=0; real Beta=0.25, Gama=0.5; real A[ND_ND],press_forces[ND_ND],viscous_forces[ND_ND]; real x_force=0.0, y_force=0.0; if ((++calls)==NUM_CALLS) last_call=1; if (last_call) { NV_S(A,=,0.0); begin_f_loop(f,t) { F_AREA(A,f,t); press_forces[0]+=F_P(f,t)*A[0]; press_forces[1]+=F_P(f,t)*A[1]; viscous_forces[0]+=F_STORAGE_R_N3V(f,t,SV_WALL_SHEAR)[0]; viscous_forces[1]+=F_STORAGE_R_N3V(f,t,SV_WALL_SHEAR)[1]; } end_f_loop(f,t) drag=press_forces[0]-viscous_forces[0]; lift=press_forces[1]-viscous_forces[1]; drag_coeff=drag/(0.5*rho*U*U*Dia); lift_coeff=lift/(0.5*rho*U*U*Dia); DeltT=dtime; a0 =1/(Beta*DeltT*DeltT); a1=Gama/(Beta*DeltT); a2=1/(Beta*DeltT); a3=1/(2*Beta)-1; a4=Gama/Beta-1; a5=(Gama/(2*Beta)-1)*DeltT; a6=(1-Gama)*DeltT; a7=DeltT*Gama; KK=k+a0*m+a1*c; PPt=lift+m*(a0*Y0+a2*YY0+a3*YYY0)+c*(a1*Y0+a4*YY0+a5*YYY0); Y=PPt/KK; YYY=(Y-Y0)*a0-a2*YY0-YYY0*a3; YY=YY0+(1-Gama)*DeltT*YYY0+Gama*DeltT*YYY; Y0=Y; YY0=YY; YYY0=YYY; Message("\n ==================================================================================="  ;Message("\n %-12s %-12s %-12s %-12s %-12s %-12s ","Flow_time","Y","lift_coeff","drag_coeff","lift","drag" ;Message("\n %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f ",time,Y,lift_coeff,drag_coeff,lift,drag); if ((fp=fopen(filename,"a" )==NULL)Message("\n Warning: Unable to open %s for writing!!!\n",filename); else { Message("\n Writing datas to %s...",filename); fprintf(fp,"%12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f %12.8f \n",time,Y,lift_coeff,drag_coeff,lift,drag,PPt,KK); fclose(fp); Message("\n Done!" ;Message("\n ===================================================================================" ;} calls=0; last_call=0; } vel[0]=0.0; vel[1]=YY; } |
3楼2016-07-08 10:00:52
4楼2018-07-18 15:01:29
