直升机滑行或滑跑状态下地面共振的分析方法技术

本申请提供了一种直升机滑行或滑跑状态下地面共振的分析方法，包括：获取直升机的轮胎的侧向刚度系数、减摆器的等效阻尼；获取横滚振型的振频，并根据横滚振型的振频，得到减摆器的有效质量、有效刚度和有效阻尼；将得到的减摆器的等效阻尼、有效质量、有效刚度和有效阻尼代入旋翼‑机体系统耦合方程，得到滑行或滑跑状态下地面共振特性。

【技术实现步骤摘要】
直升机滑行或滑跑状态下地面共振的分析方法
本申请涉及直升机动力学
，具体提供一种直升机滑行或滑跑状态下地面共振的分析方法。
技术介绍
滑行、滑跑是直升机常见的地面使用状态。在某些使用条件下，直升机必须带速度滑跑着陆或滑跑起飞，除了起落架机轮刹车速度要求外，还需从地面共振角度，给出基于该型直升机动力学稳定性分析的滑行滑跑速度限制，确保滑行滑跑状态的运转安全。直升机滑行、滑跑状态地面共振特性不同于垂直起落，一方面机轮侧向刚度会随速度增加而降低，另一方面桨叶挥舞运动会产生频率为Ω的哥氏力，也就是说旋翼系统同时承受着ωξ和Ω两个振频的简谐运动激励，此时桨叶阻尼器的阻尼特性不能用单频做功来进行处理。
技术实现思路
为了解决上述技术问题至少之一，本申请提供了一种直升机滑行或滑跑状态下地面共振的分析方法，包括：获取所述直升机的轮胎的侧向刚度系数、减摆器的等效阻尼；获取横滚振型的振频，并根据所述横滚振型的振频，得到所述减摆器的有效质量、有效刚度和有效阻尼；将得到的减摆器的等效阻尼、有效质量、有效刚度和有效阻尼代入旋翼-机体系统耦合方程，得到滑行或滑跑状态下地面共振特性。根据本申请的至少一个实施例，获取所述直升机的轮胎的侧向刚度系数，包括：按下式计算所述直升机的轮胎的侧向刚度系数：其中，Ky为轮胎侧向刚度系数，Ky0为轮胎不滚动时的侧向刚度系数，V为直升机滑行或滑跑的速度，r为轮胎未压缩时的半径，p为直升机机体振动频率。根据本申请的至少一个实施例，获取所述直升机的减摆器的等效阻尼，包括：根据单个频率振动时的阻尼特性试验结果，用函数拟合方法识别阻尼力随速度变化的函数表达式，构造阻尼力函数；计算由两个简谐运动合成情况下液压阻力在低频振动一个周期中锁座的功，获取所述减摆器的等效阻尼。根据本申请的至少一个实施例，根据单个频率振动时的阻尼特性试验结果，用函数拟合方法识别阻尼力随速度变化的函数表达式，构造阻尼力函数，包括：所述阻尼力函数如下：其中，V为活塞的相对运动速度，c为常数，ξ0为减摆器摆振幅值，bt为与振幅ξ0有关的参数，V＞0时取正号，V＜0时取负号；以单频ωξ作简谐振动u＝u0sin(ωξt)在一个周期中阻尼力所做的功为：根据同样的振幅u0和频率ωξ，从实验结果中可得出相应的阻尼力所做功的实验值AE，固定振幅u0，改变频率ωξi，并计算误差：当取极小值时，即可识别出c和bt(ξ0)的值。根据本申请的至少一个实施例，计算由两个简谐运动合成情况下液压阻力在低频振动一个周期中所做的功，获取所述减摆器的等效阻尼，包括：液压阻力在低频振动一个周期中所做的功如下：式中，V0＝lξ0ωξcos(ωξt)为减摆器活塞杆以低频ωξ振动的速度，振幅值为ξ0，V1＝lξ1Ωcos(Ωt)为减摆器活塞杆以高频Ω振动的速度，振动幅值为ξ1，l为垂直铰中心到液压减摆器活塞杆轴线的距离；因此，减摆器的等效阻尼：根据本申请的至少一个实施例，所述横滚振型的振频为：其中，Ky＝Ky1+Ky2，Ke＝Ky1e1+Ky2e2，Ky1＝2n′KYM，Ky2＝n′m′KYN，e1＝ZM，e2＝ZN，ZM为主起落架刚心到重心的Z向距离，ZN为尾起落架刚心到重心的Z向距离，m′为尾起落架机轮稳定距对尾轮侧向刚度的影响系数，n′为滑行或滑跑时，机轮侧向刚度降低系数，KYM为主起Y向刚度系数，KZM为主起Z向刚度系数，KYN为尾起Y向刚度系数，yM为主起轮距之半。根据本申请的至少一个实施例，根据所述横滚振型的振频，得到所述减摆器的有效质量、有效刚度和有效阻尼，包括：所述横滚振型为：其中，m为全机质量；通过所述横滚振型和所述横滚振型的振频，得到平面模型浆毂中心处的有效质量：有效刚度：有效阻尼：其中，Ix为全机绕x轴质量转动惯量，ZH为全机重心到桨毂中心的Z向距离，CYM为主起Y向阻尼系数，CZM为主起Z向阻尼系数，CYN为尾起Y向阻尼系数。根据本申请的至少一个实施例，将得到的减摆器的等效阻尼、有效质量、有效刚度和有效阻尼代入旋翼-机体系统耦合方程，得到滑行或滑跑状态下地面共振特性，包括：所述旋翼-机体系统耦合方程为：其中，对时间的一阶导数对时间的二阶导数直升机的总重量M＝mf+kmye，Zf(t)＝x(t)+iy(t)，mye为单片桨叶质量，Zf(t)为描述桨毂中心运动的复数坐标，为Zf(t)的共轭复数，ξ1(t)描述旋翼总重心运动的复数坐标，Icj为一片桨叶绕垂直铰的质量转动惯量，lcj为桨毂中心到垂直铰的距离，lG为桨叶重心到垂直铰的距离，Scj为一片桨叶绕摆振铰的质量静矩，Kcj为一片桨叶绕垂直铰摆振的角刚度系数，mx，my为全机在桨毂中心处x、y方向上的有效质量，Kx，Ky为全机在桨毂中心处x，y方向上的有效刚度，Cx，Cy为全机在桨毂中心处x，y方向上的有效阻尼。本申请实施例提供的直升机滑行或滑跑状态下地面共振计算方法中，通过对轮胎的侧向刚度系数计算、两个激励频率作用下的减摆器阻尼特性计算、全机振频、振型计算能够真实模拟直升机滑行滑跑状态下的地面共振特性。附图说明图1是本申请实施例提供的平面动力模型简图；图2是本申请实施例提供的某型直升机滑跑时最小重量和最大重量地面共振临界稳定转速与滑跑速度关系图；图3是本申请实施例提供的不平坦土地大速度着陆滑跑尾部典型垂向过载响应时间历程；图4是本申请实施例提供的不平坦土地着陆滑跑全程减摆器位移时间历程。具体实施方式下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释相关申请，而非对该申请的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本申请相关的部分。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。地面共振本质上是带有弹性支座的旋翼在旋转平面力的动力稳定性问题，可采用平面动力模型进行计算，如图1所示，图1是本申请实施例提供的平面动力模型简图。第一方面，本申请实施例提供的直升机滑行或滑跑状态下地面共振的计算方法包括以下步骤：步骤101，获取直升机的轮胎的侧向刚度系数、减摆器的等效阻尼。在本实施例中，直升机的轮胎的侧向刚度系数通过下式计算得到：其中，Ky为轮胎滚动时的侧向刚度系数，Ky0为轮胎不滚动时的侧向刚度系数，V为直升机滑行或滑跑的速度，r为轮胎未压缩时的半径，p为直升机机体振动频率。在本实施例的一些可选的实施例中，在进行滑行或滑跑状态下地面共振的计算时，首先在0～1之间给定一组不同的系数n′，然后按照公式：Ky＝n′Ky0计算直升机的轮胎在滚动时的侧向刚度系数。在本实施例中，当直升机在地面滑行或滑跑时，桨叶摆振运动由两个不同频率(ωξ和Ω)的简谐运动组成，在这种合成运动情况下，需要计算低频ωξ振动一周期内阻尼力做功较为困难。因此，可以在阻尼器单频阻尼特性试验的基础上，来进行计算。在本实施例的一些可选的实施例中，根据单个频率振动时的阻尼特性试验结果，用函数拟合方法识别阻尼力随速度变化的函数表达式，构造如下形式的阻尼力函数：其中，V为活塞的相对运动速度，c为常数，ξ0为减摆器摆振幅值，bt为与振幅ξ0有关的参数。可选地，V＞0时取正号，V＜0时取负号。以单频ωξ作简谐振动u＝u0si本文档来自技高网...

【技术保护点】
1.一种直升机滑行或滑跑状态下地面共振的分析方法，其特征在于，包括：获取所述直升机的轮胎的侧向刚度系数、减摆器的等效阻尼；获取横滚振型的振频，并根据所述横滚振型的振频，得到所述减摆器的有效质量、有效刚度和有效阻尼；将得到的减摆器的等效阻尼、有效质量、有效刚度和有效阻尼代入旋翼‑机体系统耦合方程，得到滑行或滑跑状态下地面共振特性。

【技术特征摘要】
1.一种直升机滑行或滑跑状态下地面共振的分析方法，其特征在于，包括：获取所述直升机的轮胎的侧向刚度系数、减摆器的等效阻尼；获取横滚振型的振频，并根据所述横滚振型的振频，得到所述减摆器的有效质量、有效刚度和有效阻尼；将得到的减摆器的等效阻尼、有效质量、有效刚度和有效阻尼代入旋翼-机体系统耦合方程，得到滑行或滑跑状态下地面共振特性。2.根据权利要求1所述的直升机滑行或滑跑状态下地面共振的分析方法，其特征在于，获取所述直升机的轮胎的侧向刚度系数，包括：按下式计算所述直升机的轮胎的侧向刚度系数：其中，Ky为轮胎侧向刚度系数，Ky0为轮胎不滚动时的侧向刚度系数，V为直升机滑行或滑跑的速度，r为轮胎未压缩时的半径，p为直升机机体振动频率。3.根据权利要求1所述的直升机滑行或滑跑状态下地面共振的分析方法，其特征在于，获取所述直升机的减摆器的等效阻尼，包括：根据单个频率振动时的阻尼特性试验结果，用函数拟合方法识别阻尼力随速度变化的函数表达式，构造阻尼力函数；计算由两个简谐运动合成情况下液压阻力在低频振动一个周期中锁座的功，获取所述减摆器的等效阻尼。4.根据权利要求3所述的直升机滑行或滑跑状态下的地面共振的分析方法，其特征在于，根据单个频率振动时的阻尼特性试验结果，用函数拟合方法识别阻尼力随速度变化的函数表达式，构造阻尼力函数，包括：所述阻尼力函数如下：其中，V为活塞的相对运动速度，c为常数，ξ0为减摆器摆振幅值，bt为与振幅ξ0有关的参数，V＞0时取正号，V＜0时取负号；以单频ωξ作简谐振动u＝u0sin(ωξt)在一个周期中阻尼力所做的功为：根据同样的振幅u0和频率ωξ，从实验结果中可得出相应的阻尼力所做功的实验值AE，固定振幅u0，改变频率ωξi，并计算误差：当取极小值时，即可识别出c和bt(ξ0)的值。5.根据权利要求4所述的直升机滑行或滑跑状态下的地面共振的分析方法，其特征在于，计算由两个简谐运动合成情况下液压阻力在低频振动一个周期中所做的功，获取所述减摆器的等效阻尼，包括：液压阻力在低频振动一个周期中所做的功如下：式中，V0＝lξ0ωξcos(ωξt)为减摆器活塞杆以低频ωξ振动的速度，振幅值为ξ0，V1＝lξ1Ωcos(...

【专利技术属性】
技术研发人员：周华，程金送，凌爱民，
申请(专利权)人：中国直升机设计研究所，
类型：发明
国别省市：江西,36