对于双组元体系,组元i的堆积密度为αi,体系堆积密度为ρ,干扰粒子的体积占整个体系的体积分数为2,粒度为r2,在干扰粒子引入前,体系密堆的那一组分其颗粒体积占整个体系的体积分数为1,粒度为r1,1 2=ρ,ηi组元i的颗粒体积占体系颗粒体积的百分数,η1 η2=1。 当r≤μ1r1时,体系为非挤塞体系,即体系为填隙机制所控制,此时体系的堆积密度可写为:ρ=α1 2(1)当μ2r1≤r2≤r1/μ2时,系体为替换机制控制,此时,体系的堆积密度不受干扰粒子的影响,恒定为一常数,即:ρ=α1=α2=α(2)当r1/μ2 r2时,ρ=α1 (1-α1)g(2,1)2(3)考虑到r2→r1/μ2时,ρ=α=α1,g(2,1)修改为:g(2,1)=1-r1/(μ2r2)(4)当μ1r1≤r2≤r1μ2时,随着r2的增大,体系从填隙机制逐渐向替换机制转变。此时,体系堆积密度为:ρ=α1/(1-η2f(r1r2))=α1/(1-η2(1-α1/αT(r1r2)(1-μ13r13/r23)))(5)其中η2为干扰粒子的固体体积分数。由(5)式可得:f(r1r2)=1-α1/αT(r1r2)(1-μ13r13/r23)(6)当r2=μ2r1时,即:1-α1/αT(r1r2)(1-μ13r13/r23)=0(7)由(7)式可得:αT(r1r2)=α(1-μ13/μ23)(8)因此,T(r1,r2)的数学模型由下式表示:T(r1,r2)=α(r1)(μ2r1-r2) α(1-μ31/μ32)(9)由于limr2/r1→μ2(1-α)g(r1,r2)=limr2/r1→μ2f(r1,r2),所以有:limr2/r1→μ2(1-α)(1-r2/(μ2r1))=limr2/r1→μ21-α/αT(r2r1)(1-μ21r31/r32)(10)考虑到(4),(9)两式,对(10)式求解后,可得:α(r1)=αr1((1-α)μ2(1-μ31μ32)-3μ31μ-42)(11)因此,f(r1,r2)修改后变为: f(r1,r2)=1-α(1-μ31r31/r32)α(r1)(μ2r1-r2) α(1-μ31/μ32)(12)α(r1)由(9)式表示。当r2 r1时,ρ=α1 (1-α1)2(13)此时,体系为非挤塞(non-crowding)体系。综上所述,体系的堆积密度随干扰粒子的增大可概括如下:当r2 r1 ρ=1 2当μ1r1≤r1=μ2r1 ρ=α1/(1-η2f(r1r2)),f(r1r2)由(12)式表示。当μ2r1≤r2≤r1/r2时,ρ=α1=α2=α当r1/μ2≤r1时,ρ=α1 (1-α1)g(2,1)2 g(2,1)由(4)式表示。当r2 r1时,ρ=α1 (1-α1)2
4 数据验证 当r2/r1=0.4时,体系同时受填隙机制和替换机制作用,由修改前后两模型计算出的堆积密度值与实验值[2]的比较如表1所示。从表1中可以看出,修改后的模型精确度得到提高。 表1 修改前后的模型计算与实验数据[2]的比较(r2/r1=0.4) 当r2/r1≥0.741时,此时体系为替换机制控制,修改后的模型与实测数据为同一恒定值,表2列出了两模型的比较。由于模型修改正是基于当粒度比大于等于0.741时,由实验观察到体系堆积密度与单组元的堆积密度相同。此时,体系的堆积密度可由单组元的堆积密度表示,因而,当体系受替换机制控制时,模型修改后的计算值与实验测量值的偏差为0。 表2 修改前后的模型计算与实验数据[2]的比较
5 结论 线性堆积密度模型是现有粉末颗粒堆积模型中,预测堆积密度准确度较高的模型。但是由于忽略了替换机制,该模型在预测粒度比大于临界值μ2的粉末体系的堆积密度缺乏精确度。改进后的模型同时考虑了替换机制和填隙机制,经过数据验证表明,模型的精确度得到较大的提高。
|