CN112784357B - 一种孔边角裂纹应力强度因子确定的参数组合方法 - Google Patents

一种孔边角裂纹应力强度因子确定的参数组合方法 Download PDF

Info

Publication number
CN112784357B
CN112784357B CN202011610909.0A CN202011610909A CN112784357B CN 112784357 B CN112784357 B CN 112784357B CN 202011610909 A CN202011610909 A CN 202011610909A CN 112784357 B CN112784357 B CN 112784357B
Authority
CN
China
Prior art keywords
crack
stress intensity
intensity factor
hole
hole corner
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
CN202011610909.0A
Other languages
English (en)
Other versions
CN112784357A (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.)
AVIC First Aircraft Institute
Original Assignee
AVIC First Aircraft Institute
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 AVIC First Aircraft Institute filed Critical AVIC First Aircraft Institute
Priority to CN202011610909.0A priority Critical patent/CN112784357B/zh
Publication of CN112784357A publication Critical patent/CN112784357A/zh
Application granted granted Critical
Publication of CN112784357B publication Critical patent/CN112784357B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

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

Abstract

本发明属于航空结构设计领域,公开了一种孔边角裂纹应力强度因子确定的参数组合方法。通过开孔疲劳试验,确定了不同尺寸角裂纹相对穿透裂纹的修正因子,结合穿透裂纹应力强度因子计算公式,得到可用于快速确定孔边角裂纹应力强度因子的方法,解决孔边角裂纹应力强度因子无理论解且分析困难,缺乏高效可靠的孔边角裂纹应力强度因子确定方法的问题。

Description

一种孔边角裂纹应力强度因子确定的参数组合方法
技术领域
本发明属于航空结构设计领域,具体涉及一种孔边角裂纹应力强度因子确定的参数组合方法。
背景技术
飞机结构中广泛可见各类耳片结构,连接耳片承受循环载荷作用且载荷仅通过连接销钉传递,磨损和高应力集中引起循环载荷作用下的裂纹萌生,继而裂纹扩展直至破坏严重影响飞机结构的安全性。应力强度因子是计算裂纹扩展速率、剩余强度和疲劳寿命的一个重要参数。因此,要研究带裂纹耳片结构的裂纹扩展特性,首先就需要对耳片结构的应力强度因子进行研究。
孔边角裂纹应力强度因子无理论解且分析困难,以往只能通过有限元分析给出。
发明内容
本发明的目的是提供一种孔边角裂纹应力强度因子确定的参数组合方法,解决孔边角裂纹应力强度因子无理论解且分析困难,缺乏高效可靠的孔边角裂纹应力强度因子确定方法的问题。
为达到上述目的,本发明采用如下技术方案予以实现。
一种孔边角裂纹应力强度因子确定的参数组合方法,所述方法包括:
S1,计算孔边穿透裂纹应力强度因子;
S2,计算孔边角裂纹应力强度因子的修正系数;
S3,根据孔边角裂纹应力强度因子的修正系数以及孔边穿透裂纹应力强度因子确定孔边角裂纹应力强度因子。
本发明技术方案的特点和进一步的改进为:
(1)S1,穿透裂纹应力强度因子计算为:
获取不同的孔边穿透裂纹长度a1下,对应的穿透裂纹应力强度因子K:
其中σ是远端应力,a1是裂纹长度,β是综合修正因子。
(2)S2之前,所述方法还包括:
在不同的孔边角裂纹长度a2和裂纹深度c组合条件下,进行开孔结构的疲劳试验,得到不同的孔边角裂纹长度和裂纹深度组合条件实际的孔边角裂纹应力强度因子Kc
(3)在得到穿透裂纹应力强度因子K和实际的孔边角裂纹应力强度因子Kc的过程中,开孔结构的结构参数和外部载荷参数相同,开孔结构的结构参数包含:耳片结构的宽度、开孔直径、耳片厚度、孔边缘到耳片边缘的距离。
(4)在S2之前,所述方法还包括:
根据不同的孔边穿透裂纹长度a1下穿透裂纹应力强度因子K,在不同的孔边角裂纹长度a2和裂纹深度c组合条件下,实际的孔边角裂纹应力强度因子Kc;将裂纹长度a1与孔边角裂纹长度a2相同的两组数据作为一组试验数据,从而得到多组试验数据,并将相同的裂纹长度a1与孔边角裂纹长度a2均记为a。
(5)S2具体为:
按照对裂纹长度a、裂纹深度c进行无量纲处理,按照/>对裂纹深度c、开孔结构的厚度t进行无量纲处理;
在多组试验数据下,根据单参数组合法,假设孔边角裂纹长度a2与裂纹深度c之间的比例关系是固定的,仅有一个裂纹参数是独立的,按照KC=KΦ的关系确定Φ的拟合表达式。
(6)孔边角裂纹应力强度因子的修正系数Φ具体为:
式中a为裂纹深度;c为表面长度,t为开孔结构的厚度。
(7)S3,确定孔边角裂纹应力强度因子Kc0,具体为:
KC0=K0Φ
式中K0为任意开孔结构上穿透裂纹应力强度因子;Φ为角裂纹应力强度因子的修正系数,KC0同一开孔结构上与穿透裂纹具有相同裂纹长度的开孔角裂纹的应力强度因子。
本发明技术方案通过开孔疲劳试验,确定了不同尺寸角裂纹相对穿透裂纹的修正因子,结合穿透裂纹应力强度因子计算公式,得到可用于快速确定孔边角裂纹应力强度因子的方法。本发明通过开孔疲劳试验,确定了不同尺寸角裂纹相对穿透裂纹的修正因子,结合穿透裂纹应力强度因子计算公式,得到可用于快速确定孔边角裂纹应力强度因子的方法。本发明所提出的孔边角裂纹应力强度因子确定的参数组合方法步骤简便,解决孔边角裂纹应力强度因子无理论解且分析困难,缺乏高效可靠的孔边角裂纹应力强度因子确定方法的问题。
附图说明
图1为耳片角裂纹示意图;
图2为参数组合法与有限元计算结果对比图。
具体实施方式
下面结合附图对本发明做进一步详细说明。
本发明实施例提供一种孔边角裂纹应力强度因子确定的参数组合方法,所述方法包括:
S1,计算孔边穿透裂纹应力强度因子;
S2,计算孔边角裂纹应力强度因子的修正系数;
S3,根据孔边角裂纹应力强度因子的修正系数以及孔边穿透裂纹应力强度因子确定孔边角裂纹应力强度因子。
进一步的:
(1)S1,穿透裂纹应力强度因子计算为:
获取不同的孔边穿透裂纹长度a1下,对应的穿透裂纹应力强度因子K:
其中σ是远端应力,a1是裂纹长度,β是综合修正因子。
(2)S2之前,所述方法还包括:
在不同的孔边角裂纹长度a2和裂纹深度c组合条件下,进行开孔结构的疲劳试验,得到不同的孔边角裂纹长度和裂纹深度组合条件实际的孔边角裂纹应力强度因子KC
(3)在得到穿透裂纹应力强度因子K和实际的孔边角裂纹应力强度因子KC的过程中,开孔结构的结构参数和外部载荷参数相同,开孔结构的结构参数包含:耳片结构的宽度、开孔直径、耳片厚度、孔边缘到耳片边缘的距离。
(4)在S2之前,所述方法还包括:
根据不同的孔边穿透裂纹长度a1下穿透裂纹应力强度因子K,在不同的孔边角裂纹长度a2和裂纹深度c组合条件下,实际的孔边角裂纹应力强度因子KC;将裂纹长度a1与孔边角裂纹长度a2相同的两组数据作为一组试验数据,从而得到多组试验数据,并将相同的裂纹长度a1与孔边角裂纹长度a2均记为a。
(5)S2具体为:
按照对裂纹长度a、裂纹深度c进行无量纲处理,按照/>对裂纹深度c、开孔结构的厚度t进行无量纲处理;
在多组试验数据下,根据单参数组合法,假设孔边角裂纹长度a2与裂纹深度c之间的比例关系是固定的,仅有一个裂纹参数是独立的,按照KC=KΦ的关系确定Φ的拟合表达式。
(6)孔边角裂纹应力强度因子的修正系数Φ具体为:
式中a为裂纹深度;c为表面长度,t为开孔结构的厚度。
(7)S3,确定孔边角裂纹应力强度因子具体为:
KC0=K0Φ
式中K0为任意开孔结构上穿透裂纹应力强度因子;Φ为角裂纹应力强度因子的修正系数,KC0同一开孔结构上与穿透裂纹具有相同裂纹长度的开孔角裂纹的应力强度因子。
具体的,如图1所示,为耳片角裂纹示意图。给出了实施例中用到的参数。
如图2所示,为参数组合法与有限元计算结果对比图。
下面以某一具体实例对本发明做进一步详细说明。
给出角裂纹两组相关参数如下:
第一组参数:施加给耳片结构的载荷P=225KN,耳片结构的宽度W=60mm,开孔直径D=24mm,耳片厚度t=25mm,孔边缘到耳片边缘的距离e=18mm,a=6mm,c=6mm
第二组参数:P=225KN,W=60mm,D=24mm,t=25mm,e=18mm,a=10mm,c=10mm
步骤一、计算穿透裂纹应力强度因子
穿透裂纹应力强度因子计算为:β可由《民机结构耐久性与损伤容限设计手册》(下册)查得:
第一组参数:β=1.60,
第二组参数:β=1.31,
步骤二、计算角裂纹应力强度因子的修正系数
按照角裂纹应力强度因子的修正系数计算公式得
第一组参数:
第二组参数:
步骤三、计算角裂纹应力强度因子
按照如下公式计算角裂纹应力强度因子
第一组参数:
第二组参数:
式中K为穿透裂纹应力强度因子;Φ为角裂纹应力强度因子的修正系数。
通过三维有限元计算耳孔角裂纹应力强度因子,得到第一组数据的应力强度因子为第二组数据的应力强度因子为/>对比图如图2所示,可以看出本技术发明给出的孔边角裂纹应力强度因子确定方法与有限元计算结果吻合很好,本发明方法步骤简便,从而解决了孔边角裂纹应力强度因子无理论解且分析困难,缺乏高效可靠的孔边角裂纹应力强度因子确定方法的问题。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (2)

1.一种孔边角裂纹应力强度因子确定的参数组合方法,其特征在于,所述方法包括:
S1,计算孔边穿透裂纹应力强度因子;
S2,计算孔边角裂纹应力强度因子的修正系数;
S3,根据孔边角裂纹应力强度因子的修正系数以及孔边穿透裂纹应力强度因子确定孔边角裂纹应力强度因子;
S2之前,所述方法还包括:
在不同的孔边角裂纹长度a2和裂纹深度c组合条件下,进行开孔结构的疲劳试验,得到不同的孔边角裂纹长度和裂纹深度组合条件实际的孔边角裂纹应力强度因子Kc
在得到穿透裂纹应力强度因子K和实际的孔边角裂纹应力强度因子Kc的过程中,开孔结构的结构参数和外部载荷参数相同,开孔结构的结构参数包含:耳片结构的宽度、开孔直径、耳片厚度、孔边缘到耳片边缘的距离;
根据不同的孔边穿透裂纹长度a1下穿透裂纹应力强度因子K,在不同的孔边角裂纹长度a2和裂纹深度c组合条件下,实际的孔边角裂纹应力强度因子Kc;将裂纹长度a1与孔边角裂纹长度a2相同的两组数据作为一组试验数据,从而得到多组试验数据,并将相同的裂纹长度a1与孔边角裂纹长度a2均记为a;
S2具体为:按照对裂纹长度a、裂纹深度c进行无量纲处理,按照/>对裂纹深度c、开孔结构的厚度t进行无量纲处理;
在多组试验数据下,根据单参数组合法,假设孔边角裂纹长度a2与裂纹深度c之间的比例关系是固定的,仅有一个裂纹参数是独立的,按照Kc=KΦ的关系确定Φ的拟合表达式;
孔边角裂纹应力强度因子的修正系数Φ具体为:
式中为c裂纹深度;a为表面裂纹长度,t为开孔结构的厚度;
S3,确定孔边角裂纹应力强度因子Kc0,具体为:
Kc0=K0Φ
式中K0为任意开孔结构上穿透裂纹应力强度因子;Φ为角裂纹应力强度因子的修正系数,Kc0同一开孔结构上与穿透裂纹具有相同裂纹长度的开孔角裂纹的应力强度因子。
2.根据权利要求1所述的一种孔边角裂纹应力强度因子确定的参数组合方法,其特征在于,S1,穿透裂纹应力强度因子计算为:
获取不同的孔边穿透裂纹长度a1下,对应的穿透裂纹应力强度因子K:
其中σ是远端应力,a1是裂纹长度,β是综合修正因子。
CN202011610909.0A 2020-12-29 2020-12-29 一种孔边角裂纹应力强度因子确定的参数组合方法 Active CN112784357B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011610909.0A CN112784357B (zh) 2020-12-29 2020-12-29 一种孔边角裂纹应力强度因子确定的参数组合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011610909.0A CN112784357B (zh) 2020-12-29 2020-12-29 一种孔边角裂纹应力强度因子确定的参数组合方法

Publications (2)

Publication Number Publication Date
CN112784357A CN112784357A (zh) 2021-05-11
CN112784357B true CN112784357B (zh) 2024-02-13

Family

ID=75753916

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011610909.0A Active CN112784357B (zh) 2020-12-29 2020-12-29 一种孔边角裂纹应力强度因子确定的参数组合方法

Country Status (1)

Country Link
CN (1) CN112784357B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006006412A1 (ja) * 2004-07-09 2006-01-19 Kyushu Tlo Company, Limited 疲労き裂成長曲線の推定法、推定プログラムおよび推定装置
JP2007016517A (ja) * 2005-07-08 2007-01-25 Sekisui House Ltd 設計gl積算システム
CN102944211A (zh) * 2012-11-12 2013-02-27 中国航空工业集团公司西安飞机设计研究所 一种飞机整体翼梁止裂筋条面积的确定方法
CN109918701A (zh) * 2018-12-29 2019-06-21 北京航空航天大学 一种基于分段权函数的涡轮盘裂纹扩展模拟方法
CN111651924A (zh) * 2020-06-04 2020-09-11 北京航空航天大学 一种基于权函数的表面裂纹应力强度因子数据库建立方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006006412A1 (ja) * 2004-07-09 2006-01-19 Kyushu Tlo Company, Limited 疲労き裂成長曲線の推定法、推定プログラムおよび推定装置
JP2007016517A (ja) * 2005-07-08 2007-01-25 Sekisui House Ltd 設計gl積算システム
CN102944211A (zh) * 2012-11-12 2013-02-27 中国航空工业集团公司西安飞机设计研究所 一种飞机整体翼梁止裂筋条面积的确定方法
CN109918701A (zh) * 2018-12-29 2019-06-21 北京航空航天大学 一种基于分段权函数的涡轮盘裂纹扩展模拟方法
CN111651924A (zh) * 2020-06-04 2020-09-11 北京航空航天大学 一种基于权函数的表面裂纹应力强度因子数据库建立方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
Continuum damage mechanics combined with the extended finite element method for the total life prediction of a metallic component;Zhixin Zhan;International Journal of Mechanical Sciences;第48-58页 *
孔边角裂纹扩展特性分析;郦正能,何犹卿;航空学报;第15卷(第03期);第372-375页 *
孔边角裂纹的应力强度因子;曹桂馨;石油化工设备(第05期);第30-36页 *
弯曲载荷作用下孔边角裂纹的应力强度因子与疲劳扩展规律;巩建鸣;沈士明;戴树和;;南京工业大学学报(自然科学版);第12卷(第03期);第51-57页 *
钢中脆硬粒子裂纹形成机理;朱亮;物理学报;第53卷(第2期);第571-576页 *
高速铁路钢桥的低温脆性断裂评估;王应良;;中国铁道科学;第33卷(第03期);第28-34页 *

Also Published As

Publication number Publication date
CN112784357A (zh) 2021-05-11

Similar Documents

Publication Publication Date Title
US10094751B2 (en) System and method for determining direct damage tolerance allowables
CN111079329B (zh) 一种基于相似结构试验的疲劳寿命评估方法
JP5720798B2 (ja) 金属部材の疲労亀裂進展抑制方法及び疲労亀裂進展抑制された金属部材
CN112784357B (zh) 一种孔边角裂纹应力强度因子确定的参数组合方法
US10989640B2 (en) Method for defining threshold stress curves utilized in fatigue and damage tolerance analysis
CN109470549B (zh) 增材制造材料p-s-n曲线表征方法及其应用
CN107101756B (zh) 一种测量工件内部残余应力的方法
US11428614B2 (en) Notch treatment methods for flaw simulation
KR20130126014A (ko) 배관의 안정성 평가용 시험장치 및 시편의 제작방법
WO2020252842A1 (zh) 基于实际载荷的齿面接触斑点齿轮的设计方法
CN110287622B (zh) 一种加工表面广义微观应力集中现象的建模及分析方法
Louks et al. Static assessment of brittle/ductile notched materials: an engineering approach based on the Theory of Critical Distances
CN106595999B (zh) 基于aml方法的复合材料填孔拉伸强度设计许用值试验方法
CN112364541B (zh) 一种构件表面裂纹打磨后应力集中系数的计算方法
JP6989941B2 (ja) 亀裂補修用治具および亀裂補修方法
CN106599496B (zh) 基于数值计算的环芯法测量残余应力的方法
CN109959555B (zh) 考虑厚度和裂纹深度影响的中心裂纹拉伸试样的孕育期预测方法
CN107478518A (zh) 高温结构中含多椭圆形埋藏缺陷的简化处理方法
KR20210150562A (ko) 가스 터빈의 날개의 보수 방법 및 가스 터빈의 날개
CN113553706B (zh) 一种锥形橡胶弹簧粘接实验拐点的判定方法
CN112446154B (zh) 一种飞机机体结构疲劳破坏概率评估方法
CN112699479B (zh) 一种飞机结构裂纹形成寿命理论值分散系数的确定方法
CN115436033B (zh) 一种考虑变幅载荷作用的榫接多通道传力结构全寿命模拟方法
KR102121195B1 (ko) 파괴 인성 시험의 좌굴방지 지그
CN115533231B (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