fluent UDF求解圆柱涡致振动问题（VIV.c）

基于Fluent软件平台，采用C语言编写UDF(User defined function)模拟圆柱振子在不同流速的风作用下的横向振动。采用Newmark-β方法求解结构运动方程，通过调节γ和β基本能得到较好结果。

#include “udf.h“
static real x_prev=0;
static real x_cur=0;
static real v_prev=0;
static real v_cur=0;
static real a_prev=0;
static real a_cur=0;
static double a[30000];
#define d_t 0.0005

DEFINE_EXECUTE_AT_END（execute_at_end）
{       
    real m c k;    
    real f_glob[ND_ND]m_glob[ND_ND]x_cg[ND_ND];    
    real time = RP_Get_Real（“flow-time“）;    
    m = 7.634;    
    c = 4.95;    
    k = 43398;    
    x_cg[0] = 0;    
    x_cg[1] = x_cur + v_cur*d_t;
    
    if （!Data_Valid_P（））        
    	return;
    	
    Domain *domain = Get_Domain（1）;  //return fluid domain    
    Thread *tf1 = Lookup_Thread（domain14）;    
    Compute_Force_And_Moment （domain tf1 x_cg f_glob m_glob TRUE）;    
    real lift = f_glob[1];  //get current lift force    
    
    real beta = 0.25;    
    real gamma = 0.25;    
    real p1 = 1/（gamma*d_t*d_t）;    
    real p2 = beta/（gamma*d_t）;