[ÇóÖú] ¡¾fluent ²¢ÐÐ UDF ÇóÖú¡¿ÎªÉ¶fluent ²¢ÐÐÕâôÂý£¿

#include <stdio.h>
#include "udf.h"
#include "math.h"
#include "mem.h"
#include "dynamesh_tools.h"
#define pi 3.141592654
#define rou 998.2
#define cord 0.2
#define h0 0.12
#define st 0.3
#define u 1.0
#define n_points 7 // number of pressure sample points
#define foil_id 3 // fluent_id of the foil
double theta0 = 28.0*pi/180.0;
double w = pi*st*u/h0;
double t = 2*h0/(st*u);
double eps = -100*pi/180 ; // -180 ~ 180
file *fc_up, *fp_up; // upstream
file *fc_dn, *fp_dn; // downstream
file *prs_foil_up, *prs_foil_dn, *fpresult; // pressure output
/**********************************************************************************************/
define_cg_motion(foil_upstream,dt,vel,omega,time,dtime)
{
real pressure = 0.0;
real f1, f3, f5;
real pl, pt, p;
real x[nd_nd] = {0.0};
real f_cg[nd_nd], m_cg[nd_nd], x_cg[nd_nd];
int n;
int i = 0;
int j = 0;
//-------------------------------------------------------------------------
#if !rp_host
face_t face ;
cell_t cell ;
domain *domain = get_domain(1);
thread *thread = dt_thread(dt);
//-------------------------------------------------------------------------
vel[0] = 0.0;
vel[1] = -h0*(cos(w*(time + dtime)) - cos(w*time)) / dtime;
omega[0] = 0.0;
omega[1] = 0.0;
omega[2] = -theta0*(sin(w*(time + dtime)) - sin(w*time)) / dtime;
//-------------------------------------------------------------------------
if (n_time % 20 == 1)
{
prs_foil_up = fopen("pressure_up.dat", "a+");
fprintf(prs_foil_up, " %.4f
", time*w / (2 * pi));
begin_f_loop(face, thread)
{
   f_centroid(x, face, thread);
   pressure = f_p(face, thread) + rp_get_float("operating-pressure"); // + absolute pressure
   fprintf(prs_foil_up, " %d %.4f %.4f %.2f
", i + 1, x[0], x[1], pressure);
    i++; }
end_f_loop(face, thread)
fclose(prs_foil_up);
}
#endif
//-------------------------------------------------------------------------
#if !rp_host if(time>8*dtime)
{
// f1 f3 f5 output
x_cg[0] = 0.0;
x_cg[1] = -h0*cos(w*time);
compute_force_and_moment (domain, thread, x_cg, f_cg, m_cg, false);
f1 = f_cg[0]/(0.5*rou*cord*u*u);
f3 = f_cg[1]/(0.5*rou*cord*u*u);
f5 = m_cg[2]/(0.25*rou*cord*cord*u*u);
pl = f_cg[1] * vel[1];
pt = m_cg[2] * omega[2];
p = pl + pt;
fc_up = fopen("force_upstream.dat", "a+");
fp_up = fopen("power_upstream.dat", "a+");
fprintf(fc_up, "%.3f %10.3f %10.3f %10.3f
", time*w / (2 * pi), f1, f3, f5);
fprintf(fp_up, "%.3f %10.3f %10.3f %10.3f
", time*w / (2 * pi), pl, pt, p );
fclose(fc_up);
fclose(fp_up);
}
#endif
return;
}