CN108169013B - 一种适用于多钉连接试件的疲劳寿命校准方法 - Google Patents

一种适用于多钉连接试件的疲劳寿命校准方法 Download PDF

Info

Publication number
CN108169013B
CN108169013B CN201711250232.2A CN201711250232A CN108169013B CN 108169013 B CN108169013 B CN 108169013B CN 201711250232 A CN201711250232 A CN 201711250232A CN 108169013 B CN108169013 B CN 108169013B
Authority
CN
China
Prior art keywords
nail
fatigue limit
value
fatigue
life
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
CN201711250232.2A
Other languages
English (en)
Other versions
CN108169013A (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 CN201711250232.2A priority Critical patent/CN108169013B/zh
Publication of CN108169013A publication Critical patent/CN108169013A/zh
Application granted granted Critical
Publication of CN108169013B publication Critical patent/CN108169013B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational 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
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0073—Fatigue

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种适用于多钉连接试件的疲劳寿命校准方法，属于飞机强度设计技术领域。步骤一、以3钉连接件中值寿命与试验载荷谱为基础，采用应力严重系数法进行寿命计算，不断调整理论疲劳极限C值，使得计算寿命与试验寿命相同，得到校准值；步骤二、按照材料疲劳极限随应力集中系数的变化关系，计算得到3钉试件的理论疲劳极限值，并获得3钉理论疲劳极限值与步骤一中的疲劳极限校准值的修正系数K；步骤三、结合2钉和4钉的理论疲劳极限值及修正系数K,分别获得2钉和4钉的疲劳极限校准值。在使用应力严重系数法对连接结构进行疲劳寿命预测时，本发明可在改变连接试件钉数的情况下，仍保证计算所得的多钉连接试件的疲劳寿命具有较高的精度。

Description

一种适用于多钉连接试件的疲劳寿命校准方法
技术领域
本发明属于飞机强度设计技术领域，具体涉及一种适用于多钉连接试件的疲劳寿命校准方法。
背景技术
对多钉连接试件进行疲劳分析时，由于连接件存在钉传载荷，钉与孔间存在配合产生的复杂非线性效应，并且连接件中钉数的随机性，导致目前尚未有一个普遍适用的高精度损伤累积模型能较好地预测该结构形式下的疲劳寿命。
但通过研究发现，钉数不同的连接结构的疲劳寿命之间往往存在一定的相关性，因此可以通过个别已完成的连接结构的疲劳试验结果，对其他钉数的连接试件的疲劳寿命进行预测，这就是校准的基本思路。基于上述疲劳寿命校准的基本思路，本发明提出一种适用于多钉连接试件的疲劳寿命校准方法，利用该方法预测不同钉数时连接试件的疲劳寿命均具有较高的计算精度。
发明内容
本发明的目的：为了解决上述问题，本发明提出了一种适用于多钉连接试件的疲劳寿命校准方法，该方法可在改变连接试件钉数的情况下，仍保证计算所得的疲劳寿命具有较高的精度。
本发明的技术方案：一种适用于多钉连接试件的疲劳寿命校准方法，基于3钉连接件的试验结果建立多钉连接试件疲劳寿命校准方法，采用该校准方法计算2钉和4钉连接件的寿命；
包括以下步骤：
步骤一、以3钉连接件中值寿命与试验载荷谱为基础，采用应力严重系数法进行寿命计算，不断调整理论疲劳极限C值，使得计算寿命与试验寿命相同，得到校准值；
a)、获取对应于关键部位应力严重系数的S—N曲线；
Figure GDA0002358455220000021
其中，C为理论疲劳极限；α，A为S-N曲线形状参数；
b)、确定名义应力法中的敏感参数；
对于同种材料，不同表面质量的试件S-N曲线的参数α、A取值基本相同；
确定理论疲劳极限C为需要优化的唯一敏感参数；
c)对理论疲劳极限C进行优化；
根据预先指定的载荷应力谱，确定其对应的中值寿命t50为目标进行优化；
设定初始值，并不断调整理论疲劳极限C的取值，反复利用名义应力法计算中值寿命，直至计算寿命等于设定的目标中值t50，此时对应的C值即为修正后的S-N曲线三参数式中的C值；
步骤二、按照材料疲劳极限随应力集中系数的变化关系，计算得到3钉试件的理论疲劳极限值，并获得3钉理论疲劳极限值与步骤一中的疲劳极限校准值的修正系数K；
步骤三、结合2钉和4钉的理论疲劳极限值及修正系数K,分别获得2钉和4钉的疲劳极限校准值。
本发明技术方案的有益效果：本发明在使用应力严重系数法对连接结构进行疲劳寿命预测时，利用本发明中的校准方法，可在改变连接试件钉数的情况下，仍保证计算所得的多钉连接试件的疲劳寿命具有较高的精度。
附图说明
图1为本发明一种适用于多钉连接试件的疲劳寿命校准方法的一优选实施例的流程示意图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行更加详细的描述。在附图中，自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例，而不是全部的实施例。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。下面结合附图对本发明的实施例进行详细说明。
如图1所示：本发明一种适用于多钉连接试件的疲劳寿命校准方法，基于3钉连接件的试验结果建立多钉连接试件疲劳寿命校准方法，采用该校准方法计算2钉和4钉连接件的寿命。由于连接件存在钉传载荷，钉与孔间存在配合产生的复杂非线性效应，因此本发明基于应力严重系数法建立多钉连接试件的寿命校准方法。
应力严重系数法是一种属于名义应力法范畴的分析方法，其核心是确定连接件的应力严重系数(SSF)及对应的细节疲劳S-N曲线参数。通过有限元细节分析，得到钉孔处的旁路载荷、传递载荷，进而求出孔边的应力严重系数，该计算结果通常是没有争议的。因此寿命校准方法的关键是校准确定与材料、紧固件、加工质量、配合情况有关的细节疲劳S-N曲线参数，其S-N曲线表达式采用三参数式，形式如下：
Figure GDA0002358455220000031
式中，C为理论疲劳极限；α，A为S-N曲线形状参数。
将上式三个参数作为敏感参数进行优化，从而校准S-N曲线提高寿命计算精度。参考以往大量的试验结果发现，对于同种材料，不同表面质量的试件S-N曲线的参数α、A取值基本相同。因此本发明中校准方法将α、A取为定值，理论疲劳极限C为需要优化的唯一敏感参数。
多钉连接试件紧固件数量对疲劳性能的影响直接表现为对应力严重系数的影响，因此首先计算得到不同钉数量对应的应力严重系数，然后确定疲劳性能参数。在确定疲劳性能参数时，假定：疲劳S-N曲线校准的实质是对S-N曲线中的疲劳强度进行修正，即在材料S-N曲线的基础上，将疲劳极限乘以考虑细节几何、孔的加工质量、钉的装配工艺等因素影响的修正系数k。对按相同材料和工艺加工得到的多钉连接试件，假定疲劳强度修正系数k与钉数量无关，即不同数量钉连接件的理论疲劳极限与实际细节疲劳极限的比值k与钉的数量无关。
以3钉连接件中值寿命与试验载荷谱为基础，采用应力严重系数法进行寿命计算，不断调整C值，使得计算寿命与试验寿命相同，得到校准值。
同时按照材料疲劳极限随应力集中系数的变化关系，计算得到3钉试件的理论疲劳极限值，并将上述两个数值作比得到k值。
根据得到的k值并结合计算得到的2钉和4钉试件的理论疲劳极限值，即可得到校准后的2钉和4钉连接件的疲劳极限，并基于该疲劳性能参数对寿命进行预测。疲劳寿命校准具体流程如图1所示。
本实施例中，细节疲劳S-N曲线的校准过程如下:
1)疲劳S-N曲线的形式
S-N曲线表达式采用三参数式，形式如下：
Figure GDA0002358455220000041
式中，C为理论疲劳极限；α，A为S－N曲线形状参数。
2)取细节S-N曲线中的α，A为与SSF较为接近的Kt下的材料S-N曲线参数α，A值；
(3)以3钉连接件中值寿命与试验载荷谱为基础，采用应力严重系数法进行寿命计算，不断调整C值，使得计算寿命与试验寿命相同，得到校准的C3值；
(4)对其他多钉连接试件的S-N曲线疲劳极限，按以下步骤进行计算：
①按材料疲劳极限随应力集中系数的变化关系CM＝a+blgKt计算得到Kt等于3钉连接件的SSF时的CM，3值，即为3钉试件的理论疲劳极限值；
②计算校准的疲劳极限C3与①中理论疲劳极限的比值k，k＝C3/CM，3；
③按CM＝a+blgKt计算得到Kt分别等于2钉和4钉连接件的SSF时的CM，2值和CM，4值，即为2钉和4钉试件的理论疲劳极限值；
④将②中计算得到的k乘以③中计算得到的CM，2值和CM，4值，即为经过校准后的2钉和4钉连接件的疲劳极限C2、C4。
根据得到的k值并结合计算得到的2钉和4钉试件的理论疲劳极限值，即可得到校准后的2钉和4钉连接件的疲劳极限，并基于该疲劳性能参数对寿命进行预测。
最后需要指出的是：以上实施例仅用以说明本发明的技术方案，而非对其限制。尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (1)

1.一种适用于多钉连接试件的疲劳寿命校准方法，其特征在于：
基于3钉连接件的试验结果建立多钉连接试件疲劳寿命校准方法，采用该校准方法计算2钉和4钉连接件的寿命；
包括以下步骤：
步骤一、以3钉连接件中值寿命与试验载荷谱为基础，采用应力严重系数法进行寿命计算，不断调整理论疲劳极限C值，使得计算寿命与试验寿命相同，得到校准值；
a)、获取对应于关键部位应力严重系数的S—N曲线；
Figure FDA0002358455210000011
其中，C为理论疲劳极限；α，A为S-N曲线形状参数；
b)、确定名义应力法中的敏感参数；
对于同种材料，不同表面质量的试件S-N曲线的参数α、A取值基本相同；
确定理论疲劳极限C为需要优化的唯一敏感参数；
c)对理论疲劳极限C进行优化；
根据预先指定的载荷应力谱，确定其对应的中值寿命t50为目标进行优化；
设定初始值，并不断调整理论疲劳极限C的取值，反复利用名义应力法计算中值寿命，直至计算寿命等于设定的目标中值t50，此时对应的C值即为修正后的S-N曲线三参数式中的C值；
步骤二、按照材料疲劳极限随应力集中系数的变化关系，计算得到3钉试件的理论疲劳极限值，并获得3钉理论疲劳极限值与步骤一中的疲劳极限校准值的修正系数K；
步骤三、结合2钉和4钉的理论疲劳极限值及修正系数K,分别获得2钉和4钉的疲劳极限校准值。
CN201711250232.2A 2017-12-01 2017-12-01 一种适用于多钉连接试件的疲劳寿命校准方法 Active CN108169013B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711250232.2A CN108169013B (zh) 2017-12-01 2017-12-01 一种适用于多钉连接试件的疲劳寿命校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711250232.2A CN108169013B (zh) 2017-12-01 2017-12-01 一种适用于多钉连接试件的疲劳寿命校准方法

Publications (2)

Publication Number Publication Date
CN108169013A CN108169013A (zh) 2018-06-15
CN108169013B true CN108169013B (zh) 2020-04-14

Family

ID=62525035

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711250232.2A Active CN108169013B (zh) 2017-12-01 2017-12-01 一种适用于多钉连接试件的疲劳寿命校准方法

Country Status (1)

Country Link
CN (1) CN108169013B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110738000B (zh) * 2019-10-10 2022-09-20 中国航空工业集团公司沈阳飞机设计研究所 一种复杂受力状态下螺栓孔的高周疲劳寿命确定方法
CN111159879B (zh) * 2019-12-25 2023-07-21 中国航空工业集团公司西安飞机设计研究所 一种集中传载结构的钉载处理方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0227382D0 (en) * 2001-12-18 2002-12-31 Visteon Global Tech Inc Fatigue sensitivity determination procedure
CN101122560A (zh) * 2007-08-23 2008-02-13 南京航空航天大学 机械结构的裂纹扩展率和裂纹扩展寿命预测方法
CN102937520A (zh) * 2012-11-12 2013-02-20 上海理工大学 一种变幅载荷下的机械零构件疲劳寿命预测方法
CN103399992A (zh) * 2013-07-22 2013-11-20 中国兵器科学研究院 一种基于可靠寿命的结构耐久性优化设计方法
US8645086B1 (en) * 2010-09-18 2014-02-04 The Boeing Company Stress function calibration method
CN103927422A (zh) * 2014-04-24 2014-07-16 北京航空航天大学 一种基于渐进损伤模型的预测复合材料螺栓连接失效的强度包线法
CN104019971A (zh) * 2014-05-20 2014-09-03 北京航空航天大学 一种预测多钉连接件谱载疲劳寿命的方法
CN104537133A (zh) * 2014-05-12 2015-04-22 中国人民解放军空军工程大学 一种基于飞机结构寿命包线理论的单机剩余寿命预测方法
CN104809299A (zh) * 2015-05-05 2015-07-29 中国航空工业集团公司沈阳飞机设计研究所 一种基于耐久性严重载荷谱计算平均安全寿命的方法
CN105260584A (zh) * 2014-07-01 2016-01-20 中国人民解放军空军工程大学 一种已服役使用过的飞机结构剩余耐久性安全寿命确定方法
CN107103121A (zh) * 2017-04-01 2017-08-29 武汉理工大学 考虑焊接残余应力的钢桥构造细节s‑n曲线的确定方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9483605B2 (en) * 2013-03-14 2016-11-01 United Technologies Corporation Probabilistic high cycle fatigue (HCF) design optimization process

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0227382D0 (en) * 2001-12-18 2002-12-31 Visteon Global Tech Inc Fatigue sensitivity determination procedure
CN101122560A (zh) * 2007-08-23 2008-02-13 南京航空航天大学 机械结构的裂纹扩展率和裂纹扩展寿命预测方法
US8645086B1 (en) * 2010-09-18 2014-02-04 The Boeing Company Stress function calibration method
CN102937520A (zh) * 2012-11-12 2013-02-20 上海理工大学 一种变幅载荷下的机械零构件疲劳寿命预测方法
CN103399992A (zh) * 2013-07-22 2013-11-20 中国兵器科学研究院 一种基于可靠寿命的结构耐久性优化设计方法
CN103927422A (zh) * 2014-04-24 2014-07-16 北京航空航天大学 一种基于渐进损伤模型的预测复合材料螺栓连接失效的强度包线法
CN104537133A (zh) * 2014-05-12 2015-04-22 中国人民解放军空军工程大学 一种基于飞机结构寿命包线理论的单机剩余寿命预测方法
CN104019971A (zh) * 2014-05-20 2014-09-03 北京航空航天大学 一种预测多钉连接件谱载疲劳寿命的方法
CN105260584A (zh) * 2014-07-01 2016-01-20 中国人民解放军空军工程大学 一种已服役使用过的飞机结构剩余耐久性安全寿命确定方法
CN104809299A (zh) * 2015-05-05 2015-07-29 中国航空工业集团公司沈阳飞机设计研究所 一种基于耐久性严重载荷谱计算平均安全寿命的方法
CN107103121A (zh) * 2017-04-01 2017-08-29 武汉理工大学 考虑焊接残余应力的钢桥构造细节s‑n曲线的确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Statistical distributions of Transition Fatigue Strength and Transition Fatigue Life in duplex S–N fatigue curves;D.S. Paolino等;《Theoretical and Applied Fracture Mechanics 》;20150808;第80卷;第31-39页 *
多钉连接件的疲劳寿命预估及优化;李宝珠等;《中国数学力学物理学高新技术交叉研究学会会议论文集》;20100801;第374-378页 *

Also Published As

Publication number Publication date
CN108169013A (zh) 2018-06-15

Similar Documents

Publication Publication Date Title
CN108108530B (zh) 一种适用于结构连接件的疲劳寿命校准方法
CN111079329B (zh) 一种基于相似结构试验的疲劳寿命评估方法
CN107038292B (zh) 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法
CN108169013B (zh) 一种适用于多钉连接试件的疲劳寿命校准方法
Miranda et al. Numerical study of the Kardar-Parisi-Zhang equation
JP2007501517A5 (zh)
CN108052717B (zh) 一种基于局部应力-应变法的疲劳寿命校准方法
CN105046079B (zh) 一种基于d-最优内表设计的田口试验设计方法
CN102385646A (zh) Mos晶体管的器件失配的修正方法
CN111881564A (zh) 一种关于机械结构变幅疲劳寿命预测方法
CN112304478A (zh) 一种基于蠕变轮廓法的残余应力测试方法
CN108062435B (zh) 一种基于名义应力法的疲劳寿命校准方法
CN108009311A (zh) 一种针对蠕变试验的蠕变本构模型参数识别方法
CN102385641B (zh) 双极型晶体管的器件失配的修正方法
KR20000074677A (ko) 피로수명 예측 방법
CN107992709B (zh) 基于中间函数的热结构模型修正方法
CN112651153A (zh) 一种确定晶体塑性有限元模型材料参数的方法
CN111610146A (zh) 一种模拟脆性固体时离散元粘结参数的自动标定方法
CN102385642B (zh) 电阻的器件失配的修正方法
CN103105778B (zh) 一种工业过程仿真数学模型参数的估计方法
Bloembergen Analytical representations of the size-of-source effect
CN108960332A (zh) 一种基于多向主元素分析法的在线监测方法
CN103337467B (zh) 接触孔之底金属层去除量的量测方法
CN109766522B (zh) 一种刮板输送机链轮的非概率可靠性灵敏度分析方法
CN110706761B (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