空气动力学模板.docxVIP

下载本文档

10
0
约2.25千字
约 5页
2020-01-11 发布于上海
举报
版权申诉

空气动力学模板.docx

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

首先根据公式，利用MATLAB计算出所需的量 % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % y1=@(x)atan((18*cot(x).*sin(x).*sin(x)-2*cot(x))./(9*cos(2*x)+14.6))-pi/18; x1=fzero(y1,0.5) x1 = 0.4779 hold on y2=@(x)atan((18*cot(x).*sin(x).*sin(x)-2*cot(x))./(9*cos(2*x)+14.6))-pi/12; x2=fzero(y2,0.5) x2 = 0.5627 Mn1=3*sin(x1) Mn1 = 1.3798 Mn2=3*sin(x2) Mn2 = 1.6004 M1=Mn1/sin(x1-pi/18) M1 = 4.6185 M2=Mn2/sin(x2-pi/12) M2 = 5.3999 % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % 然后根据③④两处流动的方向相同，能够根据几何关系得到；这一关系可以转化为有关β的方程。 ③④两处的压强相同，且根据公式,建立之间的关系。最终可以得到以下方程，，由此可以得到二元方程组。首先画出其函数图像，程序如下； % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % z=ezplot(306.7652*sin(x).^2=574.0662*sin(y).^2-0.758,[0,2*pi/3]) hold on z=ezplot((42.6611*cot(x).*sin(x).*sin(x)-2*cot(x))./(31.8628+21.3305*cos(2*x))=(19.96846+13.597*cos(2*y)+2*cot(y)-58.3178*sin(y).*sin(y).*cot(y))./(42.8225+29.1589*cos(2*y)-0.9326*cot(y)+27.194*sin(y).*sin(y).*cot(y)),[0,2*pi/3]); set(gca,XTick,0:0.1:3*pi/2) set(gca,YTick,0:0.1:3*pi/2) grid title(x与y之间的关系 ) % x=β3；y=β4 % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % 会得到如下函数图像从图像中可以看出方程的解的范围是0.4x0.5;0.3y0 之后利用MATLAB解该方程组，如下所示； % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % function F = myfun(x) F(1)=306.7652*sin(x(1)).^2-574.0662*sin(x(2)).^2+0.758; F(2)=(42.6611*cot(x(1)).*sin(x(1)).*sin(x(1))-2*cot(x(1)))./(31.8628+21.3305*cos(2*x(1)))-(19.96846+13.597*cos(2*x(2))+2*cot(x(2))-58.3178*sin(x(2)).*sin(x(2)).*cot(x(2)))./(42.8225+29.1589*cos(2*x(2))-0.9326*cot(x(2))+27.194*sin(x(2)).*sin(x(2)).*cot(x(2))); % 定义方程组F % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % 之后在命令行窗口输入； fsolve(@myfun,[0.4 0.3],optimset(fsolve)) %根据前面的判断，取初始迭代值x0=[0.4 0.3] 会得到如下的结果 ans = 0.4431 0.3210 即;；在计算相应的θ值 % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % % x1=0.4431; y1=atan((42.6611*cot(x1).*sin(x1).*sin(x1)-2*cot(x1))./(31.8628+21.3305*cos(2*x1))) y1 = 0.2650 x2=0.3210; y2=atan((58.3178*cot(x2).*sin