CN113190923A - 一种筋条整体壁板损伤容限优化方法 - Google Patents
一种筋条整体壁板损伤容限优化方法 Download PDFInfo
- Publication number
- CN113190923A CN113190923A CN202110628427.6A CN202110628427A CN113190923A CN 113190923 A CN113190923 A CN 113190923A CN 202110628427 A CN202110628427 A CN 202110628427A CN 113190923 A CN113190923 A CN 113190923A
- Authority
- CN
- China
- Prior art keywords
- rib
- crack
- parameters
- integral
- damage tolerance
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000005457 optimization Methods 0.000 title claims description 23
- 239000011159 matrix material Substances 0.000 claims description 15
- 230000002787 reinforcement Effects 0.000 claims description 10
- 238000005452 bending Methods 0.000 claims description 7
- 230000000149 penetrating effect Effects 0.000 claims description 7
- 230000035945 sensitivity Effects 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 4
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 239000011324 bead Substances 0.000 claims description 3
- 230000000644 propagated effect Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 108010063499 Sigma Factor Proteins 0.000 claims description 2
- 230000003068 static effect Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004519 manufacturing process Methods 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
- G06F2113/00—Details relating to the application field
- G06F2113/28—Fuselage, exterior or interior
-
- 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)
- 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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明属于飞机结构损伤容限设计领域,公开了一种筋条整体壁板损伤容限优化方法,包括整体壁板加筋比的确定、筋条敏感参数识别、结构尺寸优化矩阵的建立、含裂纹整体壁板有限元模型建立、应力强度因子k的计算、裂纹扩展寿命的计算,以及通过不同尺寸下裂纹扩展寿命的对比确定最优结构参数。
Description
技术领域
本发明属于飞机结构损伤容限设计领域,特别是涉及到一种筋条整体壁板损伤容限优化方法。
背景技术
整体壁板由于具有疲劳性能好、制造简单、结构重量轻的优点,已经在多个中大型飞机结构中得到应用。但是,与传统的铆接壁板相比,整体壁板缺乏天然止裂原件——长桁与蒙皮连接的铆钉,使得损伤容限性能较差,因此如何通过筋条的优化设计来提高整体壁板的损伤容限特性,是整体壁板用于飞机结构的重要设计点。在过去几年里,整体壁板在结构设计时,一般关注静强度性能,若满足要求,仅仅对壁板是否满足损伤要求进行被动式校核,而没有进行主动的损伤容限优化,这就使得壁板的裂纹扩展性能不能达到最优,没有最大限度的提高损伤容限性能。
发明内容
本发明的目的是提供一种筋条整体壁板损伤容限优化方法,解决在整体壁板优化设计上流程和方法缺乏的问题,从而最大程度上挖掘整体壁板的损伤容限性能,保证飞机使用安全。该方法也可以借鉴用于其他筋条形式的整体壁板损伤容限优化设计。
为达到上述目的,本发明采用如下技术方案予以实现。
一种筋条整体壁板损伤容限优化方法,包括以下步骤:
S1,确定筋条整体壁板的多组待优化参数,组成优化矩阵;
S2,选择所述优化矩阵中的一组筋条参数;
S3,建立至少含有两根筋条的整体壁板有限元模型,并在中心蒙皮处假设初始长度为a的穿透裂纹;
S4,在所述有限元模型中扩展所述穿透裂纹,并计算所述一组筋条参数下的裂纹扩展寿命;
S5,比较不同筋条参数下的裂纹扩展寿命,选取裂纹扩展寿命最大的筋条参数为最优结果,在所述最有结果下整体壁板损伤容限性能最好。
本发明技术方案的特点和进一步的改进为:
(1)S1,具体为:
S11,给定蒙皮厚度,桁距,根据加筋比、蒙皮厚度以及桁距计算筋条截面积;
S12,给定整体壁板的弯曲刚度、重量指标,根据弯曲刚度、重量指标确定筋条尺寸的敏感性参数;所述敏感性参数为:筋条厚度、长桁高度和缘条宽度;
S13,根据筋条截面积,在筋条厚度、长桁高度和缘条宽度的合理范围内确定出多组筋条厚度、长桁高度和缘条宽度的值;将蒙皮厚度、桁距、筋条厚度、长桁高度、缘条宽度组成优化矩阵。
(2)S4,具体为:
S41,获取应力强度因子随裂纹长度的变化曲线,且裂纹需扩展至穿过筋条;
S42,计算裂纹扩展寿命,获得该筋条参数下的裂纹扩展曲线;
S43,按照优化矩阵,选择下一组筋条参数,重复S3至S5,获得不同筋条参数下的裂纹扩展曲线及裂纹扩展寿命.
(3)S41,具体为:
(4)S42具体为:采用walker公式计算裂纹扩展寿命N:
(5)S11中,所述加筋比的确定方法为:
(6)S4中,在所述有限元模型中扩展所述穿透裂纹,采用自适应网格划分进行裂纹扩展。
(7)所述筋条整体壁板的剖面形式为“T”、“L”或“I”型。
本发明结合传统优化设计方法与断裂力学,建立了适用于整体加筋壁板损伤容限设计方法和流程,采用本方法确定出的筋条尺寸,能够在满足静力和重量指标条件下,最大限度的延缓裂纹扩展速度,提高了整体壁板的损伤容限性能,从而更能够保证整体壁板在飞机结构中的安全使用。本发明使用的理论依据正确,分析步骤清晰、简单,计算过程便于计算机自动化迭代计算。本发明的提出解决了飞机结构中采用整体壁板时是否能够满足损伤容限要求的难题。
附图说明
图1是一种筋条整体壁板损伤容限优化方的流程示意图;
图2是“Z”型筋条的整体壁板剖面示意图;
图3是整体壁板含裂纹的总体有限元模型;
图4是裂纹尖端有限元网格图;
其中,1-蒙皮,2-Z型筋条,t1-蒙皮厚度,B-长桁间距,H-筋条高度,W-筋条宽度,t2-筋条厚度。
具体实施方式
本发明实施例提供一种筋条整体壁板损伤容限优化方法,如图1所示,包括以下步骤:
S1,确定筋条整体壁板的多组待优化参数,组成优化矩阵;
S2,选择所述优化矩阵中的一组筋条参数;
S3,建立至少含有两根筋条的整体壁板有限元模型,并在中心蒙皮处假设初始长度为a的穿透裂纹;
S4,在所述有限元模型中扩展所述穿透裂纹,并计算所述一组筋条参数下的裂纹扩展寿命;
S5,比较不同筋条参数下的裂纹扩展寿命,选取裂纹扩展寿命最大的筋条参数为最优结果,在所述最有结果下整体壁板损伤容限性能最好。
进一步的:
(1)S1,具体为:
S11,给定蒙皮厚度,桁距,根据加筋比、蒙皮厚度以及桁距计算筋条截面积;
S12,给定整体壁板的弯曲刚度、重量指标,根据弯曲刚度、重量指标确定筋条尺寸的敏感性参数;所述敏感性参数为:筋条厚度、长桁高度和缘条宽度;
S13,根据筋条截面积,在筋条厚度、长桁高度和缘条宽度的合理范围内确定出多组筋条厚度、长桁高度和缘条宽度的值;将蒙皮厚度、桁距、筋条厚度、长桁高度、缘条宽度组成优化矩阵。
(2)S4,具体为:
S41,获取应力强度因子随裂纹长度的变化曲线,且裂纹需扩展至穿过筋条;
S42,计算裂纹扩展寿命,获得该筋条参数下的裂纹扩展曲线;
S43,按照优化矩阵,选择下一组筋条参数,重复S3至S5,获得不同筋条参数下的裂纹扩展曲线及裂纹扩展寿命.
(3)S41,具体为:
(4)S42具体为:采用walker公式计算裂纹扩展寿命N:
(5)S11中,所述加筋比的确定方法为:
(6)S4中,在所述有限元模型中扩展所述穿透裂纹,采用自适应网格划分进行裂纹扩展。
(7)所述筋条整体壁板的剖面形式为“T”、“L”或“I”型。
下面以某一具体实例对本发明做进一步详细说明。
一种“Z”型筋条整体壁板损伤容限优化设计的方法,包括以下步骤:
步骤二、蒙皮厚度t1=5mm,筋条间距B=160mm,根据加筋比R0=0.3计算筋条截面积Astr=325mm2;
步骤三、在给定的弯曲刚度、重量指标下,确定的敏感参数为:筋条厚度t2、、高度H和宽度w,如图2所示。
步骤四:确定筋条整体壁板的五组待优化参数,组成优化矩阵,如表1所示;
表1优化矩阵
步骤五、选择第1种筋条参数,建立至少含有2根筋条的整体壁板有限元模型,并在中心蒙皮处假设长度为a=20mm的穿透裂纹,如图3和图4所示;
步骤六、结合网格自适应软件,获取应力强度因子随裂纹长度K的变化曲线,裂纹需扩展至穿过筋条;
步骤七、利用walker计算公式,计算裂纹扩展寿命N,获得该筋条尺寸下的裂纹扩展曲线,walker公式如下:
步骤八、按照优化矩阵,依次改变筋条尺寸,重复步骤六至步骤八,获得不同筋条尺寸下的裂纹扩展寿命,如表2所示。
表2各构型下优化结果
步骤九、比较表2中不同筋条尺寸下的裂纹扩展寿命,可以看出构型3的裂纹扩展寿命最长,因此构型3的尺寸为最优结果。
本发明结合传统优化设计方法与断裂力学,建立了适用于整体加筋壁板损伤容限设计方法和流程,采用本方法确定出的筋条尺寸,能够在满足静力和重量指标条件下,最大限度的延缓裂纹扩展速度,提高了整体壁板的损伤容限性能,从而更能够保证整体壁板在飞机结构中的安全使用。本发明使用的理论依据正确,分析步骤清晰、简单,计算过程便于计算机自动化迭代计算。本发明的提出解决了飞机结构中采用整体壁板时是否能够满足损伤容限要求的难题。
Claims (8)
1.一种筋条整体壁板损伤容限优化方法,其特征在于,包括以下步骤:
S1,确定筋条整体壁板的多组待优化参数,组成优化矩阵;
S2,选择所述优化矩阵中的一组筋条参数;
S3,建立至少含有两根筋条的整体壁板有限元模型,并在中心蒙皮处假设初始长度为a的穿透裂纹;
S4,在所述有限元模型中扩展所述穿透裂纹,并计算所述一组筋条参数下的裂纹扩展寿命;
S5,比较不同筋条参数下的裂纹扩展寿命,选取裂纹扩展寿命最大的筋条参数为最优结果,在所述最有结果下整体壁板损伤容限性能最好。
2.根据权利要求1所述的一种筋条整体壁板损伤容限优化方法,其特征在于,S1,具体为:
S11,给定蒙皮厚度,桁距,根据加筋比、蒙皮厚度以及桁距计算筋条截面积;
S12,给定整体壁板的弯曲刚度、重量指标,根据弯曲刚度、重量指标确定筋条尺寸的敏感性参数;所述敏感性参数为:筋条厚度、长桁高度和缘条宽度;
S13,根据筋条截面积,在筋条厚度、长桁高度和缘条宽度的合理范围内确定出多组筋条厚度、长桁高度和缘条宽度的值;将蒙皮厚度、桁距、筋条厚度、长桁高度、缘条宽度组成优化矩阵。
3.根据权利要求1所述的一种筋条整体壁板损伤容限优化方法,其特征在于,S4,具体为:
S41,获取应力强度因子随裂纹长度的变化曲线,且裂纹需扩展至穿过筋条;
S42,计算裂纹扩展寿命,获得该筋条参数下的裂纹扩展曲线;
S43,按照优化矩阵,选择下一组筋条参数,重复S3至S5,获得不同筋条参数下的裂纹扩展曲线及裂纹扩展寿命。
7.根据权利要求1所述的一种筋条整体壁板损伤容限优化方法,其特征在于,S4中,在所述有限元模型中扩展所述穿透裂纹,采用自适应网格划分进行裂纹扩展。
8.根据权利要求1所述的一种筋条整体壁板损伤容限优化方法,其特征在于,所述筋条整体壁板的剖面形式为“T”、“L”或“I”型。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2020116129696 | 2020-12-29 | ||
CN202011612969.6A CN112733262A (zh) | 2020-12-29 | 2020-12-29 | 一种筋条整体壁板损伤容限优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113190923A true CN113190923A (zh) | 2021-07-30 |
CN113190923B CN113190923B (zh) | 2024-04-09 |
Family
ID=75611804
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011612969.6A Withdrawn CN112733262A (zh) | 2020-12-29 | 2020-12-29 | 一种筋条整体壁板损伤容限优化方法 |
CN202110628427.6A Active CN113190923B (zh) | 2020-12-29 | 2021-06-04 | 一种筋条整体壁板损伤容限优化方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011612969.6A Withdrawn CN112733262A (zh) | 2020-12-29 | 2020-12-29 | 一种筋条整体壁板损伤容限优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN112733262A (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112733262A (zh) * | 2020-12-29 | 2021-04-30 | 中国航空工业集团公司西安飞机设计研究所 | 一种筋条整体壁板损伤容限优化方法 |
CN113343529B (zh) * | 2021-06-11 | 2022-07-12 | 清华大学 | 一种整体壁板结构损伤断裂的全局控制方法和装置 |
CN113443169B (zh) * | 2021-06-11 | 2022-08-02 | 航空工业第一飞机设计研究院 | 一种整体壁板结构损伤断裂主动控制的试验方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106354934A (zh) * | 2016-08-30 | 2017-01-25 | 中航沈飞民用飞机有限责任公司 | 一种基于有限元分析求解的加筋壁板损伤容限分析方法 |
CA3037323A1 (en) * | 2017-10-31 | 2019-04-30 | China University Of Mining And Technology | Reliability robust design method for multiple failure modes of ultra-deep well hoisting container |
CN111046610A (zh) * | 2019-12-26 | 2020-04-21 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机整体翼梁无量纲应力强度因子的计算方法 |
CN111159944A (zh) * | 2019-12-26 | 2020-05-15 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机结构疲劳裂纹扩展轨迹的确定方法 |
US20200264961A1 (en) * | 2019-02-15 | 2020-08-20 | International Business Machines Corporation | Injection of simulated hardware failure(s) in a file system for establishing file system tolerance-to-storage-failure(s) |
CN112733262A (zh) * | 2020-12-29 | 2021-04-30 | 中国航空工业集团公司西安飞机设计研究所 | 一种筋条整体壁板损伤容限优化方法 |
-
2020
- 2020-12-29 CN CN202011612969.6A patent/CN112733262A/zh not_active Withdrawn
-
2021
- 2021-06-04 CN CN202110628427.6A patent/CN113190923B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106354934A (zh) * | 2016-08-30 | 2017-01-25 | 中航沈飞民用飞机有限责任公司 | 一种基于有限元分析求解的加筋壁板损伤容限分析方法 |
CA3037323A1 (en) * | 2017-10-31 | 2019-04-30 | China University Of Mining And Technology | Reliability robust design method for multiple failure modes of ultra-deep well hoisting container |
US20200264961A1 (en) * | 2019-02-15 | 2020-08-20 | International Business Machines Corporation | Injection of simulated hardware failure(s) in a file system for establishing file system tolerance-to-storage-failure(s) |
CN111046610A (zh) * | 2019-12-26 | 2020-04-21 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机整体翼梁无量纲应力强度因子的计算方法 |
CN111159944A (zh) * | 2019-12-26 | 2020-05-15 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机结构疲劳裂纹扩展轨迹的确定方法 |
CN112733262A (zh) * | 2020-12-29 | 2021-04-30 | 中国航空工业集团公司西安飞机设计研究所 | 一种筋条整体壁板损伤容限优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113190923B (zh) | 2024-04-09 |
CN112733262A (zh) | 2021-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113190923B (zh) | 一种筋条整体壁板损伤容限优化方法 | |
CN105184390A (zh) | 一种壁板结构静强度、刚度、稳定性的综合优化方法 | |
CN107506529B (zh) | 一种复合材料加筋壁板轴压稳定性计算方法 | |
CN112528540B (zh) | 一种基于力热耦合计算的结构优化方法 | |
CN111046610B (zh) | 一种飞机整体翼梁无量纲应力强度因子的计算方法 | |
Es-Haghi et al. | Evaluation of a novel Asymmetric Genetic Algorithm to optimize the structural design of 3D regular and irregular steel frames | |
CN112784359A (zh) | 薄壁梁约束扭转极限承载力迭代计算方法 | |
CN111914351B (zh) | 一种机身结构加筋壁板总体稳定性计算方法 | |
CN110489914B (zh) | 一种基于应力损伤等效的耐久性计算方法 | |
CN109614713B (zh) | 基于修正刚度方法的复合材料多钉连接钉载分配预测方法 | |
CN105184060A (zh) | 一种基于空间配面与改进傅里叶级数的层合结构振动分析方法 | |
CN113094793B (zh) | 一种减隔震结构选波及设计方法 | |
CN112733395B (zh) | 型钢混凝土组合构件从损伤到失效的综合评价分析方法 | |
CN111274670B (zh) | 一种钛合金型材组合壁板轴压承载能力试验设计方法 | |
CN105501462B (zh) | 一种机翼结构设计方法 | |
CN216380374U (zh) | 一种散热性能好的铝蜂窝板 | |
CN112711804B (zh) | 一种增升装置壁板隔波强度分析方法 | |
CN112733260B (zh) | 一种运输类飞机机身长桁通过孔的dfr疲劳计算方法 | |
CN109583090A (zh) | 一种复合材料结构的铺层优化方法 | |
US20170197170A1 (en) | Defect tolerant honeycomb structures | |
CN111027125A (zh) | 一种蜂窝等效单元及基于该蜂窝等效单元的参数计算方法 | |
CN113343529B (zh) | 一种整体壁板结构损伤断裂的全局控制方法和装置 | |
CN112800538B (zh) | 一种任意薄壁梁结构的层合化分析方法 | |
Dow et al. | Structural Evaluation of an Extruded Magnesium-Alloy T-Stiffened Panel | |
CN112699471B (zh) | 一种机身壁板轴压载荷下蒙皮有效宽度计算方法及装置 |
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 |