一种旋翼系统装试验台耦合固有频率计算方法技术方案

本发明专利技术属于直升机旋翼动力学领域，具体涉及一种旋翼系统装试验台耦合固有频率计算方法。所述方法包括：在试验台坐标系中描述试验台和旋翼桨叶的整体运动，计及试验台运动和桨叶运动的耦合作用；根据所述耦合作用得到桨叶上任意一点的矢径和速度；利用得到的矢径和速度，计算桨叶的应变能和动能；通过实验测试得到试验台的应变能和动能；利用桨叶、试验台的应变能和动能，采用Hamilton

【技术实现步骤摘要】
一种旋翼系统装试验台耦合固有频率计算方法


[0001]本专利技术属于直升机旋翼动力学领域，具体涉及一种旋翼系统装试验台耦合固有频率计算方法。

技术介绍

[0002]对于装在试验台上的直升机旋翼,在以往旋翼固有特性的计算中,没有考虑试验台与旋翼耦合对旋翼固有频率的影响,认为旋翼在旋转中心是固支的,这样使得计算结果与实际情况有比较大的差别。
[0003]孤立旋翼理论模型是基于桨毂中心坐标描述桨叶上任一点在惯性坐标系中的位形，并采用Hamilton
’
s原理推导系统的动力学方程。
[0004][0005]δU，系统应变能；δT，系统动能；δW，系统外力虚功。将其写为矩阵模式，即为孤立旋翼总体矩阵。
[0006][0007]上述方程中没有考虑试验台对旋翼系统的影响，特别是试验台旋翼轴刚度和惯量等参数对耦合系统的作用在动力学方程中无处体现，使得仿真计算值与实际测量值有很大差别，无法对旋翼/试验台耦合系统进行固有频率参数影响分析，得不到试验台各结构参数对旋翼固有频率的影响规律，具有较大的局限性。

技术实现思路

[0008]专利技术目的：提供一种旋翼系统装试验台耦合固有频率计算方法，在孤立旋翼分析方法的基础上，考虑旋翼和试验台的耦合因素，将试验台结构模态参数与旋翼结构参数进行拼装得到耦合系统动力学方程，分析试验台对旋翼固有频率的影响规律。
[0009]专利技术技术方案：提供一种旋翼系统装试验台耦合固有频率计算方法，所述方法包括：在试验台坐标系中描述试验台和旋翼桨叶的整体运动，计及试验台运动和桨叶运动的耦合作用；根据所述耦合作用得到桨叶上任意一点的矢径和速度；
[0010]利用得到的矢径和速度，计算桨叶的应变能和动能；通过实验测试得到试验台的应变能和动能；
[0011]利用桨叶、试验台的应变能和动能，采用Hamilton
’
s原理推导得到旋翼/试验台耦合系统的运动学方程；求解所述动力学方程旋翼/试验台耦合系统的固有特性。
[0012]进一步地，通过共振法获得试验台的应变能和动能，具体为：使用激振设备对试验台某点施加一个激励，再通过试验台上安装的测试仪器得出各测试点的加速度、速度和位移；根据得到的各测试点的加速度、速度和位移计算得出试验台的应变能和动能。
[0013]进一步地，在试验台坐标系中描述试验台和旋翼桨叶的整体运动包括：首先，定义
试验台/旋翼耦合系统坐标系系统，然后根据所述坐标系系统进行坐标转换得到旋翼桨叶的整体运动。
[0014]进一步地，试验台/旋翼耦合系统坐标系系统的建立包括：首先，以试验台与地面连接处某点为原点建立地面固定惯性坐标系；然后，建立试验台坐标系、桨毂不旋转坐标系、桨毂旋转坐标系、桨叶未变形坐标系和桨叶变形坐标系，及各坐标系之间的转换矩阵。
[0015]进一步地，桨叶的应变能和动能的计算方法为：根据桨叶上任意一点的矢径和速度，首先计算出桨叶上任意一点矢径和速度在地面惯性坐标系中的表达式，然后推导出桨叶上任意一点的应变能和动能变分，再通过积分求得整片桨叶的应变能δU和动能δT表达式。
[0016]进一步地，利用桨叶、试验台的应变能和动能，根据Hamilton
’
s原理，利用有限单元离散方法，推导出桨叶单元矩阵、桨叶/试验台耦合单元矩阵、试验台/桨叶耦合单元矩阵和试验台单元矩阵，并将上述单元矩阵进行组合得到对应最终的桨叶、桨叶/试验台、试验台/桨叶和试验台各自的总体阵。
[0017]进一步地，将推导得到的桨叶总体阵进行多桨叶坐标变换后得到旋翼整体矩阵，然后将旋翼整体矩阵、桨叶/试验台的总体阵、试验台/桨叶的总体阵和试验台的总体阵进行组集得到旋翼/试验台耦合系统的运动学方程。
[0018]进一步地，求解所述动力学方程旋翼/试验台耦合系统的固有特性包括：对旋翼/试验台耦合系统运动学方程中的结构运动学方程进行特征值求解，得到旋翼/试验台耦合系统的固有频率和振型。
[0019]本专利技术的技术效果：本专利基于孤立旋翼固有特性分析方法，考虑试验台旋翼进行耦合的影响，可以分析试验台各结构参数对旋翼固有频率的影响。相对传统的未考虑试验台的固有特性分析方法，本专利的试验精度更高。
附图说明
[0020]图1旋翼摆振集合型频率随旋翼轴刚度变化趋势图；
[0021]图2旋翼摆振集合型频率随旋翼轴转动惯量变化趋势图。
具体实施方式
[0022]实施例1
[0023]本实施例，提供一种计及试验台参数对旋翼影响的旋翼固有频率计算方法，解决旋翼固有频率计算与试验结果误差较大的问题。
[0024]在试验台坐标系中描述试验台和旋翼桨叶的整体运动，考虑试验台运动和桨叶运动的相互影响作用，计算桨叶和试验台的应变能和动能，并采用Hamilton
’
s原理推导旋翼/试验台耦合系统的动力学方程，进而求解方程得到旋翼/试验台耦合系统固有特性。具体包括以下步骤：
[0025]步骤1：通过共振法法获得试验台的固有特性，使用激振设备对试验台某点施加一个激励，再通过试验台上安装的加速度计等测试仪器得出各测试点的加速度、速度和位移，进而计算得出试验台的固有特性等参数(固有频率和振型等)。
[0026]步骤2：定义试验台/旋翼耦合系统坐标系，以试验台与地面连接处某点为原点建
立地面固定惯性坐标系，接着建立试验台坐标系，桨毂不旋转坐标系，桨毂旋转坐标系，桨叶未变形坐标系和桨叶变形坐标系，同时得出各坐标系的转换矩阵。
[0027]步骤3：计算桨叶上任意一点矢径和速度在地面惯性坐标系中的表达式，由此推导桨叶上任意一点应变能和动能变分，通过积分求得整片桨叶的应变能δU和动能δT表达式
[0028]步骤4：根据Hamilton
’
s原理，从应变能和动能入手，以桨叶和试验台两个角度分析，利用有限单元离散方法，推导桨叶单元矩阵、桨叶/试验台耦合单元矩阵、试验台/桨叶耦合单元矩阵和试验台单元矩阵。
[0029]步骤5：通过将局部单元自由度转换到总体自由度的形式，将单元矩阵添加到相应的位置后，考虑桨根处连接方式的影响，组成对应最终的桨叶、桨叶/试验台、试验台/桨叶和试验台总体阵。
[0030]步骤6：推导旋翼/试验台耦合系统的动力学方程，先对桨叶总体矩阵进行多桨叶坐标变换后得到旋翼整体矩阵，然后组集旋翼/试验台耦合系统各部分总体阵可以得到旋翼/试验台耦合运动方程如下。
[0031][0032][0033]将其写为矩阵模式，即为耦合系统总体矩阵。
[0034][0035]上式中：M
b
、C
b
、K
b
、F
b
分别为桨叶总质量、阻尼、刚度矩阵和载荷矢量；
[0036]M
FF
、C
FF
、K
FF
、F
F
分别为试验台总质量、阻尼、刚度矩阵和本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种旋翼系统装试验台耦合固有频率计算方法，其特征在于，所述方法包括：在试验台坐标系中描述试验台和旋翼桨叶的整体运动，计及试验台运动和桨叶运动的耦合作用；根据所述耦合作用得到桨叶上任意一点的矢径和速度；利用得到的矢径和速度，计算桨叶的应变能和动能；通过实验测试得到试验台的应变能和动能；利用桨叶、试验台的应变能和动能，采用Hamilton
’
s原理推导得到旋翼/试验台耦合系统的运动学方程；求解所述动力学方程旋翼/试验台耦合系统的固有特性。2.根据权利要求1所述的固有频率计算方法，其特征在于，通过共振法获得试验台的应变能和动能，具体为：使用激振设备对试验台某点施加一个激励，再通过试验台上安装的测试仪器得出各测试点的加速度、速度和位移；根据得到的各测试点的加速度、速度和位移计算得出试验台的应变能和动能。3.根据权利要求1所述的固有频率计算方法，其特征在于，在试验台坐标系中描述试验台和旋翼桨叶的整体运动包括：首先，定义试验台/旋翼耦合系统坐标系系统，然后根据所述坐标系系统进行坐标转换得到旋翼桨叶的整体运动。4.根据权利要求1所述的固有频率计算方法，其特征在于，试验台/旋翼耦合系统坐标系系统的建立包括：首先，以试验台与地面连接处某点为原点建立地面固定惯性坐标系；然后，建立试验台坐标系、桨毂不旋转坐标系、桨毂旋转坐标系...

【专利技术属性】
技术研发人员：程毅，赵金瑞，余智豪，周云，程起有，魏武雷，
申请(专利权)人：中国直升机设计研究所，
类型：发明
国别省市：