摘要
飞行器着陆或者坠落过程中会产生较大的冲击,为了改善乘员所受到的冲击损伤,提出了半主动磁流变库伦阻尼缓冲座椅系统。为分析缓冲座椅的力学特性,首先,建立座椅系统力学模型,利用拉格朗日方程建立座椅系统运动微分方程;其次,基于系统运动微分方程,分析系统刚度、阻尼以及不同跌落高度对缓冲座椅力学特性的影响;然后,分析了磁流变库伦阻尼器的安装方式对系统缓冲性能的影响,并利用T‑S模糊模型实现座椅系统半主动控制,将系统的极限距离缓冲效果与最优缓冲效果进行对比;最后,通过缓冲座椅原理样机的跌落试验对上述理论计算结果进行验证。研究发现,相比于传统的线性缓冲器,所提出的缓冲座椅能够适应不同的跌落冲击强度,具有良好的缓冲性能。
直升飞机着陆、航空器坠落过程中与地面碰撞时,所产生的冲击经过结构框架后传递给座椅,最后作用在人体上,从而对人体产生损伤作
剪式支撑结构具有结构简单紧凑的优点,方便安装弹性元件和阻尼元件,目前在座椅悬架中主要应用在隔振领
在剪式支撑结构上安装可变阻尼元件能够实现半主动控制。半主动控制系统能够针对不同的冲击强度来设置阻尼元件参数,从而适应不同冲击环境。磁流变阻尼器具有响应快、阻尼可变的特性,因此被应用在飞行器起落架、悬架座椅上完成结构隔振或缓
传统的磁流变阻尼器输出阻尼力受速度影响较大,在冲击作用下输出阻尼力十分不稳定。库伦阻尼器输出阻尼力与加载速度相关性较小,能够在冲击作用下稳定输出阻尼力。笔者基于一种具有库伦阻尼特性的磁流变库伦阻尼器,提出一种半主动磁流变库伦阻尼缓冲座椅系统,以提高缓冲座椅在不同冲击强度下的缓冲效果。
半主动缓冲座椅系统原理如

图1 半主动缓冲座椅系统原理
Fig.1 Schematic of buffing seat with semi-active control
将

图2 半主动缓冲座椅系统简化原理图
Fig.2 Structural schematic of scissor support structure
通常情况下,缓冲座椅在实际跌落过程中质量块动能不断增加,在与地面碰撞之后,剪式支撑结构才开始缓冲。在缓冲过程中,质量块的动能转化为弹簧的弹性势能和阻尼器消耗产生的热能。为说明缓冲座椅在冲击过程中的受力状态,对剪式支撑结构运动进行分析,如

图3 剪式支撑结构运动分析
Fig.3 Motion analysis of scissor support structure
在剪式支撑结构发生变形过程中,假设杆长为2l,质量块相对基础的变形为,那么可以获得剪式支架铰链的转角关于垂向相对变形的关系式为
(1) |
因此,水平移动距离关于垂直相对变形的关系式为
(2) |
假设质量块向上运动和弹簧压缩变形为正向,那么当质量块向下运动时,弹簧处于拉长状态,所以符号为负。在整个缓冲过程中,假设剪式支撑结构的2根杆为轻质杆,不考虑其质量及重力作用,那么系统的动能主要为质量块的动能
(3) |
在缓冲过程中系统的势能主要表现为弹性元件的弹性势能
(4) |
为方便计算,将跌落过程中质量块所产生的速度设置为系统的初始条件。根据拉格朗日方程,有
(5) |
其中:;为垂向库伦阻尼力。
由于缓冲座椅中的磁流变库伦阻尼器在实际试验过程中呈现出与库伦阻尼相同的特性,因此在建模过程中将该阻尼等效为库伦阻尼进行计算。这2个阻尼均表现为库伦阻尼特征,即
(6) |
其中:为库伦阻尼器阻尼力;表示质量块运动方向;为坐标变换关系。
(7) |
其中
(8) |
由
(9) |
(10) |
当质量块向下变形时,恢复力向上,其符号为负。根据
(11) |
根据

图4 座椅垂向恢复力及刚度
Fig.4 Seat vertical restoring force and stiffness
系统在垂直方向上的等效阻尼力绝对值为
(12) |
由

图5 等效阻尼力与垂向变形的关系
Fig.5 The relationship between equivalent damping force and vertical deformation
由
笔者通过四阶龙格库塔数值求解方法计算运动微分方程(9),从而获得缓冲座椅冲击响应解。弹簧刚度对系统冲击响应的影响如


图6 系统刚度对冲击响应的影响
Fig.6 The influence of system stiffness on impact response
计算2种刚度条件下系统的冲击响应,分别为低刚度和高刚度。当刚度值为时,剪式支撑结构在不同跌落高度冲击下的系统冲击响应如
系统的负向加速度响应随着跌落高度的增加而增加,但是正向加速度响应峰值却保持在一定的恒定值,这是由于系统刚度为时,负向缓冲距离较长并大于临界位移值,因此最大加速度响应等于系统临界位移值对应的加速度。由
由于在不同跌落高度下系统响应规律基本一致,这里以跌落高度为0.2 m的情况进行说明。当剪式支撑结构仅水平安装库伦阻尼器时,零刚度情况下水平库伦阻尼对缓冲效果的影响如

图7 零刚度情况下水平库伦阻尼对缓冲效果的影响
Fig.7 Influence of horizontal Coulomb damping on cushioning effect under zero stiffness
为了合理利用剪式支撑结构的非线性特征,将库伦阻尼器水平安装在剪式支撑结构上,根据

图8 水平库伦阻尼对系统冲击响应的影响
Fig.8 Influence of horizontal Coulomb damping on system shock response
由
在不同的跌落高度下,需要通过计算获得最优阻尼值使得系统具有最优缓冲效果。在跌落过程中,乘员的耐冲极限可以用加速度响应来衡量,因此需要以系统的最大加速度响应为优化目标,库伦阻尼值为参变量,将系统的最大缓冲距离设定为边界条件。
一般结构的安装空间尺寸是有限的,缓冲结构都有一个极限缓冲距离,当最大位移响应超过该缓冲距离,系统就会发生二次刚性碰撞,从而导致人体二次碰撞,该冲击损伤可能比缓冲损伤更加严重。因此,需要将系统的最大缓冲距离设为最大可移动距离的2/3,留出一定的阈值避免发生二次碰撞。另外,由于剪式支撑结构的恢复力随着变形增加而降低,为了保证系统在极限距离处的恢复能够支撑系统重力,需要将刚度设为65 kN/m,质量块质量为60 kg,杆长为0.3 m,剪式座椅初始角为45°,极限位置恢复力为1 252 N。为了简化计算,将跌落高度所产生的冲击速度作为质量块的初始速度。系统的运动微分方程可以简化为
(13) |
本算例中最大缓冲距离为0.14 m,当跌落高度超过0.4 m时,无阻尼系统的最大位移响应就会超过最大缓冲距离。为保证系统具有足够大的位移响应,这里仅计算0.4,0.5和0.6 m这3种跌落高度下系统的冲击响应,水平库伦阻尼系统响应如

图9 水平库伦阻尼系统响应
Fig.9 Horizontal Coulomb damping system response
由
为进一步了解最优缓冲效果与极限缓冲效果的具体过程,分别计算了3种不同高度条件下时域缓冲曲线,水平库伦阻尼系统最优加速度响应和极限缓冲加速度响应分别如图

图10 水平库伦阻尼系统最优加速度响应
Fig.10 Optimal acceleration response of horizontal Coulomb damping system

图11 水平库伦阻尼系统极限缓冲加速度响应
Fig.11 Limit buffering acceleration response of horizontal Coulomb damping system
水平安装阻尼具有可变阻尼特征,也为系统的变刚度、变阻尼特性提供了可能性。本研究通过计算装有不同缓冲器座椅的冲击响应,对水平阻尼安装、垂向阻尼安装和传统库伦阻尼器进行对比说明。首先,对安装有垂直库伦阻尼器的剪式支撑座椅进行计算,当系统的刚度相同时,在不同高度下系统的阻尼值、位移响应以及加速度响应如12所示。由图可以看出,垂向安装阻尼器的剪式支撑缓冲座椅与水平安装阻尼器的响应规律相类似,在极限缓冲距离内均存在一个最优缓冲效果。

图12 垂直库伦阻尼系统响应
Fig.12 Vertical Coulomb damping system response
为了进一步了解最优缓冲效果与极限缓冲效果的具体过程,分别计算了3种不同高度条件下的时域缓冲曲线,垂直库伦阻尼系统最优加速度响应和极限缓冲加速度响应分别如图

图13 垂直库伦阻尼系统最优加速度响应
Fig.13 Optimal acceleration response of vertically Coulomb damping system

图14 垂直库伦阻尼系统极限缓冲加速度响应
Fig.14 Limit buffering acceleration response of vertically Coulomb damping system
从上述计算结果可以发现,安装有垂直库伦阻尼器的剪式支撑结构的冲击响应与水平库伦阻尼器规律基本相同,但垂直库伦阻尼器所设置的阻尼值要小于水平库伦阻尼器,这是因为水平库伦阻尼器在垂直方向上的等效阻尼力值随着变形增加而降低,在库伦阻尼器输出阻尼力相同的情况下,水平库伦阻尼器的有效输出值要低于垂直库伦阻尼器。
为了进行对比分析,还计算了无刚度水平库伦阻尼器以及经典线性库伦阻尼器系统的冲击响应。不同结构条件下系统响应对比如

图15 各类型阻尼器系统响应对比
Fig.15 System response comparison of various types of dampers
在本研究所计算的工况中,水平库伦阻尼器和垂直库伦阻尼器的最优缓冲效果均优于极限缓冲效果、无刚度水平阻尼器和经典刚度库伦阻尼系统。由
T‑S模糊控制能够根据实际经验值或者理论计算数据库,完成冲击载荷与输出阻尼力之间的参数化建模,能够针对特定冲击载荷输出最优阻尼力。

图16 T-S模糊控制模型
Fig.16 T-S fuzzy control model
基于上述参数,根据冲击输入载荷就能够完成输出阻尼力D1u控制。针对缓冲座椅水平安装阻尼器的情况,基于最优阻尼响应计算结果,建立冲击输入载荷与输出阻尼之间的参数化模型,如

图17 参数化模型
Fig.17 Parametric modeling
通过将T‑S模糊控制模型写入控制器中,在特定跌落冲击强度作用下,缓冲座椅所需要的最优缓冲阻尼力可以按照
与无刚度水平库伦阻尼器相比,具有刚度特性的水平库伦阻尼器系统缓冲优势更明显。当系统具有最优阻尼时,则有更小的加速度响应且缓冲距离更短。为了验证上述结果,对水平库伦阻尼系统进行跌落试验测试。磁流变库伦阻尼器扭矩测试曲线如

图18 磁流变库伦阻尼器扭矩测试曲线
Fig.18 Torque test curve of magnetorheological damper
本次试验中,假人体质量为60 kg,磁流变库伦阻尼器的阻尼力可调范围为500~4 000 N。以0.5和0.6 m的跌落高度为例,结合理论计算结果进行对比。缓冲座椅跌落试验和跌落冲击响应曲线分别如图

图19 缓冲座椅跌落试验
Fig.19 Drop test of the buffing seat

图20 缓冲座椅跌落冲击响应曲线
Fig.20 Shock response of buffing seat under drop test
在理论计算过程中,当跌落高度为0.5 m时,无弹簧系统和有弹簧系统对应的库伦阻尼力分别为3 800 和1 920 N。当跌落高度为0.6 m时,无弹簧系统和有弹簧系统对应的库伦阻尼力分别为4 600和2 676 N。
与试验相比,理论计算所使用的库伦阻尼力要大于试验中所设置的库伦阻尼力,这是因为在原理样机中,除了磁流变液阻尼器的阻尼力作用,还存在系统自身结构之间的摩擦阻尼作用,使实际最优阻尼力与理论最优阻尼力产生了差别。
当系统安装弹簧时,弹簧在缓冲过程中可以将一部分冲击能量存储为弹性势能,并且在释放弹性势能的过程中还可以通过阻尼器消耗能量。由于无刚度水平安装阻尼器系统无弹性元件,因此整个缓冲过程仅靠阻尼器消耗能量。为了保证系统能够在极限距离处完成缓冲,必须将阻尼力提升到一定值。通过观察可以发现:当系统无弹簧安装时,试验数据与理论数据响应趋势基本一致,如
1) 缓冲座椅的剪式支撑结构具有较强的结构非线性,系统体现出显著的变刚度、变阻尼特性。当缓冲座椅负向变形时,系统的刚度阻尼特性均呈现出“软特性”,从而保证系统负向变形过程中加速度响应维持在极限值范围内。
2) 当跌落高度低于0.5 m时,安装水平库伦阻尼器的系统具有最优缓冲效果;当跌落高度为0.6 m时,水平库伦阻尼器系统的最优缓冲效果与垂直安装阻尼器相接近,但是其缓冲效果仍然优于其他系统。
3) 半主动阻尼控制可以使系统方便获得最优阻尼力,从而实现不同跌落冲击强度的最优缓冲效果。通过跌落试验与理论计算结果对比可知,冲击响应曲线趋势基本吻合。
参考文献
闫振. 垂直起降可重复使用运载器的支腿缓冲器设计及分析[D]. 哈尔滨:哈尔滨工业大学,2019. [百度学术]
谢燕. 飞船返回舱着陆冲击缓冲座椅系统改进设计研究[D].长沙:国防科学技术大学,2010. [百度学术]
熊鹏宇,周云波,胡文海,等.车辆防护座椅性能及缓冲吸能器性能分析[J].科学技术与工程,2019,19(27):381-386. [百度学术]
XIONG Pengyu, ZHOU Yunbo, HU Wenhai,et al. Research on performance of vehicle anti-shock seat and parameters of energy absorber[J]. Science Technology and Engineering, 2019, 19(27):381-386.(in Chinese) [百度学术]
邓先来, 王翔. 航空座椅抗坠毁环境模拟技术分析[J].直升机技术,2017(4):23-28. [百度学术]
DENG Xianlai, WANG Xiang. Study the simulation technique of anti-crash environment about airline seat [J]. Helicopter Technique, 2017(4):23-28. (in Chinese) [百度学术]
晏涛, 薛萌萌, 肖光润, 等. X型结构航空座椅吸能参数确定方法[J].装备制造技术,2018(6):139-141. [百度学术]
YAN Tao, XUE Mengmeng, XIAO Guangrun,et al. Method for determining energy absorption parameters of X-type structure aviation seat[J]. Equipment Manufacturing Technology, 2018(6):139-141. (in Chinese) [百度学术]
江浩斌,胡贝,聂佳梅,等.含惯容器的剪式座椅悬架动态特性研究[J].车辆与动力技术,2013(4):11-16. [百度学术]
JIANG Haobin, HU Bei, Nie Jiamei,et al. Research on dynamic characteristic of ISD scissors seat suspension [J]. Vehicle & Power Technology, 2013(4):11-16. (in Chinese) [百度学术]
朱思洪,王家胜,王敏娜.带附加气室空气悬架剪式座椅振动特性试验研究[J].振动与冲击,2009,28(11):104-106,206. [百度学术]
ZHU Sihong, WANG Jiasheng, WANG Minna. Experimental study on vibration characteristics of scissors linkage seat with air suspension attache[J]. Journal of Vibration and Shock, 2009,28(11):104-106,206. (in Chinese) [百度学术]
徐晓美,朱思洪.基于MBS的剪式座椅静态和动态特性仿真研究[J].机械科学与技术,2006(6):684-686,731. [百度学术]
XU Xiaomei, ZHU Sihong. Simulation research on the static and dynamic characteristics of a scissors linkage seat based on MBS [J]. Mechanical Science and Technology, 2006(6):684-686,731. (in Chinese) [百度学术]
WU Z, JING X, SUN B, et al. A 6DOF passive vibration isolator using X-shape supporting structures[J]. Journal of Sound and Vibration, 2016, 380:90-111. [百度学术]
董小闵,李军礼,于建强,等.月面低空飞行器着陆缓冲机构设计与仿真分析[J].载人航天,2019,25(6):779-782,798. [百度学术]
DONG Xiaomin, LI Junli, YU Jianqiang,et al. Design and simulation analysis of landing gear for lunar low altitude aircraft [J]. Manned Spaceflight,2019,25(6):779-782,798. (in Chinese) [百度学术]
董小闵,余淼,廖昌荣,等.飞行器磁流变自适应半主动冲击缓冲器[J].西南交通大学学报,2009,44(5):748-752. [百度学术]
DONG Xiaomin, YU Miao, LIAO Changrong,et al. Adaptive semiactive bumper for aircraft[J]. Journal of Southwest Jiaotong University, 2009,44(5):748-752. (in Chinese) [百度学术]
ZHU X X, ZHU S H. A theoretical model for calculating vibration characteristics of a kind of driver seat with air spring and mr damper[J]. Applied Mechanics and Materials, 2012, 141:8-14. [百度学术]
LIU C, JING X, LI F. Vibration isolation using a hybrid lever-type isolation system with an X-shape supporting structure[J]. International Journal of Mechanical Sciences, 2015, 98:169-177. [百度学术]
NING D, SUN S, DU H, et al. Vibration control of an energy regenerative seat suspension with variable external resistance[J]. Mechanical Systems and Signal Processing, 2018, 106:94-113. [百度学术]
水奕洁, SUBHASH R, 上官文斌.剪式悬架座椅等效刚度阻尼的计算与分析[J].振动与冲击,2016,35(19):38-44. [百度学术]
SHUI Yijie, SUBHASH R, SHANGGUAN Wenbin. Calculation and analysis for equivalent stiffness and damping of a scissor type seat-seat suspension[J]. Journal of Vibration and Shock, 2016,35(19):38-44. (in Chinese) [百度学术]
王利娟, 闫建国, 侯占峰, 等. 剪式座椅结构参数变化对其减振性能的影响分析[J]. 机电工程, 2019, 36(4):368-373. [百度学术]
WANG Lijuan, YAN Jianguo, HOU Zhanfeng, et al. Influence analysis for structural parameters of scissor type seat on seat′s vibration characteristic[J]. Journal of Mechanical & Electrical Engineering, 2019,36(4):368-373. (in Chinese) [百度学术]
SUN X, JING X, XU J, et al. Vibration isolation via a scissor-like structured platform[J]. Journal of Sound and Vibration, 2014, 333(9):2404-2420. [百度学术]