ROCKET PROPULSION ELEMENTS——DEFINITIONS AND FUNDAMENTALS笔记
火箭推进理论包括力学、热力学、化学
单位换算
磅力
l
b
f
lbf
lbf
磅,磅质量
l
b
m
lbm
lbm,
l
b
lb
lb,
p
o
u
n
d
s
pounds
pounds,
l
b
s
lbs
lbs
1
s
l
u
g
=
1
l
b
f
1
f
t
/
s
2
1 slug = \frac{1 lbf}{1 ft/s^2}
1slug=1ft/s21lbf
1
k
g
=
1
N
1
m
/
s
2
1 kg = \frac{1 N}{1 m/s^2}
1kg=1m/s21N
1
l
b
f
=
4.45
N
1 lbf = 4.45 N
1lbf=4.45N
1
f
t
=
30.48
c
m
1 ft = 30.48 cm
1ft=30.48cm
1
s
l
u
g
=
14.593904
k
g
1 slug = 14.593904kg
1slug=14.593904kg
1
c
a
l
=
4.184
J
1 cal = 4.184 J
1cal=4.184J
1
e
r
g
=
1
0
−
7
J
1 erg = 10^{-7} J
1erg=10−7J
1
e
V
=
1.60218
×
1
0
−
19
J
1 eV = 1.60218\times 10^{-19} J
1eV=1.60218×10−19J
1
B
t
u
=
1055.05585
J
1 Btu = 1055.05585 J
1Btu=1055.05585J
1
l
b
=
0.454
k
g
1 lb = 0.454kg
1lb=0.454kg
1
k
g
=
2.20462262
l
b
1kg = 2.20462262 lb
1kg=2.20462262lb
定义
总冲(量) I t = ∫ 0 t F d t I_t=\int_{0}^{t} F dt It=∫0tFdt
质量比冲
I
s
=
∫
0
t
F
d
t
∫
m
˙
d
t
I_s = \frac{\int_{0}^{t} F dt}{\int \dot{m} dt}
Is=∫m˙dt∫0tFdt,单位
m
/
s
m/s
m/s
重量比冲
I
s
=
∫
0
t
F
d
t
∫
m
˙
g
d
t
I_s = \frac{\int_{0}^{t} F dt}{\int \dot{m} g dt}
Is=∫m˙gdt∫0tFdt,单位
s
s
s
书中的比冲指重量比冲,对于固发,重量比冲往往根据总冲 与 发动机初始与最终重量之差 来计算
推进剂总流量
m
˙
\dot{m}
m˙
等效推进剂总重量
w
=
I
t
/
I
s
w = I_t / I_s
w=It/Is
推进剂重量流量
w
˙
=
m
˙
g
\dot{w} = \dot{m} g
w˙=m˙g
体积流量是指流经管道横截面的流体的体积与该体积通过截面所用的时间之比。$m^3/t$
质量流量是指流经管道横截面的流体的质量与该体积通过截面所用的时间之比。$kg/s$
等效排气速度 c = F / m ˙ c = F/\dot{m} c=F/m˙
推进剂比耗量为比冲倒数
质量比为最终质量与初始质量之比
M
R
=
m
f
/
m
o
MR = m_f / m_o
MR=mf/mo
总的质量比为各级质量比之积
举例:三级火箭质量分别是
m
0
,
m
1
,
m
2
m_0, m_1, m_2
m0,m1,m2
各级质量比为:
(
m
1
+
m
2
)
/
(
m
0
+
m
1
+
m
2
)
(m_1+m_2)/(m_0+m_1+m_2)
(m1+m2)/(m0+m1+m2)、
m
2
/
(
m
1
+
m
2
)
m_2/(m_1+m_2)
m2/(m1+m2)
总质量比为:
m
2
/
(
m
0
+
m
1
+
m
2
)
m_2/(m_0+m_1+m_2)
m2/(m0+m1+m2)
推进剂质量分数 为推进剂质量与初始质量之比
ζ
=
(
m
o
−
m
f
)
/
m
0
=
m
p
/
m
o
\zeta = (m_o - m_f)/m_0 = m_p / m_o
ζ=(mo−mf)/m0=mp/mo
推进剂质量
m
p
m_p
mp
推进系统死重包括:贮存推进剂、燃烧推进剂所需的硬件
冲重比为总冲与初始重量之比
推重比为推力与初始重量之比
初始重量
w
0
=
(
m
f
+
m
p
)
g
w_0 = (m_f+m_p)g
w0=(mf+mp)g
补充
重量(weight):由于引力(Gravity )而作用到物体上的力
W
=
m
g
W=mg
W=mg,
g
g
g为重力(Gravity )加速度
测量 measurements
标准 standard
国际单位制 International System of Units (SI)
公制 metric system;SI只是公制中的一种
英制 British Imperial Units,美国采用英制
视重,表观重量(apparent weight):体重计的读出的重量(由于体重计与人整体作变速运动,所以体重计对人的支撑力是变化)
人静止时,
F
N
=
m
g
F_N = mg
FN=mg
人在下蹲时,经历先加速后减速最后静止,三个阶段
(以向下为正)
向下加速,速度向下,加速度方向向下,
m
g
−
F
N
=
m
a
mg - F_N = ma
mg−FN=ma,
F
N
<
m
g
F_N < mg
FN<mg,失重weightlessness
向下减速,速度向下,加速度方向向上,
m
g
−
F
N
=
−
m
a
mg - F_N = -ma
mg−FN=−ma,
F
N
>
m
g
F_N > mg
FN>mg,超重overweight
地球上物体所受的重力只是引力的一个分量
质量有三个属性:惯性质量、主动引力质量、被动引力质量
惯性质量(inertial mass): F = m a F=ma F=ma
引力质量(gravitational mass)主要可以分为主动引力质量和被动引力质量两种,前者决定物体产生的引力场的强弱,后者决定物体在处于其他引力场时受到的引力大小
F
=
G
M
1
M
2
r
2
F=G\frac{M_1 M_2}{r^2}
F=Gr2M1M2
狭义相对论提出,物体质量随其运动速度的增加而增加
m
=
m
0
1
−
(
v
/
c
)
2
m = \frac{m_0}{\sqrt{1-(v/c)^2}}
m=1−(v/c)2m0,静止质量
m
0
m_0
m0
狭义相对论提出,质能方程
E
=
m
c
2
E=mc^2
E=mc2
由上面两个方程可以得出
m
0
c
2
=
E
2
−
p
2
c
2
m_0 c^2 = \sqrt{E^2 - p^2 c^2}
m0c2=E2−p2c2,
p
=
m
v
p=mv
p=mv
广义相对论中的等效原理:惯性质量和引力质量是等价
推力
假设推力与质量流量恒定,气体出口速度均匀并指向轴向,总推力为:
F
=
m
˙
v
2
+
(
p
2
−
p
3
)
A
2
F= \dot{m} v_2 + (p_2 - p_3)A_2
F=m˙v2+(p2−p3)A2
当
p
2
=
p
3
p_2=p_3
p2=p3,则称发动机喷管具有最佳膨胀比
第一项为动量推力,是推进剂质量流量与排气相对飞行器的速度之积
第二项为压差推力,喷管出口面积乘以排气压力与环境压力之差
燃烧室内表面的压力
p
1
p_1
p1是最高的
喷管喉部压力
p
t
p_t
pt
喷管出口压力
p
2
p_2
p2
大气压
p
3
p_3
p3随高度增加而降低
补充—拉瓦尔喷管
拉瓦尔喷管,渐缩渐阔喷管(de Laval nozzle,convergent-divergent nozzle,CD nozzle,con-di nozzle)
拉瓦尔喷管三段中的流体为:亚音速、音速、超音速
等熵流动方程
d
p
p
=
γ
d
ρ
ρ
\frac{dp}{p} = \gamma \frac{d \rho}{\rho}
pdp=γρdρ
d
p
=
c
2
d
ρ
dp = c^2 d\rho
dp=c2dρ
其中,
γ
=
C
p
/
C
v
,
c
=
γ
p
ρ
=
γ
R
T
\gamma = C_p / C_v,c = \sqrt{\gamma \frac{p}{\rho}}=\sqrt{\gamma RT}
γ=Cp/Cv,c=γρp=γRT,
c
c
c为音速
动量守恒方程
u
′
=
u
+
(
Δ
v
Δ
x
)
Δ
x
u' = u + (\frac{\Delta v}{\Delta x})\Delta x
u′=u+(ΔxΔv)Δx
p
′
=
p
+
(
Δ
p
Δ
x
)
Δ
x
p' = p + (\frac{\Delta p}{\Delta x})\Delta x
p′=p+(ΔxΔp)Δx
A
′
=
A
A' = A
A′=A
Δ
x
Δ
t
=
v
\frac{\Delta x}{\Delta t} = v
ΔtΔx=v
m
=
ρ
Δ
x
A
m = \rho \Delta x A
m=ρΔxA
F
=
m
a
=
m
Δ
v
Δ
t
F = ma = m \frac{\Delta v}{\Delta t}
F=ma=mΔtΔv
p
A
−
p
′
A
′
=
m
v
′
−
v
Δ
t
pA - p'A' = m \frac{v' - v}{\Delta t}
pA−p′A′=mΔtv′−v,将上面前三个方程代入得
p
A
−
(
p
+
(
Δ
p
Δ
x
)
Δ
x
)
A
=
m
(
Δ
v
Δ
x
)
Δ
x
Δ
t
pA - (p + (\frac{\Delta p}{\Delta x})\Delta x) A = m \frac{(\frac{\Delta v}{\Delta x})\Delta x}{\Delta t}
pA−(p+(ΔxΔp)Δx)A=mΔt(ΔxΔv)Δx,再代入后两个方程得
−
Δ
p
Δ
x
=
ρ
v
Δ
v
Δ
x
- \frac{\Delta p}{\Delta x} = \rho v \frac{\Delta v}{\Delta x}
−ΔxΔp=ρvΔxΔv
−
d
p
d
x
=
ρ
v
d
v
d
x
- \frac{d p}{d x} = \rho v \frac{d v}{d x}
−dxdp=ρvdxdv
ρ
v
d
v
=
−
d
p
\rho v dv = - dp
ρvdv=−dp
连续性方程(质量守恒方程)
m
˙
=
ρ
v
A
=
C
\dot{m} = \rho v A = C
m˙=ρvA=C,
m
˙
\dot{m}
m˙为质量流量
l
n
ρ
+
l
n
v
+
l
n
A
=
l
n
C
ln \rho + ln v + ln A = ln C
lnρ+lnv+lnA=lnC
d
(
l
n
ρ
+
l
n
v
+
l
n
A
)
=
0
d(ln \rho + ln v + ln A) = 0
d(lnρ+lnv+lnA)=0
d
ρ
ρ
+
d
v
v
+
d
A
A
=
0
\frac{d \rho}{\rho} + \frac{d v}{v} + \frac{d A}{A} = 0
ρdρ+vdv+AdA=0
由动量守恒方程、等熵流动方程联立得
d
p
=
c
2
d
ρ
,
ρ
v
d
v
=
−
d
p
dp = c^2 d\rho, \rho v dv = - dp
dp=c2dρ,ρvdv=−dp
−
ρ
v
d
v
=
c
2
d
ρ
-\rho v dv = c^2 d\rho
−ρvdv=c2dρ
−
ρ
M
2
d
v
=
v
d
ρ
-\rho M^2 dv = v d\rho
−ρM2dv=vdρ
−
M
2
d
v
/
v
=
d
ρ
/
ρ
-M^2 dv / v = d \rho / \rho
−M2dv/v=dρ/ρ
再将上面方程代入质量守恒方程得
−
M
2
d
v
/
v
=
d
ρ
/
ρ
,
d
ρ
ρ
+
d
v
v
+
d
A
A
=
0
-M^2 dv / v = d \rho / \rho, \frac{d \rho}{\rho} + \frac{d v}{v} + \frac{d A}{A} = 0
−M2dv/v=dρ/ρ,ρdρ+vdv+AdA=0
−
M
2
d
v
/
v
+
d
v
v
+
d
A
A
=
0
-M^2 dv / v + \frac{d v}{v} + \frac{d A}{A} = 0
−M2dv/v+vdv+AdA=0
(
1
−
M
2
)
d
v
/
v
=
−
d
A
/
A
(1 - M^2) dv/v = -dA/A
(1−M2)dv/v=−dA/A
(
1
−
M
2
)
d
v
/
v
=
−
d
A
/
A
(1 - M^2) dv/v = -dA/A
(1−M2)dv/v=−dA/A
如果流动为超音速(
M
>
1
M>1
M>1),若截面积增大(
d
A
>
0
dA>0
dA>0),则
d
v
>
0
dv>0
dv>0,即速度增加
如果流动为超音速(
M
>
1
M>1
M>1),若截面积减小(
d
A
<
0
dA<0
dA<0),则
d
v
<
0
dv<0
dv<0,即速度减小
如果流动为亚音速(
M
<
1
M<1
M<1),若截面积增大(
d
A
>
0
dA>0
dA>0),则
d
v
<
0
dv<0
dv<0,即速度减小
如果流动为亚音速(
M
<
1
M<1
M<1),若截面积减小(
d
A
<
0
dA<0
dA<0),则
d
v
>
0
dv>0
dv>0,即速度增加
动量守恒
ρ
v
d
v
=
−
d
p
\rho v dv = - dp
ρvdv=−dp,由于
ρ
v
>
0
\rho v > 0
ρv>0,因此
d
v
dv
dv与
d
p
dp
dp反号,即速度增加,压力减小;速度减小,压力增加
亚音速流可视为不可压缩流,密度恒定,而质量守恒 ρ v A = C \rho v A = C ρvA=C,因此面积减小时,速度增大;面积增大时,速度减小。
超音速流为可压缩流,密度变化,质量守恒
ρ
v
A
=
C
\rho v A = C
ρvA=C,因此面积减小时,
ρ
v
\rho v
ρv增大;面积增大时,
ρ
v
\rho v
ρv减小。
如果流动为超音速(
M
>
1
M>1
M>1)
−
M
2
d
v
/
v
=
d
ρ
/
ρ
-M^2 dv / v = d \rho / \rho
−M2dv/v=dρ/ρ,则
−
ρ
/
v
>
d
ρ
/
d
v
- \rho / v > d \rho / dv
−ρ/v>dρ/dv,
0
>
d
ρ
/
d
v
0 > d \rho / dv
0>dρ/dv,因此
d
v
dv
dv与
d
ρ
d\rho
dρ反号。
如果流动为超音速(
M
>
1
M>1
M>1),
−
M
2
d
v
/
v
=
d
ρ
/
ρ
-M^2 dv / v = d \rho / \rho
−M2dv/v=dρ/ρ,则
∣
d
(
l
n
ρ
)
∣
>
∣
d
(
l
n
v
)
∣
|d(ln\rho)| > |d(lnv)|
∣d(lnρ)∣>∣d(lnv)∣,即速度变化量远小于密度变化量,
ρ
v
\rho v
ρv的变化主要取决于
d
ρ
d \rho
dρ。
因此,面积增大,
ρ
v
\rho v
ρv增大,则
d
ρ
>
0
,
d
v
<
0
d \rho >0, dv < 0
dρ>0,dv<0;面积减小,
ρ
v
\rho v
ρv减小,则
d
ρ
<
0
,
d
v
>
0
d \rho < 0, dv > 0
dρ<0,dv>0。
https://www.grc.nasa.gov/WWW/k-12/airplane/conmo.html
https://www.grc.nasa.gov/WWW/k-12/airplane/isentrop.html
https://www.grc.nasa.gov/WWW/k-12/airplane/tunnozt.html
排气速度
等效排气速度 c = v 2 + ( p 2 − p 3 ) A 2 / m ˙ c = v_2 + (p_2 - p_3)A_2/\dot{m} c=v2+(p2−p3)A2/m˙
特征排气速度 c ∗ = p 1 A t / m ˙ c^* = p_1 A_t / \dot{m} c∗=p1At/m˙
能量与效率
射流功率为单位时间内输出的喷射物质动能
对于恒定的气体喷射速度
v
v
v,射流功率
P
j
e
t
=
m
˙
v
2
2
P_{jet} = \frac{\dot{m} v^2}{2}
Pjet=2m˙v2
比功率为射流功率除以推进剂系统装填质量
化学火箭发动机输入功率
P
c
h
e
m
=
m
˙
Q
R
J
P_{chem} = \dot{m} Q_R J
Pchem=m˙QRJ
单位质量化学推进剂能获得的最大能量为其燃烧反应热
Q
R
Q_R
QR
化学火箭发动机燃烧效率是单位质量推进剂的实际反应热与理论反应热之比
燃烧效率乘以
P
c
h
e
m
P_{chem}
Pchem即为可用功率,可用功率被转换为排气射流的动能
化学火箭发动机内效率 η i n t = P j e t / ( η c o m b P c h e m ) \eta_{int} = P_{jet}/(\eta_{comb} P_{chem}) ηint=Pjet/(ηcombPchem)
推进剂热量 = 燃烧室可用能量 + 燃烧损失(不良混合、不完全燃烧)
燃烧室可用能量 = 排气射流总能量 + 壁面热量损失
排气射流总能量 = 排气射流动能 + 排气射流不可用热能
排气射流动能 = 可用于推进飞行器的能量 + 排气残余动能
飞行器功率
P
v
e
h
i
c
l
e
=
F
v
P_{vehicle}=Fv
Pvehicle=Fv
推进效率
η
p
=
P
v
e
h
i
c
l
e
/
(
P
v
e
h
i
c
l
e
+
(
w
˙
/
g
0
)
(
c
−
u
)
2
/
2
)
\eta_p = P_{vehicle}/(P_{vehicle} + (\dot{w}/g_0)(c-u)^2/2)
ηp=Pvehicle/(Pvehicle+(w˙/g0)(c−u)2/2)