CN113536456B - 一种基于动力学响应辨识的大柔性结构降阶模型建模方法 - Google Patents

一种基于动力学响应辨识的大柔性结构降阶模型建模方法 Download PDF

Info

Publication number
CN113536456B
CN113536456B CN202110746615.9A CN202110746615A CN113536456B CN 113536456 B CN113536456 B CN 113536456B CN 202110746615 A CN202110746615 A CN 202110746615A CN 113536456 B CN113536456 B CN 113536456B
Authority
CN
China
Prior art keywords
structural
discrete
order
mode
dynamic response
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.)
Active
Application number
CN202110746615.9A
Other languages
English (en)
Other versions
CN113536456A (zh
Inventor
谢长川
张铎耀
安朝
杨超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202110746615.9A priority Critical patent/CN113536456B/zh
Publication of CN113536456A publication Critical patent/CN113536456A/zh
Application granted granted Critical
Publication of CN113536456B publication Critical patent/CN113536456B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • G06F17/141—Discrete Fourier transforms
    • 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)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Operations Research (AREA)
  • Discrete Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Feedback Control In General (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种基于动力学响应辨识的大柔性结构降阶模型建模方法，将非线性结构动力学方程以等时间间隔进行离散化，形成广义坐标离散序列，然后对其进行离散傅里叶变换，形成频域变换序列，交换移动方程部分项顺序构成线性方程组，利用模态组合乘比例系数的方式构建测试载荷加载于结构上计算动力学响应，提前动响应数据后求方程组的最小二乘解即可给出非线性刚度系数，建立非线性结构降阶模型。本发明构建的基于动力学响应辨识的大柔性结构降阶模型，与传统的基于静力学响应数据辨识的大柔性结构降阶模型相比，利用结构动力学响应数据样本进行非线性刚度系数的辨识，可以提供更准确的结构动力学特性信息。

Description

一种基于动力学响应辨识的大柔性结构降阶模型建模方法
技术领域
本发明涉及技术领域，尤其涉及一种基于动力学响应辨识的大柔性结构降阶模型建模方法。
背景技术
现代飞行器设计对于飞行性能的要求不断提高，以高空长航时无人机、飞翼布局无人机为代表的大柔性飞行器，自20世纪90年代以来逐渐成为各航空大国的研究重点。该类飞行器能够以更灵活的方式不间断地执行诸如军事侦察、民用科学探索及通信中继等任务，具有远超传统飞行器的飞行高度及续航能力，具有广阔的应用前景。由于其飞行性能的要求，大柔性飞行器通常利用较大的机翼展弦比来提高升阻比，同时在结构设计汇总增加复合材料比例，这些因素导致这类飞行器相对结构重量低、柔性大，气动力加载引起的结构变形较大，带来结构几何非线性问题。
结构建模是气动弹性分析中的核心问题，准确的大变形结构建模方法是几何非线性气动弹性分析的基础。大变形梁理论与非线性有限元模型是传统几何非线性气动弹性系统中经常使用的两类大柔性结构建模方法。气动弹性理论分析中经常使用大变形梁理论进行结构建模，其计算效率高，系统阶数低，但模型适用性差，只能用于简单模型的分析。非线性有限元模型是工程分析中经常使用的大柔性结构建模方法，其模型适用性强，可以应用于复杂结构模型的分析计算，但模型自由度数大，非线性阶数高，收敛性差计算效率，且难以对结构动力学特性进行机理性研究。非线性结构降阶模型是一种能够较好的平衡求解效率、计算精度及复杂模型适用性的大柔性结构建模方法。非线性结构降阶模型中，将非线性刚度表达为广义坐标多项式形式，利用回归分析方法求解非线性刚度系数给出结构动力学方程。
在“Static aeroelastic analysis including geometric nonlinearitesbased on reduced order model”Xie Changchuan等，《Chinese Journal ofAeronautics》,第638-650页，20170216中，将非线性结构降阶模型引入几何非线性静气动弹性分析，构建静气动弹性分析框架。在“Flutter and gust response analysis of awing model including geometric nonlinearities based on a modified structuralROM”An Chao等，《Chinese Journal of Aeronautics》,第48-63页，20190731以及公开号为CN112580241A的专利“一种基于结构降阶模型的非线性气动弹性动响应分析方法”中，将非线性结构降阶模型引入颤振与阵风分析。在公开号为CN112580240A的专利“一种适用于复杂大柔性飞机建模的非线性子结构方法”中，将非线性结构降阶模型引入子结构方法中，形成适用于复杂全机的大柔性结构建模方法。上述研究中，在计算非线性结构降阶模型中的非线性刚度系数前，在非线性有限元软件中计算一定数量的静力学载荷变形样本，在这些样本基础上，通过回归分析方法求解非线性刚度系数。虽然非线性静力学分析在得到样本时计算效率高，但是基于静力学数据的非线性刚度系数辨识使得系统动力学特性有一定程度的损失，不能提供完整准确的结构动力学特性信息。
发明内容
有鉴于此，本发明提供了一种基于动力学响应辨识的大柔性结构降阶模型建模方法，用以克服现有大柔性结构降阶模型建模中只用到静力学响应数据做辨识样本导致结构动力学特征无法被准确捕捉的不足，建立适用于大柔性结构变形分析的结构降阶模型。
本发明提供的一种基于动力学响应辨识的大柔性结构降阶模型建模方法，包括如下步骤：
S1：建立大柔性结构的有限元模型，划分有限元模型的结构有限元网格，定义结构有限元网格的单元属性和材料特性，设置计算边界条件；
S2：基于大柔性结构的有限元模型，计算大柔性结构在垂直弯曲方向、扭转方向及水平弯曲方向的线性模态；
S3：考虑包含几何非线性大变形作用的结构动力学方程，将结构动力学方程中的广义坐标和模态力向量以等时间间隔进行离散，得到离散形式的结构动力学方程；
S4：对离散形式的结构动力学方程进行离散傅里叶变换；
S5：选择有限元模型在垂直弯曲方向和扭转方向的线性模态，得到垂直弯曲方向和扭转方向的线性模态组合后的测试载荷，作为外载加载于大柔性结构上进行动力学响应计算，得到结构时域动响应位移数据；
S6：对测试载荷和结构时域动响应位移数据进行傅里叶变换，得到频域离散序列后代入经离散傅里叶变换后的结构动力学方程，利用最小二乘方法计算得到非线性刚度系数后代入步骤S3中未经离散的结构动力学方程，得到完整的非线性结构动力学方程。
在一种可能的实现方式中，在本发明提供的上述基于动力学响应辨识的大柔性结构降阶模型建模方法中，步骤S3，具体包括：
考虑包含几何非线性大变形作用的结构动力学方程如下：
Figure BDA0003144536260000031
方程(1)对应第m阶结构模态，并应用Einstein求和约定；其中，m＝1，…，N，N表示参与降阶的结构模态数目；Mm表示第m阶结构模态的广义质量项；Km表示第m阶结构模态的广义刚度项；
Figure BDA0003144536260000032
表示第m阶结构模态下的第n阶结构模态广义坐标与第l阶结构模态广义坐标的乘积多项式的非线性刚度项，
Figure BDA0003144536260000033
表示第m阶结构模态下的第n阶结构模态广义坐标、第l阶结构模态广义坐标和第p阶结构模态广义坐标的乘积多项式的非线性刚度系数；qm表示第m阶结构模态的广义坐标，qn表示第n阶结构模态的广义坐标，ql表示第l阶结构模态的广义坐标，qp表示第p阶结构模态的广义坐标；Fm表示第m阶结构模态的模态力向量；
将方程(1)中的广义坐标qm、qn、ql、qp及模态力向量Fm以等时间间隔进行离散，得到离散形式的结构动力学方程：
Figure BDA0003144536260000041
其中，i＝0，1，…，NT-1，表示离散时间点序列；NT表示离散时间点的总数量；qm(i)表示第m阶结构模态的广义坐标离散序列，qn(i)表示第n阶结构模态的广义坐标离散序列，ql(i)表示第l阶结构模态的广义坐标离散序列，qp(i)表示第p阶结构模态的广义坐标离散序列，Fm(i)表示第m阶结构模态的模态力向量离散序列。
在一种可能的实现方式中，在本发明提供的上述基于动力学响应辨识的大柔性结构降阶模型建模方法中，步骤S4，具体包括：
对离散形式的结构动力学方程中的广义坐标离散序列qm(i)、qn(i)、ql(i)和qp(i)进行离散傅里叶变换，变换后分别表示为Qm(k)、Qn(k)、Ql(k)和Qp(k)，变换关系为：
Figure BDA0003144536260000042
Figure BDA0003144536260000043
Figure BDA0003144536260000044
Figure BDA0003144536260000045
对离散形式的结构动力学方程中的广义坐标二阶项qnql和广义坐标三阶项qnqlqp进行离散傅里叶变换，变换后分别表示为
Figure BDA0003144536260000046
和
Figure BDA0003144536260000047
变换关系为：
Figure BDA0003144536260000051
Figure BDA0003144536260000052
对离散形式的结构动力学方程中的模态力向量Fm(i)进行离散傅里叶变换，变换后表示为fm(k)，变换关系为：
Figure BDA0003144536260000053
将方程(3)～(9)代入方程(2)中，得到经过傅里叶变换后的离散形式的结构动力学方程为：
Figure BDA0003144536260000054
在方程(3)～(10)中取NK个不同的频域k值{k1 k2 … kNK}，得到NK个方程组成的方程组：
Figure BDA0003144536260000055
在一种可能的实现方式中，在本发明提供的上述基于动力学响应辨识的大柔性结构降阶模型建模方法中，步骤S5中，得到的垂直弯曲方向和扭转方向的线性模态组合后的测试载荷为：
Ftest＝∑piφib+∑pjφjt  (12)
其中，pi和pj表示人为给定的比例系数，φib表示选定的线性弯曲模态振型，φjt表示选定的线性扭转模态振型。
在一种可能的实现方式中，在本发明提供的上述基于动力学响应辨识的大柔性结构降阶模型建模方法中，步骤S6中，对测试载荷和结构时域动响应位移数据进行傅里叶变换，得到频域离散序列后代入经离散傅里叶变换后的结构动力学方程，利用最小二乘方法计算得到非线性刚度系数，具体包括：
对测试载荷和结构时域动响应位移数据进行傅里叶变换，得到五个频域离散序列
Figure BDA0003144536260000061
Qm(k)、
Figure BDA0003144536260000062
和fm(k)，代入经离散傅里叶变换后的结构动力学方程组(11)后，利用最小二乘方法计算得到非线性刚度系数
Figure BDA0003144536260000063
和
Figure BDA0003144536260000064
本发明提供的上述基于动力学响应辨识的大柔性结构降阶模型建模方法，针对大柔性结构建模问题，基于动力学响应数据进行结构非线性刚度系数的辨识，提供一套准确可靠适用于大柔性结构变形分析的几何非线性结构降阶模型。将非线性结构动力学方程以等时间间隔进行离散化，形成广义坐标离散序列，然后对其进行离散傅里叶变换，形成频域变换序列，交换移动方程部分项顺序构成线性方程组，利用模态组合乘比例系数的方式构建测试载荷加载于结构上计算动力学响应，提前动响应数据后求方程组的最小二乘解即可给出非线性刚度系数，建立非线性结构降阶模型。本发明构建的基于动力学响应辨识的大柔性结构降阶模型，与传统的基于静力学响应数据辨识的大柔性结构降阶模型相比，利用结构动力学响应数据样本进行非线性刚度系数的辨识，可以提供更准确的结构动力学特性信息。
附图说明
图1为本发明提供的一种基于动力学响应辨识的大柔性结构降阶模型建模方法的流程图；
图2为本发明实施例1中机翼模型示意图；
图3为本发明实施例1中机翼模型有限元模型示意图；
图4为本发明实施例1中机翼线性模态振型图。
具体实施方式
下面将结合本发明实施方式中的附图，对本发明实施方式中的技术方案进行清楚、完整的描述，显然，所描述的实施方式仅仅是作为例示，并非用于限制本发明。
本发明提供的一种基于动力学响应辨识的大柔性结构降阶模型建模方法，如图1所示，包括如下步骤：
第一步：计算初始化。
建立大柔性结构的有限元模型，划分有限元模型的结构有限元网格，定义结构有限元网格的单元属性和材料特性，设置计算边界条件。
第二步：计算结构模态。
基于大柔性结构的有限元模型，计算大柔性结构在垂直弯曲方向、扭转方向及水平弯曲方向的线性模态。
第三步：结构降阶方程。
(1)结构动力学方程形式。
考虑包含几何非线性大变形作用的结构动力学方程如下：
Figure BDA0003144536260000071
方程(1)对应第m阶结构模态，并应用Einstein求和约定；其中，m＝1，…，N，N表示参与降阶的结构模态数目；Mm表示第m阶结构模态的广义质量项；Km表示第m阶结构模态的广义刚度项；
Figure BDA0003144536260000072
为非线性刚度项，上标表示对应的广义坐标的多项式阶数，下标对应的广义坐标的模态阶次，即
Figure BDA0003144536260000073
表示第m阶结构模态下的第n阶结构模态广义坐标与第l阶结构模态广义坐标的乘积多项式的非线性刚度项，
Figure BDA0003144536260000074
表示第m阶结构模态下的第n阶结构模态广义坐标、第l阶结构模态广义坐标和第p阶结构模态广义坐标的乘积多项式的非线性刚度系数；qm表示第m阶结构模态的广义坐标，qn表示第n阶结构模态的广义坐标，ql表示第l阶结构模态的广义坐标，qp表示第p阶结构模态的广义坐标；Fm表示第m阶结构模态的模态力向量；
(2)等时间间隔离散。
将方程(1)中的广义坐标qm、qn、ql、qp及模态力向量Fm以等时间间隔进行离散，得到离散形式的结构动力学方程：
Figure BDA0003144536260000081
其中，i＝0，1，…，NT-1，表示离散时间点序列；NT表示离散时间点的总数量；qm(i)表示第m阶结构模态的广义坐标离散序列，qn(i)表示第n阶结构模态的广义坐标离散序列，ql(i)表示第l阶结构模态的广义坐标离散序列，qp(i)表示第p阶结构模态的广义坐标离散序列，Fm(i)表示第m阶结构模态的模态力向量离散序列。
第四步：结构动力学方程转换为频域序列。
对离散形式的结构动力学方程中的广义坐标离散序列qm(i)、qn(i)、ql(i)和qp(i)进行离散傅里叶变换，变换后分别表示为Qm(k)、Qn(k)、Ql(k)和Qp(k)，变换关系为：
Figure BDA0003144536260000082
Figure BDA0003144536260000083
Figure BDA0003144536260000084
Figure BDA0003144536260000085
同理，对离散形式的结构动力学方程中的广义坐标二阶项qnql和广义坐标三阶项qnqlqp进行离散傅里叶变换，变换后分别表示为
Figure BDA0003144536260000086
和
Figure BDA0003144536260000087
变换关系为：
Figure BDA0003144536260000088
Figure BDA0003144536260000091
同理，对离散形式的结构动力学方程中的模态力向量Fm(i)进行离散傅里叶变换，变换后表示为fm(k)，变换关系为：
Figure BDA0003144536260000092
将方程(3)～(9)代入方程(2)中，得到经过傅里叶变换后的离散形式的结构动力学方程为：
Figure BDA0003144536260000093
在方程(3)～(10)中取NK个不同的频域k值{k1 k2 … kNK}，得到NK个方程组成的方程组：
Figure BDA0003144536260000094
第五步：确定测试载荷及变形样本。
测试载荷的给定是降阶模型建立的关键部分，测试载荷合适与否直接影响回归分析的准确性及降阶模型的合理性。
选择有限元模型在垂直弯曲方向和扭转方向的线性模态，得到垂直弯曲方向和扭转方向的线性模态组合后的测试载荷为：
Ftest＝∑piφib+∑pjφjt  (12)
其中，pi和pj表示人为给定的比例系数，φib表示选定的线性弯曲模态振型，φjt表示选定的线性扭转模态振型。将Ftest作为外载加载于大柔性结构上进行动力学响应计算，得到结构时域动响应位移数据。
第六步：求解非线性刚度系数，建立非线性结构降阶模型。
对测试载荷和结构时域动响应位移数据进行傅里叶变换，得到五个频域离散序列
Figure BDA0003144536260000101
Qm(k)、
Figure BDA0003144536260000102
和fm(k)，代入经离散傅里叶变换后的结构动力学方程组(11)后，利用最小二乘方法计算得到非线性刚度系数
Figure BDA0003144536260000103
和
Figure BDA0003144536260000104
将非线性刚度系数
Figure BDA0003144536260000105
和
Figure BDA0003144536260000106
代入步骤S3中未经离散的结构动力学方程，即可得到完整的非线性结构动力学方程。
下面通过一个具体的实施例对本发明提供的上述基于动力学响应辨识的大柔性结构降阶模型建模方法的具体实施进行详细说明。
实施例1：
采用大柔性机翼，翼型选用NACA0015对称翼型，翼根完全固支。机翼参数如下：翼根弦长100mm，翼尖弦长100mm，展长1000mm，后掠角0Deg，无扭转角。机翼主梁采用矩形截面，截面长35mm，截面宽1.5mm，位于机翼50％弦长处。机翼外形由8个轻木盒段来维持，每个盒段由轻木骨架与热缩膜蒙皮组成。为了减小由于弯曲变形使得盒段相互挤压而产生对机翼刚度的影响，盒段与主梁单点连接并且相邻盒段之间留有一定间隙。为了调整机翼的动力学特性和结构变形，在翼尖设置了配重杆，配重杆总长200mm。模型总重0.613kg。该机翼模型示意图如图2所示，包括配重杆1、机翼翼梁2和维形框3，配重杆1分为配重杆前锥1a、配重杆中段1b和配重杆后锥1c，机翼翼梁2的一端2a与配重杆1固定连接，另一端2b为固支端。
本发明实施例1提供的基于动力学响应辨识的大柔性结构降阶模型建模方法的步骤如下：
第一步：计算初始化。
建立机翼的有限元模型，利用有限元软件MSC.Nasrtan完成有限元建模。有限元模型采用CBEAM单元模拟主梁及翼肋，采用CONM2单元模拟集中质量，FORCE1添加随动力结构载荷，根部六自由度固支。有限元模型信息如表1所示。有限元模型共有节点213个，梁单元212个，集中质量单元65个，在51个节点上加载随动力，有限元模型示意图如图3所示，图3中的坐标系沿来流方向为x轴，沿展向为y轴，满足右手定则。
表1有限元模型信息
模型节点数 213
模型梁单元数 212
模型集中质量数 65
随动力数 51
第二步：计算结构模态。
有限元模型建立完成后，利用Nastran软件SOL1 03线性模态求解器计算结构模态，前5阶模态及第8阶模态二阶扭转模态的信息如表2所示。振型信息如图4所示，图4中的坐标系沿来流方向为x轴，沿展向为y轴，满足右手定则。第1阶模态为垂直一弯模态(如图4中的(a)所示)，固有频率1.179Hz，模态刚度5.488E+01，固有频率较低，模型柔度较大；第2阶模态为垂直二弯模态(如图4中的(b)所示)，固有频率7.724Hz，模态刚度2.355E+03；第3阶模态为垂直三弯模态(如图4中的(c)所示)，固有频率22.19Hz，模态刚度1.944E+04；第4阶模态为一阶扭转模态(如图4中的(d)所示)，固有频率22.95Hz，模态刚度2.080E+04；第5阶模态为水平一弯模态(如图4中的(e)所示)，固有频率27.47Hz，模态刚度2.978E+04；第8阶模态为二阶扭转模态(如图4中的(f)所示)，固有频率95.72Hz，模态刚度3.617E+05。
表2线性模态信息
Figure BDA0003144536260000111
Figure BDA0003144536260000121
第三步：结构降阶方程。
结构动力学方程形如方程(1)所示，本实施例1中采用前5阶线性模态，方程(1)中已知项有如下数值：
Figure BDA0003144536260000122
Figure BDA0003144536260000123
以等时间间隔进行离散后的动力学方程为：
Figure BDA0003144536260000124
第四步：结构动力学方程转换为频域序列。
按照方程(3)～(9)的形式进行傅里叶变换，经过变换后的方程为：
Figure BDA0003144536260000125
本实施例1中取0～100Hz中等间距50个不同的频率作为k值。
第五步：确定测试载荷及变形样本。
本实施例1中，测试载荷选定为在由一阶弯曲模态和一阶扭转模态组合出的作用力：
Ftest＝∑piφib+∑pjφjt  (17)
模态加权系数pi取值150，模态加权系数pj取值150。将测试载荷加载在有限元模型上并利用Nastran软件SOL400几何非线性静力学分析求解器计算，得到5s的测试变形动响应数据。翼端最大垂直变形271mm，接近机翼翼展的27％。
第六步：求解非线性刚度系数，建立非线性结构降阶模型。
将第五步中得到的位移数据及力数据进行傅里叶变换，得到频域离散序列
Figure BDA0003144536260000131
Qm(k)，
Figure BDA0003144536260000132
fm(k)。
NK个不同的频域k值取为{k1 k2 … k50}，最小二乘问题可表达为：
Figure BDA0003144536260000133
求解未知非线性刚度系数后带入方程(1)中即可得到完整的非线性结构动力学方程。
本发明实施例1的计算结果与基于静力学数据辨识的传统非线性结构降阶模型相比，能够提供更加准确的结构动力学特性信息，在实际应用中具有重要意义。
本发明提供的上述基于动力学响应辨识的大柔性结构降阶模型建模方法，针对大柔性结构建模问题，基于动力学响应数据进行结构非线性刚度系数的辨识，提供一套准确可靠适用于大柔性结构变形分析的几何非线性结构降阶模型。将非线性结构动力学方程以等时间间隔进行离散化，形成广义坐标离散序列，然后对其进行离散傅里叶变换，形成频域变换序列，交换移动方程部分项顺序构成线性方程组，利用模态组合乘比例系数的方式构建测试载荷加载于结构上计算动力学响应，提前动响应数据后求方程组的最小二乘解即可给出非线性刚度系数，建立非线性结构降阶模型。本发明构建的基于动力学响应辨识的大柔性结构降阶模型，与传统的基于静力学响应数据辨识的大柔性结构降阶模型相比，利用结构动力学响应数据样本进行非线性刚度系数的辨识，可以提供更准确的结构动力学特性信息。
显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (5)

1.一种基于动力学响应辨识的大柔性结构降阶模型建模方法，其特征在于，包括如下步骤：
S1：建立大柔性结构的有限元模型，划分有限元模型的结构有限元网格，定义结构有限元网格的单元属性和材料特性，设置计算边界条件；
S2：基于大柔性结构的有限元模型，计算大柔性结构在垂直弯曲方向、扭转方向及水平弯曲方向的线性模态；
S3：考虑包含几何非线性大变形作用的结构动力学方程，将结构动力学方程中的广义坐标和模态力向量以等时间间隔进行离散，得到离散形式的结构动力学方程；
S4：对离散形式的结构动力学方程进行离散傅里叶变换；
S5：选择有限元模型在垂直弯曲方向和扭转方向的线性模态，得到垂直弯曲方向和扭转方向的线性模态组合后的测试载荷，作为外载加载于大柔性结构上进行动力学响应计算，得到结构时域动响应位移数据；
S6：对测试载荷和结构时域动响应位移数据进行傅里叶变换，得到频域离散序列后代入经离散傅里叶变换后的结构动力学方程，利用最小二乘方法计算得到非线性刚度系数后代入步骤S3中未经离散的结构动力学方程，得到完整的非线性结构动力学方程。
2.如权利要求1所述的基于动力学响应辨识的大柔性结构降阶模型建模方法，其特征在于，步骤S3，具体包括：
考虑包含几何非线性大变形作用的结构动力学方程如下：
Figure FDA0003144536250000011
方程(1)对应第m阶结构模态，并应用Einstein求和约定；其中，m＝1,…,N，N表示参与降阶的结构模态数目；Mm表示第m阶结构模态的广义质量项；Km表示第m阶结构模态的广义刚度项；
Figure FDA0003144536250000012
表示第m阶结构模态下的第n阶结构模态广义坐标与第l阶结构模态广义坐标的乘积多项式的非线性刚度项，
Figure FDA0003144536250000021
表示第m阶结构模态下的第n阶结构模态广义坐标、第l阶结构模态广义坐标和第p阶结构模态广义坐标的乘积多项式的非线性刚度系数；qm表示第m阶结构模态的广义坐标，qn表示第n阶结构模态的广义坐标，ql表示第l阶结构模态的广义坐标，qp表示第p阶结构模态的广义坐标；Fm表示第m阶结构模态的模态力向量；
将方程(1)中的广义坐标qm、qn、ql、qp及模态力向量Fm以等时间间隔进行离散，得到离散形式的结构动力学方程：
Figure FDA0003144536250000022
其中，i＝0,1,…,NT-1，表示离散时间点序列；NT表示离散时间点的总数量；qm(i)表示第m阶结构模态的广义坐标离散序列，qn(i)表示第n阶结构模态的广义坐标离散序列，ql(i)表示第l阶结构模态的广义坐标离散序列，qp(i)表示第p阶结构模态的广义坐标离散序列，Fm(i)表示第m阶结构模态的模态力向量离散序列。
3.如权利要求2所述的基于动力学响应辨识的大柔性结构降阶模型建模方法，其特征在于，步骤S4，具体包括：
对离散形式的结构动力学方程中的广义坐标离散序列qm(i)、qn(i)、ql(i)和qp(i)进行离散傅里叶变换，变换后分别表示为Qm(k)、Qn(k)、Ql(k)和Qp(k)，变换关系为：
Figure FDA0003144536250000023
Figure FDA0003144536250000024
Figure FDA0003144536250000025
Figure FDA0003144536250000031
对离散形式的结构动力学方程中的广义坐标二阶项qnql和广义坐标三阶项qnqlqp进行离散傅里叶变换，变换后分别表示为
Figure FDA0003144536250000032
和
Figure FDA0003144536250000033
变换关系为：
Figure FDA0003144536250000034
Figure FDA0003144536250000035
对离散形式的结构动力学方程中的模态力向量Fm(i)进行离散傅里叶变换，变换后表示为fm(k)，变换关系为：
Figure FDA0003144536250000036
将方程(3)～(9)代入方程(2)中，得到经过傅里叶变换后的离散形式的结构动力学方程为：
Figure FDA0003144536250000037
在方程(3)～(10)中取NK个不同的频域k值{k1k2…kNK}，得到NK个方程组成的方程组：
Figure FDA0003144536250000038
4.如权利要求3所述的基于动力学响应辨识的大柔性结构降阶模型建模方法，其特征在于，步骤S5中，得到的垂直弯曲方向和扭转方向的线性模态组合后的测试载荷为：
Ftest＝∑piφib+∑pjφjt            (12)
其中，pi和pj表示人为给定的比例系数，φib表示选定的线性弯曲模态振型，φjt表示选定的线性扭转模态振型。
5.如权利要求4所述的基于动力学响应辨识的大柔性结构降阶模型建模方法，其特征在于，步骤S6中，对测试载荷和结构时域动响应位移数据进行傅里叶变换，得到频域离散序列后代入经离散傅里叶变换后的结构动力学方程，利用最小二乘方法计算得到非线性刚度系数，具体包括：
对测试载荷和结构时域动响应位移数据进行傅里叶变换，得到五个频域离散序列
Figure FDA0003144536250000041
Qm(k)、
Figure FDA0003144536250000042
和fm(k)，代入经离散傅里叶变换后的结构动力学方程组(11)后，利用最小二乘方法计算得到非线性刚度系数
Figure FDA0003144536250000043
和
Figure FDA0003144536250000044
CN202110746615.9A 2021-07-02 2021-07-02 一种基于动力学响应辨识的大柔性结构降阶模型建模方法 Active CN113536456B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110746615.9A CN113536456B (zh) 2021-07-02 2021-07-02 一种基于动力学响应辨识的大柔性结构降阶模型建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110746615.9A CN113536456B (zh) 2021-07-02 2021-07-02 一种基于动力学响应辨识的大柔性结构降阶模型建模方法

Publications (2)

Publication Number Publication Date
CN113536456A CN113536456A (zh) 2021-10-22
CN113536456B true CN113536456B (zh) 2023-05-16

Family

ID=78097564

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110746615.9A Active CN113536456B (zh) 2021-07-02 2021-07-02 一种基于动力学响应辨识的大柔性结构降阶模型建模方法

Country Status (1)

Country Link
CN (1) CN113536456B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2953846A1 (en) * 2013-02-05 2015-12-16 MacGregor Pusnes AS. System and method for transfer of hydrocarbon containing fluids
CN108052772A (zh) * 2017-12-30 2018-05-18 北京航空航天大学 一种基于结构降阶模型的几何非线性静气动弹性分析方法
CN108363843A (zh) * 2018-01-25 2018-08-03 北京航空航天大学 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法
CN108427322A (zh) * 2018-03-16 2018-08-21 南京航空航天大学 一种大柔性飞行器基于在线辨识的建模方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2953846A1 (en) * 2013-02-05 2015-12-16 MacGregor Pusnes AS. System and method for transfer of hydrocarbon containing fluids
CN108052772A (zh) * 2017-12-30 2018-05-18 北京航空航天大学 一种基于结构降阶模型的几何非线性静气动弹性分析方法
CN108363843A (zh) * 2018-01-25 2018-08-03 北京航空航天大学 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法
CN108427322A (zh) * 2018-03-16 2018-08-21 南京航空航天大学 一种大柔性飞行器基于在线辨识的建模方法

Also Published As

Publication number Publication date
CN113536456A (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
CN102012953B (zh) Cfd/csd耦合求解非线性气动弹性仿真方法
Raveh CFD-based models of aerodynamic gust response
CN112580241B (zh) 一种基于结构降阶模型的非线性气动弹性动响应分析方法
Murua 'Flexible aircraft dynamics with a geometrically-nonlinear description of the unsteady aerodynamics'
Dang et al. Accelerated loosely-coupled CFD/CSD method for nonlinear static aeroelasticity analysis
CN110162823B (zh) 考虑气动面曲面效应和法向运动的非定常气动力计算方法
CN110162822B (zh) 耦合结构模态的时域快速非定常气动力计算方法
Komarov et al. Aircraft design using a variable density model
CN113723027A (zh) 一种针对弹性飞机的静气动弹性计算方法
CN109858133B (zh) 一种基于应力映射的点阵结构设计与优化方法
CN113886967A (zh) 多巡航工况的大型飞机机翼气动弹性优化方法
Valente et al. A doublet-lattice method correction approach for high fidelity gust loads analysis
Ritter et al. Comparison of nonlinear aeroelastic methods for maneuver simulation of very flexible aircraft
Zhang et al. A morphing wing with cellular structure of non-uniform density
Ritter et al. Dynamic aeroelastic simulations of the Pazy wing by UVLM with nonlinear viscous corrections
CN112580240B (zh) 一种适用于复杂大柔性飞机建模的非线性子结构方法
Jennings et al. Effect of stabilizer dihedral and static lift on T-tail flutter
CN113536456B (zh) 一种基于动力学响应辨识的大柔性结构降阶模型建模方法
Valente et al. An optimized doublet-lattice method correction approach for a large civil aircraft
Opgenoord et al. Aeroelastic tailoring using additively manufactured lattice structures
CN115795643A (zh) 一种大柔性结构降阶模型快速构建方法
Huo et al. Composite wing elastic axis for aeroelasticity optimization design
Nguyen et al. Multi-point jig twist optimization of mach 0.745 transonic truss-braced wing aircraft and high-fidelity cfd validation
Castells Marin et al. Simulation of Flexible Aircraft Response to Gust and Turbulence for Flight Dynamics Investigations
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