CN109214104A - 一种预测变应力下铝铜合金蠕变时效行为的方法 - Google Patents
一种预测变应力下铝铜合金蠕变时效行为的方法 Download PDFInfo
- Publication number
- CN109214104A CN109214104A CN201811123118.8A CN201811123118A CN109214104A CN 109214104 A CN109214104 A CN 109214104A CN 201811123118 A CN201811123118 A CN 201811123118A CN 109214104 A CN109214104 A CN 109214104A
- Authority
- CN
- China
- Prior art keywords
- stress
- creep
- dislocation
- equation
- theta
- 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
- 230000035882 stress Effects 0.000 title claims abstract description 102
- 230000032683 aging Effects 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 45
- JRBRVDCKNXZZGH-UHFFFAOYSA-N alumane;copper Chemical compound [AlH3].[Cu] JRBRVDCKNXZZGH-UHFFFAOYSA-N 0.000 title abstract 3
- 239000000463 material Substances 0.000 claims abstract description 34
- 238000005728 strengthening Methods 0.000 claims abstract description 25
- 239000002245 particle Substances 0.000 claims abstract description 16
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000012360 testing method Methods 0.000 claims abstract description 14
- 230000005540 biological transmission Effects 0.000 claims abstract description 4
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 4
- 238000005457 optimization Methods 0.000 claims abstract description 4
- 230000008859 change Effects 0.000 claims description 23
- 238000001556 precipitation Methods 0.000 claims description 20
- 229910000881 Cu alloy Inorganic materials 0.000 claims description 18
- WPPDFTBPZNZZRP-UHFFFAOYSA-N aluminum copper Chemical compound [Al].[Cu] WPPDFTBPZNZZRP-UHFFFAOYSA-N 0.000 claims description 18
- 230000001419 dependent effect Effects 0.000 claims description 15
- 230000007246 mechanism Effects 0.000 claims description 15
- 238000004881 precipitation hardening Methods 0.000 claims description 9
- 238000010008 shearing Methods 0.000 claims description 9
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 7
- 229910052782 aluminium Inorganic materials 0.000 claims description 7
- 230000033001 locomotion Effects 0.000 claims description 7
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical group [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 claims description 6
- 230000002860 competitive effect Effects 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 239000006104 solid solution Substances 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 238000004090 dissolution Methods 0.000 claims description 4
- 238000011084 recovery Methods 0.000 claims description 4
- 230000003014 reinforcing effect Effects 0.000 claims description 4
- 235000015220 hamburgers Nutrition 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000002787 reinforcement Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000009864 tensile test Methods 0.000 claims description 2
- 238000001816 cooling Methods 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract description 4
- 238000002474 experimental method Methods 0.000 abstract description 3
- 239000000470 constituent Substances 0.000 abstract 1
- 230000004069 differentiation Effects 0.000 abstract 1
- 239000004576 sand Substances 0.000 abstract 1
- 239000004575 stone Substances 0.000 abstract 1
- 239000012071 phase Substances 0.000 description 34
- 238000003483 aging Methods 0.000 description 5
- 229910000838 Al alloy Inorganic materials 0.000 description 4
- 229910045601 alloy Inorganic materials 0.000 description 2
- 239000000956 alloy Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000006355 external stress Effects 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000002244 precipitate Substances 0.000 description 2
- 230000001376 precipitating effect Effects 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000000930 thermomechanical effect Effects 0.000 description 2
- 229910000906 Bronze Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000006911 nucleation Effects 0.000 description 1
- 238000010899 nucleation Methods 0.000 description 1
- 239000007790 solid phase Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- 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)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Hardware Design (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Geometry (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating And Analyzing Materials By Characteristic Methods (AREA)
Abstract
一种预测变应力下铝铜合金蠕变时效行为的方法,在蠕变试验机上进行单轴拉伸应力时效处理,在特定温度不同的应力水平下进行蠕变实验,获得不同应力下的蠕变曲线,然后将实验空冷到室温后进行拉伸实验,获得不同应力时效时屈服强度的变化;用X射线方法计算步骤一蠕变过程位错密度的演变,通过透射电镜观测得到蠕变过程主要强化相θ′和θ″尺寸和体积分数的变化;建立预测铝铜合金变应力下蠕变和时效强化的宏微观统一本构模型;确定材料常数,采用粒子群优化算法得到模型中的参数值。本发明解决了现有方法只能在较窄应力范围下使用的弊端,对准确预测实用材料和构件复杂应力下的蠕变量和性能演化具有重要的指导意义,可以应用在构件制造的数值模拟中。
Description
技术领域
本发明属于铝合金热加工和成形技术领域,特别涉及一种预测变应力下铝铜合金蠕变时效行为的方法。
背景技术
蠕变时效成形(CAF)通常用于制造航空航天工业中大型铝合金整体部件。通过蠕变/应力松弛和时效硬化,该工艺可以同时实现板材的成形和强化。弹性变形通过蠕变/应力松弛转化为塑性变形,而合金强度主要通过固相析出和少量蠕变时效成形过程位错增殖来提高。而带加强筋的板材与模具表面接触时将受到复杂非均匀应力。例如,加筋部位的应力要高的多,甚至在初始加载阶段超过了铝合金的屈服强度。因此,了解不同应力水平下的变形和强化行为对CAF的准确预测至关重要,特别是对于具有复杂结构的构件。
铝铜合金广泛应用于制造轻质航空航天部件,在铝铜合金中不同晶体结构和大小的相共存导致不同强化机制,在实际的CAF过程中,需要建立能准确量化各机制综合效应的模型来描述析出强化。在蠕变时效过程中施加的应力产生了热-机械耦合情况,增加了析出的复杂性。同时,在高应力位置发生的塑性变形会引起一定程度的位错,对析出相的成核和生长产生很大影响。因此,在模拟蠕变时效过程中屈服强度的演化时,需要调整析出强化方程。析出相倾向于以异质形核的方式在预先存在的位错上形成,析出路径放生变化。但在蠕变时效过程中,描述这种行为的模型很少。
目前尽管存在一些描述蠕变时效行为的唯象模型,但不同应力下蠕变时效的物理机制和预测仍然很难解决。一个能预测不同应力下蠕变变形和析出强化共存的本构模型对CAF过程的模拟和设计是至关重要的。目前的模型主要考虑了位错硬化和回复,以及溶质析出和溶解。它们在描述主导机制发生变化的蠕变时效行为方面存在很大局限性,所以不足以预测大应力范围下的蠕变时效行为。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种预测变应力下铝铜合金蠕变时效行为的方法,在系统研究应力对铝铜合金蠕变时效行为影响的基础上,分析了蠕变时效应力水平相关性的机理,构建的本构模型来预测所得到的强度和蠕变应变随外加应力的变化,模型预测数据与实验结果吻合。本发明解决了现有方法只能在较窄应力范围下使用的弊端,对准确预测实用材料和构件复杂应力下的蠕变量和性能演化具有重要的指导意义,可以应用在构件制造的数值模拟中。
为了实现上述目的,本发明采用的技术方案是:
一种预测变应力下铝铜合金蠕变时效行为的方法,包括:
步骤一:在蠕变试验机上进行单轴拉伸应力时效处理(蠕变时效),在特定温度(120-200℃)不同的应力水平(应力值从低到高于屈服强度)下进行蠕变实验,获得不同应力下的蠕变曲线,然后将实验空冷到室温后进行拉伸实验,获得不同应力时效时屈服强度的变化;
步骤二:用X射线方法计算步骤一蠕变过程位错密度的演变,通过透射电镜观测得到蠕变过程主要强化相θ′和θ″尺寸和体积分数的变化;
步骤三:建立预测铝铜合金变应力下蠕变和时效强化的宏微观统一本构模型,建立的模型可以描述在较大应力范围下的复杂蠕变时效行为,蠕变时效过程中产生的内部微结构变量以速率形式给出,将材料属性与这些变量关联,模型特别考虑与蠕变时效过程有关的位错和析出的相互作用,以提高预测能力。使用无量纲时间和应力来对方程进行无量纲化。模型具体包含描述屈服强度σy演变的强度方程和描述蠕变应变与应力关系即蠕变应变率的变形方程,过程如下:
1),计算加载阶段应力导致的位错,当施加的应力超过材料的屈服强度时,在加载阶段可以产生明显的位错,得到加载阶段位错密度的方程:
ρI0=c1·(σ-σIY)2+ρq
其中ρI0是加载阶段位错密度,c1是材料相关常数,σ是施加应力,σIY是加载前材料的初始屈服强度,ρq是假定的淬火位错密度,为0.05×1014m-2;
2),建立应力时效过程位错演化的方程,应力时效过程中的位错演变分为两部分:在加载阶段积累位错的静态回复和与蠕变有关的位错演变;
其中是标准位错密度率,c2和c3是材料相关常数,n1是常数,是标准位错密度;
3),拟合通过XRD测试测量的位错密度数据,获得标准化位错密度曲线;
4),高应力时效时θ′/θ″竞争析出,建立如下析出相体积分数演化的方程描述其竞争关系:
其中,ftotal是析出相体积分数,是ftotal随时间变化的速率,是fθ′随时间变化的速率,是fθ″随时间变化的速率,fθ′是θ′的相对体积分数,fθ″是θ″的相对体积分数,c4、c5、c6、n2、n3是材料相关常数,t是时间,是位错密度,c7是铝基体中铜原子转变为θ″相的百分比,c8是θ″到θ′时铜原子的溶解析出率;
5),建立描述析出相尺寸演化的方程,析出相尺寸是决定析出硬化潜力的重要参数,在等温应力时效过程中,析出相主要通过体积扩散机制生长,其中析出相大小与时效时间的平方根成比例,满足菲克第二扩散定律,应力诱发的位错可以显着影响析出相的生长,因此,方程中涉及应力项和位错项,它们以直径D和T厚度的生长速率的形式写出,如下:
其中,是Dθ″随时间变化的速率,是Dθ′随时间变化的速率,是Tθ″随时间变化的速率,是Tθ′随时间变化的速率,也即成长率,Dθ″是θ″的平均直径,Dθ′是θ′的平均直径,Tθ″是θ″的平均厚度,Tθ′是θ′的平均厚度,c9、c10、c11、c12、c13、c14、c15、c16、c17、c18、n4、n5、n6、n7是材料相关常数;
6),建立可剪切θ″相的强化方程:
其中Γ是位错线张力,Lp是平均中心到中心粒子间距,F是颗粒对位错剪切的抵抗力的量度,M是泰勒因子,b是伯格斯矢量;
考虑析出相的应力位向效应Lpm=Lp/β,其中Lpm是改进的平均中心到中心粒子间距,β是应力位向常数,因此考虑到θ″粒子的优先取向,强化方程变为:
Γ=(Gb2/2π)·ln(ro/ri)
其中γi是θ″和铝基体间的界面能量,r0和ri分别是外部和内部的截止距离,是特定相θm厚度的平方;
7),建立不可剪切θ'相的强化方程:
其中ν是泊松比,G是剪切模量,fma是最大体积分数,用于归一化;
8),建立析出硬化的方程,考虑剪切和非剪切两种析出强化机制,整体析出硬化方程为:
(σppt)1.4=(σshearable)1.4+(σshear-resistant)1.4
由此得到描述屈服强度σy演变的强度方程为:
即,屈服强度σy由混合规则给出,强化组成:内在强度σi、位错强度σp、析出强化σppt和固溶强化σss;
9),描述蠕变应变与应力关系的变形方程基于经典的奥罗万方程式,是位错密度和平均位错速度的乘积,后者受到驱动位错运动的外部应力和阻碍位错运动的位错或析出相等障碍物的影响:
即,蠕变时效过程中产生的内部微结构变量以速率形式给出,且材料属性与这些变量关联;
其中,v代表平均位错速度,c20、c21和n8是材料相关常数,k1和k2是加权因子代表不同强化机制对位错运动的抵抗力,σρ是位错贡献的强度;
步骤四:确定材料常数,采用粒子群优化算法得到所述模型中的参数值。大量参数使得必须逐步解耦模型并拟合参数,过程如下:
首先,将蠕变应变曲线与多项式分别拟合,以获得函数形式的蠕变应变率而不是散射点;
其次,用多项式拟合位错演化方程,相对体积分数,析出相大小演变和固溶体硬化贡献;
最后,拟合屈服强度以及内在强度和应力定向系数。
与现有技术相比,本发明的有益效果是:
材料及其构件在蠕变时经历应力时效处理,实际过程热机械条件较复杂,各种应力下的蠕变同时发生。然而,现有的方法无法描述高应力水平下蠕变应变的剧烈变化以及应力对时效硬化变化的影响。这源于对蠕变时效行为随应力增加转变机制的理解不准确有关。本方法考虑高应力引起的位错及其对析出行为的影响,开发的统一本构模型可以描述在大应力范围下的复杂蠕变时效行为。本预测方法对于蠕变和蠕变时效成形过程的模拟至关重要。
附图说明
图1表示预测和实验铝铜合金屈服强度曲线。
图2表示预测和实验铝铜合金蠕变曲线。
具体实施方式
以下结合附图和具体实施方式对本发明方法进一步阐释并进行数据说明。实施例所用样品为商用2219铝合金。拉伸试验采用标准为:GB/T 228-2002。高温蠕变时效试验在在三思泰捷公司生产的RMT-D10型电子式高温蠕变持久强度试验机上进行,试验机的控温精度为±2℃,载荷精度为±3N。
本发明提供一种精确预测铝铜合金变应力下蠕变时效行为的方法,下面以2219合金(典型铝铜合金)的蠕变和时效强化行为的预测为例,详细介绍本发明提供的预测方法和具体实施细节,涉及的方法如下:
步骤一:在蠕变试验机上进行了单轴拉伸应力时效处理(蠕变时效),在165℃不同的应力水平下(60MPa,90MPa,120MPa,150MPa,180MPa,210MPa)进行1-48h的蠕变实验,获得不同应力下的蠕变曲线。然后将实验空冷到室温后进行拉伸实验,获得不同应力时效时屈服强度的变化。
步骤二:用X射线方法计算上述蠕变过程位错密度的演变。通过透射电镜观测得到上述蠕变过程主要强化相θ′和θ″尺寸和体积分数的变化。
步骤三:建立预测铝铜合金变应力下蠕变和时效强化的宏微观统一本构模型。建立的模型可以描述在较大应力范围下的复杂蠕变时效行为。蠕变时效过程中产生的内部微结构变量以速率形式给出,将材料属性与这些变量关联,开发的模型包含描述屈服强度演变的强度方程和描述蠕变应变与应力关系的变形方程。屈服强度σy由混合规则给出,强化组成:内在强度σi、位错强度σp、析出强化σppt和固溶强化σss;。
模型特别考虑与蠕变时效过程有关的位错和析出的相互作用,以提高预测能力。使用无量纲时间和应力来对方程进行无量纲化。
计算加载阶段应力导致的位错,当施加的应力超过材料的屈服强度时,在加载阶段可以产生明显的位错:
ρIO=c1·(σ-σIY)2+ρq
其中ρI0是加载阶段位错密度,c1是材料相关常数,σ是施加应力,σIY是加载前材料的初始屈服强度,ρq是假定的淬火位错密度,为0.05×1014m-2。
建立应力时效过程位错演化的方程,应力时效过程中的位错演变分为两部分:在加载阶段积累位错的静态回复和与蠕变有关的位错演变。
其中是标准位错密度率,c2和c3是材料相关常数,n1是常数,是标准位错密度;
拟合通过XRD测试测量的位错密度数据,获得标准化位错密度曲线。
建立析出硬化的方程,考虑剪切和非剪切两种析出强化机制,整体析出硬化方程为:
(σppt)1.4=(σshearable)1.4+(σshear-resistant)1.4
建立描述析出相体积分数演化的方程,高应力时效时θ′/θ″竞争析出,建立如下模型描述其竞争关系:
其中,ftotal是析出相体积分数,是ftotal随时间变化的速率,是fθ′随时间变化的速率,是fθ″随时间变化的速率,fθ′是θ′的相对体积分数,fθ″是θ″的相对体积分数,c4、c5、c6、n2、n3是材料相关常数,t是时间,是位错密度,c7是铝基体中铜原子转变为θ″相的百分比,c8是θ″到θ′时铜原子的溶解析出率。
建立描述析出相尺寸演化的方程,析出相尺寸是决定析出硬化潜力的重要参数。在等温应力时效过程中,析出相主要通过体积扩散机制生长,其中析出相大小与时效时间的平方根成比例,满足菲克第二扩散定律。应力诱发的位错可以显着影响析出相的生长,因此,方程中涉及应力项和位错项,它们以D
(直径)和T(厚度)的生长速率的形式写出,如下:
其中,是Dθ″随时间变化的速率,是Dθ′随时间变化的速率,是Tθ″随时间变化的速率,是Tθ′随时间变化的速率,也即成长率,Dθ″是θ″的平均直径,Dθ′是θ′的平均直径,Tθ″是θ″的平均厚度,Tθ′是θ′的平均厚度,c9、c10、c11、c12、c13、c14、c15、c16、c17、c18、n4、n5、n6、n7是材料相关常数;方程以生长速率的形式给出,其直接受位错密度的影响。
建立可剪切θ“相的强化方程:
其中Γ是位错线张力,Lp是平均中心到中心粒子间距,F是颗粒对位错剪切的抵抗力的量度,M是泰勒因子,b是伯格斯矢量。考虑析出相的应力位向效应Lpm=Lp/β,其中Lpm是改进的平均中心到中心粒子间距,β是应力位向常数,
因此考虑到θ″粒子的优先取向,强化方程变为:
Γ=(Gb2/2π)·ln(ro/ri)
其中γi是θ″和铝基体间的界面能量,r0和ri分别是外部和内部的截止距离,是特定相θm厚度的平方;
建立不可剪切θ'相的强化方程:
其中ν是泊松比,G是剪切模量,fma是最大体积分数,用于归一化;
建立蠕变应变-应力关系的方程,基于经典的奥罗万方程式,蠕变应变率是位错密度和平均位错速度的乘积。后者受到驱动位错运动的外部应力和阻碍位错运动的位错或析出相等障碍物的影响:
其中,v代表平均位错速度,c20、c21和n8是材料相关常数,k1和k2是加权因子代表不同强化机制对位错运动的抵抗力,σρ是位错贡献的强度。
步骤四:确定材料常数采用粒子群优化算法得到所提模型中的参数值。大量参数使得必须逐步解耦模型并拟合参数。首先,将蠕变应变曲线与多项式分别拟合,以获得函数形式的蠕变应变率而不是散射点。其次,多项式用于拟合位错演化方程,然后是相对体积分数,析出相大小演变和固溶体硬化贡献。最后,拟合屈服强度,以及内在强度和应力定向系数。最后,蠕变应变曲线统一到蠕变和应力关系的方程,并获得所有参量。本构模型中的参数列于下表1中。
表1:拟合得到的材料常数
本发明提供的方法使宏观变形和微观结构演化耦合,可用于预测在纯弹性区域内的应力范围内的蠕变时效行为以及产生初始塑性变形的应变范围内的蠕变时效行为。通过考虑不同应力下的位错密度和平均位错速度,可以反映蠕变行为的应力相关性。将屈服强度和蠕变的预测结果与实验数据进行比较,如图1和2所示,可看出预测值和实验值匹配非常好。
本发明方法可以准确地描述铝铜合金在各种应力水平下的蠕变时效行为,特别是高应力下蠕变应变的显著增加。计算结果与低应力和高应力时效铝铜合金的总屈服强度演变的实验数据也非常吻合,实现了形/性演化行为的高精度预测。
Claims (4)
1.一种预测变应力下铝铜合金蠕变时效行为的方法,其特征在于,包括:
步骤一:在蠕变试验机上进行单轴拉伸应力时效处理,在特定温度不同的应力水平下进行蠕变实验,获得不同应力下的蠕变曲线,然后将实验空冷到室温后进行拉伸实验,获得不同应力时效时屈服强度的变化;
步骤二:用X射线方法计算步骤一蠕变过程位错密度的演变,通过透射电镜观测得到蠕变过程主要强化相θ′和θ″尺寸和体积分数的变化;
步骤三:建立预测铝铜合金变应力下蠕变和时效强化的宏微观统一本构模型,包含描述屈服强度σy演变的强度方程和描述蠕变应变与应力关系即蠕变应变率的变形方程,过程如下:
1),计算加载阶段应力导致的位错,当施加的应力超过材料的屈服强度时,在加载阶段可以产生明显的位错,得到加载阶段位错密度的方程:
ρ1O=c1·(σ-σIY)2+ρq
其中ρIO是加载阶段位错密度,c1是材料相关常数,σ是施加应力,σIY是加载前材料的初始屈服强度,ρq是假定的淬火位错密度,为0.05×1014m-2;
2),建立应力时效过程位错演化的方程,应力时效过程中的位错演变分为两部分:在加载阶段积累位错的静态回复和与蠕变有关的位错演变;
其中是标准位错密度率,c2和c3是材料相关常数,n1是常数,是标准位错密度;
3),拟合通过XRD测试测量的位错密度数据,获得标准化位错密度曲线;
4),高应力时效时θ′/θ″竞争析出,建立如下析出相体积分数演化的方程描述其竞争关系:
其中,ftota1是析出相体积分数,是ftotal随时间变化的速率,是f′θ随时间变化的速率,是f″θ随时间变化的速率,f′θ是θ′的相对体积分数,f″θ是θ″的相对体积分数,c4、c5、c6、n2、n3是材料相关常数,t是时间,是位错密度,c7是铝基体中铜原子转变为θ″相的百分比,c8是θ″到θ′时铜原子的溶解析出率;
5),建立描述析出相尺寸演化的方程,析出相尺寸是决定析出硬化潜力的重要参数,在等温应力时效过程中,析出相主要通过体积扩散机制生长,其中析出相大小与时效时间的平方根成比例,满足菲克第二扩散定律,应力诱发的位错可以显着影响析出相的生长,因此,方程中涉及应力项和位错项,它们以直径D和T厚度的生长速率的形式写出,如下:
其中,是Dθ″随时间变化的速率,是Dθ′随时间变化的速率,是Tθ″随时间变化的速率,是Tθ′随时间变化的速率,也即成长率,Dθ″是θ″的平均直径,Dθ′是θ′的平均直径,Tθ″是θ″的平均厚度,Tθ′是θ′的平均厚度,c9、c10、c11、c12、c13、c14、c15、c16、c17、c18、n4、n5、n6、n7是材料相关常数;
6),建立可剪切θ″相的强化方程:
其中Γ是位错线张力,Lp是平均中心到中心粒子间距,F是颗粒对位错剪切的抵抗力的量度,M是泰勒因子,b是伯格斯矢量;
7),建立不可剪切θ′相的强化方程:
其中v是泊松比,G是剪切模量,fma是最大体积分数,用于归一化;
8),建立析出硬化的方程,考虑剪切和非剪切两种析出强化机制,整体析出硬化方程为:
(σppt)1.4=(σshearable)1.4+(σshear-resistant)1.4
由此得到描述屈服强度σy演变的强度方程为:
即,屈服强度σy由混合规则给出,强化组成:内在强度σi、位错强度σp、析出强化σppt和固溶强化σss;
9),描述蠕变应变与应力关系的变形方程基于经典的奥罗万方程式,是位错密度和平均位错速度的乘积,为:
即,蠕变时效过程中产生的内部微结构变量以速率形式给出,且材料属性与这些变量关联;
其中,v代表平均位错速度,c20、c21和n8是材料相关常数,k1和k2是加权因子代表不同强化机制对位错运动的抵抗力,σρ是位错贡献的强度;
步骤四:确定材料常数,采用粒子群优化算法得到所述模型中的参数值。
2.根据权利要求1所述预测变应力下铝铜合金蠕变时效行为的方法,其特征在于,所述步骤一中,所述特定温度的取值范围为120-200℃。
3.根据权利要求1所述预测变应力下铝铜合金蠕变时效行为的方法,其特征在于,所述步骤三中,考虑析出相的应力位向效应Lpm=Lp/β,其中Lpm是改进的平均中心到中心粒子间距,β是应力位向常数,因此考虑到θ″粒子的优先取向,强化方程变为:
Γ=(Gb2/2π)·ln(ro/ri)
其中γi是θ″和铝基体间的界面能量,r0和ri分别是外部和内部的截止距离,是特定相θm厚度的平方。
4.根据权利要求1所述预测变应力下铝铜合金蠕变时效行为的方法,其特征在于,所述步骤四中,逐步解耦模型并拟合参数,过程如下:
首先,将蠕变应变曲线与多项式分别拟合,以获得函数形式的蠕变应变率而不是散射点;
其次,用多项式拟合位错演化方程,相对体积分数,析出相大小演变和固溶体硬化贡献;
最后,拟合屈服强度以及内在强度和应力定向系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811123118.8A CN109214104B (zh) | 2018-09-26 | 2018-09-26 | 一种预测变应力下铝铜合金蠕变时效行为的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811123118.8A CN109214104B (zh) | 2018-09-26 | 2018-09-26 | 一种预测变应力下铝铜合金蠕变时效行为的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109214104A true CN109214104A (zh) | 2019-01-15 |
CN109214104B CN109214104B (zh) | 2022-12-20 |
Family
ID=64981671
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811123118.8A Active CN109214104B (zh) | 2018-09-26 | 2018-09-26 | 一种预测变应力下铝铜合金蠕变时效行为的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109214104B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110232223A (zh) * | 2019-05-27 | 2019-09-13 | 西北工业大学 | 镍基单晶高温合金微结构筏化类型的预测方法 |
CN111145842A (zh) * | 2020-02-27 | 2020-05-12 | 西安交通大学 | 一种面向高速变形过程的多相金属材料本构行为构建方法 |
CN111914438A (zh) * | 2020-08-27 | 2020-11-10 | 浙江华电器材检测研究所有限公司 | 老化态纯铝线或铝合金线屈服强度的预测公式和预测方法 |
CN112485113A (zh) * | 2020-11-17 | 2021-03-12 | 核工业西南物理研究院 | 一种小尺寸样品的材料拉伸性能测试方法及装置 |
CN112733330A (zh) * | 2020-12-22 | 2021-04-30 | 武汉理工大学 | 铝合金热变形行为和成形后屈服强度的预测方法 |
CN113076627A (zh) * | 2021-03-17 | 2021-07-06 | 湖南大学 | 一种镍基高温合金宽温域蠕变强度预测方法 |
CN113976707A (zh) * | 2021-11-03 | 2022-01-28 | 西北工业大学 | 一种铝合金大曲率复杂壁板构件复合成形方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103422035A (zh) * | 2012-09-20 | 2013-12-04 | 中南大学 | 一种用于Al-Cu-Mg系合金板材的蠕变时效成形方法 |
CN105004617A (zh) * | 2015-07-20 | 2015-10-28 | 沈阳工业大学 | 一种描述金属材料蠕变曲线的方法 |
CN105784508A (zh) * | 2016-04-11 | 2016-07-20 | 沈阳工业大学 | 一种表征单晶Ni基合金蠕变性能的方法 |
-
2018
- 2018-09-26 CN CN201811123118.8A patent/CN109214104B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103422035A (zh) * | 2012-09-20 | 2013-12-04 | 中南大学 | 一种用于Al-Cu-Mg系合金板材的蠕变时效成形方法 |
CN105004617A (zh) * | 2015-07-20 | 2015-10-28 | 沈阳工业大学 | 一种描述金属材料蠕变曲线的方法 |
CN105784508A (zh) * | 2016-04-11 | 2016-07-20 | 沈阳工业大学 | 一种表征单晶Ni基合金蠕变性能的方法 |
Non-Patent Citations (2)
Title |
---|
"汽车工业用力学性能和抗热蠕变性能均好的Al-Cu合金", 《铝加工》 * |
李喜财等: "基于蠕变时效交互作用机理的2219铝合金统一本构建模", 《中南大学学报(自然科学版)》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110232223A (zh) * | 2019-05-27 | 2019-09-13 | 西北工业大学 | 镍基单晶高温合金微结构筏化类型的预测方法 |
CN110232223B (zh) * | 2019-05-27 | 2022-09-27 | 西北工业大学 | 镍基单晶高温合金微结构筏化类型的预测方法 |
CN111145842A (zh) * | 2020-02-27 | 2020-05-12 | 西安交通大学 | 一种面向高速变形过程的多相金属材料本构行为构建方法 |
CN111145842B (zh) * | 2020-02-27 | 2022-07-12 | 西安交通大学 | 一种面向高速变形过程的多相金属材料本构行为构建方法 |
CN111914438A (zh) * | 2020-08-27 | 2020-11-10 | 浙江华电器材检测研究所有限公司 | 老化态纯铝线或铝合金线屈服强度的预测公式和预测方法 |
CN112485113A (zh) * | 2020-11-17 | 2021-03-12 | 核工业西南物理研究院 | 一种小尺寸样品的材料拉伸性能测试方法及装置 |
CN112485113B (zh) * | 2020-11-17 | 2023-04-21 | 核工业西南物理研究院 | 一种小尺寸样品的材料拉伸性能测试方法及装置 |
CN112733330A (zh) * | 2020-12-22 | 2021-04-30 | 武汉理工大学 | 铝合金热变形行为和成形后屈服强度的预测方法 |
CN112733330B (zh) * | 2020-12-22 | 2022-09-23 | 武汉理工大学 | 铝合金热变形行为和成形后屈服强度的预测方法 |
CN113076627A (zh) * | 2021-03-17 | 2021-07-06 | 湖南大学 | 一种镍基高温合金宽温域蠕变强度预测方法 |
CN113976707A (zh) * | 2021-11-03 | 2022-01-28 | 西北工业大学 | 一种铝合金大曲率复杂壁板构件复合成形方法 |
CN113976707B (zh) * | 2021-11-03 | 2023-09-01 | 西北工业大学 | 一种铝合金大曲率复杂壁板构件复合成形方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109214104B (zh) | 2022-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109214104B (zh) | 一种预测变应力下铝铜合金蠕变时效行为的方法 | |
Li et al. | A unified constitutive model for asymmetric tension and compression creep-ageing behaviour of naturally aged Al-Cu-Li alloy | |
Chen et al. | Stress relaxation behavior of an Al–Zn–Mg–Cu alloy in simulated age-forming process | |
Zhang et al. | Constitutive modeling for creep age forming of heat-treatable strengthening aluminum alloys containing plate or rod shaped precipitates | |
Xu et al. | Experimental research on creep aging behavior of Al-Cu-Mg alloy with tensile and compressive stresses | |
Wang et al. | Anomalous behaviors of a single-crystal Nickel-base superalloy over a wide range of temperatures and strain rates | |
CN106503397B (zh) | 一种获得金属材料再结晶组织演化晶界可动性参数的方法 | |
Zheng et al. | Experimental investigation of multi-step stress-relaxation-ageing of 7050 aluminium alloy for different pre-strained conditions | |
Zhan et al. | Effects of process parameters on mechanical properties and microstructures of creep aged 2124 aluminum alloy | |
CN108256179B (zh) | 一种预测材料蠕变曲线的方法 | |
Whittaker et al. | An analysis of modern creep lifing methodologies in the titanium alloy Ti6-4 | |
Wang et al. | Prediction on hot deformation behavior of spray formed ultra-high strength aluminum alloy—A comparative study using constitutive models | |
Yang et al. | Tension-compression asymmetry of stress-relaxation ageing behavior of AA2219 alloy over a wide range of stress levels | |
Herasymchuk | Microstructurally-dependent model for predicting the kinetics of physically small and long fatigue crack growth | |
Zhang et al. | Dynamic impact constitutive relation of 6008-T6 aluminum alloy based on dislocation density and second-phase particle strengthening effects | |
Chandra et al. | Deformation behavior of Nickel-based superalloy Su-263: Experimental characterization and crystal plasticity finite element modeling | |
Selyutina et al. | Structural–temporal peculiarities of dynamic deformation of nanostructured and nanoscaled metals | |
Fan et al. | Fatigue life evaluation for notched single-crystal Ni-based superalloys considering inhomogeneous rafting microstructure | |
Xia et al. | Effect of different pre-treatments on anisotropy in stress relaxation aging behavior of AA2195 alloy | |
CN113008677B (zh) | 一种镍基高温合金的蠕变持久性能预测方法 | |
Wei et al. | Model for texture evolution in cold rolling of 2.4 wt.-% Si non-oriented electrical steel | |
Dancette et al. | Multisite model prediction of texture induced anisotropy in brass | |
CN112733330B (zh) | 铝合金热变形行为和成形后屈服强度的预测方法 | |
Ma et al. | Parametric study of cyclic plasticity behavior in a directionally solidified superalloy with partial recrystallization by crystal plasticity finite element simulation | |
Friis et al. | Work hardening behaviour of heat-treatable Al-Mg-Si-alloys |
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 |