CN108062435B - 一种基于名义应力法的疲劳寿命校准方法 - Google Patents

一种基于名义应力法的疲劳寿命校准方法 Download PDF

Info

Publication number
CN108062435B
CN108062435B CN201711255528.3A CN201711255528A CN108062435B CN 108062435 B CN108062435 B CN 108062435B CN 201711255528 A CN201711255528 A CN 201711255528A CN 108062435 B CN108062435 B CN 108062435B
Authority
CN
China
Prior art keywords
life
nominal stress
median
value
fatigue limit
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
Application number
CN201711255528.3A
Other languages
English (en)
Other versions
CN108062435A (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.)
Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
Original Assignee
Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
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 Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC filed Critical Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
Priority to CN201711255528.3A priority Critical patent/CN108062435B/zh
Publication of CN108062435A publication Critical patent/CN108062435A/zh
Application granted granted Critical
Publication of CN108062435B publication Critical patent/CN108062435B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04—Ageing analysis or optimisation against ageing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种基于名义应力法的疲劳寿命校准方法，属于飞机强度设计技术领域，步骤一、对于同种材料，不同表面质量的试件S‑N曲线的参数α、A取值基本相同；确定理论疲劳极限C为需要优化的唯一敏感参数；步骤二、根据预先指定的载荷应力谱，确定其对应的中值寿命t50为目标进行优化；设定初始值，并不断调整理论疲劳极限C的取值，反复利用名义应力法计算中值寿命，直至计算寿命等于设定的目标中值t50，此时对应的C值即为修正后的S‑N曲线三参数式中的C值；步骤三、对确定的名义应力法及其对应的理论疲劳极限C进行验证；经验证采用优化得到的参数和名义应力法进行其他载荷谱下的寿命预测，使得其他载荷谱下的寿命计算结果具有较高的精度。

Description

一种基于名义应力法的疲劳寿命校准方法
技术领域
本发明属于飞机强度设计技术领域，具体涉及一种基于名义应力法的疲劳寿命校准方法。
背景技术
当前，对飞机疲劳关键部位进行随机谱下的寿命分析时，由于随机谱中存在复杂的载荷交互作用，因此目前尚未有一个普遍适用的高精度损伤累积模型能较好地预测其疲劳寿命。
但通过研究发现，系列化飞机往往在结构和载荷谱上具有很大的继承性，因此可以采用已有型号的疲劳试验结果，对后续型号的疲劳寿命进行预测。
发明内容
本发明的目的：为了解决上述问题，本发明提出了一种基于名义应力法的疲劳寿命校准方法，该方法可在变换随机载荷谱的情况下，仍保证计算所得的疲劳寿命具有较高的精度。
本发明的技术方案：一种基于名义应力法的疲劳寿命校准方法，所述名义应力法S-N曲线表达式为：
Figure GDA0003086174540000011
其中，C为理论疲劳极限；α，A为S-N曲线形状参数，S、N分别表示应力水平与失效寿命；
包括以下步骤：
步骤一、确定名义应力法中的敏感参数；
对于同种材料，不同表面质量的试件S-N曲线的参数α、A取值相同；
确定理论疲劳极限C为需要优化的唯一敏感参数；
步骤二、对理论疲劳极限C进行优化；
a)根据预先指定的载荷应力谱，确定其对应的中值寿命t50为目标进行优化；
b)设定初始值，并不断调整理论疲劳极限C的取值，反复利用名义应力法计算中值寿命，直至计算的寿命等于设定的目标中值t50，此时对应的C值即为修正后的S-N曲线三参数式中的C值；
步骤三、对确定的名义应力法及其对应的理论疲劳极限C进行验证；
利用优化得到的参数和名义应力法进行载荷谱下的寿命预测并将预测得到的中值寿命t’与所述载荷谱下的实测寿命中值比较误差；
若误差小于20％,则确定该校准方法。
本发明技术方案的有益技术效果：本发明在使用名义应力法进行疲劳寿命预测时，利用本发明中的校准方法，可在变换随机载荷谱的情况下，仍保证计算所得的疲劳寿命具有较高的精度。
附图说明
图1为本发明基于名义应力法的疲劳寿命校准方法的一优选实施例的流程示意图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行更加详细的描述。在附图中，自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例，而不是全部的实施例。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。下面结合附图对本发明的实施例进行详细说明。
一种基于名义应力法的疲劳寿命校准方法涉及到疲劳寿命校准方法及名义应力法；
疲劳寿命校准方法的核心是：以已有的试验结果为基础，采用合理的疲劳分析方法，确定疲劳分析的敏感参数，通过调整参数使得计算的寿命和试验寿命相同，然后采用该参数和相同的寿命分析方法进行其他载荷谱下的寿命预测，使得其他载荷谱下的寿命计算结果具有较高的精度。
名义应力法是依据结构疲劳关键部位的应力谱、关键部位的应力-寿命曲线(S-N曲线)与等寿命曲线，按照疲劳累积损伤理论估算结构的疲劳寿命。而疲劳累积损伤理论通常采用Miner线性累积损伤理论，该理论是理论界和工程界最普遍接受的，因此提高疲劳寿命预测结果精度的核心在于细节疲劳S-N曲线的确定。S-N曲线表达式采用三参数式，形式如下：
Figure GDA0003086174540000031
式中，C为理论疲劳极限；α，A为S-N曲线形状参数，S、N分别表示应力水平与失效寿命。将上式三个参数作为敏感参数进行优化，从而校准S-N曲线提高寿命计算精度。
在传统的名义应力法中，影响疲劳寿命分析结果准确性的主要因素包括累积损伤理论和细节的S-N曲线；上述两个因素中，由于目前缺乏普遍适用的非线性累积损伤理论，Miner线性累积损伤被理论界和工程界普遍接受。因此提高疲劳寿命预测结果精度的核心在于细节疲劳S-N曲线的确定，它应充分考虑材料、细节形式、细节几何、表面加工质量等因素的影响，建立充分反映细节疲劳特性的S-N曲线。
该曲线的建立是以材料S-N曲线参数为基础，考虑细节与材料的差别，对该S-N曲线参数进行修正，获得能够更好描述真实疲劳结构细节的S-N曲线。
修正方法为：以典型细节模拟试件在某实际载荷谱下的成组疲劳试验结果为依据，从相关材料手册中查找与结构关键部位应力集中系数相同的材料S-N曲线参数，作为细节S-N曲线参数的初始值，采用名义应力法进行该谱下的疲劳寿命分析，以计算的寿命与试验寿命相同为目标，优化得到细节S-N曲线参数。
如图1所示：一种基于名义应力法的疲劳寿命校准方法，包括以下步骤：
步骤一、确定名义应力法中的敏感参数；
优化参数时，参数越多，需要的参考样本数量就越多，且优化结果的稳健性越差。因此，应对优化参数进行适当约束。参考以往大量的试验结果发现，对于同种材料，不同表面质量的试件S-N曲线的参数α、A取值基本相同。因此本发明中校准方法将α、A取为定值，理论疲劳极限C为需要优化的唯一敏感参数。
步骤二、对理论疲劳极限C在某随机谱下试验寿命反推确定；
a)根据预先指定的载荷应力谱，确定其对应的中值寿命t50为目标进行优化；
b)设定初始值，并不断调整理论疲劳极限C的取值，反复利用名义应力法计算中值寿命，直至计算的寿命等于设定的目标中值t50，此时对应的C值即为修正后的S-N曲线三参数式中的C值；
本实施例中，可通过某随机谱下试验寿命反推确定。反推方法为：依据事先指定的载荷应力谱，并以其对应的中值寿命t50为目标，假定C0的初始值，并以一定步长不断调整C值，反复用名义应力法计算中值寿命，直至计算的寿命等于t50，所对应的C值即为结构关键部位S-N曲线三参数式中的C值。
步骤三、对确定的名义应力法及其对应的理论疲劳极限C进行验证；
利用优化得到的参数和名义应力法进行载荷谱下的寿命预测并将预测得到的中值寿命t’与所述载荷谱下的实测寿命中值比较误差，若误差小于20％,则确定该校准方法。
经验证采用优化得到的参数和名义应力法进行其他载荷谱下的寿命预测，使得其他载荷谱下的寿命计算结果具有较高的精度。
最后需要指出的是：以上实施例仅用以说明本发明的技术方案，而非对其限制。尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (1)

1.一种基于名义应力法的疲劳寿命校准方法，所述名义应力法S-N曲线表达式为：
Figure FDA0003005012570000011
其中，C为理论疲劳极限；α，A为S-N曲线形状参数，S、N分别表示应力水平与失效寿命；
其特征在于，包括以下步骤：
步骤一、确定名义应力法中的敏感参数；
对于同种材料，不同表面质量的试件S-N曲线的参数α、A取值相同；
确定理论疲劳极限C为需要优化的唯一敏感参数；
步骤二、对理论疲劳极限C进行优化；
a)根据预先指定的载荷应力谱，确定其对应的中值寿命t50为目标进行优化；
b)设定初始值，并不断调整理论疲劳极限C的取值，反复利用名义应力法计算中值寿命，直至计算的寿命等于设定的目标中值t50，此时对应的C值即为修正后的S-N曲线三参数式中的C值；
步骤三、对确定的名义应力法及其对应的理论疲劳极限C进行验证；
利用优化得到的参数和名义应力法进行载荷谱下的寿命预测并将预测得到的中值寿命t’与所述载荷谱下的实测寿命中值比较误差；
若误差小于20％,则确定该校准方法。
CN201711255528.3A 2017-12-01 2017-12-01 一种基于名义应力法的疲劳寿命校准方法 Active CN108062435B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711255528.3A CN108062435B (zh) 2017-12-01 2017-12-01 一种基于名义应力法的疲劳寿命校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711255528.3A CN108062435B (zh) 2017-12-01 2017-12-01 一种基于名义应力法的疲劳寿命校准方法

Publications (2)

Publication Number Publication Date
CN108062435A CN108062435A (zh) 2018-05-22
CN108062435B true CN108062435B (zh) 2021-07-09

Family

ID=62135203

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711255528.3A Active CN108062435B (zh) 2017-12-01 2017-12-01 一种基于名义应力法的疲劳寿命校准方法

Country Status (1)

Country Link
CN (1) CN108062435B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112229654B (zh) * 2020-09-30 2023-01-31 湖南行必达网联科技有限公司 板簧试验方法、装置、电子设备和计算机可读介质
CN113702062B (zh) * 2021-07-27 2023-10-13 岚图汽车科技有限公司 基于实际构件的疲劳特性测试方法及系统

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101718651A (zh) * 2009-11-18 2010-06-02 上海理工大学 强化和损伤共同作用下剩余强度和寿命的评价方法
US20120259593A1 (en) * 2011-04-07 2012-10-11 El-Zein Mohamad S Method for the prediction of fatigue life for welded structures
CN103344514B (zh) * 2013-07-05 2016-01-27 北京航空航天大学 一种基于名义应力法的高周疲劳和低强度冲击耦合的损伤计算方法
CN103838975A (zh) * 2014-03-18 2014-06-04 唐山轨道客车有限责任公司 车体焊缝疲劳寿命计算方法及装置
CN107103121B (zh) * 2017-04-01 2020-05-05 武汉理工大学 考虑焊接残余应力的钢桥构造细节s-n曲线的确定方法

Also Published As

Publication number Publication date
CN108062435A (zh) 2018-05-22

Similar Documents

Publication Publication Date Title
CN108108530B (zh) 一种适用于结构连接件的疲劳寿命校准方法
CN106886663B (zh) 齿轮弯曲疲劳寿命预测方法及装置
US9280620B2 (en) Method and system for probabilistic fatigue crack life estimation
CN108052717B (zh) 一种基于局部应力-应变法的疲劳寿命校准方法
CN109635385B (zh) 一种综合考虑疲劳强度影响因素的零部件寿命预测方法
CN107545110B (zh) 一种动态应力加速寿命试验剖面编制方法
JP2014071053A (ja) 高温部材のクリープ損傷評価方法および損傷評価システム
KR101526313B1 (ko) 피로수명 예측방법
CN108918261B (zh) 一种用少量试验测定材料构件疲劳寿命规律的方法
CN106844901B (zh) 一种基于多因素融合修正的结构件剩余强度评估方法
CN108062435B (zh) 一种基于名义应力法的疲劳寿命校准方法
WO2011026874A3 (en) String failure monitoring
CN112100713B (zh) 一种基于变间距栅格的索力自动识别方法
CN104990893A (zh) 一种基于相似判别的汽油辛烷值检测方法
CN111222095A (zh) 一种大坝变形监测中的粗差判别方法、装置及系统
CN103632035A (zh) 一种确定机械零部件的寿命概率分布与平均寿命的方法
CN116519206A (zh) 基于深度学习的压力传感器校准系统及方法
KR101428720B1 (ko) 비행자료를 이용한 응력 스펙트럼 생성 방법
CN108169013B (zh) 一种适用于多钉连接试件的疲劳寿命校准方法
CN104458630B (zh) 一种紫外差分气体分析仪的数据处理方法及系统
CN116805220A (zh) 一种基于熵理论的标准效能分析模型与评估方法
CN111022124A (zh) 一种桥隧工程短、长期变形的超前预警方法
Long et al. The effects of SHM system parameters on the value of damage detection information
JP2013044666A (ja) 多軸疲労寿命評価方法
CN113609675B (zh) 一种疲劳短裂纹扩展速率预测模型的确定方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant