CN109359330A - 一种描述高温材料蠕变变形三阶段的方法及模型 - Google Patents
一种描述高温材料蠕变变形三阶段的方法及模型 Download PDFInfo
- Publication number
- CN109359330A CN109359330A CN201811032610.4A CN201811032610A CN109359330A CN 109359330 A CN109359330 A CN 109359330A CN 201811032610 A CN201811032610 A CN 201811032610A CN 109359330 A CN109359330 A CN 109359330A
- Authority
- CN
- China
- Prior art keywords
- creep
- exp
- stress
- temperature
- test
- 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
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]
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
本发明公开了一种描述高温材料蠕变变形三阶段的方法及模型，该方法的步骤为：(1)通过试验获得材料的工程蠕变变形‑时间曲线，并将其转换为真实蠕变应变‑时间曲线；(2)将真实蠕变应变‑时间曲线进行归一化；(3)将蠕变模型中的参数k、η、α表示为应力和温度的函数；(4)利用试验数据拟合蠕变模型参数；(5)将蠕变模型嵌入有限元软件，实现对实际结构的蠕变分析。本发明通过参数k、η、α能够很好的描述蠕变的第一阶段和第二阶段材料内部位错的运动变化导致的蠕变速率逐渐减小和趋于稳定的现象，此外通过断裂寿命归一化后的时间项，能够反应蠕变空洞、裂纹等微观损伤从而表征蠕变第三阶段蠕变速率迅速增大的现象。
Description
技术领域
本发明涉及一种能够完整描述高温材料蠕变变形三阶段的方法及模型，属于高温结构强度技术领域。
背景技术
蠕变通常是指在温度、载荷/应力不变的条件下，变形随时间的增加而增大的现象。其主要特点在于当结构的温度达到一定时，即使其内部的应力远小于材料的屈服强度仍然能产生明显的蠕变。航空发动机热端部件特别是转动件承受着高温和大应力，在使用中会发生明显的蠕变现象。
目前常用的蠕变模型有时间硬化模型和应变硬化模型等，但这类模型只能描述蠕变第一阶段或者前两个阶段。然而，在实际的蠕变问题中第三阶段的蠕变变形占据了整个寿命期内蠕变变形的70％以上，在应力水平较小时甚至占到总变形的80％以上，且第三阶段寿命占到总寿命的50％左右，因此，对蠕变第三阶段的变形描述是分析高温结构/材料失效过程、预测构件变形或寿命的基础。1985年Evans和Wilshire在其著作中提出一种能描述完整的蠕变变形的模型——θ映射法，虽然θ映射法可以模拟蠕变变形的3个阶段，但其模型中的参数并不具有很好的物理基础，此外，在2015年由王延荣等人提出的归一化参数蠕变模型与θ映射法类似，均能描述蠕变三个阶段，但是蠕变参数缺少物理意义。而对基于纯蠕变物理过程建立的模型而言，确定参数的方法太过繁琐，且计算过程相当复杂，不利于工程应用。因此，有必要发展一种既能反应蠕变物理过程又能够更加便于应用的蠕变模型。
发明内容
为发展一种既能更好反应蠕变物理过程又能便于工程应用的蠕变模型，解决目前商用有限元软件中自带模型不能描述蠕变三阶段的问题，使得模型既能描述蠕变变形整体过程又能反应蠕变的微观物理过程，本发明的目的是提供一种描述高温材料蠕变变形三阶段的方法及模型。
为实现上述目的，本发明采用的技术方案为：
一种描述高温材料蠕变变形三阶段的方法，包括以下步骤：
(1)通过试验获得材料的工程蠕变变形-时间曲线，并将其转换为真实蠕变应变-时间曲线；
(2)将步骤(1)得到的真实蠕变应变-时间曲线进行归一化；
(3)将蠕变模型中的参数k、η、α表示为应力和温度的函数；其中，上式中：εc为蠕变应变，σ为应力，t、tf分别为蠕变时间和材料在指定温度、应力下的断裂寿命；
(4)利用试验数据拟合蠕变模型参数；
(5)将蠕变模型嵌入有限元软件，实现对实际结构的蠕变分析。
所述步骤(1)中，利用体积不变假设，将试验获得的工程应力和工程蠕变应变转换为真实应力和真实蠕变应变，并与试验时间一起构成蠕变模型参数拟合的初始数据；真实应力和真实蠕变应变转换公式为：
式中：σT、σE分别为真实蠕变应变、真实应力、工程蠕变应变以及工程应力，εe,εp分别为工程弹性应变和工程塑性应变。
所述步骤(2)中，分别运用各自试验载荷、温度下的断裂寿命将蠕变试验时间除以断裂寿命，获得归一化时间、真实应力以及真实蠕变应变三列数据，作为模型参数拟合的输入数据。
所述步骤(3)中，参数k、η、α的温度和应力函数表达式为：
k＝c1+c2σ+c3T+c4σT，式中σ为施加应力，T为试验温度，c1-c4为试验拟合常数；
η＝c5+c6σ+c7T+c8σT，式中σ为施加应力，T为试验温度，c5-c8为试验拟合常数；
α＝c9+c10σ+c11T+c12σT，式中σ为施加应力，T为试验温度，c9-c12为试验拟合常数；
当各试验条件下的温度相同时上述表达式简化为k＝c1+c2σ、η＝c3+c4σ、α＝c5+c6σ。
所述步骤(4)中，利用1stOpt软件中的参数拟合功能，采用麦夸特法和通用全局优化法，拟合得到参数ci的值，下标i为1至12，当各试验温度相同时得到简化后的参数ci的值，下标i为1至6。
所述步骤(5)中，编写有限元蠕变子程序usercreep时需要提供蠕变应变增量delcr、蠕变应变增量对等效应力的导数dcrda(1)和蠕变应变增量对蠕变应变的导数dcrda(2)；
若各试验温度相同，则其表达式分别为：
delcr＝((σ·exp(-k/η·t/tf)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)-(σ·exp(-k/η·t/tf)·k·(exp((α+k/η)·t/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α)))·Δt，
dcrda(1)＝((exp(-k/η·t/tf)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)+(σ·exp(-k/η·t/tf)·((c6·(t/tf)^(α-1))/tf+(exp((t·(α+k/η))/tf)·(c6+c2/η-(c4·k)/(η)^2))/tf-(c6·(t/tf)^(1/α))/(tf·(α)^2)+(c6·log(t/tf)·(t/tf)^(α-1)·α)/tf+(t·exp((t·(α+k/η))/tf)·(α+k/η)·(c6+c2/η-(c4·k)/(η)^2))/tf^2-(c6·log(t/tf)·(t/tf)^(1/α)·(1/α+1))/(tf·(α)^2)))/(k+η·α)-(σ·exp(-(t·k)/(tf·η))·((c2·t)/(tf·η)-(c4·t·k)/(tf·(η)^2))·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)-(σ·exp(-(t·k)/(tf·η))·(c2+c6·η+c4·α)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)^2-(exp(-(t·k)/(tf·η))·k·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α))-(c2·σ·exp(-(t·k)/(tf·η))·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α))-(σ·exp(-(t·k)/(tf·η))·k·(c6·log(t/tf)·(t/tf)^α+(t·exp((t·(α+k/η))/tf)·(c6+c2/η-(c4·k)/(η)^2))/tf-(c6·log(t/tf)·(t/tf)^(1/α+1))/(α)^2))/(tf·η·(k+η·α))+(σ·exp(-k/η·t/tf)·k·(c2+c6·η+c4·α)·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α)^2)+(c4·σ·exp(-k/η·t/tf)·a1·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·(η)^2·(k+η·α))+(σ·exp(-(t·k)/(tf·η))·((c2·t)/(tf·η)-(c4·t·k)/(tf·(η)^2))·k·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α)))·Δt，
dcrda(2)＝((σ·exp(-(t·k)/(tf·η))·((exp((t·(α+k/η))/tf)·(α+k/η)^2)/tf^2+((t/tf)^(1/α-1)·(1/α+1))/(tf^2·α)+((t/tf)^(α-2)·α·(α-1))/tf^2))/(k+η·α)+(σ·exp(-(t·k)/(tf·η))·k^2·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf^2·η^2·(k+η·α))-(2·σ·exp(-(t·k)/(tf·η))·k·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(tf·η·(k+η·α)))·(1/((σ·exp(-k/η·t/tf)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)-(σ·exp(-k/η·t/tf)·k·(exp((α+k/η)·t/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α))))·Δt
上式中，Δt为软件计算中的时间增量。
一种描述高温材料蠕变变形三阶段的模型，该模型的表达式为：
式中：εc为蠕变应变，σ为应力，t、tf分别为蠕变时间和材料在指定温度、应力下的断裂寿命，k、η、α均表示为应力和温度的函数，为反应材料内部蠕变过程中位错变化的参数。
有益效果：本发明提出的方法及模型既反应了由材料内部位错变化引起的蠕变第一阶段和第二阶段的变形特征，又包含了蠕变空洞等损伤引起的第三阶段应变速率随时间增加迅速增大的现象，因此能够更加准确的描述材料的蠕变变形全过程。通过编写有限元子程序，可以将本发明的模型应用于实际结构的蠕变分析中。
附图说明
图1为本发明的实施流程图；
图2为典型蠕变曲线及其三个阶段示意图；
图3为TC11材料500℃下4种应力条件的工程蠕变应变图；
图4为TC11材料500℃下4种应力条件的真实蠕变应变图；
图5为TC11材料500℃下4种应力条件归一化时间的真实蠕变应变图；
图6为TC11材料500℃下4种应力条件试验的归一化蠕变曲线与模型拟合曲线对比图；
图7为TC11材料500℃下4种应力条件试验的蠕变曲线与模型拟合曲线对比图。
具体实施方式
本发明提供一种描述高温材料蠕变变形三阶段的方法及模型，包括能够反应材料内部因应变硬化而导致蠕变速率逐渐下降的蠕变第一阶段、硬化过程与恢复过程达到平衡而产生的蠕变速率趋于恒定的蠕变第二阶段以及因材料内部产生空洞和裂纹而导致蠕变速率快速增大的蠕变第三阶段。并且通过编写有限元软件ANSYS的usercreep子程序，将该方法应用于实际结构的蠕变变形分析中。蠕变应变的表达式为：
式中：εc为蠕变应变，σ为应力，t、tf分别为蠕变时间和材料在指定温度、应力下的断裂寿命，k、η、α均表示为应力和温度的函数，该模型的主要特点在于：k、η、α为反应材料内部蠕变过程中位错变化的参数，能够很好的描述蠕变的第一阶段和第二阶段。此外，通过断裂寿命对时间进行归一化，引入第三阶段的蠕变损伤来描述第三阶段的变形规律，即式(1)中等号右边中括号里的第二项。同时式(1)中等号右边中括号里的第三项对蠕变变形的第一阶段进行修正，最终达到能够更好描述材料蠕变变形规律的目的。
下面结合实际应用实例和附图对本发明作更进一步的说明：
本发明涉及的一种描述高温材料蠕变变形三阶段的方法，其实施流程如附图1所示，典型的材料蠕变曲线如附图2所示，其在TC11钛合金材料500℃下的蠕变变形分析中的运用包括以下步骤：
(1)通过试验获得TC11材料500℃下的工程蠕变应变-时间曲线，如附图3所示，将其转换为真实蠕变应变-时间曲线，如附图4所示；
光滑试件蠕变试验获得材料的工程蠕变-时间曲线(εc-t)，由如下变换公式转化为真实蠕变-时间曲线和真实应力；
式中：σT、σE分别为真实蠕变应变、真实应力、工程蠕变应变以及工程应力，εe,εp分别为工程弹性应变和工程塑性应变；
(2)将TC11材料500℃下的真实蠕变应变-时间曲线进行归一化：
分别运用各自试验载荷、温度下的断裂寿命tf将时间进行归一化，TC11材料在500℃、4种工程应力537MPa、558MPa、580MPa和590MPa作用下的断裂寿命tf分别为645.256h、392.444h、313.728h和215.203h，由此获得归一化时间、真实应力以及真实蠕变应变三列数据，作为模型参数拟合的输入数据，归一化时间的真实蠕变应变曲线如附图5所示；
(3)将参数k、η、α表示为应力和温度的函数；
参数k、η、α的温度和应力函数表达式为：
k＝c1+c2σ+c3T+c4σT，式中σ为施加应力，T为试验温度，c1-c4为试验拟合常数；
η＝c5+c6σ+c7T+c8σT，式中σ为施加应力，T为试验温度，c5-c8为试验拟合常数；
α＝c9+c10σ+c11T+c12σT，式中σ为施加应力，T为试验温度，c9-c12为试验拟合常数；
当各试验条件下的温度相同时上述表达式简化为k＝c1+c2σ、η＝c3+c4σ、α＝c5+c6σ。
(4)利用试验数据拟合蠕变模型参数，拟合结果与试验数据对比如附图6和图7所示；
利用1stOpt软件中的参数拟合功能，采用麦夸特法+通用全局优化法，拟合得到参数ci的值，下标i为1至12，当各试验温度相同时可得到简化后的参数ci的值，下标i为1至6，拟合TC11材料在500℃、5种应力下的蠕变试验数据得到的参数c1～c6分别为：
c1＝267782.6665,
c2＝-349.6828，
c3＝48451.6357，
c4＝-175.1019，
c5＝0.55695，
c6＝-0.00044296。
(5)将该模型嵌入有限元软件，实现对实际结构的蠕变分析：
编写有限元蠕变子程序usercreep时需要提供蠕变应变增量delcr、蠕变应变增量对等效应力的导数dcrda(1)和蠕变应变增量对蠕变应变的导数dcrda(2)；
以各试验温度相同的条件为例，其表达式分别为：
delcr＝((σ·exp(-k/η·t/tf)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)-(σ·exp(-k/η·t/tf)·k·(exp((α+k/η)·t/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α)))·Δt，
dcrda(1)＝((exp(-k/η·t/tf)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)+(σ·exp(-k/η·t/tf)·((c6·(t/tf)^(α-1))/tf+(exp((t·(α+k/η))/tf)·(c6+c2/η-(c4·k)/(η)^2))/tf-(c6·(t/tf)^(1/α))/(tf·(α)^2)+(c6·log(t/tf)·(t/tf)^(α-1)·α)/tf+(t·exp((t·(α+k/η))/tf)·(α+k/η)·(c6+c2/η-(c4·k)/(η)^2))/tf^2-(c6·log(t/tf)·(t/tf)^(1/α)·(1/α+1))/(tf·(α)^2)))/(k+η·α)-(σ·exp(-(t·k)/(tf·η))·((c2·t)/(tf·η)-(c4·t·k)/(tf·(η)^2))·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)-(σ·exp(-(t·k)/(tf·η))·(c2+c6·η+c4·α)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)^2-(exp(-(t·k)/(tf·η))·k·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α))-(c2·σ·exp(-(t·k)/(tf·η))·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α))-(σ·exp(-(t·k)/(tf·η))·k·(c6·log(t/tf)·(t/tf)^α+(t·exp((t·(α+k/η))/tf)·(c6+c2/η-(c4·k)/(η)^2))/tf-(c6·log(t/tf)·(t/tf)^(1/α+1))/(α)^2))/(tf·η·(k+η·α))+(σ·exp(-k/η·t/tf)·k·(c2+c6·η+c4·α)·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α)^2)+(c4·σ·exp(-k/η·t/tf)·a1·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·(η)^2·(k+η·α))+(σ·exp(-(t·k)/(tf·η))·((c2·t)/(tf·η)-(c4·t·k)/(tf·(η)^2))·k·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α)))·Δt，
dcrda(2)＝((σ·exp(-(t·k)/(tf·η))·((exp((t·(α+k/η))/tf)·(α+k/η)^2)/tf^2+((t/tf)^(1/α-1)·(1/α+1))/(tf^2·α)+((t/tf)^(α-2)·α·(α-1))/tf^2))/(k+η·α)+(σ·exp(-(t·k)/(tf·η))·k^2·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf^2·η^2·(k+η·α))-(2·σ·exp(-(t·k)/(tf·η))·k·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(tf·η·(k+η·α)))·(1/((σ·exp(-k/η·t/tf)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)-(σ·exp(-k/η·t/tf)·k·(exp((α+k/η)·t/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α))))·Δt。
上式中，Δt为软件计算中的时间增量。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (7)
1.一种描述高温材料蠕变变形三阶段的方法，其特征在于：包括以下步骤：
(1)通过试验获得材料的工程蠕变变形-时间曲线，并将其转换为真实蠕变应变-时间曲线；
(2)将步骤(1)得到的真实蠕变应变-时间曲线进行归一化；
(3)将蠕变模型中的参数k、η、α表示为应力和温度的函数；其中，上式中：εc为蠕变应变，σ为应力，t、tf分别为蠕变时间和材料在指定温度、应力下的断裂寿命；
(4)利用试验数据拟合蠕变模型参数；
(5)将蠕变模型嵌入有限元软件，实现对实际结构的蠕变分析。
2.根据权利要求1所述的描述高温材料蠕变变形三阶段的方法，其特征在于：所述步骤(1)中，利用体积不变假设，将试验获得的工程应力和工程蠕变应变转换为真实应力和真实蠕变应变，并与试验时间一起构成蠕变模型参数拟合的初始数据；真实应力和真实蠕变应变转换公式为：
式中：σT、σE分别为真实蠕变应变、真实应力、工程蠕变应变以及工程应力，εe,εp分别为工程弹性应变和工程塑性应变。
3.根据权利要求1所述的描述高温材料蠕变变形三阶段的方法，其特征在于：所述步骤(2)中，分别运用各自试验载荷、温度下的断裂寿命将蠕变试验时间除以断裂寿命，获得归一化时间、真实应力以及真实蠕变应变三列数据，作为模型参数拟合的输入数据。
4.根据权利要求1所述的描述高温材料蠕变变形三阶段的方法，其特征在于：所述步骤(3)中，参数k、η、α的温度和应力函数表达式为：
k＝c1+c2σ+c3T+c4σT，式中σ为施加应力，T为试验温度，c1-c4为试验拟合常数；
η＝c5+c6σ+c7T+c8σT，式中σ为施加应力，T为试验温度，c5-c8为试验拟合常数；
α＝c9+c10σ+c11T+c12σT，式中σ为施加应力，T为试验温度，c9-c12为试验拟合常数；
当各试验条件下的温度相同时上述表达式简化为k＝c1+c2σ、η＝c3+c4σ、α＝c5+c6σ。
5.根据权利要求1所述的描述高温材料蠕变变形三阶段的方法，其特征在于：所述步骤(4)中，利用1stOpt软件中的参数拟合功能，采用麦夸特法和通用全局优化法，拟合得到参数ci的值，下标i为1至12，当各试验温度相同时得到简化后的参数ci的值，下标i为1至6。
6.根据权利要求1所述的描述高温材料蠕变变形三阶段的方法，其特征在于：所述步骤(5)中，编写有限元蠕变子程序usercreep时需要提供蠕变应变增量delcr、蠕变应变增量对等效应力的导数dcrda(1)和蠕变应变增量对蠕变应变的导数dcrda(2)；
若各试验温度相同，则其表达式分别为：
delcr＝((σ·exp(-k/η·t/tf)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)-(σ·exp(-k/η·t/tf)·k·(exp((α+k/η)·t/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α)))·Δt，
dcrda(1)＝((exp(-k/η·t/tf)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)+(σ·exp(-k/η·t/tf)·((c6·(t/tf)^(α-1))/tf+(exp((t·(α+k/η))/tf)·(c6+c2/η-(c4·k)/(η)^2))/tf-(c6·(t/tf)^(1/α))/(tf·(α)^2)+(c6·log(t/tf)·(t/tf)^(α-1)·α)/tf+(t·exp((t·(α+k/η))/tf)·(α+k/η)·(c6+c2/η-(c4·k)/(η)^2))/tf^2-(c6·log(t/tf)·(t/tf)^(1/α)·(1/α+1))/(tf·(α)^2)))/(k+η·α)-(σ·exp(-(t·k)/(tf·η))·((c2·t)/(tf·η)-(c4·t·k)/(tf·(η)^2))·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)-(σ·exp(-(t·k)/(tf·η))·(c2+c6·η+c4·α)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)^2-(exp(-(t·k)/(tf·η))·k·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α))-(c2·σ·exp(-(t·k)/(tf·η))·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α))-(σ·exp(-(t·k)/(tf·η))·k·(c6·log(t/tf)·(t/tf)^α+(t·exp((t·(α+k/η))/tf)·(c6+c2/η-(c4·k)/(η)^2))/tf-(c6·log(t/tf)·(t/tf)^(1/α+1))/(α)^2))/(tf·η·(k+η·α))+(σ·exp(-k/η·t/tf)·k·(c2+c6·η+c4·α)·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α)^2)+(c4·σ·exp(-k/η·t/tf)·a1·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·(η)^2·(k+η·α))+(σ·exp(-(t·k)/(tf·η))·((c2·t)/(tf·η)-(c4·t·k)/(tf·(η)^2))·k·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α)))·Δt，
dcrda(2)＝((σ·exp(-(t·k)/(tf·η))·((exp((t·(α+k/η))/tf)·(α+k/η)^2)/tf^2+((t/tf)^(1/α-1)·(1/α+1))/(tf^2·α)+((t/tf)^(α-2)·α·(α-1))/tf^2))/(k+η·α)+(σ·exp(-(t·k)/(tf·η))·k^2·(exp((t·(α+k/η))/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf^2·η^2·(k+η·α))-(2·σ·exp(-(t·k)/(tf·η))·k·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(tf·η·(k+η·α)))·(1/((σ·exp(-k/η·t/tf)·(((t/tf)^(α-1)·α)/tf+((t/tf)^(1/α)·(1/α+1))/tf+(exp((t·(α+k/η))/tf)·(α+k/η))/tf))/(k+η·α)-(σ·exp(-k/η·t/tf)·k·(exp((α+k/η)·t/tf)+(t/tf)^(1/α+1)+(t/tf)^α-1))/(tf·η·(k+η·α))))·Δt
上式中Δt为软件计算中的时间增量。
7.一种描述高温材料蠕变变形三阶段的模型，其特征在于：该模型的表达式为：
式中：εc为蠕变应变，σ为应力，t、tf分别为蠕变时间和材料在指定温度、应力下的断裂寿命，k、η、α均表示为应力和温度的函数，为反应材料内部蠕变过程中位错变化的参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811032610.4A CN109359330B (zh) | 2018-09-05 | 2018-09-05 | 一种描述高温材料蠕变变形三阶段的方法及模型 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811032610.4A CN109359330B (zh) | 2018-09-05 | 2018-09-05 | 一种描述高温材料蠕变变形三阶段的方法及模型 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109359330A true CN109359330A (zh) | 2019-02-19 |
| CN109359330B CN109359330B (zh) | 2021-02-26 |
Family
ID=65350309
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811032610.4A Active CN109359330B (zh) | 2018-09-05 | 2018-09-05 | 一种描述高温材料蠕变变形三阶段的方法及模型 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109359330B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110688788A (zh) * | 2019-08-28 | 2020-01-14 | 南京航空航天大学 | 一种高温材料蠕变变形和寿命预测方法及模型 |
| CN112926230A (zh) * | 2019-12-06 | 2021-06-08 | 核工业理化工程研究院 | 结构件的材料蠕变性能指标确立方法 |
| CN113125275A (zh) * | 2021-04-06 | 2021-07-16 | 西北工业大学 | 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080204451A1 (en) * | 2007-02-26 | 2008-08-28 | Samsung Electronics Co., Ltd | Geometry processor using a post-vertex cache and method thereof |
| CN104809273A (zh) * | 2015-04-03 | 2015-07-29 | 北京航空航天大学 | 一种描述蠕变变形的方法 |
| CN106529017A (zh) * | 2016-11-08 | 2017-03-22 | 华南理工大学 | 一种高铬钢构件高温蠕变变形预测与蠕变损伤分析方法 |
-
2018
- 2018-09-05 CN CN201811032610.4A patent/CN109359330B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080204451A1 (en) * | 2007-02-26 | 2008-08-28 | Samsung Electronics Co., Ltd | Geometry processor using a post-vertex cache and method thereof |
| CN104809273A (zh) * | 2015-04-03 | 2015-07-29 | 北京航空航天大学 | 一种描述蠕变变形的方法 |
| CN106529017A (zh) * | 2016-11-08 | 2017-03-22 | 华南理工大学 | 一种高铬钢构件高温蠕变变形预测与蠕变损伤分析方法 |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110688788A (zh) * | 2019-08-28 | 2020-01-14 | 南京航空航天大学 | 一种高温材料蠕变变形和寿命预测方法及模型 |
| CN110688788B (zh) * | 2019-08-28 | 2021-06-22 | 南京航空航天大学 | 一种高温材料蠕变变形和寿命预测方法 |
| CN112926230A (zh) * | 2019-12-06 | 2021-06-08 | 核工业理化工程研究院 | 结构件的材料蠕变性能指标确立方法 |
| CN113125275A (zh) * | 2021-04-06 | 2021-07-16 | 西北工业大学 | 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法 |
| CN113125275B (zh) * | 2021-04-06 | 2024-05-03 | 西北工业大学 | 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109359330B (zh) | 2021-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109359330A (zh) | 一种描述高温材料蠕变变形三阶段的方法及模型 | |
| CN109142083B (zh) | 一种变载荷历程下的蠕变损伤计算方法 | |
| WO2021114994A1 (zh) | 一种基于纳米压痕实验的晶体塑性材料参数反演识别的实现方法 | |
| CN110688788B (zh) | 一种高温材料蠕变变形和寿命预测方法 | |
| CN109598079B (zh) | 一种气缸盖分区疲劳寿命预估方法 | |
| Whittaker et al. | An analysis of modern creep lifing methodologies in the titanium alloy Ti6-4 | |
| CN108256179B (zh) | 一种预测材料蠕变曲线的方法 | |
| CN102331377B (zh) | 一种评估t/p92钢的蠕变性能的方法 | |
| CN110220805B (zh) | 一种基于蠕变疲劳损伤累积的变幅多轴热机疲劳寿命预测方法 | |
| CN115017745B (zh) | 一种用于火电机组结构件的蠕变疲劳寿命预测方法及系统 | |
| CN109117604B (zh) | 一种镍基单晶合金基于滑移系的Walker模型材料常数的获取方法 | |
| CN108254250B (zh) | 一种考虑动态应变时效的热机多轴应力应变关系确定方法 | |
| CN111400876B (zh) | 一种基于不同约束条件下的金属构件热机械疲劳寿命预测方法 | |
| CN111222267A (zh) | 一种冲压发动机热端部件寿命分析方法 | |
| CN108182327B (zh) | 一种基于线性损伤累积的多轴热机械疲劳寿命预测方法 | |
| Banaszkiewicz et al. | The effect of strength differential on material effort and lifetime of steam turbine rotors under thermo-mechanical load | |
| CN113125275A (zh) | 镍基单晶高温合金蠕变模型参数确定与蠕变寿命预测方法 | |
| CN112630045A (zh) | 基于真实叶片样品的镍基单晶合金的蠕变寿命预测方法 | |
| Lissenden et al. | Experimental investigation of cyclic and time-dependent deformation of titanium alloy at elevated temperature | |
| CN110411863B (zh) | 一种基于蠕变延性的高温蠕变寿命预测方法 | |
| Molla et al. | Transient deformational properties of high temperature alloys used in solid oxide fuel cell stacks | |
| CN109101750B (zh) | 一种电子器件湿热环境下的加速寿命预测方法 | |
| CN116359060A (zh) | 一种金属材料基于应变速率影响的热机械疲劳寿命预测方法及存储介质 | |
| CN112149334A (zh) | 一种基于反演优化算法的蠕变本构参数获取方法 | |
| CN110967245B (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 |