[ÇóÖú] Fluent ÈçºÎ»ñÈ¡ »úÒí ±íÃæѹÁ¦·Ö²¼£¿

DEFINE_CG_MOTION(foil,dt,vel,omega,time,dtime)
{
   real x_A=0.0,y_A=0.0,x_N=0.0,y_N=0.0;
   real theta_A,theta_N;
   real omegap = 0.0;
/////////////////////////////////////////////////////////
    face_t f;
    cell_t c;
        real f_glob[ND_ND],m_glob[ND_ND],x_cg[ND_ND];
    real f1,f3,f5;
    Domain *domain= Get_Domain (1);
    Thread *tf1 = Lookup_Thread (domain, 38);
/////////////////////////////////////////////////////////
   theta_A = -theta0*cos(w*CURRENT_TIME);
   theta_N = -theta0*cos(w*(CURRENT_TIME+CURRENT_TIMESTEP));
   omegap = (theta_N-theta_A)/CURRENT_TIMESTEP;

   x_A= 0.0;
   x_N= 0.0;
   y_A=H0*sin(w*CURRENT_TIME);
   y_N=H0*sin(w*(CURRENT_TIME+CURRENT_TIMESTEP));
  
   vel[0] =0.0;
   vel[1] =(y_N - y_A) / CURRENT_TIMESTEP;
   omega[0] =0.0;
   omega[1] =0.0;
   omega[2] = omegap;
   // f1  f3  f5 output //////////////////////////////////
    x_cg[0]=DT_CG(dt)[0];
    x_cg[1]=DT_CG(dt)[1];

    Compute_Force_And_Moment (domain, tf1, x_cg, f_glob, m_glob,1);
   
    f1 = f_glob[0]/(0.5*998.1*H0*U*U);
    f3 = f_glob[1]/(0.5*998.1*H0*U*U);
    f5 = m_glob[2]/(0.25*998.1*H0*H0*U*U);

    fp=fopen("force.dat", "a+");
    fp1=fopen("point.dat", "a+");
    fprintf(fp,"%.6f %.5f %.5f %.5f\n",CURRENT_TIME*w/(2*PI),f1,f3,f5);
    fprintf(fp1,"%.6f %.6f\n",x_cg[0],x_cg[1]);
        fclose(fp);
        fclose(fp1);
   //////////////////// //////////////////////////////////
  return;
}

3Â¥: Originally posted by xuwenhua625 at 2016-03-18 16:55:37
ÎÒÖªµÀ Ôõô×öÁË

   // presure ///////////////////////////////////////////
   pre=fopen("presure.dat", "a+";
     begin_f_loop(f, tf1 )
    {
        F_CENTROID(x,f,tf1 );
...

4Â¥: Originally posted by hitzhx at 2016-05-17 10:22:15
Â¥Ö÷д¹ýÒíÐͱíÃæ»ý±ùµÄ¶¯Íø¸ñô£¬Ñ§Ï°Ò»Ï¡£¡£...