涡格法代码及解释_物理_自然科学_专业资料

《涡格法代码及解释_物理_自然科学_专业资料》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、#include"iostream.h"#includc"stdio.h"#include"math.h"#definePI3.1415926classAIRFOIL〃用来存放翼型的信息{public:doubleL,Bg,S;doubleXo,Xc;doubleY,Cy;AIRFOIL(){Y=0.0f,S=0.0f,L=0.0f,Bg=0.0f,Xo=0.0f,Xc=0.0f;}}；classGIRD〃网格信息{public:doublex1,z1,x2,z2;〃左右自由涡的坐标doublex3,z3,x4,z4;〃3/4弦线处的坐标doub

2、lex,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;}}；doublevec(doublex,doublez,doublex1,doublezl,doublex2,doublez2){doublea,b,c,d,c;a=l/((x2-x)*(zl-z)-(x1-x)*(z2-z));b=((x2-xl)*(xI-x)+(z2-zl)*(zl-z))/sqrt(pow((x1・x),2)+p

3、ow((z1-z),2));c=((x2-xI)*(x2-x)+(z2-z1)*(z2・z))/sqrt(pow((x2・x),2)+pow((z2・z),2));d=(l-(x1-x)/sqrt(pow((x1-x),2)+pow((z1-z),2)))/(z1-z);c=(l-(x2-x)/sqrt(pow((x2-x),2)+pow((z2-z),2)))/(z2-z);return(a*(b-c)+d-e)/4/PI;}voidGaussscidcl(intn,double*M,doubledouble*x,double*b)〃高斯■■塞

4、得尔迭带法{intt二0,i,j;〃迭代次数while(t<20)//次数限制，精度要求，此处可修改，是迭带开关{for(i=0;i

5、j];}}x[i]=(b[i]-M[i])/a[i][i];//迭代)cout«++t;for(i=0;i

6、;doubleY=0.0,M,a,ep=1e-10,p=1.22505,Cy=0.0;//p为海平面空气密度comvv“这是一个用涡格法计算机翼升力的程序!"«endl;cout«niH输入翼型个参数:展长L,根弦Bg,前缘后掠角Xo,丿匸缘后掠角XcH«endl;whilc(l){cin»airfoil.L»airfoil.Bg»airfoil.Xo»airfoil.Xc;if(airfoil.Bg-airfoil.L*(tan(airfoil.Xo*PI/180)+tan(airfoil.Xc*PI/180))/2>0){cout«airfo

7、il.L«Hu«airfoil.Bg«HH«airfoil.Xo«u“vvairfoil.Xcvv“"«endl;break;}else{cout«**翼型的稍弦为0!请重新输入翼型数据“vvendl;}}cout«H请输入来流马赫数和攻角n«endl;cin»M»a;a=a*PI/180;cout«M«,t,«a«endl;coutvv”请输入根弦上的节点数，前缘上的节点数:”vvendl;cin»Ng»Nq;cout«Ng«Mn«Nq«HH«endl;Nq・・;Ng・・;〃变成分多少块double*baseq=newdoublefNq+l]

8、;double*bascB=ncwdouble[Nq+l];double*result=newdouble[2*Nq*Ng];double*b=newdouble[2*Nq*N創;double*Ml=newdouble[2*Nq*Ng];GIRD**girdleft,**girdright;//左半边机翼，右半边机翼girdlcft=ncwGIRD*[Ng];for(i=0;i

9、]=newGIRD[Nq];}doublewidth二airfoil.L/Nq/2;〃展长每个分块的长度〃前缘节点的x坐标cout«n前