文章编号：1004-0609(2014)06-1449-10

基于均匀化方法的钨丝增强锆基块体非晶复合材料等效弹性常数预测

廖光开，李乡安，邹  萍，陈舒敏，龙志林

(湘潭大学 土木工程与力学学院，湘潭 411105)

摘  要：基于钨丝增强锆基块体非晶复合材料结构特点及均匀化方法，建立计算该复合材料等效弹性常数的理论模型；结合有限单元法对钨丝增强锆基块体非晶复合材料等效弹性常数进行数值计算；分析增强相的体积分数或尺寸对钨丝增强锆基块体非晶复合材料有效性能的影响。研究结果表明：将均匀化理论与有限元方法相结合能有效地预测具有周期性细观结构的块体非晶复合材料的有效弹性性能，为合理设计该类材料奠定基础。

关键词：块体非晶复合材料；均匀化理论；有限元分析；增强相；弹性常数

中图分类号：TG139           　   　     文献标志码：A

Homogenization-based approach for predicting equivalent elastic constants of tungsten fiber reinforced Zr-based bulk metallic glass matrix composites

LIAO Guang-kai, LI Xiang-an, ZOU Ping, CHEN Shu-min, LONG Zhi-lin

(College of Civil Engineering and Mechanics, Xiangtan University, Xiangtan 411105, China)

Abstract: Based on the structural signature of tungsten-fiber reinforced Zr-based bulk metallic glass matrix composites (BMGMCs) and homogenization theory, a theoretical model for calculating their equivalent elastic constants was established. The equivalent elastic constants of these composites were predicted by the proposed model combined with the finite element analysis. The effects of the volume fraction and size of the reinforced phases on the effective properties of tungsten-fiber reinforced BMGMCs were further studied. The results demonstrate that a combined approach of the homogenization theory and the finite element formulation can effectively predict the elastic properties of tungsten-fiber reinforced BMGMCs with periodic  meso-structure and thus, providing bases for the rational design of this kind of composite.

Key words: bulk metallic-glass matrix composite; homogenization theory; finite element analysis; reinforced phase; elastic constant

钨丝增强锆基块体非晶复合材料作为一种新型高级的复合材料，具有高强、高弹、耐磨和耐蚀等特点^[1-3]。由于钨丝的引入使锆基块体非晶合金中局部软化的剪切带^[4-5]得到加固，从而极大地提高了这种复合材料的塑性，并具有良好的自锐性^[6-8]。因此，钨丝增强锆基块体非晶复合材料在军工、航空、体育等领域得到广泛应用。近年来，国内外学者对该种新型材料的力学性能进行了大量实验研究(如刚度、强度、断裂、损伤等)^[1]，但现阶段对其弹性常数的预测还鲜有文献报道。建立合适的等效弹性计算模型也许有益于理性设计钨丝增强锆基块体非晶复合材料。

钨丝增强锆基块体非晶复合材料由块体非晶基体和钨纤维组合而成，其几何细观结构呈周期性，对其力学性能的研究可以尝试选择具有代表性的细观单元进行分析，预测其等效弹性性能。目前，对复合材料等效弹性常数的预测方法和模型很多，如自洽理论^[9-10]、Mori-Tanaka方法^[11]、自治模型^[12]等。这几种模型在分析之前都进行了较多的假设，其分析结果和材料的实际情况存在较大的偏差，具有一定的局限性，不能精确地描述材料的应力应变场。随着国内外学者多年来研究的深入，综合考察了纤维取向、界面和增强相的几何结构等多方面因素对弹性性能的影响，HASSANI等^[13-14]成功地运用均匀化理论模型预测或计算了具有细观周期性结构复合材料的弹性常数。同时均匀化方法是一种具有严格理论基础的数学方法。根据复合材料细观周期性的特点，通过多尺度小参数渐近展开的摄动方法，建立有关位移、应变和应力场的摄动方程，通过求解这些方程可分别得到宏观和细观两个尺度的相关力学性能参数，既能从细观尺度分析材料的等效模量和变形，又能从宏观尺度分析结构的响应^[15]，所以均匀化方法被广泛应用于复合材料力学性能分析^[16-17]。

本文作者基于文献[18]中报道的钨丝增强锆基块体非晶复合材料(W/Zr_41.2Ti_13.8Cu_12.5Ni₁₀Be_22.5)的结构特点，利用均匀化方法建立计算该复合材料等效弹性常数的理论模型，并结合有限单元法建立了计算钨丝增强锆基块体非晶复合材料等效弹性常数三维有限格式和完整的宏观、细观分析步骤；用Visual C⁺⁺编写均匀化理论的计算程序Homo3D，研究增强相的体积分数及其尺寸对钨丝增强锆基块体非晶复合材料有效性能的影响。

1  均匀化过程及其有限元格式

1.1  渐近均匀化过程

如图1(a)所示，设钨丝增强锆基块体非晶合金材料在三维空间中占据区域为，其细观结构可看成是非均质单元在空间的周期性重复排列。设该物体受体积力f、力边界Γ_t上的表面力t以及位移边界Γ_v上的给定位移的作用。当该物体处于静力平衡状态时，应满足弹性力学的基本方程和边界条件。由于复合材料具有非均匀结构(即钨丝增强相分散于基体中)，当涉及到钨丝增强相的细微结构时，用传统方法直接求解此边值问题异常困难。

在宏观某一点处的细观结构可看成是非均质单胞在空间中周期性重复排列，因此，宏观结构的性能参数是微结构的平均值。实际的非均质材料细观结构具有高度非均匀性和周期性，当宏观结构受外部载荷作用时，使得结构场变量(如位移和应力)在宏观位置x的非常小的邻域内也会有很大变化，并且呈周期性，即

           (1)

式中：x为宏观尺寸坐标；y为细观尺寸的坐标；为两种尺度之比；N为一整数；Y为周期。假设宏观的波动(x的波动)对场变量的影响较细观的(y的波动)要小得多，通过坐标变换，可以将宏观区域中的一点变换为一个具有周期性结构的单胞Y(见图1(b))。将单胞进行等效均匀化，当趋近于零时，等效均匀化后的单胞弹性力学常数反映了原复合材料等效均质的结果，由此可求得钨丝增强锆基块体非晶合金复合材料的等效弹性力学常数。

图1  具有细观周期性结构的钨丝增强锆基块体非晶复合材料结构及其单胞

Fig. 1  Illustration of structure of tungsten fiber reinforced Zr-based bulk metallic glass matrix composites with periodic microstructures

假设单胞为线弹性介质，在中，满足下列基本方程和边界条件：

平衡方程，

                                  (2)

几何方程，

                            (3)

本构方程，

                                (4)

力学边界条件，

                                    (5)

位移边界条件，

                                    (6)

式中：为材料的弹性张量；为位移张量；为应力张量；为应变张量；为力边界上的单位法向量；上标表示该函数具有两尺度的特征。

将宏观尺度下的位移场展开成关于小参数的渐近展开式：

                 (7)

式(7)中展开的各项表征了位移局部振荡的精细程度，随着展开项数的增加，式(7)逐渐趋近于真实位移场在微细观尺度下的振荡情况。根据线弹性范围的虚位移原理有^[19]

     (8)

式中：v_i为满足位移边界条件的虚位移。由及链式法得

                  (9)

将式(7)和(9)代入式(8)，得

                      (10)

根据Y周期函数的性质^[20]得

             (11)

结合式(10)和式(11)得

        (12)

                     (13)

            (14)

式(13)中：p_i为作用在增强相边界上的力。

根据周期性边界条件，对式(12)和(13)进行相关数学变换可求得

                             (15)

        (16)

式中：和分别是式(17)和(18)的解

      (17)

         (18)

设，代入式(14)得

     (19)

将式(16)代入式(19)得

      (20)

根据式(20)和虚位移原理得胞元宏观等效弹性张量：

           (21)

式中：为单胞域上的位移场；k、l、p和q为相应的张量指标符号。令

                    (22)

                          (23)

将式(20)改写为

          (24)

式中：和b_i分别为胞元内的平均残余应力和平均体力张量。式(24)即为均质化问题在宏观坐标系统中构建的整体虚位移方程。

通过求解可将细观非均质问题转化为均质问题，借助数值算法可求得解析解，根据式(20)和(21)，结合有限元算法，推导适合于钨丝增强锆基块体非晶复合材料的均匀化法的三维有限元格式。

1.2  均匀化理论的三维有限元计算格式

图2(a)所示为钨丝增强锆基块体非晶复合材料周期性结构的几何示意图。图中方框所示区域表示所选择的单胞结构。根据均匀化理论，结构的宏观等效弹性张量，必须通过求解式(17)单胞结构的位移得到，即利用周期性边界条件，通过求解六次方程组可分别求得：，，，，，，然后通过回代到式(21)中可得。采用六面体8节点等参单元对单胞区域进行离散，并利用相应的插值函数、及进行插值有：

单元内任意一点的坐标变换公式，

                       (25)

单元内任意一点的虚位移，

                       (26)

单元内任意一点的位移，

                    (27)

其中，节点m的形函数由局部坐标给出，，为母单元中的节点坐标，m为哑指标可取1, 2, …, 8，i为自由指标，可取1, 2, 3；表示单元中结点m的第个i坐标分量；表示单元中结点m的第i个虚位移分量；表示单元中结点m的第i个位移分量。根据偏微分法则及式(25)可知：

                     (28)

                   (29)

式中：称为Jacobian矩阵。定义应变矩阵为、单位阵及初始应变矩阵。

均质化问题式(24)可表示成矩阵形式：

                (30)

由虚位移v的任意性可得

                 (31)

即

                           (32)

式(32)即为求解均匀化有限单元附加参数的方程。[K]和分别为结构总体刚度矩阵和总体载荷矩阵，可分别通过单元刚度和单元载荷矩阵来集成。

(33)

 (34)

式中：为作用在单元边界上的力矩阵。

同理等效弹性模量式(21)的矩阵形式为

(35)

式中：为所有单元体积之和。式(33)~(35)的积分可通过三维的8节点高斯积分公式求得。因此，细观均匀化问题可化为具有初应力作用及周期性边界条件的空间有限元问题来求解，针对单胞上不同的初应变受载模式，采用相应的固定边界条件。

2  单胞的周期性边界条件

渐近均匀化理论要求细观单胞呈周期性排列，周期性边界条件施加在单胞相互平行的边界上，HASSANI等^[13-14]给出了周期性边界条件施加方式。在进行细观单胞有限元分析求解时，细观单胞内广义位移需要满足周期性边界条件，即

      (36)

这里，i=1, 2分别表示二维坐标的y₁和y₂两个方向。也就是说对于单胞边界上对称的两点A和B(见图2(a))，有。对于二维问题，细观均匀化问题(35)需解三次方程(对于三维问题需解六次方程)。该方法在处理位移和力边界条件方面有着很大的优越性，即只用到广义位移的周期性边界条件。对于一个具体的标准单胞来说，上述边界条件的施加将会根据单胞变形的形式而不同。单胞的变形方式有以下3种况：1) 在正方形单胞水平方向上作用单位初应变，如图2(b)所示，其周期性边界条件的施加图3(a)为取对称单胞1/4的状态；2) 在正方形单胞垂直方向上作用单位初应变，如图2(c)所示，其周期性边界条件的施加与第一种情况相同；3) 在正方形单胞上作用单位初剪应变，如图2(d)所示，其周期性边界条件的施加如图3(b)所示，取对称单胞1/4的状态。

3  分析与讨论

利用Visual C⁺⁺语言编写Homo3D均匀化程序来计算钨丝增强锆基块体非晶复合材料等效弹性常数。利用图4所示的有限元模型进行均匀化理论计算。在如图4(a)所示的正方形单胞中，基体和增强相都呈对称分布^[18]。图4(b)所示的有限元模型中圆柱形增强相的半径为R，正方体基体的边长为L。设钨丝为横观各向同性材料，弹性常数为E_W=410 GPa，泊松比0.28，块体非晶合金基体为各向同性材料，弹性常数为E_BMG=96 GPa，泊松比0.36。单元网格划分如图4(b)所示，考察增强相的体积分数(V_f)和微结构几何尺寸比值R/L对宏观力学性能的影响。

图2  不同条件下周期性单胞的变形方式

Fig. 2  Different deformation modes for periodic unit cell under different conditions

图3  对应图2(b)、(c)和(d)单胞施加的边界条件

Fig. 3  Boundary conditions corresponding to unit cells in Figs. 2(b), (c) and (d)

将不同体积分数下钨丝增强锆基块体非晶复合材料弹性模量的预测值(本文中建立的理论模型计算值，用Homo FEM表示)与Voigt模型和Reuss模型^[21]的计算结果(分别用Voigt和Reuss表示)以及文献[22-26]报道的实验值(Experimental results)进行比较，如表1所列。Voigt模型和Reuss模型公式^[21]如下：

Voigt:               (37)

Reuss:              (38)

图4  钨丝增强锆基块体非晶复合材料的周期性细观结构及其三维单胞网格模型

Fig. 4  Periodic microstructure (a) and three-dimensional grid cell model (b) of tungsten fiber reinforced Zr-based bulk metallic glass matrix composites

式中：E^H和V_f分别表示复合材料的等效弹性模量和钨丝体积分数。将计算结果作无量纲化处理即E^H/E_W、μ^H/μ_BMG和/E₀，其中。钨丝增强锆基块体非晶复合材料的等效弹性模量与钨丝体积分数关系如图5所示。在代表性体积单元内，单个纤维增强情况下不同R/L(不同体积分数)对整体等效弹性常数的影响如图6和表2所示。一定体积分数下不同R/L对整体等效弹性常数的影响如图7所示。

表1  不同体积分数钨丝增强锆基块体非晶复合材料弹性模量的预测值(Homo FEM)与文献报道的实验值及Voigt模型和Reuss模型的计算结果的比较

Table 1  Comparison of predicted elastic moduli (Homo FEM) of tungsten fiber reinforced Zr-based bulk metallic glass matrix composites with different volume fractions of tungsten fiber with their experimental values from literature and calculated values based on Voigt and Reuss models

从表1和图5可以看出，钨丝增强锆基块体非晶复合材料的弹性模量随着钨丝的体积分数的增加而增大，用均匀化方法预测的结果处在Voigt模型和Reuss模型的结果之间, 且与Voigt模型近乎重合。Voigt模型采用均匀应变场假设，得到的等效弹性模量为真实等效弹性模量的上限，Reuss模型采用均匀应力场假设，得到的等效弹性模量为真实等效弹性模量的下限。

表2  不同R/L下单相钨丝增强锆基块体非晶复合材料等效弹性常数的预测值

Table 2  Predicted equivalent elastic constants of tungsten fiber reinforced Zr-based bulk metallic glass matrix composites with different R/L ratio values

图5  复合材料等效弹性模量与钨丝体积分数关系

Fig. 5  Relationship between equivalent elastic modulus of composite and volume fraction of tungsten fiber

图6  单相钨丝增强下复合材料等效弹性参数与微结构尺寸比R/L的关系

Fig. 6  Relationships among equivalent elastic parameters and R/L ratio value under condition of single-phase tungsten enhanced

图7  多相钨丝增强且定钨丝体积分数(60%)时等效弹性参数与微结构尺寸比R/L的关系

Fig. 7  Relationship between equivalent elastic parameters of composites and R/L ratio value under condition of multi-phase tungsten enhanced and given volume fraction of tungsten fiber of 60%

说明使用均匀化方法预测非晶复合材料的等效弹性模量的合理性和正确性。预测结果与文献报道的实验数据基本吻合，但还存在一定差异，其原因可能是：钨丝增强锆基块体非晶复合材料并非完全线弹性材料，如基体锆基块体非晶合金像聚合物材料一样，具有黏弹性特征；有限元模型与实体模型存在差异；计算程序中解方程的方法存在误差。为减小与实验数据的差异，在进一步的计算中，尝试将基体的黏弹性能对整体弹性常数的影响考虑进去；选用高配置计算机，进行有限元建模，进一步细化单元网格；选用合适的解方程方法。

从图6和表2可以看出，在代表性体积单元内单相钨丝增强锆基块体非晶复合材料的等效弹性常数(其中和分别等效于拉梅常数λ和G)随着微结构几何尺寸比值R/L增加而增大；而等效泊松比随着R/L的增加而减小。图7显示，代表性体积单元内钨丝体积分数一定(60%)而增强相数量不同时复合材料等效弹性常数随着微结构几何尺寸比值R/L的增加呈略微减小的趋势，但趋势不明显，近似于有限元单胞中含单相钨丝的计算值；等效泊松比随着R/L的增加略微增大，同样近似于有限元单胞中含单相钨丝的计算结果。因此，在给定钨丝体积分数时，单胞中钨丝相数量对整体弹性常数影响很小，可忽略不计。在实际工程计算中，只考虑有限元单胞中含单相钨丝情况下，不同微结构几何尺寸比值R/L对复合材料整体等效弹性常数的影响即可。

4  结论

1) 利用均匀化方法的三维有限元格式对钨丝增强锆基块体非晶复合材料的等效弹性常数进行预测，预测结果与实验数据相吻合，且均匀化方法与其他方法相近，是一种具有严格数学依据、高效而实用的方法，在预测钨丝增强锆基块体非晶复合材料等效弹性常数方面具有相对较强的优势。

2) 通过计算，研究了增强相的体积分数对钨丝增强锆基块体非晶复合材料宏观力学性能的影响，得到了合理的数值结果，并探讨了单相增强和多相增强情况下增强相的尺寸大小对复合材料的宏观力学性能的影响，得到了有意义的结论，为进一步确定钨丝增强锆基块体非晶合金的成分设计和微观应力，奠定了有力的理论基础。

REFERENCES

[1] 张晓立, 王金相, 孙宇新, 付艳恕, 刘家骢. 钨丝增强块体非晶复合材料的研究进展[J]. 稀有金属材料与工程, 2008, 37(8): 1323-1328.

ZHANG Xiao-li, WANG Jin-xiang, SUN Yu-xin, FU Yan-shu, LIU Jia-cong. Study progress of tungsten-fiber reinforced bulk metallic-glass matrix composites[J]. Rare Metal Materials and Engineering, 2008, 37(8): 1323-1328.

[2] 柳 林, 孙 民, 谌 祺, 刘 兵, 邱春雷. Zr-Cu-Ni-Al-Nb大块非晶合金的晶化行为、力学性能及电化学腐蚀行为的研究[J]. 物理学报, 2006, 55(4): 1930-1935.

LIU Lin, SUN Min, CHEN Qi, LIU Bing, QIU Chun-lei. Crystallization, mechanical and corrosion properties of Zr-Cu-Ni-Al-Nb bulk glassy alloys[J]. Acta Physica Sinica, 2006, 55(4): 1930-1935.

[3] PEKER A, JOHNSON W L. A highly processable metallic glass Zr_41.2Ti_13.8Cu_12.5Ni_10.0Be_22.5[J]. Applied Physics Letters, 1993, 63(17): 2342-2344.

[4] 许 福, 龙志林, 彭 建, 张 平. 块体非晶合金剪切带的原子力纳米压痕行为[J]. 中国有色金属学报, 2011, 21(6): 1444-1449.

XU Fu, LONG Zhi-lin, PENG Jian, ZHANG Ping. Atomic force microscope nanoindentation behavior of shear bands of bulk metallic glasses[J]. The Chinese Journal of Nonferrous Metals, 2011, 21(6): 1444-1449.

[5] 杨元政, 董振江, 仇在宏, 陈小祝, 谢致嶶, 白晓军. 块体非晶合金Cu₅₈Zr₂₀Ti₂₀Mo₂的形成与力学性能[J]. 中国有色金属学报, 2007, 17(7): 1090-1095.

YANG Yuan-zheng, DONG Zhen-Jiang, QIU Zai-hong, CHEN Xiao-zhu, XIE Zhi-wei, BAI Xiao-jun. Formation and mechanical properties of Cu₅₈Zr₂₀Ti₂₀Mo₂ bulk metallic glass[J]. The Chinese Journal of Nonferrous Metals, 2007, 17(7): 1090-1095.

[6] BRUCK H A, CHRISTMAN T, ROSAKIS A J, JOHNSON W L. Quasi-static constitutive behavior of Zr_41.25Ti_13.75 Ni₁₀ Cu_12.5Be_22.5 bulk amorphous alloys[J]. ScriptaMaterialia, 1994, 30(4): 429-434.

[7] LONG Z L, SHAO Y, XIE G Q, ZHANG P, SHEN B L, INOUE A. Enhanced soft-magnetic and corrosion properties of Fe-based bulk glassy alloys with improved plasticity through the addition of Cr[J]. Journal of Alloys and Compounds, 2008, 462: 52-59.

[8] 马卫锋, 寇宏超, 李金山, 陈春生, 杜三明, 周 廉, 付恒志. 钨丝增强Zr基非晶复合材料动态力行为及断裂特性[J]. 中国有色金属学报, 2008, 18(6): 1045-1050.

MA Wei-feng, KOU Hong-chao, LI Jin-shan, CHEN Chun-sheng, DU San-ming, ZHOU Lian, FU Heng-zhi. Dynamic mechanical behaviors and fracture characteristic of tungsten fiber reinforced Zr-based metallic glass matrix composites[J]. The Chinese Journal of Nonferrous Metals, 2008, 18(6): 1045-1050.

[9] HILL R. A self-consistent mechanics of composite materials[J]. Journal of the Mechanics and Physics of Solids, 1965, 13(4): 213-222.

[10] BUDIANSKY B. On the elastic moduli of some heterogeneous materials[J]. Journal of the Mechanics and Physics of Solids, 1965, 13(4): 223-227.

[11] MORI T, TANAKA K. Average stress in matrix and average elastic energy of material with misfitting inclusions[J]. Acta Metallurgica, 1973, 21(5): 571-574.

[12] 杨庆生, 陈浩然. 夹杂问题中的自洽有限元法和复合材料的平均弹性性能[J]. 复合材料学报, 1992, 9(1): 79-84.

YANG Qing-sheng, CHEN Hao-ran. Self-consistent finite element method for the problems of inclusion and the average elastic properties of composite materials[J]. Acta Materiae Compositae Sinica, 1992, 9(1): 79-84.

[13] HASSANI B, HINTON E. A review of homogenization and topology optimization I-homogenization theory for media with periodic structure[J]. Computer & Structures, 1998, 69: 707-717.

[14] HASSANI B, HINTON E. A review of homogenization and topology optimization II-analytical and numerical solution of homogenization equations[J]. Computers & Structures, 1998, 69: 719-738.

[15] 董纪伟, 孙良新, 洪 平. 基于均匀化方法的三维编织复合材料等效弹性性能预测[J]. 宇航学报, 2005, 26(4): 482-486.

DONG Ji-wei, SUN Liang-xin, HONG Ping. Homogenization- based method for predicting effective elastic properties of three-dimensional braided composites[J]. Journal of Astronautics, 2005, 26(4): 482-486.

[16] 付海雄, 马石城, 李晓全. 用渐近均匀化方法计算水泥土的等效弹性常数[J]. 湘潭大学自然科学学报, 2006, 28(1): 117-120.

FU Hai-xiong, MA Shi-cheng, LI Xiao-quan. Application of asymptotic homogenization theory to calculate the equivalent elastic constants of cement-soil[J]. Natural Science Journal of Xiangtan University, 2006, 28(1): 117-120.

[17] 闫晓鹏, 牛卫晶, 王志华, 马宏伟. 基于均匀化理论的混凝土等效弹性模量的数值模拟[J]. 太原理工大学学报, 2011, 42(2): 212-214.

YAN Xiao-peng, NIU Wei-jin, WANG Zhi-hua, MA Hong-wei. Numerical simulation of equivalent modulus of concrete based on homogenization theory[J]. Journal of Taiyuan University of Technology, 2011, 42(2): 212-214.

[18] DRAGOI D, E, CLAUSEN B, BOURKE M A M. Investigation of thermal residual stresses in tungsten-fiber/ bulk metallic glass matrix composites[J]. Scripta Materialia, 2001, 45: 245-252.

[19] 王勖成. 有限单元法[M]. 北京: 清华大学出版社, 2003: 46-47.

WANG Xu-cheng. Finite element method[M]. Beijing: Tsinghua University Press, 2003: 46-47.

[20] GUEDES J, KIKUCHI N. Preprocessing and post-processing for materials based on the homogenization method with adaptive finite element methods[J]. Computer Methods in Applied Mechanics and Engineering, 1990, 83: 143-198.

[21] SHEN T K, HING P. Ultrasonic through-transmission method of evaluation the modulus of elasticity of Al₂O₃-ZrO₂ composite[J]. Journal of Materials Science, 1997, 32: 6633-6638.

[22] CONNER R D, DANDLIKER R B, JOHNSON W L, Mechanical properties of tungsten and steel fiber reinforced Zr_41.25Ti_13.75Cu_12.5Ni₁₀Be_22.5 metallic glass matrix composites[J]. Acta Materialia, 1998, 46(17): 6089-6102.

[23] ZHANG H, ZHANG Z F, WANG Z G, QIU K Q, ZHANG H F, ZANG Q S, HU Z Q. Fatigue damage and fracture behavior of tungsten fiber reinforced Zr-based metallic glassy composite[J]. Materials Science and Engineering A, 2006, 418: 146-154.

[24] ZHANG H, ZHANG Z F, WANG Z G, ZHANG H F. Deformation and damage evolution of tungsten fiber reinforced metallic glass matrix composite induced by compression[J]. Materials Science and Engineering A, 2008, 483/484: 164-167.

[25] ZHANG B Y, CHEN X H, WANG S S, LIN D Y, HUI X D. High strength tungsten wire reinforced Zr-based bulk metallic glass matrix composites prepared by continuous infiltration process[J]. Materials Letters, 2013, 93: 210-214.

[26] WANG G, CHEN D M, SHEN J, STACHURSKI Z H, QIN Q H, SUN J F, ZHOU B D. Deformation behaviors of a tungsten-wire/bulk metallic glass matrix composite in a wide strain rate range[J]. Journal of Non-Crystalline Solids, 2006, 352: 3872-3878.

(编辑  陈卫萍)

基金项目：国家自然科学基金资助项目(51071134)；湖南省自然科学基金重点项目(12JJ2024)

收稿日期：2013-10-08；修订日期：2014-01-23

通信作者：龙志林，教授，博士；电话：0731-58298287；E-mail: longzl@xtu.edu.cn