資源描述：
《繞圓柱無環(huán)量流動和有環(huán)量流動流線分布圖.doc》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、繞圓柱無環(huán)量流動和有環(huán)量流動流線分布圖一、繞圓柱無環(huán)量流動由均勻直線流動和偶極流的流函數(shù)可疊加得繞圓柱無環(huán)量流動的流函數(shù)方程為：Ψ=V∞r(nóng)sin(θ)-Msin(θ)2πr相應(yīng)流動的流線方程為：V∞r(nóng)sin(θ)-Msin(θ)2πr=C（常量）不同C值對應(yīng)不同的流函數(shù)曲線，將極坐標(biāo)形式變成直角坐標(biāo)形式：V∞y-M2πryx2+y2=C（常量）1.V∞=10，M=500,C=10,-10,0,20,-20ezplot('10*y-500*y/(2*pi*(x^2+y^2))=10')holdonezplot('10*y-500*y/(2*pi*(x^2+y^2))=-10')h
2、oldonezplot('10*y-500*y/(2*pi*(x^2+y^2))=0')holdonezplot('10*y-500*y/(2*pi*(x^2+y^2))=20')holdonezplot('10*y-500*y/(2*pi*(x^2+y^2))=-20')holdon2.V∞=20，M=520,C=20,-20,0,40,-40ezplot('20*y-520*y/(2*pi*(x^2+y^2))=20')holdonezplot('20*y-520*y/(2*pi*(x^2+y^2))=-20')holdonezplot('20*y-520*y/(2*pi*
3、(x^2+y^2))=0')holdonezplot('20*y-520*y/(2*pi*(x^2+y^2))=40')holdonezplot('20*y-520*y/(2*pi*(x^2+y^2))=-40')holdon3.V∞=30，M=540,C=30,-30,0,60,-60ezplot('30*y-540*y/(2*pi*(x^2+y^2))=30')holdonezplot('30*y-540*y/(2*pi*(x^2+y^2))=-30')holdonezplot('30*y-540*y/(2*pi*(x^2+y^2))=0')holdonezplot('30
4、*y-540*y/(2*pi*(x^2+y^2))=60')holdonezplot('30*y-540*y/(2*pi*(x^2+y^2))=-60')holdon4.V∞=100，M=1000,C=30,-30,0,60,-60ezplot('100*y-1000*y/(2*pi*(x^2+y^2))=30')holdonezplot('100*y-1000*y/(2*pi*(x^2+y^2))=-30')holdonezplot('100*y-1000*y/(2*pi*(x^2+y^2))=0')holdonezplot('100*y-1000*y/(2*pi*(x^2+
5、y^2))=60')holdonezplot('100*y-1000*y/(2*pi*(x^2+y^2))=-60')holdon5.V∞=10，M=100,C=30,-30,0,60,-60ezplot('10*y-100*y/(2*pi*(x^2+y^2))=30')holdonezplot('10*y-100*y/(2*pi*(x^2+y^2))=-30')holdonezplot('10*y-100*y/(2*pi*(x^2+y^2))=0')holdonezplot('10*y-100*y/(2*pi*(x^2+y^2))=60')holdonezplot('10*y
6、-100*y/(2*pi*(x^2+y^2))=-60')holdon二、繞圓柱有環(huán)量流動由均勻直線流動、偶極流、點(diǎn)渦三者疊加后的復(fù)合流動，流函數(shù)方程為：Ψ=V∞（r-R2r）sin(θ)+Γ2πl(wèi)n(r)相應(yīng)流動的流線方程為：V∞（r-R2r）sin(θ)+Γ2πl(wèi)n(r)=C（常量）不同C值對應(yīng)不同的流函數(shù)曲線，將極坐標(biāo)形式變成直角坐標(biāo)形式：V∞(x2+y2-R2x2+y2)+Γ2πl(wèi)nx2+y2=C（常量）1.V∞=10,R=3,Γ=0.1,C=0,10,20,30,40ezplot('10*y*(1-9/(x^2+y^2))+0.1*exp(sqrt(x^2+y^2))
7、/(2*pi)=0')holdonezplot('10*y*(1-9/(x^2+y^2))+0.1*exp(sqrt(x^2+y^2))/(2*pi)=10')holdonezplot('10*y*(1-9/(x^2+y^2))+0.1*exp(sqrt(x^2+y^2))/(2*pi)=20')holdonezplot('10*y*(1-9/(x^2+y^2))+0.1*exp(sqrt(x^2+y^2))/(2*pi)=30')holdonezplot('10*y*(1-9/(x^2+y
