3 计算
将KKV的表面分成5个流动区域,从前到后分别为1区,2、3区,4、5区,从上到下分别为2、4区,1区和3、5区。明智的导引法能够尽可能的使用RCS和TVC提供过载,因此假设AOA为0。假设THAAD拦截弹速度为8马赫,假设空气比热比γ为1.4,假设空气密度ρ为1.225kg/m^3。THAAD拦截弹的阻力表达为
D=p2sinδS2+p3sinδS2+p4sin(δ-θ)S4+p5sin(δ-θ)S5-p6S6
由于KKV速度为8Ma,因此Ma^2>>1。在高超音速条件下,圆锥激波角近似表示为
β=γ+1/γ+3[1+(1+2(γ+3)/(γ+1)^2Ma^2δ^2)^0.5]δ
β=6/11*[1+(1+8.8/(2.4^2)*64*(π/20)^2)]*π/20
β=0.09
由等熵条件可知,2、3区压力系数为
Cp2=Cp3=4δ^2/γ+1(β^2/δ^2-1/Ma^2δ^2)+[2(γ+1)Ma^2δ^2]/2δ^2/[β^2+(δ-1)Ma^2δ^2]*(β/δ-1)^2
Cp2=Cp3=4*(π/20)^2/2.4*[0.09^2/(π/20)^2-1/64*(π/20)^2)+4.8*64*(π/20)^2/[2*(π/20)^2/0.09^2+0.4*64*(π/20)^2]*(0.09/π/20-1)^2
Cp2=Cp3=0.021
P2=P3=ρV^2Cp2/2
P2=P3=1.225*(5000)^2*0.021/2
P2=P3=321562.5
根据斜激波公式,V2^2及V3^2为
Ma2^2=Ma3^2=【[Ma1+2/(γ-1])/(2γ/γ-1)*Ma1^2sin^2β-1】+Ma1^2cos^2β/(γ-1/2)Ma1^2sin^2β+1
Ma2^2=Ma3^2=64+(2/0.4)/(2.8/0.4)*64sinβ-1+64cos^2β/(0.2*64*sin^2β+1)
Ma2^2=Ma3^2=74.53
2、3区空气流至4、5区产生普朗特-迈耶流动,由P-M关系得
V(Ma4)-v(Ma2)=λ
其中,v(Man)=1+1/2*(γ-1)Ma^2
V(Ma2)=1+0.2*74.53
V(Ma2)=16.9
v(Ma4)=v(Ma5)=0.06+16.9
V(Ma4)=v(Ma5)=16.96
由等熵条件可知
P4=p5=[v(Ma2)/v(Ma4)]^(γ/γ-1)*p2
P4=p5=(16.9/16.96)^3.5*321562.5
P4=p5=317060.625
现在,阻力表达式中所有参数已知
D=p2sinδS2+p3sinδS3+p4sin(δ-θ)S4+p5sin(δ-θ)S5-p6S6
D=2p2sinδS2+2p4sin(δ-θ)S4-p6S6
D=p2So1+p4(So2-So1)-p6So2
其中So1为前锥体底面积,So2为后锥体底面积
D=2*321562.5*0.1*0.1*3.14+2*317060.625*(0.18*0.18-0.1*0.1)*3.14-101300*0.18*0.18*3.14
D=20194.125+22300.776-10305.86
D=32189.041N
下面计算HTV2的阻力。HTV2飞行时的AOA未知,因此难以使用斜激波公式计算,但可以使用其它参数估算。设升力为L,重力为W。已知HTV-2的升阻比在3.5-4.0之间,假设HTV-2突防时过载为20g。在升力系数最大时升阻比肯定不为最大值,即HTV-2在过载为20g时升阻比小于4。
L=20W
D>L/4
D>5W
D>5*9000=45000N
此外,HTV-2的长、宽、高都大于KKV,浸润面积更大,因此摩擦阻力也更大。我认为,HTV-2阻力大于THAAD拦截弹。