涡格法代码及解释解读.docVIP

#include iostream.h #include stdio.h #include math.h #define PI 3.1415926 class AIRFOIL //用来存放翼型的信息 { public: double L,Bg,S; double Xo,Xc; double Y,Cy; AIRFOIL(){Y=0.0f,S=0.0f,L=0.0f,Bg=0.0f,Xo=0.0f,Xc=0.0f;} }; class GIRD //网格信息 { public: double x1,z1,x2,z2;//左右自由涡的坐标 double x3,z3,x4,z4;//3/4弦线处的坐标 double x,z;//控制点的坐标,3/4弦线中点 GIRD(){x1=0.0f,x2=0.0f,z1=0.0f,z2=0.0f,x3=0.0f,x4=0.0f,z3=0.0f,z4=0.0f,x=0.0f,z=0.0f;} }; double vec(double x,double z,double x1,double z1,double x2,double z2 ) { double a,b,c,d,e; a=1/((x2-x)*(z1-z)-(x1-x)*(z2-z)); b=((x2-x1)*(x1-x)+(z2-z1)*(z1-z))/