周岸峰,郭 振,李道奎,周仕明,姜人伟1,
(1. 国防科技大学 空天科学学院, 湖南 长沙 410073; 2. 空天任务智能规划与仿真湖南省重点实验室, 湖南 长沙 410073;3. 中国运载火箭技术研究院 北京宇航系统工程研究所, 北京 100076)
导弹发射井是指供陆基战略弹道导弹垂直贮存、准备和实施发射的地下工程设施[1]。无论是在地震等严酷自然环境中,还是常规武器轰炸、核打击等恶劣战场条件下,发射井只有具备充足的防护能力,才能有效提高战略导弹的生存和反击能力[2-3]。但随着大当量、高精度、分导式多弹头核武器的快速发展[4-5],特别是其命中精度的日益提高,导弹发射井受到的威胁也越来越大,发射井的加固和超级加固势在必行[6],如此才能确保己方在遭受核打击之后能够给对方致命的反击。然而,当发射井遭受核打击时,井体抗压能力越高,其“感受”的地冲击震动也越强烈[1],使得井内导弹及其发射装置也面临着强烈的冲击震动,大位移导致的失稳、侧翻和碰撞等潜在危险[7]。因此在加固井体强度的同时,也必须对井内的导弹及发射装置进行科学合理的减震设计。
目前,关于井下发射导弹的减震设计的研究报道相对较少。湘言[8]深入分析了美国井下发射的民兵导弹的摆式减震系统,建立了相应的二维分析模型并进行了数值计算与实验验证。陈钦佩[7]总结了美国导弹发射井的加固措施,并详细介绍了民兵导弹的摆式减震系统以及改进的悬吊式和斜吊式减震系统。张蕴华[9]简要介绍了美国MX导弹与发射筒之间、发射筒与发射井壁之间的隔震元件。赵国柱等[10]则简单描述了SS-19战略导弹的摆式液压气动减震系统。Wood等[11]在其研制的导弹发射筒悬吊系统中,利用径向对称布置的空气弹簧和底部环形空气弹簧分别抵抗横向和纵向的冲击。Kendall等[12]提出了一种吊笼式的井基导弹悬挂系统,能有效减小横向和垂直方向的冲击震动。Larson等[13]设计了一种同心筒式导弹垂直发射装置,其垂直减震器像倒置的摆式减震系统一样布置,只承受压力载荷。
然而上述工作多集中于调研总结,详细深入的研究较少。Chobotov[14]将使用摆式悬挂系统的井基导弹视为刚体,推导了导弹的运动学方程,获得了导弹受到激励作用后的稳定性判据。Jiang等[15]建立了舰载导弹四联发射筒的精细有限元模型,对发射筒结构的抗冲击能力进行了评估。张筱等[16]基于俄罗斯井下冷发射的SS-18导弹发射平台建立了动力学仿真模型,并分析了减震装置的刚度对导弹发射初始扰动的影响。荣吉利等[17]分析了不同减震系统对发射筒位移与倾斜角度的影响,弹体与发射筒之间相对倾斜的影响以及弹体自身速度、加速度的影响,为井下导弹悬挂系统选型提供了重要参考。但荣吉利等在研究过程中将发射筒与弹体均等效为刚体,忽略了发射筒和弹体自身柔性对结构动力响应的影响,不利于减震系统设计参数的确定。
本文进一步考虑弹-筒系统结构柔性的影响,分析爆炸地冲击载荷作用下不同减震系统中弹体和发射筒关键部位的动力响应特性,并以悬吊式减震系统为例讨论减震器刚度、阻尼对减震性能的影响,为井下冷发射导弹的减震设计与优化提供有效参考。
对于地下防护工程,导弹等武器触地爆炸时引起的水平地冲击震动是主要的危害,在爆炸近区地冲击加速度波形[18-19]可以表示为:
(1)
其中:τ、e1、e2、s和c的计算如式(2)所示;T为加速度时程持续时间,计算如式(3)所示;Am为加速度峰值,可根据美国陆军技术手册《常规武器防护原理》[20]计算,其表达式为如式(4)所示。
(2)
(3)
式中:λ为常系数,取值范围为0.656 48~1.312 96,本文取平均值0.984 72;R为爆炸距离;cP为纵波波速。
(4)
式中,f为爆炸能量耦合系数,α为衰减系数,Q为爆炸当量。
以爆炸当量为2×105t TNT、爆炸距离为1 000 m时的触地爆作为计算条件,根据文献[20]选取爆炸能量耦合系数为0.6、纵波波速为1 500 m/s、衰减系数为2.375。从而,按式(3)和式(4)计算得到地冲击过程的加速度时程持续时间和加速度峰值分别为0.656 5 s、14.00 m/s2,最终得到如图1所示的加速度时程曲线。
图1 爆炸地冲击加速度时程曲线Fig.1 Acceleration history of blast-induced ground shock
冷发射导弹的三种典型减震系统[17]如图2所示。在三种减震系统中,弹体都储存在发射筒内,弹体与筒壁之间安装适配器以保持弹体稳定,并起到一定的缓冲作用。悬吊式减震系统中,发射筒与井壁之间用垂直减震器悬挂于发射井内,径向布设若干个液压弹簧和阻尼器,如俄罗斯的SS-18导弹。下支承式减震系统则没有悬挂装置,由支承底座提供垂直减震缓冲功能,由填充于发射筒与井壁之间的特制发泡塑料实现水平减震(如美国MX导弹),也可使用弹簧阻尼装置实现水平减震。斜吊式减震系统通过斜置的减震器将发射筒悬挂于井内,其水平震动与垂直震动相互耦合。
图2 典型减震系统Fig.2 Typical shock absorption system
由于水平地冲击作用对弹-筒系统的危害最大,故主要考虑水平减震措施的设计,借助MSC.Patran前处理软件进行建模,将导弹等效为梁模型,发射筒和适配器均等效为均质弹性体,减震器等效为弹簧阻尼单元,减震器与井壁的连接点通过多点约束(multi-point constraint, MPC)单元与一参考点耦合,用来模拟刚性井壁,对参考点施加图1所示的加速度激励。有限元模型主要参数见表1[17, 21-22]。
表1 有限元模型主要参数Tab.1 Major parameters of finite element models
为合理进行有限元网格划分,对有限元模型的网格参数进行了多次预估,表2给出了最后两次估计的网格参数,图3是不同网格参数对应的模型的仿真结果。由图可知,网格加密前后的结果一致性良好,相对误差小于0.6%,由此说明表2所示的网格参数设置是合理有效的。为兼顾计算效率和计算精度,选取表2中“加密前”对应的模型网格。
表2 有限元网格参数Tab.2 Parameters of finite element mesh
图3 有限元网格划分的合理性验证Fig.3 Verification of the rationality of finite element mesh division
考虑弹体上敏感设备对加速度的耐受性、弹体与发射筒碰撞以及发射筒与井壁碰撞的可能性,将弹体纵向各关键部位的加速度、弹-筒之间的相对位移以及发射筒的位移作为评价减震效果的影响因素,并引入传递系数作为评价减震效果的最直接指标[12]。传递系数η定义为传递到弹体上的最大加速度max{a(x,t)}与基础运动的最大加速度Am之比,表达式为:
(5)
为分析结构柔性对弹-筒系统动力响应的影响,以悬吊式为例,将发射筒和弹体的弹性模量提高1 000倍以模拟刚体,并与柔性体模型进行对比,相关结果如图4~6所示。由图可知,与柔性体模型相比,刚体模型中弹体大部分位置的加速度响应峰值和弹-筒之间的相对位移响应峰值均大大降低,而发射筒大部分位置的位移响应峰值略有升高。因此,若将发射筒和弹体等效为刚体进行建模,将导致减震性能设计不足,从而严重影响导弹的使用安全。为合理确定减震系统的设计参数必须进一步考虑结构柔性的影响。
图4 弹体加速度响应峰值Fig.4 Peak acceleration response of missile
图5 弹-筒相对位移响应峰值Fig.5 Peak relative displacement response between missile and canister
图6 发射筒位移响应峰值Fig.6 Peak displacement response of canister
在爆炸地冲击作用下,不同减震系统中弹-筒系统动力响应峰值沿弹体纵向的分布规律如图7~9所示,动力响应的传递系数如表3所示。可以看出:①悬吊式和下支承式减震系统的减震效果比较接近,前者略好。这是因为支承台也传递了地冲击载荷,导致系统动力响应稍有增大,但支承台水平方向的刚度仅为水平减震器刚度的1.15%,传递作用相对很小。因此,为避免下支承台过多地传递地冲击载荷,应适当降低其水平方向的刚度。②与悬吊式和下支承式相比,整体上斜吊式的减震效果要差。这与斜吊式水平和垂直震动相互耦合导致结构局部振动效应增强或减弱有关,如图6所示,弹体中部加速度响应峰值增大而尾部加速度响应峰值减小。另一方面,由于斜吊式减震系统水平和垂直隔震不解耦,在重力作用下,头部的减震器与水平线的夹角变大,减震器刚度的水平分量变小;而尾部的减震器与水平线的夹角变小,减震器刚度的水平分量变大,从而出现图9所示的发射筒头部位移变化增大而尾部位移变化减小的现象。③由图7和图8可知,弹体纵向23 m附近的加速度响应峰值以及弹-筒相对位移响应峰值明显大于其他位置,原因在于弹体通过适配器约束在发射筒内,相当于多点简支梁,适配器之间的间隔距离越长,简支梁的跨度越大,因此简支梁中间部位的动力响应和变形都比较大。
图7 各减震系统中弹体加速度响应峰值Fig.7 Peak acceleration response of missile in different shock absorption systems
图8 各减震系统中弹-筒相对位移响应峰值Fig.8 Peak relative displacement response between missile and canister in different shock absorption systems
图9 各减震系统中发射筒位移响应峰值Fig.9 Peak displacement response of canister in different shock absorption systems
表3 各减震系统的传递系数Tab.3 Transfer coefficients of different shock absorption systems
以悬吊式减震系统为例,当水平减震器阻尼保持在4×105N·s/m时,不同减震器刚度下系统的动力响应分布如图10~12所示,传递系数如表4所示。可以看出,当减震器刚度从2×107N/m逐渐降至2×105N/m时,在相同激励作用下减震器由于刚度降低而行程增大,导致发射筒位移响应增大;同时使得通过适配器传递给弹体的载荷减弱,因而弹体弯矩、加速度响应峰值以及弹-筒之间相对位移响应峰值均大大减小。而以悬吊减震系统为例,当水平减震器刚度保持在2×106N/m时,不同减震器阻尼下系统的动力响应分布如图13~15所示,传递系数如表5所示。可以看出,当减震器阻尼从4×105N·s/m逐渐降至5×104N·s/m时,弹体加速度响应峰值和弹-筒之间相对位移响应峰值均随之减小,发射筒位移响应峰值随之增大。原因在于地冲击加速度激励是一个幅值较大而作用时间极短的冲击过程,速度变化剧烈,阻尼器在冲击瞬间被“刚化”,还未来得及耗散冲击能量就已经提供了巨大的阻尼力,虽然有利于发射筒恢复稳定,但同时也导致传递至弹体的载荷偏高,且阻尼越大,这种“刚化”效应越强,因此出现弹体动力响应峰值随减震器阻尼减小而减小的现象。
图10 不同减震器刚度下弹体加速度响应峰值Fig.10 Peak acceleration response of missile under different shock absorber stiffness
图11 不同减震器刚度下弹-筒相对位移响应峰值Fig.11 Peak relative displacement response between missile and canister under different shock absorber stiffness
图12 不同减震器刚度下发射筒位移响应峰值Fig.12 Peak displacement response of canister under different shock absorber stiffness
表4 不同减震器刚度下的传递系数Tab.4 Transfer coefficients under different shock absorber stiffness
表5 不同减震器阻尼下的传递系数Tab.5 Transfer coefficients under different shock absorber damping
图13 不同减震器阻尼下弹体加速度响应峰值Fig.13 Peak acceleration response of missile under different shock absorber damping
图14 不同减震器阻尼下弹-筒相对位移响应峰值Fig.14 Peak relative displacement response between missile and canister under different shock absorber damping
图15 不同减震器阻尼下发射筒位移响应峰值Fig.15 Peak displacement response of canister under different shock absorber damping
当减震器刚度由2×105N/m降至2×104N/m时,系统动力响应峰值略有增大,则与弹-筒耦合振动增强有关。
减震器阻尼与弹-筒系统恢复稳定的时间息息相关,直接决定了导弹反击发射准备时长,因此其也是影响导弹综合作战能力的重要因素。图16给出了不同减震器阻尼下弹体尾部的加速度响应时程曲线。由图可知,减小减震器阻尼虽然能降低弹体加速度响应峰值,但不利于抑制弹体的晃动;对比来看,减震器阻尼取2×105N·s/m为宜。
图16 不同减震器阻尼下弹体尾部的加速度响应Fig.16 Acceleration response of missile tail under different shock absorber damping
本文开展了弹-筒系统抗爆减震设计研究,比较了悬吊式、下支承式、斜吊式等多种减震系统的减震效果,并分析了减震器刚度和阻尼对减震效果的影响,得到以下结论:
1)将发射筒和弹体等效为刚体进行建模,将导致减震性能设计不足,从而严重影响导弹的使用安全,为合理确定减震系统的设计参数必须进一步考虑结构柔性的影响。
2)悬吊式减震系统的减震效果略好于下支承式,斜吊式的减震效果最差;使用下支承式减震系统时,应适当降低支承台的水平刚度。
3)对于悬吊式减震系统,弹体加速度响应偏大的部位一般为头部、中部和尾部,弹体与发射筒之间最容易发生碰撞的部位为中部。
4)在一定范围内,悬吊式减震系统的减震效果随减震器刚度的降低或减震器阻尼的减小而显著提升;但刚度过小容易导致缓冲行程增大,而阻尼过小会使弹体趋于稳定的时间更长,设计时应综合考虑。