涡格法代码及解释.pdf资源-CSDN文库

版权申诉

5星 · 超过95%的资源 100 浏览量 2023-09-30 15:43:39 上传评论 1 收藏 165KB PDF 举报

资源推荐

资源详情

资源评论

#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))/sqrt(pow((x1-x),2)+pow((z1-z),2));

c=((x2-x1)*(x2-x)+(z2-z1)*(z2-z))/sqrt(pow((x2-x),2)+pow((z2-z),2));

d=(1-(x1-x)/sqrt(pow((x1-x),2)+pow((z1-z),2)))/(z1-z);

e=(1-(x2-x)/sqrt(pow((x2-x),2)+pow((z2-z),2)))/(z2-z);

return (a*(b-c)+d-e)/4/PI;

}

void Gaussseidel(int n,double *M,double **a,double *x,double *b)//高斯--塞得尔迭带法

{

int t=0,i,j;//迭代次数

while(t<20)//次数限制,精度要求,此处可修改,是迭带开关

{

for(i=0;i<n;i++)

{

M[i] = 0;

本内容试读结束，登录后可阅读更多

下载后可阅读完整内容，剩余8页未读，立即下载

内容反馈

版权申诉

睡大街的蒲牢

2023-11-06

这个资源对我启发很大，受益匪浅，学到了很多，谢谢分享~

hhappy0123456789

粉丝: 61
资源: 5万+

最新资源

资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！点击此处反馈

feedback-tip