CN113673023A - 一种飞机复杂块体结构疲劳载荷谱确定方法 - Google Patents
一种飞机复杂块体结构疲劳载荷谱确定方法 Download PDFInfo
- Publication number
- CN113673023A CN113673023A CN202010412386.2A CN202010412386A CN113673023A CN 113673023 A CN113673023 A CN 113673023A CN 202010412386 A CN202010412386 A CN 202010412386A CN 113673023 A CN113673023 A CN 113673023A
- Authority
- CN
- China
- Prior art keywords
- block structure
- complex block
- airplane
- finite element
- element model
- 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
- 238000001228 spectrum Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000013461 design Methods 0.000 claims abstract description 18
- 239000000463 material Substances 0.000 claims description 3
- 206010016256 fatigue Diseases 0.000 description 35
- 238000004458 analytical method Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012407 engineering method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000013585 weight reducing agent Substances 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)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Toys (AREA)
Abstract
一种飞机复杂块体结构疲劳载荷谱确定方法,在复杂块体结构三维实体细节有限元模型的表面,用位于上述棱边上的三维实体单元节点创建多个连续的杆单元;将复杂块体结构细节有限元模型嵌入飞机总体有限元模型,分别向飞机总体有限元模型施加飞机各使用状态下的设计载荷,计算获得飞机各使用状态下的复杂块体结构表面棱边上所有杆单元应力,形成该复杂块体结构的疲劳载荷谱。
Description
技术领域
本发明属于飞机设计中的疲劳分析相关技术,进一步是一种飞机复杂块体结构疲劳危险部位及其疲劳载荷谱的确定。
背景技术
在飞机结构设计中,为满足承载要求,需要采用一些复杂块体接头结构,其结构特点:1)通过多个耳片与周围结构进行连接,除耳片外,其他部位都为非典型结构,无法通过工程方法获取细节应力,需要建立细节有限元三维实体模型进行应力分析;2)疲劳危险部位较多,包括:台阶、开口、过渡圆角等,,分析工作量大;3)对于其中任意一个疲劳危险部位,在不同飞机使用情况下,例如:地面滑跑情况、爬升情况、巡航情况、着陆撞击情况等,有限元三维实体单元最大主应力方向不一致。不便于组成飞机在整个飞行过程中,疲劳危险部位的载荷时间历程,也称作载荷谱。例如:机身与机翼连接接头结构、机身与起落架连接接头结构都是一种复杂块体结构。
众所周知,在不断变化的载荷作用下,结构会发生疲劳失效破坏。在飞机设计过程中,为确定结构疲劳寿命,需要对所有疲劳危险部位、疲劳薄弱部位,主要为台阶、开口、过渡圆角等,进行疲劳分析。而疲劳分析核心工作为:建立疲劳危险部位载荷谱,也就是对飞机在整个飞行过程中,所受的载荷时间历程的确定。
对于复杂块体结构,耳片疲劳分析方法较为成熟。而对于耳片除外的其他非典型结构,任意一个疲劳危险部位,在不同使用情况下,有限元三维实体单元应力最大主应力方向不一致,无法通过直接提取各个使用情况下疲劳危险部位有限元三维实体单元最大主应力进行组谱分析。现用解决方法为:1)首先提取所有使用情况下,该疲劳危险部位有限元三维实体单元最大主应力,比较确定最大值,记录对应应力方向(定义为:参考方向);2)将所有使用情况下疲劳危险部位有限元三维实体单元最大主应力向参考方向进行投影,获得在参考方向上的应力值;3)将各个使用情况下投影后的单元最大主应力进行组谱分析,形成该危险部位载荷谱。
现有解决方法存在问题:将所有使用情况下疲劳危险部位有限元三维实体单元应力向参考方向进行投影,需要建立大量局部坐标,涉及大量坐标变换工作,工作量大,不便于进行批量编程处理。该方法不适合飞机设计过程中,多轮次结构优化减重分析工作。
发明内容
本申请的目的在于提供一种飞机复杂块体结构疲劳疲劳载荷谱的确定方法,便于进行批量编程处理,以满足型号研制时该类结构疲劳强度评估的需要。
一种飞机复杂块体结构疲劳载荷谱确定方法,已知飞机的设计数模,该设计数模包含所述复杂块体结构的设计数模、飞机总体有限元模型和飞机各使用状态下的设计载荷,所述的复杂块体结构表面含有许多棱边,其特征在于包含以下内容:1)依据飞机复杂块体结构的设计数模建立复杂块体结构三维实体细节有限元模型;2)在复杂块体结构三维实体细节有限元模型的表面,用位于上述棱边上的三维实体单元节点创建多个连续的杆单元;3)将2)中复杂块体结构细节有限元模型嵌入飞机总体有限元模型,分别向飞机总体有限元模型施加飞机各使用状态下的设计载荷,计算获得飞机各使用状态下的复杂块体结构表面棱边上所有杆单元应力;4)重复步骤2)和步骤3),计算获得复杂块体结构的每一个棱边上所有杆单元应力,形成该复杂块体结构的疲劳载荷谱。
本发明的优点是:1)在复杂块体结构表面棱边借用已有三维实体单元节点增加杆单元,由于杆单元仅承受轴力,在飞机各使用状态下,应力方向一致,便于快速形成应力谱。2)引入杆单元方法简单操作方便,可借用已有三维实体单元节点进行设置。3)杆单元属性设置时,材料与复杂块体结构相同,设置截面面积不大于0.001mm2,杆单元的引入,不改变复杂块体结构刚度,不影响真实复杂块体结构传载路径。
以下结合实施例附图对本申请做进一步详细描述。
附图说明
图1是一个飞机复杂块体结构示意图。
图2是基于有限元模型在复杂块体结构的棱边引入杆单元示意图。
图3含杆单元的复杂块体结构细节有限元模型与飞机总体有限元模型示意图。
图4是复杂块体结构耳片根部棱边区疲劳载荷示意图。
图中编号说明:1复杂块体结构、2棱边、3复杂块体结构细节有限元模型、4节点、5杆单元、6飞机总体有限元模型、7复杂块体结构耳片根部弧形区棱边疲劳载荷曲线。
具体实施方式
参见附图,图1提供一个飞机复杂块体结构示意,该复杂块体结构1是飞机机身与起落架的连接接头,如图所示,在复杂块体结构1上有多个需要进行疲劳载荷分析的棱边2,疲劳分析中,复杂块体结构疲劳危险部位全部位于其表面棱边处,本实施例以复杂块体结构耳片根部弧形区棱边2为例依据本申请提供的方法进行疲劳载荷谱的确定。
首先,依据飞机复杂块体结构的设计数模建立该复杂块体结构2的三维实体细节有限元模型3,然后在复杂块体结构三维实体细节有限元模型3的表面,用位于上述棱边2上的三维实体单元节点4创建多个连续的杆单元5,两个节点4之间为一个杆单元5,设定该杆单元5的材料与复杂块体结构2的材料相同,且杆单元的截面面积不大于0.001mm2,这样做的有益效果在于,疲劳分析时,该杆单元不改变复杂块体的结构刚度,不影响复杂块体结构1的真实传载路径;再将上述复杂块体结构细节有限元模型3嵌入飞机总体有限元模型6中,如图3所示,分别向飞机总体有限元模型施加飞机各使用状态下的设计载荷,计算获得飞机各使用状态下的复杂块体结构耳片根本弧形区表面棱边2上所有杆单元5的应力,组建应力谱,形成以时间历程为横轴,以应力值为纵轴的复杂块体结构耳片根部弧形区棱边疲劳载荷曲线7,如图4所示;同理,重复上述步骤,计算获得该复杂块体结构的每一个棱边上所有杆单元应力,形成该复杂块体结构的疲劳载荷谱。
Claims (3)
1.一种飞机复杂块体结构疲劳载荷谱确定方法,已知飞机的设计数模,该设计数模包含所述复杂块体结构的设计数模、飞机总体有限元模型和飞机各使用状态下的设计载荷,所述的复杂块体结构表面含有许多棱边,其特征在于包含以下内容:1)依据飞机复杂块体结构的设计数模建立复杂块体结构三维实体细节有限元模型;2)在复杂块体结构三维实体细节有限元模型的表面,用位于上述棱边上的三维实体单元节点创建多个连续的杆单元;3)将2)中复杂块体结构细节有限元模型嵌入飞机总体有限元模型,分别向飞机总体有限元模型施加飞机各使用状态下的设计载荷,计算获得飞机各使用状态下的复杂块体结构表面棱边上所有杆单元应力;4)重复步骤2)和步骤3),计算获得复杂块体结构的每一个棱边上所有杆单元应力,形成该复杂块体结构的疲劳载荷谱。
2.如权利要求1所述的飞机复杂块体结构疲劳载荷谱确定方法,其特征在于,在步骤2)中,上述杆单元的材料与复杂块体结构相同。
3.如权利要求1或2所述的飞机复杂块体结构疲劳载荷谱确定方法,其特征在于,在步骤2)中,设置杆单元的截面面积不大于0.001mm2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010412386.2A CN113673023B (zh) | 2020-05-15 | 2020-05-15 | 一种飞机复杂块体结构疲劳载荷谱确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010412386.2A CN113673023B (zh) | 2020-05-15 | 2020-05-15 | 一种飞机复杂块体结构疲劳载荷谱确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113673023A true CN113673023A (zh) | 2021-11-19 |
CN113673023B CN113673023B (zh) | 2024-05-24 |
Family
ID=78537643
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010412386.2A Active CN113673023B (zh) | 2020-05-15 | 2020-05-15 | 一种飞机复杂块体结构疲劳载荷谱确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113673023B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114329768A (zh) * | 2021-12-06 | 2022-04-12 | 中航飞机起落架有限责任公司 | 起落架疲劳应力计算方法、系统、设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573172A (zh) * | 2014-11-19 | 2015-04-29 | 新疆金风科技股份有限公司 | 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置 |
CN109710982A (zh) * | 2018-12-04 | 2019-05-03 | 中国特种飞行器研究所 | 基于有限元模拟分析的部件级模拟试验件设计验证方法 |
CN109726436A (zh) * | 2018-12-04 | 2019-05-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种机翼盒段刚度确定方法 |
-
2020
- 2020-05-15 CN CN202010412386.2A patent/CN113673023B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104573172A (zh) * | 2014-11-19 | 2015-04-29 | 新疆金风科技股份有限公司 | 风力发电机组中的结构件的疲劳分析方法和疲劳分析装置 |
CN109710982A (zh) * | 2018-12-04 | 2019-05-03 | 中国特种飞行器研究所 | 基于有限元模拟分析的部件级模拟试验件设计验证方法 |
CN109726436A (zh) * | 2018-12-04 | 2019-05-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种机翼盒段刚度确定方法 |
Non-Patent Citations (5)
Title |
---|
SLOBODANKA BOLJANOVIĆ等: ""Computational modeling of aircraft lugs failure under fatigue loading"", 《INTERNATIONAL JOURNAL OF FATIGUE》, 22 May 2018 (2018-05-22), pages 252 - 261, XP085528765, DOI: 10.1016/j.ijfatigue.2018.05.022 * |
Y. GARBATOV等: ""Fatigue assessment of marine structures"", 《MARINE TECHNOLOGY AND ENGINEERING – GUEDES SOARES ET AL》, 31 January 2011 (2011-01-31), pages 864 - 888 * |
万一品等: "《装载机工作装置有限元分析与疲劳强度评估》", 《起重运输机械》, 20 September 2016 (2016-09-20), pages 30 - 34 * |
李毅等: ""操纵面非线性气动弹性响应引起的结构疲劳损伤分析方法研究"", 《振动与冲击》, vol. 34, no. 11, 15 June 2015 (2015-06-15), pages 41 - 45 * |
龚思楚等: "《某型飞机机翼壁板对接区细节分析及 寿命估算》", 《TRAINER》, 15 March 2020 (2020-03-15), pages 61 - 64 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114329768A (zh) * | 2021-12-06 | 2022-04-12 | 中航飞机起落架有限责任公司 | 起落架疲劳应力计算方法、系统、设备及存储介质 |
CN114329768B (zh) * | 2021-12-06 | 2024-05-07 | 中航飞机起落架有限责任公司 | 起落架疲劳应力计算方法、系统、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113673023B (zh) | 2024-05-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106156449B (zh) | 一种复合材料机翼壁板优化设计方法 | |
CN109918694B (zh) | 一种基于离散化载荷模型的飞机载荷计算方法 | |
CN101604353B (zh) | 翼组件的设计方法 | |
CN111274648B (zh) | 一种民用飞机前缘襟翼的分布式飞行载荷设计方法 | |
Zakuan et al. | Modelling and structural analysis of three-dimensional wing | |
CN110929444A (zh) | 一种滑橇式起落架仿真设计方法 | |
CN113673023A (zh) | 一种飞机复杂块体结构疲劳载荷谱确定方法 | |
CN107273638B (zh) | 一种基于平尾载荷的机身与尾翼连接铰点载荷分配方法 | |
CN109726435B (zh) | 一种接头钉群载荷计算的方法 | |
Schuhmacher et al. | Multidisciplinary design optimization of a regional aircraft wing box | |
CN105528478A (zh) | 一种飞机翼面结构传力特点的图示分析方法 | |
CN103366094B (zh) | 一种基于力载荷和力边界条件的确定复合材料结构强度的方法 | |
US11479336B2 (en) | Diagonal pressure deck | |
CN108052711A (zh) | 一种并联组合发动机安装结构的热力耦合分析方法 | |
CN108090260B (zh) | 一种接头约束载荷的分析方法 | |
Meink | Composite grid vs. composite sandwich: a comparison based on payload shroud requirements | |
CN109388868B (zh) | 一种考虑钉孔接触及钉弯曲的钉载计算方法 | |
CN111008498A (zh) | 超静定吊挂与机翼连接单元结构尺寸的分析方法 | |
CN109827723A (zh) | 固定翼飞机全机强度及模态试验柔性支架系统 | |
Ahaus et al. | Assessment of CFD/CSD analytical tools for improved rotor loads | |
CN110795871A (zh) | 一种水陆两栖飞机船底框强度计算和试验方法 | |
Zhou et al. | Characteristics analysis and optimization of flying-wing vehicle structure | |
CN113449375B (zh) | 一种复合材料桨叶疲劳寿命的半解析计算方法 | |
CN113483989A (zh) | 一种飞机发动机模型设计方法、制造方法及飞机发动机模型 | |
Sobieczky | Aircraft surface generation |
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 |