6xxx本构关系讲解.docx
《6xxx本构关系讲解.docx》由会员分享,可在线阅读,更多相关《6xxx本构关系讲解.docx(14页珍藏版)》请在冰豆网上搜索。
![6xxx本构关系讲解.docx](https://file1.bdocx.com/fileroot1/2022-10/18/ef6115a7-71ad-4728-83b3-683aeeb1e3f1/ef6115a7-71ad-4728-83b3-683aeeb1e3f11.gif)
6xxx本构关系讲解
6061、6063、6082铝合金高温变形行为研究及本构方程
一.试验过程
1)试验前压缩试样加丄成两端带有凹槽(09mmXO・2mm)的试样(GlOmmX
2)在试样上焊上两根用来测量温度的金属丝,这两根金属丝是不同的,其
中一根有磁性,而另一根则无,金属丝在不相互接触的条件下应尽可能的接近。
3)凹槽内填充润滑剂(石墨乳),变形时,封闭在腔体内的润滑剂可以减小
平面压头与试样接触面的摩擦,从而减少样品的不均匀变形。
开动气动阀使压头夹紧试样,要注意对中;并将两根金属丝接在相应地接头上,需要注意,有磁性的金属丝和无磁性的金属纟幺所接位置不同;测量应变的玻璃仪器贴着试样安放好,如图1・2所示。
Load
Load
4)在计算机上设置控制参数,并调节与检査好仪器,准确无误后即可启动
计算机的程序开始模拟压缩实验。
5)所有试样均利用自身电阻进行加热,采用Ni-NiAl热电偶直接焊在试样
中部连续测温,升温速率10(rC/min,达到所设定的温度后,保温3nun后进行恒温恒应变速率的压缩试验。
6)变形后立即对试样进行水淬,以冻结变形组织,用于金相组织分析,水
淬延迟时间大约为0・5s。
7)取出压缩后的试样,山Gleeble・1500系统的计算机自动采集真应力、真
应变、压力、温度、时间等数据。
为获得较为恒定的应变速率,压头位移速度匕按下式进行设定。
式中g为应变速率,力为样品瞬时高度,&为真应变,每隔0.1真应变值分段控
二.本构方程的建立2J材料模型
对于不同材料髙温塑性变形的研究发现,材料变形时的应力水平和应变速率.温度之间满足指数关系:
式中b为一定温度和应变条件下的流变应力,名为真实应变,cTjn和m(r)为与温度有关的常数,小应变条件下,这些常数随应变发生变化,一旦进入稳态流变阶段,则一定温度下它们保持恒定。
热变形过程中,材料在任何应变或稳态下的高温流变应力<7强烈地取决于变形温度r和应变速率对不同热加工数据的仔细研究表明,低应力水平下稳态流变应力<7和应变速率g之间的关系可用指数关系进行描述:
(1-2)
言=A*"
式中料|为与温度无关的常数。
而在高应力水平下稳态流变应力<7和应变速率g之间的关系可用幕指数关系来加以描述:
(1-3)
式中0也是与温度无关的常数。
这些关系描述了应变硬化和动态软化过程之间的动态平衡,与稳态蠕变变形对应的关系非常相似。
根据这种相似性,
Sellars和Tegart于1966年提出了一种包含变形激活能Q和温度T的双曲正弦形式的修正Arrhenius关系来描述这种热激活稳态变形行为:
£=Asinh(ab)]"exp
‘2
、而丿
(1-4)
式(1・4)中A、a.H为于温度无关的常数,/?
为气体常数,r为绝对温度。
比
较式(1・2)、式(1・3)和式(1・4)可以发现,在低应力水平下(abvO・8),式(14)接近式(1・2)的指数关系,高应力水平(ab>l・2)时则接近(1・3)式的幕指数关系,常数a、尸和畀之间满足Q=创5因此,0和《可山低应力水平下的
实验数据求解。
热加工变形时的应变速率通常比蠕变时的应变速率大儿个数量级,但山于蠕变和热加工均属于热激活过程,热加工可视作蠕变在大应变速率和较高应力水平条件下的一种外延,两者的变形机制和软化机制都非常相似,因此它们都可以用热激活的Arrhenius式(1・4)进行描述。
本实验中铝合金热压缩变形就属于这种1W况,故可用该式来进行描述。
Zener和HoHomon在1944年提出并实验了一种确定钢高速拉伸实验应力一
应变关系的方法。
在室温和低于室温变形时,钢的应力一应变关系取决于应变速
率g和温度几农和T的关系可用一项参数Z表示,BIJ:
b=b(Z,f)
(1-5)
该参数Z包含激活能2项:
Z=gexp
(2)。
变形激活能2通常和激活熔AH相
RT
等,它提供了速率控制机制中原子重排难易程度的有关信息,山于高温塑性变形
存在热激活过程,也是Zener和Hollomon提及的条件,据此可将式(1・5)写成:
Z=Z(6g)
(1-6)
它依赖于流变应力<7而与温度无关。
Z打<7之间遵从下述关系:
Z=Asinh(ab)
(1-7)
研究表明,该式在较宽应变速率和温度范围内与实验数据吻合得较好。
这样,
我们可以得到所谓“温度补偿应变速率”,即Zener-Hollomon参数Z值的
定义:
Z=£exp
"刈"=A[sinh(\RT
(1-8)
式中4,Zba和A//均属于材料常数。
一般来说,A在高应力水平时为与速率控制机制中热激活位置成正比的结构因子。
当应力降低时,A,“和a的物理意义也发生变化,常数a为温度补偿应变速率和流变应力之间的相关性从指数关系变化到探函数时对应的流变应力的倒数,畀为温度补偿应变速率敏感性的倒数,A为与变形材料内部激活位置密度、空位浓度.位错上割阶的平均间距、位错柏氏矢量、原子配数、跃迁频率以及激活爛有关的函数。
为了更好的研究材料在变形时的力学行为,还应了解与应变速率和温度有关
(1-9)
的流变应力<7的变化规律。
从式(1・8)可以推出:
smh(ab)=1—
\A>
根据双曲正弦函数的定义,应有:
由此可以将流变应力表述成应变速率和温度的函数,亦可表达成Zener-Hollo参
数Z值的函数:
1/2
<7=—In'
a
只要知道A,小a、等材料常数,便可求得材料在任意变形条件下的流变应力值。
为了增加式(MI)对各种变形条件的普遍适用性,还应考虑材料变形时流变应力的应变敏感性。
但对高温塑性变形过程来说,若材料的软化足以抵消硬化的作用,可以认为材料在稳态变形阶段流变应力是应变不敬感的,因而可以忽略应变大小对流变应力的影响。
2・2流变应力曲线
①6063铝合金
0.2
0.6
0.4Tn^Strain
oooooooo
64208642
11111111
0.6
0.20.4
Tn^Strain
0.0
0.0
sans«2H
TrueStrain
图36063铝合金不同变形条件下的真应力一真应变曲线(a〜d)
表16063铝合金压缩变形时实测峰值应力(MPa)
^度co
应变速率
300
350
400
450
5
79.27
63.99
58.20
47.63
10
80.27
67.59
59.63
54.79
20
99.44
82.30
80.94
78.90
30
101.18
86.74
76.86
72.24
6063铝合金流变应力模型
利用以上模型,建立6063铝合金的流变应力模型,对峰值应力其应变速率
为:
(M2)
g=AF(b)cxp(—Q/R7*)
式中,F(b)是应力的函数,可以表示为以下三种形式:
F(b)=b"
当acr<0.8
(M3)
F(b)=exp(0b)
当07>1・2
(M4)
e为变形激活能,R为气体常数,r为绝对温度,a、0、H和A为材料
常数。
同时,式(1-12)可以很方便地表示为温度补偿应变速率参数,Zener-Hollomon
参数乙
Z=£exp(0/RT)
(1-17)
在高应力和低应力下,式(1・1)可分别表示为:
£=BL
(M8)
s=B'exp(Qb)
(M9)
对式(M8)和武(1・19)两边分别求对数得:
In(^)=ln(B)+nln(b)
(1-20)
根据实验结果绘制的峰值应力与变形速率、变形温度之间的关系曲线,如图
4所示。
从图中可看出,稳态流变应力和应变速率的双对数关系、流变应力的双
曲线正弦对数项和温度的倒数之间皆较好的满足线性关系。
山此可以认为6063
铝合金压缩变形时应力一应变速率关系满足双曲正弦形式,流变应力与变形温度
满足Arrhenius关系,即可以用包含Arrhenius项的Z参数描述6063铝合金变形
时的流变行为。
«值和0值可以通过式(1・20)和式(1-21)分别利用图l-3(a)(b)求ln(^)-ln(0=O」42MpG,«=8・47,此时对应的a=0.0168O
对所有应力状态有:
£=A[sinh(ab)rc)^(—Q/RT)
(1-22)
O-=丄ln{(彳)m+[(彳+1]"2}对式(1■⑵求导得:
(1-23)
(1-24)
求lnsinh(ac7)-(r-*)的斜率,即:
d[lnsinh(aJ)]/*厂|)=2469。
将此值和R、
«值代入(1・24)式,得:
Q=\l\.^^KJhnol
将e值代入(1・12)武,两边求对数得:
(1-25)
ln(£)=}x\A—QfRT+z/ln[sinh(cir<7)]
作图l・3(d)川Msinh(如)]Tn(g),可知两者截距即2//er-ln(A)的值,将
T值代入即可得到A值:
A=L904X10*V
6063铝合金峰值应力下的各参数值如下:
6
8642Q£
444443
2
0
A
1⑷
表26063铝合金峰值应力下的各参数值
n
ez/MPa'
0/MPaJ
Q/(KJ•ino|-')
A/s・*
8.47
0.0168
0.142
171.68
1.904X10'^
-ro
ri1,213141,51・6151・8
(C)
-6■
•2
-1
0
2345
Inf
(d)
2OS642024
11--
图46063铝合金峰值应力与温度、应变速率关系
(c)3对应5s'\10s-\30s-'
(a)(b)(d)I-4对应30(rC,35(rC,400€,450€
利用上述所得的数据,可以用包含Arrhenius项的Z参数描述6063铝合金变形时的流变行为,6063铝合金变形激活能e为171.68其Z参数可表述
(1-26)
为:
Z=^exp(17L68/Rr)s・i
流变应力、应变速率与温度的关系可用Z参数表示为:
b=59・5In{[Z/(l・904xl0j]""+{[(z/i・904xl0H)]2g7+l}"2}MPa(1-27)
从上面可以看出:
经典的双曲正弦本构方程可以精确地描述6063铝合金变形时的流变行为。
流变应力被描述成应变速率及变形温度的函数,在不同变形条件下,上述参数值都可认为在保持峰值应力时是不变的。
在进行低温快速挤压数值模拟时,只需把上式(1-27)输入到有限元中的材料本构关系模型中去,DEFORM软件利用内插和外插法,就能把材料在低温快速压缩试验中所得的应力应变曲线准确地描述出来。
@6061铝合金
160
140
120
100
80
60
40
20
0
0.00.20-40.6
TrueStrain
01g
160
140
120
100
80
60
40
20
0
0.00-20.6
TrueStem
6061铝合金不同变形条件下的真应力一真应变曲线(a〜d)
度CC)
应变速率
300
350
400
450
5
82.48
73.15
62.42
52.34
10
94.66
80.91
64.73
60.59
20
110.92
97