DOI:10.19476/j.ysxb.1004.0609.2000.02.016
六方结构金属板材弹性的计算法
王超群 商顺利 王宁
北京有色金属研究总院!北京100088
摘 要：
具有六方结构的金属板材呈现较强的物理和力学性能各向异性 , 尤其是钛及其合金呈现的弹性各向异性 , 其轧向与横向的弹性模量差异较大。根据Bunge等人的三维取向分布函数 (ODF) , 由单晶弹性系数和少量有限的织构系数如直至四阶的C系数 , 简便地计算了低对称六方结构的板材的弹性。钛及其合金板材弹性的理论计算结果与实测数值相吻合 , 钛合金板材的弹性各向异性是由织构产生的 , 具有 { 1 2 1 5}〈1 0 1 0〉织构的钛板材的横向弹性模量为 1 1 6GPa, 大于轧向的弹性模量 (1 0 8GPa)
关键词：
六方金属板材 ;弹性 ;织构 ;极图 ;ODF ;
中图分类号： TG113
收稿日期： 1999-01-19
Calculation of elastic properties of metals and alloys with a hexagonal structure
Abstract：
The presence of texture was known to cause anisotropy of elastic properties of polycrystalline materials with a hexagonal structure. By using texture coefficients and elastic constant of single crystal, the elastic property of polycrystalline plate can be calculated. The present calculation method has been applied to predicting elastic modulus for a rolled titanium sheet and the results are consistent well with experimental ones.
Keyword：
hexagonal metallic plate; elasticity; texture; pole figure; ODFs;
Received： 1999-01-19
具有低对称六方结构的金属晶体, 其物理和力学性能呈强的各向异性, 如弹性模量绕C 轴呈对称分布, C 轴方向的弹性模量 (143.3 GPa) 远大于底面各方向的弹性模量 (104.4 GPa) 。 因此只要建立弹性模量E 的方向分布与C 轴织构的方向分布之间的关系, 通过板材的织构测定, 就可以计算合金板材的弹性模量, 并可通过合金化、 加工与热处理工艺的调整, 有目的地控制织构与弹性各向异性, 以满足实际使用的要求。 人们一般采用反极图或极图法, 由单晶弹性数据, 结合极密度加权平均, 计算多晶金属板材的弹性
[1 ]
, 然而反极图与极图法均为二维的, 且在织构测量时需要叠片制样, 既费时又不精确。 为此, Bunge
[2 ]
, Morris
[3 ]
和ZUO
[4 ]
等人发展了一种以三维取向分布函数 (ODF) 为基础的新方法, 根据单晶弹性数据和少量有限的织构系数 (至四阶的C 系数) , 便可简捷地计算六方织构板材的弹性向各异性。 其中根据Voigt-Reuss-Hill等的平均法, 结合织构加权函数, 对含磷深冲钢板以及六方Zn板的弹性计算已取得很好的结果, 但对钛合金板材尚未见报道。 为此本文将主要采用ZUO
[4 ]
的方法计算织构钛合金板材的弹性模量, 并讨论了几种方法的一致性。
1金属板材弹性各向异性的计算方法
1.1多晶材料的弹性性质
上面已扼要介绍了在一定的假设条件下, 通过相应的平均方法, 能够从单晶的数据计算多晶材料的弹性性质。 Voigt平均法系假设多晶体内的所有晶粒的应变增量相等, 而Reuss平均法则假设在整个试样中应力均匀分布。 Hill
[5 ,6 ]
认为这两种模型均不真实, 因为在恒应变模型中, 相邻晶粒之间的应力难以平衡, 而在恒应力模型中, 晶粒的畸变在晶界上不能连续。 Hill指出根据Voigt研究的弹性模量要大于根据Reuss所得的结果, 正确的数值应处于两者之间。 Bunge
[7 ]
也指出在单轴拉伸下, Voigt近似对于铅笔型晶粒组成的多晶体结构是正确的, 而Reuss近似则对于片状晶组成的多晶体是合理的。 因此真实的计算模型应是Hill近似, 即由晶粒形状因子m 作为加权值, 对于等轴晶的多晶体m =0.5。 Bunge将Voigt-Reuss-Hill近似法推广应用于具有立方晶系数的织构多晶材料。 Morris
[3 ]
则发展了一种利用织构加权函数的平均四阶张量法, 该法适用于正交晶系和更高的晶系 (四方、 六方和立方系) 多晶板材的弹性计算。 近年来Bunge
[2 ]
, ZUO
[4 ]
和LI
[8 ]
等人也就六方结构金属的弹性进行了计算, 其中ZUO, Morris等人对多晶弹性的Voigt-Reuss-Hill近似法给出了普遍解
[4 ,9 ,10 ]
。 本文将主要介绍这种方法并将其用于钛及其合金板材的弹性测算。
1.2六方结构金属板材的弹性计算方法原理 [4]
作为一级近似, 多晶材料弹性的平均值
E
?
i
j
k
l
=
∮
E
i
j
k
l
?
f
(
g
)
d
g
=
Τ
?
i
j
k
l
o
p
q
r
·E opqr (1)
式中
Τ
?
i
j
k
l
o
p
q
r
为织构因子, E opqr 为四阶弹性张量分量, f (g ) 为晶体取向分布函数, 其中E ijkl
对于六方晶系, 考虑到晶体对称性以及应力、 应变张量的对称性, E opqr 只有5个独立分量即E 1111 , E 1122 , E 1133 , E 2323 和E 3333 , 因此
E
?
i
j
k
l
=
Τ
?
i
j
k
l
1
1
1
1
?
E
1
1
1
1
+
Τ
?
i
j
k
l
1
1
2
2
?
E
1
1
2
2
+
Τ
?
i
j
k
l
1
1
3
3
?
E
1
1
3
3
+
Τ
?
i
j
k
l
3
3
3
3
?
E
3
3
3
3
+
Τ
?
i
j
k
l
2
3
2
3
?
E
2
3
2
3
?
?
?
(
2
)
式中
Τ
?
i
j
k
l
o
p
q
r
是对四阶弹性张量E opqr 为独立的T ijkl opqr
Τ
?
i
j
k
l
o
p
q
r
独立, 它们的系数e
L
′
i
j
k
l
o
p
q
r
n
′
=a L ′n ′ijkl opqrib L ′n ′ijkl opqrT uvst
Τ
?
i
j
k
l
o
p
q
r
, 线性独立T uvst
[
4 ]
。
对于板材试样, 当拉伸坐标系 (O -xyz ) 相对于试样坐标系 (通常取轧向-横向-板面法向) 的取向为{α , 0, 0}, 则
Τ
?
u
v
s
t
{
α
?
0
?
0
}
=
∑
l
=
0
(
2
)
4
∑
n
=
-
L
+
L
[
4
π
(
2
L
+
1
)
]
-
1
/
2
?
C
L
1
n
e
L
*
n
(
u
v
s
t
)
?
exp
(
i
n
α
)
?
?
?
(
3
)
由式 (3) 可见, 织构因子的计算和平均弹性张量的计算只需要测算至四阶的织构系数C 1n L e
L
*
n
是纯数学量, aL′n′ , bL′n′ 可查表1。
实际计算如下:
例如, 根据Reuss 平均法,
S
?
1
1
1
1
R
的计算如式 (4) 。
S
?
1
1
1
1
R
=
Τ
?
1
1
1
1
S
1
1
1
1
+
Τ
?
1
1
1
2
S
1
1
2
2
+
Τ
?
1
1
1
3
S
1
1
3
3
+
Τ
?
3
3
3
3
S
3
3
3
3
+
Τ
?
1
1
4
4
S
2
3
2
3
?
?
?
(
4
)
根据与T uvst
[4 ]
, 有
Τ
1
1
1
2
=
0
,
Τ
1
1
1
3
=
2
Τ
1
2
1
1
+
2
Τ
1
3
1
1
?
Τ
1
1
3
3
=
1
-
Τ
1
1
1
1
-
2
Τ
1
2
1
1
-
2
Τ
1
3
1
1
?
Τ
1
1
4
4
=
4
Τ
1
2
1
1
+
4
Τ
1
3
1
1
}
?
?
?
(
5
)
因此式 (4) 变为
S
?
1
1
1
1
R
=
Τ
?
1
1
1
1
S
1
1
1
1
+
Τ
?
1
1
1
3
S
1
1
3
3
+
Τ
?
1
1
3
3
S
3
3
3
3
+
Τ
?
1
1
4
4
S
2
3
2
3
?
?
?
(
6
)
其中只有织构因子
Τ
?
1
1
1
1
?
Τ
?
1
2
1
1
及
Τ
?
1
3
1
1
待算。 它们的计算式可写成
Τ
?
u
v
s
t
{
α
?
0
?
0
}
=
a
0
0
+
1
5
π
(
a
2
0
C
2
1
0
+
a
2
2
C
2
1
2
cos
?
2
α
)
+
1
6
π
(
a
4
0
C
4
1
0
+
a
4
2
C
4
1
2
cos
2α+a44 C
4
1
4
cos 4α) (7)
式中
a
0
0
=
(
2
π
)
-
1
?
a
0
0
=
8
/
1
5
(对于
Τ
?
1
1
1
1
)
, 而对于其它的织构因子, T 1211 及T 1311 的a 00 , a 20 , a 22 , a 40 , a 42 和a 44 可查表1。
表1 六方晶体对称的级数展开系数 [4]
Table 1 Series expansion coefficients e
L
′
n
′
(uvst ) =α L ′n ′uvst ) +ib L ′n ′uvst ) for hexagonal crystal symmetry
L ′n ′uvst ) a L ′n ′uvst ) b L ′n ′uvst )
B 0 ×20B 1 ×
21B 2 ×
22B 2 ×
40B 0 ×
41B 1 ×
42B 3 ×
43B 4 ×
44B 5 ×
B 2 ×22B 2 ×
41B 1 ×
42B 3 ×
43B 4 ×
44B 5 ×
56
16
-8
3
-2
1
7
-4
1
-1
35
28
7
2
1
-4
2
1
-1
1
4
3
-1
4
1
-1
56
16
8
3
2
1
7
2
-1
-4
-2
-4
-3
-1
-1
1
1
4
1
1
-4
4
4
-4
-1
-2
B 0 =2 (π) 1/2 /105; B 1 = (5π) 1/2 /105; B 2 = (30π) 1/2 /105; B 3 = (10π) 1/2 /105; B 4 = (35π) 1/2 /105; B 5 = (70π) 1/2 /105
对于板材平面里的弹性模量
E
?
(
α
)
的变化可按下式计算:
E
?
(
α
)
=
1
S
?
(
1
1
1
1
)
R
{
α
,
0
,
0
}
?
?
?
(
8
)
对于Voigt平均法也可以进行类似的计算, 但E opqr 用劲度张量C opqr 替换。 而Hill平均法可用下式计算:
C
?
i
j
k
l
Η
=
[
C
?
i
j
k
l
V
+
(
S
?
i
j
k
l
R
)
-
1
]
/
2
S
?
i
j
k
l
R
=
[
(
C
?
i
j
k
l
V
)
-
1
+
S
?
i
j
k
l
R
]
/
2
?
?
?
(
9
)
式中
(
S
?
i
j
k
l
R
)
-
1
,
(
C
?
i
j
k
l
V
)
-
1
分别是
S
?
i
j
k
l
R
?
C
?
i
j
k
l
V
的逆阵。
2实验方法
选用纯钛和高弹钛合金板材。 在X′pert MRD衍射仪上, 用Cu靶, 40 kV, 40 mA, 石墨单色器和正比计数管, 测量了 (010) , (002) , (110) , (011) 和 (012) 极图, 步进Δα =Δβ =5 °, α 扫测范围0 °～85 °, β 范围为0 °～360 °, ODF计算的展开度L =22。
计算弹性模量时选用C
2
1
0
, C
2
1
2
, C
4
1
0
, C
4
1
2
和C
4
1
4
五个C 系数, 而钛单晶的弹性常数用C 1111 =162.4 GPa, C 1122 =92.0 GPa, C 1133 =69 GPa, C 3333 =180.7 GPa和C 2323 =46.7 GPa; 并由式 (4) , (7) , (8) 和 (9) 计算Voigt-Reuss-Hill平均法的E (α ) 值。
3结果与讨论
图1是热轧钛计算的ODF。 将由ODF获得的C 1n 4 和钛单晶的弹性常数代入式 (4) , (7) , (8) 和 (9) , 计算了Voigt, Reuss和Hill平均法的E (α ) 值, 其结果如图2所示。 由图2可见, 横向的弹性模量 (≈116 GPa) 大于轧向的弹性模量 (≈108 GPa) , 这与图1/图2的极图数据和ODF分析结果相吻合。 由极图分析可知该板材的主要织构为 (0001) ±30
°
Ν
D
-
Τ
D
?
1
0
1
?
0
?
, 而由ODF分析可知主要的织构类型为
{
1
2
1
?
5
}
?
1
0
1
?
0
?
, 因此轧向的杨氏模量偏低是可以理解的。 轧向与横向弹性模量实际测量结果分别为107.7 GPa和116 GPa, 实验测试结果与理论预测相吻合。
无论是ZUO
[4 ]
, 还是Bunge
[2 ]
的方法, 其计算方法是一致的。 例如Bunge
[2 ]
导出的Reuss平均法的表达式为
E
?
R
(
α
)
=
[
S
?
R
(
α
)
]
-
1
?
?
?
(
1
0
)
S
?
R
(
α
)
=
S
0
/
1
5
+
S
2
?
F
2
(
α
)
+
S
4
?
F
4
(
α
)
?
?
?
(
1
1
)
S
0
=
-
8
S
1
1
+
3
S
3
3
+
4
S
1
3
+
2
S
4
4
?
?
?
(
1
2
)
S
2
=
-
8
S
1
1
+
6
S
3
3
+
2
S
1
3
+
S
4
4
?
?
?
(
1
3
)
S
4
=
S
1
1
+
S
3
3
-
S
1
3
-
S
4
4
?
?
?
(
1
4
)
将式 (12) , (13) 和 (14) 代入式 (11) 可得 (用4位数字标记)
S
?
R
(
α
)
=
8
1
5
S
1
1
1
1
+
3
1
5
S
3
3
3
3
+
4
1
5
S
1
1
3
3
+
8
1
5
S
1
3
1
3
+
(
-
8
S
1
1
1
1
+
6
S
3
3
3
3
+
2
S
1
1
3
3
+
4
S
1
3
1
3
)
?
F
2
(
α
)
+
(
S
1
1
1
1
+
S
3
3
3
3
-
S
1
1
3
3
-
4
S
1
3
1
3
)
F
4
(
α
)
?
?
?
(
1
5
)
而由式 (4) ZUO
[4 ]
的结果, 即可以得
S
?
R
(
α
)
=
8
1
5
S
1
1
1
1
+
3
1
5
S
3
3
3
3
+
4
1
5
S
1
1
3
3
+
8
1
5
S
1
3
1
3
+
?
?
?
?
(
1
6
)
因此两者是一致的, 而由Bunge的F 2 及F 4 计算式:
F
2
(
α
)
=
4
5
2
5
2
(
1
2
C
2
1
1
Ρ
?
2
0
+
C
2
1
2
Ρ
?
2
2
cos
2α) (17)
图1 热轧工业钛板的ODFs (φ2截面, L=22)
Fig.1 ODFs of hot-rolled titanium plate (φ 2 section, L =22)
图2 热轧工业钛板弹性模量的变化 (按照Voigt, Reuss和Hill平均法计算)
Fig.2 Elastic modulus of hot rolled titanium plate (calculated according to V-R-S average) ○—Experimental values
?
?
F
4
(
α
)
=
1
6
9
4
5
[
1
2
C
2
1
1
Ρ
?
4
0
+
C
2
1
2
Ρ
?
4
2
cos
?
2
α
+
C
2
1
4
Ρ
?
4
4
cos
?
4
α
)
?
?
?
(
1
8
)
式中 C
4
1
n
为C 系数,
Ρ
?
2
n
为连带勒让德多项式的傅氏系数, 可以进一步同ZUO的方法比较, 得出:
Ρ
?
2
0
=
-
1
.
2
5
?
Ρ
?
2
2
=
2
.
1
6
5
1
?
Ρ
?
4
0
=
0
.
7
9
5
5
?
Ρ
?
4
2
=
-
1
.
1
8
5
9和
Ρ
?
4
4
=
1
.
5
6
8
7
。
此外还可以进一步比较, 证明ZUO
[4 ]
, LI
[8 ]
和Bunge方法之间的一致性, 这里不再一一赘述。
参考文献
[1 ] CHENHong sun (陈洪荪 ) .ElasticanisotropyofMetalMaterials (金属的弹性各向异性 ) [M] .Beijing :Metal lurgicalIndustryPress .1 996 .
[] [2 ] BungeHJandParkNJ.ProgramSystemforPhysicalPropertiesofTexturedMaterials [M] .CuvillierVerlag , G ttingen , 1 993 .
[3] MorrisPR .Averagingfourth rankstensorswithweightfunction [J] .JApplPhys, 1 969, 40 (2 ) :447.
[4] ZUOLiang (左 良 ) , HunbertMandEslingC .ElasticpropertiesofpolycrystalsintheVoigt Reuss Hillapproxi mation [J] .JApplCryst, 1 992 , 2 5 :751 .
[5] HillR .Theelasticbehaviourofacrystallineaggregate[J] .ProcPhysSoc , 1 952 , A65 :349.
[6] HuHsun (胡 甸阝 ) .高法向塑性的各向异性冷轧和退火磷钢板材的弹性性质 [J] .ActaMetallSinica, 1 981 , 1 7 (6) :595 .
[7] BungeHJ.Partialtextureanalysis [J] .TexturesandMicostructures, 1 990 , 1 2 :47.
[8] LIYandThompsonRB .Relationsbetweenelasticcons tantsCijandtextureparametersforhexagonalmaterials[J] .JApplPhys, 1 990 , 67 (5) :2 663 .
[9] MorrisPR .Elasticconstantsofpolycrystals [J] .IntJEngSci, 1 970 , 8:49.
[1 0 ] ZUOLiang (左 良 ) , XUJia zen (徐家桢 ) andLIANGZhi de (梁志德 ) .Averagingfourth rankselas tictensorsfortexturedpolycrystallineaggregateswith outphysicalsymmetry [J] .JApplPhys, 1 989, 66 (6) :2 338
