CN105716935B - 预测试样尺寸对疲劳寿命影响的方法 - Google Patents
预测试样尺寸对疲劳寿命影响的方法 Download PDFInfo
- Publication number
- CN105716935B CN105716935B CN201610059873.9A CN201610059873A CN105716935B CN 105716935 B CN105716935 B CN 105716935B CN 201610059873 A CN201610059873 A CN 201610059873A CN 105716935 B CN105716935 B CN 105716935B
- Authority
- CN
- China
- Prior art keywords
- fatigue life
- stress
- specimen
- life
- bulk specimen
- 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
Links
- 0 CCC*C*NC Chemical compound CCC*C*NC 0.000 description 2
- UXXWCYCBIWTDQB-UHFFFAOYSA-N CC1C=CC2C1C2 Chemical compound CC1C=CC2C1C2 UXXWCYCBIWTDQB-UHFFFAOYSA-N 0.000 description 1
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
-
- 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)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开一种预测试样尺寸对疲劳寿命影响的方法,包括通过寿命与应力之间的函数关系,将不同应力水平下的疲劳寿命转化为同一应力水平下的寿命;对同一应力下的疲劳寿命进行统计分析,确定其最服从的概率分布形式;计算小试样和大试样的控制体积,确定大试样与小试样控制体积之比;大试样的疲劳寿命NL由同应力下小试样的疲劳寿命N的n次Monte Carlo模拟获得,即NL=min{N1,N2,…,Nn};获得大试样在同一应力下多个疲劳寿命数据;获得大试样在不同应力下的疲劳寿命数据;对获得的大试样疲劳寿命数据进行统计分析,得出P‑S‑N曲线。有效克服大试样疲劳实验费时、费力的缺点。
Description
技术领域
本发明涉及疲劳寿命预测理论方法领域,尤其是一种预测试样尺寸对疲劳寿命影响的方法。
背景技术
疲劳失效分析与寿命预测是工程结构完整性的关键力学问题之一。大量事实表明,大尺寸试样往往比相同加载条件下的小试样具有较低的疲劳寿命。疲劳实验耗时、费力,特别是大尺寸试样长寿命疲劳实验。因此,发展通过小试样的疲劳实验数据对大试样的疲劳寿命和P-S-N(失效概率-应力-寿命)曲线进行预测具有十分重要的应用价值。
发明内容
本发明提供一种预测试样尺寸对疲劳寿命影响的方法,用于克服现有技术中的缺陷,方便、快捷获得不同尺寸试样的疲劳寿命。
本发明提供一种预测试样尺寸对疲劳寿命影响的方法,包括以下步骤:
步骤1,通过寿命与应力之间的函数关系,将不同应力水平下的疲劳寿命转化为同一应力水平下的寿命;
步骤2,对同一应力下的疲劳寿命进行统计分析,确定其最服从的概率分布形式;
步骤3,计算小试样和大试样的控制体积,确定大试样与小试样控制体积之比;
步骤4,大试样的疲劳寿命NL由同应力下小试样的疲劳寿命N的n次Monte Carlo模拟获得,即NL=min{N1,N2,…,Nn};
步骤5,重复步骤4,获得大试样在同一应力下多个疲劳寿命数据;
步骤6,重复5,获得大试样在不同应力下的疲劳寿命数据;
步骤7,对获得的大试样疲劳寿命数据进行统计分析,得出P-S-N曲线。
本发明提供的预测试样尺寸对疲劳寿命影响的方法,能够通过少数小试样的疲劳实验数据获得大试样的疲劳寿命和P-S-N曲线结果,有效克服大试样(特别是大尺寸构件)疲劳实验高时间成本、高经费成本的缺点,特别是超高周范畴(107周次及以上寿命)的疲劳实验,为方便、快捷获得不同尺寸试样的疲劳寿命提供了方法。
附图说明
图1a是本发明实施例提供的预测试样尺寸对疲劳寿命影响的方法中弹簧钢小试样的形状及尺寸示意图;
图1b是本发明实施例提供的预测试样尺寸对疲劳寿命影响的方法中弹簧钢大试样的形状及尺寸示意图;
图2是图1a中弹簧钢小试样S-N数据和P-S-N曲线预测结果示意图;
图3是图2中弹簧钢大试样S-N数据与Monte Carlo预测结果比较示意图;
图4是图2中弹簧钢大试样S-N数据与P-S-N曲线预测结果比较示意图。
具体实施方式
实施例
本发明提供一种预测试样尺寸对疲劳寿命影响的方法,并选用已有文献中的一组小试样和大试样实验数据,表明本发明方法的一般性和有效性。主要包括以下步骤:
步骤1,首先确定小试样的疲劳寿命概率分布函数,本实施例中对弹簧钢的大试样(参见图1b)和同材质的弹簧钢小试样(图1a)为例进行说明和具体实验,具体步骤如下:
步骤11,通过寿命N与应力σ之间的函数N(σ),如N(σ)=Aiσa(Ai为除 应力之外其他因素的影响;a为常数,可通过最小二乘法获得),将不同应力水平下的疲劳寿命转化为同一应力水平下的寿命。寿命与应力之间的函数跟试样的材质有直接关系,可根据具体的材质,通过试验及理论计算获得寿命与应力之间的函数关系。
作为实际应用的例子,选取文献中一种弹簧钢小试样(参见图1a)及其寿命与应力的试验数据(参见图2中的方块■),寿命与应力关联的函数取为:N(σ)=Aiσa则任意应力水平σk(k=1,2,...,n)下的寿命Nk均可转化为某一给定应力σ′下的寿命,即
或
通过最小二乘法得到a=-30.11;
通过该步骤将任意应力下的寿命转换为任意的同一应力下的寿命数据;
若同一应力下已有多个疲劳寿命数据供统计分析,则不需要此步骤,直接对上述某一给定应力下的多个疲劳寿命数据进行统计分析,确定这些数据最服从的概率分布形式;如果在对小试样进行试验时,选择同一应力(例如800MPa)下进行,可以通过上述步骤11获得其他同一应力(例如700MPa)下的疲劳寿命数据。
步骤12,根据转化后的某一给定应力σ′下的疲劳寿命选择一种概率分布形式,如正态分布形式或Weibull分布形式;
本实施例中对弹簧钢小试样在同一应力下的对数疲劳寿命服从双参数Weibull分布时的P-S-N曲线进行预测,参见图2,其中包含概率分别是50%、95%、99%及99.9%的P-S-N曲线,其中方块(■)是对弹簧钢小试样进行的疲劳寿命实验数据,通过比较可以看出:实验数据基本分布在概率50%的P-S-N预测曲线周围;
步骤2,计算弹簧钢小试样和弹簧钢大试样的控制体积(这里是指疲劳裂纹可能萌生的高应力区域,通常取应力≥90%最大应力的区域,也可以根据预测要求取不同的值),确定大试样与小试样控制体积之比:n=VL/VS(n取整数);弹簧钢大试样的示意图参见图1b,采用有限元计算得到弹簧钢小试样的控制体积为VS=33mm3,大试样的控制体积为VL=912mm3,求出控 制体积之比n=VL/VS=28;如果试样的形状规则,还可采用理论计算方法求出控制体积;无论试样的形状如何,均可采用有限元计算方法求出控制体积,适用范围广。
步骤3,由弹簧钢小试样的疲劳寿命数据,获得不同应力下弹簧钢大试样的疲劳寿命数据,对其进行统计分析获得P-S-N曲线;
步骤31,弹簧钢大试样的疲劳寿命NL由同应力下弹簧钢小试样的疲劳寿命N的n次Monte Carlo模拟获得,即NL=min{N1,N2,…,Nn};
步骤32,重复步骤31,获得弹簧钢大试样在同一应力下多个疲劳寿命数据;图3给出了采用弹簧钢小试样实验数据和Monte Carlo模拟对弹簧钢大试样(VL=912mm3)S-N数据的预测结果与实验结果的比较示意图;其中方块(■)是弹簧钢小试样的疲劳寿命实验数据,三角(△)是弹簧钢大试样的疲劳寿命实验数据,叉号(×)是弹簧钢大试样疲劳寿命的预测数据;经比较可见预测结果与实验结果吻合。
步骤33,重复步骤32,获得弹簧钢大试样在不同应力下的疲劳寿命数据;
步骤34,对获得的大试样疲劳寿命数据进行统计分析,得出P-S-N曲线预测结果,其中包含概率分别是50%、95%、99%及99.9%的P-S-N曲线,其中三角(△)是弹簧钢大试样的疲劳寿命的实验数据。经比较可以看出,实验数据分布在概率50%的P-S-N曲线周围。
上述方法,其中大试样和小试样需要同种材质。(注:所选用的小试样和大试样实验数据引自“Furuya,Y.Notable size effects on very high cycle fatigueproperties of high-strength steel.Mater.Sci.Eng.A,2011,528:5234–5240.) 。
Claims (5)
1.一种预测试样尺寸对疲劳寿命影响的方法,其特征在于,包括以下步骤:
步骤1,通过寿命与应力之间的函数关系,将不同应力水平下的疲劳寿命转化为同一应力水平下的寿命;
步骤2,对同一应力下的疲劳寿命进行统计分析,确定其最服从的概率分布形式;
步骤3,计算小试样和大试样的控制体积,确定大试样与小试样控制体积之比;
步骤4,大试样的疲劳寿命NL由同应力下小试样的疲劳寿命N的n次Monte Carlo模拟获得,即NL=min{N1,N2,…,Nn};其中,n为大试样与小试样的控制体积之比且n取整数;
步骤5,重复步骤4,获得大试样在同一应力下多个疲劳寿命数据;
步骤6,重复步骤5,获得大试样在不同应力下的疲劳寿命数据;
步骤7,对获得的大试样疲劳寿命数据进行统计分析,得出P-S-N曲线。
2.根据权利要求1所述的预测试样尺寸对疲劳寿命影响的方法,其特征在于,所述小试样和大试样的材质相同,均是金属材料;其中步骤1中,寿命与应力之间的函数为:N(σ)=Aiσa;
其中,寿命为N,应力为σ,Ai为除应力之外的其他影响因素;a为常数,可通过最小二乘法获得。
3.根据权利要求2所述的预测试样尺寸对疲劳寿命影响的方法,其特征在于,其中步骤2中,概率分布形式为正态分布或双参数Weibull分布。
4.根据权利要求1所述的预测试样尺寸对疲劳寿命影响的方法,其特征在于,其中步骤3中,采用有限元计算得到小试样和大试样的控制体积。
5.根据权利要求1-4任一所述的预测试样尺寸对疲劳寿命影响的方法,其特征在于,所述小试样和大试样的控制体积为疲劳裂纹可能萌生的区域,所述区域的所受应力为大于或等于最大应力的90%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610059873.9A CN105716935B (zh) | 2016-01-28 | 2016-01-28 | 预测试样尺寸对疲劳寿命影响的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610059873.9A CN105716935B (zh) | 2016-01-28 | 2016-01-28 | 预测试样尺寸对疲劳寿命影响的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105716935A CN105716935A (zh) | 2016-06-29 |
CN105716935B true CN105716935B (zh) | 2018-02-16 |
Family
ID=56154327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610059873.9A Active CN105716935B (zh) | 2016-01-28 | 2016-01-28 | 预测试样尺寸对疲劳寿命影响的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105716935B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106294953B (zh) * | 2016-08-02 | 2019-07-23 | 浙江大学 | 一种考虑应力水平的水泥基材料疲劳失效概率和疲劳寿命概率模型的建立方法 |
CN107145641B (zh) * | 2017-04-13 | 2020-04-24 | 南京航空航天大学 | 叶片振动疲劳概率寿命预测方法 |
CN108645706B (zh) * | 2018-04-28 | 2020-12-29 | 中国科学院金属研究所 | 一种通过硬化强度和抗拉强度预测金属材料疲劳强度的方法 |
CN108871981B (zh) * | 2018-05-25 | 2021-04-20 | 中国科学院金属研究所 | 一种金属材料疲劳性能优劣的评价方法 |
CN109086541B (zh) * | 2018-08-14 | 2021-03-26 | 中国科学院力学研究所 | 预测试样几何形状和加载方式对疲劳强度影响的方法 |
CN109145426B (zh) * | 2018-08-14 | 2021-04-09 | 中国科学院力学研究所 | 预测试样几何形状和加载方式对疲劳寿命影响的方法及装置 |
CN109470549B (zh) * | 2018-09-07 | 2020-07-28 | 北京航空航天大学 | 增材制造材料p-s-n曲线表征方法及其应用 |
CN113627049B (zh) * | 2021-07-16 | 2023-08-18 | 武汉钢铁有限公司 | 超高强钢大试样超长寿命疲劳强度评估方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100533119C (zh) * | 2007-08-23 | 2009-08-26 | 南京航空航天大学 | 机械结构的裂纹扩展率和裂纹扩展寿命预测方法 |
US8155940B2 (en) * | 2008-07-30 | 2012-04-10 | GM Global Technology Operations LLC | Methods and systems for predicting very high cycle fatigue properties in metal alloys |
CN102384856B (zh) * | 2011-08-15 | 2014-01-22 | 东南大学 | 基于概率有限元法的钢桥疲劳可靠度评估方法 |
CN103698222B (zh) * | 2014-01-20 | 2015-09-30 | 哈尔滨工程大学 | 一种利用锅炉壁温、应力测量装置的锅炉疲劳寿命测定方法 |
CN104034599A (zh) * | 2014-05-20 | 2014-09-10 | 北京航空航天大学 | 一种测试发动机材料极高温疲劳s-n曲线的试验系统 |
CN104951660A (zh) * | 2015-07-06 | 2015-09-30 | 北京航空航天大学 | 一种确定飞机结构首翻期和检修周期的方法 |
-
2016
- 2016-01-28 CN CN201610059873.9A patent/CN105716935B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105716935A (zh) | 2016-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105716935B (zh) | 预测试样尺寸对疲劳寿命影响的方法 | |
Cao | Models for ductile damage and fracture prediction in cold bulk metal forming processes: a review | |
Bagherifard et al. | On the shot peening surface coverage and its assessment by means of finite element simulation: A critical review and some original developments | |
Schmitz et al. | A probabilistic model for LCF | |
Li et al. | Quantification of the interaction within defect populations on fatigue behavior in an aluminum alloy | |
Agius et al. | Sensitivity and optimisation of the Chaboche plasticity model parameters in strain-life fatigue predictions | |
CN111339700B (zh) | 核电汽轮机叶片的疲劳损伤评估方法、装置和存储介质 | |
CN109933925B (zh) | 一种金属板材的冲压成形性能预测方法 | |
WO2013179877A1 (ja) | 部品の破壊評価装置、部品の破壊評価方法、及びコンピュータプログラム | |
Wan et al. | Effect of micro-porosities on fatigue behavior in aluminum die castings by 3D X-ray tomography inspection | |
Yan et al. | Constraining evolution in the halo model using galaxy redshift surveys | |
Yuan et al. | A modified constitutive equation for elevated temperature flow behavior of Ti–6Al–4V alloy based on double multiple nonlinear regression | |
CN107506544A (zh) | 一种金属材料疲劳裂纹扩展预测方法 | |
An et al. | Modified Arrhenius-type constitutive model and artificial neural network-based model for constitutive relationship of 316LN stainless steel during hot deformation | |
Fujii et al. | Strain-based approach to investigate intergranular stress corrosion crack initiation on a smooth surface of austenitic stainless steel | |
CN103308372B (zh) | 基于透射电子显微镜的t91钢老化评级方法 | |
CN103439473A (zh) | 一种12Cr1MoV钢受热面状态评估方法 | |
Laz et al. | A probabilistic total fatigue life model incorporating material inhomogeneities, stress level and fracture mechanics | |
CN109948216B (zh) | 总应变能密度修正的缺口件低周疲劳预测方法 | |
CN102788748B (zh) | 82b盘条索氏体含量分析用标准样品的制备方法 | |
Halilovič et al. | Prediction of elastic strain recovery of a formed steel sheet considering stiffness degradation | |
Yan et al. | Comparison between impression creep and uni-axial tensile creep performed on nickel-based single crystal superalloys | |
CN113358678B (zh) | α钛变形过程介观应力和织构的半定量预测及可视化方法 | |
Luo et al. | Statistical aspects of fatigue behaviour in a PM Ni-base superalloy Udimet 720 | |
Nioaţă | Researches regarding the optimization of thermal treatment depending on hardness for maraging 300 steel |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |