绕圆柱无环量流动和有环量流动流线分布图_第1页
绕圆柱无环量流动和有环量流动流线分布图_第2页
绕圆柱无环量流动和有环量流动流线分布图_第3页
绕圆柱无环量流动和有环量流动流线分布图_第4页
绕圆柱无环量流动和有环量流动流线分布图_第5页
已阅读5页，还剩6页未读， 继续免费阅读

付费下载

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、绕圆柱无环量流动和有环量流动流线分布图一、绕圆柱无环量流动由均匀直线流动和偶极流的流函数可叠加得绕圆柱无环量流动的流函数方程为：=V r sin()-M sin()2r 相应流动的流线方程为：V r sin()- M sin()2r =C（常量）不同C值对应不同的流函数曲线，将极坐标形式变成直角坐标形式：V y-M2ryx2+y2 =C（常量）1. V=10，M=500, C=10,-10,0,20,-20 ezplot(10*y-500*y/(2*pi*(x2+y2)=10) hold on ezplot(10*y-500*y/(2*pi*(x2+y2)=-10) hold on ezplo

2、t(10*y-500*y/(2*pi*(x2+y2)=0) hold on ezplot(10*y-500*y/(2*pi*(x2+y2)=20) hold on ezplot(10*y-500*y/(2*pi*(x2+y2)=-20) hold on2. V=20，M=520, C=20,-20,0,40,-40 ezplot(20*y-520*y/(2*pi*(x2+y2)=20) hold on ezplot(20*y-520*y/(2*pi*(x2+y2)=-20) hold on ezplot(20*y-520*y/(2*pi*(x2+y2)=0) hold on ezplot(20

3、*y-520*y/(2*pi*(x2+y2)=40) hold on ezplot(20*y-520*y/(2*pi*(x2+y2)=-40) hold on3. V=30，M=540, C=30,-30,0,60,-60 ezplot(30*y-540*y/(2*pi*(x2+y2)=30) hold on ezplot(30*y-540*y/(2*pi*(x2+y2)=-30) hold on ezplot(30*y-540*y/(2*pi*(x2+y2)=0) hold on ezplot(30*y-540*y/(2*pi*(x2+y2)=60) hold on ezplot(30*y-

4、540*y/(2*pi*(x2+y2)=-60) hold on4. V=100，M=1000, C=30,-30,0,60,-60 ezplot(100*y-1000*y/(2*pi*(x2+y2)=30) hold on ezplot(100*y-1000*y/(2*pi*(x2+y2)=-30) hold on ezplot(100*y-1000*y/(2*pi*(x2+y2)=0) hold on ezplot(100*y-1000*y/(2*pi*(x2+y2)=60) hold on ezplot(100*y-1000*y/(2*pi*(x2+y2)=-60) hold on5.

5、V=10，M=100, C=30,-30,0,60,-60ezplot(10*y-100*y/(2*pi*(x2+y2)=30) hold on ezplot(10*y-100*y/(2*pi*(x2+y2)=-30) hold on ezplot(10*y-100*y/(2*pi*(x2+y2)=0) hold on ezplot(10*y-100*y/(2*pi*(x2+y2)=60) hold on ezplot(10*y-100*y/(2*pi*(x2+y2)=-60) hold on绕圆柱有环量流动由均匀直线流动、偶极流、点涡三者叠加后的复合流动，流函数方程为：=V（r-R2r）si

6、n()+2 ln(r)相应流动的流线方程为：V（r-R2r）sin()+2 ln(r)=C（常量）不同C值对应不同的流函数曲线，将极坐标形式变成直角坐标形式：V (x2+y2-R2x2+y2) +2 lnx2+y2 =C（常量）1.V=10, R=3, =0.1, C=0,10,20,30,40ezplot(10*y*(1-9/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=0) hold on ezplot(10*y*(1-9/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=10) hold on ezplot(10*y*(1-9/(x2+y2)+

7、0.1*exp(sqrt(x2+y2)/(2*pi)=20) hold onezplot(10*y*(1-9/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=30) hold onezplot(10*y*(1-9/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=40) hold on2.V=20, R=2, =0.1, C=0,20,-20,40,-40ezplot(20*y*(1-4/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=0) hold on ezplot(20*y*(1-4/(x2+y2)+0.1*exp(sqrt

8、(x2+y2)/(2*pi)=20) hold on ezplot(20*y*(1-4/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=-20) hold onezplot(20*y*(1-4/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=40) hold onezplot(20*y*(1-4/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=-40) hold on3.V=50, R=4, =0.1, C=0,10,-10,20,-20ezplot(50*y*(1-16/(x2+y2)+0.1*exp(sqrt(x2+y2)/

9、(2*pi)=0) hold on ezplot(50*y*(1-16/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=10) hold on ezplot(50*y*(1-16/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=-10) hold onezplot(50*y*(1-16/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=20) hold onezplot(50*y*(1-16/(x2+y2)+0.1*exp(sqrt(x2+y2)/(2*pi)=-20) hold on4.V=50, R=3, =-0.3, C=0

10、,20,-20ezplot(50*y*(1-9/(x2+y2)-0.3*exp(sqrt(x2+y2)/(2*pi)=0,-6 6,-5 5) hold on ezplot(50*y*(1-9/(x2+y2)-0.3*exp(sqrt(x2+y2)/(2*pi)=-20,-6 6,-5 5) hold on ezplot(50*y*(1-9/(x2+y2)-0.3*exp(sqrt(x2+y2)/(2*pi)=20,-6 6,-5 5) hold on5.V=1000, R=5, =-1, C=0,20,-20,40，-40ezplot(1000*y*(1-25/(x2+y2)-1*exp(sqrt(x2+y2)/(2*pi)=0,-6 6,-5 5) hold on ezplot(1000*y*(1-25/(x2+y2)-1*exp(sqrt(x2+y2)/(2*pi)=-20,-6 6,-5 5) hold on ezplot(1000*y*(1-25/(x2+y2)-1*exp(sqrt(x2+y2)/(2*pi)=20,-6

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论