CN110147643B - 车钩钩体剩余寿命确定方法和装置 - Google Patents

车钩钩体剩余寿命确定方法和装置 Download PDF

Info

Publication number
CN110147643B
CN110147643B CN201910506291.4A CN201910506291A CN110147643B CN 110147643 B CN110147643 B CN 110147643B CN 201910506291 A CN201910506291 A CN 201910506291A CN 110147643 B CN110147643 B CN 110147643B
Authority
CN
China
Prior art keywords
crack
determining
node
load spectrum
stress intensity
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
CN201910506291.4A
Other languages
English (en)
Other versions
CN110147643A (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.)
China Shenhua Energy Co Ltd
Shenhua Rail and Freight Wagons Transport Co Ltd
Original Assignee
China Shenhua Energy Co Ltd
Shenhua Rail and Freight Wagons Transport Co Ltd
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 China Shenhua Energy Co Ltd, Shenhua Rail and Freight Wagons Transport Co Ltd filed Critical China Shenhua Energy Co Ltd
Priority to CN201910506291.4A priority Critical patent/CN110147643B/zh
Publication of CN110147643A publication Critical patent/CN110147643A/zh
Application granted granted Critical
Publication of CN110147643B publication Critical patent/CN110147643B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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/06—Power analysis or power optimisation

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)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明实施例提供一种车钩钩体剩余寿命确定方法和装置，属于铁路运输领域。该方法包括：检测车钩钩体的初始裂纹尺寸；确定每级载荷谱的裂纹前沿各节点的应力强度因子范围；确定裂纹前沿各节点的平均裂纹扩展速率；确定单次裂纹扩展所需的循环加载次数；确定所述扩展点以外的节点的扩展量；确定新的裂纹尺寸，并使用所述新的裂纹尺寸作为所述初始裂纹尺寸重复上述步骤，直到裂纹前沿各节点的最大应力强度因子有任一者等于断裂韧性；根据每次重复步骤得到的单次裂纹扩展所需的循环加载次数、载荷谱的频次以及载荷谱对应的服役公里确定所述车钩钩体的剩余寿命。该车钩钩体剩余寿命确定方法和装置可以准确评估车钩钩体产生裂纹后的剩余寿命。

Description

车钩钩体剩余寿命确定方法和装置
技术领域
本发明涉及铁路运输，具体地涉及车钩钩体剩余寿命确定方法和装置。
背景技术
随着人们对结构疲劳寿命的深入研究，损伤容限法逐渐被应用到工程结构的剩余寿命预测中。车钩钩体作为铁路货车关键承载部件，对其剩余寿命的准确把握对优化车钩维修策略，提高结构使用寿命具有较大的工程实用价值。当前对车钩钩体出现裂纹后结构剩余寿命预测的方法较少，缺少一套完整的车钩剩余寿命预测理论和方法。目前因为无法准确了解车钩钩体出现裂纹后的剩余服役寿命，对车钩裂纹的检查相对严格，一旦车钩主体结构出现裂纹，即视为报废，导致了车钩使用寿命较短。
发明内容
本发明实施例的目的是提供一种车钩钩体剩余寿命确定方法和装置，该车钩钩体剩余寿命确定方法和装置可以准确评估车钩钩体产生裂纹后的剩余寿命。
为了实现上述目的，本发明提供一种车钩钩体剩余寿命确定方法，该方法包括：检测所述车钩钩体的初始裂纹尺寸；根据每级载荷谱的载荷以及所述车钩钩体的初始裂纹尺寸，确定每级载荷谱的裂纹前沿各节点的应力强度因子范围；根据所述每级载荷谱的裂纹前沿各节点的应力强度因子范围、所述每级载荷谱的载荷应力比以及所述载荷谱总循环次数确定裂纹前沿各节点的平均裂纹扩展速率；根据扩展点的平均裂纹扩展速率和设定的裂纹扩展点的扩展量确定单次裂纹扩展所需的循环加载次数；根据所述单次裂纹扩展所需的循环加载次数和所述扩展点以外的节点的平均裂纹扩展速率确定所述扩展点以外的节点的扩展量；根据所述扩展点以外的节点的扩展量和初始裂纹尺寸确定新的裂纹尺寸，并使用所述新的裂纹尺寸作为所述初始裂纹尺寸重复上述步骤，直到所述裂纹前沿各节点的最大应力强度因子有任一者等于断裂韧性；根据每次重复步骤得到的单次裂纹扩展所需的循环加载次数、载荷谱的频次以及载荷谱对应的服役公里确定所述车钩钩体的剩余寿命。
优选地，所述扩展点为所述裂纹前沿各节点的最大和最小的应力强度因子的中间值所对应的点。
优选地，每级载荷谱的裂纹前沿各节点的应力强度因子范围通过以下公式得到：
其中，ΔK为所述每级载荷谱的裂纹前沿各节点的应力强度因子范围，Kmax为每级载荷谱的裂纹前沿各节点的最大应力强度因子，Kmin为每级载荷谱的裂纹前沿各节点的最小应力强度因子。
优选地，所述裂纹前沿各节点的平均裂纹扩展速率通过以下公式得到：
其中，为所述裂纹前沿各节点的平均裂纹扩展速率，/>为每级载荷谱对应的裂纹前沿各节点的裂纹扩展速率，n为载荷谱总循环次数。
优选地，所述单次裂纹扩展所需的循环加载次数通过以下公式得到：
其中，ΔN为所述单次裂纹扩展所需的循环加载次数，Δa为设定的裂纹扩展点的扩展量，为所述扩展点的平均裂纹扩展速率。
优选地，所述扩展点以外的节点的扩展量通过以下公式得到：
其中，Δaj为所述扩展点以外的节点的扩展量，为所述扩展点以外的节点的平均裂纹扩展速率，ΔN为所述单次裂纹扩展所需的循环加载次数。
优选地，所述车钩钩体的剩余寿命通过以下公式得到：
其中，Lresidual为所述车钩钩体的剩余寿命，为每次重复步骤得到的单次裂纹扩展所需的循环加载次数之和，ΔL为载荷谱对应的服役公里，n为载荷谱总循环次数。
本发明提供一种车钩钩体剩余寿命确定装置，该装置包括：检测单元以及处理单元，其中，所述检测单元用于检测所述车钩钩体的初始裂纹尺寸；所述处理单元用于：根据每级载荷谱的载荷以及所述车钩钩体的初始裂纹尺寸，确定每级载荷谱的裂纹前沿各节点的应力强度因子范围；根据所述每级载荷谱的裂纹前沿各节点的应力强度因子范围、所述每级载荷谱的载荷应力比以及所述载荷谱总循环次数确定裂纹前沿各节点的平均裂纹扩展速率；根据扩展点的平均裂纹扩展速率和设定的裂纹扩展点的扩展量确定单次裂纹扩展所需的循环加载次数；根据所述单次裂纹扩展所需的循环加载次数和所述扩展点以外的节点的平均裂纹扩展速率确定所述扩展点以外的节点的扩展量；根据所述扩展点以外的节点的扩展量和初始裂纹尺寸确定新的裂纹尺寸，并使用所述新的裂纹尺寸作为所述初始裂纹尺寸重复上述步骤，直到所述裂纹前沿各节点的最大应力强度因子有任一者等于断裂韧性；根据每次重复步骤得到的单次裂纹扩展所需的循环加载次数、载荷谱的频次以及载荷谱对应的服役公里确定所述车钩钩体的剩余寿命。
优选地，所述扩展点为所述裂纹前沿各节点的最大和最小的应力强度因子的中间值所对应的点。
优选地，所述车钩钩体的剩余寿命通过以下公式得到：
其中，Lresidual为所述车钩钩体的剩余寿命，为每次重复步骤得到的单次裂纹扩展所需的循环加载次数之和，ΔL为载荷谱对应的服役公里，n载荷谱总循环次数。
通过上述技术方案，采用本发明提供的车钩钩体剩余寿命确定方法，首先检测所述车钩钩体的初始裂纹尺寸，然后确定每级载荷谱的裂纹前沿各节点的应力强度因子范围，接着确定裂纹前沿各节点的平均裂纹扩展速率，然后确定单次裂纹扩展所需的循环加载次数，再确定扩展点以外的节点的扩展量，最后裂纹前沿各节点的最大应力强度因子有任一者等于断裂韧性时确定车钩钩体的剩余寿命。通过仿真计算的方法模拟车钩钩体裂纹扩展过程，得出结构的裂纹扩展速率方程，即车钩钩体特定部位裂纹尺寸与服役里程的关系，基于裂纹扩展速率方程评估车钩钩体产生裂纹后的剩余寿命。
本发明实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施例，但并不构成对本发明实施例的限制。在附图中：
图1是本发明一实施例提供的车钩钩体剩余寿命确定方法的流程图；
图2是本发明一实施例提供的载荷谱的示意图；
图3是本发明一实施例提供的车钩钩体剩余寿命确定装置的结构示意图。
附图标记说明
1 检测单元 2 处理单元。
具体实施方式
以下结合附图对本发明实施例的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明实施例，并不用于限制本发明实施例。
图1是本发明一实施例提供的车钩钩体剩余寿命确定方法的流程图。如图1所示，该方法包括：
步骤S11，检测所述车钩钩体的初始裂纹尺寸；
步骤S12，根据每级载荷谱的载荷以及所述车钩钩体的初始裂纹尺寸，确定每级载荷谱的裂纹前沿各节点的应力强度因子范围；
步骤S13，根据所述每级载荷谱的裂纹前沿各节点的应力强度因子范围、所述每级载荷谱的载荷应力比以及所述载荷谱总循环次数确定裂纹前沿各节点的平均裂纹扩展速率；
步骤S14，根据扩展点的平均裂纹扩展速率和设定的裂纹扩展点的扩展量确定单次裂纹扩展所需的循环加载次数；
步骤S15，根据所述单次裂纹扩展所需的循环加载次数和所述扩展点以外的节点的平均裂纹扩展速率确定所述扩展点以外的节点的扩展量；
步骤S16，根据所述扩展点以外的节点的扩展量和初始裂纹尺寸确定新的裂纹尺寸，并使用所述新的裂纹尺寸作为所述初始裂纹尺寸重复上述步骤，直到所述裂纹前沿各节点的最大应力强度因子有任一者等于断裂韧性；
步骤S17，根据每次重复步骤得到的单次裂纹扩展所需的循环加载次数、载荷谱的频次以及载荷谱对应的服役公里确定所述车钩钩体的剩余寿命。
本发明实施例提供以下方式实现车钩钩体的剩余寿命的计算：
首先，检查车钩钩体初始裂纹位置和尺寸，确定初始裂纹参数，车钩裂纹大多为表面裂纹，裂纹深度尺寸用aA,0表示，裂纹表面长度尺寸用2×aC,0表示，aC,0为裂纹表面长度尺寸的1/2。
接着，根据车钩钩体的尺寸参数建立车钩的三维几何模型。
接着，取车钩钩体表面材料加工材料试件进行相应的疲劳和断裂实验，得到材料的疲劳极限、S-N曲线、裂纹扩展门槛值、裂纹扩展速率方程、断裂韧性。
接着，根据车钩实测载荷进行载荷谱编制块，用于车钩疲劳寿命仿真计算，如图2的载荷谱示意图所示。(载荷谱编制方法有相关论文介绍，在此不再赘述)。
接着，建立车钩钩体有限元模型，根据车钩钩体载荷谱施加最危险的载荷工况，进行静力学分析，校核车钩构体静强度是否满足列车运行要求。
接着，参考实际情况在车钩钩体有限元模型中引入初始裂纹，由网格自适应划分技术自动离散网格，为准确计算裂纹前沿的应力状态，在裂纹前沿形成对称环形网格区域。内圈单元类型为15节点楔形单元，将特定位置节点移动到1/4节点位置，靠近裂纹前沿应力强度因子计算点，即形成奇异单元；外圈单元类型为20节点六面体单元。
接着，根据M积分法(也叫交互积分法)计算裂纹前沿各个节点在特定载荷下的应力强度因子值。具体说就是先应用有限元计算软件计算出裂纹附近的应力应变场，在根据M积分法计算出裂纹前沿的应力前度因子值。
接着，确定裂纹前沿所有节点最大和最小的应力强度因子值，计算最大和最小应力强度因子值的中间值，并将裂纹前沿应力强度因子值等于或最接近中间值的节点作为扩展点。
接着，针对裂纹前沿各节点的裂纹扩展速率的计算方法相同，以下以扩展点为例，裂纹扩展过程是在载荷谱块的作用下进行的，但有限元模型采用了单载荷形式，就必须采用一种换算方法。扩展点在单载荷下有一个应力强度因子值KF，根据线性换算扩展点在其他载荷作用下的应力强度因子值KFi，计算方程如下：
载荷谱块每一个载荷循环都有一个最大值和一个最小值载荷，也就有两个应力强度因子值(应力强度因子与载荷以及裂纹尺寸相关)，最大与最小应力强度因子的差值就是这个载荷循环下的应力强度因子范围ΔK。计算方程如下：
其中，ΔK为所述每级载荷谱的扩展点(或裂纹前沿各节点)的应力强度因子范围，Kmax为每级载荷谱的扩展点(或裂纹前沿各节点)的最大应力强度因子，Kmin为每级载荷谱的扩展点(或裂纹前沿各节点)的最小应力强度因子。
根据材料裂纹扩展速率方程，每一个应力强度因子范围值都对应一个裂纹扩展速率。裂纹扩展速率方程如下：
其中，n为载荷谱总循环次数。
载荷谱块有一定的循环次数，本实施例使用的是18级载荷谱(但本发明不限于此)，循环次数为1240。也就是说有1240个裂纹扩展速率，我们将所有的裂纹扩展速率进行平均，得到扩展点(或裂纹前沿各节点)的平均裂纹扩展速率。计算方程如下：
其中，为所述扩展点(或裂纹前沿各节点)的平均裂纹扩展速率，/>为每级载荷谱对应的扩展点(或裂纹前沿各节点)的裂纹扩展速率，n为载荷谱总循环次数。
接着，人为设定单次裂纹扩展量Δa，即扩展点的扩展量为Δa，通过扩展点平均裂纹扩展速率方程得到单次扩展所需的循环加载次数。计算方程如下：
其中，ΔN为所述单次裂纹扩展所需的循环加载次数，Δa为设定的裂纹扩展点的扩展量，为所述扩展点的平均裂纹扩展速率。
裂纹前沿除扩展点外的节点扩展量Δaj计算方程如下：
其中，Δaj为所述扩展点以外的节点的扩展量，为所述扩展点以外的节点的平均裂纹扩展速率，ΔN为所述单次裂纹扩展所需的循环加载次数。
接着，扩展方向按照最大张开型裂纹应力强度因子理论确定。理论公式为式如下：
max(KI)
接着，裂纹扩展参量计算完成，进行重新进行网格自适应划分，形成新的裂纹前沿有限元模型。也就确定了新的裂纹深度aA,1和裂纹长度2×aC,1，以及进行的循环次数ΔN1。
最后，反复进行上述的有限元建模和后续步骤，直至裂纹前沿节点的最大应力强度因子值达到材料断裂韧性值。确定此时的aA,m和2×aC,m为裂纹临界允许尺寸，为裂纹扩展过程在的载荷加载次数，/>为载荷谱块的加载次数。m为扩展步数，n1为载荷谱块的频次。/>乘上载荷谱块对应的服役公里数就是车钩出现裂纹的剩余服役寿命，计算式如下：
其中，Lresidual为所述车钩钩体的剩余寿命，为每次重复步骤得到的单次裂纹扩展所需的循环加载次数之和，ΔL为载荷谱对应的服役公里，n为载荷谱总循环次数。
图3是本发明一实施例提供的车钩钩体剩余寿命确定装置的结构示意图。如图3所示，该装置包括：检测单元1以及处理单元2，其中，所述检测单元1用于检测所述车钩钩体的初始裂纹尺寸；所述处理单元2用于：根据每级载荷谱的载荷以及所述车钩钩体的初始裂纹尺寸，确定每级载荷谱的裂纹前沿各节点的应力强度因子范围；根据所述每级载荷谱的裂纹前沿各节点的应力强度因子范围、所述每级载荷谱的载荷应力比以及所述载荷谱总循环次数确定裂纹前沿各节点的平均裂纹扩展速率；根据扩展点的平均裂纹扩展速率和设定的裂纹扩展点的扩展量确定单次裂纹扩展所需的循环加载次数；根据所述单次裂纹扩展所需的循环加载次数和所述扩展点以外的节点的平均裂纹扩展速率确定所述扩展点以外的节点的扩展量；根据所述扩展点以外的节点的扩展量和初始裂纹尺寸确定新的裂纹尺寸，并使用所述新的裂纹尺寸作为所述初始裂纹尺寸重复上述步骤，直到所述裂纹前沿各节点的最大应力强度因子有任一者等于断裂韧性；根据每次重复步骤得到的单次裂纹扩展所需的循环加载次数、载荷谱的频次以及载荷谱对应的服役公里确定所述车钩钩体的剩余寿命。
优选地，所述扩展点为所述裂纹前沿各节点的最大和最小的应力强度因子的中间值所对应的点。
优选地，所述车钩钩体的剩余寿命通过以下公式得到：
其中，Lresidual为所述车钩钩体的剩余寿命，为每次重复步骤得到的单次裂纹扩展所需的循环加载次数之和，ΔL为载荷谱对应的服役公里，n为载荷谱总循环次数。
上文所述的装置的实施例与上文所述的方法的实施例类似，在此不再赘述。
通过上述技术方案，采用本发明提供的车钩钩体剩余寿命确定方法，首先检测所述车钩钩体的初始裂纹尺寸，然后确定每级载荷谱的裂纹前沿各节点的应力强度因子范围，接着确定裂纹前沿各节点的平均裂纹扩展速率，然后确定单次裂纹扩展所需的循环加载次数，再确定扩展点以外的节点的扩展量，最后裂纹前沿各节点的最大应力强度因子有任一者等于断裂韧性时确定车钩钩体的剩余寿命。通过仿真计算的方法模拟车钩钩体裂纹扩展过程，得出结构的裂纹扩展速率方程，即车钩钩体特定部位裂纹尺寸与服役里程的关系，基于裂纹扩展速率方程评估车钩钩体产生裂纹后的剩余寿命。
以上结合附图详细描述了本发明实施例的可选实施方式，但是，本发明实施例并不限于上述实施方式中的具体细节，在本发明实施例的技术构思范围内，可以对本发明实施例的技术方案进行多种简单变型，这些简单变型均属于本发明实施例的保护范围。
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明实施例对各种可能的组合方式不再另行说明。
本领域技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序存储在一个存储介质中，包括若干指令用以使得单片机、芯片或处理器(processor)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
此外，本发明实施例的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明实施例的思想，其同样应当视为本发明实施例所公开的内容。

Claims (6)

1.一种车钩钩体剩余寿命确定方法，其特征在于，该方法包括：
检测所述车钩钩体的初始裂纹尺寸；
根据每级载荷谱的载荷以及所述车钩钩体的初始裂纹尺寸，确定每级载荷谱的裂纹前沿各节点的应力强度因子范围，其中，所述每级载荷谱的裂纹前沿各节点的应力强度因子范围通过以下公式得到：
其中，ΔK为所述每级载荷谱的裂纹前沿各节点的应力强度因子范围，Kmax为每级载荷谱的裂纹前沿各节点的最大应力强度因子，Kmin为每级载荷谱的裂纹前沿各节点的最小应力强度因子，
所述最大应力强度因子Kmax和所述最小应力强度因子Kmin的确定方法包括：
应用有限元计算软件计算出所述车钩钩体的应力应变场，其中，所述有限元计算软件的网格单元划分包括：将内圈单元类型设置为15节点楔形单元，将外圈单元类型设置为20节点六面体单元，以及
根据M积分法以及所述应力应变场计算出所述每级载荷谱的裂纹前沿的应力强度因子值，在所述每级载荷谱的裂纹前沿的应力强度因子值中确定最大应力强度因子Kmax和最小应力强度因子Kmin；
根据所述每级载荷谱的裂纹前沿各节点的应力强度因子范围、所述每级载荷谱的载荷应力比以及所述载荷谱总循环次数确定裂纹前沿各节点的平均裂纹扩展速率，其中，所述裂纹前沿各节点的平均裂纹扩展速率通过以下公式得到：
其中，为所述裂纹前沿各节点的平均裂纹扩展速率，/>为每级载荷谱对应的裂纹前沿各节点的裂纹扩展速率，n为载荷谱总循环次数；
根据扩展点的平均裂纹扩展速率和设定的裂纹扩展点的扩展量确定单次裂纹扩展所需的循环加载次数，其中，所述扩展点为所述裂纹前沿各节点的最大和最小的应力强度因子的中间值所对应的点；
根据所述单次裂纹扩展所需的循环加载次数和所述扩展点以外的节点的平均裂纹扩展速率确定所述扩展点以外的节点的扩展量；
根据所述扩展点以外的节点的扩展量和初始裂纹尺寸确定新的裂纹尺寸，并使用所述新的裂纹尺寸作为所述初始裂纹尺寸重复上述步骤，直到所述裂纹前沿各节点的最大应力强度因子有任一者等于断裂韧性；
根据每次重复步骤得到的单次裂纹扩展所需的循环加载次数、载荷谱的频次以及载荷谱对应的服役公里确定所述车钩钩体的剩余寿命。
2.根据权利要求1所述的车钩钩体剩余寿命确定方法，其特征在于，所述单次裂纹扩展所需的循环加载次数通过以下公式得到：
其中，ΔN为所述单次裂纹扩展所需的循环加载次数，Δa为设定的裂纹扩展点的扩展量，为所述扩展点的平均裂纹扩展速率。
3.根据权利要求1所述的车钩钩体剩余寿命确定方法，其特征在于，所述扩展点以外的节点的扩展量通过以下公式得到：
其中，Δaj为所述扩展点以外的节点的扩展量，为所述扩展点以外的节点的平均裂纹扩展速率，ΔN为所述单次裂纹扩展所需的循环加载次数。
4.根据权利要求1所述的车钩钩体剩余寿命确定方法，其特征在于，所述车钩钩体的剩余寿命通过以下公式得到：
其中，Lresidual为所述车钩钩体的剩余寿命，为每次重复步骤得到的单次裂纹扩展所需的循环加载次数之和，ΔL为载荷谱对应的服役公里，n为载荷谱总循环次数。
5.一种车钩钩体剩余寿命确定装置，其特征在于，该装置包括：
检测单元以及处理单元，其中，
所述检测单元用于检测所述车钩钩体的初始裂纹尺寸；
所述处理单元用于：
根据每级载荷谱的载荷以及所述车钩钩体的初始裂纹尺寸，确定每级载荷谱的裂纹前沿各节点的应力强度因子范围，其中，所述每级载荷谱的裂纹前沿各节点的应力强度因子范围通过以下公式得到：
其中，ΔK为所述每级载荷谱的裂纹前沿各节点的应力强度因子范围，Kmax为每级载荷谱的裂纹前沿各节点的最大应力强度因子，Kmin为每级载荷谱的裂纹前沿各节点的最小应力强度因子，
所述最大应力强度因子Kmax和所述最小应力强度因子Kmin的确定方法包括：
应用有限元计算软件计算出所述车钩钩体的应力应变场，其中，所述有限元计算软件的网格单元划分包括：将内圈单元类型设置为15节点楔形单元，将外圈单元类型设置为20节点六面体单元，以及
根据M积分法以及所述应力应变场计算出所述每级载荷谱的裂纹前沿的应力强度因子值，在所述每级载荷谱的裂纹前沿的应力强度因子值中确定最大应力强度因子Kmax和最小应力强度因子Kmin；
根据所述每级载荷谱的裂纹前沿各节点的应力强度因子范围、所述每级载荷谱的载荷应力比以及所述载荷谱总循环次数确定裂纹前沿各节点的平均裂纹扩展速率，其中，所述裂纹前沿各节点的平均裂纹扩展速率通过以下公式得到：
其中，为所述裂纹前沿各节点的平均裂纹扩展速率，/>为每级载荷谱对应的裂纹前沿各节点的裂纹扩展速率，n为载荷谱总循环次数；
根据扩展点的平均裂纹扩展速率和设定的裂纹扩展点的扩展量确定单次裂纹扩展所需的循环加载次数，其中，所述扩展点为所述裂纹前沿各节点的最大和最小的应力强度因子的中间值所对应的点；
根据所述单次裂纹扩展所需的循环加载次数和所述扩展点以外的节点的平均裂纹扩展速率确定所述扩展点以外的节点的扩展量；
根据所述扩展点以外的节点的扩展量和初始裂纹尺寸确定新的裂纹尺寸，并使用所述新的裂纹尺寸作为所述初始裂纹尺寸重复上述步骤，直到所述裂纹前沿各节点的最大应力强度因子有任一者等于断裂韧性；
根据每次重复步骤得到的单次裂纹扩展所需的循环加载次数、载荷谱的频次以及载荷谱对应的服役公里确定所述车钩钩体的剩余寿命。
6.根据权利要求5所述的车钩钩体剩余寿命确定装置，其特征在于，所述车钩钩体的剩余寿命通过以下公式得到：
其中，Lresidual为所述车钩钩体的剩余寿命，为每次重复步骤得到的单次裂纹扩展所需的循环加载次数之和，ΔL为载荷谱对应的服役公里，n载荷谱总循环次数。
CN201910506291.4A 2019-06-12 2019-06-12 车钩钩体剩余寿命确定方法和装置 Active CN110147643B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910506291.4A CN110147643B (zh) 2019-06-12 2019-06-12 车钩钩体剩余寿命确定方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910506291.4A CN110147643B (zh) 2019-06-12 2019-06-12 车钩钩体剩余寿命确定方法和装置

Publications (2)

Publication Number Publication Date
CN110147643A CN110147643A (zh) 2019-08-20
CN110147643B true CN110147643B (zh) 2023-11-14

Family

ID=67591173

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910506291.4A Active CN110147643B (zh) 2019-06-12 2019-06-12 车钩钩体剩余寿命确定方法和装置

Country Status (1)

Country Link
CN (1) CN110147643B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111027208B (zh) * 2019-12-10 2022-02-01 中国航空综合技术研究所 在役机械装备关键结构元件疲劳载荷作用下的定延寿方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221473A (zh) * 2010-04-14 2011-10-19 广州市特种机电设备检测研究院 一种起重机主金属结构剩余疲劳寿命的估算方法
JP2013050377A (ja) * 2011-08-31 2013-03-14 National Institute Of Advanced Industrial & Technology 寿命・余寿命評価機能付き非破壊検査装置
CN106290559A (zh) * 2016-07-30 2017-01-04 董超超 发动机剩余寿命预测系统
JP2017009330A (ja) * 2015-06-18 2017-01-12 株式会社東芝 き裂進展解析装置およびき裂進展解析方法
WO2017220923A1 (fr) * 2016-06-20 2017-12-28 Safran Aircraft Engines Procédé d'estimation du facteur d'intensité des contraintes et procédé de calcul de durée de vie associé

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221473A (zh) * 2010-04-14 2011-10-19 广州市特种机电设备检测研究院 一种起重机主金属结构剩余疲劳寿命的估算方法
JP2013050377A (ja) * 2011-08-31 2013-03-14 National Institute Of Advanced Industrial & Technology 寿命・余寿命評価機能付き非破壊検査装置
JP2017009330A (ja) * 2015-06-18 2017-01-12 株式会社東芝 き裂進展解析装置およびき裂進展解析方法
WO2017220923A1 (fr) * 2016-06-20 2017-12-28 Safran Aircraft Engines Procédé d'estimation du facteur d'intensité des contraintes et procédé de calcul de durée de vie associé
CN106290559A (zh) * 2016-07-30 2017-01-04 董超超 发动机剩余寿命预测系统

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
三维结构中片状疲劳裂纹扩展行为的数值模拟方法;李有堂等;《兰州理工大学学报》;20160425;第42卷(第2期);第2.3节 *
含广布腐蚀坑结构寿命评估方法;董妍等;《力学与实践》;20110608(第03期);全文 *
基于子模型的车钩裂纹扩展特性研究;薛海等;《机械强度》;20181010(第05期);全文 *
破损舰体结构疲劳裂纹扩展及剩余寿命研究;刘燕红等;《船舶力学》;20101015(第10期);全文 *
重载列车E级钢钩舌疲劳裂纹扩展条件及寿命预测;李晓慧等;《北京交通大学学报》;20060815(第04期);全文 *

Also Published As

Publication number Publication date
CN110147643A (zh) 2019-08-20

Similar Documents

Publication Publication Date Title
US20210027001A1 (en) Method for Predicting Sealing Reliability of Soft Packing Lithium Ion Battery
CN105718633B (zh) 一种底盘件的载荷分析方法
CN105117536B (zh) 一种rpv含裂纹类缺陷的简化弹塑性断裂力学分析方法
CN104820781B (zh) 考虑温循载荷顺序加载影响的bga焊点热疲劳寿命预测方法
CN103616179A (zh) 一种基于缺陷建模的传动齿轮疲劳寿命评估方法
CN107545110B (zh) 一种动态应力加速寿命试验剖面编制方法
CN103712865A (zh) 高温部件的蠕变损伤评价方法以及损伤评价系统
CN103424654A (zh) 一种敏感设备电压暂降敏感度的评估方法
CN106354898A (zh) 一种基于总应变能量密度的焊缝疲劳寿命计算方法
CN105260574A (zh) 一种基于临界面法疲劳破坏准则的高周多轴疲劳寿命预测方法
KR101817343B1 (ko) 수정된 h-m 모델에 기반한 파괴곡선 산출 방법
CN110263401B (zh) 钢箱梁顶板-纵肋焊接细节残余应力松弛效应的评估方法
CN110147643B (zh) 车钩钩体剩余寿命确定方法和装置
CN111881564B (zh) 一种关于机械结构变幅疲劳寿命预测方法
CN110987676A (zh) 随机多轴载荷下考虑裂纹闭合效应的全寿命预测方法
CN111797749A (zh) 电动卡车控制器支架结构振动激励载荷计算方法及系统
CN113609584B (zh) 重复使用火箭发动机推力室可靠性优化设计方法
CN113607580B (zh) 一种金属构件疲劳试验方法及剩余寿命预测方法
CN110955933B (zh) 一种基于响应面法的机械结构模糊疲劳可靠度计算方法
CN111506977A (zh) 一种动力电池建模方法
CN114282423B (zh) 一种岸桥结构的疲劳监测方法和寿命预测方法
CN108334666A (zh) 一种车轴抗疲劳断裂分析方法
Hectors et al. Fracture mechanics and hot spot stress‐based fatigue life calculation: case study for a crane runway girder
CN110736595A (zh) 轨道交通舱室内电子设备寿命评估方法
Gao et al. A novel probabilistic fatigue life prediction method for welded structures based on gPC

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