一种基于疲劳强度的无人机结构设计方法与流程

本发明属于无人机结构设计领域，具体涉及一种减少计算量，提高设计效率的基于疲劳强度的无人机结构设计方法。

背景技术：

疲劳损毁是无人机结构失效的主要模式之一，对于无人机机体结构中出现的疲劳裂纹，一旦没有及时维修，导致裂纹扩展将导致无人机机体的毁灭性灾难。目前世界上主要的无人机结构设计标准都规定了疲劳度校核的相关指标。事实上，无人机结构的初始设计并非均能满足结构疲劳强度设计使用寿命的要求，在评估实际的机体节点结构疲劳强度时，需要对无法满足设计标准和要求的结构进行优化，通常的做法是根据结构优化方案不断修改机体结构的有限元模型，再重复进行谱分析计算。然而，由于无人机机体结构疲劳损伤程度大的位置往往是结构很复杂的部位，而且谱分析计算过程复杂、计算量十分庞大，这为无人机的节点结构节点优化设计带来很大的困难。

技术实现要素：

本发明的目的在于提出一种减少计算量，提高设计效率的基于疲劳强度的无人机结构设计方法。

一种基于疲劳强度的无人机结构设计方法，包括如下步骤：

(1)输入无人机机体结构设计模型、待计算的风向u和风向的频率ω，将无人机机体结构设计模型转化为有限元分析模型；

(2)评估无人机机体结构的有限元分析模型的疲劳强度，确认无人机机体结构的易疲劳热点；

(3)向无人机机体结构的易疲劳热点施加m种不同角度的风向，确定每种风向条件下易疲劳热点的疲劳累积损伤度bi，i∈(0,m)，i为风向标号；选择bi值最大的风向，设置为最危险风向；最危险风向对易疲劳热点的疲劳累积损伤度为bmax，；

(4)计算bmax与总疲劳累积损伤度b的确定性关系系数q；

(5)将最危险风向持续施加在无人机机体结构的易疲劳热点上，连续调整最危险风向的频率，采集并记录易疲劳热点的疲劳累积损伤度，选择对应易疲劳热点受到的最大应力的频率，设置为易疲劳热点的最危险频率，并获取该易疲劳热点应力的最大值τmax；

(6)计算易疲劳热点的应力与疲劳累积损伤度的关系p；

其中，n为设定的常数；

(7)计算对易疲劳热点施加频率为ω的风向u时，易疲劳热点处的应力τhot；

(8)计算对易疲劳热点施加频率为v的风向u时，易疲劳热点处的总疲劳累积损伤度bu；

所述的疲劳累积损伤度的确定方法为：

(3.1)采集无人机航率系数α；

(3.2)确认无人机疲劳计算的回复期tr，tr＝4.7×10⁸s；

(3.3)根据无人机结构疲劳标准选取盈利疲劳s-n曲线参数和n；

(3.4)根据n构建伽玛函数：

(3.5)确定在第n1个装载、风向i、航向j的情况下的应力响应谱的零阶矩m0ijn1和二阶矩m2ijn1；

(3.6)确定第n1个装载的时间分配系数pn1、第i个风向出现的概率pi和第j个航向出现的概率pj；

(3.7)确定在第n1个装载、风向i、航向j的情况下的应力响应过零率νijn1

(3.8)计算无人机遭遇频率ωek下航向角为θ时易疲劳热点的波能谱gηη(ωek,θ)

(3.9)计算疲劳累积损伤度：

所述应力的计算方法包括：

根据易疲劳热点周边的应力梯度，根据网格模型提取距离热点3t/2和t/2处的表面应力，t为风向周期，根据线性关系外插得到易疲劳热点处的应力，提取不同频率的具有单位波幅的规则波产生的各项载荷的实部τc和虚部τs，得到无人机机体结构在对应频率下的合成应力τ＝τc+iτs，选取易疲劳热点附近的四个有限元网格单元，计算得到的各单元应力τkx，τky，k＝1,2,3,4，为有限元网格标号，热点处的应力τhotx、τhoty为：

τhot的单位为n/mm²。

所述的计算无人机遭遇频率下航向角为时易疲劳热点的波能谱，包括如下步骤：

(3.8.1)确定平均空气波压力hz和平均跨零周期tz；

(3.8.2)测量无人机航速u；

(3.8.3)计算风波谱

(3.8.4)计算波能谱

所述的易疲劳热点的平均跨零周期tz为：

所述的平均空气波压力hz通过当地航空资料查得。

本发明的有益效果在于：

本发明通过无人机机体的典型节点结构在某一风向的特定计算频率下，通过确定无人机结构的热点应力与此风向下结构疲劳累积损伤度的关系，以及风向疲劳累积损伤度与疲劳累积总损伤度的关系，提供了一种基于疲劳强度的无人机机体结构设计方法。本发明可以有效解决由于机体易疲劳位置结构复杂，因谱分析方法计算量庞大，从而导致降低设计容错率的技术问题。本发明为通过谱分析方法开展无人机结构疲劳度设计的方法带来了极大减少设计工作量、提高设计效率的技术效果。

附图说明

图1为无人机机体结构谱分析流程图；

图2为无人机机体疲劳优化简化方法流程图；

图3为无人机机体易疲劳热点的热点应力插值方式；

图4为易疲劳热点180度风向的曲线图。

具体实施方式

下面结合附图对本发明做进一步描述。

本发明涉及无人机机体结构优化设计领域，特别是涉及基于疲劳强度评估结果的无人机机体的优化设计方法。本方法使无人机机体的典型节点在某一风向的特定频率下，通过结构热点应力与风向的疲劳累积损伤度的关系，风向疲劳累积损伤度与疲劳累积总损伤度的关系，给出一种基于疲劳强度的无人机结构优化简化计算方法。该方法有效的解决了由于无人机机体结构疲劳损伤较大的复杂部位计算过程复杂、任务繁重的问题，极大的减少了设计工作量，提高了设计效率，为无人机机体结构的优化设计带来极大便利。

实施例1

研究对象：目标无人机具有四旋翼结构，且其在机翼与本体连接处附近焊接，造成无人机翼内端部与主体连接处存在明显的应力集中。在对目标无人机的典型节点结构进行疲劳强度评估时，发现图3标注位置的疲劳热点的疲劳累积损伤度过大，同时，通过图4研究发现：该结构的疲劳累积损伤度与疲劳累积总损伤度存在确定关系。因此基于疲劳强度的无人机结构设计方法，包括如下步骤：

(1)输入无人机机体结构设计模型、待计算的风向u和风向的频率ω，将无人机机体结构设计模型转化为有限元分析模型；

(2)评估无人机机体结构的有限元分析模型的疲劳强度，确认无人机机体结构的易疲劳热点；

(3)向无人机机体结构的易疲劳热点施加m种不同角度的风向，确定每种风向条件下易疲劳热点的疲劳累积损伤度bi，i∈(0,m)，i为风向标号；选择bi值最大的风向，设置为最危险风向；最危险风向对易疲劳热点的疲劳累积损伤度为bmax；

所述的疲劳累积损伤度的确定方法为：

(3.1)采集无人机航率系数α；

(3.2)确认无人机疲劳计算的回复期tr，tr＝4.7×10⁸s；

(3.3)根据无人机结构疲劳标准选取盈利疲劳s-n曲线参数和n；

(3.4)根据n构建伽玛函数：

(3.5)确定在第n1个装载、风向i、航向j的情况下的应力响应谱的零阶矩m0ijn1和二阶矩m2ijn1；

(3.6)确定第n1个装载的时间分配系数pn1、第i个风向出现的概率pi和第j个航向出现的概率pj；

(3.7)确定在第n1个装载、风向i、航向j的情况下的应力响应过零率νijn1

(3.8)计算无人机遭遇频率ωek下航向角为θ时易疲劳热点的波能谱gηη(ωek,θ)；

所述的计算无人机遭遇频率下航向角为时易疲劳热点的波能谱，包括如下步骤：

(3.8.1)确定平均空气波压力hz和平均跨零周期tz；

(3.8.2)测量无人机航速u；

(3.8.3)计算风波谱

所述的易疲劳热点的平均跨零周期tz为：

所述的平均空气波压力hz通过当地航空资料查得。

(3.8.4)计算波能谱

(3.9)计算疲劳累积损伤度：

(4)计算bmax与总疲劳累积损伤度b的确定性关系系数q；

(5)将最危险风向持续施加在无人机机体结构的易疲劳热点上，连续调整最危险风向的频率，采集并记录易疲劳热点的疲劳累积损伤度，选择对应易疲劳热点受到的最大应力的频率，设置为易疲劳热点的最危险频率，并获取该易疲劳热点应力的最大值τmax；所述应力的计算方法包括：

根据易疲劳热点周边的应力梯度，根据网格模型提取距离热点3t/2和t/2处的表面应力，t为风向周期，根据线性关系外插得到易疲劳热点处的应力，提取不同频率的具有单位波幅的规则波产生的各项载荷的实部τc和虚部τs，得到无人机机体结构在对应频率下的合成应力τ＝τc+iτs，选取易疲劳热点附近的四个有限元网格单元，计算得到的各单元应力τkx，τky，k＝1,2,3,4，为有限元网格标号，热点处的应力τhotx、τhoty为：

τhot的单位为n/mm²。

(6)计算易疲劳热点的应力与疲劳累积损伤度的关系p；

其中，n为设定的常数；

(7)计算对易疲劳热点施加频率为ω的风向u时，易疲劳热点处的应力τhot；

(8)计算对易疲劳热点施加频率为v的风向u时，易疲劳热点处的总疲劳累积损伤度bu；

本发明的特点如下：

在执行无人机机体结构易疲劳热点的结构优化设计过程中，只需计算易疲劳热点在某一计算风向的某个计算频率下的结构热点应力，即可计算得到该优化方案下的典型节点结构的疲劳累积损伤度，而不需要进行全风向全频率的结构疲劳谱分析计算。

该简化计算方法的计算过程简单，可靠性高。经过大量数据的分析研究发现，该简化计算方法计算得到的易疲劳热点的疲劳累积损伤度与采用全风向全频率的谱分析计算方法的误差基本在1％范围内，最高误差可控制在5％左右。该优化方法可以极大减少计算工作量，提高设计效率。设有n2个风向和m2个计算频率，那么该方法完成疲劳热点的优化设计工作的时间最少可为原来全风向全频率的谱分析计算时间的此外，通过合理地取定各类系数值，可保证该简化计算方法的计算结果与采用全风向全频率的谱分析方法计算结果的误差在可接受的范围，进而能充分保障无人机结构的疲劳强度安全。对于任何风向及在某一风向的任意频率下，该简化计算方法均成立；可确保该方法的安全性更好。