Sound Velocity and Shock Response Behavior of Cu/PMMA Composites
-
摘要: 基于熔融共混法,制备了一系列不同配比且随机分散的Cu/PMMA复合材料,重点研究了Cu颗粒含量对PMMA基体声速与冲击压缩行为的影响。超声测试结果表明,随着Cu颗粒含量的增加,声波的衰减使材料的横、纵波声速皆呈缓慢下降趋势,由此使体积声速亦呈减小趋势。基于平板撞击实验,获得了冲击压力在1.1~6.0 GPa范围内各复合材料的冲击波速度-粒子速度方程。Cu/PMMA复合材料声阻抗的升高使Hugoniot参数λ逐渐增大,而零压体积声速减小,与常压体积声速所表现出的变化趋势一致。结合已有的压力-粒子速度关系模型,对各材料的压力-粒子速度曲线进行了讨论。在此基础上,归纳出一种基于上述模型的用于预测金属粒子填充聚合物基复合材料压力-密度关系的可靠方法。
-
关键词:
- 冲击压缩 /
- 声速 /
- 冲击波速度-粒子速度方程 /
- 压力-密度关系
Abstract: In our work, a series of Cu/PMMA composites with different components were prepared using the melting blending method, in which particles are randomly dispersed in PMMA matrix without agglomeration. Then again study was conducted on the Cu particles content’s influence on the sound velocity and impact compression behavior of PMMA matrix. The ultrasonic test results show that with Cu particles content increasing, the sound wave attenuation makes a slow decreasing tendency of transversal and longitudinal sound velocities in material, which in turn decreases its bulk sound velocity. Based on the plate impact test, the shock wave velocity-particle velocity (D-u) equations of Cu/PMMA composites in the impact pressure range of 1.1–6.0 GPa were obtained. Owing to the increase of the acoustic impedance of Cu/PMMA composites, Hugoniot parameter shows an increase while the fitted zero-pressure sound velocity tends to decrease, which turns to be consistent with the variation of bulk sound velocity at atmospheric pressure. In addition, pressure-particle velocity (p-u) curves of the composites were discussed on the basis of the p-u model. And a reliable method was proposed to predict pressure-density (p-ρ ) relationship of polymer matrix composites filled with metal particles. -
球体的冲击破碎是颗粒物质力学的基本问题。球体结构材料广泛存在于日常生活和工程应用中,例如煤粒破碎[1-2]、落石灾害[3]、化工产品[4]等。研究球体材料在冲击载荷下的破碎,对于提高颗粒物质的加工技术、节约能源、促进经济发展、有效监控及预防自然灾害等具有重要价值。
早期国内外研究人员对单颗粒的冲击破碎开展了大量的实验研究,研究的重点是颗粒的临界破碎速度、颗粒的破碎模式、碎片形貌以及碎块的统计特征等。Andrews等[5-6]对钠钙玻璃球和陶瓷球进行了研究,在临界速度下,球体从刚性靶板回弹;超过临界速度后,钠钙玻璃球碎裂为大量尺寸不均的细小碎片,而陶瓷球则碎裂为若干大小均匀形似“月牙状”的碎块,并将钠钙玻璃球体破碎归因于球体压缩端的剪切应力,而将陶瓷球体碎裂归因于球体中部的冲击拉伸应力。Salman等[7-9]对不同类型材料和大小的脆性球体进行试验,结果表明大部分脆性材料的破坏模式是相似的,即:在中低速下,球体的局部裂纹扩展导致局部破碎,在撞击端形成压缩平台;在高速下,球体碎裂为大量难以识别的粉状碎片。Chau[10]、Wu[11]等根据冲击动能的转换分析了脆性球体破碎后碎块的统计特征。
在数值计算方面,易洪昇等[12]使用ABAQUS软件模拟了冲击过程中的球体内部应力情况,由于颗粒破碎涉及大变形,采用有限元方法会出现网格畸变问题,因此不适合模拟大变形和损伤现象。离散元方法(Discrete element method, DEM)能够满足岩体工程中破裂和裂纹发展等问题的研究,适合于从本质上研究固体介质的力学特性[13]。Shen等[14]利用离散元软件分析了球形岩石在不同应变率情况下碎块特征、碎块数量、损伤率等的变化情况。Carmona等[15]采用离散元方法研究了球体碎裂过程中裂纹的扩展过程,发现碎裂的大块碎片质量符合两参数的Weibull分布。
以石英玻璃球作为典型的脆性球体,开展了较广泛冲击速度下石英玻璃球撞击刚性壁的破坏现象研究,利用离散元数值软件PFC3D分析了石英玻璃球在不同撞击速度下的撞击碎裂过程。
1. 实验装置及试件
实验采用高速气枪作为石英玻璃球的加载发射装置,撞击区域使用碳化钨材料作为近似刚性靶板,碳化钨的弹性模量为720 GPa,密度为15.6 g/cm3,泊松比为0.21。碳化钨周围附有厚壁钢块,如图1所示。采用Kirana超高速摄像机记录球体的破碎过程,拍摄幅频为2×105帧每秒,高速相机采用断路方式触发,触发线设置在气枪的出口位置。在球体撞击区域的外围设计了一个有机玻璃防护装置用于收集碎片。
高速气枪出口距离靶板250 mm,通过圆柱形长杆调节高速气枪和靶板的垂直度,从而实现玻璃球的正撞破坏。通过水平放置和位于球体上方放置的高速相机,实验观察石英玻璃球的飞行姿态,结果表明玻璃球体未发生明显偏转,撞击靶板位置为初始预设位置。根据超高速相机拍摄可获得玻璃球体的飞行时间,从而确定玻璃球体飞行的平均速度,近似等效为球体撞击速度。
实验试件采用昊天玻璃制品有限公司生产的超纯石英玻璃球,其SiO2的质量分数大于99.99%。球体直径为(10.00±0.02)mm,质量为(1.150±0.015)g,具有良好的透光性,且表面无明显缺陷。
2. 碎裂过程及结果分析
为了研究玻璃球的临界冲击破碎速度,在低速下开展了大量的实验。实验结果表明:当石英玻璃球以11.87 m/s撞击刚性靶板时,碰撞为近似弹性碰撞,球体以略低于原速(回弹速度10.25 m/s)从刚性壁回弹,除却撞击点应力集中破坏外,球体内部没有产生宏观可见裂纹。石英玻璃球撞击前和回弹后如图2所示。
图3是冲击速度为78 m/s的冲击破碎过程。从图3可以看出,球体的破碎可以分为3个阶段。在初始阶段,撞击端首先发生破坏,形成经过子午面的宏观裂纹,同时形成细条状碎片,这些碎片以撞击点为中心向四周扩散,而在远离撞击端部分还保持相对完整;随着撞击的持续进行,在撞击端的碎片继续向四周运动,在远离撞击端的部分产生“龟甲状”的裂纹,并最终碎裂为若干“圆锥状”的碎块。
当冲击速度为35 m/s时,球体从刚性壁回弹,在撞击端产生了一些细小的碎片,如图4(a)所示。图5(a)表明,实验回收后的碎片存在3个破坏区域,分别为压缩破碎区、表面剥落区以及剪切破坏区。其中压缩破碎区是一个圆锥状的区域,这块区域的石英玻璃已不具备透光性,说明这一区域的石英玻璃发生了剧烈的压缩破碎,但是这些碎块还有部分镶嵌在球体上并没有从球体上脱落。在压缩破碎区的周围是一个环形的表面剥落区,受到表面波作用一些片状的碎片从球体上剥落,产生层片状的碎片。在球体内部,石英玻璃的折光性发生了改变,说明在剪切力的作用下石英玻璃内部产生剪切裂纹带,剪切裂纹带由压缩破碎区开始起裂,沿着子午面由撞击端向试件内部扩展。
当冲击速度为50 m/s时,从球体的透光性可以看出,球体完全碎裂为若干碎块,如图4(b)所示。在撞击端同样会出现类似于35 m/s时的压缩破碎现象,并且涉及的范围更大。由于冲击动能的增加,剪切破碎区完全贯穿了球体的子午面,最终导致石英玻璃球碎裂成若干类似于“月牙状”的碎块,如图5(b)所示。
当冲击速度为78 m/s时,与上述几种破碎情况完全不同,石英玻璃球不再从刚性壁回弹,而是发生了“坍塌式”的破碎,如图4(c)所示。在撞击端,产生大量细小的碎块,在球体的另一端表面有许多类似于“经线形”和“纬线形”的宏观裂纹,最终产生大量类似于“龟甲状”的碎片,内部也产生大量细小碎片,如图5(c)所示。
当冲击速度为135 m/s时,发生“蘑菇云”似的破碎现象,如图4(d)所示。产生的碎片更加细小,在远离撞击端的地方发生明显的“层裂”现象。球体在高速撞击过程中碎裂为极细小的碎片,透明球体不易观察裂纹发展,将石英玻璃球表面涂黑,有利于观测球体表面裂纹扩展。由图5(d)可知,随着速度的增加,球体内部破碎更为剧烈,产生大量粉末状细小碎片。
综上所述,动态冲击下石英玻璃球的破坏形式主要有4种特征:(1)局部破碎,在球体上形成3个破碎区域;(2)中心劈裂,沿着冲击轴线碎裂为类似于“月牙状”的碎片;(3)坍塌破坏,石英玻璃球整体破碎,碎裂为若干碎块;(4)层裂破坏,在远离撞击端的地方产生“蘑菇云”似的“层裂”破坏。
3. 离散元数值计算
3.1 离散元模型的建立
选择目前离散元中使用广泛的平行黏结模型(Linear parallel bond mode)。平行黏结模型本质上是用一层线性的黏结体将相互重叠的颗粒绑定起来,这些黏结的颗粒能很好地模拟岩石、玻璃等脆性材料的压缩破坏过程[16]。
为了使模型具有较高的密实度,将模型中的孔隙率设定为0.2,最大颗粒与最小颗粒的粒径比为3∶2。颗粒间的摩擦系数对材料的宏观参数影响不大[17],故取经验值0.577。
根据文献[18-19]的标定方法,通过单轴压缩标定了模型的弹性模量、泊松比和抗压强度,通过单轴拉伸标定抗拉强度,由三点弯标定弯曲强度,带缺口的三点弯标定断裂韧性和颗粒半径。最终得到的主要微观参数如表1所示。标定的宏观参数如表2所示,得到石英玻璃的宏观参数与制造商提供的材料参数基本一致,因此可以用于模拟石英玻璃的冲击破碎。
表 1 石英玻璃离散元模型的主要微观参数Table 1. Main microscopic parameters of discrete element mode of quartz glassEffective modulus of linear contact/GPa Normal to shear stiffness ratio of linear contact Porosity Minimum radius of particles/mm Size ratio of maximum and minimum particles Tensile strength of contact/MPa Shear strength of contact/MPa 55 2.9 0.2 0.1 1.5 300 600 表 2 石英玻璃在常态下的物理参数Table 2. Physical properties of quartz glass under ambient conditionMethod Equivalent density/(kg·m–3) Elastic modulus/GPa Poisson’s ratio Compressive strength/MPa Tensile strength/MPa Bending strength/MPa Fracture toughness/(N·m–3/2) Manufacture provide 2.203 77.8 0.170 860 50 67.0 0.78 DEM numerical simulation 2.203 78.0 0.172 798 50 67.4 0.85 根据上述标定的离散元微观参数建立与实验同尺寸的几何模型。采用刚性墙面模拟实验中的刚性壁,在模拟初始时刻,给所有颗粒施加垂直撞击刚性壁的速度,使其匀速撞击刚性壁。
3.2 石英玻璃球冲击刚性壁破碎过程
图6是球体以78 m/s撞击刚性壁的过程。模拟中不同颜色代表的是碰撞过程中产生的不同碎块。对比图3和图6可知,模拟与实验吻合较好,模拟中玻璃球的破坏过程存在明显的时序性。在初始时刻,在撞击端发生局部的破碎,产生细小的碎片,同时裂纹向远离撞击端的地方传播;随着撞击的进行,撞击端产生的碎片向外运动,远离撞击端的地方碎裂为较大的锥形碎块,当撞击端处的细小碎片运动到更远处时,这些较大的锥形碎块以一定的速度继续撞击刚性壁。
球体以78 m/s撞击刚性壁的过程中,球体速度(左端粒子平均速度)、球体内部裂纹数量(粒子之间的力键)以及球体撞击力的时程曲线如图7所示。石英玻璃球撞击破碎过程大致可分为3个阶段:第1阶段,可近似为弹性压缩过程,球体所受到的冲击载荷也近似线性增加,大约经历
3μs 后撞击力达到最大值,在弹性波还未传播至左端时,球体速度保持不变,之后速度下降,在球体撞击端裂纹发展并产生少量碎片,球体整体裂纹发展较缓慢;第2阶段,为球体破碎阶段,撞击力在碎裂过程中发生明显振荡,同时球体速度也存在显著振荡,这是由于该阶段大量密集裂纹的生成,使弹性卸载波在球体中传播并来回反射,最终球体撞击端(右端粒子)几乎完全破坏,球体裂纹数整体稳定,撞击力也卸载为零附近;第3阶段,为球体二次撞击阶段,远离撞击端(左端粒子)在碎裂过程中还保留残余速度,在球体整体破碎后仍以较低的速度再次撞击刚性壁,但对石英玻璃球整体的撞击力和裂纹增长影响甚微。图8为在不同速度下石英玻璃球速度随时间变化的曲线。图8(a)显示了撞击速度低于临界破碎速度时球体速度的变化过程。可以看出,球体撞击靶板后匀速下降,最终几乎以原速从靶板回弹。图8(b)显示了撞击速度高于临界破碎速度时球体速度的变化过程。可以看出,球体左端粒子在撞击刚性壁后平均速度急剧下降并保留部分残余速度,撞击速度越大,速度下降越快,残余速度也越大。球体左端粒子以残余速度自由飞行,石英玻璃球撞击端(右端粒子)完全破碎,左端粒子在向右运动过程中没有任何阻力,将球体匀速运动部分对时间积分,发现球体的运动距离均约为球体直径的1/2。
图9给出了0~300 m/s撞击速度下石英玻璃球撞击力随撞击速度的变化规律。从图9可以看出,撞击速度为20 m/s左右时为撞击力增长规律的临界转折点,此前球体未产生明显破碎,此后球体破碎使得撞击力上涨趋势大幅衰减。
在临界转折速度前,石英玻璃球撞击刚性靶板过程可近似认为是弹性碰撞,通过Hertz接触理论可得弹性碰撞接触力[20]
Pm=(53πρ)35(34KI)−25V65R2 (1) 式中:
KI=1−μ12E1+1−μ22E2 ,其中E1、μ1 和E2、μ2 分别为球体材料和刚性靶板的弹性模量和泊松比;ρ 为球体材料的密度;R为球体的半径。图10给出了临界转折速度附近石英玻璃球撞击力随撞击速度的变化规律。从图10可以看出:当速度低于临界转折速度时,球体不破碎,内部几乎不产生裂纹,球体整体近似弹性碰撞,Hertz接触理论可以较好地预测石英玻璃球的撞击力;而速度高于临界转折速度时,球体发生破碎,不再具有球体整体性,破碎过程将卸载撞击力,并消耗部分能量,飞散的碎片也将带走大量动能;随着冲击速度的增加,石英玻璃球的极限撞击将等效于流体撞击。
4. 结 论
冲击速度在临界破坏速度下为弹性碰撞,球体以略低于原速从刚性壁回弹,除却撞击点应力集中破坏外,球体内部没有产生宏观可见裂纹。低速碰撞时,球体呈现“压缩破碎区-表面剥落区-剪切破坏区”的破坏结构,剪切裂纹沿着子午面由撞击端向试件内部扩展,并还不足以形成碎片;中速碰撞时,由于剪切破坏区充分扩展,致使球体碎裂为若干“月牙状”的碎块;高速碰撞时,石英玻璃球发生坍塌式破碎,表面产生大量类似于“龟甲状”的碎片,内部也产生大量细小碎片;当速度更高时,在远离撞击端有明显的层裂现象产生,球体破碎产生大量粉末状碎片。
离散元软件PFC3D再现了石英玻璃球撞击刚性壁的过程,球体在高速碰撞下破碎可以分为弹性压缩、整体破碎和二次撞击3个阶段。石英玻璃球撞击力随撞击速度呈两段式增长,在球体碎裂前,撞击过程近似为弹性碰撞过程,Hertz接触理论可以较好地描述其撞击力,在球体碎裂后,破碎过程将卸载撞击力,并消耗部分能量,飞散的碎片也将带走大量动能,撞击力远低于Hertz接触理论值,并且随着撞击速度增大,偏差越大。
-
表 1 Cu/PMMA复合材料的弹性力学参数
Table 1. Elastic mechanical parameters of Cu/PMMA composites
wCu/% K/GPa G/GPa E/GPa 10 5.66 2.33 6.17 25 5.80 2.43 6.40 40 6.24 2.71 7.09 60 7.03 3.19 8.32 表 2 Cu/PMMA复合材料的击靶参数
Table 2. Target parameters of Cu/PMMA composites
wCu/% W/(m·s−1) d/mm D/(km·s−1) u/(km·s−1) 10 328 1.814 2.988 0.296 462 1.802 3.046 0.416 621 1.814 3.161 0.558 902 1.803 3.725 0.797 40 313 1.754 2.662 0.276 462 1.742 2.702 0.407 625 1.751 2.891 0.546 901 1.742 3.380 0.773 60 307 1.913 2.437 0.264 446 1.893 2.553 0.381 640 1.891 2.836 0.539 897 1.918 3.252 0.741 表 3 Cu/PMMA复合材料的Li模型参数
Table 3. Impact parameters corresponding to Li model for Cu/PMMA composites
wCu/% W/(m·s−1) p1/GPa p2/GPa p/GPa 10 328 1.05 7.89 1.14 462 1.51 11.31 1.64 621 2.10 15.74 2.28 902 3.53 26.49 3.81 40 313 0.87 6.56 1.33 462 1.31 9.81 1.99 625 1.88 14.09 2.86 901 3.11 23.32 4.73 60 307 0.77 5.74 1.59 446 1.16 8.68 2.40 640 1.82 13.64 3.77 897 2.87 21.50 5.94 -
[1] MEYERS M A. Dynamic behavior of materials [M]. New York: Wiley, 1994: 1−5. [2] ZHANG Y N, ZHENG L X, SUN G Z, et al. Failure mechanisms of carbon nanotube fibers under different strain rates [J]. Carbon, 2012, 50(8): 2887–2893. doi: 10.1016/j.carbon.2012.02.057 [3] BIE B X, HAN J H, LU L, et al. Dynamic fracture of carbon nanotube/epoxy composites under high strain-rate loading [J]. Composites Part A: Applied Science and Manufacturing, 2015, 68: 282–288. doi: 10.1016/j.compositesa.2014.10.001 [4] YASHIRO S, OGI K, NAKAMURA T, et al. Characterization of high-velocity impact damage in CFRP laminates: Part I–experiment [J]. Composites Part A: Applied Science and Manufacturing, 2013, 48: 93–100. doi: 10.1016/j.compositesa.2012.12.015 [5] YASHIRO S, OGI K, YOSHIMURA A, et al. Characterization of high-velocity impact damage in CFRP laminates: Part Ⅱ - prediction by smoothed particle hydrodynamics [J]. Composites Part A: Applied Science and Manufacturing, 2014, 56: 308–318. doi: 10.1016/j.compositesa.2013.04.012 [6] XIE W B, ZHANG W, KUANG N H, et al. Experimental investigation of normal and oblique impacts on CFRPs by high velocity steel sphere [J]. Composites Part B: Engineering, 2016, 99: 483–493. doi: 10.1016/j.compositesb.2016.06.020 [7] GAY E, BERTHE L, BOUSTIE M, et al. Study of the response of CFRP composite laminates to a laser-induced shock [J]. Composites Part B: Engineering, 2014, 64: 108–115. doi: 10.1016/j.compositesb.2014.04.004 [8] GIANNAROS E, KOTZAKOLIOS A, KOSTOPOULOS V, et al. Hypervelocity impact response of CFRP laminates using smoothed particle hydrodynamics method: implementation and validation [J]. International Journal of Impact Engineering, 2019, 123: 56–69. doi: 10.1016/j.ijimpeng.2018.09.016 [9] CHEN X, LI Y L, ZHI Z, et al. The compressive and tensile behavior of a 0/90 C fiber woven composite at high strain rates [J]. Carbon, 2013, 61: 97–104. doi: 10.1016/j.carbon.2013.04.073 [10] LONG X J, LI B, WANG L, et al. Shock response of Cu/graphene nanolayered composites [J]. Carbon, 2016, 103: 457–463. doi: 10.1016/j.carbon.2016.03.039 [11] XIE W B, ZHANG W, GUO L C, et al. The shock and spallation behavior of a carbon fiber reinforced polymer composite [J]. Composites Part B: Engineering, 2018, 153: 176–183. doi: 10.1016/j.compositesb.2018.07.047 [12] DANDEKAR D P, HALL C A, CHHABILDAS L C, et al. Shock response of a glass-fiber-reinforced polymer composite [J]. Composite Structures, 2003, 61(1/2): 51–59. doi: 10.1016/S0263-8223(03)00031-X [13] JIAN W R, LONG X J, TANG M X, et al. Deformation and spallation of shock-loaded graphene: effects of orientation and grain boundary [J]. Carbon, 2018, 132: 520–528. doi: 10.1016/j.carbon.2018.02.070 [14] MENG Z X, HAN J L, QIN X, et al. Spalling-like failure by cylindrical projectiles deteriorates the ballistic performance of multi-layer graphene plates [J]. Carbon, 2018, 126: 611–619. doi: 10.1016/j.carbon.2017.10.068 [15] TIAN Y, ZHANG H, ZHAO J, et al. High strain rate compression of epoxy based nanocomposites [J]. Composites Part A: Applied Science and Manufacturing, 2016, 90: 62–70. doi: 10.1016/j.compositesa.2016.06.008 [16] REN S Y, ZHANG Q M, WU Q, et al. A debris cloud model for hypervelocity impact of the spherical projectile on reactive material bumper composed of polytetrafluoroethylene and aluminum [J]. International Journal of Impact Engineering, 2019, 130: 124–137. doi: 10.1016/j.ijimpeng.2019.04.011 [17] RAULS M B, RAVICHANDRAN G. Structure of shock waves in particulate composites [J]. Journal of Applied Physics, 2020, 127(6): 065902. [18] LI J B, LI W B, WANG X M, et al. Shock response and prediction model of equation of state for aluminum powder/rubber matrix composites [J]. Materials and Design, 2020, 191: 108632. doi: 10.1016/j.matdes.2020.108632 [19] BEK Y K, HAMDIA K M, RABCZUK T, et al. Micromechanical model for polymeric nano-composites material based on SBFEM [J]. Composite Structures, 2018, 194: 516–526. doi: 10.1016/j.compstruct.2018.03.064 [20] REN H L, LI W, NING J G, et al. Effect of temperature on the impact ignition behavior of the aluminum/polytetrafluoroethylene reactive material under multiple pulse loading [J]. Materials and Design, 2020, 189: 108522. [21] REN S Y, ZHANG Q M, WU Q, et al. Influence of impact-induced reaction characteristics of reactive composites on hypervelocity impact resistance [J]. Materials and Design, 2020, 192: 108722. [22] 胡昌明, 李雪梅, 彭建祥, 等. 冲击载荷下K9玻璃的光学特性 [J]. 高压物理学报, 2017, 31(5): 573–578. doi: 10.11858/gywlxb.2017.05.010HU C M, LI X M, PENG J X, et al. Optical properties of K9 glass under shock loading [J]. Chinese Journal of High Pressure Physics, 2017, 31(5): 573–578. doi: 10.11858/gywlxb.2017.05.010 [23] 谭华. 实验冲击波物理[M]. 北京: 国防工业出版社, 2018: 281−304.TAN H. Experimental shock wave physics [M]. Beijing: National Defense Industry Press, 2018: 281−304. [24] CARTER W J, MARSH S P. Hugoniot equation of state of polymers [R]. New Mexico: Los Alamos National Laboratory, 1995. [25] JORDAN J L, CASEM D, ZELLNER M. Shock response of polymethylmethacrylate [J]. Journal of Dynamic Behavior of Materials, 2016, 2(3): 372–378. doi: 10.1007/s40870-016-0071-5 [26] TSOU F K, CHOU P C. The control-volume approach to hugoniot of macroscopically homogeneous composites [J]. Journal of Composite Materials, 1970, 4(4): 526–537. doi: 10.1177/002199837000400408 [27] TORVIK P J. Shock propagation in a composite material [J]. Journal of Composite Materials, 1970, 4(3): 296–309. doi: 10.1177/002199837000400302 [28] 经福谦. 实验物态方程导引[M]. 2版. 北京: 科学出版社, 1999: 95−99.JING F Q. Introduction to experimental equation of state [M]. 2nd ed. Beijing: Science Press, 1999: 95−99. -