CN114186334A - 一种引入全弹性边界条件的旋翼气弹耦合计算方法 - Google Patents
一种引入全弹性边界条件的旋翼气弹耦合计算方法 Download PDFInfo
- Publication number
- CN114186334A CN114186334A CN202111382566.1A CN202111382566A CN114186334A CN 114186334 A CN114186334 A CN 114186334A CN 202111382566 A CN202111382566 A CN 202111382566A CN 114186334 A CN114186334 A CN 114186334A
- Authority
- CN
- China
- Prior art keywords
- blade
- rotor
- section
- model
- deformation
- 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
- 230000008878 coupling Effects 0.000 title claims abstract description 42
- 238000010168 coupling process Methods 0.000 title claims abstract description 42
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 42
- 238000004364 calculation method Methods 0.000 title claims abstract description 31
- 230000005489 elastic deformation Effects 0.000 claims description 42
- 239000011159 matrix material Substances 0.000 claims description 18
- 238000000034 method Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 239000006185 dispersion Substances 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 4
- 230000006698 induction Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 1
- 238000002789 length control Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明属于直升机旋翼设计与理论建模技术,涉及一种用于旋翼气动弹性耦合仿真分析的计算方法。针对旋翼气动弹性耦合问题,在含自由尾迹的旋翼载荷计算模型中引入弹性桨叶旋翼动力学模型,替代刚体挥舞动力学模型,更真实地反映桨叶气动弹性耦合的物理现象,具有更高的旋翼振动载荷预估精度。
Description
技术领域
本发明属于直升机旋翼设计与理论建模技术,涉及一种用于旋翼气动弹性耦合仿真分析的计算方法。
背景技术
旋翼振动载荷预估是贯穿直升机设计阶段的重要环节之一,是一个涉及旋翼气动、结构、数值计算等多学科耦合的综合分析,也是最具挑战的难点之一。细长梁结构的桨叶在旋转工作时受气动载荷影响而产生结构变形,而桨叶结构变形又引起旋翼流场变化进而改变桨叶气动载荷分布,二者存在复杂的耦合关系。尤其在大前进比、变转速等工作条件下以及存在强烈干扰的多旋翼构型直升机当中,气动与结构的耦合关系影响更为显著。
在传统含自由尾迹的旋翼载荷计算模型中,将桨叶假设为刚体且只采用刚体挥舞模型,简化了桨叶挥舞运动,忽略桨叶摆振与扭转运动以及弹性变形,即未反映实际桨叶真实的挥/摆/扭运动以及弹性变形引起的旋翼尾迹变化以及气动载荷变化,在大前进比飞行状态,变转速飞行状态以及存在强烈干扰的多旋翼构型的直升机前飞状态中,桨叶结构变形不可忽略,对旋翼振动载荷结果影响不容忽略,而未考虑弹性变形的旋翼载荷计算模型精度明显不足,无法真实反映气动弹性耦合问题。
发明内容
发明目的:本发明提出一种引入全弹性边界条件的旋翼气弹耦合计算方法,针对旋翼气动弹性耦合问题,在含自由尾迹的旋翼载荷计算模型中引入弹性桨叶旋翼动力学模型,替代刚体挥舞动力学模型,更真实地反映桨叶气动弹性耦合的物理现象,具有更高的旋翼振动载荷预估精度。
发明技术方案:提供一种引入全弹性边界条件的旋翼气弹耦合计算方法,所述方法包括:
建立旋翼动力学模型。首先,建立含桨叶挥舞、摆振、扭转的非线性弹性变形量的运动方程,所述运动方程中包含有非线性弹性变形量、所述非线性弹性变形量之间的非线性耦合关系和剖面运动速度;然后,将所述运动方程和非线性耦合关系引入旋翼动力学模型中;
建立旋翼气动模型。首先,采用非定常L-B模型建立桨叶气动力模型;然后,采用涡流理论,通过线涡离散方式对旋翼桨尖涡进行空间离散,并采用粘性涡核经验模型反映桨尖涡的畸变与耗散,建立自由尾迹模型;将所述桨叶气动力模型和自由尾迹模型集成为旋翼气动模型;
设定初始的旋翼工作状态,通过自由尾迹模型计算得到旋翼的流场分布,所述流场分布包括桨叶各剖面3/4弦点位置处的合速度;将所述桨叶各剖面3/4弦点位置处的合速度作为输入传递给桨叶气动力模型计算得到桨叶各剖面气动力,所述各剖面气动力包括升力、阻力和力矩;将桨叶各剖面气动力作为输入传递给旋翼动力学模型,求解有限元节点处的1/4弦点位置弹性变形及剖面运动;通过有限元形函数差值和线性差值方法得到整片桨叶任意剖面3/4弦点位置的弹性变形及剖面运动;
将所述整片桨叶任意剖面3/4弦点位置的弹性变形及剖面运动作为输入传递给自由尾迹模型,并计算得到旋翼气弹耦合之后的桨叶各剖面3/4弦点位置处的合速度;将旋翼气弹耦合之后的桨叶各剖面3/4弦点位置处的合速度作为输入传递给桨叶气动力模型,并计算得到旋翼气弹耦合之后的桨叶各剖面气动力;将旋翼气弹耦合之后的桨叶各剖面气动力作为输入传递给旋翼动力学模型,求解有限元节点处的1/4弦点位置弹性变形及剖面运动;然后桨叶旋翼气弹耦合之后的各剖面的弹性变形及剖面运动再次作为输入传递给自由尾迹模型,进行迭代循环,当整片桨叶各剖面的弹性变形趋于收敛后,得到最终整片桨叶各剖面的弹性变形及剖面运动。
进一步地,确定桨叶变形之后的空间坐标关系,根据所述空间坐标关系建立运动方程。
进一步地,确定桨叶未变形坐标系、桨叶变形坐标系和桨毂中心固定坐标系;建立桨叶未变形坐标系和桨毂中心固定坐标系之间的转换矩阵TBH、桨叶变形坐标系和未变形坐标系转换矩阵TEB;通过转换矩阵TEB得到桨叶在未变形坐标下的失径根据转换矩阵TBH得到桨毂中心固定坐标系下桨叶各剖面的非线性弹性变形量及其非线性弹性变形量之间的非线性耦合关系。
进一步地,所述非线性弹性变形量包括桨叶拉伸变形u、摆振变形v、挥舞变形w和扭转变形p。
进一步地,初始的旋翼工作状态包括旋翼的初始状态参数、旋翼气动参数和结构参数。
进一步地,桨叶各剖面3/4弦点位置处的合速度的计算公式为:
附图标记说明
图1是本发明涉及的旋翼动力学模型中的弹性桨叶变形示意图;
图2是本发明涉及的某时刻下桨叶挥舞和摆振位移结果对比示意图,其中点划线为传统的刚体挥舞动力学模型计算结果,实线为本发明涉及的旋翼动力学方程建模中的弹性桨叶变形计算结果;(a)图为挥舞位移(b)图为摆振位移;
图3是本发明涉及的旋翼旋转周期内桨叶运动对比示意图,其中点划线为刚体桨叶模型计算结果,实线为弹性桨叶模型计算结果;
图4是本发明涉及的引入全弹性边界条件后的尾迹几何结果示意图,其中粗实线表示桨叶形状,细实线为桨尖涡曲线(尾迹曲线)。
具体实施方式
本实施例,提供一种引入全弹性边界条件的旋翼气弹耦合计算方法包括:
(1)建立旋翼动力学模型,建立含桨叶挥舞、摆振、扭转的非线性弹性变形量的运动方程,包含有非线性弹性变形量、所述非线性弹性变形量之间的非线性耦合关系和剖面运动速度,进而准确描述弹性桨叶的非线性运动,弹性桨叶变形如图1,在桨叶剖面e处通过坐标转换形式得到桨叶变形坐标系E与桨叶未变形坐标系B下的几何关系(包括轴向变形u,摆振变形v,和挥舞变形w)。并将系统能量分为应变能U,动能T,外载荷W,最后采用Hamilton原理建立旋翼动力学方程(其中δ表示变分)。
∫(δU-δT-δW)dt=0
(2)建立含自由尾迹的旋翼气动模型,采用非定常L-B模型建立桨叶气动力模型;采用涡流理论,通过线涡离散方式对旋翼桨尖涡进行空间离散,并采用粘性涡核经验模型反映桨尖涡的畸变与耗散,建立自由尾迹模型并最终集成旋翼气动模型。
(3)在自由尾迹模型的输入中引入全弹性桨叶变形的边界条件,将旋翼动力学模型中求解的整片桨叶有限元节点剖面1/4弦点位置的弹性变形及剖面运动通过有限元形函数差值、线性差值以及坐标系矩阵转换方式得到整片桨叶任意剖面3/4弦点位置的弹性变形及剖面运动,包含轴向变形u、摆振变形v、和挥舞变形w、扭转变形p以及剖面运动参数最后作为输入传递给自由尾迹模型,进行气弹耦合计算。
具体地,包括如下步骤:
第一步:引入弹性变形的弹性桨叶控制点空间坐标计算。
设定初始的旋翼工作状态,自由尾迹模型中采用的W-L升力面理论是将桨叶剖面信息集中反映至3/4弦点位置的控制点处,而旋翼动力学模型求解出的桨叶弹性变形以及剖面运动参数结果是1/4弦长有限元节点处,首先需通过差值方法得到整片桨叶任意1/4弦点处的桨叶变形及运动信息。
然后利用空间坐标系关系(桨叶未变形坐标系和桨毂中心固定坐标系之间的转换矩阵TBH,桨叶变形坐标系和未变形坐标系转换矩阵TEB),通过坐标矩阵转化得到桨叶在未变形坐标下的失径(表达式如下),其中x为剖面位置,u,v,w为当前剖面的拉伸,摆振和挥舞变形,为变形坐标系下的剖面上任意一点坐标,即可以得到3/4弦长控制点处的失径。最后,通过TBH坐标系转换与线性差值方法,最终得到桨毂中心固定坐标系下的桨叶任意3/4弦长控制点处的失径,即包含了桨叶拉伸变形u、摆振变形v、挥舞变形w和扭转变形p。
转换矩阵TBH与TEB表达式如下:
其中δ为旋翼轴倾角,ψ为旋转角,βP为预锥角,[I J K]为基矢量,下标H表示桨毂中心固定坐标系,B表示桨叶未变形坐标系。
第二步:包含弹性桨叶运动的控制点处合速度求解。
在自由尾迹模型中利用桨叶控制点处(3/4弦点)的气流不穿透条件求解该点处的合速度,表达式如下:
第三步:气弹耦合计算
通过自由尾迹模型计算得到旋翼桨叶各剖面3/4弦点位置处的合速度;并将桨叶各剖面3/4弦点位置处的合速度作为输入传递给桨叶气动力模型,并计算得到桨叶各剖面气动力(包括升力、阻力和力矩);将桨叶各剖面气动力作为输入传递给旋翼动力学模型,求解有限元节点处的1/4弦点位置弹性变形及剖面运动;通过有限元形函数差值和线性差值方法得到整片桨叶任意剖面3/4弦点位置的弹性变形及剖面运动;然后将整片桨叶任意剖面3/4弦点位置的弹性变形及剖面运动再次作为输入传递给自由尾迹模型,进行迭代循环,当整片桨叶各剖面的弹性变形趋于收敛后,得到最终整片桨叶各剖面的弹性变形及剖面运动。
Claims (8)
1.一种引入全弹性边界条件的旋翼气弹耦合计算方法,其特征在于,所述方法包括:
建立旋翼动力学模型:首先,建立含桨叶挥舞、摆振、扭转的非线性弹性变形的运动方程,所述运动方程中包含有非线性弹性变形量、所述非线性弹性变形量之间的非线性耦合关系和剖面运动速度;然后,将所述运动方程和非线性耦合关系引入旋翼动力学模型中;
建立旋翼气动模型:首先,采用非定常L-B模型建立桨叶气动力模型;然后,采用涡流理论,通过线涡离散方式对旋翼桨尖涡进行空间离散,并采用粘性涡核经验模型反映桨尖涡的畸变与耗散,建立自由尾迹模型;将所述桨叶气动力模型和自由尾迹模型集成为旋翼气动模型;
设定初始的旋翼工作状态,通过自由尾迹模型计算得到旋翼的流场分布,所述流场分布包括桨叶各剖面3/4弦点位置处的合速度;将所述桨叶各剖面3/4弦点位置处的合速度作为输入传递给桨叶气动力模型计算得到桨叶各剖面气动力,所述各剖面气动力包括升力、阻力和力矩;将桨叶各剖面气动力作为输入传递给旋翼动力学模型,求解有限元节点处的1/4弦点位置弹性变形及剖面运动;通过有限元形函数差值和线性差值方法得到整片桨叶任意剖面3/4弦点位置的弹性变形及剖面运动;
将所述整片桨叶任意剖面3/4弦点位置的弹性变形及剖面运动作为输入传递给自由尾迹模型,并计算得到旋翼气弹耦合之后的桨叶各剖面3/4弦点位置处的合速度;将旋翼气弹耦合之后的桨叶各剖面3/4弦点位置处的合速度作为输入传递给桨叶气动力模型计算得到旋翼气弹耦合之后的桨叶各剖面气动力;将旋翼气弹耦合之后的桨叶各剖面气动力作为输入传递给旋翼动力学模型,求解有限元节点处的1/4弦点位置弹性变形及剖面运动;然后桨叶旋翼气弹耦合之后的各剖面的弹性变形及剖面运动再次输入自由尾迹模型,依次迭代,当整片桨叶各剖面的弹性变形趋于收敛后,得到最终整片桨叶各剖面的弹性变形及剖面运动。
2.根据权利要求1所述的旋翼气弹耦合计算方法,其特征在于,确定桨叶变形之后的空间坐标关系,根据所述空间坐标关系建立运动方程。
4.根据权利要求1所述的旋翼气弹耦合计算方法,其特征在于,所述非线性弹性变形量包括桨叶拉伸变形u、摆振变形v、挥舞变形w和扭转变形p。
7.根据权利要求1所述的旋翼气弹耦合计算方法,其特征在于,初始的旋翼工作状态包括旋翼的初始状态参数、旋翼气动参数和结构参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111382566.1A CN114186334B (zh) | 2021-11-19 | 2021-11-19 | 一种引入全弹性边界条件的旋翼气弹耦合计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111382566.1A CN114186334B (zh) | 2021-11-19 | 2021-11-19 | 一种引入全弹性边界条件的旋翼气弹耦合计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114186334A true CN114186334A (zh) | 2022-03-15 |
CN114186334B CN114186334B (zh) | 2024-09-24 |
Family
ID=80541097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111382566.1A Active CN114186334B (zh) | 2021-11-19 | 2021-11-19 | 一种引入全弹性边界条件的旋翼气弹耦合计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114186334B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118094780A (zh) * | 2024-04-23 | 2024-05-28 | 华东交通大学 | 基于相邻涵道中螺旋桨气动干扰的动力学机理建模方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009022033A1 (es) * | 2007-08-16 | 2009-02-19 | Indra Sistemas, S.A. | Procedimiento de simulación en tiempo real de un rotor de helicóptero |
KR20140140279A (ko) * | 2013-05-29 | 2014-12-09 | 한국항공우주산업 주식회사 | 회전익 항공기 미부 공력진동 예측시스템 및 그 제어방법 |
CN104317980A (zh) * | 2014-08-26 | 2015-01-28 | 中国直升机设计研究所 | 一种共轴刚性旋翼气动弹性响应方法 |
CN112528408A (zh) * | 2020-12-11 | 2021-03-19 | 中国直升机设计研究所 | 一种直升机旋翼与机身耦合稳定性建模方法 |
-
2021
- 2021-11-19 CN CN202111382566.1A patent/CN114186334B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009022033A1 (es) * | 2007-08-16 | 2009-02-19 | Indra Sistemas, S.A. | Procedimiento de simulación en tiempo real de un rotor de helicóptero |
KR20140140279A (ko) * | 2013-05-29 | 2014-12-09 | 한국항공우주산업 주식회사 | 회전익 항공기 미부 공력진동 예측시스템 및 그 제어방법 |
CN104317980A (zh) * | 2014-08-26 | 2015-01-28 | 中国直升机设计研究所 | 一种共轴刚性旋翼气动弹性响应方法 |
CN112528408A (zh) * | 2020-12-11 | 2021-03-19 | 中国直升机设计研究所 | 一种直升机旋翼与机身耦合稳定性建模方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118094780A (zh) * | 2024-04-23 | 2024-05-28 | 华东交通大学 | 基于相邻涵道中螺旋桨气动干扰的动力学机理建模方法 |
CN118094780B (zh) * | 2024-04-23 | 2024-07-02 | 华东交通大学 | 基于相邻涵道中螺旋桨气动干扰的动力学机理建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114186334B (zh) | 2024-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109614633B (zh) | 一种复合式旋翼飞行器非线性建模及线性化配平方法 | |
Hui et al. | Aerodynamic performance investigation on a morphing unmanned aerial vehicle with bio-inspired discrete wing structures | |
Gamboa et al. | Optimization of a morphing wing based on coupled aerodynamic and structural constraints | |
Bae et al. | Aerodynamic and static aeroelastic characteristics of a variable-span morphing wing | |
CN113806871B (zh) | 一种考虑结构非线性的柔性飞行动力学建模方法 | |
CN104317980A (zh) | 一种共轴刚性旋翼气动弹性响应方法 | |
CN111581784B (zh) | 一种基于数据驱动自适应准稳态模型的扑翼运动参数优化方法 | |
Öhrle et al. | Compound helicopter X 3 in high-speed flight: correlation of simulation and flight test | |
CN113051662B (zh) | 一种基于cfd和datcom的折叠翼尖变体飞行器的气动建模及其性能评估方法 | |
CN117108445B (zh) | 一种串列式双风轮风电机组数字孪生仿真方法 | |
Bosch et al. | Nonlinear aeroelasticity, flight dynamics and control of a flexible membrane traction kite | |
Fasel et al. | Aeroservoelastic optimization of morphing airborne wind energy wings | |
Fengqian et al. | Continuous morphing trailing-edge wing concept based on multi-stable nanomaterial | |
CN114186334B (zh) | 一种引入全弹性边界条件的旋翼气弹耦合计算方法 | |
Cong et al. | Design and analysis of passive variable-pitch propeller for VTOL UAVs | |
Yang et al. | Aeroelastic trim and flight loads analysis of flexible aircraft with large deformations | |
CN113221237B (zh) | 一种基于降阶建模的大迎角颤振分析方法 | |
Candade et al. | Aero-structural design of composite wings for airborne wind energy applications | |
CN117034451A (zh) | 一种旋翼/螺旋桨气动干扰计算方法 | |
Cooper et al. | Optimization of a scaled SensorCraft model with passive gust alleviation | |
Anusonti-Inthra et al. | Design of a Comfortable Rotor Airfoil Using Distributed Piezoelectric Actuators | |
CN115982837A (zh) | 一种多桨倾转机翼旋翼飞行器多系统耦合动力学建模方法 | |
CN115729264A (zh) | 一种基于柔性自适应翼梢小翼的变稳隐身飞机控制方法 | |
Tchatchueng Kammegne et al. | Control validation of a morphing wing in an open loop architecture | |
Wang et al. | Static aeroelastic analysis of flexible aircraft with large deformations |
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 |