CN112989659B - 一种基于点权函数法的表面裂纹强度因子数据库建立方法 - Google Patents

一种基于点权函数法的表面裂纹强度因子数据库建立方法 Download PDF

Info

Publication number
CN112989659B
CN112989659B CN202110259966.7A CN202110259966A CN112989659B CN 112989659 B CN112989659 B CN 112989659B CN 202110259966 A CN202110259966 A CN 202110259966A CN 112989659 B CN112989659 B CN 112989659B
Authority
CN
China
Prior art keywords
crack
point
stress
weight function
formula
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
CN202110259966.7A
Other languages
English (en)
Other versions
CN112989659A (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.)
Beihang University
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN202110259966.7A priority Critical patent/CN112989659B/zh
Publication of CN112989659A publication Critical patent/CN112989659A/zh
Application granted granted Critical
Publication of CN112989659B publication Critical patent/CN112989659B/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/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
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/08—Probabilistic or stochastic CAD
    • 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)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种基于点权函数的表面裂纹强度因子数据库建立方法，通过对表面裂纹面上所有点的应力和权函数进行积分，高速计算表面裂纹在复杂应力梯度下的应力强度因子，从而为表面概率损伤容限评估中的线弹性断裂力学分析提供基础。具体地，采用有限元参数化建模，建立立方体结构构型，基于J积分应力强度因子理论，计算给定应力梯度下的参考应力强度因子，计算精度不受裂纹尖端网格数量影响，计算结果准确；点权函数法相较于通用权函数法，克服了只考虑应力在裂纹深度单向变化的局限性，适用于更复杂的应力情况，弥补了应力强度因子手册的局限性，为航空发动机限寿件的概率损伤容限评估提供技术支撑，具有重要的工程意义和实际价值。

Description

一种基于点权函数法的表面裂纹强度因子数据库建立方法
技术领域
本发明属于航空发动机概率损伤容限评估技术领域，尤其涉及一种基于点权函数的表面裂纹强度因子数据库建立方法。
背景技术
航空发动机限寿件是航空发动机的关键安全件。在加工及使用维护过程中，不可避免地会在限寿件上引入小概率缺陷，这会对航空发动机的安全性带来极大的威胁。因而美国适航规章FAR-33部及我国适航规章CCAR-33部R2版第33.70条款提出一套加强的寿命管理流程，用概率损伤容限评估方法来考虑限寿件中出现概率极小的缺陷对安全性的影响。该评估方法将表面缺陷视为结构体中的表面裂纹，以线弹性断裂力学为基础，以应力强度因子作为裂纹应力集中的表征量，并结合缺陷、载荷、无损检测等随机因素，处理限寿件在带初始裂纹体的条件下低循环疲劳失效问题，从而计算航空发动机限寿件的失效风险，评估限寿件安全性。
高效的裂纹扩展计算方法是概率损伤容限评估的基础。航空发动机限寿件工作应力复杂，以航空发动机轮盘转子为例，其承受的工作载荷包括叶片拉力、离心力、热应力以及加工残余应力等，所承受的载荷在裂纹平面上具有任意应力分布。因而需要发展可计算任意应力分布下的裂纹应力强度因子的方法。工程常用的应力强度因子计算方法为解析法和数值计算法。解析法的基本原理是求解满足几何、平衡和物理方程与边界条件的应力函数，该方法只对简单结构体上具有常规应力分布的裂纹存在闭合解，因而具有一定的局限性。数值计算法主要采用有限元方法计算应力强度因子，可计算任意结构体中裂纹在任意应力分布下的应力强度因子，然而其效率受到建模和分网过程的影响，不适用于计算量大的概率损伤容限评估。
在此基础上，基于Buckner和Rice提出的权函数理论，可以结合解析法与数值计算法，建立表面裂纹应力强度因子数据库。其中，对于二维平面上应力分布主要随裂纹深度变化的裂纹，Glinka和Shen提出了通用权函数公式，从而可计算裂纹在裂纹深度方向上任意应力梯度下的应力强度因子。然而，通用权函数方法只适用于应力在裂纹深度单方向上变化的情况，而对于二维平面上应力分布在裂纹深度和长度两个维度上都不均匀的裂纹，通用权函数法并不能适用。实际部件中存在多处二维应力梯度，为解决具有复杂应力分布的部件裂纹扩展计算问题，Orynyak提出了点权函数公式，可高效计算具有复杂应力分布的裂纹的应力强度因子。点权函数法与通用权函数法相比，描述的应力分布形式从一维增加到二维，适用于存在复杂应力梯度情况的航空发动机限寿件，从而为限寿件概率损伤容限评估提供高效准确的线弹性断裂力学分析方法。
发明内容
有鉴于此，本发明提出一种基于点权函数的表面裂纹强度因子数据库建立方法，用以提高缺陷裂纹扩展计算效率，本发明的具体技术方案如下：
一种基于点权函数的表面裂纹强度因子数据库建立方法，包括以下步骤：
S1：针对具有表面裂纹的结构体，建立有限元模型，并在裂纹面上给定三种垂直于裂纹面的应力，作为有限元模型的边界条件；
S2：求解步骤S1建立的表面裂纹在三种应力边界条件下的参考应力强度因子；
S3：结合点权函数公式，利用三种应力边界条件与参考应力强度因子，求解点权函数公式中的权重系数；
S4：采用响应面分析法，将步骤S3获得的权重系数转化为权重系数代理模型；
S5：结合步骤S4得到的权重系数代理模型，对点权函数公式进行二维积分运算，快速获得二维任意载荷下的应力强度因子，建立表面裂纹应力强度因子数据库。
进一步地，所述步骤S1的具体过程为：
有限元模型涉及的四个参数为半圆形裂纹半径a、结构体厚度T、结构体宽度2L以及结构体半高H，无量纲化后简化为三个参数
Figure GDA0003457933140000021
和
Figure GDA0003457933140000022
取值范围分别为：
Figure GDA0003457933140000023
其中，
Figure GDA0003457933140000024
和
Figure GDA0003457933140000025
按照公式(1)的范围等间距取值，建立含有表面裂纹的结构体的有限元模型，并在每个有限元模型的裂纹面上给定两种垂直于裂纹面向外的应力，作为裂纹扩展驱动力，应力梯度分别为：
Figure GDA0003457933140000031
其中，σ1、σ2和σ3分别表示裂纹面上给定的三种应力梯度；σ0为任意正数，反映应力大小；x,y分别表示裂纹面内任意点的横坐标值和纵坐标值。
进一步地，所述步骤S2的具体过程为：
采用有限元软件ABAQUS计算裂纹尖端J积分：
Figure GDA0003457933140000032
其中，V表示位移矢量分量；ds表示积分线路上微分弧长；Ti和ui分别表示积分线路上作用于ds积分单位上任意指定方向i方向的应力分量与位移分量；Γ表示裂纹下表面某点到裂纹上表面某点的积分回路；x'和y'分别表示裂纹尖端积分区域的横坐标值和纵坐标值；
对于单一的I型裂纹，对有限元模型求解得到的参考应力强度因子Kr为：
Figure GDA0003457933140000033
对于平面应力：E′＝E；对于平面应变：
Figure GDA0003457933140000034
其中，E表示弹性模量，ν表示泊松比；Kr包含半圆形裂纹边界上所有点的参考应力强度因子，采用三种应力梯度σ1、σ2和σ3作为边界条件，通过公式(4)计算得到半圆形裂纹的最深点B和表面点A在两种应力边界条件下的参考应力强度因子
Figure GDA0003457933140000035
和
Figure GDA0003457933140000036
进一步地，所述步骤S3的具体过程为：
点权函数理论公式为：
Figure GDA0003457933140000037
Figure GDA0003457933140000041
其中，l＝1,2,3；k＝A,B，A是半圆形裂纹表面点，B是半圆形裂纹最深点；
Figure GDA0003457933140000042
表示三种应力梯度下的参考应力强度因子；σl是公式(2)中给定的位于裂纹面上的三种参考应力梯度；Q是裂纹面上任意一点，r代表Q到裂纹圆心的距离；Qy是Q点关于y轴对称点；Q'是裂纹尖端任意一点，Q'y是Q'的纵坐标值，sign(Q'y)是该纵坐标值的符号函数；
Figure GDA0003457933140000043
是裂纹面上Q点与裂纹尖端Q'的距离的平方；
Figure GDA0003457933140000044
是裂纹关于y轴对称点Qy与裂纹尖端Q'的距离的平方；
Figure GDA0003457933140000045
是半圆形裂纹表面点A的权重系数，
Figure GDA0003457933140000046
是半圆形裂纹最深点B的权重系数；
将点权函数公式进行柱坐标系变换，得到：
Figure GDA0003457933140000047
其中，
Figure GDA0003457933140000048
分别是裂纹面上任意点与y轴负方向的夹角，以及裂纹面任意点与裂纹圆心的距离，θ是裂纹尖端点Q'与y轴负方向的夹角；
简化点权函数公式：
Figure GDA0003457933140000049
Figure GDA00034579331400000410
Figure GDA0003457933140000051
l通过将公式(2)以及参考应力强度因子
Figure GDA0003457933140000052
和
Figure GDA0003457933140000053
代入公式(9)和(10)，通过求解方程组(8)，获取点权函数公式中的半圆形裂纹表面点A的权重系数
Figure GDA0003457933140000054
以及半圆形裂纹最深点B的权重系数
Figure GDA0003457933140000055
进一步地，所述步骤S4的具体过程为：
采用响应面分析法，获取参考应力强度因子代理模型：
Figure GDA0003457933140000056
其中k＝A,B；n＝α,β,γ；
用含有交叉项的线性多项式表达：
Figure GDA0003457933140000057
其中，i＝1,2,3；j＝1,2,3；po、pj和pij分别表示待定系数，ei和ej均表示无量纲几何参数，e1、e2和e3分别代表
Figure GDA0003457933140000058
和
Figure GDA0003457933140000059
ε表示误差项；
通过回归方法求解出公式(7)中的待定系数，记为
Figure GDA00034579331400000510
和
Figure GDA00034579331400000511
则公式(7)表达为：
Figure GDA00034579331400000512
进一步地，所述步骤S5的具体过程为：
权重系数
Figure GDA00034579331400000513
表示为与结构体几何参数
Figure GDA00034579331400000514
和
Figure GDA00034579331400000515
相关的函数形式：
Figure GDA00034579331400000516
当结构体几何参数
Figure GDA0003457933140000061
和
Figure GDA0003457933140000062
取值在公式(1)范围内时，由公式(14)得到权重系数
Figure GDA0003457933140000063
给定作用于裂纹面的任意应力梯度σ，利用点权函数公式(14)和公式(15)，通过积分计算，快速获取半圆形裂纹的最深点B和表面点A的应力强度因子：
Figure GDA0003457933140000064
得到结构体上的表面裂纹在任意应力梯度下的应力强度因子，建立不同结构体的表面裂纹应力强度因子数据库。
本发明的有益效果在于：
1.本发明的方法通过对表面裂纹面上所有点的应力和权函数进行积分，高速计算表面裂纹在复杂应力梯度下的应力强度因子，从而为表面概率损伤容限评估中的线弹性断裂力学分析提供基础。
2.本发明计算给定应力梯度下的参考应力强度因子，计算精度不受裂纹尖端网格数量影响，计算结果准确。
3.本发明基于点权函数理论，建立表面裂纹在复杂应力梯度下的应力强度因子高效计算方法，弥补了应力强度因子手册的局限性；与通用权函数法相比，通用权函数法可以获取一维任意应力梯度下的裂纹应力强度因子，点权函数法可以获得二维任意应力梯度下的裂纹应力强度因子，因而点权函数法适用的范围更广。
4.本发明点权函数数据库在适用的结构体构型中，对于一次二次分布应力梯度计算误差在5％以内。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，通过参考附图会更加清楚的理解本发明的特征和优点，附图是示意性的而不应理解为对本发明进行任何限制，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，可以根据这些附图获得其他的附图。其中：
图1为本发明的方法流程图；
图2为具有表面裂纹的结构体示意图；
图3为具有表面裂纹面示意图；
图4为本发明实施的流程示意图；
图5为本发明实施例1的结构体半圆形裂纹面应力分布图，其中(a)为应力为σv1时的裂纹面上应力分布，(b)为应力为σv2时的裂纹面上应力分布；
图6为本发明实施例1的结构体半圆形裂纹面通过本发明的方法与有限元方法计算结果对比，其中(a)为应力为σv1时的计算结果，(b)为应力为σv2时的计算结果。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本发明的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
本发明提出一种基于点权函数表面裂纹强度因子数据库建立方法，通过对表面裂纹面上所有点的应力和权函数进行积分，高速计算表面裂纹在复杂应力梯度下的应力强度因子，从而为表面概率损伤容限评估中的线弹性断裂力学分析提供基础。
具体地，如图1-4所示，本发明采用有限元参数化建模，建立立方体结构构型，基于J积分应力强度因子理论，计算给定应力梯度下的参考应力强度因子，计算精度不受裂纹尖端网格数量影响，计算结果准确；基于点权函数理论，建立表面裂纹在复杂应力梯度下的应力强度因子高效计算方法，弥补了应力强度因子手册的局限性；与通用权函数法相比，通用权函数法可以获取一维任意应力梯度下的裂纹应力强度因子，点权函数法可以获得二维任意应力梯度下的裂纹应力强度因子，因而点权函数法适用的范围更广。
具体地，一种基于点权函数的表面裂纹强度因子数据库建立方法，包括以下步骤：
S1：针对具有表面裂纹的结构体，建立有限元模型，并在裂纹面上给定三种垂直于裂纹面的应力，作为有限元模型的边界条件；
S2：求解步骤S1建立的表面裂纹在三种应力边界条件下的参考应力强度因子；
S3：结合点权函数公式，利用三种应力边界条件与参考应力强度因子，求解点权函数公式中的权重系数；
S4：采用响应面分析法，将步骤S3获得的权重系数转化为权重系数代理模型；
S5：结合步骤S4得到的权重系数代理模型，对点权函数公式进行二维积分运算，快速获得二维任意载荷下的应力强度因子，建立表面裂纹应力强度因子数据库。
步骤S1的具体过程为：
有限元模型涉及的四个参数为半圆形裂纹半径a、结构体厚度T、结构体宽度2L以及结构体半高H，无量纲化后简化为三个参数
Figure GDA0003457933140000081
和
Figure GDA0003457933140000082
取值范围分别为：
Figure GDA0003457933140000083
其中，
Figure GDA0003457933140000084
和
Figure GDA0003457933140000085
按照公式(1)的范围等间距取值，建立含有表面裂纹的结构体的有限元模型，并在每个有限元模型的裂纹面上给定两种垂直于裂纹面向外的应力，作为裂纹扩展驱动力，应力梯度分别为：
Figure GDA0003457933140000086
其中，σ1、σ2和σ3分别表示裂纹面上给定的三种应力梯度；σ0为任意正数，反映应力大小；x,y分别表示裂纹面内任意点的横坐标值和纵坐标值。
步骤S2的具体过程为：
采用有限元软件ABAQUS计算裂纹尖端J积分：
Figure GDA0003457933140000087
其中，V表示位移矢量分量；ds表示积分线路上微分弧长；Ti和ui分别表示积分线路上作用于ds积分单位上任意指定方向i方向的应力分量与位移分量；Γ表示裂纹下表面某点到裂纹上表面某点的积分回路；x'和y'分别表示裂纹尖端积分区域的横坐标值和纵坐标值；
对于单一的I型裂纹，对有限元模型求解得到的参考应力强度因子Kr为：
Figure GDA0003457933140000091
对于平面应力：E′＝E；对于平面应变：
Figure GDA0003457933140000092
其中，E表示弹性模量，ν表示泊松比；Kr包含半圆形裂纹边界上所有点的参考应力强度因子，采用三种应力梯度σ1、σ2和σ3作为边界条件，通过公式(4)计算得到半圆形裂纹的最深点B和表面点A在两种应力边界条件下的参考应力强度因子
Figure GDA0003457933140000093
和
Figure GDA0003457933140000094
步骤S3的具体过程为：
点权函数理论公式为：
Figure GDA0003457933140000095
Figure GDA0003457933140000096
其中，l＝1,2,3；k＝A,B，A是半圆形裂纹表面点，B是半圆形裂纹最深点；
Figure GDA00034579331400000911
表示三种应力梯度下的参考应力强度因子；σl是公式(2)中给定的位于裂纹面上的三种参考应力梯度；Q是裂纹面上任意一点，r代表Q到裂纹圆心的距离；Qy是Q点关于y轴对称点；Q'是裂纹尖端任意一点，Q'y是Q'的纵坐标值，sign(Q'y)是该纵坐标值的符号函数；
Figure GDA0003457933140000097
是裂纹面上Q点与裂纹尖端Q'的距离的平方；
Figure GDA0003457933140000098
是裂纹关于y轴对称点Qy与裂纹尖端Q'的距离的平方；
Figure GDA0003457933140000099
是半圆形裂纹表面点A的权重系数，
Figure GDA00034579331400000910
是半圆形裂纹最深点B的权重系数；
将点权函数公式进行柱坐标系变换，得到：
Figure GDA0003457933140000101
其中，
Figure GDA0003457933140000102
分别是裂纹面上任意点与y轴负方向的夹角，以及裂纹面任意点与裂纹圆心的距离，θ是裂纹尖端点Q'与y轴负方向的夹角；
简化点权函数公式：
Figure GDA0003457933140000103
Figure GDA0003457933140000104
Figure GDA0003457933140000105
l通过将公式(2)以及参考应力强度因子
Figure GDA0003457933140000106
和
Figure GDA0003457933140000107
代入公式(9)和(10)，通过求解方程组(8)，获取点权函数公式中的半圆形裂纹表面点A的权重系数
Figure GDA0003457933140000108
以及半圆形裂纹最深点B的权重系数
Figure GDA0003457933140000109
步骤S4的具体过程为：
采用响应面分析法，获取参考应力强度因子代理模型：
Figure GDA00034579331400001010
其中k＝A,B；n＝α,β,γ；
用含有交叉项的线性多项式表达：
Figure GDA0003457933140000111
其中，i＝1,2,3；j＝1,2,3；po、pj和pij分别表示待定系数，ei和ej均表示无量纲几何参数，e1、e2和e3分别代表
Figure GDA0003457933140000112
和
Figure GDA0003457933140000113
ε表示误差项；
通过回归方法求解出公式(7)中的待定系数，记为
Figure GDA0003457933140000114
和
Figure GDA0003457933140000115
则公式(7)表达为：
Figure GDA0003457933140000116
步骤S5的具体过程为：
权重系数
Figure GDA0003457933140000117
表示为与结构体几何参数
Figure GDA0003457933140000118
和
Figure GDA0003457933140000119
相关的函数形式：
Figure GDA00034579331400001110
当结构体几何参数
Figure GDA00034579331400001111
和
Figure GDA00034579331400001112
取值在公式(1)范围内时，由公式(14)得到权重系数
Figure GDA00034579331400001113
给定作用于裂纹面的任意应力梯度σ，利用点权函数公式(14)和公式(15)，通过积分计算，快速获取半圆形裂纹的最深点B和表面点A的应力强度因子：
Figure GDA00034579331400001114
得到结构体上的表面裂纹在任意应力梯度下的应力强度因子，建立不同结构体的表面裂纹应力强度因子数据库。
为了方便理解本发明的上述技术方案，以下通过具体实施例对本发明的上述技术方案进行详细说明。
实施例1
对表面半圆形裂纹给定不同应力梯度，如图5所示，其中，(a)为应力为σv1时的裂纹面上应力分布，(b)为应力为σv2时的裂纹面上应力分布，
Figure GDA0003457933140000121
Figure GDA0003457933140000122
其中，σv1及σv2为裂纹面上应力分布，为以半圆裂纹圆心为原点，结构体厚度方向为X轴，结构体宽度方向为Y轴，x,y分别为裂纹面任意点的横坐标、纵坐标值。
根据本发明的方法，可获取表面裂纹在不同的半径a以及不同结构体宽度及厚度L、T下的应力强度因子。此外，根据有限元建模方法，可以获取应力强度因子精确值，如图6(a)和图6(b)所示，可见本发明的方法具有较高的计算精度，在适用尺寸范围中，对于一次二次多项式裂纹，其误差在5％以内。
综上，本发明的点权函数数据库在适用的结构体构型中，对于一次二次分布应力梯度计算误差在5％以内。本发明提供的上述基于权函数的表面裂纹应力强度因子数据库建立方法，为航空发动机限寿件的概率损伤容限评估提供技术支撑，具有重要的工程意义和实际价值。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于点权函数的表面裂纹强度因子数据库建立方法，其特征在于，包括以下步骤：
S1：针对具有表面裂纹的结构体，建立有限元模型，并在裂纹面上给定三种垂直于裂纹面的应力，作为有限元模型的边界条件；具体过程为：
有限元模型涉及的四个参数为半圆形裂纹半径a、结构体厚度T、结构体宽度2L以及结构体半高H，无量纲化后简化为三个参数
Figure FDA0003464346280000011
和
Figure FDA0003464346280000012
取值范围分别为：
Figure FDA0003464346280000013
其中，
Figure FDA0003464346280000014
和
Figure FDA0003464346280000015
按照公式(1)的范围等间距取值，建立含有表面裂纹的结构体的有限元模型，并在每个有限元模型的裂纹面上给定两种垂直于裂纹面向外的应力，作为裂纹扩展驱动力，应力梯度分别为：
Figure FDA0003464346280000016
其中，σ1、σ2和σ3分别表示裂纹面上给定的三种应力梯度；σ0为任意正数，反映应力大小；x,y分别表示裂纹面内任意点的横坐标值和纵坐标值；
S2：求解步骤S1建立的表面裂纹在三种应力边界条件下的参考应力强度因子；具体过程为：
采用有限元软件ABAQUS计算裂纹尖端J积分：
Figure FDA0003464346280000017
其中，V表示位移矢量分量；ds表示积分线路上微分弧长；Ti和ui分别表示积分线路上作用于ds积分单位上任意指定方向i方向的应力分量与位移分量；Γ表示裂纹下表面某点到裂纹上表面某点的积分回路；x'和y'分别表示裂纹尖端积分区域的横坐标值和纵坐标值；
对于单一的I型裂纹，对有限元模型求解得到的参考应力强度因子Kr为：
Figure FDA0003464346280000018
对于平面应力：E′＝E；对于平面应变：
Figure FDA0003464346280000021
其中，E表示弹性模量，ν表示泊松比；Kr包含半圆形裂纹边界上所有点的参考应力强度因子，采用三种应力梯度σ1、σ2和σ3作为边界条件，通过公式(4)计算得到半圆形裂纹的最深点B和表面点A在两种应力边界条件下的参考应力强度因子
Figure FDA0003464346280000022
和
Figure FDA0003464346280000023
S3：结合点权函数公式，利用三种应力边界条件与参考应力强度因子，求解点权函数公式中的权重系数；具体过程为：
点权函数理论公式为：
Figure FDA0003464346280000024
Figure FDA0003464346280000025
其中，l＝1,2,3；k＝A,B，A是半圆形裂纹表面点，B是半圆形裂纹最深点；
Figure FDA00034643462800000212
表示三种应力梯度下的参考应力强度因子；σl是公式(2)中给定的位于裂纹面上的三种参考应力梯度；Q是裂纹面上任意一点，r代表Q到裂纹圆心的距离；Qy是Q点关于y轴对称点；Q'是裂纹尖端任意一点，Q'y是Q'的纵坐标值，sign(Q'y)是该纵坐标值的符号函数；
Figure FDA0003464346280000026
是裂纹面上Q点与裂纹尖端Q'的距离的平方；
Figure FDA0003464346280000027
是裂纹关于y轴对称点Qy与裂纹尖端Q'的距离的平方；
Figure FDA0003464346280000028
是半圆形裂纹表面点A的权重系数，
Figure FDA0003464346280000029
是半圆形裂纹最深点B的权重系数；
将点权函数公式进行柱坐标系变换，得到：
Figure FDA00034643462800000210
其中，
Figure FDA00034643462800000211
分别是裂纹面上任意点与y轴负方向的夹角，以及裂纹面任意点与裂纹圆心的距离，θ是裂纹尖端点Q'与y轴负方向的夹角；
简化点权函数公式：
Figure FDA0003464346280000031
Figure FDA0003464346280000032
Figure FDA0003464346280000033
l通过将公式(2)以及参考应力强度因子
Figure FDA0003464346280000034
和
Figure FDA0003464346280000035
代入公式(9)和(10)，通过求解方程组(8)，获取点权函数公式中的半圆形裂纹表面点A的权重系数
Figure FDA0003464346280000036
Figure FDA0003464346280000037
以及半圆形裂纹最深点B的权重系数
Figure FDA0003464346280000038
S4：采用响应面分析法，将步骤S3获得的权重系数转化为权重系数代理模型；
S5：结合步骤S4得到的权重系数代理模型，对点权函数公式进行二维积分运算，快速获得二维任意载荷下的应力强度因子，建立表面裂纹应力强度因子数据库。
2.根据权利要求1所述的一种基于点权函数的表面裂纹强度因子数据库建立方法，其特征在于，所述步骤S4的具体过程为：
采用响应面分析法，获取参考应力强度因子代理模型：
Figure FDA0003464346280000039
其中k＝A,B；n＝α,β,γ；
用含有交叉项的线性多项式表达：
Figure FDA0003464346280000041
其中，i＝1,2,3；j＝1,2,3；po、pj和pij分别表示待定系数，ei和ej均表示无量纲几何参数，e1、e2和e3分别代表
Figure FDA0003464346280000042
和
Figure FDA0003464346280000043
ε表示误差项；
通过回归方法求解出公式(7)中的待定系数，记为
Figure FDA0003464346280000044
和
Figure FDA0003464346280000045
则公式(7)表达为：
Figure FDA0003464346280000046
3.根据权利要求1或2所述的一种基于点权函数的表面裂纹强度因子数据库建立方法，其特征在于，所述步骤S5的具体过程为：
权重系数
Figure FDA0003464346280000047
表示为与结构体几何参数
Figure FDA0003464346280000048
和
Figure FDA0003464346280000049
相关的函数形式：
Figure FDA00034643462800000410
当结构体几何参数
Figure FDA00034643462800000411
和
Figure FDA00034643462800000412
取值在公式(1)范围内时，由公式(14)得到权重系数
Figure FDA00034643462800000413
给定作用于裂纹面的任意应力梯度σ，利用点权函数公式(14)和公式(15)，通过积分计算，快速获取半圆形裂纹的最深点B和表面点A的应力强度因子：
Figure FDA00034643462800000414
得到结构体上的表面裂纹在任意应力梯度下的应力强度因子，建立不同结构体的表面裂纹应力强度因子数据库。
CN202110259966.7A 2021-03-10 2021-03-10 一种基于点权函数法的表面裂纹强度因子数据库建立方法 Active CN112989659B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110259966.7A CN112989659B (zh) 2021-03-10 2021-03-10 一种基于点权函数法的表面裂纹强度因子数据库建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110259966.7A CN112989659B (zh) 2021-03-10 2021-03-10 一种基于点权函数法的表面裂纹强度因子数据库建立方法

Publications (2)

Publication Number Publication Date
CN112989659A CN112989659A (zh) 2021-06-18
CN112989659B true CN112989659B (zh) 2022-03-04

Family

ID=76336258

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110259966.7A Active CN112989659B (zh) 2021-03-10 2021-03-10 一种基于点权函数法的表面裂纹强度因子数据库建立方法

Country Status (1)

Country Link
CN (1) CN112989659B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113591355B (zh) * 2021-08-06 2022-03-11 中山政数大数据科技有限公司 基于大数据的桥梁拉索钢丝腐蚀程度智能自动测量平台
CN113610812B (zh) * 2021-08-10 2024-02-02 哈尔滨理工大学 超声振动磨削光学玻璃亚表面裂纹形式判别与影响因素的权重解算方法
CN114511534B (zh) * 2022-01-28 2023-05-05 江苏泰和木业有限公司 一种基于图像处理的pc板裂纹判断方法及系统
CN117057166B (zh) * 2023-10-11 2023-12-26 合肥通用机械研究院有限公司 应力集中部位裂纹自由表面处应力强度因子的计算方法
CN117057167B (zh) * 2023-10-11 2024-01-09 合肥通用机械研究院有限公司 一种应力集中部位裂纹最深点处应力强度因子的计算方法
CN117236069B (zh) * 2023-11-08 2024-02-02 合肥通用机械研究院有限公司 任意应力分布下裂纹自由表面处应力强度因子的计算方法
CN117195608B (zh) * 2023-11-08 2024-01-12 合肥通用机械研究院有限公司 任意应力分布下裂纹最深点处应力强度因子的计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111191391A (zh) * 2019-12-10 2020-05-22 中国航空工业集团公司成都飞机设计研究所 一种基于局部应力分布加权系数的模拟件设计方法
CN111651924A (zh) * 2020-06-04 2020-09-11 北京航空航天大学 一种基于权函数的表面裂纹应力强度因子数据库建立方法
CN111783243A (zh) * 2020-06-18 2020-10-16 东南大学 一种基于滤波算法的金属结构疲劳裂纹扩展寿命预测方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107208842A (zh) * 2014-10-10 2017-09-26 红叶资源公司 气体围闭系统
CN109918701B (zh) * 2018-12-29 2021-09-10 北京航空航天大学 一种基于分段权函数的涡轮盘裂纹扩展模拟方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111191391A (zh) * 2019-12-10 2020-05-22 中国航空工业集团公司成都飞机设计研究所 一种基于局部应力分布加权系数的模拟件设计方法
CN111651924A (zh) * 2020-06-04 2020-09-11 北京航空航天大学 一种基于权函数的表面裂纹应力强度因子数据库建立方法
CN111783243A (zh) * 2020-06-18 2020-10-16 东南大学 一种基于滤波算法的金属结构疲劳裂纹扩展寿命预测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"Engineering Method For Stress Intensity Factor Caculation And Their Practical Application";I.V. Orynyak;《17th European Conference on Fracture》;20080905;1911-1923 *
"Universal Weight Function Method on the Probabilistic Surface Damage Tolerance Assessment of Aeroengine Rotors";Huimin Zhou .etc;《 ASME 2020 International Mechanical Engineering Congress and Exposition》;20210216;1-10 *
"大深长比表面裂纹应力强度因子的二维权函数";袁奎霖等;《船舶力学》;20201130;第24卷(第11期);1453-1461 *
"航空发动机限寿件概率失效风险评估的等效应力转化";丁水汀等;《燃气涡轮试验与研究》;20190228;第32卷(第1期);1-7,19 *

Also Published As

Publication number Publication date
CN112989659A (zh) 2021-06-18

Similar Documents

Publication Publication Date Title
CN112989659B (zh) 一种基于点权函数法的表面裂纹强度因子数据库建立方法
CN111651924B (zh) 一种基于权函数的表面裂纹应力强度因子数据库建立方法
CN109614715B (zh) 一种多轴载荷作用下考虑缺口效应的能量场强法及其应用
CN105653794B (zh) 一种含初始缺陷矩形板结构的时变可靠性设计方法
CN104316388A (zh) 一种对各向异性材料结构件进行疲劳寿命测定的方法
US10989640B2 (en) Method for defining threshold stress curves utilized in fatigue and damage tolerance analysis
US11428614B2 (en) Notch treatment methods for flaw simulation
CN109948216B (zh) 总应变能密度修正的缺口件低周疲劳预测方法
Banerjee et al. Advanced elastic and inelastic three‐dimensional analysis of gas turbine engine structures by BEM
Dey et al. Multiple‐crack identification in a channel section steel beam using a combined response surface methodology and genetic algorithm
Mathiyazhagan et al. Finite element analysis on curved beams of various sections
Montebello Analysis of the stress gradient effect in Fretting-Fatigue through a description based on nonlocal intensity factors
CN114065423B (zh) 快速评估航空发动机风扇叶片颤振的方法
Hoang Computational investigation of variation in wing aerodynamic load under effect of aeroelastic deformations
Liu et al. Numerical limit analysis of cylindrical shells with part-through slots
CN107066727B (zh) 三维空间矢量应力场强法
CN114139276A (zh) 一种盘轴一体式整体叶盘结构疲劳寿命分析方法
CN109242335B (zh) 一种基于概率密度演化理论的概率失效风险评估高效计算方法
Lai et al. Effects of local thickness defects on the buckling of micro-beam
Zhou et al. Non-uniform rational B-spline based free vibration analysis of axially functionally graded tapered Timoshenko curved beams
CN112906255B (zh) 一种裂纹尖端应力强度因子的测量方法
Chakshu et al. Natural frequencies of pre-twisted airfoil blades
Li et al. Stress intensity factors of corner cracks in three dimensional finite plates based on universal weight function method
Li et al. Stress Intensity Factors for Corner Cracks under Arbitrary Stress Fields in the Finite Plates Based on the Point Weight Function Method
Ball The Impact of Bulk Residual Stress on the Qualification of Large Aluminum Forgings

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Li Guo

Inventor after: Ding Shuiting

Inventor after: Zhou Huimin

Inventor after: Liu Junbo

Inventor after: Liu Xiaojing

Inventor before: Ding Shuiting

Inventor before: Li Guo

Inventor before: Zhou Huimin

Inventor before: Liu Junbo

Inventor before: Liu Xiaojing

GR01 Patent grant
GR01 Patent grant