船舶流体力学第7章打印.docx
《船舶流体力学第7章打印.docx》由会员分享,可在线阅读,更多相关《船舶流体力学第7章打印.docx(20页珍藏版)》请在冰豆网上搜索。
船舶流体力学第7章打印
第七章势流理论
(二)
本章主要讨论:
轴对称有势流动和机翼绕流的有关理论
§7.1轴对称流动
一条曲线绕轴旋转一周形成的物体形状称为旋成体。
当来流沿旋成体中轴线方向绕流旋成体时,通过中轴线的各子午面上的流动均相同,这种流动称为
轴对称流动。
比如,均匀流绕圆球的流动。
r
―►V
X
—►►
轴对称轴
对于无旋轴对称流动,存在速度势函数0和流函数。
但,速度势函数0是调和函数,流函数不是调和函数。
采用柱坐标(r,,x),设x轴为对称轴,流动参数不随变化。
vrVr(r,x,t)VxVx(r,x,t)
比如:
不可压缩流体的轴对称势流应该满足:
求解不可压缩流体轴对称势流问题的主要任务就是寻求满足以上方程组和边界条件的速度矢量。
有两种数学求解途经:
途径一:
控制方程:
物面无穿透条件:
无穷远处来流:
这里:
Vr,Vx
rx
速度势函数0是调和函数,可以采用叠加法求解。
途径二:
控制方程:
D2
物面无穿透条件:
Vn
0无穷远处来流:
S
V
—*■
V
这里:
vr1,
rx
1
Vx
rr
2x
2r
rr
流函数函数屮不是调和函数,称为斯托克斯函数。
但它是线性的,也可采用叠加法求解。
••基本的轴对称势流:
1.均匀直线流:
>V
x
轴对称轴
Vr0,VxV
V0
Vr0,
VxV
Vx
r
x
又Vr0,
1
Vx
V
1—Vr
rx
rr
2
2.空间点源(汇)流:
(0,0)处有一点源Q:
Q4
R2vr
又:
如图,有:
4R24r2x2
Vr
—vrvRsin
r
Qr
4r2x2—厂x2
-Vx
x
vRcos
rvx
rxQ
4r2
Q
4r2x2
且:
rvr
r2Q
Qx
4.r2x2
即:
当点源在X。
点(轴对称轴上),速度势函数和流函数为:
1
r
xo
x
Q
1
Q
IP
xxo
4F
2
xX。
4
Jr
22
xX。
3.空间偶极子流:
令:
limQxM0
x0
Q
令:
1
1
lim
x0Q
Qx
22.rx
x
2x
M
1
M
x
4
x
4
x、r2x2
4
3
22?
rx
亦可得:
当偶极子在X0点(轴对称轴上),速度势函数和流函数为:
XXo
3
22
rxXo
M
2r
3
4
2
2t
r
xXo
二•均匀来流绕圆球体的流动:
xRcos
采用球坐标(R入)。
柱坐标与球坐标的关系为:
rRsin
均匀流:
偶极子流:
叠加后得到
Vx1Vr2
2
求出速度:
VRcos
2cos
4R2
」VR2sin2
2
M—sin
M
1
2
3COS
R3
M1.
3sin
4R3
在球表面Vr=0,故:
M2R3V
Ro
3匹
\2V
2RCOS0
"*VR寻Si『0
相应地:
R)3
VrV1COS0
R3
球表面速度分布:
VR0V0
设无穷远处压强为
P,由伯努利方程,有:
veV1
R03斎
sin0
3、,
V
sin0
2
V2
V2
VR2v2
P
P
Pc
2
2
2
于是,得到球表面的压强分布:
2
V彳9•2A
pp1sin0
24
球表面的压强系数分布:
CP1討20
例1:
x=d,点汇
-Q;x=-d,点源
流体作用在球体上的阻力和升力均为零。
解:
叠加三个基本势流的流函数,得到:
x
r2
令=o,得到零流线方程
-Vr2—
24
r2xd2
代数方程给出了两条曲线,一条是与轴重合的直线,另一条是卵形封闭曲线。
显然,流函数屮=C.给出了均匀直线流绕流卵形回转体所形成的势流流场的流线。
这类卵形回转体也称为兰金(Rankine)体。
§7.3有限翼展机翼
对机翼理论的研究是流体力学中最引人注目的应用课题之一。
舰船上的舵、水翼、减摇鳍等本身就是机翼,螺旋桨、透平机械的叶片、水泵的叶片等都是利用机翼的原理工作的。
我们可以利用机翼原理来产生升力(例如飞机、风筝等)或推力(例如螺旋桨等),因此机翼理论
的研究对船舶工程有重要意义。
.机翼的几何参数:
翼型:
翼型是机翼剖面的基本形状。
翼型具有产生的升力与阻力之比(升阻比」可能大的体形,整体上是优良流线形,使流体能顺着其表面尽可能无分离地向尖后缘流去。
如图所示为翼型无分离地绕流。
前缘或导边(leadingedge):
迎流的一端。
后缘或随边(trailingedge):
攻角a(angleofattack):
来流与弦之间的夹角。
竝音連tPLtTLM
工程实际中应用的一些翼型的基本形状:
舵u江寸林a
后缘总是尖的(产生环量)圆前缘:
减小形状阻力。
尖前缘:
减小压缩性所引起的激波阻力或自由表面所引起的兴波阻力。
中线(centerline):
翼型内各圆弧中点的连线。
翼弦(chord):
中线两端的连线,常作为翼型基线。
对称翼型:
中线与弦线重合的翼型。
厚度(thicheness)t:
翼弦的垂线与翼型上下表面交点之间的最大距离。
相对厚度:
翼厚与弦长之比。
二、机翼的平面图形
机翼的常见平面图形:
般来说,翼型的厚度与翼弦相比要小得多,许多实用场合中翼展比翼弦大得多。
展弦比:
入=翼展的平方/翼面积S
入无限翼展机翼,即为二元机翼。
二.有限翼展机翼:
实际上机翼的展弦比均为有限值,故流动是三维的。
对于无限翼展机翼,可近似用一根无限长的涡线(涡线有r)来代替,称附着涡。
而对于有限翼展机翼,去卩不能用有限长附着涡来代替机翼,因为这样旋涡会在流体内终止。
对于有限翼展机翼,由于下翼面压力大于上翼面:
卜、、、$
乙於八、\、、1
<二-一
*4-+亠
上翼面流线向中间偏移,下翼面流线相反。
上下压差作用下产生自由涡。
—Hr
*+44丿
上翼面
下翼面
上
下
自由涡与附着涡联成n形涡。
由海姆霍兹定理已知n形涡r=常数。
图片
三.下洗和诱导阻力:
如图,对于矩形机翼上任一点A,坐标为y,用半无穷直线涡公式得左自由涡在该点所诱导的速度:
左自由涡产生的沿翼展的平均诱导速度为:
因左右对称,整个机翼下面的平均诱导速度为:
Vz
wi
Wi2elVzdy
1e
2e1
W—
l
el
—3—4
2ley2le
左、右翼端涡在机翼下面产生的平均诱导速度,
方向向下,称为下洗速度,或称为下滑速度。
来流速度与下洗速两速度矢相加:
Vw式中V为实际(有效)来流速度。
i
e
式中a为有效攻角,a为下洗角或下滑角。
V的方向与翼弦的夹角为:
下洗角可由下式计算:
i
tan1
Wi
V
因为Wi
向下故为负值。
库塔一儒柯夫斯基力为:
LV
力L'在升力和阻力方向的投影分别为:
L
Vcosi
RiVsini
一般地,下洗速度Wi很小,
即ai很小,
故有:
isin
itani
这时:
r-Wi
下洗角:
ii,
升力:
L
V,
诱导阻力:
Ri
Li
wi
如果在翼端装上当板,限制绕流,可减小诱导阻力,如图所示:
§7.4升力线理论
一.有限翼展机翼的升力模型:
实际有限翼展机翼沿翼展方向的剖面的形状,安装角度有变化,各个截面环量也变化。
如图,用n形涡系代替单一的n形涡,附着涡在翼展上迭合在一起形成升力线,n形涡系的自由涡连成一整体而形成涡面。
虽然每根n形涡环量不变,但沿翼展不同截面有数目不同的n形涡,所以沿翼展环量是变化的。
二.有限翼展机翼的升力线理论:
入>2:
大展弦比机翼。
入V2:
小展弦比机翼或短翼。
入>2时机翼的附着涡系可用一根涡丝来代替,这根涡丝通常称为升力线(liftline)。
升力线理论:
以升力线为理想模型的计算机翼动力特性的理论。
弓I入两点假定:
(1)自由涡面是平面,延伸至无穷远而不翻卷成两股大涡,自由涡面旋涡角速度矢量平行来流。
(2)翼面上横向流动很小,任一剖面处可作平面流动处理,三元效应仅考虑各翼剖面处下洗速度和
下洗角的不同。
这就是“简单的切片理论”方法:
处强度为d^d()d的涡丝在升力线上y点产生的下洗速度为:
d
dw1()d
UVVj
4y
沿展向积分得整个自由涡在y处的诱导速度:
1;2()d
\A/
VV|
4iy
对于小攻角,下洗角ai为小量,有:
i
wi
V
宽度为dy的一段机翼的二维升力为:
dL
V
(y)dy
按定义升力垂直于来流:
dLdL
COSi
V
(y)dy
诱导阻力:
dRdLtanjw}
(y)dy
整个机翼的升力和诱导阻力:
!
2
l2
12J2zx|
LV(y)dyRi
Wi(y)
(y)dy
—(y)』一dy
l2
412
41212y
.升力系数和诱导阻力系数:
LV(y)dy
而:
(y)2lVAnsinn式中A为待定常数。
n1
22
LVlAnsinnsind
n10
升力系数:
CLL一-A1A
1V2Ss
2
四.具有最小诱导阻力的机翼平面形状—椭圆机翼:
显然,当6=0时,阻力最小。
对应的机翼环量分布为:
()2VlA1sin
即:
0sin或:
—sin(a)
0
其中:
0—2VlA1且:
Icos(b)
2l2
22
(a),(b)两式两边平方后相加得:
一1
0少2
即,最小诱导阻力系数的机翼的环量分布为椭圆形状。
相应的下洗角为:
Wi
V
Ai
诱导阻力系数为:
CR
显然,
Cl
时,Cr0。
即,无限翼展机翼没有诱导阻力。
对于其它无扭转的非椭圆机翼,其下洗角和阻力系数修正为:
Cl
CRi
ci
(i
实际中常采用梯形
机翼。
五.展弦比换算:
在进行机翼设计,比如船用舵的设计时,常采用展弦比换算方法。
设两机翼平面形状,翼型及弦长都相同,例如矩形机翼1、2,展弦比分别为入1和入2
下洗角沿翼展的分布为:
i(y)
沿翼展下洗角的平均值:
或:
1
「°
k
nAn
n1sin
sinnI
sin
i(y)dy
k
A2n1
n1
所以:
i^1
k
A2n1
n2
A
展弦比换算步骤如下:
el
e2
Cl
⑴
1
ii
1
1
Cl
⑵
由相似原理知:
兔1=
2
在其上任取一点E,
Cl
1
Cl
所对应的升力系数为C
设翼1的C
l〜a'曲线已知,
l,求出几何攻角之差:
CM
I1
重复上面步骤得一系列翼2上的点,连接它便是
2曲线。
1
例1:
一飞机自重21582N,机翼面积为20m-翼展11m,若水平方向飞行速度为280km/h,流体密
度卩=1.226kg/m3。
求:
1)升力系数,展弦比,环量;
2)设机翼平面形状为矩形,求诱导阻力系数。
解:
展弦比:
l26.05。
因飞行水平升力与飞机自重平衡,则升力系数:
Cl
0.29
环量为:
-20.57m2/s
Vl
查表7.4.1(参见P133)得机翼平面形状为矩形时:
-
(1)
0.335。
诱导阻力系数:
CRi
0.0047。
例2:
一机翼弦长2m,展长10m,两端为零,环量沿翼展呈椭园型分布。
以360km/h的速度在大气中飞行,
设机翼中部的环量r0=20m2/s,
求:
升力系数及诱导阻力系数(p
=1.2kg/m3)。
解:
环量分布为:
即:
220.52y2
10
4.52y2
升力:
(y)dy
5
V452y2dy18.85(KN)
5
升力系数:
Cl
又:
故,诱导阻力系数:
0.157
V2S
102
210
C2
CRL0.00157。