CN105447269B - 一种飞机在定常过载状态下的非结构质量配平方法 - Google Patents

一种飞机在定常过载状态下的非结构质量配平方法 Download PDF

Info

Publication number
CN105447269B
CN105447269B CN201510939210.1A CN201510939210A CN105447269B CN 105447269 B CN105447269 B CN 105447269B CN 201510939210 A CN201510939210 A CN 201510939210A CN 105447269 B CN105447269 B CN 105447269B
Authority
CN
China
Prior art keywords
collection
load
trim
mass
structural
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.)
Expired - Fee Related
Application number
CN201510939210.1A
Other languages
English (en)
Other versions
CN105447269A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical 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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201510939210.1A priority Critical patent/CN105447269B/zh
Publication of CN105447269A publication Critical patent/CN105447269A/zh
Application granted granted Critical
Publication of CN105447269B publication Critical patent/CN105447269B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • 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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种飞机在定常过载状态下的非结构质量配平方法，通过引入非结构质量单元的增量使飞机气动载荷与惯性载荷配平，从而使由无位移约束结构有限元模型分析得到的结果较位移约束方法和惯性释放方法更为合理。本发明通过的缩聚技术得到缩聚后的载荷列阵，从而精确得到飞机结构相对于参考点的非平衡力系；再进一步由飞机过载得到配平载荷，并将其等效地转化为配平后非结构质量单元的增量属性。由此，利用RBE3刚性单元的连接，将这些质量载荷增量传递到结构单元上，即将配平后的质量单元惯性载荷传递到油箱、货仓或设备舱的结构模型承力点上，从而精确得到飞机结构的受力特性。本发明对于提高飞机结构的设计精确性具有重要的应用价值及意义。

Description

一种飞机在定常过载状态下的非结构质量配平方法
技术领域
本发明涉及飞机结构设计领域，具体是一种飞机在定常过载状态下的非结构质量配平方法。
背景技术
飞机在飞行过程中受到各种载荷的作用，影响飞机结构强度的载荷主要有飞行中的空气动力、空气阻力、发动机推力以及结构的惯性载荷等。采用有限元离散化数值分析方法对飞机结构进行设计分析时，会遇到气动力与飞机惯性载荷不平衡的情况。这种不平衡状态可能来自于初始气动载荷与惯性载荷建模的不平衡，也可能是尽管气动载荷与结构惯性载荷在初始条件下是为平衡的，但在飞机结构的优化设计阶段中导致气动载荷与优化后的结构惯性载荷不平衡。飞机结构的优化设计是在满足多种工程约束，如位移约束、应力应变约束或屈曲极限载荷约束等的条件下，通过改变结构部件的尺寸参数，如飞机的蒙皮厚度、梁缘条的截面尺寸参数等，以达到飞机结构质量最轻的设计目的。飞机结构部件的尺寸参数在优化迭代过程中不断变化，则会相应地引起飞机惯性载荷与气动载荷的不平衡。合理有效的配平载荷力系，对于飞机结构的优化设计有着至关重要的影响。
针对上述飞机设计过程中的飞行载荷不平衡问题，现行飞机结构设计中常用的方法主要有两种：一是位移约束方法，另一是惯性释放方法。位移约束方法通常是在机头机尾或翼根处等远离载荷校核区施加固支或简支约束，将原奇异方程转化为非奇异方程，从而得到飞机结构的位移及内力分布。然而这种设置位移约束的方法会将未配平的载荷通过约束节点传递到基础上，造成不恰当的传力路线，并导致约束节点相关单元的应力集中以及变形的不正确，降低了设计分析结果的计算精度。文献1“MSC Nastran Linear StaticAnalysis User’s Guide”公开了一种惯性释放技术。该技术根据达朗伯原理,通过结构载荷列阵的缩聚方法，对飞机结构有限元模型施加与加速度相应的修正惯性载荷，从而配平由气动载荷与初始惯性载荷的差异所引入的非平衡力系；再通过对刚体运动自由度施加约束，进一步得到静态问题的解，从而解决约束不足问题。但该方法与飞机设计的载荷状态相悖，因为飞机的设计状态确定的过载由飞行包线指定，其数值并不会随着不平衡力系的变化而变化。
本发明针对位移约束方法以及惯性释放方法在飞机结构设计应用中如上所述的不足，提出了一种非结构质量配平方法，该方法通过对结构设计过程中的非结构质量载荷配平，达到消除数值模型求解的奇异性问题。数值算例表明该技术更为可靠，不对结构载荷的传力路线产生任何影响，且能够给出结构质量变化时所带来的非结构质量增益，如燃油或设备的改变量。
发明内容
为克服在飞机结构设计中，尤其在结构优化设计过程中，对于飞行状态下因飞机气动载荷与惯性载荷不平衡而导致无位移约束有限元模型的无解问题，本发明提出了一种飞机在定常过载状态下的非结构质量配平方法。
本发明的具体过程是：
步骤1，有限元非结构网格的划分：根据用户给定的飞机CAD数据，采用Hypermesh软件对飞机CAD几何模型进行有限元非结构网格划分。
步骤2，在Hypermesh软件中输入飞机各部位材料的参数。所述的飞机各部位材料的参数包括：杨氏模量E，单位为MPa；剪切模量G，单位为MPa；泊松比μ和密度ρ，单位为Kg/mm3；以及复合材料铺层单向带材料性能；所述的复合材料铺层单向带材料性能包括1-方向杨氏模量E1、2-方向杨氏模量E2、1-2平面泊松比μ12、1-2平面剪切模量G12、1-3平面剪切模量G13、2-3平面剪切模量G23和密度ρ。
步骤3，将建立好以后的机翼有限元模型文件导入Patran软件。
步骤4，建立集中质量点：在油箱承力节点处建立了两个集中质量点，以模拟油箱燃油载荷对飞机结构的作用。在所述飞机中，油箱位于机身段。在飞机有限元设计中，燃油采用集中质量单元代替，通过腹板的作用将其载荷传递至机身。
步骤5，建立非结构配平质量单元：在机身对称面的下端建立一个质量单元。该质量单元x-向位置与油箱承力节点的x-向坐标一致，记为非结构配平质量单元。在有限元初始分析中非结构配平质量单元的质量设置为零，并通过Nastran软件中的RBE3单元连接至步骤4中建立的两个集中质量单元节点上。
步骤6，获取初始刚度矩阵K、初始质量阵M和外部作用的气动载荷列阵{F}：在Patran软件中输入用户给定的飞机的气动载荷，最终生成Nastran软件可执行BDF模型文件，并提交至Nastran软件计算并导出结构的初始刚度矩阵K、初始质量阵M和外部作用的气动载荷列阵{F}。
步骤7，求解飞机有限元结构模型的节点载荷列阵{P}：根据Nastran软件导出的初始质量阵M和外部作用的气动载荷列阵{F}，求解飞机有限元结构模型的节点载荷列阵{P}。节点载荷列阵{P}包括外部作用的气动载荷列阵{F}及结构质量惯性载荷{GI}，结构质量惯性载荷{GI}由设计人员根据当前飞行状态下的载荷系数ny提出。令飞机有限元模型中的节点总数为n，则各节点的加速度向量表示为：
式中，ai＝[0 ay 0 0 0 0]T，ay＝nyg，ay为飞机升力方向的加速度值，g为重力加速度。当飞机匀速平飞时，ny＝1时，飞机的结构质量惯性载荷表示为：
{GI}＝Ma
由上式得到飞机有限元结构模型的节点载荷列阵{P}为：
{P}＝{F}-{GI}
步骤8，选取相对参考点：在机身对称面下端处选取一个有限元节点作为相对参考点，记为r-集。相对参考点x-向位置与油箱承力节点的x-向坐标一致，坐标为(1781.0,-449.964,0.0)。设相对参考点的平动位移及转角位移约束为零。
步骤9，对刚度矩阵进行分块处理：定义飞机有限元模型的完全节点集为a-集，参考节点为r-集，并将a-集中不包含r-集的节点集记为l-集，称为剩余节点集。则飞机结构的刚度矩阵分块表示为
所述Kaa是a-集相关的刚度矩阵，Krr是r-集相关的刚度矩阵，Kll是l-集相关的刚度矩阵，Klr是l-集与r-集的耦合刚度矩阵，Krl是r-集与l-集的耦合刚度矩阵。
步骤10，确定飞机结构模型中a-集与r-集间的约束关系并聚缩载荷列阵：按Guyan聚缩原理，确定飞机结构模型中a-集与r-集间的约束关系并聚缩载荷列阵，得到a-集与r-集间的约束关系[Gar]：
所述ua是a-集内的有限元节点位移，ur是r-集内的有限元节点位移，ul是l-集r-集内的有限元节点位移，Irr是r-集的惯性矩，Grl是r-集与l-集的耦合约束关系；
在确定所述飞机结构模型中a-集与r-集间的约束关系[Gar]时，通过求解所述a-集的结构有限元方程：
得到：
[Klr]{ur}+[Kll]{ul}＝0
{ul}＝-[Kll]-1[Klr]{ur}＝[Glr]{ur}
公式中，所述[Glr]为l-集与r-集间的约束关系矩阵；ur是参考节点r-集的位移向量，ul是参考节点l-集的位移向量。
步骤11，建立减缩载荷列阵{Pr}：由得到的r-集和a-集刚***移约束[Gar]，得到减缩载荷列阵{Pr}：
{Pr}＝[Gar]T{Pa}＝[Gar]T{P} (2)
飞机在重力方向为定常过载状态以及重力方向上的载荷不平衡的飞行状态时，该飞机x-轴方向和z-轴方向受力为零，得到聚缩后刚体飞机的非平衡载荷{Pr}具体形式为：
{Pr}＝[0 Pn 0 M1 0 M3]T (3)
公式中，Pn为飞机沿重力方向载荷，M1为飞机沿x-轴所受弯矩，M3为飞机沿z-轴所受弯矩。
步骤12，由r-集广义载荷确定r-集配平质量阵：令r-集的配平质量阵为r-集节点加速度ar＝[0 ay 0 0 0 0]T，由步骤11中的减缩载荷列阵{Pr}得到r-集的刚体运动方程为
故
由公式(5)给出r-集配平质量阵为：
该质量阵即能够保证飞机在飞行载荷作用下过载状态的不变性，且表明配平质量的增量应为-Pn/ay。
步骤13，确定r-集和非结构配平质量单元节点之间的刚体约束关系：
记相对参考点为r点，非结构配平质量单元节点为s点。相对参考点与非结构配平质量单元节点的转角分别为θr和θs，平移为dr和ds。则非结构配平质量单元节点的平移及转角与相对参考点平移和转角的矩阵表达关系为：
其中，ur＝[dr,θr]T，us＝[ds,θs]T，dx、dy、dz为矢量rsr在x、y、z方向的分量，Grs为相对参考点和非结构配平质量单元节点之间的刚体约束关系矩阵。通过公式(7)能够确定r-集和非结构配平质量单元节点之间的刚体约束关系。
骤14，确定非结构配平质量单元的配平质量阵：由步骤12中的r-集广义配平质量阵等效变换至步骤5中非结构配平质量单元的配平质量阵。令非结构配平质量单元的位移和加速度分别为us和as，配平质量阵为Ms，
步骤15：赋予非结构配平质量单元质量属性，并完成有限元分析：将Ms赋予非结构配平质量单元，通过Nastran求解器即能够进一步完成相应的节点位移及应力分布分析。至此，完成了飞机在定常过载状态下非结构质量的配平。
本发明以飞机载荷典型飞行包线中的恒定过载状态为基本依据，针对结构设计过程中的结构质量不确定而提出的一种非结构质量配平方法，能够实现飞行状态下飞机结构有限元模型的合理分析，从而为飞机结构设计提供有力的技术支持。
本发明的基本特征在于：通过引入非结构质量单元的增量使飞机气动载荷与惯性载荷配平，从而使由无位移约束结构有限元模型分析得到的结果较位移约束方法和惯性释放方法更为合理。其目的在于满足飞机结构优化设计减轻飞机结构质量、增加有效质量载荷的基本目的。本发明中，非结构质量是指飞机中的燃油质量、设备质量或商载质量等。在方法上主要通过载荷列阵的缩聚技术，得到缩聚后的载荷列阵，从而精确得到飞机结构相对于参考点的非平衡力系；再进一步由飞机过载得到配平载荷，并将其等效地转化为配平后非结构质量单元的增量属性。由此，利用RBE3刚性单元的连接，将这些质量载荷增量传递到结构单元上，也即将配平后的质量单元惯性载荷传递到油箱、货仓或设备舱的结构模型承力点上，从而精确得到飞机结构的受力特性。本发明技术对于提高飞机结构的设计精确性具有重要的应用价值及意义。
本发明针对飞机设计过程中的飞行载荷不平衡问题，首先要将其转化为载荷平衡问题，进而才能通过有限元方法进行数值模拟分析。载荷处理方法现常用的方法主要有两种：一是位移约束方法，另一是惯性释放方法。位移约束方法通常是在机头机尾或翼根处等远离载荷校核区施加固支或简支约束，将原奇异方程转化为非奇异方程，从而得到飞机结构的位移及内力分布。然而这种设置位移约束的方法会将未配平的载荷通过约束节点传递到基础上，造成不恰当的传力路线，并导致约束节点相关单元的应力集中以及变形的不正确，降低了设计分析结果的计算精度。文献1“MSC Nastran Linear Static AnalysisUser’s Guide”公开了一种惯性释放技术。该技术根据达朗伯原理,通过结构载荷列阵的缩聚方法，对飞机结构有限元模型施加与加速度相应的修正惯性载荷，从而配平由气动载荷与初始惯性载荷的差异所引入的非平衡力系；再通过对刚体运动自由度施加约束，进一步得到静态问题的解，从而解决约束不足问题。但该方法与飞机设计的载荷状态相悖，因为飞机的设计状态确定的过载由飞行包线指定，其数值并不会随着不平衡力系的变化而变化。
本发明针对前述两种方法的弊端，其目的在于两点：一是不会对结构造成不恰当的传力路线；二是不会改变飞机的过载状态。本发明通过引入非结构配平质量配平单元，并由非结构配平质量单元在给定过载状态下的惯性载荷对飞机结构设计过程中的非平衡惯性载荷进行配平，使得非结构配平质量单元的惯性载荷与飞机的惯性载荷关于相对参考点的作用之和为零。此技术保证了飞机在特定外部载荷的作用下过载状态的不变性，这与飞机设计过程中过载状态不变的要求相符合。
数值算例表明该技术能够保证飞机的过载状态在非平衡结构载荷作用下不发生改变，且不对结构载荷的传力路线产生任何影响。此外，由非结构配平质量单元的质量能够给出结构质量变化时所带来的非结构质量增益，如燃油或设备的改变量等，这对飞机结构设计有着重要的意义。
本发明适用于飞行载荷状态下，飞机结构设计过程中质量变化时的惯性载荷配平计算问题，更具体地，涉及飞机在重力方向为定常过载状态时，结构质量设计迭代过程中有限元模型的惯性质量载荷配平计算，实现了飞行状态下飞机结构有限元模型的合理分析，从而为飞机结构设计提供有力的技术支持。
附图说明
图1是无人机模型外形图；
图2是无人机模型骨架及机体坐标系示意图；
图3是无人机油箱承力节点位置示意图；
图4是无人机非结构配平质量单元模型图；
图5是无人机相对参考点示意图；
图6时本发明的流程图。
附图中，1.油箱承力节点，2.油箱承力节点位置，3.相对参考点。
具体实施方式
本实施例通过对某型无人机在定常过载状态下的非结构质量配平过程对本发明的技术方案进行详细描述。
本实施例具体是针对某用户给定的无人机CAD数据、过载系数及气动载荷，通过非结构质量配平方法完成该无人机在自由飞行状态下的受力分析设计过程。具体过程包括以下步骤：
步骤1，有限元非结构网格的划分。根据用户给定的无人机CAD数据，采用Hypermesh软件对无人机CAD几何模型进行有限元非结构网格划分。无人机几何外形如图1所示，内部骨架如图2所示。该无人机分为机身、内翼、外翼三段。机身长为3.893m，内翼的主要承力构件包括3个腹板和5个翼肋；外翼的主要承力构件包括3个腹板和16个翼肋。机翼主要参数如下：展长为15.000m，根弦长为2.791m，翼尖弦长为0.900m，前缘后掠角为28°，后缘后掠角为24°，机翼无扭转角，蒙皮厚度为0.002m。所述的无人机几何模型中，所有的梁橼条、肋橼条、加强筋和长桁均使用两节点均直梁单元模拟，共有12334个单元。其中梁橼条、肋橼条为“L”型梁单元，加强筋与长桁使用矩形梁单元。所有蒙皮、梁腹板、肋腹板和框腹板均使用壳单元模拟，壳单元共有63864个。在无人机各个接头处采用Nastran软件中的RBE2单元模拟。整个无人机共有节点59802个，单元76199个。
步骤2，在Hypermesh软件中输入用户指定的无人机材料参数。该无人机共采用5种材料：发动机喷口处使用耐高温材料30CrMnSi；辅助梁橼条与腹板采用7050铝合金材料；内翼蒙皮、外翼蒙皮、梁腹板以及机身蒙皮使用对称铺层复合材料；其余部分材料均使用2024铝合金。各材料属性如表1，2所示。表2中，E1为1-方向杨氏模量，E2为2-方向杨氏模量，μ12为1-2平面泊松比，G12为1-2平面剪切模量，G13为1-3平面剪切模量，G23为2-3平面剪切模量，ρ为密度。
表1模型中金属材料性能
表2复合材料铺层单向带材料性能
步骤3，将建立的机翼有限元模型文件导入Patran软件。
步骤4，建立集中质量点。在油箱承力节点1处建立了两个集中质量点。这两个集中质量点用于模拟油箱燃油载荷对飞机结构的作用。在所述无人机中，油箱位于机身段。在无人机有限元设计中，燃油通常采用集中质量单元代替，通过腹板的作用将所述燃油载荷传递至机身。本实例中，根据设计，两个油箱承力节点位置2分别为(1781.0,43.7799,400.0)和(1781.0,43.7799,-400.0)，如图3所示。
步骤5，建立非结构配平质量单元。在机身对称面的下端建立一个质量单元，该质量单元x-向位置与油箱承力节点1的x-向坐标一致，并把该质量单元记为非结构配平质量单元。通过确定该非结构配平质量单元的质量特性，以达到在无人机过载状态不改变的条件下无人机结构载荷的自平衡。本实例中，所述非结构配平质量单元的坐标为(1781.0,-643.78,0.0)。在有限元初始分析中，所述非结构配平质量单元质量为零，并通过Nastran软件中的RBE3单元连接至步骤4中建立的两个集中质量单元节点上，也就是油箱的承力接头处，如图4所示。所述RBE3单元是Nastran软件中的连接单元，作用是实现载荷从非结构配平质量单元到油箱承力节点1的传递。
步骤6，获取初始刚度矩阵K、初始质量阵M和气动载荷列阵{F}。在Patran软件中输入用户给定的无人机的气动载荷，并提交至Nastran软件计算并导出结构的初始刚度矩阵K、初始质量阵M和气动载荷列阵{F}。
步骤7，求解无人机有限元结构模型的节点载荷列阵{P}。根据Nastran软件导出的初始质量阵M和用户给定的气动载荷{F}，求解无人机有限元结构模型的节点载荷列阵{P}。节点载荷列阵{P}由外部作用的气动载荷{F}及结构质量惯性载荷{G}两项构成，质量惯性载荷{G}由设计人员根据当前飞行状态下的载荷系数ny提出。令无人机有限元模型中的节点总数为n，则各节点的加速度向量表示为：
式中，ai＝[0 ay 0 0 0 0]T，ay＝nyg，ay为无人机升力方向的加速度值，g为重力加速度，即g＝-9800mm/s2。本实例中，无人机为匀速平飞，ny＝1。于是，无人机的惯性载荷表示为：
G＝Ma
由上式得到无人机有限元结构模型的节点载荷列阵{P}为：
{P}＝{F}-{G}
步骤8，选取相对参考点。在机身对称面下端处选取一个有限元节点作为相对参考点3，相对参考点x-向位置与油箱承力节点的x-向坐标一致，记为r-集。所述相对参考点3的坐标为(1781.0,-449.964,0.0)，位置如图5所示。设相对参考点的平动位移及转角位移约束为零。相对参考点3表示结构模型中其他节点的变形位移均相对于该相对参考点，且该相对参考点的形变位移为零。由于无人机在自由飞行状态下刚度矩阵奇异，故无人机在结构载荷作用下该节点的位移包含3个平动刚***移和3个转动刚***移。参考点的选取不同并不会对节点间的相对位移及飞行器受力特性产生影响。
步骤9：对刚度矩阵进行分块处理。定义无人机有限元模型的完全节点集为a-集、相对参考点为r-集，并记a-集中不包含r-集的节点集为l-集，称为剩余节点集。则无人机结构的刚度矩阵K的分块表示为
这里Kaa是a-集相关的刚度矩阵，Krr是r-集相关的刚度矩阵，Kll是l-集相关的刚度矩阵，Klr是l-集与r-集的耦合刚度矩阵，Krl是r-集与l-集的耦合刚度矩阵。
步骤10：确定无人机结构模型中a-集与r-集间的约束关系并聚缩载荷列阵。按Guyan聚缩原理，确定无人机结构模型中a-集与r-集间的约束关系并聚缩载荷列阵。获此关系的关键是求解一个无载荷作用的有限元方程，故a-集的结构有限元方程写成如下分块形式：
式中，ur是参考节点r-集的位移向量，ul是参考节点l-集的位移向量。分块求解上式的第二块方程组得到：
[Klr]{ur}+[Kll]{ul}＝0
{ul}＝-[Kll]-1[Klr]{ur}＝[Glr]{ur}
公式中，[Glr]为l-集与r-集间的约束关系矩阵。于是，得到a-集与r-集间的约束关系[Gar]：
步骤11，建立减缩载荷列阵{Pr}。由得到的r-集和a-集刚***移约束[Gar]，得到减缩载荷列阵{Pr}：
{Pr}＝[Gar]T{Pa}＝[Gar]T{P} (2)
本发明只针对无人机在重力方向为定常过载状态以及重力方向上的载荷不平衡情况。针对这一常见的飞行状态，无人机在x-向方向和z-向方向受力为零，且由于无人机仅受重力方向的载荷作用，因此刚体无人机绕竖直方向转动的扭矩为零。故得到聚缩后刚体无人机的非平衡载荷{Pr}具体形式为：
{Pr}＝[0 Pn 0 M1 0 M3]T (3)
这里，Pn为无人机沿重力方向载荷，M1为无人机沿x-轴所受弯矩，M3为无人机沿z-轴所受弯矩。本实例中，将无人机气动载荷代入(2)式，得到减缩载荷列阵{Pr}为
{Pr}＝[0 1.86×105 0 7.30×105 0 8.79×106]T (4)
步骤12：由r-集广义载荷确定r-集配平质量阵。令r-集的配平质量阵为r-集节点加速度ar＝[0 ay 0 0 0 0]T，为保证结构的过载状态不变，即无人机载荷关于相对参考点作用为零。则r-集配平质量的惯性载荷与无人机结构的不平衡载荷{Pr}关于r-集作用为零，由步骤11中的减缩载荷列阵{Pr}能够得到r-集的刚体运动方程为
故
由上式即能够给出一个可行的r-集配平质量阵解：
广义质量增量阵即能够保证无人机在飞行载荷作用下过载状态的不变性，且表明配平质量的增量应为-Pn/ay。
将(4)式和ay代入(6)式能够得关于r-集配平质量阵
需特别指出，在本发明方法中，相对参考点3的节点位置是固定的，从而导致非配平质量阵是广义的。事实上在空间中存在着一点，在该点上采用传统的质量阵，即仅在对角线元素中的前三个有值，且大小为-Pn/ay，即能够实现为M1及M3平衡。令该点位置相对于相对参考点r的x向距离为rx，z向距离为rz，则
rx＝M3/Pn，rz＝M1/Pn
所述的rx与rz主要用于结构优化后对于油箱等货载结构进行空间位置调整，通过对其重心的调整，以消除优化过程中由于尺寸更新带来的无人机重心位置变化的影响。根据得知无人机在当前构型下，需增加载荷质量为18.98Kg，且该质量重心相对于相对参考点x向距离为47.25mm，z向距离为3.92mm。
步骤13：确定相对参考点r-集和非结构配平质量单元节点之间的刚体约束关系。记相对参考点为r点，非结构配平质量单元节点为s点。相对参考点与非结构配平质量单元节点的转角分别为θr和θs，位移分别为dr和ds。非结构配平质量单元节点的位移和转角用相对参考点的位移和转角表示为：
dr＝ds+θs×rsr
θr＝θs
式中，rsr为由节点s到节点r的矢量。则非结构配平质量单元节点的位移及转角与相对参考点位移和转角的矩阵表达关系为：
其中，ur＝[dr,θr]T，us＝[ds,θs]T，dx、dy、dz为矢量rsr在x、y、z方向的分量，Grs为相对参考点和非结构配平质量单元节点之间的刚体约束关系矩阵。将相对参考点坐标(1781.0,-643.78,0.0)和非结构配平质量单元节点坐标(1781.0,-448.96,0.0)代入(7)式得到相对参考点和非结构配平质量单元节点之间的刚体约束关系Grs
骤14，确定非结构配平质量单元的配平质量阵。由步骤12中的r-集配平质量阵等效变换至步骤5中非结构配平质量单元的配平质量阵。令非结构配平质量单元的位移和加速度分别为us和as，配平质量阵为Ms。则非结构配平质量单元惯性力在相应的位移上所做的功Ws为
同理能够得r-集配平质量阵的惯性载荷在相应的位移上所做的功Wr为
由于等效转换两者能量一致，即Ws＝Wr。故得
将步骤13得到的r-集和非结构配平质量单元节点之间的刚体约束关系Grs，代入(8)式即能够得到非结构配平质量单元质量特性Ms
上式中的Ms即为结构模型中所加非结构质量点的配平质量阵，应当注意的是，该质量配平增量阵是广义的，不仅包含质量信息，还包含了转动惯量信息。Ms反映了为保证无人机在飞行载荷作用下过载状态不变而需增加的惯性力及力矩，这是因为设计过程中，结构构形参数变化所导致的质心变化是不定的，而非结构配平质量点的设置是确定的，即惯性载荷的配平不仅只是质量力的配平，也需力矩的配平。
步骤15：赋予非结构配平质量单元质量属性，并完成有限元分析。赋予非结构配平质量单元质量属性，并完成有限元分析。将由步骤14计算所得的Ms赋予非结构配平质量单元，通过NASTRAN求解器即能够进一步完成相应的节点位移及应力分布分析。
通过该非结构配平质量单元的作用，能够得到一个自平衡的无人机结构有限元模型，即在相对参考点处支反力为零。本发明在不改变过载状态的情况下，通过非结构质量的配平得到一个自平衡的无人机结构，这对于无人机的有限元数值模拟是至关重要的。如果采用传统的位移约束法或惯性释放法则会使得有限元模拟得到不合理的数值解。
本实施例中，步骤7～14主要完成不平衡载荷及非结构配平质量单元质量特性分析，通过非结构配平质量单元的引入，整体保持无人机在设计过程中的质量特性的不变，包括无人机的质量大小及无人机的重心位置不变，从而保证在特定外部载荷的作用下，无人机的过载不变。该技术满足了工程设计需求，且实施过程形式简单，能够方便的通过编写计算机程序自动实现非结构配平质量单元的质量特性分析，从而大大加快无人机有限元的优化设计分析。该技术克服了位移约束方法及惯性释放方法引入的不合理分析结果，为高效的无人机结构设计，尤其是高效的无人机结构优化设计，提供良好的技术支持。

Claims (2)

1.一种飞机在定常过载状态下的非结构质量配平方法，其特征在于，具体过程是：
步骤1，有限元非结构网格划分：根据用户给定的无人机CAD数据，采用Hypermesh软件对无人机CAD几何模型进行有限元非结构网格划分；
步骤2，在Hypermesh软件中输入无人机各部位材料的参数：所述的无人机各部位材料的参数包括：杨氏模量E，单位为MPa；剪切模量G，单位为MPa；泊松比μ和密度ρ，单位为Kg/mm3；以及复合材料铺层单向带材料性能；所述的复合材料铺层单向带材料性能包括1-方向杨氏模量E1、2-方向杨氏模量E2、1-2平面泊松比μ12、1-2平面剪切模量G12、1-3平面剪切模量G13、2-3平面剪切模量G23和密度ρ；
步骤3，将建立好以后的机翼有限元模型文件导入Patran软件；
步骤4，建立集中质量点：在油箱承力节点处建立了两个集中质量点，以模拟油箱燃油载荷对飞机结构的作用；在所述无人机中，油箱位于机身段；在飞机有限元设计中，燃油采用集中质量单元代替，通过腹板的作用将其载荷传递至机身；
步骤5，建立质量单元：在机身对称面的下端建立一个质量单元；该质量单元x-向位置与油箱承力节点的x-向坐标一致，记为非结构配平质量单元；在有限元初始分析中非结构配平质量单元的质量设置为零，并通过Nastran软件中的RBE3单元连接至步骤4中建立的两个集中质量单元节点上；
步骤6，获取初始刚度矩阵K、初始质量阵M和外部作用的气动载荷列阵{F}：在Patran软件中输入用户给定的无人机的气动载荷，最终生成Nastran软件可执行BDF模型文件，并提交至Nastran软件计算并导出结构的初始刚度矩阵K、初始质量阵M和外部作用的气动载荷列阵{F}；
步骤7，求解飞机有限元结构模型的节点载荷列阵{P}：根据Nastran软件导出的初始质量阵M和外部作用的气动载荷列阵{F}，求解飞机有限元结构模型的节点载荷列阵{P}；节点载荷列阵{P}包括外部作用的气动载荷列阵{F}及结构质量惯性载荷{GI}，结构质量惯性载荷{GI}由设计人员根据当前飞行状态下的载荷系数ny提出；令飞机有限元模型中的节点总数为n，则各节点的加速度向量表示为：
式中，ai＝[0 ay 0 0 0 0]T，ay＝nyg，ay为飞机升力方向的加速度值，g为重力加速度；当飞机匀速平飞时，ny＝1时，结构质量惯性载荷表示为：
{GI}＝Ma
由上式得到飞机有限元结构模型的节点载荷列阵{P}为：
{P}＝{F}-{GI}
步骤8，选取相对参考点：在机身对称面下端处选取一个有限元节点作为相对参考点，记为r-集；相对参考点x-向位置与油箱承力节点的x-向坐标一致，坐标为(1781.0,-449.964,0.0)；设相对参考点的平动位移及转角位移约束为零；
步骤9，对刚度矩阵进行分块处理：定义飞机有限元模型的完全节点集为a-集，参考节点为r-集，并将a-集中不包含r-集的节点集记为l-集，称为剩余节点集；则飞机结构的刚度矩阵分块表示为
所述Kaa是a-集相关的刚度矩阵，Krr是r-集相关的刚度矩阵，Kll是l-集相关的刚度矩阵，Klr是l-集与r-集的耦合刚度矩阵，Krl是r-集与l-集的耦合刚度矩阵；
步骤10，确定飞机结构模型中a-集与r-集间的约束关系并聚缩载荷列阵：按Guyan聚缩原理，确定飞机结构模型中a-集与r-集间的约束关系并聚缩载荷列阵，得到
a-集与r-集间的约束关系[Gar]：
所述ua是a-集内的有限元节点位移，ur是r-集内的有限元节点位移，ul是l-集r-集内的有限元节点位移，Irr是r-集的惯性矩，Grl是r-集与l-集的耦合约束关系；
步骤11，建立聚缩后的非平衡载荷列阵{Pr}：由得到的r-集和a-集刚***移约束[Gar]，得到聚缩后的非平衡载荷列阵{Pr}：
{Pr}＝[Gar]T{Pa}＝[Gar]T{P} (2)
飞机在重力方向为定常过载状态以及重力方向上的载荷不平衡的飞行状态时，该飞机x-轴方向和z-轴方向受力为零，得到聚缩后的非平衡载荷列阵{Pr}具体形式为：
{Pr}＝[0 Pn 0 M1 0 M3]T (3)
公式中，Pn为飞机沿重力方向载荷，M1为飞机沿x-轴所受弯矩，M3为飞机沿z-轴所受弯矩；
步骤12，由r-集广义载荷确定r-集配平质量阵：令r-集的配平质量阵为r-集节点加速度ar＝[0 ay 0 0 0 0]T，由步骤11中的减缩载荷列阵{Pr}得到r-集的刚体运动方程为
故
由公式(5)给出r-集配平质量阵为：
该质量阵即能够保证飞机在飞行载荷作用下过载状态的不变性，且表明配平质量的增量应为-Pn/ay；
步骤13，确定r-集和非结构配平质量单元节点之间的刚体约束关系：
记相对参考点为r点，非结构配平质量单元节点为s点；相对参考点与非结构配平质量单元节点的转角分别为θr和θs，平移为dr和ds；则非结构配平质量单元节点的平移及转角与相对参考点平移和转角的矩阵表达关系为：
其中，ur＝[dr,θr]T，us＝[ds,θs]T，dx、dy、dz为矢量rsr在x、y、z方向的分量，Grs为相对参考点和非结构配平质量单元节点之间的刚体约束关系矩阵；通过公式(7)能够确定r-集和非结构配平质量单元节点之间的刚体约束关系；
步骤14，确定非结构配平质量单元的配平质量阵：由步骤12中的r-集广义配平质量阵等效变换至步骤5中非结构配平质量单元的配平质量阵；令非结构配平质量单元的位移和加速度分别为us和as，配平质量阵为Ms，
步骤15，赋予非结构配平质量单元质量属性，并完成有限元分析：将Ms赋予非结构配平质量单元，通过Nastran求解器即能够进一步完成相应的节点位移及应力分布分析；至此，完成了飞机在定常过载状态下非结构质量的配平。
2.如权利要求1所述飞机在定常过载状态下的非结构质量配平方法，其特征在于，在确定所述飞机结构模型中a-集与r-集间的约束关系[Gar]时，通过求解所述a-集的结构有限元方程：
得到：
[Klr]{ur}+[Kll]{ul}＝0
{ul}＝-[Kll]-1[Klr]{ur}＝[Glr]{ur}
公式中，所述[Glr]为l-集与r-集间的约束关系矩阵；ur是参考节点r-集的位移向量，ul是参考节点l-集的位移向量。
CN201510939210.1A 2015-12-15 2015-12-15 一种飞机在定常过载状态下的非结构质量配平方法 Expired - Fee Related CN105447269B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510939210.1A CN105447269B (zh) 2015-12-15 2015-12-15 一种飞机在定常过载状态下的非结构质量配平方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510939210.1A CN105447269B (zh) 2015-12-15 2015-12-15 一种飞机在定常过载状态下的非结构质量配平方法

Publications (2)

Publication Number Publication Date
CN105447269A CN105447269A (zh) 2016-03-30
CN105447269B true CN105447269B (zh) 2019-01-11

Family

ID=55557440

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510939210.1A Expired - Fee Related CN105447269B (zh) 2015-12-15 2015-12-15 一种飞机在定常过载状态下的非结构质量配平方法

Country Status (1)

Country Link
CN (1) CN105447269B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106874572B (zh) * 2017-01-21 2020-06-19 西安交通大学 一种考虑串油特性的飞机油箱承载结构轻量化设计方法
CN107145646B (zh) * 2017-04-19 2021-07-09 中国航空工业集团公司西安飞机设计研究所 一种计算飞机尾撬支反力的方法
CN107844124B (zh) * 2017-12-01 2020-07-14 吉林大学 一种四旋翼飞行器携带不平衡负载稳定飞行的控制方法
CN110197033B (zh) * 2019-05-30 2023-12-01 上海海事大学 一种含连接单元的柔性仿生翼有限元建模方法
CN110941920B (zh) * 2019-09-10 2022-09-27 厦门大学 一种用于无人机飞行载荷数据计算与后处理的方法
CN110704950B (zh) * 2019-09-27 2020-07-31 西北工业大学 一种消除自由飞配平载荷下飞机变形中刚性位移的方法
CN113704881B (zh) * 2021-08-07 2023-08-04 中国航空工业集团公司沈阳飞机设计研究所 一种结构有限元模型惯性载荷施加方法
CN114065394B (zh) * 2021-11-19 2024-07-19 中国直升机设计研究所 一种直升机机体主承力结构应力分析方法
CN118153211B (zh) * 2024-05-10 2024-07-16 中国航发四川燃气涡轮研究院 航空发动机极限载荷确定方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2936882A1 (fr) * 2008-10-07 2010-04-09 Thales Sa Aeronef pourvu d'un manche de commande, d'une gouverne de profondeur et d'un systeme de pilotage automatique
CN103983463A (zh) * 2014-04-17 2014-08-13 中国航空工业集团公司沈阳飞机设计研究所 一种飞机机体与起落架联合加载的验证试验方法
CN104408231A (zh) * 2014-11-04 2015-03-11 江西洪都航空工业集团有限责任公司 一种小飞机全机载荷配平方法
CN104931250A (zh) * 2015-06-29 2015-09-23 中国航空工业集团公司西安飞机设计研究所 一种高升力***全机加载动态试验方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2936882A1 (fr) * 2008-10-07 2010-04-09 Thales Sa Aeronef pourvu d'un manche de commande, d'une gouverne de profondeur et d'un systeme de pilotage automatique
CN103983463A (zh) * 2014-04-17 2014-08-13 中国航空工业集团公司沈阳飞机设计研究所 一种飞机机体与起落架联合加载的验证试验方法
CN104408231A (zh) * 2014-11-04 2015-03-11 江西洪都航空工业集团有限责任公司 一种小飞机全机载荷配平方法
CN104931250A (zh) * 2015-06-29 2015-09-23 中国航空工业集团公司西安飞机设计研究所 一种高升力***全机加载动态试验方法

Also Published As

Publication number Publication date
CN105447269A (zh) 2016-03-30

Similar Documents

Publication Publication Date Title
CN105447269B (zh) 一种飞机在定常过载状态下的非结构质量配平方法
Kenway et al. Aerostructural optimization of the Common Research Model configuration
Murua et al. Assessment of wake-tail interference effects on the dynamics of flexible aircraft
Cavagna et al. NeoCASS: an integrated tool for structural sizing, aeroelastic analysis and MDO at conceptual design level
Cavagna et al. Structural sizing, aeroelastic analysis, and optimization in aircraft conceptual design
CN102867097B (zh) 计及静弹性变形影响的光固化快速成型风洞模型设计方法
Meddaikar et al. Aircraft aeroservoelastic modelling of the FLEXOP unmanned flying demonstrator
Murua Flexible aircraft dynamics with a geometrically-nonlinear description of the unsteady aerodynamics
Xie et al. Static aeroelastic analysis including geometric nonlinearities based on reduced order model
Dussart et al. Flight dynamic modelling and simulation of large flexible aircraft
Cho et al. Combined co-rotational beam/shell elements for fluid–structure interaction analysis of insect-like flapping wing
Ritter et al. Comparison of nonlinear aeroelastic methods for maneuver simulation of very flexible aircraft
CN103853890A (zh) 一种高超声速飞行器气动弹性剪裁方法
Fugate et al. Aero-Structural Modeling of the Truss-Braced Wing Aircraft Using Potential Method with Correction Methods for Transonic Viscous Flow and Wing-Strut Interference Aerodynamics
Farsadi et al. High aspect ratio composite wings: geometrically nonlinear aeroelasticity, multi-disciplinary design optimization, manufacturing, and experimental testing
Klimmek et al. Aircraft loads-an important task from pre-design to loads flight testing
Guimarães et al. Flexibility Assessment of the Aeroelastic-flight-dynamic Behavior for Supersonic Aircraft
Laughlin et al. A parametric and physics-based approach to structural weight estimation of the hybrid wing body aircraft
Reed et al. Coaxial Rotor Wake and Prop Induction Impact on a Horizontal Tail Using HPCMP CREATE™-AV Helios
Zhao et al. Fuel Weight Minimization for Large N+ 3 Composite Transports with Multiple Control Surfaces
Corcione et al. Impact at aircraft level of elastic efficiency of a forward-swept tailplane
Nguyen et al. Multi-point jig twist optimization of mach 0.745 transonic truss-braced wing aircraft and high-fidelity cfd validation
Ramos Construction and analysis of a lightweight UAV wing prototype
Gupta Flight Dynamics Model of a Small Flexible Aircraft
Xiong et al. Jig Twist Optimization of Mach 0.745 Transonic Truss Braced Wing Aircraft and High-Fidelity CFD Validation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190111

Termination date: 20211215

CF01 Termination of patent right due to non-payment of annual fee