UDF_fluent造波UDF_边界造波_fluent造波_fluent波_udf_

#include "udf.h" static real v_prev = 0.0; /*定义阀体速度*/ static real loc_prev=0.0; /*定义弹簧变形量*/ #define k 200000 /*定义弹簧弹性系数*/ #define m 0.005 /*定义阀体质量*/ #define Lk 0.003 /*定义弹簧行程*/ DEFINE_CG_MOTION(value, dt, vel, omega, time, dtime) { Thread *t; face_t f; real NV_VEC (A); /*定义流体作用力，合外力，弹簧弹力，加速度*/ real f1,f2,f_s, dv; NV_S (vel, =, 0.0); NV_S (omega, =, 0.0); if (!Data_Valid_P ()) return; t = DT_THREAD (dt); /*求解流体作用力*/ f1= 0.0; begin_f_loop (f, t) { F_AREA (A, f, t); f1 = f1 + F_P (f, t) * A[0]; } end_f_loop (f, t) /*求解弹簧弹力*/ f_s=k*loc_prev; /*求解合外力*/ f2=f1-f_s; /*求解加速度*/ dv = dtime * f2 / m; /*判断弹簧是否运动至上止点*/ /*如果弹簧位移大于总行程且合外力仍然大于0，则处于上止点,阀体速度为0*/ if ((loc_prev>Lk)&&(f2>0)) { v_prev=0; } /*否则阀体继续运动*/ else { v_prev = v_prev + dv; loc_prev=loc_prev+v_prev*dtime; } /*显示一些变量*/ Message("\n\ntime = %f, x_vel = %f, f2 = %f, f_s = %f, loc_prev=%f, dv=%f\n", time, v_prev, f2, f_s ,loc_prev, dv); vel[0] = v_prev; } DEFINE_CG_MOTION(value_1, dt, vel, omega, time, dtime) { vel[0] = v_prev; }