CN106886663A - 齿轮弯曲疲劳寿命预测方法及装置 - Google Patents
齿轮弯曲疲劳寿命预测方法及装置 Download PDFInfo
- Publication number
- CN106886663A CN106886663A CN201710197972.8A CN201710197972A CN106886663A CN 106886663 A CN106886663 A CN 106886663A CN 201710197972 A CN201710197972 A CN 201710197972A CN 106886663 A CN106886663 A CN 106886663A
- Authority
- CN
- China
- Prior art keywords
- crack
- fatigue
- model
- gear
- 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.)
- Granted
Links
- 238000005452 bending Methods 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000000977 initiatory effect Effects 0.000 claims abstract description 61
- 239000000463 material Substances 0.000 claims abstract description 20
- 238000012937 correction Methods 0.000 claims abstract description 16
- 230000007246 mechanism Effects 0.000 claims abstract description 10
- 230000003746 surface roughness Effects 0.000 claims abstract description 9
- 238000006073 displacement reaction Methods 0.000 claims description 15
- 238000009826 distribution Methods 0.000 claims description 13
- 208000037656 Respiratory Sounds Diseases 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 6
- 150000001875 compounds Chemical class 0.000 claims 1
- 230000006870 function Effects 0.000 description 17
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 7
- 238000003754 machining Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 229910000851 Alloy steel Inorganic materials 0.000 description 1
- 230000002929 anti-fatigue Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000009661 fatigue test Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
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/04—Ageing analysis or optimisation against ageing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种齿轮弯曲疲劳寿命预测方法及装置，该齿轮弯曲疲劳寿命预测方法包括：基于表面粗糙度建立疲劳极限的修正模型，并根据所述修正模型对材料的疲劳极限进行修正，得到修正后疲劳极限；根据门槛应力强度因子幅度及所述修正后疲劳极限确定门槛裂纹长度；创建疲劳裂纹尺寸萌生模型，并基于所述疲劳裂纹尺寸萌生模型预测齿轮疲劳裂纹萌生寿命；基于线弹性断裂力学准则预测齿轮疲劳裂纹扩展寿命；根据所述齿轮疲劳裂纹萌生寿命和齿轮疲劳裂纹扩展寿命，建立齿轮弯曲疲劳计算模型，计算齿轮弯曲疲劳寿命。
Description
技术领域
本发明涉及技术齿轮弯曲疲劳寿命预测技术领域，特别涉及一种齿轮弯曲疲劳寿命预测方法及装置。
背景技术
齿轮是传动系统的关键部件，轮齿弯曲疲劳造成的齿轮折断是齿轮最普遍的一种失效模式。明确齿轮真实应力分布，预测齿轮弯曲载荷下工作寿命，已成为齿轮抗疲劳设计的重要依据。
传统预测齿轮弯曲疲劳寿命的方法是依据大量齿轮弯曲疲劳试验获得齿轮S-N曲线，在其基础上进行强度计算设计，进而预测齿轮的弯曲疲劳寿命。然而这种传统方法并未考虑表面加工状况、齿轮结构的几何特征、应力梯度和平均应力等因素的影响。很难准确地预测齿轮弯曲疲劳寿命，不能揭示齿轮弯曲疲劳失效的机理，此外这种预测方法主要建立在大量试验的基础上，成本较高且周期较长。
发明内容
本发明提供一种齿轮弯曲疲劳寿命预测方法及装置，以准确地预测齿轮弯曲疲劳寿命的方法，减小对齿轮材料、结构尺寸、工艺参数、试验量等因素的依赖性。
为了实现上述目的，本发明实施例提供了一种齿轮弯曲疲劳寿命预测方法，该齿轮弯曲疲劳寿命预测方法包括：
基于表面粗糙度建立疲劳极限的修正模型，并根据所述修正模型对材料的疲劳极限进行修正，得到修正后疲劳极限；
根据门槛应力强度因子幅度及所述修正后疲劳极限确定门槛裂纹长度；
创建疲劳裂纹尺寸萌生模型，并基于所述疲劳裂纹尺寸萌生模型预测齿轮疲劳裂纹萌生寿命；
基于线弹性断裂力学准则预测齿轮疲劳裂纹扩展寿命；
根据所述齿轮疲劳裂纹萌生寿命和齿轮疲劳裂纹扩展寿命，建立齿轮弯曲疲劳计算模型，计算齿轮弯曲疲劳寿命。
一实施例中，该齿轮弯曲疲劳寿命预测方法还包括：
步骤1：根据包含模数、齿数、压力角的齿轮基本参数绘制齿轮齿根二维几何模型；
步骤2：基于所述齿轮齿根二维几何模型，划分网格，施加边界约束，确定承载工况，建立齿轮齿根二维有限元模型；
步骤3：确定平面应变下裂纹尖端处张开型应力强度因子及滑开型应力强度因子与节点位移的拟合函数关系；
步骤4：根据所述拟合函数关系建立平面应变下裂纹尖端复合应力强度因子方程；
步骤5：基于所述齿轮齿根二维有限元模型建立裂纹扩展角计算模型和裂纹扩展增量计算模型，并预测裂纹扩展路径。
一实施例中，所述创建疲劳裂纹尺寸萌生模型，并基于所述疲劳裂纹尺寸萌生模型预测齿轮疲劳裂纹萌生寿命，包括：
根据应力梯度下的裂纹萌生模型，建立齿轮局部应力分布关系；
基于所述齿轮局部应力分布关系，确定作用在裂纹上的平均应力范围；
根据所述平均应力范围建立所述疲劳裂纹尺寸萌生模型；
根据所述疲劳裂纹尺寸萌生模型、修正后疲劳极限及门槛裂纹长度建立裂纹萌生寿命预测模型；
基于所述裂纹萌生寿命预测模型计算齿轮疲劳裂纹萌生寿命。
一实施例中，所述基于线弹性断裂力学准则预测齿轮疲劳裂纹扩展寿命，包括：
根据应力强度因子范围计算模型计算应力强度因子范围；
根据门槛应力强度因子及断裂韧度求解裂纹扩展应力强度因子范围；
基于平均应力对长裂纹扩展的影响，建立裂纹扩展速率修正Paris公式；
基于修正后的所述Paris公式、门槛裂纹长度、临界裂纹尺寸建立齿轮疲劳裂纹扩展寿命模型；
根据所述齿轮疲劳裂纹扩展寿命模型计算所述齿轮疲劳裂纹扩展寿命。
一实施例中，所述基于所述齿轮齿根二维有限元模型建立裂纹扩展角计算模型和裂纹扩展增量计算模型，并预测裂纹扩展路径，包括：
基于最大切向应力准则，建立裂纹扩展角计算模型；
根据所述应力强度因子范围及修正后的所述Paris公式建立裂纹扩展增量计算模型；
重复执行上述步骤3至步骤5，直至应力强度因子达到临界应力强度因子；
基于应力强度因子达到临界应力强度因子过程中所得到的裂纹扩展角及裂纹扩展增量预测试件发生失效时的裂纹扩展路径。
为了实现上述目的，本发明实施例还提供了一种齿轮弯曲疲劳寿命预测装置，该齿轮弯曲疲劳寿命预测装置包括：
疲劳极限修正单元，用于基于表面粗糙度建立疲劳极限的修正模型，并根据所述修正模型对材料的疲劳极限进行修正，得到修正后疲劳极限；
门槛裂纹长度确定单元，用于根据门槛应力强度因子幅度及所述修正后疲劳极限确定门槛裂纹长度；
裂纹萌生寿命预测单元，用于创建疲劳裂纹尺寸萌生模型，并基于所述疲劳裂纹尺寸萌生模型预测齿轮疲劳裂纹萌生寿命；
裂纹扩展寿命预测单元，用于基于线弹性断裂力学准则预测齿轮疲劳裂纹扩展寿命；
齿轮弯曲疲劳寿命计算单元，用于根据所述齿轮疲劳裂纹萌生寿命和齿轮疲劳裂纹扩展寿命，建立齿轮弯曲疲劳计算模型，计算齿轮弯曲疲劳寿命。
一实施例中，该齿轮弯曲疲劳寿命预测装置还包括：
几何模型绘制单元，用于根据包含模数、齿数、压力角的齿轮基本参数绘制齿轮齿根二维几何模型；
有限元模型创建单元，用于基于所述齿轮齿根二维几何模型，划分网格，施加边界约束，确定承载工况，建立齿轮齿根二维有限元模型；
拟合函数关系确定单元，用于确定平面应变下裂纹尖端处张开型应力强度因子及滑开型应力强度因子与节点位移的拟合函数关系；
应力强度方程建立单元，用于根据所述拟合函数关系建立平面应变下裂纹尖端复合应力强度因子方程；
裂纹扩展路径预测单元，用于基于所述齿轮齿根二维有限元模型建立裂纹扩展角计算模型和裂纹扩展增量计算模型，并预测裂纹扩展路径。
一实施例中，所述裂纹萌生寿命预测单元包括：
应力分布建立模块，用于根据应力梯度下的裂纹萌生模型，建立齿轮局部应力分布关系；
平均应力范围确定模块，用于基于所述齿轮局部应力分布关系，确定作用在裂纹上的平均应力范围；
裂纹尺寸萌生模型创建模块，用于根据所述平均应力范围建立所述疲劳裂纹尺寸萌生模型；
裂纹萌生寿命预测模型创建模块，用于根据所述疲劳裂纹尺寸萌生模型、修正后疲劳极限及门槛裂纹长度建立裂纹萌生寿命预测模型；
裂纹萌生寿命预测模块，用于基于所述裂纹萌生寿命预测模型计算齿轮疲劳裂纹萌生寿命。
一实施例中，所述裂纹扩展寿命预测单元包括：
应力强度因子范围计算模块，用于根据应力强度因子范围计算模型计算应力强度因子范围；
裂纹扩展应力强度因子范围计算模块，用于根据门槛应力强度因子及断裂韧度求解裂纹扩展应力强度因子范围；
公式修正模块，用于基于平均应力对长裂纹扩展的影响，建立裂纹扩展速率修正Paris公式；
齿轮疲劳裂纹扩展寿命模型创建模块，用于基于修正后的所述Paris公式、门槛裂纹长度、临界裂纹尺寸建立齿轮疲劳裂纹扩展寿命模型；
齿轮疲劳裂纹扩展寿命计算模块，用于根据所述齿轮疲劳裂纹扩展寿命模型计算所述齿轮疲劳裂纹扩展寿命。
一实施例中，所述裂纹扩展路径预测单元包括：
裂纹扩展角计算模块，用于基于最大切向应力准则，建立裂纹扩展角计算模型；
裂纹扩展增量计算模块，用于根据所述应力强度因子范围及修正后的所述Paris公式建立裂纹扩展增量计算模型；
裂纹扩展路径预测模块，用于基于应力强度因子达到临界应力强度因子过程中所得到的裂纹扩展角及裂纹扩展增量预测试件发生失效时的裂纹扩展路径。
本发明基于裂纹萌生和扩展机理，虑及表面加工状况、结构几何特征，应力梯度以及平均应力的影响，建立了裂纹萌生寿命和扩展寿命计算模型，简化齿轮工作寿命预测流程，可以方便快捷并且准确地预测齿轮弯曲疲劳寿命的方法，大大减少了试验成本，并且减小了对齿轮材料、结构尺寸、工艺参数、试验量等因素的依赖性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例的齿轮弯曲疲劳寿命预测方法流程图；
图2为本发明实施例预测齿轮疲劳裂纹萌生寿命的方法流程图；
图3为本发明实施例预测齿轮疲劳裂纹扩展寿命的方法流程图；
图4为本发明实施例的裂纹扩展路径预测方法流程图；
图5为本实施例的齿轮弯曲疲劳寿命预测装置的结构框图一；
图6为本实施例的裂纹萌生寿命预测单元的结构框图；
图7为本实施例的裂纹扩展寿命预测单元的结构框图；
图8为本实施例的齿轮弯曲疲劳寿命预测装置的结构框图二；
图9本实施例的裂纹扩展路径预测单元单元的结构框图；
图10为本发明实施例的齿轮齿根二维有限元模型示意图；
图11为本发明实施例的表面粗糙度Ra、抗拉强度Rm与表面光洁度修正系数ks函数关系图；
图12为本发明实施例的裂纹尖端奇异节点位移模型；
图13为本发明实施例的平面应变下裂纹尖端处张开型和滑开型应力强度因子与节点位移的拟合函数关系图；
图14为本发明实施例的平面应变下裂纹尖端处复合应力强度因子与节点位移的拟合函数关系图；
图15为本发明实施例的预测裂纹扩展路径示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明实施例的齿轮弯曲疲劳寿命预测方法流程图，如图1所示，该齿轮弯曲疲劳寿命预测方法包括：
S101：基于表面粗糙度建立疲劳极限的修正模型，并根据所述修正模型对材料的疲劳极限进行修正，得到修正后疲劳极限；
S102：根据门槛应力强度因子幅度及所述修正后疲劳极限确定门槛裂纹长度；
S103：创建疲劳裂纹尺寸萌生模型，并基于所述疲劳裂纹尺寸萌生模型预测齿轮疲劳裂纹萌生寿命；
S104：基于线弹性断裂力学准则预测齿轮疲劳裂纹扩展寿命；
S105：根据所述齿轮疲劳裂纹萌生寿命和齿轮疲劳裂纹扩展寿命，建立齿轮弯曲疲劳计算模型，计算齿轮弯曲疲劳寿命。
由图1所示的流程可知，本实施例首先对疲劳极限进行修正，然后根据门槛应力强度因子幅度及修正后疲劳极限确定门槛裂纹长度，然后预测齿轮疲劳裂纹萌生寿命及齿轮疲劳裂纹扩展寿命，最后基于预测对齿轮疲劳裂纹萌生寿命及齿轮疲劳裂纹扩展寿命，计算齿轮弯曲疲劳寿命。利用该方法，可以方便快捷并且准确地预测齿轮弯曲疲劳寿命的方法，大大减少了试验成本。
S101中，疲劳极限修正模型如下：
σfr＝ks·σf (1)
公式(1)中，σfr是修正后疲劳极限，σf是疲劳极限，ks是表面系数。
根据公式(1)所示的疲劳极限修正模型对材料的疲劳极限进行修正，即可得到修正后疲劳极限σfr。
一实施例中，疲劳极限σf＝650MPa，抗拉极限Rm＝1080MPa，粗糙度Ra＝6.4，基于表面粗糙度Ra、抗拉强度Rm与表面系数ks函数关系，如图11，得到此实施例的表面光洁度参数ks＝0.65，据此可由(1)式得到修正后疲劳强度σfr＝423MPa。
S102具体实施时，需要根据门槛应力强度因子范围ΔKth及修正后疲劳极限σfr建立门槛裂纹长度ath的求解方程：
根据门槛应力强度因子幅度ΔKth＝237MPa√mm和S101得到的真实疲劳极限σfr＝423MPa,通过(2)式计算出门槛裂纹长度ath＝0.1mm。
S103具体实施时，如图2所示，包括如下步骤：
S201：根据应力梯度下的裂纹萌生模型，建立齿轮局部应力分布关系：
公式(3)中，σ是齿根局部应力，ρ是齿根半径，kt是缺口弹性应力集中因子，Δσ是应力范围，x是裂纹尖端与齿面距离。
S202：基于齿轮局部应力分布关系，确定作用在裂纹上的平均应力范围。平均应力范围的计算公式如下：
公式(4)中，a为裂纹长度。
积分可得：
根据二项式定理：
(5)式进一步简化为：
S203：根据所述平均应力范围建立所述疲劳裂纹尺寸萌生模型：
公式(8)中，α是萌生指数，M是Taylor因数，μ剪切模量，υ是泊松比，h是滑移带宽度，d是材料颗粒尺寸，常数λ一般取0.005。
S204：根据所述疲劳裂纹尺寸萌生模型、修正后疲劳极限及门槛裂纹长度建立裂纹萌生寿命预测模型。
将(6)式代入(7)式，基于修正后疲劳极限σfr和上述确定的门槛裂纹长度ath，建立裂纹萌生寿命预测模型：
S205：基于所述裂纹萌生寿命预测模型计算齿轮疲劳裂纹萌生寿命。
根据下述公式确定裂纹萌生寿命Ni：
一实施例中，选取参数ρ＝0.25mm，α＝0.5，μ＝7.76×104MPa、h＝1.5×10-3μm、λ＝0.005、d＝1μm、kt＝5、Δσ＝609MPa。
基于S101中σfr＝423MPa和S102中ath＝0.1mm，根据(10)式可以求得裂纹萌生寿命Ni＝2.327×106。
S104具体实施时，如图3所示，包括如下步骤：
S301：根据应力强度因子范围计算模型计算应力强度因子范围；
应力强度因子范围计算模型如下：
ΔK＝Kmax-Kmin (11)
公式(11)中，ΔK为应力强度因子范围，Kmax为最大应力强度因子，Kmin为最小应力强度因子。
S302：根据门槛应力强度因子及断裂韧度求解裂纹扩展应力强度因子范围。
具体地，首先根据材料属性的门槛应力强度因子Kth和断裂韧度Kc，求解裂纹扩展应力强度因子范围ΔKp：
ΔKp＝Kc-Kth (12)
然后求解应力比R：
公式(13)中，σmin为最小应力，σmax为最大应力，σm为平均应力，σa为应力幅值。
基于Paris公式：
其中，C和m是材料参数。
S303：基于平均应力对长裂纹扩展的影响，建立裂纹扩展速率修正Paris公式。
考虑到平均应力对长裂纹扩展的影响，可以建立裂纹扩展速率修正Paris公式：
S304：基于修正后的所述Paris公式、门槛裂纹长度、临界裂纹尺寸建立齿轮疲劳裂纹扩展寿命模型。
根据上述求得的门槛裂纹长度ath以及材料的临界裂纹长度ac，建立裂纹扩展阶段寿命Np计算模型：
S305：根据所述齿轮疲劳裂纹扩展寿命模型计算所述齿轮疲劳裂纹扩展寿命。
具体实施为，门槛应力强度因子Kth≈269MPa√mm，断裂韧度Kc＝2620MPa√mm。可由(12)式求解应力强度因子范围ΔKp＝2351MPa√mm。
根据上述(13)式可以求解应力比R＝0。
选取材料参数C＝3.31×10-17mm/cycl/(MPa√mm)m，m＝4.16，临界裂纹长度ac＝8.6mm以及S102中ath＝0.1mm。根据上述(16)式求解裂纹扩展寿命Np＝4.372×105。
在图1所示方法基础上，本发明还可以预测裂纹扩展路径，图4为本实施例的裂纹扩展路径预测流程图，如图4所示，裂纹扩展路径预测方法包括：
S401：根据包含模数、齿数、压力角的齿轮基本参数绘制齿轮齿根二维几何模型。
一实施例中，模数mn＝4.5mm、齿数z＝39、压力角αn＝24°。网格单元为自由三角形网格，单元采用四节点双线性平面应力单元(GPS4R)。在此实施例中，边界约束为齿根下边界和两侧边界固定、加载载荷为F＝1000N/mm，如图10所示。
S402：基于所述齿轮齿根二维几何模型，划分网格，施加边界约束，确定承载工况，建立齿轮齿根二维有限元模型。
本实施例中，齿轮材料为高强度合金钢42CrMo4,表面处理为完全硬化热处理，其材料参数包括：弹性模量E＝2.1×105MPa、泊松比υ＝0.3。
S403：确定平面应变下裂纹尖端处张开型应力强度因子及滑开型应力强度因子与节点位移的拟合函数关系。
基于S402建立的二维有限元模型，根据裂纹尖端1/4奇异节点位移模型，分别求解张开型应力强度因子KⅠ和滑开型应力强度因子KⅡ，确定平面应变下裂纹尖端处应力强度因子KⅠ、KⅡ与节点位移的拟合函数关系：
其中，G是材料的剪切模量，υ是泊松比，L是有限元网格长度，v、u是b、c、d、e四个节点分别是法向和切向上的位移。
S404：根据所述拟合函数关系建立平面应变下裂纹尖端复合应力强度因子方程。
建立平面应变下裂纹尖端复合应力强度因子方程：
一实施例中，建立的裂纹尖端1/4奇异节点位移模型，如图12所示。基于(17)式求解张开型和滑开型应力强度因子，分别建立张开型和滑开型应力强度因子与节点位移的拟合函数关系，如图13所示。
基于(18)式建立的平面应变下裂纹尖端处应力强度因子与节点位移的拟合函数关系如图14所示。
S405：基于所述齿轮齿根二维有限元模型建立裂纹扩展角计算模型和裂纹扩展增量计算模型，并预测裂纹扩展路径。
基于最大切向应力准则，建立裂纹扩展角计算模型：
根据所述应力强度因子范围ΔK及修正后的所述Paris公式建立裂纹扩展增量计算模型：
其中，ΔN是裂纹扩展增量Δa所需的循环次数，ΔK为裂纹扩展增量Δa对应的应力强度因子范围。
重复执行S403至S405，直至应力强度因子K达到临界应力强度因子Kc，此时发生失效，基于应力强度因子达到临界应力强度因子过程中所得到的裂纹扩展角及裂纹扩展增量预测试件发生失效时的裂纹扩展路径。从而预测裂纹扩展路径，如图15所示。
S105中，根据公式(10)及公式(16)，就可以计算出齿轮弯曲疲劳寿命：
基于上述得到的裂纹萌生寿命Ni＝2.327×106和裂纹扩展寿命Np＝4.372×105，根据上述(21)式，可以得到齿轮弯曲疲劳寿命N＝2.7642×106。
本发明基于裂纹萌生和扩展机理，虑及表面加工状况、结构几何特征，应力梯度以及平均应力的影响，建立了裂纹萌生寿命和扩展寿命计算模型，简化齿轮工作寿命预测流程，可以方便快捷并且准确地预测齿轮弯曲疲劳寿命的方法，大大减少了试验成本，并且减小了对齿轮材料、结构尺寸、工艺参数、试验量等因素的依赖性。
基于与上述齿轮弯曲疲劳寿命预测方法相同的发明构思，本申请提供一种齿轮弯曲疲劳寿命预测装置，如下面实施例所述。由于该齿轮弯曲疲劳寿命预测装置解决问题的原理与齿轮弯曲疲劳寿命预测方法相似，因此该齿轮弯曲疲劳寿命预测装置终端的实施可以参见齿轮弯曲疲劳寿命预测方法的实施，重复之处不再赘述。
图5为本实施例的齿轮弯曲疲劳寿命预测装置的结构框图，该齿轮弯曲疲劳寿命预测装置包括：疲劳极限修正单元501，门槛裂纹长度确定单元502，裂纹萌生寿命预测单元503，裂纹扩展寿命预测单元504及齿轮弯曲疲劳寿命计算单元505。
疲劳极限修正单元501基于表面粗糙度建立疲劳极限的修正模型，并根据所述修正模型对材料的疲劳极限进行修正，得到修正后疲劳极限；
门槛裂纹长度确定单元502，用于根据门槛应力强度因子幅度及所述修正后疲劳极限确定门槛裂纹长度；
裂纹萌生寿命预测单元503，用于创建疲劳裂纹尺寸萌生模型，并基于所述疲劳裂纹尺寸萌生模型预测齿轮疲劳裂纹萌生寿命；
裂纹扩展寿命预测单元504，用于基于线弹性断裂力学准则预测齿轮疲劳裂纹扩展寿命；
齿轮弯曲疲劳寿命计算单元505，用于根据所述齿轮疲劳裂纹萌生寿命和齿轮疲劳裂纹扩展寿命，建立齿轮弯曲疲劳计算模型，计算齿轮弯曲疲劳寿命。
一实施例中，如图6所示，裂纹萌生寿命预测单元503包括：
应力分布建立模块601，用于根据应力梯度下的裂纹萌生模型，建立齿轮局部应力分布关系；
平均应力范围确定模块602，503用于基于所述齿轮局部应力分布关系，确定作用在裂纹上的平均应力范围；
裂纹尺寸萌生模型创建模块603，用于根据所述平均应力范围建立所述疲劳裂纹尺寸萌生模型；
裂纹萌生寿命预测模型创建模块604，用于根据所述疲劳裂纹尺寸萌生模型、修正后疲劳极限及门槛裂纹长度建立裂纹萌生寿命预测模型；
裂纹萌生寿命预测模块605，用于基于所述裂纹萌生寿命预测模型计算齿轮疲劳裂纹萌生寿命。
一实施例中，如图7所示，裂纹扩展寿命预测单元504包括：
应力强度因子范围计算模块701，用于根据应力强度因子范围计算模型计算应力强度因子范围；
裂纹扩展应力强度因子范围计算模块702，用于根据门槛应力强度因子及断裂韧度求解裂纹扩展应力强度因子范围；
公式修正模块703，用于基于平均应力对长裂纹扩展的影响，建立裂纹扩展速率修正Paris公式；
齿轮疲劳裂纹扩展寿命模型创建模块704，用于基于修正后的所述Paris公式、门槛裂纹长度、临界裂纹尺寸建立齿轮疲劳裂纹扩展寿命模型；
齿轮疲劳裂纹扩展寿命计算模块705，用于根据所述齿轮疲劳裂纹扩展寿命模型计算所述齿轮疲劳裂纹扩展寿命。
一实施例中，如图8所示，该齿轮弯曲疲劳寿命预测装置还包括：
几何模型绘制单元801，用于根据包含模数、齿数、压力角的齿轮基本参数绘制齿轮齿根二维几何模型；
有限元模型创建单元802，用于基于所述齿轮齿根二维几何模型，划分网格，施加边界约束，确定承载工况，建立齿轮齿根二维有限元模型。
拟合函数关系确定单元803，用于确定平面应变下裂纹尖端处张开型应力强度因子及滑开型应力强度因子与节点位移的拟合函数关系；
应力强度方程建立单元804，用于根据所述拟合函数关系建立平面应变下裂纹尖端复合应力强度因子方程；
裂纹扩展路径预测单元805，用于基于所述齿轮齿根二维有限元模型建立裂纹扩展角计算模型和裂纹扩展增量计算模型，并预测裂纹扩展路径。
一实施例中，如图9所示，裂纹扩展路径预测单元805包括：
裂纹扩展角计算模块901，用于基于最大切向应力准则，建立裂纹扩展角计算模型；
裂纹扩展增量计算模块902，用于根据所述应力强度因子范围及修正后的所述Paris公式建立裂纹扩展增量计算模型；
裂纹扩展路径预测模块903，用于基于应力强度因子达到临界应力强度因子过程中所得到的裂纹扩展角及裂纹扩展增量预测试件发生失效时的裂纹扩展路径。
本发明基于裂纹萌生和扩展机理，虑及表面加工状况、结构几何特征，应力梯度以及平均应力的影响，建立了裂纹萌生寿命和扩展寿命计算模型，简化齿轮工作寿命预测流程，可以方便快捷并且准确地预测齿轮弯曲疲劳寿命的方法，大大减少了试验成本，并且减小了对齿轮材料、结构尺寸、工艺参数、试验量等因素的依赖性。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (9)
1.一种齿轮弯曲疲劳寿命预测方法，其特征在于，包括：
基于表面粗糙度建立疲劳极限的修正模型，并根据所述修正模型对材料的疲劳极限进行修正，得到修正后疲劳极限；
根据门槛应力强度因子幅度及所述修正后疲劳极限确定门槛裂纹长度；
创建疲劳裂纹尺寸萌生模型，并基于所述疲劳裂纹尺寸萌生模型预测齿轮疲劳裂纹萌生寿命；
基于线弹性断裂力学准则预测齿轮疲劳裂纹扩展寿命；
根据所述齿轮疲劳裂纹萌生寿命和齿轮疲劳裂纹扩展寿命，建立齿轮弯曲疲劳计算模型，计算齿轮弯曲疲劳寿命。
2.根据权利要求1所述的齿轮弯曲疲劳寿命预测方法，其特征在于，还包括：
步骤1：根据包含模数、齿数、压力角的齿轮基本参数绘制齿轮齿根二维几何模型；
步骤2：基于所述齿轮齿根二维几何模型，划分网格，施加边界约束，确定承载工况，建立齿轮齿根二维有限元模型；
步骤3：确定平面应变下裂纹尖端处张开型应力强度因子及滑开型应力强度因子与节点位移的拟合函数关系；
步骤4：根据所述拟合函数关系建立平面应变下裂纹尖端复合应力强度因子方程；
步骤5：基于所述齿轮齿根二维有限元模型建立裂纹扩展角计算模型和裂纹扩展增量计算模型，并预测裂纹扩展路径。
3.根据权利要求2所述的齿轮弯曲疲劳寿命预测方法，其特征在于，所述基于线弹性断裂力学准则预测齿轮疲劳裂纹扩展寿命，包括：
根据应力强度因子范围计算模型计算应力强度因子范围；
根据门槛应力强度因子及断裂韧度求解裂纹扩展应力强度因子范围；
基于平均应力对长裂纹扩展的影响，建立裂纹扩展速率修正Paris公式；
基于修正后的所述Paris公式、门槛裂纹长度、临界裂纹尺寸建立齿轮疲劳裂纹扩展寿命模型；
根据所述齿轮疲劳裂纹扩展寿命模型计算所述齿轮疲劳裂纹扩展寿命。
4.根据权利要求3所述的齿轮弯曲疲劳寿命预测方法，其特征在于，基于所述齿轮齿根二维有限元模型建立裂纹扩展角计算模型和裂纹扩展增量计算模型，并预测裂纹扩展路径，包括：
基于最大切向应力准则，建立裂纹扩展角计算模型；
根据所述应力强度因子范围及修正后的所述Paris公式建立裂纹扩展增量计算模型；
重复执行所述步骤3至步骤5，直至应力强度因子达到临界应力强度因子；
基于应力强度因子达到临界应力强度因子过程中所得到的裂纹扩展角及裂纹扩展增量预测试件发生失效时的裂纹扩展路径。
5.一种齿轮弯曲疲劳寿命预测装置，其特征在于，包括：
疲劳极限修正单元，基于表面粗糙度建立疲劳极限的修正模型，并根据所述修正模型对材料的疲劳极限进行修正，得到修正后疲劳极限；
门槛裂纹长度确定单元，用于根据门槛应力强度因子幅度及所述修正后疲劳极限确定门槛裂纹长度；
裂纹萌生寿命预测单元，用于创建疲劳裂纹尺寸萌生模型，并基于所述疲劳裂纹尺寸萌生模型预测齿轮疲劳裂纹萌生寿命；
裂纹扩展寿命预测单元，用于基于线弹性断裂力学准则预测齿轮疲劳裂纹扩展寿命；
齿轮弯曲疲劳寿命计算单元，用于根据所述齿轮疲劳裂纹萌生寿命和齿轮疲劳裂纹扩展寿命，建立齿轮弯曲疲劳计算模型，计算齿轮弯曲疲劳寿命。
6.根据权利要求5所述的齿轮弯曲疲劳寿命预测装置，其特征在于，还包括：
几何模型绘制单元，用于根据包含模数、齿数、压力角的齿轮基本参数绘制齿轮齿根二维几何模型；
有限元模型创建单元，用于基于所述齿轮齿根二维几何模型，划分网格，施加边界约束，确定承载工况，建立齿轮齿根二维有限元模型；
拟合函数关系确定单元，用于确定平面应变下裂纹尖端处张开型应力强度因子及滑开型应力强度因子与节点位移的拟合函数关系；
应力强度方程建立单元，用于根据所述拟合函数关系建立平面应变下裂纹尖端复合应力强度因子方程；
裂纹扩展路径预测单元，用于基于所述齿轮齿根二维有限元模型建立裂纹扩展角计算模型和裂纹扩展增量计算模型，并预测裂纹扩展路径。
7.根据权利要求5所述的齿轮弯曲疲劳寿命预测装置，其特征在于，所述裂纹萌生寿命预测单元包括：
应力分布建立模块，用于根据应力梯度下的裂纹萌生模型，建立齿轮局部应力分布关系；
平均应力范围确定模块，用于基于所述齿轮局部应力分布关系，确定作用在裂纹上的平均应力范围；
裂纹尺寸萌生模型创建模块，用于根据所述平均应力范围建立所述疲劳裂纹尺寸萌生模型；
裂纹萌生寿命预测模型创建模块，用于根据所述疲劳裂纹尺寸萌生模型、修正后疲劳极限及门槛裂纹长度建立裂纹萌生寿命预测模型；
裂纹萌生寿命预测模块，用于基于所述裂纹萌生寿命预测模型计算齿轮疲劳裂纹萌生寿命。
8.根据权利要求5所述的齿轮弯曲疲劳寿命预测装置，其特征在于，所述裂纹扩展寿命预测单元包括：
应力强度因子范围计算模块，用于根据应力强度因子范围计算模型计算应力强度因子范围；
裂纹扩展应力强度因子范围计算模块，用于根据门槛应力强度因子及断裂韧度求解裂纹扩展应力强度因子范围；
公式修正模块，用于基于平均应力对长裂纹扩展的影响，建立裂纹扩展速率修正Paris公式；
齿轮疲劳裂纹扩展寿命模型创建模块，用于基于修正后的所述Paris公式、门槛裂纹长度、临界裂纹尺寸建立齿轮疲劳裂纹扩展寿命模型；
齿轮疲劳裂纹扩展寿命计算模块，用于根据所述齿轮疲劳裂纹扩展寿命模型计算所述齿轮疲劳裂纹扩展寿命。
9.根据权利要求8所述的齿轮弯曲疲劳寿命预测装置，其特征在于，所述裂纹扩展路径预测单元包括：
裂纹扩展角计算模块，用于基于最大切向应力准则，建立裂纹扩展角计算模型；
裂纹扩展增量计算模块，用于根据所述应力强度因子范围及修正后的所述Paris公式建立裂纹扩展增量计算模型；
裂纹扩展路径预测模块，用于基于应力强度因子达到临界应力强度因子过程中所得到的裂纹扩展角及裂纹扩展增量预测试件发生失效时的裂纹扩展路径。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710197972.8A CN106886663B (zh) | 2017-03-29 | 2017-03-29 | 齿轮弯曲疲劳寿命预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710197972.8A CN106886663B (zh) | 2017-03-29 | 2017-03-29 | 齿轮弯曲疲劳寿命预测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106886663A true CN106886663A (zh) | 2017-06-23 |
| CN106886663B CN106886663B (zh) | 2020-04-03 |
Family
ID=59182497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710197972.8A Expired - Fee Related CN106886663B (zh) | 2017-03-29 | 2017-03-29 | 齿轮弯曲疲劳寿命预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106886663B (zh) |
Cited By (30)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107389471A (zh) * | 2017-07-13 | 2017-11-24 | 电子科技大学 | 裂纹失效模式判定方法及基于该方法的疲劳寿命预测方法 |
| CN107449894A (zh) * | 2017-08-23 | 2017-12-08 | 中国人民解放军海军工程大学 | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 |
| CN107577902A (zh) * | 2017-10-23 | 2018-01-12 | 哈尔滨工业大学 | 一种基于ukf的飞机疲劳结构剩余寿命预测方法 |
| CN108444697A (zh) * | 2017-12-04 | 2018-08-24 | 重庆大学 | 渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法 |
| CN108627328A (zh) * | 2018-05-24 | 2018-10-09 | 电子科技大学 | 一种耦合应力梯度修正的缺口件多轴疲劳寿命预测方法 |
| CN109165407A (zh) * | 2018-07-18 | 2019-01-08 | 上海工程技术大学 | 一种用于机械构件疲劳裂纹寿命的预估方法 |
| CN109241147A (zh) * | 2018-10-22 | 2019-01-18 | 上海工程技术大学 | 一种用于评估统计值变异性的方法 |
| CN109948216A (zh) * | 2019-03-12 | 2019-06-28 | 华东理工大学 | 总应变能密度修正的缺口件低周疲劳预测方法 |
| CN110059432A (zh) * | 2019-04-26 | 2019-07-26 | 上海工程技术大学 | 一种估算裂纹萌生尺寸数值的方法 |
| CN110222439A (zh) * | 2019-06-12 | 2019-09-10 | 四川大学 | 基于Abaqus平台疲劳损伤与寿命评估方法 |
| CN111027208A (zh) * | 2019-12-10 | 2020-04-17 | 中国航空综合技术研究所 | 在役机械装备关键结构元件疲劳载荷作用下的定延寿方法 |
| CN111209677A (zh) * | 2020-01-13 | 2020-05-29 | 上海工程技术大学 | 一种基于快速系数的铝合金疲劳寿命计算方法 |
| CN111324989A (zh) * | 2020-03-19 | 2020-06-23 | 重庆大学 | 一种基于ga-bp神经网络的齿轮接触疲劳寿命预测方法 |
| CN111737861A (zh) * | 2020-06-09 | 2020-10-02 | 西北工业大学 | 基于不确定性预测试件安全断裂疲劳寿命的方法 |
| CN111751199A (zh) * | 2020-06-09 | 2020-10-09 | 西北工业大学 | 基于eifs分布的疲劳寿命预测方法 |
| CN111783243A (zh) * | 2020-06-18 | 2020-10-16 | 东南大学 | 一种基于滤波算法的金属结构疲劳裂纹扩展寿命预测方法 |
| CN112182807A (zh) * | 2020-10-26 | 2021-01-05 | 东北大学 | 基于有限元仿真的直齿轮齿根裂纹扩展规律分析方法 |
| CN112179765A (zh) * | 2020-06-09 | 2021-01-05 | 西北工业大学 | 基于ct试件的疲劳寿命预测方法 |
| CN112329166A (zh) * | 2020-10-28 | 2021-02-05 | 西安热工研究院有限公司 | 一种风电机组主轴疲劳寿命反推方法 |
| CN112800624A (zh) * | 2021-02-08 | 2021-05-14 | 上海工程技术大学 | 基于路径规划算法的复合材料寿命预测方法 |
| CN112881167A (zh) * | 2020-09-10 | 2021-06-01 | 北京理工大学 | 一种微动疲劳裂纹萌生位置判定方法 |
| CN112903435A (zh) * | 2021-01-19 | 2021-06-04 | 南京航空航天大学 | 基于单个缺陷的多轴疲劳寿命预测方法 |
| CN113158496A (zh) * | 2021-05-21 | 2021-07-23 | 北部湾大学 | 储液罐疲劳寿命预测方法 |
| JP6944736B1 (ja) * | 2021-01-08 | 2021-10-06 | 株式会社寿ホールディングス | 疲労寿命予測方法、疲労寿命予測装置、疲労寿命予測プログラムおよび記憶媒体 |
| CN113591268A (zh) * | 2021-06-24 | 2021-11-02 | 内蒙古工业大学 | 变幅载荷下齿轮接触疲劳寿命可靠性评估方法及装置 |
| CN113704907A (zh) * | 2021-08-13 | 2021-11-26 | 湖南磐钴传动科技有限公司 | 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 |
| CN114216803A (zh) * | 2021-09-30 | 2022-03-22 | 中国航发北京航空材料研究院 | 一种金属材料的高周疲劳全寿命预测方法 |
| CN115169694A (zh) * | 2022-07-06 | 2022-10-11 | 天津大学 | 一种实现亚临界裂纹动态扩展与寿命预测的方法 |
| CN117057166A (zh) * | 2023-10-11 | 2023-11-14 | 合肥通用机械研究院有限公司 | 应力集中部位裂纹自由表面处应力强度因子的计算方法 |
| CN113704907B (zh) * | 2021-08-13 | 2024-06-11 | 湖南磐钴传动科技有限公司 | 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103616179A (zh) * | 2013-12-05 | 2014-03-05 | 广西大学 | 一种基于缺陷建模的传动齿轮疲劳寿命评估方法 |
| CN104281782A (zh) * | 2014-10-13 | 2015-01-14 | 北京理工大学 | 基于缺口试件的啮合齿轮弯曲疲劳极限评估方法及装置 |
-
2017
- 2017-03-29 CN CN201710197972.8A patent/CN106886663B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103616179A (zh) * | 2013-12-05 | 2014-03-05 | 广西大学 | 一种基于缺陷建模的传动齿轮疲劳寿命评估方法 |
| CN104281782A (zh) * | 2014-10-13 | 2015-01-14 | 北京理工大学 | 基于缺口试件的啮合齿轮弯曲疲劳极限评估方法及装置 |
Non-Patent Citations (4)
| Title |
|---|
| HAILONG DENG 等: "Very High Cycle Fatigue Failure Analysis and Life Prediction of Cr-Ni-W Gear Steel Based on Crack Initiation and Growth Behaviors", 《MATERIALS》 * |
| 林腾蛟 等: "圆柱齿轮齿根三维裂纹扩展分析及寿命预测", 《重庆大学学报》 * |
| 沈士明: "《含缺陷承压设备安全分析技术》", 31 March 2011, 北京：中国石化出版社 * |
| 王延忠 等: "航空重载面齿轮三维裂纹分析与疲劳寿命预测", 《北京航空航天大学学报》 * |
Cited By (44)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107389471A (zh) * | 2017-07-13 | 2017-11-24 | 电子科技大学 | 裂纹失效模式判定方法及基于该方法的疲劳寿命预测方法 |
| CN107449894B (zh) * | 2017-08-23 | 2019-08-20 | 中国人民解放军海军工程大学 | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 |
| CN107449894A (zh) * | 2017-08-23 | 2017-12-08 | 中国人民解放军海军工程大学 | 氮化硅结合碳化硅耐火陶瓷的疲劳寿命预测方法 |
| CN107577902A (zh) * | 2017-10-23 | 2018-01-12 | 哈尔滨工业大学 | 一种基于ukf的飞机疲劳结构剩余寿命预测方法 |
| CN107577902B (zh) * | 2017-10-23 | 2020-11-13 | 哈尔滨工业大学 | 一种基于ukf的飞机疲劳结构剩余寿命预测方法 |
| CN108444697A (zh) * | 2017-12-04 | 2018-08-24 | 重庆大学 | 渗碳硬化齿轮弹塑性接触疲劳点蚀与剥落风险预测方法 |
| CN108627328A (zh) * | 2018-05-24 | 2018-10-09 | 电子科技大学 | 一种耦合应力梯度修正的缺口件多轴疲劳寿命预测方法 |
| CN109165407B (zh) * | 2018-07-18 | 2023-08-08 | 上海工程技术大学 | 一种用于机械构件疲劳裂纹寿命的预估方法 |
| CN109165407A (zh) * | 2018-07-18 | 2019-01-08 | 上海工程技术大学 | 一种用于机械构件疲劳裂纹寿命的预估方法 |
| CN109241147A (zh) * | 2018-10-22 | 2019-01-18 | 上海工程技术大学 | 一种用于评估统计值变异性的方法 |
| CN109241147B (zh) * | 2018-10-22 | 2022-07-01 | 上海工程技术大学 | 一种用于评估统计值变异性的方法 |
| CN109948216A (zh) * | 2019-03-12 | 2019-06-28 | 华东理工大学 | 总应变能密度修正的缺口件低周疲劳预测方法 |
| CN109948216B (zh) * | 2019-03-12 | 2023-01-03 | 华东理工大学 | 总应变能密度修正的缺口件低周疲劳预测方法 |
| CN110059432A (zh) * | 2019-04-26 | 2019-07-26 | 上海工程技术大学 | 一种估算裂纹萌生尺寸数值的方法 |
| CN110222439A (zh) * | 2019-06-12 | 2019-09-10 | 四川大学 | 基于Abaqus平台疲劳损伤与寿命评估方法 |
| CN111027208A (zh) * | 2019-12-10 | 2020-04-17 | 中国航空综合技术研究所 | 在役机械装备关键结构元件疲劳载荷作用下的定延寿方法 |
| CN111027208B (zh) * | 2019-12-10 | 2022-02-01 | 中国航空综合技术研究所 | 在役机械装备关键结构元件疲劳载荷作用下的定延寿方法 |
| CN111209677B (zh) * | 2020-01-13 | 2022-03-25 | 上海工程技术大学 | 一种基于快速系数的铝合金疲劳寿命计算方法 |
| CN111209677A (zh) * | 2020-01-13 | 2020-05-29 | 上海工程技术大学 | 一种基于快速系数的铝合金疲劳寿命计算方法 |
| CN111324989B (zh) * | 2020-03-19 | 2024-01-30 | 重庆大学 | 一种基于ga-bp神经网络的齿轮接触疲劳寿命预测方法 |
| CN111324989A (zh) * | 2020-03-19 | 2020-06-23 | 重庆大学 | 一种基于ga-bp神经网络的齿轮接触疲劳寿命预测方法 |
| CN112179765A (zh) * | 2020-06-09 | 2021-01-05 | 西北工业大学 | 基于ct试件的疲劳寿命预测方法 |
| CN111737861A (zh) * | 2020-06-09 | 2020-10-02 | 西北工业大学 | 基于不确定性预测试件安全断裂疲劳寿命的方法 |
| CN111751199A (zh) * | 2020-06-09 | 2020-10-09 | 西北工业大学 | 基于eifs分布的疲劳寿命预测方法 |
| CN111751199B (zh) * | 2020-06-09 | 2021-07-09 | 西北工业大学 | 基于eifs分布的疲劳寿命预测方法 |
| CN111737861B (zh) * | 2020-06-09 | 2022-02-01 | 西北工业大学 | 基于不确定性预测试件安全断裂疲劳寿命的方法 |
| CN111783243A (zh) * | 2020-06-18 | 2020-10-16 | 东南大学 | 一种基于滤波算法的金属结构疲劳裂纹扩展寿命预测方法 |
| CN111783243B (zh) * | 2020-06-18 | 2021-02-26 | 东南大学 | 一种基于滤波算法的金属结构疲劳裂纹扩展寿命预测方法 |
| CN112881167A (zh) * | 2020-09-10 | 2021-06-01 | 北京理工大学 | 一种微动疲劳裂纹萌生位置判定方法 |
| CN112182807A (zh) * | 2020-10-26 | 2021-01-05 | 东北大学 | 基于有限元仿真的直齿轮齿根裂纹扩展规律分析方法 |
| CN112329166A (zh) * | 2020-10-28 | 2021-02-05 | 西安热工研究院有限公司 | 一种风电机组主轴疲劳寿命反推方法 |
| WO2022149297A1 (ja) * | 2021-01-08 | 2022-07-14 | 株式会社寿ホールディングス | 疲労寿命予測方法、疲労寿命予測装置、疲労寿命予測プログラムおよび記憶媒体 |
| JP6944736B1 (ja) * | 2021-01-08 | 2021-10-06 | 株式会社寿ホールディングス | 疲労寿命予測方法、疲労寿命予測装置、疲労寿命予測プログラムおよび記憶媒体 |
| CN112903435A (zh) * | 2021-01-19 | 2021-06-04 | 南京航空航天大学 | 基于单个缺陷的多轴疲劳寿命预测方法 |
| CN112800624B (zh) * | 2021-02-08 | 2022-09-13 | 上海工程技术大学 | 基于路径规划算法的复合材料寿命预测方法 |
| CN112800624A (zh) * | 2021-02-08 | 2021-05-14 | 上海工程技术大学 | 基于路径规划算法的复合材料寿命预测方法 |
| CN113158496A (zh) * | 2021-05-21 | 2021-07-23 | 北部湾大学 | 储液罐疲劳寿命预测方法 |
| CN113591268A (zh) * | 2021-06-24 | 2021-11-02 | 内蒙古工业大学 | 变幅载荷下齿轮接触疲劳寿命可靠性评估方法及装置 |
| CN113704907A (zh) * | 2021-08-13 | 2021-11-26 | 湖南磐钴传动科技有限公司 | 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 |
| CN113704907B (zh) * | 2021-08-13 | 2024-06-11 | 湖南磐钴传动科技有限公司 | 一种基于齿面应力场的齿轮接触疲劳寿命预测方法 |
| CN114216803A (zh) * | 2021-09-30 | 2022-03-22 | 中国航发北京航空材料研究院 | 一种金属材料的高周疲劳全寿命预测方法 |
| CN115169694A (zh) * | 2022-07-06 | 2022-10-11 | 天津大学 | 一种实现亚临界裂纹动态扩展与寿命预测的方法 |
| CN117057166A (zh) * | 2023-10-11 | 2023-11-14 | 合肥通用机械研究院有限公司 | 应力集中部位裂纹自由表面处应力强度因子的计算方法 |
| CN117057166B (zh) * | 2023-10-11 | 2023-12-26 | 合肥通用机械研究院有限公司 | 应力集中部位裂纹自由表面处应力强度因子的计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106886663B (zh) | 2020-04-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106886663A (zh) | 齿轮弯曲疲劳寿命预测方法及装置 | |
| McCartney | Mechanics of matrix cracking in brittle-matrix fibre-reinforced composites | |
| CN107884290B (zh) | 一种考虑磨损影响下的微动疲劳裂纹扩展寿命预测方法 | |
| Gheitasi et al. | Slenderness classification of unstiffened metal plates under shear loading | |
| CN105803941B (zh) | 一种斜拉桥双索同时张拉无应力状态施工控制方法 | |
| Yan et al. | Ductile fracture simulation of constructional steels based on yield-to-fracture stress–strain relationship and micromechanism-based fracture criterion | |
| CN104809311A (zh) | 一种基于多因素融合修正的结构件剩余寿命预测方法 | |
| Kadkhodayan et al. | A new fictitious time for the dynamic relaxation (DXDR) method | |
| Deng et al. | Nonlinear analysis of reinforced concrete beam bending failure experimentation based on ABAQUS | |
| Mendes et al. | A simplified reinforcing steel model suitable for cyclic loading including ultra-low-cycle fatigue effects | |
| CN107014541A (zh) | 基于线性模型修正的等代铰接梁模型索力测量方法 | |
| Wu et al. | Random fatigue damage accumulation analysis of composite thin-wall structures based on residual stiffness method | |
| Muñoz-Reja et al. | A numerical implementation of the Coupled Criterion of Finite Fracture Mechanics for elastic interfaces | |
| Malíková et al. | Significance of higher-order terms of the Williams expansion for plastic zone extent estimation demonstrated on a mixed-mode geometry | |
| JP2015132564A (ja) | 熱変形解析方法、熱変形解析プログラム、および熱変形解析装置 | |
| Yang et al. | Numerical study on steady-state cracking of composites | |
| Wen et al. | A numerical algorithm for multiple cracks propagation in concrete structure | |
| CN102733525B (zh) | 一种基于承载全过程开口式整体张拉结构多阶段设计方法 | |
| CN110501127B (zh) | 一种基于损伤状态倾角斜率的等截面梁损伤识别方法 | |
| CN111429981B (zh) | 针对单晶材料热-机疲劳本构模型的损伤计算方法 | |
| CN112307547B (zh) | 隧道掌子面支护压力设计方法 | |
| Sandovic et al. | Experimental and analytical investigation of deformations and stress distribution in steel bands of a two-span stress-ribbon pedestrian bridge | |
| CN110501177B (zh) | 基于自由端倾角影响线曲率的悬臂梁损伤识别方法 | |
| He et al. | Damage constitutive model and numerical simulation of in-plane tear failure for architectural fabric membranes | |
| CN109522570B (zh) | 一种基于指数Weibull方程的混凝土疲劳变形演化模型 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200403 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |