• 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=107J
    1 e V = 1.60218 × 1 0 − 19 J 1 eV = 1.60218\times 10^{-19} J 1eV=1.60218×1019J
    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˙dt0tFdt,单位 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˙gdt0tFdt,单位 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$
    
    • 1
    • 2

    等效排气速度 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 ζ=(momf)/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 mgFN=ma F N < m g F_N < mg FN<mg,失重weightlessness
    向下减速,速度向下,加速度方向向上, m g − F N = − m a mg - F_N = -ma mgFN=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)2 m0,静止质量 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=E2p2c2 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+(p2p3)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/Cvc=γρ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} pApA=mΔtvv,将上面前三个方程代入得
    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 (1M2)dv/v=dA/A

    ( 1 − M 2 ) d v / v = − d A / A (1 - M^2) dv/v = -dA/A (1M2)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
    
    • 1
    • 2
    • 3

    排气速度

    等效排气速度 c = v 2 + ( p 2 − p 3 ) A 2 / m ˙ c = v_2 + (p_2 - p_3)A_2/\dot{m} c=v2+(p2p3)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)(cu)2/2)

  • 相关阅读:
    【华为机试真题 JAVA】报数游戏-100
    mysql连接池DataSource,DruidDataSource的理解及其使用
    文件共享服务器
    【tensorflow】制作自己的数据集
    非零基础自学Java (老师:韩顺平) 第4章 运算符 4.3 关系运算符(比较运算符)
    使命担当 守护安全 | 中睿天下获全国海关信息中心感谢信
    【Python 实战基础】Pandas如何将表格的前几行生成html
    CPU密集型、IO密集型
    Spring-JdbcTemplate环境搭配及基本功能
    无向图三元环计数(根号算法)
  • 原文地址:https://blog.csdn.net/qq_37083038/article/details/126331088