CN107273649B - 一种脆性材料在高温蠕变状态下失效概率的预测方法 - Google Patents
一种脆性材料在高温蠕变状态下失效概率的预测方法 Download PDFInfo
- Publication number
- CN107273649B CN107273649B CN201710699545.XA CN201710699545A CN107273649B CN 107273649 B CN107273649 B CN 107273649B CN 201710699545 A CN201710699545 A CN 201710699545A CN 107273649 B CN107273649 B CN 107273649B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msup
- creep
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- 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
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种脆性材料在高温蠕变状态下失效概率的预测方法，其在现有技术的基础上，结合脆性材料内部缺陷随机分布的自然属性，假定单轴蠕变失效应变服从威布尔分布，利用单轴蠕变试验获得单轴蠕变失效应变的概率密度分布曲线，通过单轴与多轴蠕变失效应变的转化关系，得到多轴蠕变失效应变的概率密度函数，进而积分得到失效概率计算模型；在此基础上，结合蠕变‑损伤本构方程，利用Fortran语言编写子程序并嵌入到有限元软件中，进而得到脆性材料在高温蠕变状态下的失效概率的预测结果。本发明解决了现有技术不能进行脆性材料在高温蠕变状态下的可靠性预测的技术问题，所获得的预测结果，真实、准确、合理和可靠。
Description
技术领域
本发明涉及可靠性工程技术领域，具体涉及一种脆性材料在高温蠕变状态下失效概率的预测方法。
背景技术
当今国内外的失效评定工作主要采用确定性断裂力学方法的“合乎使用”原则，该方法取结构、缺陷、材料性能等参数的某一个给定值，配合一定的安全系数进行分析，给出安全或不安全的评定结果。
然而，在实际工程中，脆性材料内部缺陷随机分布，其结构尺寸、材料性能参数、载荷等也都有不确定性，可视为具有一定分布的随机变量。
因此，确定性断裂力学将所有参量都作为单值确定量的处理方法，会使评定结构与实际情况产生较大偏差甚至得到错误的评定结果。
为了研究各种不确定性因素对结构失效的影响，定量评估含缺陷结构的安全性，出现了概率断裂力学评定方法。
概率断裂力学将不确定性变量视作服从一定分布的随机变量，采用失效概率表示危险程度，为工程应用中评价构件安全程度提供精确的定量指标，并可以应用这种理论和方法指导可靠性设计和寿命预测。
现有的威布尔分布(即Weibull分布，也称作韦伯分布或韦氏分布)失效概率计算表达式以应力为基础，然而，脆性材料在高温蠕变状态下，不可避免地会发生应力松弛效应，应力迅速减小，接近于零。此时，若采用基于应力的计算表达式来计算失效概率，将会产生很大的偏差，甚至会出现相反的结论。
因此，现有技术的威布尔失效概率计算表达式并不适合评价脆性材料在高温蠕变状态下的可靠性，需要建立新型的失效概率计算模型。
发明内容
为了解决现有技术中威布尔失效概率计算表达式存在的不足，本发明旨在根据威布尔理论以及脆性材料单轴蠕变失效应变呈概率分布的自然属性，以获得新型的失效概率预测公式，从而更准确地预测脆性材料在高温蠕变状态下的失效概率。
本发明为解决上述技术问题所采用的技术方案是，一种脆性材料在高温蠕变状态下失效概率的预测方法，其特征在于，包括以下步骤：
第一步，根据脆性材料内部缺陷随机分布的自然属性，假定反应脆性材料属性的单轴蠕变失效应变εf服从威布尔分布；则单轴蠕变失效应变的概率密度函数f(εf)满足下式(1)：
上式(1)中：
η为变量的尺度参数，η＞0；
β为变量的形状参数，β＞0；
第二步，根据下式(2)(此式为现有技术中本领域公知的公式)所示的单轴与多轴蠕变失效应变εf *的转化关系，根据数学转换关系，得到如下式(3)所示的多轴蠕变失效应变的概率密度分布函数f(εf *)：
上式(2)中：
σm是指材料所承受的静水应力；
σeq为米塞斯应力(即von Mises应力)；
n表示蠕变指数；
为与单轴蠕变失效应变无关的系数；得出多轴蠕变失效应变εf *服从威布尔分布，多轴蠕变失效应变的概率密度分布函数的数学表达式(3)为：
第三步，依据结构失效的条件为等效蠕变应变值εe大于多轴蠕变失效应变值εf *的原则，对多轴蠕变失效应变的概率密度分布函数的数学表达式(3)进行积分，即得到如下式(4)所示的失效概率的计算表达式：
在此基础上，考虑到材料内部缺陷的不同，对于体积为V的脆性材料试样，考虑到体积效应，相应的失效概率表达式为下式(5)：
上式(5)中：
V0为特征体积；
第四步，在相同的试验条件下，对若干组体积为V0的试样在相同应力水平下进行单轴蠕变断裂试验，记录每个断裂蠕变应变值，并以蠕变断裂应变为横坐标，在某一个蠕变断裂应变区间的断裂试样数量为纵坐标，绘制出单轴蠕变失效应变值累积分布直方图；
第五步，根据绘制出的单轴蠕变失效应变值累积分布直方图，用每个蠕变断裂应变区间断裂试样的数量除以总的断裂试样数量，即为体积V0的试样在该区间内的断裂概率值PF0，将V0和PF0带入上述的失效概率计算公式(4)并两边取两次对数，得到：
ln[-ln(1-PF0)]＝βlnεe-lnηβ (6)
根据各试样在相同应力水平下进行单轴蠕变断裂的试验结果，做出ln[-ln(1-PF0)]与lnεe的曲线，并进行线性回归，所得到的直线的斜率即为参数β，根据所得到的直线与y轴的截距能够得到参数η；
第六步，根据上式(5)，结合蠕变-损伤本构方程，利用Fortran语言，编写子程序并嵌入到有限元软件ABAQUS中，即得到脆性材料在高温蠕变状态下失效概率的预测结果。
其中，蠕变-损伤本构方程如下所示：
式中，为蠕变应变，σI为最大主应力B为蠕变第二阶段的常数，β0是与应力相关的函数，ρ是微裂纹损伤参数，ω为蠕变损伤量。
优选地，上述的脆性材料在高温蠕变状态下失效概率的预测方法，其第四步中所述的若干组，优选为10～20组。
上述技术方案直接带来的技术效果是，为更好地理解本发明的技术特点，下面简要说明本发明的技术原理和理论依据。
上述技术方案的理论依据为，单轴蠕变失效应变是反应脆性材料自身蠕变性能的参数，由于脆性材料内部的缺陷分布具有随机性，那么单轴蠕变试验得到的单轴蠕变失效应变值也具有不确定性，而威布尔分布具有较强的拟合能力，在含缺陷结构的可靠性分析领域具有很强的适应性。
因此，可以假设单轴蠕变失效应变服从威布尔分布，并且概率密度分布函数中的尺度参数η的大小表征了分布分散程度的大小，形状参数β取不同的值，可分别得到正、负偏差及对称的概率密度函数。
由于威布尔理论需要考虑最弱链假设，即结构在恒定的单轴载荷下，认为它类似于拉伸的N链，每条链都有不同的失效强度，当最弱链失效时，整个结构失效，因此链的强度与最弱链相关。每条链的失效强度不同，取决于试样内部的缺陷不同，即“体积效应”：试样体积越大，内部缺陷越大，对应的产生较大的应力强度。
因此体积为V的试样，对应的失效概率表达式为：
亦即，上述技术方案的失效概率预测计算模型对现有技术的基于应力的失效概率计算模型进行了科学合理的校正。
优选地，上述的脆性材料在高温蠕变状态下失效概率的预测方法，其第四步中所述的若干组，优选为10-20组。该优选技术方案直接带来的技术效果是，我们的经验表明，兼顾结果的可靠性与工作效率，在相同的试验条件下，对10-20组体积为V0的试样在相同应力水平下进行单轴蠕变断裂试验即可获得比较理想的预测结果。
实践表明，本发明相对于现有技术具有如下有益效果：
1、本发明较好地解决了现有技术不能进行脆性材料在高温蠕变状态下的可靠性预测的技术问题。
2、所获得的预测结果真实、准确、合理和可靠。
附图说明
图1为基于应变的失效概率预测方法流程图。
图2为试样体积大小与缺陷大小的关系示意图。
图3为单轴蠕变失效应变累积分布直方图。
图4为实施例1中600℃下玻璃陶瓷GC-9材料等效蠕变应变和米塞斯应力(Mises应力)随蠕变时间变化的关系曲线。
图5为实施例1中600℃下玻璃陶瓷GC-9材料在本发明的失效概率计算模型下获得的失效概率与现有技术的基于应力的失效概率计算模型下获得的失效概率的对比曲线图。
图6为实施例2中650℃下陶瓷材料YSZ等效蠕变应变和米塞斯应力(Mises应力)随蠕变时间变化的关系曲线。
图7为实施例2中650℃下陶瓷材料YSZ在本发明的失效概率计算模型下获得的失效概率与现有技术的基于应力的失效概率计算模型下获得的失效概率的对比曲线图。
具体实施方式
下面结合附图和实施例，对本发明进行详细说明。
实施例1：
预测玻璃陶瓷GC-9材料在600℃下蠕变50000h的失效概率。
玻璃陶瓷GC-9材料在600℃下蠕变50000h的失效概率预测过程，按照如图1所示的流程进行。
实施例2：
预测陶瓷材料YSZ在650℃下蠕变50000h的失效概率。
陶瓷材料YSZ在650℃下蠕变50000h的失效概率预测过程，按照如图1所示的流程进行。
实施例1和实施例2在计算过程所用参数见表1所示：
表1
威布尔理论需要考虑最弱链假设，即结构在恒定的单轴载荷下，认为它类似于拉伸的N链，每条链都有不同的失效强度，当最弱链失效时，整个结构失效。因此，链的强度与最弱链相关。每条链的失效强度不同，取决于试样内部的缺陷不同，即“体积效应”。
图2为本发明的试样体积大小与缺陷大小的关系示意图，如图2所示，试样体积越大，内部缺陷越大，对应的产生较大的应力强度。
图3为本发明的单轴蠕变失效应变累积分布直方图，如图3所示，在相同的试验条件下，对20组体积为V0的试样在相同应力水平下进行单轴蠕变断裂试验，记录每个断裂蠕变应变值，并以蠕变断裂应变为横坐标，在某一个蠕变断裂应变区间的断裂试样数量为纵坐标，所绘制出的单轴蠕变失效应变值累积分布直方图。
图4为实施例1中等效蠕变应变和米塞斯应力(Mises应力)随蠕变时间的变化曲线；图5为600℃下玻璃陶瓷GC-9材料在本发明的失效概率计算模型下获得的失效概率与现有技术的基于应力的失效概率计算模型下获得的失效概率的对比曲线图。
图6为实施例2中陶瓷材料YSZ等效蠕变应变和米塞斯应力(Mises应力)随蠕变时间的变化曲线；图7为的650℃下陶瓷材料YSZ在本发明的失效概率计算模型下获得的失效概率与现有技术的基于应力的失效概率计算模型下获得的失效概率的对比曲线图。
从图5和图7中可以看出，采用本发明中所提出的基于应变的失效概率计算模型的失效概率随时间的增加而增大，这与工程实际相符合。因为，脆性材料结构在高温长时间服役过程中，其蠕变变形(如图4和图6所示)和损伤逐渐增大，破坏的可能性也逐渐增加，结构可靠性性越来越差，故失效概率在逐渐增加。
而现有技术的基于应力的失效概率计算模型所得到的失效概率随时间的增大而减小，因为蠕变过程中会有应力松弛现象，应力逐渐减小(如图4和图6所示)，所以失效概率逐渐减小。但这与实际不符，因此，现有技术的基于应力的失效概率计算模型不能用于计算高温蠕变状态下的失效概率。
图5和图7中的对比结果进一步证明了上述结论。
图5和图7中的对比结果清楚地表明：本发明的基于应变的失效概率计算模型所获得的脆性材料在高温蠕变状态下的失效概率的预测结果，相对于现有技术，更真实、准确、合理和可靠。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。
Claims (2)
1.一种脆性材料在高温蠕变状态下失效概率的预测方法，其特征在于，包括以下步骤：
第一步，根据脆性材料内部缺陷随机分布的自然属性，假定反应脆性材料属性的单轴蠕变失效应变εf服从威布尔分布；则单轴蠕变失效应变的概率密度函数f(εf)满足下式(1)：
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>&beta;</mi>
<mi>&eta;</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>&beta;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>&beta;</mi>
</msup>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
上式(1)中：
η为变量的尺度参数，η＞0；
β为变量的形状参数，β＞0；
第二步，根据下式(2)所示的单轴与多轴蠕变失效应变εf *的转化关系，根据数学转换关系，得到如下式(3)所示的多轴蠕变失效应变的概率密度分布函数f(εf *)：
<mrow>
<msup>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mo>*</mo>
</msup>
<mo>=</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mfrac>
<mn>2</mn>
<mn>3</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>0.5</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>/</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>0.5</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mi>m</mi>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
上式(2)中：
σm是指材料所承受的静水应力；
σeq为米塞斯应力；
n表示蠕变指数；
为与单轴蠕变失效应变无关的系数；得出多轴蠕变失效应变εf *服从威布尔分布，多轴蠕变失效应变的概率密度分布函数的数学表达式(3)为：
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mi>&beta;</mi>
<mi>&eta;</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mo>*</mo>
</msup>
</mrow>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>&beta;</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mo>*</mo>
</msup>
</mrow>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>&beta;</mi>
</msup>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
第三步，依据结构失效的条件为等效蠕变应变值εe大于多轴蠕变失效应变值εf *的原则，对多轴蠕变失效应变的概率密度分布函数的数学表达式(3)进行积分，即得到如下式(4)所示的失效概率的计算表达式：
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>F</mi>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<msub>
<mi>&epsiv;</mi>
<mi>e</mi>
</msub>
</msubsup>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&epsiv;</mi>
<mi>f</mi>
</msub>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<msup>
<msub>
<mi>d&epsiv;</mi>
<mi>f</mi>
</msub>
<mo>*</mo>
</msup>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mi>e</mi>
</msub>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>&beta;</mi>
</msup>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
在此基础上，考虑到材料内部缺陷的不同，对于体积为V的脆性材料试样，考虑到体积效应，相应的失效概率表达式为下式(5)：
<mrow>
<mi>P</mi>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&epsiv;</mi>
<mi>e</mi>
</msub>
<mi>&eta;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mi>&beta;</mi>
</msup>
<mfrac>
<mi>V</mi>
<msub>
<mi>V</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
上式(5)中：
V0为特征体积；
第四步，在相同的试验条件下，对若干组体积为V0的试样在相同应力水平下进行单轴蠕变断裂试验，记录每个断裂蠕变应变值，并以蠕变断裂应变为横坐标，在某一个蠕变断裂应变区间的断裂试样数量为纵坐标，绘制出单轴蠕变失效应变值累积分布直方图；
第五步，根据绘制出的单轴蠕变失效应变值累积分布直方图，用每个蠕变断裂应变区间断裂试样的数量除以总的断裂试样数量，即为体积V0的试样在该区间内的断裂概率值PF0，将V0和PF0带入上述的失效概率计算公式(4)并两边取两次对数，得到：
ln[-ln(1-PF0)]＝βlnεe-lnηβ (6)
根据各试样在相同应力水平下进行单轴蠕变断裂的试验结果，做出ln[-ln(1-PF0)]与lnεe的曲线，并进行线性回归，所得到的直线的斜率即为参数β，根据所得到的直线与y轴的截距得到参数η；
第六步，根据上式(5)，结合蠕变-损伤本构方程，利用Fortran语言，编写子程序并嵌入到有限元软件ABAQUS中，即得到脆性材料在高温蠕变状态下失效概率的预测结果；
其中，蠕变-损伤本构方程如下所示：
<mrow>
<msub>
<mover>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>3</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>B&sigma;</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msub>
<mi>S</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mi>I</mi>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</mfrac>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&beta;</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mi>&rho;</mi>
</mrow>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>n</mi>
<mo>+</mo>
<mn>3</mn>
<mo>)</mo>
<msup>
<mi>&rho;</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mi>n</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>3</mn>
<mo>)</mo>
<msup>
<mi>&rho;</mi>
<mn>3</mn>
</msup>
</mrow>
<mrow>
<mn>9</mn>
<mi>n</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>3</mn>
<mo>)</mo>
<msup>
<mi>&rho;</mi>
<mn>4</mn>
</msup>
</mrow>
<mrow>
<mn>108</mn>
<mi>n</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mn>4</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&pi;</mi>
<msqrt>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mn>3</mn>
<mo>/</mo>
<mi>n</mi>
</mrow>
</msqrt>
</mrow>
</mfrac>
<msup>
<mi>&omega;</mi>
<mrow>
<mn>3</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&omega;</mi>
<mo>=</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<mfrac>
<msub>
<mover>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>e</mi>
</msub>
<msubsup>
<mi>&epsiv;</mi>
<mi>f</mi>
<mo>*</mo>
</msubsup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，为蠕变应变，σI为最大主应力，B为蠕变第二阶段的常数，β0是与应力相关的函数，ρ是微裂纹损伤参数，ω为蠕变损伤量。
2.根据权利要求1所述的一种脆性材料在高温蠕变状态下失效概率的预测方法，其特征在于，第四步中所述的若干组为10～20组。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710699545.XA CN107273649B (zh) | 2017-08-16 | 2017-08-16 | 一种脆性材料在高温蠕变状态下失效概率的预测方法 |
| US16/476,004 US20190354889A1 (en) | 2017-08-16 | 2017-12-20 | Method of predicting failure probability of brittle material in high temperature creep state |
| PCT/CN2017/117405 WO2019033668A1 (zh) | 2017-08-16 | 2017-12-20 | 一种脆性材料在高温蠕变状态下失效概率的预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710699545.XA CN107273649B (zh) | 2017-08-16 | 2017-08-16 | 一种脆性材料在高温蠕变状态下失效概率的预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107273649A CN107273649A (zh) | 2017-10-20 |
| CN107273649B true CN107273649B (zh) | 2018-05-04 |
Family
ID=60077282
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710699545.XA Active CN107273649B (zh) | 2017-08-16 | 2017-08-16 | 一种脆性材料在高温蠕变状态下失效概率的预测方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US20190354889A1 (zh) |
| CN (1) | CN107273649B (zh) |
| WO (1) | WO2019033668A1 (zh) |
Families Citing this family (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109060549B (zh) * | 2018-10-15 | 2019-07-16 | 西安建筑科技大学 | 基于s形函数建立岩石单轴压缩全过程蠕变模型的方法 |
| CN110991109B (zh) * | 2019-11-22 | 2023-04-21 | 西安航天动力技术研究所 | 一种适用于柔性接头摆动密封可靠性分析方法 |
| CN111177889B (zh) * | 2019-12-10 | 2021-11-23 | 北京航空航天大学 | 一种基于弹用电点火头的密封火工品贮存寿命评估方法 |
| CN111144016B (zh) * | 2019-12-30 | 2021-03-02 | 河海大学 | 一种多轴应力状态混凝土变四参数损伤本构 |
| CN111832124B (zh) * | 2020-05-28 | 2022-05-31 | 西北工业大学 | 元模型重要抽样结合空间分割的涡轮叶片重要性分析方法 |
| CN111859686B (zh) * | 2020-07-27 | 2023-10-10 | 北京卫星环境工程研究所 | 一种基于退化过程模型的机械构件寿命预测方法 |
| CN112100806B (zh) * | 2020-07-31 | 2023-12-15 | 东南大学 | 一种表征复合材料本构曲线的获取方法 |
| CN112630044B (zh) * | 2020-11-19 | 2024-06-18 | 西北工业大学 | 基于晶体取向的镍基单晶合金的蠕变寿命预测方法 |
| CN112730092B (zh) * | 2020-12-24 | 2022-07-12 | 南京航空航天大学 | 一种基于失效模式一致的蠕变载荷等效加速方法 |
| CN113076627A (zh) * | 2021-03-17 | 2021-07-06 | 湖南大学 | 一种镍基高温合金宽温域蠕变强度预测方法 |
| CN113111410B (zh) * | 2021-03-19 | 2023-04-28 | 西安建筑科技大学 | 混杂纤维混凝土动态劈拉强度标度律模型的建立方法 |
| CN113128026B (zh) * | 2021-03-22 | 2022-06-24 | 长沙理工大学 | 再生沥青混合料高温蠕变和低温松弛性能平衡设计方法 |
| CN113514343B (zh) * | 2021-07-13 | 2022-05-06 | 华东理工大学 | 一种多轴蠕变性能参数的测试方法 |
| CN114329930A (zh) * | 2021-12-15 | 2022-04-12 | 华东理工大学 | 一种基于骨点应力概念的断裂寿命确定方法和系统 |
| CN114091315B (zh) * | 2022-01-24 | 2022-04-15 | 西北工业大学 | 高温合金应力及损伤演化方法、装置、存储介质及电子设备 |
| CN114547753B (zh) * | 2022-03-03 | 2024-03-22 | 中交基础设施养护集团有限公司 | 考虑地层空间变异性的长条形地下结构分段长度确定方法 |
| CN114692468B (zh) * | 2022-04-24 | 2023-09-19 | 西南交通大学 | 一种连续纤维增强陶瓷基复合材料跨尺度损伤预测方法 |
| CN115017641B (zh) * | 2022-05-31 | 2023-04-07 | 西安交通大学 | 基于大变形理论的反应堆压力容器失效评估方法及系统 |
| CN115032156B (zh) * | 2022-06-07 | 2024-07-26 | 安徽农业大学 | 一种热处理木材脆性的预测方法及预测模型 |
| CN117933471B (zh) * | 2024-01-24 | 2024-07-05 | 北京建筑大学 | 一种浸水脆性固体材料的长期压缩寿命预测评价方法 |
| CN118446094B (zh) * | 2024-05-07 | 2025-03-11 | 西南石油大学 | 一种基于时空贝叶斯的埋地管道失效预测方法 |
| CN119203452B (zh) * | 2024-09-19 | 2025-05-20 | 南京林业大学 | 一种基于改进Kelvin-Voight模型的重组竹材压缩蠕变应变预测方法 |
| CN120012445B (zh) * | 2025-04-17 | 2025-06-24 | 中国科学院力学研究所 | 深海耐压圆柱壳结构破坏尺度预测方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0421229B1 (en) * | 1989-10-04 | 1995-03-01 | General Electric Company | Creep, stress rupture and hold-time fatigue crack resistant alloys |
| JP2003114177A (ja) * | 2001-10-02 | 2003-04-18 | Mitsubishi Heavy Ind Ltd | 熱サイクルクリープ疲労の試験片、その試験装置、及び、その試験方法 |
| CN101710053A (zh) * | 2009-11-06 | 2010-05-19 | 上海师范大学 | 高温材料的蠕变寿命预测方法 |
| CN101994908A (zh) * | 2010-08-12 | 2011-03-30 | 华东理工大学 | 高温管线系统实现可靠性维修规划的方法 |
| CN103323343A (zh) * | 2013-06-17 | 2013-09-25 | 四川大学 | 聚合物材料蠕变失效寿命的判定方法及预测方法 |
| CN106557630A (zh) * | 2016-11-21 | 2017-04-05 | 中国石油大学(华东) | 一种材料在多轴应力状态下的蠕变‑损伤寿命预测方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7979237B2 (en) * | 2008-05-07 | 2011-07-12 | A. M. Fitzgerald & Associates, Llc | Fracture prediction for crystalline microstructures |
| US20170124448A1 (en) * | 2015-10-30 | 2017-05-04 | Northrop Grumman Systems Corporation | Concurrent uncertainty management system |
-
2017
- 2017-08-16 CN CN201710699545.XA patent/CN107273649B/zh active Active
- 2017-12-20 WO PCT/CN2017/117405 patent/WO2019033668A1/zh active Application Filing
- 2017-12-20 US US16/476,004 patent/US20190354889A1/en not_active Abandoned
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0421229B1 (en) * | 1989-10-04 | 1995-03-01 | General Electric Company | Creep, stress rupture and hold-time fatigue crack resistant alloys |
| JP2003114177A (ja) * | 2001-10-02 | 2003-04-18 | Mitsubishi Heavy Ind Ltd | 熱サイクルクリープ疲労の試験片、その試験装置、及び、その試験方法 |
| CN101710053A (zh) * | 2009-11-06 | 2010-05-19 | 上海师范大学 | 高温材料的蠕变寿命预测方法 |
| CN101994908A (zh) * | 2010-08-12 | 2011-03-30 | 华东理工大学 | 高温管线系统实现可靠性维修规划的方法 |
| CN103323343A (zh) * | 2013-06-17 | 2013-09-25 | 四川大学 | 聚合物材料蠕变失效寿命的判定方法及预测方法 |
| CN106557630A (zh) * | 2016-11-21 | 2017-04-05 | 中国石油大学(华东) | 一种材料在多轴应力状态下的蠕变‑损伤寿命预测方法 |
Non-Patent Citations (4)
| Title |
|---|
| Effect of Temperature Fluctuation on Creep and Failure Probability for Planar Solid Oxide Fuel Cell;Wenchun Jiang 等;《Journal of Fuel Cell Science and Technology》;20151031;第12卷(第5期);第051004-1--051004-10页 * |
| Using short-time creep relaxation effect to decrease the residual stress in the bonded compliant seal of planar solid oxide fuel cell e A-finite element simulation;Yu-Cai Zhang 等;《Journal of Power Sources》;20140103;第255卷;第108-115页 * |
| 基于连续损伤理论的多轴蠕变设计;姚华堂 等;《中国机械工程》;20070630;第18卷(第12期);第1438-1443页 * |
| 高温蠕变与断裂评价的若干关键问题;涂善东 等;《金属学报》;20090731;第45卷(第7期);第781-787页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| US20190354889A1 (en) | 2019-11-21 |
| WO2019033668A1 (zh) | 2019-02-21 |
| CN107273649A (zh) | 2017-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107273649B (zh) | 一种脆性材料在高温蠕变状态下失效概率的预测方法 | |
| Mishnaevsky Jr et al. | Materials of large wind turbine blades: recent results in testing and modeling | |
| CN105117536B (zh) | 一种rpv含裂纹类缺陷的简化弹塑性断裂力学分析方法 | |
| CN104462757B (zh) | 基于监测数据的Weibull分布可靠性序贯验证试验方法 | |
| CN107563054B (zh) | 一种基于SWT参数的Weakest-Link方法的涡轮盘概率寿命分析方法 | |
| CN103778276B (zh) | 基于疲劳寿命分布的复合材料可靠度预测方法 | |
| CN111751199A (zh) | 基于eifs分布的疲劳寿命预测方法 | |
| CN103970999B (zh) | 飞机结构疲劳裂纹安全损伤扩展周期确定方法 | |
| CN106383959A (zh) | 一种基于最大熵模型的材料疲劳寿命的预测方法 | |
| CN103940662B (zh) | 高温材料应力松弛剩余应力和损伤的预测方法 | |
| CN117114240A (zh) | 一种加高重力坝结构的安全评价方法、系统、设备及介质 | |
| Wang et al. | Fatigue reliability updating evaluation of existing steel bridges | |
| CN1945269A (zh) | 高温构件材料的寿命及寿命消耗的预测方法 | |
| Zhao et al. | A multi‐axial low‐cycle fatigue life prediction model considering effects of additional hardening | |
| CN106290775A (zh) | 一种电站锅炉受热面sa210c钢材料状态评估方法 | |
| CN114936498A (zh) | 考虑材料性能退化的蠕变疲劳损伤等级确定方法及系统 | |
| CN106653099B (zh) | 一种聚变堆偏滤器的疲劳寿命分析评判方法 | |
| CN111881564A (zh) | 一种关于机械结构变幅疲劳寿命预测方法 | |
| CN117877646B (zh) | 基于不同评价指标的轮轨材料滚动接触疲劳损伤预测方法 | |
| CN106372275A (zh) | 用于蠕变裂纹扩展数值模拟结果的修正方法 | |
| CN104281750B (zh) | 一种故障率为浴盆曲线的寿命分布方法 | |
| CN117993204A (zh) | 一种长时蠕变寿命预测方法、系统及电子设备 | |
| CN106933780A (zh) | 一种风力发电机叶片疲劳寿命的计算方法 | |
| Wei et al. | Fatigue life assessment of welded structures with the linear traction stress analysis approach | |
| CN101482464A (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 |