CN104034576B - 一种飞机结构裂纹扩展试验用载荷谱简化方法 - Google Patents
一种飞机结构裂纹扩展试验用载荷谱简化方法 Download PDFInfo
- Publication number
- CN104034576B CN104034576B CN201410266931.6A CN201410266931A CN104034576B CN 104034576 B CN104034576 B CN 104034576B CN 201410266931 A CN201410266931 A CN 201410266931A CN 104034576 B CN104034576 B CN 104034576B
- Authority
- CN
- China
- Prior art keywords
- crack
- test
- loading spectrum
- intensity factor
- stress
- 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
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种飞机结构裂纹扩展试验用载荷谱简化方法,其包括下列步骤:(1)确定裂纹扩展试验件的结构形式及含裂纹部分的材料,并确定材料的疲劳裂纹扩展门槛值ΔKth;(2)再确定结构裂纹尖端的应力强度因子;(3)对裂纹扩展试验中的临界裂纹长度进行估算,并进行分段处理;(4)根据分段情况,对载荷谱中低于疲劳裂纹扩展门槛值的应力循环进行删除,形成各段裂纹长度对应的载荷谱,最终应用于裂纹扩展试验。本发明确定含裂纹结构的裂尖应力强度因了值,结合材料的疲劳裂纹扩展门槛值ΔKth;对裂纹扩展试验用载荷谱进行分段简化,形成了一系列简化的载荷谱,可直接用于飞机结构裂纹扩展试验,方法简单易行,且能有效缩短试验周期。
Description
技术领域
本发明属于飞机载荷谱加载技术,涉及一种飞机结构裂纹扩展试验用载荷谱简化方法。
背景技术
国内航空界蓬勃发展,越来越多的民/军机需要进行疲劳/裂纹扩展试验,而试验周期一般都比较长,严重制约了飞机的研制时间,同时对人力、财力也是一种浪费。
一直以来,通过简化载荷谱来缩短试验周期的思路被延续着,在此基础上,形成了低载截除、高载截取技术、简化为常幅谱等。其中低载截除技术的思想就是对于载荷谱中对裂纹扩展没有做出贡献的载荷级进行截除,从而减小一块载荷谱中的载荷级数以达到减少试验时间的目的。其中一项截除标准就是裂纹扩展门槛值ΔKth,其与裂纹尖端应力强度因子相关。
正是基于裂纹尖端应力强度因子在裂纹扩展寿命计算中的重要地位,国内外对其的研究已经很深入,并形成了成熟的技术,如有限元素法(用的比较多)、工程法及试验方法。这些方法可以直接用来计算裂尖应力强度因子,从而用于裂纹扩展寿命的计算。
对于一般航空金属材料,其裂纹扩展寿命曲线分为三个部分:裂纹萌生阶段、裂纹扩展阶段、快速扩展阶段,这是已被航空界所接受的理论。然而在现有技术航空疲劳/裂纹扩展试验中,由于航空金属材料裂纹扩展的复杂性,使得疲劳试验周期一般都比较长,严重制约了飞机的研制时间,增加了研制成本。
发明内容
本发明的目的:通过对试验用载荷谱进行分段处理,提供一种飞机结构裂纹扩展试验用载荷谱简化方法,有效缩短了裂纹扩展试验周期。
本发明的技术方案:一种飞机结构裂纹扩展试验用载荷谱简化方法,先确定裂纹扩展试验件的结构形式及含裂纹部分的材料,并确定材料的疲劳裂纹扩展门槛值ΔKth;再确定结构裂纹尖端的应力强度因子;然后对裂纹扩展试验中的临界裂纹长度进行估算,并进行分段处理,对裂纹扩展的分段处理时,借鉴同类材料及结构的裂纹扩展曲线;最后根据分段情况,对载荷谱中低于疲劳裂纹扩展门槛值的应力循环进行删除,形成各段裂纹长度对应的载荷谱。
所述的飞机结构裂纹扩展试验用载荷谱简化方法,其具体过程如下:步骤1:确定裂纹扩展试验件的结构形式及含裂纹部分的材料,并针对不同的材料,根据试片试验测得疲劳裂纹扩展门槛值ΔKth;
步骤2:确定结构裂纹尖端的应力强度因子
其中,应力强度因子的通常函数表达式如下:
K=f(β,σ,a)
式中:
σ—加载应力;
a—为半裂纹长度;
β—构型参数,与结构构型参数、裂纹长度有关;
其中,构型参数β根据具体结构查应力强度因子手册求得,从而求得应力强度因子或者采用有限元方法间接求得应力强度因子;
步骤3:对裂纹扩展试验中的临界裂纹长度进行估算;
步骤4:对裂纹长度进行分段处理,裂纹长度的划分与后续试验的精度关系对应,对裂纹扩展的分段处理时,借鉴同类材料及结构的裂纹扩展曲线,裂纹萌生阶段及稳态扩展阶段,多分段;而对裂纹快速扩展阶段,分成一段;
步骤5:对载荷谱进行截除,形成一系列新的载荷谱
根据步骤4中的分段情况,并结合步骤2中的应力强度因子计算公式,计算每段裂纹长度下对应的最大应力强度因子值,采用疲劳裂纹扩展门槛值对载荷谱中的应力值进行截除,最终形成一系列试验用载荷谱。
裂纹扩展试验过程中,根据实际裂纹扩展的长度选取不同的载荷谱进行加载。
裂纹扩展试验中,临界裂纹长度的确定有净截面屈服方法、断裂韧性判据、R曲线判据及COD判据,选择其中的一种或者几种估算得到对具体结构的临界裂纹长度ac。
所述应力强度因子的函数表达式为
分段的方法根据截除载荷的标准而定,先确定截除的载荷大小σ,然后通过公式(5)计算相应的半裂纹长度a,那么对应a时,删除载荷谱谱中所有小于σ的应力值从而简化试验载荷谱,达到缩短试验时间的目的,
其中,ΔKth为疲劳裂纹扩展门槛值,w为板宽。
本发明的有益效果:本发明飞机结构裂纹扩展试验用载荷谱简化方法通过分段处理简化试验用载荷谱,能达到缩短试验周期的效果,其处理过程简单,可操作性强,有效性强并方法简单易行,非常适用于飞机结构(金属)的裂纹扩展试验,具有较大的工程应用价值。
附图说明
图1是本发明一种飞机结构裂纹扩展试验用载荷谱简化方法—典型裂纹扩展曲线的分段处理关系图;
图2是本发明一种飞机结构裂纹扩展试验用载荷谱简化方法的流程图。
具体实施方式
下面通过具体实施方式对本发明作进一步的详细说明:
请参阅图1,其是本发明飞机结构裂纹扩展试验用载荷谱简化方法—典型裂纹扩展曲线的分段处理关系图。本发明飞机结构裂纹扩展试验用载荷谱简化方法的主要特征参数包括裂纹扩展门槛值ΔKth、裂尖应力强度因子K、半裂纹长度a、载荷谱。
在整个裂纹扩展过程中,基于低载截除思想并结合DKth、K及a对载荷谱进行分段截除,由图1可以看出,每个阶段加载的载荷谱是不同的。
请参阅图2,其是本发明飞机结构裂纹扩展试验用载荷谱简化方法的流程图,本发明的具体流程步骤为:
步骤1:确定裂纹扩展试验件的结构形式及含裂纹部分的材料,并确定材料的疲劳裂纹扩展门槛值ΔKth;
具体实施方式中,针对不同的材料,可根据小试片试验测得疲劳裂纹扩展门槛值ΔKth,对于常用材料,可通过查1993年科学出版社出版的《应力强度因子手册》得到,而对于新兴材料,则需要通过小试片试验测得。
步骤2:确定结构裂纹尖端的应力强度因子;
对含裂纹结构裂纹尖端应力强度因子的求解是本文中至关重要的一步,其直接影响了后续试验的精度和误差。应力强度因子的通常函数表达式如公式(1),而典型结构的应力强度因子表达式为公式(2):
K=f(β,σ,a)(1)
式中:
σ—加载应力;
a—为半裂纹长度;
β—与结构构型参数、裂纹长度有关的参数,简称构型参数;
需要指出的是,其中构型参数β的确定是重中之重,通常的方法是根据具体结构查《应力强度因子手册》得到;而对于手册中没有的或者是结构形式比较复杂的,通常会采用有限元方法进行解决,有限元方法往往采用间接法求得应力强度因子,如计算J积分、裂尖开张位移量等。对于具体的问题需要进行具体分析。例如,对于含中心孔孔边开裂的有限宽板而言,其构型因子β的表达式为公式(3);而对于典型机身壁板结构,若其框开裂,构型因子β的表达式为公式(4);而对于复杂结构,如翼身接头而言,若其表面出现裂纹,其应力强度因子则需要通过有限元软件计算得到,如Abaqus中的J积分得到。
β=β1β2β3… (4)
式中:
w—板宽;
β1、β2、β3—为与结构有关的参数,详细可参考《结构耐久性损伤容限设计手册》;
通过步骤2求出的裂尖应力强度因子K值是为了与疲劳裂纹扩展门槛值进行比较,而由于K值是与裂纹长度有关的,故此处的K是动态的。
步骤3:对裂纹扩展试验中的临界裂纹长度进行估算:
对于一般金属结构而言,裂纹扩展试验中,临界裂纹长度的确定有净截面屈服方法、断裂韧性判据、R曲线判据及COD判据等,可以选择其中的一种或者几种估算得到对具体结构的临界裂纹长度ac。
步骤4:对裂纹长度进行分段处理,裂纹长度的划分与后续试验的精度关系比较密切。因此,对裂纹扩展的分段处理时,应该借鉴同类材料及结构的裂纹扩展曲线(a-N曲线);此处采用的原则是:
a)裂纹萌生阶段及稳态扩展阶段,多分段;
b)快速扩展阶段,分成一段;
具体分段数根据实际结构而定,分段越多,裂纹扩展过程中节省的时间越少,但是,载荷谱调用的时间会偏多,应该予以权衡。此处以含中心孔孔边开裂的有限宽板为例进行简要说明。结合公式(2)和公式(3),并带入疲劳裂纹扩展门槛值,可以得到公式(5),可以看出,公式(5)的、板宽w和疲劳裂纹扩展门槛值ΔKth均为常数,则半裂纹长度a与截除的应力值σ成反比,可以得到重要结论:分段越多,截除的小载荷也就越多。
分段的方法可以根据截除载荷的标准而定,如先确定截除的载荷大小,即σ,然后通过公式(5)计算相应的半裂纹长度a,那么对应a时,载荷谱谱中所有小于σ的应力值即可被删除,这样就简化试验载荷谱,从而达到缩短试验时间的目的。
步骤5:对载荷谱进行截除,形成一系列新的载荷谱;
根据步骤4中的分段情况,并结合步骤2中的应力强度因子计算公式,可以计算每段裂纹长度下对应的最大应力强度因子值,采用疲劳裂纹扩展门槛值对载荷谱中的应力值进行截除,具体方法如步骤4中描述,确定相应的半裂纹长度a,通过公式(5)计算截除应力σ,那么载荷谱中小于应力σ的值即可被删除,载荷谱谱中的其他应力大小和位置均不变。最终形成一系列试验用载荷谱,在试验过程中,及时观察结构的裂纹长度,当裂纹长度达到之前计算时用到的长度a时,选用此处对应的载荷谱进行试验。
本发明飞机结构裂纹扩展试验用载荷谱简化方法通过载荷谱分段处理技术提出了一种试验加载方法,能达到缩短试验周期的效果,其处理过程简单,可操作性强,有效性强并方法简单易行,非常适用于飞机结构(金属)的裂纹扩展试验,具有较大的工程应用价值。
Claims (6)
1.一种飞机结构裂纹扩展试验用载荷谱简化方法,其特征在于,先确定裂纹扩展试验件的结构形式及含裂纹部分的材料,并确定材料的疲劳裂纹扩展门槛值ΔKth;再确定结构裂纹尖端的应力强度因子;然后对裂纹扩展试验中的临界裂纹长度进行估算,并进行分段处理,对裂纹扩展的分段处理时,借鉴同类材料及结构的裂纹扩展曲线;最后根据分段情况,对载荷谱中低于疲劳裂纹扩展门槛值的应力循环进行删除,形成各段裂纹长度对应的载荷谱。
2.根据权利要求1所述的飞机结构裂纹扩展试验用载荷谱简化方法,其特征在于,具体过程如下:
步骤1:确定裂纹扩展试验件的结构形式及含裂纹部分的材料,并针对不同的材料,根据试片试验测得疲劳裂纹扩展门槛值ΔKth;
步骤2:确定结构裂纹尖端的应力强度因子
其中,应力强度因子的函数表达式如下:
K=f(β,σ,a)
式中:
σ—加载应力;
a—为半裂纹长度;
β—构型参数,与结构构型参数、裂纹长度有关;
其中,构型参数β根据具体结构查应力强度因子手册求得,从而求得应力强度因子或者采用有限元方法间接求得应力强度因子;
步骤3:对裂纹扩展试验中的临界裂纹长度进行估算;
步骤4:对裂纹长度进行分段处理,裂纹长度的划分与后续试验的精度关系对应,对裂纹扩展的分段处理时,借鉴同类材料及结构的裂纹扩展曲线,裂纹萌生阶段及稳态扩展阶段,多分段;而对裂纹快速扩展阶段,分成一段;
步骤5:对载荷谱进行截除,形成一系列新的载荷谱
根据步骤4中的分段情况,并结合步骤2中的应力强度因子计算公式,计算每段裂纹长度下对应的最大应力强度因子值,采用疲劳裂纹扩展门槛值对载荷谱中的应力值进行截除,最终形成一系列试验用载荷谱。
3.根据权利要求2所述的飞机结构裂纹扩展试验用载荷谱简化方法,其特征在于:裂纹扩展试验过程中,根据实际裂纹扩展的长度选取不同的载荷谱进行加载。
4.根据权利要求3所述的飞机结构裂纹扩展试验用载荷谱简化方法,其特征在于:裂纹扩展试验中,临界裂纹长度的确定有净截面屈服方法、断裂韧性判据、R曲线判据及COD判据,选择其中的一种或者几种估算得到对具体结构的临界裂纹长度ac。
5.根据权利要求4所述的飞机结构裂纹扩展试验用载荷谱简化方法,其特征在于:应力强度因子的函数表达式为
6.根据权利要求4所述的飞机结构裂纹扩展试验用载荷谱简化方法,其特征在于:分段的方法根据截除载荷的标准而定,先确定截除的载荷大小σ,然后通过公式(5)计算相应的半裂纹长度a,那么对应a时,删除载荷谱谱中所有小于σ的应力值从而简化试验载荷谱,达到缩短试验时间的目的,
其中,ΔKth为疲劳裂纹扩展门槛值,w为板宽。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410266931.6A CN104034576B (zh) | 2014-06-16 | 2014-06-16 | 一种飞机结构裂纹扩展试验用载荷谱简化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410266931.6A CN104034576B (zh) | 2014-06-16 | 2014-06-16 | 一种飞机结构裂纹扩展试验用载荷谱简化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104034576A CN104034576A (zh) | 2014-09-10 |
CN104034576B true CN104034576B (zh) | 2016-11-02 |
Family
ID=51465431
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410266931.6A Active CN104034576B (zh) | 2014-06-16 | 2014-06-16 | 一种飞机结构裂纹扩展试验用载荷谱简化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104034576B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104330253B (zh) * | 2014-11-28 | 2017-03-08 | 中航沈飞民用飞机有限责任公司 | 一种加筋壁板损伤容限特性分析方法 |
CN105716934B (zh) * | 2014-12-03 | 2018-04-10 | 中国飞机强度研究所 | 一种结构缺陷尺寸的试验检测方法 |
CN104614227B (zh) * | 2015-02-15 | 2017-04-19 | 东南大学 | 超高周疲劳裂纹萌生时间计算方法 |
CN106706275A (zh) * | 2015-07-29 | 2017-05-24 | 哈尔滨飞机工业集团有限责任公司 | 一种确定直升机复合材料涵道垂尾寿命的试验方法 |
CN106092785B (zh) * | 2016-06-17 | 2019-02-12 | 北京航空航天大学 | 采用不对称裂纹测试疲劳裂纹扩展速率的方法 |
CN106354934B (zh) * | 2016-08-30 | 2019-07-02 | 中航沈飞民用飞机有限责任公司 | 一种基于有限元分析求解的加筋壁板损伤容限分析方法 |
CN106596255A (zh) * | 2016-11-18 | 2017-04-26 | 北京航空航天大学 | 一种基于光纤光谱图像分析的孔边裂纹诊断方法 |
CN107389460B (zh) * | 2017-07-21 | 2019-12-24 | 中国工程物理研究院化工材料研究所 | 液压加载试验机的材料水平断裂检测方法 |
CN107884290B (zh) * | 2017-10-16 | 2019-11-19 | 北京理工大学 | 一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法 |
CN108982205B (zh) * | 2018-03-24 | 2019-04-05 | 北京工业大学 | 基于飞机实际结构件裂纹尖端局部实测应变的模拟件载荷转化方法 |
CN108984866B (zh) * | 2018-06-28 | 2022-09-30 | 中国铁道科学研究院集团有限公司金属及化学研究所 | 一种试验载荷谱的编制方法 |
CN111191391B (zh) * | 2019-12-10 | 2022-08-23 | 中国航空工业集团公司成都飞机设计研究所 | 一种基于局部应力分布加权系数的模拟件设计方法 |
CN112213090B (zh) * | 2020-09-25 | 2022-11-18 | 中国直升机设计研究所 | 一种直升机动部件损伤容限简化谱编制方法 |
CN114509260B (zh) * | 2021-11-29 | 2023-12-15 | 中国航发沈阳发动机研究所 | 航空发动机高速内外圈同转滚子轴承的加速等效试验方法 |
CN114778051B (zh) * | 2022-06-27 | 2022-09-02 | 中国飞机强度研究所 | 飞机垂向振动测试的试验载荷谱确定方法及应用 |
CN117332622B (zh) * | 2023-12-01 | 2024-02-13 | 无锡华天燃气轮机有限公司 | 一种裂纹扩展寿命确定方法、装置、设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129512A (zh) * | 2011-02-24 | 2011-07-20 | 西北工业大学 | 基于Paris公式的疲劳寿命分析方法 |
CN103544402A (zh) * | 2013-11-05 | 2014-01-29 | 中国航空工业集团公司西安飞机设计研究所 | 一种疲劳开裂结构等效分析谱的构造方法 |
CN103678256A (zh) * | 2012-09-20 | 2014-03-26 | 中国人民解放军装甲兵工程学院 | 车用发动机载荷谱编制方法 |
-
2014
- 2014-06-16 CN CN201410266931.6A patent/CN104034576B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102129512A (zh) * | 2011-02-24 | 2011-07-20 | 西北工业大学 | 基于Paris公式的疲劳寿命分析方法 |
CN103678256A (zh) * | 2012-09-20 | 2014-03-26 | 中国人民解放军装甲兵工程学院 | 车用发动机载荷谱编制方法 |
CN103544402A (zh) * | 2013-11-05 | 2014-01-29 | 中国航空工业集团公司西安飞机设计研究所 | 一种疲劳开裂结构等效分析谱的构造方法 |
Non-Patent Citations (2)
Title |
---|
含MSD 共线多孔平板应力强度因子有限元分析;郁大照等;《海军航空工程学院学报》;20060930;561-565 * |
载荷谱简化对广布疲劳损伤裂纹扩展影响研究;董登科,张文东;《2013 年首届中国航空科学技术大会论文集》;20131231;1-7 * |
Also Published As
Publication number | Publication date |
---|---|
CN104034576A (zh) | 2014-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104034576B (zh) | 一种飞机结构裂纹扩展试验用载荷谱简化方法 | |
CN107145641B (zh) | 叶片振动疲劳概率寿命预测方法 | |
Kishi et al. | Nonlinear finite element analysis of top-and seat-angle with double web-angle connections | |
CN103970999B (zh) | 飞机结构疲劳裂纹安全损伤扩展周期确定方法 | |
CN104020045B (zh) | 一种测定飞机/直升机动部件高周疲劳特性的步进加速试验方法 | |
US10989640B2 (en) | Method for defining threshold stress curves utilized in fatigue and damage tolerance analysis | |
CN101976311A (zh) | 基于漂移布朗运动模型的加速退化试验贝叶斯评估方法 | |
Metkar et al. | Evaluation of FEM based fracture mechanics technique to estimate life of an automotive forged steel crankshaft of a single cylinder diesel engine | |
CN111159943A (zh) | 一种动翼面封严结构的屈曲处理方法 | |
CN102944211B (zh) | 一种飞机整体翼梁止裂筋条面积的确定方法 | |
Lin et al. | Reliability analysis of aero-engine blades considering nonlinear strength degeneration | |
CN103544402A (zh) | 一种疲劳开裂结构等效分析谱的构造方法 | |
CN106372275A (zh) | 用于蠕变裂纹扩展数值模拟结果的修正方法 | |
Kato et al. | Data assimilation for turbulent flows | |
CN103323272A (zh) | 一种dfr值上限的确定方法 | |
CN104715112A (zh) | 批量化钻孔的方法 | |
CN112213090B (zh) | 一种直升机动部件损伤容限简化谱编制方法 | |
Kumar et al. | Pin shear–material testing and validation using finite element analysis (ANSYS) | |
CN104764421B (zh) | 一种重离子加速器单粒子试验降能片厚度快速确定方法 | |
Doherty | Transonic airfoil study using sonic plateau, optimization and off-design performance maps | |
Ju et al. | Dynamic fracture analysis technique of aircraft fuselage containing damage subjected to blast | |
Cho | Near-tip grid refinement for the effective and reliable natural element crack analysis | |
Dompierre et al. | Crack propagation methodology under complex loadings | |
Oliveira et al. | The influence of ratios and combined ratios on the distribution of the product of two independent gaussian random variables | |
Zaplana et al. | Fatigue and crack propagation life improvement quantification for coldworked holes outside nominal conditions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |