CN113955142A - 一种旋翼1ω谐波振动实时调整控制方法 - Google Patents
一种旋翼1ω谐波振动实时调整控制方法 Download PDFInfo
- Publication number
- CN113955142A CN113955142A CN202111176139.8A CN202111176139A CN113955142A CN 113955142 A CN113955142 A CN 113955142A CN 202111176139 A CN202111176139 A CN 202111176139A CN 113955142 A CN113955142 A CN 113955142A
- Authority
- CN
- China
- Prior art keywords
- vibration
- pull rod
- omega
- data
- disturbance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/40—Maintaining or repairing aircraft
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C27/00—Rotorcraft; Rotors peculiar thereto
- B64C27/51—Damping of blade movements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- Control Of Turbines (AREA)
Abstract
本申请属于旋翼振动控制技术领域,公开了一种旋翼1Ω谐波振动实时调整控制方法,包括以下步骤:步骤一:采集数据;步骤二:对采集的数据进行处理;得到振动数据的1Ω振动幅值;步骤三:重复N次步骤一、二,得到N组1Ω振动幅值、N组拉杆扰动位移和N组桨叶桨尖高度;步骤四:根据N组的1Ω振动幅值和对应的N组拉杆扰动位移通过最小二乘辨识算法得到拉杆扰动位移和1Ω振动幅值之间的转换矩阵T;步骤五:根据N组桨叶桨尖高度和对应的N组拉杆扰动位移通过最小二乘辨识算法得到拉杆扰动位移和桨叶桨尖高度之间的转换矩阵A;步骤六:建立拉杆扰动位移和1Ω振动幅值的二次型性能指标函数;求解所述二次型性能指标函数。
Description
技术领域
本发明属于旋翼振动控制技术领域,公开了一种旋翼1Ω谐波振动实时调整控制方法。
背景技术
旋翼1Ω谐波振动调整是直升机的一个常规维护项目,目前调整采用人工调整方式,通过改变旋翼变距拉杆长度或者在桨毂支臂上加减配重片进行调整。实际操作时,机械师需要根据每次飞行后振动数据,人工计算拉杆调整量去进行地面手动调整,然后再飞行验证。如果效果不行,需重复在地面进行拉杆调整后飞行验证。
当前人工调整的方法,无法对任一飞行状态进行旋翼1Ω谐波振动实时自动调整,且人工调整方式,因人而异,不同的机械师可能给出不同的调整量,使整个旋翼1Ω谐波振动调整过程耗费大量时间、人力、物力、财力。因此可以通过旋翼智能变距拉杆自动调整的方式,在直升机飞行过程实时对旋翼1Ω谐波振动进行自动调整,能节省调整过程的人力、物力、财力。本文提出了一种基于旋翼智能变距拉杆自动调整方式的旋翼1Ω谐波振动实时调整控制方法,能够应用在旋翼智能拉杆上进行直升机振动主动控制,大大提升了直升机的使用维护性和经济性。
发明内容
本发明提出一种旋翼1Ω谐波振动实时调整控制方法,包括以下步骤:
步骤一:在飞行过程中采集旋翼转速数据、主减速器振动数据、拉杆扰动位移数据、各片桨叶桨尖高度数据;
步骤二:对采集的旋翼转速数据和主减速器振动数据进行处理;得到振动数据的1Ω振动幅值;
步骤三:重复N次步骤一、二,得到N组1Ω振动幅值、N组拉杆扰动位移和N组桨叶桨尖高度;
步骤四:根据N组的1Ω振动幅值和对应的N组拉杆扰动位移通过最小二乘辨识算法得到拉杆扰动位移和1Ω振动幅值之间的转换矩阵T;
步骤五:根据N组桨叶桨尖高度和对应的N组拉杆扰动位移通过最小二乘辨识算法得到拉杆扰动位移和桨叶桨尖高度之间的转换矩阵A;
步骤六:建立拉杆扰动位移和1Ω振动幅值的二次型性能指标函数;求解所述二次型性能指标函数得到的拉杆扰动位移作为最优拉杆扰动位移;并输出所述最优拉杆扰动位移给智能拉杆进行旋翼1Ω谐波振动实时调整。
进一步,所述步骤二中,数据处理过程如下:对主减速器振动数据进行傅里叶变化,提取1Ω振动频率对应的振动幅值。
进一步,所述旋翼转速数据为周期数据;所述步骤二中,数据处理过程还包括:截取至少一个完整周期旋翼转速数据所对应的主减速器振动数据进行傅里叶变化。
进一步,所述步骤四中,还包括:根据拉杆扰动位移和1Ω振动幅值之间的转换矩阵T建立拉杆扰动位移和1Ω振动幅值之间的线性关系:Z=Z0+T·u;
其中,Z为预估的主减速器振动值,u为拉杆扰动位移,Z0为主减速器初始振动值。
进一步,所述步骤四中,还包括:根据拉杆扰动位移和桨叶桨尖高度的转换矩阵A建立拉杆扰动位移和桨叶桨尖高度之间线性关系:G=G0+A·u;
其中,G为预估的桨叶桨尖高度,u为拉杆扰动位移,G0为桨叶桨尖高度初始值。
进一步,所述步骤六中,拉杆扰动位移和1Ω振动幅值的二次型性能指标函数把桨叶桨尖高度最大允许差值、拉杆扰动位移允许的最小和最大位移作为约束条件,二次型性能指标函数表示如下:
min.J=ZT·WZ·Z+uT·Wu·u
s.t.Z=Z0+T·u
A·u=G-G0≤b
umin≤u≤umax
其中,Wz,Wu为权矩阵,取单位对角矩阵,
b为桨叶桨尖高度最大允许差值,
uminumax为拉杆扰动位移允许的最小和最大位移。
进一步,所述步骤六中,二次型性能指标函数求解过程如下:将二次型性能指标函数需转化为标准二次规划问题;
标准二次规划问题数据表达式如下:
H=2·(TT·WZ·T+Wu)
q=(2·Z0 T·WZ·T)T
s.t.A·u=G-G0≤b
umin≤u≤umax
进一步,所述步骤六中,采用内点法求解标准二次规划问题。
有益效果:
本专利提出的旋翼1Ω谐波振动实时调整控制方法,能够应用在旋翼智能拉杆上,在直升机飞行过程实时对旋翼1Ω谐波振动进行自动调整,大大提升了直升机的使用维护性和经济性。
附图说明
图1是旋翼1Ω谐波振动实时调整控制方法流程示意图。
具体实施方式
下面通过具体实施方式对本发明做进一步说明:
本发明的一种旋翼1Ω谐波振动实时调整控制方法主要包含数据采集、数据处理与存储、数据辨识和数据求解输出四大步骤。
1.数据采集:为实现对旋翼1Ω谐波振动的实时调整,首先需要在飞行过程中获取与控制量和控制目标相关的原始数据,采集的数据包括:各片桨叶桨尖高度、主减速器振动、拉杆扰动位移信号、旋翼转速数据,;
2.数据处理:本方法主要是调整旋翼1Ω谐波振动,因此需要从采集的数据中提取振动数据的1Ω振动幅值。为确保提取的1Ω振动幅值的准确性,需结合转速信号,截取转速整周期时间段相对应的振动数据段,然后对截取数据段采用傅里叶变换方法提取1Ω振动幅值;
3.数据存储:重复N次1)、2)过程,并把N次的1Ω振动幅值、拉杆扰动位移、桨叶桨尖高度进行存储,获得N组1Ω振动幅值、拉杆扰动位移、桨叶桨尖高度数据,以便下一步通过对N组数据的辨识获得拉杆扰动位移与1Ω振动幅值以及桨叶桨尖高度的线性关系;
4.数据辨识:采用最小二乘法辨识方法,辨识得到拉杆扰动位移与1Ω振动幅值转换矩阵T以及拉杆扰动位移与桨叶桨尖高度的转换矩阵A,从而可得到拉杆扰动位移与1Ω振动幅值以及桨叶桨尖高度的线性关系。
转换矩阵的关系表达式如下:
Z=Z0+T·u
G=G0+A·u
其中,Z为预估的1Ω振动幅值,G为预估的桨叶桨尖高度,T、A为转换矩阵,u为拉杆扰动位移,Z0为初始振动值,G0为初始桨叶桨尖高度。
5.数据求解与输出
在直升机旋翼1Ω谐波振动调整的要求中,除了要求1Ω振动幅值降低到允许范围内,还需要满足各片桨叶的桨尖高度互差不超过限制值的要求,以及智能拉杆的扰动位移需在允许作动范围内。
因此首先需要以桨叶桨尖高度最大允许差值、拉杆扰动位移允许的最小和最大位移作为约束条件,建立拉杆扰动位移和1Ω振动幅值的二次型性能指标函数J。
二次型性能指标函数如下:
min.J=ZT·WZ·Z+uT·Wu·u
s.t.Z=Z0+T·u
A·u=G-G0≤b
umin≤u≤umax
其中,Wz,Wu为权矩阵,取单位对角矩阵,
b为桨叶桨尖高度最大允许差值,
uminumax为拉杆扰动位移允许的最小和最大位移
然后将性能指标函数转化为标准二次规划问题,并使用内点法等优化算法求解带有约束条件的标准二次规划问题的最小值,即可求解得到最优控制量。
标准二次规划问题数学表达式如下:
H=2·(TT·WZ·T+Wu)
q=(2·Z0 T·WZ·T)T
s.t.A·u=G-G0≤b
umin≤u≤umax
最后将求解得到的最优控制量(拉杆扰动位移)输出给智能拉杆,从而实现对旋翼1Ω谐波振动的调整控制。
以上所述,仅为本发明的具体实施例,对本发明进行详细描述,未详尽部分为常规技术。但本发明的保护范围不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。本发明的保护范围应以所述权利要求的保护范围为准。
Claims (8)
1.一种旋翼1Ω谐波振动实时调整控制方法,其特征在于:所述方法包括以下步骤:
步骤一:在飞行过程中采集旋翼转速数据、主减速器振动数据、拉杆扰动位移数据、各片桨叶桨尖高度数据;
步骤二:对采集的旋翼转速数据和主减速器振动数据进行处理;得到振动数据的1Ω振动幅值;
步骤三:重复N次步骤一、二,得到N组1Ω振动幅值、N组拉杆扰动位移和N组桨叶桨尖高度;
步骤四:根据N组的1Ω振动幅值和对应的N组拉杆扰动位移通过最小二乘辨识算法得到拉杆扰动位移和1Ω振动幅值之间的转换矩阵T;
步骤五:根据N组桨叶桨尖高度和对应的N组拉杆扰动位移通过最小二乘辨识算法得到拉杆扰动位移和桨叶桨尖高度之间的转换矩阵A;
步骤六:建立拉杆扰动位移和1Ω振动幅值的二次型性能指标函数;求解所述二次型性能指标函数得到的拉杆扰动位移作为最优拉杆扰动位移;并输出所述最优拉杆扰动位移给智能拉杆进行旋翼1Ω谐波振动实时调整。
2.根据权利要求1所述的一种旋翼1Ω谐波振动实时调整控制方法,其特征在于:所述步骤二中,数据处理过程如下:对主减速器振动数据进行傅里叶变化,提取1Ω振动频率对应的振动幅值。
3.根据权利要求2所述的一种旋翼1Ω谐波振动实时调整控制方法,其特征在于:所述旋翼转速数据为周期数据;所述步骤二中,数据处理过程还包括:截取至少一个完整周期旋翼转速数据所对应的主减速器振动数据进行傅里叶变化。
4.根据权利要求3所述的一种旋翼1Ω谐波振动实时调整控制方法,其特征在于:所述步骤四中,还包括:根据拉杆扰动位移和1Ω振动幅值之间的转换矩阵T建立拉杆扰动位移和1Ω振动幅值之间的线性关系:Z=Z0+T·u;
其中,Z为预估的主减速器振动值,u为拉杆扰动位移,Z0为主减速器初始振动值。
5.根据权利要求4所述的一种旋翼1Ω谐波振动实时调整控制方法,其特征在于:所述步骤四中,还包括:根据拉杆扰动位移和桨叶桨尖高度的转换矩阵A建立拉杆扰动位移和桨叶桨尖高度之间线性关系:G=G0+A·u;
其中,G为预估的桨叶桨尖高度,u为拉杆扰动位移,G0为桨叶桨尖高度初始值。
6.根据权利要求5所述的一种旋翼1Ω谐波振动实时调整控制方法,其特征在于:所述步骤六中,拉杆扰动位移和1Ω振动幅值的二次型性能指标函数把桨叶桨尖高度最大允许差值、拉杆扰动位移允许的最小和最大位移作为约束条件,二次型性能指标函数表示如下:
min.J=ZT·Wz·Z+uT·Wu·u
s.t.Z=Z0+T·u
A·u=G-G0≤b
umin≤u≤umax
其中,Wz,Wu为权矩阵,取单位对角矩阵,
b为桨叶桨尖高度最大允许差值,
umin umax为拉杆扰动位移允许的最小和最大位移。
8.根据权利要求7所述的一种旋翼1Ω谐波振动实时调整控制方法,其特征在于:所述步骤六中,采用内点法求解标准二次规划问题。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111176139.8A CN113955142B (zh) | 2021-10-09 | 2021-10-09 | 一种旋翼1ω谐波振动实时调整控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111176139.8A CN113955142B (zh) | 2021-10-09 | 2021-10-09 | 一种旋翼1ω谐波振动实时调整控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113955142A true CN113955142A (zh) | 2022-01-21 |
CN113955142B CN113955142B (zh) | 2023-09-05 |
Family
ID=79463201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111176139.8A Active CN113955142B (zh) | 2021-10-09 | 2021-10-09 | 一种旋翼1ω谐波振动实时调整控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113955142B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GR3023986T3 (en) * | 1992-12-11 | 1997-10-31 | Dynamic Eng Inc | System for controlling higher harmonic vibrations in helicopter rotor blades. |
CN106814609A (zh) * | 2017-01-06 | 2017-06-09 | 西安交通大学 | 一种频谱塑形主动控制方法及主动控制系统 |
CN208544421U (zh) * | 2018-06-14 | 2019-02-26 | 南京航空航天大学 | 一种独立变桨距旋翼系统 |
CN112731814A (zh) * | 2021-01-14 | 2021-04-30 | 南京航空航天大学 | 基于自适应谐波识别频响修正的直升机振动主动控制方法 |
CN112859589A (zh) * | 2021-01-14 | 2021-05-28 | 南京航空航天大学 | 旋翼变转速直升机振动主动控制的混合控制方法 |
-
2021
- 2021-10-09 CN CN202111176139.8A patent/CN113955142B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GR3023986T3 (en) * | 1992-12-11 | 1997-10-31 | Dynamic Eng Inc | System for controlling higher harmonic vibrations in helicopter rotor blades. |
CN106814609A (zh) * | 2017-01-06 | 2017-06-09 | 西安交通大学 | 一种频谱塑形主动控制方法及主动控制系统 |
CN208544421U (zh) * | 2018-06-14 | 2019-02-26 | 南京航空航天大学 | 一种独立变桨距旋翼系统 |
CN112731814A (zh) * | 2021-01-14 | 2021-04-30 | 南京航空航天大学 | 基于自适应谐波识别频响修正的直升机振动主动控制方法 |
CN112859589A (zh) * | 2021-01-14 | 2021-05-28 | 南京航空航天大学 | 旋翼变转速直升机振动主动控制的混合控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113955142B (zh) | 2023-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101400892B (zh) | 用于减小受到风轮面不对称加载的风力涡轮机的部件的疲劳负载的方法与控制系统 | |
EP3791060B1 (en) | Wind turbine control method | |
Rašuo et al. | Harmonization of new wind turbine rotor blades development process: A review | |
CN112412700A (zh) | 一种风电机组塔架振动控制系统 | |
CN106795857A (zh) | 与确定风力涡轮机中的转子失衡有关的改善 | |
Molland et al. | Wind tunnel investigation of the influence of propeller loading on ship rudder performance | |
CN110761945B (zh) | 风力发电机组的叶片失速控制方法及装置 | |
CN110454329B (zh) | 一种风电机组桨距角控制方法 | |
DE102010017777A1 (de) | Verfahren und Vorrichtung zum Steuern der Umfangsgeschwindigkeit an Rotorflügelspitzen | |
WO2020011323A1 (en) | Method and system for controlling a wind turbine to reduce nacelle vibration | |
CN107782443B (zh) | 一种风力发电机叶片固有频率自动提取方法 | |
CN110714881B (zh) | 一种风力机性能预测控制方法及装置 | |
CN113657662A (zh) | 一种基于数据融合的降尺度风电功率预测方法 | |
CN111878308A (zh) | 风力机集群预测群控方法 | |
CN113283024A (zh) | 一种基于有限元建模的拓扑优化和参数优化的风机轮毂结构优化方法 | |
CN107480324B (zh) | 一种风电机组叶片状态监测系统构建方法 | |
CN113955142A (zh) | 一种旋翼1ω谐波振动实时调整控制方法 | |
DE102010027229A1 (de) | Verfahren und Vorrichtung zur Bereitstellung eines Abstellwinkel-Korrektursignals für ein vorbestimmtes Rotorblatt eier Windkraftanlage | |
CN110457800B (zh) | 考虑机械惯性的水平轴风机风速出力折算方法 | |
Arcidiacono | Theoretical performance of helicopters having second and higher harmonic feathering control | |
Furlanetto et al. | Design optimization of tapered steel wind turbine towers by QPSO algorithm | |
CN111622893B (zh) | 非对称翼型垂直轴风力机的变桨控制方法及系统 | |
Schmaus et al. | Experimental Study of Rotor Performance in Deep Ground Effect with Application to a Human-Powered Helicopter | |
EP3536949B1 (en) | Wind farm, and operation method, control device, and operation control program for wind farm | |
Xu et al. | Development and application of a dynamic stall model for rotating wind turbine blades |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |