CN112052529A - 一种提高大前进比旋翼配平收敛性的计算方法 - Google Patents
一种提高大前进比旋翼配平收敛性的计算方法 Download PDFInfo
- Publication number
- CN112052529A CN112052529A CN202011028637.3A CN202011028637A CN112052529A CN 112052529 A CN112052529 A CN 112052529A CN 202011028637 A CN202011028637 A CN 202011028637A CN 112052529 A CN112052529 A CN 112052529A
- Authority
- CN
- China
- Prior art keywords
- iteration
- convergence
- value
- rotor
- trim
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/20—Design reuse, reusability analysis or reusability optimisation
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Testing Of Balance (AREA)
Abstract
本发明属于直升机旋翼设计与理论建模技术,公开了一种提高大前进比旋翼配平收敛性的计算方法,采用逐步逼近法方法和残差监控方法改进旋翼配平计算方法,区别于传统的配平迭代计算方法,一种提高大前进比旋翼配平收敛性的计算方法对前进比参数进行加密划分,逐步逼近法计算和实时监控迭代循环中的残差变化,该方法有效地拓展了旋翼配平的计算能力且保证配平求解精度,在大前进比状态中有效避免了配平迭代中发散问题。
Description
技术领域
本发明属于直升机旋翼设计与理论建模技术,涉及一种提高大前进比旋翼配平收敛性的计算方法。
背景技术
旋翼配平是直升机分析设计阶段的重要环节之一,通过设定稳定飞行状态下目标来确定直升机周期稳态特性和驾驶员操纵,在直升机气动、动力学和飞行力学分析中都有重要意义。在大前进比变转速旋翼气弹分析中,桨盘气流的不对称性加剧,气动与结构耦合效应增加,微弱的操纵变化会带来强烈的气动、结构载荷变化,使得配平计算困难性增加,难以收敛,因此要建立适用于大前进比的旋翼配平计算方法。
目前国内外针对旋翼配平主要采用牛顿迭代法。传统旋翼配平迭代中只在初始时刻计算Jaocbian矩阵,即整个过程只利用初始输入生成一次Jacobian矩阵,未考虑实际结果与初始输入的偏差对Jacobian矩阵的影响。这导致在大前进比配平中Jacobian矩阵质量差,且大前进比状态下的气弹计算复杂性增加,对初始输入的敏感性增加,都容易引发数值发生问题。
发明内容
本发明要解决的技术问题:提出一种提高大前进比旋翼配平收敛性的迭代计算方法,用于大前进比变转速旋翼气弹配平计算,解决大前进比变转速旋翼配平计算困难甚至发生的问题,能在设计阶段通过仿真分析确定直升机稳态飞行中的驾驶员操纵和旋翼特性。
本发明的技术方案:
一种提高大前进比旋翼配平收敛性的计算方法,包括以下步骤:
第一步:划分目标前进比参数;
第二步:设定配平目标,所述配平目标为桨毂六力,建立配平平衡方程组;采用手动输入或者内部计算两种方式给出第一个前进比参数下的初始操纵输入;
第三步:根据初始操纵输入和配平目标生成Jacobian矩阵;
第四步:进入迭代循环进行旋翼气弹计算;循环中采用牛顿迭代法调用Jacobian矩阵,直至迭代收敛得到最新的操纵值;
第五步:将迭代收敛时的操纵值作为下一个前进比参数下的操纵输入,重复进行第三、四步,直到完成所有前进比参数的计算,完成大前进比旋翼配平求解,得到配平目标配平后的操纵值。
进一步,所述第四步中,还包括:
设定收敛准则设定值;每一次迭代得到新的操纵值,并结合上一次迭代结果计算出残差值,迭代循环中引入实时监控残差方法,根据残差值与收敛准则设定值的关系判断迭代的收敛性。
进一步,若残差值高于收敛准则设定值,则迭代尚未收敛,进入残差监控环节;若残差值小于上一次迭代的残差值,则继续进行迭代。
进一步,若残差值大于与上一次迭代的残差值,则跳出迭代循环并根据此时的操纵值与配平目标重新生成Jacobian矩阵,重新开始新一轮的迭代循环。
进一步,若残差低于收敛准则设定值,则迭代收敛,当前前进比参数下的配平计算结果作为最新的操纵值。
进一步,所述第一步中的前进比参数划分包括:采用正弦加密将前进比参数划分为一组数值;划分结果如下:
其中,Ns为划分的前进比参数个数。
进一步,所述第一步中,还包括:建立旋翼动力学模型和旋翼载荷计算模型。
进一步,所述第四步中,还包括:根据旋翼动力学模型和旋翼载荷计算模型进行旋翼气弹计算。
进一步,所述牛顿迭代法包括:通过以下公式计算操纵值;
λn=λn-1-[J]-1·f(λn-1)
其中:λn为当前迭代得到的操纵值,λn-1为上一次迭代的操纵值,J为Jacobian矩阵,f为配平平衡方程组
本发明的有益效果:一种提高大前进比旋翼配平收敛性的迭代计算方法采用逐步逼近法,采用逐步逼近法方法和残差监控方法改进旋翼配平计算方法,区别于传统的配平迭代计算方法,一种提高大前进比旋翼配平收敛性的计算方法对前进比参数进行加密划分,逐步逼近法计算和实时监控迭代循环中的残差变化,该方法有效地拓展了旋翼配平的计算能力且保证配平求解精度,在大前进比状态中有效避免了配平迭代中发散问题。
附图说明
图1是本发明涉及的一种提高大前进比旋翼配平收敛性的迭代计算方法;
图2是本发明涉及的旋翼动力学方程建模中的弹性桨叶变形示意图
图3是本发明涉及的旋翼桨毂固定坐标系
图4是两种配平计算方法的迭代效果图。
具体实施方式
下面结合附图对本发明所涉及的提高大前进比旋翼配平收敛性的计算方法做进一步详细说明。
一种提高大前进比旋翼配平收敛性的迭代计算方法,如图1所示,包括以下步骤:
(1)建立旋翼动力学模型。建立旋翼气动力模型,考虑桨叶挥舞、摆振、扭转的非线性弹性变形以及相互之间的非线性耦合关系,准确描述弹性桨叶非线性运动,弹性桨叶坐标系如图2所示,通过坐标转换形式得到桨叶变形坐标系E与桨叶未变形坐标系B下的几何关系,包括轴向变形u,摆振变形v,和挥舞变形w。并将系统能量分为应变能U,动能T,外载荷W,最后采用Hamilton原理建立旋翼动力学方程,其中δ表示变分:
∫(δU-δT-δW)dt=0
(2)建立旋翼载荷计算模型。采用力积分方法建立桨叶剖面载荷计算模型,将剖面以外的作用力,包括气动力,惯性力和阻尼器等外部作用力。得到桨叶任意剖面的六力,包括剖面展向拉力fx,剖面弦向剪力剖面弦向剪力fy,剖面垂直剪力fz,剖面挥舞弯矩mx,剖面摆振弯矩mz,剖面扭矩my。通过从桨根到桨尖的积分计算得到桨叶桨根载荷,并通过坐标系转换,桨根载荷转换到旋翼桨毂固定坐标系中如图3所示,叠加得到桨毂六力,包括阻力Fx,侧向力Fy,拉力Fz,滚转力矩Mx,俯仰力矩My,扭矩Mz。
(3)集成旋翼动力学模型和旋翼载荷计算模型,最终得到旋翼气弹计算方法。
(4)建立旋翼配平分析计算方法。以桨毂原点为中心,采用牛顿迭代法建立桨毂六力平衡方程,即旋翼配平方程,根据初始输入条件与配平目标生成Jacobian矩阵。牛顿迭代法公式如下:
(5)建立提高大前进比旋翼配平收敛性的迭代计算方法。采用逐步逼近法方法和残差监控方法改进旋翼配平计算方法,在配平循环流程中对前进比参数进行加密划分,逐步配平,逼近设定的前进比目标值,迭代循环过程中实时监控迭代循环中的残差变化,及时自主更新Jacobian矩阵,控制迭代计算中的参数数值,提高收敛性,解决大前进比状态下配平计算发散问题。
一种提高大前进比旋翼配平收敛性的计算方法具体实施方式如下,如图1所示:
第一步:建立旋翼动力学模型和旋翼载荷计算模型。针对旋翼前进比参数进行加密划分,采用正弦加密法将前进比参数划分为一组数值,加密划分表示如下,Ns为划分的前进比参数个数。对每一个前进比参数进行配平,当前前进比参数下配平收敛后进行下一个前进比参数配平,直至配平至目标前进比参数。即逐步逼近方式进行配平求解。
第二步:设定配平目标,配平目标为桨毂六力,建立配平平衡方程组,进入配平循环,采用手动输入或者内部计算两种方式给出第一个前进比参数下的初始操纵输入,后续的配平循环采用上一个前进比参数配平操纵结果作为初始操纵输入;
第三步:根据初始操纵输入和配平目标生成Jacobian矩阵;
第四步:设置收敛准则设定值,进入迭代循环,进行旋翼气弹计算;循环中采用牛顿迭代法调用Jacobian矩阵,直至迭代收敛得到最新的操纵值;期间通过每一步的迭代得到更新的操纵值,并结合上一步迭代结果计算当前迭代步的残差值error,并引入实时监控残差方法判断迭代循环的收敛性,根据残差值error的大小以及变化会出现以下几种选择:
1)若残差高于收敛准则设定值,则迭代还未收敛,则进入残差监控环节,若残差值小于上一次迭代的残差值,则继续进行迭代。
2)若残差高于收敛准则设定值,则迭代还未收敛,则进入残差监控环节,若残差值大于上一次迭代的残差值,则跳出迭代循环并根据此时的操纵值与配平目标重新生成Jacobian矩阵,重新开始新一轮的迭代循环。
3)若残差低于收敛准则设定值,则迭代收敛,当前前进比参数下的配平计算结果作为最新的操纵值,进行下一轮配平。
第五步:将迭代收敛时的操纵值作为下一个前进比参数下的操纵输入,重复进行第三、四步,直到完成所有前进比参数的计算,完成大前进比旋翼配平求解,得到配平目标配平后的操纵值。
一种提高大前进比旋翼配平收敛性的迭代计算方法采用逐步逼近法,采用逐步逼近法方法和残差监控方法改进旋翼配平计算方法,区别于传统的配平迭代计算方法,一种提高大前进比旋翼配平收敛性的计算方法对前进比参数进行加密划分,逐步逼近法计算和实时监控迭代循环中的残差变化,该方法有效地拓展了旋翼配平的计算能力且保证配平求解精度,在大前进比状态中有效避免了配平迭代中发散问题,如图4所示。
Claims (9)
1.一种提高大前进比旋翼配平收敛性的计算方法,其特征在于:所述方法包括:
第一步:获取目标前进比参数的正常范围,并在目标前进比参数的正常范围内将目标前进比参数划分为一组数值;
第二步:设定配平目标,所述配平目标为桨毂六力,建立配平平衡方程组;采用手动输入或者内部计算两种方式给出第一个前进比参数下的初始操纵输入;
第三步:根据初始操纵输入和配平目标生成Jacobian矩阵;
第四步:进入迭代循环进行旋翼气弹计算;循环中采用牛顿迭代法调用Jacobian矩阵,直至迭代收敛得到最新的操纵值,将初始操纵输入替换为最新的操纵值;
第五步:将迭代收敛时的操纵值作为下一个前进比参数下的操纵输入,重复进行第三、四步,直到完成所有前进比参数的计算,完成大前进比旋翼配平求解,得到配平目标配平后的操纵值。
2.根据权利要求1所述的一种提高大前进比旋翼配平收敛性的计算方法,其特征在于:所述第四步中,还包括:设定收敛准则设定值;每一次迭代得到新的操纵值,并结合上一次迭代结果计算出残差值,迭代循环中引入实时监控残差方法,根据残差值与收敛准则设定值的关系判断迭代的收敛性。
3.根据权利要求2所述的一种提高大前进比旋翼配平收敛性的计算方法,其特征在于:若残差值高于收敛准则设定值,则迭代尚未收敛,进入残差监控环节;若残差值小于上一次迭代的残差值,则继续进行迭代。
4.根据权利要求3所述的一种提高大前进比旋翼配平收敛性的计算方法,其特征在于:若残差值大于与上一次迭代的残差值,则跳出迭代循环并根据此时的操纵值与配平目标重新生成Jacobian矩阵,重新开始新一轮的迭代循环。
5.根据权利要求4所述的一种提高大前进比旋翼配平收敛性的计算方法,其特征在于:若残差低于收敛准则设定值,则迭代收敛,当前前进比参数下的配平计算结果作为最新的操纵值。
7.根据权利要求1所述的一种提高大前进比旋翼配平收敛性的计算方法,其特征在于:所述第一步中,还包括:建立旋翼动力学模型和旋翼载荷计算模型。
8.根据权利要求7所述的一种提高大前进比旋翼配平收敛性的计算方法,其特征在于:所述第四步中,还包括:根据旋翼动力学模型和旋翼载荷计算模型进行旋翼气弹计算。
9.根据权利要求1所述的一种提高大前进比旋翼配平收敛性的计算方法,其特征在于:所述牛顿迭代法包括:通过以下公式计算操纵值;
λn=λn-1-[J]-1·f(λn-1)
其中:λn为当前迭代得到的操纵值,λn-1为上一次迭代的操纵值,J为Jacobian矩阵,f为配平平衡方程组。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011028637.3A CN112052529B (zh) | 2020-09-25 | 2020-09-25 | 一种提高大前进比旋翼配平收敛性的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011028637.3A CN112052529B (zh) | 2020-09-25 | 2020-09-25 | 一种提高大前进比旋翼配平收敛性的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112052529A true CN112052529A (zh) | 2020-12-08 |
CN112052529B CN112052529B (zh) | 2023-03-17 |
Family
ID=73604945
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011028637.3A Active CN112052529B (zh) | 2020-09-25 | 2020-09-25 | 一种提高大前进比旋翼配平收敛性的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112052529B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060047737A1 (en) * | 2004-08-30 | 2006-03-02 | Synopsys, Inc. | Simulation of processes, devices and circuits by a modified newton method |
AU2006225121A1 (en) * | 2005-03-15 | 2006-09-21 | Chevron U.S.A. Inc. | Stable method and apparatus for solving S-shaped non -linear functions utilizing modified Newton-Raphson algorithms |
CN102520620A (zh) * | 2011-11-18 | 2012-06-27 | 南京航空航天大学 | 单旋翼直升机/涡轴发动机通用综合模型的构建方法 |
CN104317980A (zh) * | 2014-08-26 | 2015-01-28 | 中国直升机设计研究所 | 一种共轴刚性旋翼气动弹性响应方法 |
CN104615086A (zh) * | 2014-12-25 | 2015-05-13 | 江苏高精机电装备有限公司 | 一种基于流函数的螺旋桨桨叶的数控加工方法 |
CN105260492A (zh) * | 2015-04-27 | 2016-01-20 | 中国直升机设计研究所 | 一种旋翼与机体耦合动力学模态综合建模方法 |
CN107330152A (zh) * | 2017-06-01 | 2017-11-07 | 南京航空航天大学 | 一种适用于旋翼飞行器的高效气动配平方法 |
CN108628803A (zh) * | 2017-03-20 | 2018-10-09 | 北京佳士乐动漫科技有限公司 | 一种基于牛顿迭代的线性方程求解方法 |
CN109614633A (zh) * | 2018-10-25 | 2019-04-12 | 南京航空航天大学 | 一种复合式旋翼飞行器非线性建模方法及配平方法 |
CN110435877A (zh) * | 2019-08-30 | 2019-11-12 | 西安倾云无人机技术有限公司 | 一种自适应气动变距螺旋桨 |
CN112199784A (zh) * | 2020-11-04 | 2021-01-08 | 南京航空航天大学 | 一种共轴刚性双旋翼气动配平方法及系统 |
-
2020
- 2020-09-25 CN CN202011028637.3A patent/CN112052529B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060047737A1 (en) * | 2004-08-30 | 2006-03-02 | Synopsys, Inc. | Simulation of processes, devices and circuits by a modified newton method |
AU2006225121A1 (en) * | 2005-03-15 | 2006-09-21 | Chevron U.S.A. Inc. | Stable method and apparatus for solving S-shaped non -linear functions utilizing modified Newton-Raphson algorithms |
CN102520620A (zh) * | 2011-11-18 | 2012-06-27 | 南京航空航天大学 | 单旋翼直升机/涡轴发动机通用综合模型的构建方法 |
CN104317980A (zh) * | 2014-08-26 | 2015-01-28 | 中国直升机设计研究所 | 一种共轴刚性旋翼气动弹性响应方法 |
CN104615086A (zh) * | 2014-12-25 | 2015-05-13 | 江苏高精机电装备有限公司 | 一种基于流函数的螺旋桨桨叶的数控加工方法 |
CN105260492A (zh) * | 2015-04-27 | 2016-01-20 | 中国直升机设计研究所 | 一种旋翼与机体耦合动力学模态综合建模方法 |
CN108628803A (zh) * | 2017-03-20 | 2018-10-09 | 北京佳士乐动漫科技有限公司 | 一种基于牛顿迭代的线性方程求解方法 |
CN107330152A (zh) * | 2017-06-01 | 2017-11-07 | 南京航空航天大学 | 一种适用于旋翼飞行器的高效气动配平方法 |
CN109614633A (zh) * | 2018-10-25 | 2019-04-12 | 南京航空航天大学 | 一种复合式旋翼飞行器非线性建模方法及配平方法 |
CN110435877A (zh) * | 2019-08-30 | 2019-11-12 | 西安倾云无人机技术有限公司 | 一种自适应气动变距螺旋桨 |
CN112199784A (zh) * | 2020-11-04 | 2021-01-08 | 南京航空航天大学 | 一种共轴刚性双旋翼气动配平方法及系统 |
Non-Patent Citations (2)
Title |
---|
余智豪等: "大前进比变转速旋翼气弹动力学建模与载荷特性分析", 《振动与冲击》 * |
吴杰: "直升机旋翼气动弹性振动载荷研究", 《中国优秀博硕士学位论文全文数据库(博士)工程科技Ⅱ辑(月刊)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112052529B (zh) | 2023-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102520620B (zh) | 单旋翼直升机/涡轴发动机通用综合模型的构建方法 | |
CN111523267B (zh) | 一种基于参数化有限元模型的风机主轴结构优化方法 | |
CN107665278B (zh) | 一种大型直驱风电机组的心轴强度优化设计方法 | |
CN102653315B (zh) | 一种提高无轴承旋翼气弹稳定性的装置及其设计方法 | |
CN107762730B (zh) | 一种带有尾缘襟翼的大型变桨风力机控制系统及控制方法 | |
CN106483405B (zh) | 基于隐马尔可夫模型的npc光伏逆变器的故障诊断方法 | |
CN105888971A (zh) | 一种大型风力机叶片主动降载控制系统及方法 | |
CN113283024A (zh) | 一种基于有限元建模的拓扑优化和参数优化的风机轮毂结构优化方法 | |
CN115525988B (zh) | 一种风电机组自主载荷仿真计算与修正系统 | |
CN112031996A (zh) | 一种优化风电机组偏航运动时机舱振动超限的方法及系统 | |
CN117108445B (zh) | 一种串列式双风轮风电机组数字孪生仿真方法 | |
CN112052529B (zh) | 一种提高大前进比旋翼配平收敛性的计算方法 | |
Fasel et al. | Aeroservoelastic optimization of morphing airborne wind energy wings | |
CN113221251B (zh) | 一种高超声速飞行器整机结构设计方法及系统 | |
CN112199784B (zh) | 一种共轴刚性双旋翼气动配平方法及系统 | |
CN114065514A (zh) | 一种串列式双风轮风电机组仿真的bem数据映射方法及系统 | |
CN116305589B (zh) | 一种直升机桨叶结构降阶分析方法、系统、设备及介质 | |
Anusonti-Inthra | Full vehicle simulations for a coaxial rotorcraft using high-fidelity CFD/CSD coupling | |
CN116628854A (zh) | 一种翼型气动特性预测方法、系统、电子设备及存储介质 | |
Cocco et al. | Tiltrotor whirl flutter analysis in support of NGCTR aeroelastic wind tunnel model design | |
CN105260498B (zh) | 一种大型民机机翼变弯度设计方法 | |
CN114186334A (zh) | 一种引入全弹性边界条件的旋翼气弹耦合计算方法 | |
CN109726407A (zh) | 一种风力发电机组变桨轴承连接螺栓极限强度计算方法 | |
CN112069588A (zh) | 一种风电机组塔架屈曲稳定性的预测方法及系统 | |
CN117077403A (zh) | 一种基于主动扭转旋翼的直升机飞行性能优化方法 |
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 |