中国有色金属学报 2003,(06),1338-1342 DOI:10.19476/j.ysxb.1004.0609.2003.06.003
高能球磨Ti/Al粉末挤压固结致密过程数值分析
李小强 李元元 龙雁 胡连喜 王尔德
华南理工大学机械工程学院,华南理工大学机械工程学院,华南理工大学机械工程学院,哈尔滨工业大学材料科学与工程学院,哈尔滨工业大学材料科学与工程学院 广州510640 ,广州510640 ,广州510640 ,哈尔滨150001 ,哈尔滨150001
摘 要:
基于高能球磨过程中粉末的大塑性变形行为的分析 ,推导了球磨粉体的本构方程 ,并针对球磨Ti/Al粉末的挤压过程进行了有限元模拟。分析结果表明 ,随着球磨时间的延长 ,粉末固结致密难度增加 ,所需挤压力显著增大。该趋势和模拟数值均与挤压实验结果基本一致 ,验证了所推导的球磨粉体本构方程的合理性。
关键词:
Ti/Al粉末 ;高能球磨 ;本构关系 ;挤压 ;数值分析 ;
中图分类号: TF124
作者简介: 李小强(1972),男,博士;电话:02087113832;Email:lxqgyy@163.com;
收稿日期: 2003-01-08
基金: 国家自然科学基金资助项目 (5 9875 0 15 ); 华南理工大学自然科学基金资助项目 (E5 3 0 5 2 93 );
Numerical simulation of extrusion of high energy milled Ti/Al powders
Abstract:
On the basis of analysis of deformation behavior of powders during high energy milling, the plastic constitutive equation of milled powders was established. And using the equation, the extrusion of milled Ti/Al powders was further simulated by finite element method. The simulation results show that, with increasing the milling time, the densification of powders becomes more and more difficult, and the needed extrusion pressure substantially rises. Both the trend and the simulation data agree well with the experimental results, which affirms that the established plastic constitutive equation of milled powders is reasonable.
Keyword:
Ti/Al powder; high energy ball milling; plastic constitutive equation; extrusion; numerical simulation;
Received: 2003-01-08
自Benjamin等
[1 ]
提出高能球磨工艺以来, 由于它在细化粉末晶粒、 实现强制固溶以及制备非晶材料等方面的优越性, 因而受到了广泛关注并得到了飞速发展
[2 ,3 ,4 ,5 ,6 ]
。 然而,由于高能球磨工艺的固有特性, 它获得的只是粉末, 为了拓宽其应用领域, 往往需对球磨粉末进行有效固结成形处理。 由于球磨加工硬化和细晶强化等作用, 经球磨后粉体的塑性加工难度显著增加, 从而对固结致密工艺提出了更高的要求
[7 ,8 ]
。 为了优化致密参数和减少实验工作量, 有必要对球磨粉末固结致密过程进行数值模拟研究。 尽管目前对元素和预合金粉末致密过程的数值模拟已进行了较多的研究工作
[9 ,10 ,11 ]
, 同时高能球磨机理和球磨模型的研究也得到了进一步发展
[12 ,13 ,14 ]
, 然而, 关于球磨粉末本构方程的研究仍鲜见报道, 更不用说对在此基础上的球磨粉末固结致密成形过程进行数值模拟。 这在一定程度上阻碍了对球磨粉末致密行为的研究。 本文作者基于对高能球磨过程中粉末的大塑性变形行为的分析, 推导出了球磨粉体的本构方程, 并通过数值模拟方法研究了Ti/Al二元球磨粉末挤压致密成形过程, 对所推导的球磨粉末本构方程的合理性进行了验证。
1球磨粉体本构方程
1.1球磨粉末变形分析
根据Maurice
[12 ]
等的研究, 高能球磨过程中两个磨球每次碰撞捕获的粉末量V c 、 被捕获粉末变形前的初始高度h 0 、 粉末的变形量ε max 及碰撞接触时间2 τ 可分别描述为
V c = f p ? V s = ρ b / ( ρ p C R ) V t / V b ? 1 ? L f ? 1 3 π r 2 h ? ? ? ( 1 ) h 0 = 1 3 L f ? ρ b / ( ρ p C R ) V t / V b ? 1 ? ? ? ( 2 ) ε max = ln ( h 0 h 0 ? v τ / 2 ) ? ? ? ( 3 ) 2 τ = 5 . 5 7 4 4 v ? 0 . 2 ( ρ b / E b ) 0 . 4 R ? ? ? ( 4 )
V
c
=
f
p
?
V
s
=
ρ
b
/
(
ρ
p
C
R
)
V
t
/
V
b
-
1
?
L
f
?
1
3
π
r
h
2
?
?
?
(
1
)
h
0
=
1
3
L
f
?
ρ
b
/
(
ρ
p
C
R
)
V
t
/
V
b
-
1
?
?
?
(
2
)
ε
max
=
ln
(
h
0
h
0
-
v
τ
/
2
)
?
?
?
(
3
)
2
τ
=
5
.
5
7
4
4
v
-
0
.
2
(
ρ
b
/
E
b
)
0
.
4
R
?
?
?
(
4
)
式中 f p 为球磨桶中粉末相对体积含量, %;ρ p 为粉末密度; V t 和V b 分别为球磨桶容积和磨球的体积; C R 为球料质量比; h 为被捕获粉末变形后的高度; t 为变形时间; r h 为磨球的弹性变形区半径; ρ b 和E b 分别为磨球密度及其弹性模量; v 为磨球线速度; R 为球磨机搅拌杆长度; L f 为磨球空间运动平均自由程, 它是一个只与装球量、 磨球半径相关的参量。
当对摩尔比为1∶1(即体积比为16∶15)的Ti、 Al元素混合粉进行球磨时, 如果将复合粉中Ti、 Al看成是交替均匀层叠排列, 则Ti、 Al层厚度之比应为16∶15, 相邻两层间距可表示为
l = l 0 exp ( ? ∑ ε max ) = ? ? l 0 exp [ ? v m b 2 ( L f + 2 τ v ) V b ρ b ? V p V c ? t ? ? ? ln ( h 0 h 0 ? 1 . 3 9 3 7 v 0 . 8 ( ρ b / E b ) 0 . 4 R ) ] ? ? ? ( 5 )
l
=
l
0
exp
(
-
∑
ε
max
)
=
?
?
l
0
exp
[
-
v
m
b
2
(
L
f
+
2
τ
v
)
V
b
ρ
b
?
V
p
V
c
?
t
?
?
?
ln
(
h
0
h
0
-
1
.
3
9
3
7
v
0
.
8
(
ρ
b
/
E
b
)
0
.
4
R
)
]
?
?
?
(
5
)
式中 m b 为单个磨球的质量; l 0 和l 分别为起始和球磨过程中某一时刻粉末中组元层片的间距。 在球料比、 装料量、 装球量等磨球条件一定时, 式(5)可进一步写成
ε t = ? ln ( l / l 0 ) = A t B v ? 1 + C ? ln ( k k ? D v 0 . 8 )
ε
t
=
-
ln
(
l
/
l
0
)
=
A
t
B
v
-
1
+
C
?
ln
(
k
k
-
D
v
0
.
8
)
;
A =m b V p /V c ; B =2L f V b ρ b ; C =4τV b ρ b ;
D =1.393 7R (ρ b /E b )0.4 ; k =h 0 (6)
式中 A 、 B 、 C 、 D 和k 均为与球磨环境相关的正常数, ε t 为球磨时间为t 时粉末的总压缩应变量的绝对值。
1.2本构方程推导
从细观角度看, 可将球磨二元粉末看作是一种层片复合材料。 尽管单个粉末具有各向异性, 但由于粉坯是由众多颗粒无序取向构成的, 从宏观上仍可视为均匀的各向同性材料来处理。 因而, 可在粉坯上取包含某点的一个非无限小的微元体, 无论二元粉末是否经过机械球磨, 该微元体所包含的粉末颗粒数均远大于1, 并基于此处微元体中组元间的弱界面结合, 可假设单元内组元间的应力传递为无损传递。 当组元材质在一定温度下的流动应力方程为
σ i = C i ε ? n i i ε ? m i i + y i , 0
σ
i
=
C
i
ε
?
i
n
i
ε
?
i
m
i
+
y
i
,
0
时, 则根据复合材料力学理论, 全致密微元体的应力
σ ? c
σ
?
c
可写成
σ ? c = Σ i ( C i ε ? n i i ε ? m i i + y i , 0 ) φ i ? ? ? = σ ? i = C i ε ? n i i ε ? m i i + y i , 0 ? ? ? ( 7 ) y i , 0 = C i ε ? n i i ε ˙ m i i + y ′ i , 0 ? ? ? ( 8 )
σ
?
c
=
Σ
i
(
C
i
ε
?
i
n
i
ε
?
i
m
i
+
y
i
,
0
)
φ
i
?
?
?
=
σ
?
i
=
C
i
ε
?
i
n
i
ε
?
i
m
i
+
y
i
,
0
?
?
?
(
7
)
y
i
,
0
=
C
i
ε
?
i
n
i
ε
˙
i
m
i
+
y
′
i
,
0
?
?
?
(
8
)
式中 符号下标c 和i分别表示微元体和不同组元; C、 n和m为常数;y′i,0 为组元粉末未球磨时的屈服应力; yi, 0 为组元经不同时间球磨后的屈服应力; φi 为组元体积分数, %。 值得注意的是, 式(7)与式(8)中
ε ? i
ε
?
i
描述的对象并不相同, 前者指球磨后续挤压变形时产生的应变; 而后者指球磨过程中组元累积应变, 是一个与球磨时间历史和球磨变形历史相关的量, 由式(6)计算。
对于二元体系, 且组元变形时应变速率对应力的影响可以忽略时, 即mi =0, 球磨粉体本构方程则可表述为
y 1 , 0 ≤ σ ? c ≤ y 2 , 0
y
1
,
0
≤
σ
?
c
≤
y
2
,
0
时,
ε ? c = φ 1 ( σ ? c ? y 1 , 0 ) 1 n 1 C 1 ? ? ? ( 9 ) σ ? c ≥ y 2 , 0 ≥ y 1 , 0 时 ?
ε
?
c
=
φ
1
(
σ
?
c
-
y
1
,
0
)
1
n
1
C
1
?
?
?
(
9
)
σ
?
c
≥
y
2
,
0
≥
y
1
,
0
时
?
ε ? c = φ 1 ( σ ? c ? y 1 , 0 ) 1 n 1 C 1 + ( 1 ? φ 1 ) ( σ ? c ? y 2 , 0 ) C 2 ) 1 n 2 ? ? ? ( 1 0 )
ε
?
c
=
φ
1
(
σ
?
c
-
y
1
,
0
)
1
n
1
C
1
+
(
1
-
φ
1
)
(
σ
?
c
-
y
2
,
0
)
C
2
)
1
n
2
?
?
?
(
1
0
)
这是一个与组元粉体力学性能、 组元体积分数以及组元变形历史相关的全致密微元体应力—应变表达式。
为了客观反映二元复合粉体挤压过程, 在借助有限元软件分析时, 尚需为其提供全致密粉体材料传热系数和比热容等材料物性参数。 这些参数均可按复合原则处理, 球磨粉体传热系数κ c 和热容C c 计算方法如下:
κ c = ∑ i λ i φ i ? ? ? ( 1 1 ) C c = ∑ i C i φ i ? ? ? ( 1 2 )
κ
c
=
∑
i
λ
i
φ
i
?
?
?
(
1
1
)
C
c
=
∑
i
C
i
φ
i
?
?
?
(
1
2
)
上述推导的球磨粉末致密微元体的本构方程, 不仅适用于二元球磨粉末, 也可推广应用于一元和多元球磨粉末。 但由于该本构方程的推导是以粉末变形为基础的, 因此对于长时间球磨形成强烈固溶和生成新相的情况不宜选用。
2Ti/Al球磨粉末挤压有限元分析
2.1有限元模型
图1所示为模拟d 48 mm×24 mm球磨Ti/Al(摩尔比为1∶1)粉末坯挤压过程的有限元模型。 模具与工件之间的摩擦属剪切摩擦, 即摩擦力F =μτ , 式中τ 为粉体剪切屈服强度, μ 为摩擦因子(取0.3)。 粉末坯料的平均传热系数和热容分别为κ c ρ r 和C c ρ r , 均与实时相对密度ρ r 成正比。 坯料与模具之间的传热系数为50.0 kW/(m2 ·℃)。 粉坯和模具的初始温度都为250 ℃。 球磨不同时间后粉坯(室温,1 000 MPa下单向压制而成)挤压起始相对密度如表1所示。 球磨粉致密体的应力—应变关系由前面所推导的本构方程计算确定。 球磨条件为: 磨球材质和直径分别为GCr15和8 mm, 装球量为0.5, 球料质量比20∶1, 球磨机搅拌杆长100 mm, 球磨转速400 r/min。 本研究仅对球磨时间不超过3 h的粉末挤压过程进行数值分析。 因为更长时间的球磨会引发粉末强烈固溶, 甚至反应生成金属间化合物
[15 ]
, 已不宜选用上述本构方程。
图1 粉末坯料与模具的初始有限元网格(λ=9)
Fig.1 Finite element model of powder billet and die set
表1 球磨粉坯挤压起始相对密度
Table1 Raw relative densities of milled powder billets before extrusion
Milling time/h
Relative density/%
0
88.25
0.5
87.70
1
84.70
2
81.42
3
75.41
有限元数值模拟软件采用非线性有限元软件DEFORM2D。
2.2挤压过程数值分析
图2所示为对球磨1 h后粉末坯挤压过程(挤压比λ =9)模拟所获得的坯料内部等效应力分布图。 粉坯在挤压过程中, 最大等效应力出现在凹模口附近; 在挤出坯料的前端存在一个低应力; 与凸模端部接触部位的应力相对较低,只是当凸模行程D 达到最大值(48 mm)附近时, 与凸模端部接触部位的应力才达到较大值。
对球磨不同时间粉末坯料挤压(λ =9和16)过程的数值模拟结果显示, 经不同时间球磨粉末坯料在挤压过程中的应力响应规律以及挤压载荷随凸模压下量的变化规律与球磨1 h挤压比为9时的相似。
图2 球磨1 h粉末坯料挤压等效应力分布(λ=9)
Fig.2 Effective stress distribution during extrusion of billets of powder milled for 1 h
(a)—D=10 mm; (b)—D=20 mm; (c)—D=30 mm; (d)—D=40 mm
图3所示为数值模拟获得的球磨不同时间粉末坯料挤压过程中最大挤压载荷及其拟合曲线。 在同一挤压条件下, 随球磨时间的延长, 所需挤压力显著升高; 对于球磨相同时间粉末坯料, 随挤压比的增加, 挤压载荷亦明显升高。
图3 球磨不同时间粉末坯料挤压力及其拟合曲线
Fig.3 Extrusion pressures for billets of powder milled for various time and their fit curves
3实验验证
挤压实验选择与前述模拟初始条件相同的预热温度、 模具形状和尺寸, 并采用石墨和黄油混合物作润滑介质。
图4所示为经不同时间球磨粉末坯料在不同挤压比下的挤压力。 可以看出: 随球磨时间的延长, 由于球磨累积应变增加, 故挤压力增加; 同时随挤压比的增加, 由于粉体挤压变形程度增加, 挤压力也相应提高。 结合表1还可推断, 随球磨时间的延长, 粉末坯料的固结致密难度相应增大。
图4 球磨时间对粉末坯料挤压力的影响
Fig.4 Extrusion pressure as function of milling time
综上所述, 实验获得的挤压规律与有限元数值模拟结果是一致的, 在数值上两者是基本吻合的, 只是实验挤压力略偏高。 分析认为这主要是由以下原因所造成: 1) 模拟时选用的材质为纯Ti和纯Al, 而实际上粉末因含杂质和表面氧化等因素, 球磨前的真实屈服应力要偏高; 2) 模拟时未考虑实际球磨粉末组元间因存在一定程度的固溶而引起的强化效应; 3) 模拟选用的Ti、 Al应力—应变公式外推至大塑性变形时的精确性尚有待进一步证实。
4结论
1) 通过对机械球磨过程数值描述, 推导出了复合粉末组元层片厚度计算公式, 并结合细观复合力学理论给出了球磨粉末连续体的本构方程, 为球磨粉末塑性变形数值分析提供了基础。
3) 高能球磨Ti/Al二元粉末坯料挤压过程的数值分析表明, 随着球磨时间的延长, 粉末固结致密难度增加, 挤压力显著增大。 该趋势和模拟数值均与挤压实验结果基本一致, 证实了所推导的球磨粉末本构方程是合理的。
参考文献
[1] BenjaminJS .Dispersionstrengthenedsuperalloysbyme chanicalalloying[J].MetTrans,1970,8(1):29432951.
[2] 朱心昆,林秋实,陈铁力,等.机械合金化的研究及进展[J].粉末冶金技术,1999,17(4):291296.ZHUXin kun,LINQiu shi,CHENTie li,etal.Re searchand progressinmechanicalalloying[J].PowderMetallurgyTechnology,1999,17(4):291296.
[3] 唐振方,高 勇,杨元政,等.WC粉体的高能球磨细化[J].中国有色金属学报,2000,10(S1):246249.TANGZhen fang,GAOYong,YANGYuan zheng,etal.Triturationofmicrogradedtungstencarbidebyhighenergyballmilling[J].TheChineseJournalofNonfer rousMetals,2000,10(S1):246249.
[4] ZouY ,SajiS ,KusabirakiK .SolidstatereactionsandmicrostructuralevolutionofAl90-xFexNi10 mixedpow dersduringhigh energyballmilling[J].InternationalJournalofMaterialsandProductTechnology,2001(2):682687.
[5] HwangS ,NishimuraC ,McCormickPG .Mechanicalmillingofmagnesium powder[J].MaterialsScienceandEngineeringA ,2001(318):2233.
[6] 张厚安,刘心宇,陈 平.Al/MoSi2复合粉末材料的机械合金化合成[J].粉末冶金技术,2002,20(3):166168.ZHANGHou an,LIUXin yu,CHENPing.SynthesisofAl/MoSi2 powdersbymechanicalalloying[J].PowderMetallurgyTechnology,2002,20(3):166168.
[7] LiangGX ,LiZM ,WangED .Hothydrostaticextru sionofmechanicallyalloyedAl4.9Fe4.9Nipowders[J].JournalofMaterialsScienceandTechnology,1994,10(4):285288.
[8] AhnIS ,KimBS ,MoonLH .Densificationbehaviorofhighenergyball milledW 20%Cu[J].AdvancesinPow derMetallurgyandParticulateMaterials,1995,5(1):209217.
[9] SmithLN ,MidhaPS .Computersimulationofmor phologyandpackingbehaviourofirregularparticles,forpredictingapparentpowderdensities[J].ComputationalMaterialsScience,1997,7(4):377383.
[10] LiuL ,YuanY .Dynamicsimulationofpowdercompactbyrandompackingofmonosizedandpolydisperseparti cles[J].JournalofMaterialsScienceLetters,2000,19(10):841843.
[11] HuangCC ,ChengJH .Forgingsimulationofsinteredpowdercompactsundervariousfrictionalconditions[J].InternationalJournalofMechanicalSciences,2002,44(3):489507.
[12] MauriceDR ,CourtneyTH .Thephysicsofmechanicalalloying:afirstreport[J].MetTransA ,1990,21:289295.
[13] 杨君友,吴建生,曾振鹏.机械合金化过程中粉末的变形及其能量的转化[J].金属学报,1998,34(10):10611067.YANGJun you,WUJian sheng,ZENGZhen peng.Investigationonthedeformationandenergytransferofpowdersduringmechanicalalloying[J].ActaMetallur gicaSinica,1998,34(10):10611067.
[14] WatanabeR ,HashimotoH ,LeeGG .Computersimu lationofmillingballmotioninmechanicalalloying[J].MaterialsTransactionsJIM ,1995,36(2):102109.
[15] 李小强,胡连喜,王尔德.TiAl二元粉末机械球磨过程中组织结构的变化[J].稀有金属材料与工程,2001,30(3):183186.LIXiao qiang,HULian xi,WANGEr de.EvolutionofmicrostructuresofTi/Albinarypowdersduringme chanicalalloyingprocess[J].RareMetalMaterialsandEngineering,2001,30(3):183186.