网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

水泥基材料高速纳米压痕聚类分析与矿相识别  PDF

  • 单钰涵 1,2
  • 田东诚 3
  • 赵建胜 4
  • 陈兵 1,2
  • 侯东伟 1,2
1. 上海交通大学 船舶海洋与建筑工程学院,上海 200240; 2. 上海交通大学 上海市公共建筑和础设施数字化运维重点实验室,上海 200240; 3. 中国石油集团大庆钻探工程公司 钻井生产技术服务一公司,黑龙江 大庆 163358; 4. 中石化中原石油工程有限公司 固井公司,河南 濮阳 457000

中图分类号: TU528.01

最近更新:2023-05-30

DOI:10.3969/j.issn.1007-9629.2023.05.015

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

通过高速纳米压痕技术测绘了3种水灰比(0.2、0.4和0.6)水泥净浆的弹性模量与硬度云图,引进K‑medoids聚类分析方法划分水泥基材料的微观矿相,根据K‑medoids聚类结果将力学性能图转化为矿相分布图,并对矿相的尺寸分布进行了统计分析.结果表明:水泥基材料主要包括多孔复合相(PP)、低密度/高密度C‑S‑H相(LD/HD C‑S‑H)、氢氧化钙相(CH)和超高性能相(UHP)等;随着水灰比的增加,PP、LD C‑S‑H的含量逐渐增加,HD C‑S‑H、UHP的含量逐渐减少;水泥石微观矿相的尺寸符合对数正态分布.

水泥基材料具有非均质性和多尺度的特

1‑2.在微观尺度上,硅相(C3S、C2S)、铝相(C3A)和铁相(C4AF)的水化产物发生聚集、胶凝,形成了微观结构.各矿相的力学性能和分布形态决定了水泥基材料的宏观性能,同时也是对其进行多尺度描述的基础参3‑4.

纳米压痕技

5‑7是近年来发展的一种先进的微观测量技术,可以在纳米和亚微米尺度测量材料的弹性模量、硬度和蠕变性能等.Oliver8‑9系统建立了基于纳米压痕曲线的弹性模量、硬度分析理论.Constantinides10采用点阵压痕方法,结合统计分析确定了水泥基多相材料微观各相的力学性能和含量.此外,采用纳米压痕对微观矿相的性能差异进行表征,还提供了一种分析材料矿相构成与分布的新方法.通过纳米压痕绘图技术,亦称为3Dmapping技11,可以得到复合材料中特定区域的微观力学性能及其位置分布之间的联系,已经在水12、岩13、合14‑18以及生物材19的微观组成分析中得到了应用.

基于纳米压痕的矿相分析通常采用解卷积方

10进行数据处理,但存在初值影响大、计算复杂程度高等问题.在人工智能领域,聚类分析是一种高效的数据统计方法和智能识别技术,已广泛应用于图像分割、信息检索和数据挖掘20.Kmedoids聚类分21方法通过对数据的聚合特征进行数学刻画,可以统计分析3Dmapping云图的矿相性能,同时得到矿相的几何分布图.

本文首先通过高速纳米压痕绘图技术,测绘了3种水灰比(质量比,文中涉及的水灰比、组成等除特别说明外均为质量比或质量分数)水泥净浆的弹性模量和硬度云图;然后,采用Kmedoids聚类分析识别水泥净浆中的主要微观矿相,获得矿相分布图;最后,通过灰度统计方法分析各相的尺寸分布.

1 试验

1.1 试验材料

由普通P·O 42.5硅酸盐水泥制备水泥净浆试块,水灰比(mW/mC)分别为0.2、0.4和0.6.在标准养护条件 ((20±2) ℃、RH≥95%)下养护28 d后进行纳米压痕试验.P·O 42.5水泥的主要化学组成如表1所示,粒径分布如图1所示.3个水灰比水泥净浆试样的28 d抗压强度分别为(97.8±3.7)、(53.0±1.6)、(22.0±1.1) MPa.

表1  P·O 42.5水泥的主要化学组成
Table 1  Main chemical composition of P·O 42.5 cement ( w/% )
SiO2CaOSO3Al2O3Fe2O3MgONa2OK2OTotal
22.75 50.56 0.79 10.38 10.37 1.02 0.73 1.31 97.91

图1  P·O 42.5水泥的粒径分布

Fig.1  Particle size distribution of P·O 42.5 cement

1.2 纳米压痕试验

1.2.1 计算原理

纳米压痕试验通过特定形状的压头压入材料表面,测量压入荷载-压入深度(Fh)曲线,基于接触力学分析材料的微观力学性能,计算其接触约化模量(Er)和压痕硬度(H

8

Er=π2βAcS (1)
H=FAc (2)
S=dFdhhmax (3)

式中:β是与压头形状有关的无量纲修正系数,当使用Berkovich压头时,β =1.05;S是压痕卸载刚度, N/m;Ac为最大压入深度(hmax)处的接触面积,nm2,是hmax处接触深度(hc)的函

9.

对于各向同性的均质材料,Er可根据式(4)计算.

1Er=1-v2E+1-vi2Ei (4)

式中:v为被测材料的泊松比;E为被测材料的弹性模量,GPa;vi为压头的泊松比;Ei为压头的弹性模量,GPa.

1.2.2 试验方法

首先,将养护至28 d的水泥净浆试样切割成10 mm的小块,放入异丙醇中浸泡2 d以终止其水化;然后,用环氧树脂进行镶嵌.高速高分辨率的纳米压痕力学绘图对试件表面平整度有很高的要求,因此需要精细的打磨抛光处理.首先,在自动磨抛机上用800#、1 200#、1 500#、2 000#、4 000#的砂纸上依次打磨2.5、5.0、10.0、30.0、30.0 min;然后,在高分子合成革抛光布上,采用磨料粒度为3.00、1.00、0.25、0.05 μm的金刚石悬浮液进行抛光处理,总抛光时间2.5 h.

纳米压痕测试采用iMicro Nanoindenter型纳米压痕仪,荷载和位移分辨率分别为6 nN和0.04 nm.试验采用Berkovich金刚石压头,其弹性模量和泊松比分别为1 141 GPa和0.07.采用NanoBlitz3D的高速纳米压痕方法,在该模式下,包含接触检测、加载、卸载在内的单个压痕测试过程所需时间不到1 s,从而可以在短时间内测得大量数据,实现材料的高速、高分辨率力学性能测绘.选取100×100点阵进行纳米压痕试验,点间距1 μm,压入荷载0.4 mN.每种水灰比测绘了3张空间分辨率为1 μm的100 μm×100 μm力学性能图.

1.3 聚类分析

聚类分析是一种无监督学习方法,常见的基于欧氏距离的聚类方法有Kmeans算法和Kmedoids算法.Kmedoids聚

22的计算原理是将数据划分成k个组,并选择组内一点为聚类中心:首先,初始化聚类中心,并将数据点分配到最近的聚类中心得到各聚类簇;然后,通过迭代不断更新聚类中心和聚类簇.设定E为所有数据组内点距离各自中心点的平方和,如式(5)所示.

E=i=1kxCidistx,oi2 (5)

式中:dist(xixj)为拥有p个维度的样本点xixj之间的欧式距离,即:

distxi,xj=a=1pxia-xja2 (6)

E达到最小时,得到k个簇C1C2,…Ck及其聚类中心o1o2,…ok.

Kmedoids聚类过程如图2所示.

图2  Kmedoids聚类过程示意图

Fig.2  Schematic diagram of Kmedoids clustering

在Kmedoids算法中,每个组的聚类中心是组内的数据点,因此比起将组内数据的均值当作聚类中心的Kmeans算法,Kmedoids算法的鲁棒性更好,不易受异常值的影响.与传统的解卷积方

10相比,聚类法还具有初值依赖性小、计算效率高等优势,适合处理大规模数据.此外,聚类算法可以获得各数据点的分类标签,便于建立数据点分相结果与位置坐标的关系,将力学性能图转化为矿相的几何分布图.

2 结果与分析

水泥基胶凝材料中的矿相十分复杂,包含水化硅酸钙(CSH)、氢氧化钙(CH)、水化硫铝酸钙(AFm和AFt)等水化产物,同时还有未水化水泥颗粒中的硅酸三钙(C3S)、硅酸二钙(C2S)、铝酸三钙(C3A)、铁铝酸四钙(C4AF)以及石膏

23.由于点阵纳米压痕试验方法具有统计性,因此可以对水泥基材料中性能相近的矿相进行简并处理,一般可划分为多孔复合相(PP)、低密度CSH相(LD CSH)、高密度CSH相(HD CSH)、CH相和超高性能相(UHP)1024‑25.其中:PP相是水泥石中毛细孔附近的疏松区域,UHP相主要包括未水化水泥颗粒的有关矿相,可简化认为UHP相代表了未水化的水泥颗粒相.

2.1 纳米压痕3D‑mapping测试结果

34分别为3个水灰比水泥净浆试样弹性模量和硬度的3Dmapping云图.其中:不同灰度代表不同的弹性模量、硬度值,归属于不同的矿相.由图34可见:弹性模量与硬度表征的矿相分布一致;随着水灰比的增加,代表力学性能较低的浅灰色区域面积逐渐增大,力学性能较高的深灰色区域面积逐渐减小,这从矿物组成的角度解释了高水灰比试块整体弹性模量和硬度减小的原因.

图3  3个水灰比水泥净浆试样弹性模量的3Dmapping云图

Fig.3  3Dmapping cloud map of elastic modulus for cement paste samples with three watercement ratios

图4  3个水灰比水泥净浆试样硬度的3Dmapping云图

Fig.4  3Dmapping cloud map of hardness for cement paste samples with three watercement ratios

2.2 矿相聚类

高速纳米压痕测试速度快,单个测试可在1 s内完成.同时,为了提高分辨率,点与点之间通常设置为较小间距.因为较高的应变率和较小的点间距,导致测量的材料力学参数的准确性受到一定影

12,因此一般不作为力学性能表征的准确值,而更多关注其力学性能的点阵分布及其所表达的矿相分布.Sudharshan26采用试验和有限元模拟研究压痕间距(d)与压入深度(h)的比值(d/h)对不同材料测试结果的影响,结果发现:当d/h>10时,材料的硬度保持稳定;当d/h=5时,材料的硬度变化范围在10%以内;当d/h>8时,SiO2、CaCl2、聚碳酸酯、Al、Ni、Cu、Mo的弹性模量和硬度偏差均小于10%.Němeček12采用有限元方法模拟应变率对材料压痕结果的影响,发现高速纳米压痕条件下(1.542 s-1)测得的弹性模量比标准模式(0.067 s-1)提高了约10%.在本文中,除部分PP相外,其余矿相的压痕深度都在150 nm以内,d/h值约为6.7.试验数据与真实值有一定的偏差,但不同相间的力学差异可以区分,所以矿相的划分是合理的.

按照Chen

21的计算步骤,对各水灰比水泥净浆试样压痕点的弹性模量和硬度进行Kmedoids聚类,包括PP相、LD CSH相、HD CSH相、CH相和UHP相,结果如图5所示.将聚类分析方法与传统解卷积方10的计算结果进行对比(见表2),发现聚类法具有与解卷积方法一致的准确度.在以往的研究中,Vandamme27对水灰比在0.15~0.40的5种水灰比试块进行了纳米压痕试验,解卷积结果表明,随着水灰比的增大,未水化水泥颗粒逐渐减少,LD CSH相的含量逐渐增加,当水灰比大于0.2时,增大水灰比会使HD CSH相的含量减少.魏亚28通过点阵压痕测试0.3、0.4、0.5这3种水灰比水泥净浆试块,解卷积结果表明水泥各矿相的力学性质不随水灰比变化.这与本文中Kmedoids聚类结果一致,体现了聚类结果的准确性.

图5  3个水灰比试样的Kmedoids聚类结果

Fig.5  Kmedoids clustering results of cement samples with three watercement ratios

表2  不同水灰比水泥净浆主要微观矿相解卷积结果和聚类结果
Table 2  Deconvolution results and clustering results of main mineral phases in cement paste with different w/c ratios
mW/mCPhaseDeconvolutionKmedoids clustering
E/GPaH/GPaVolume fractionE/GPaH/GPaVolume fraction
0.2 PP 13.69±2.04 0.35±0.21 0.06 13.71±2.57 0.46±0.14 0.10
LD CSH 20.02±2.47 0.69± 0.13 0.19 18.78±1.67 0.68± 0.13 0.17
HD CSH 28.49±5.99 1.15±0.33 0.47 27.94±6.00 1.15±0.37 0.48
CH 48.08±6.75 2.39±0.44 0.14 47.02±8.70 2.45±0.57 0.12
UHP 73.51±18.60 4.64±1.14 0.14 74.53±18.73 4.68±1.23 0.13
0.4 PP 12.14±3.22 0.45±0.15 0.13 12.22±2.99 0.47±0.19 0.15
LD CSH 19.32±2.82 0.78± 0.14 0.26 19.31±2.01 0.79± 0.16 0.25
HD CSH 26.16±3.95 1.17±0.24 0.36 27.35±4.28 1.25±0.32 0.39
CH 40.52±8.10 2.45±0.65 0.15 43.34±8.54 2.86±0.74 0.13
UHP 70.00±17.62 4.80±0.82 0.10 71.32±16.53 5.11±1.12 0.08
0.6 PP 11.90±3.29 0.48±0.19 0.28 11.49±2.92 0.47±0.18 0.27
LD CSH 18.51±3.19 0.89±0.21 0.30 18.76±2.92 0.80±0.21 0.31
HD CSH 25.91±3.81 1.43±0.29 0.24 26.23±4.21 1.42±0.27 0.24
CH 36.34±6.16 2.29±0.51 0.16 37.25±7.83 2.40±0.58 0.17
UHP 76.62±11.01 4.90±0.68 0.02 74.23±15.14 4.75±1.07 0.01

此外,聚类分析的另一优势在于其结果可直接将各个测点的几何坐标与其所归属的矿相相对应,从而得到矿相的几何分布,如图6所示.由图6可见:在0.2水灰比的水泥净浆试样中,UHP相代表的未水化水泥颗粒含量较高,低强度PP相的含量少;在0.6水灰比试样中,未水化水泥颗粒几乎消失,低强度PP相的含量较高;从几何分布的位置关系上,CH相、HD CSH相分布在未水化水泥颗粒相周围,而LD CSH相与PP相邻近;含量较高的LD和HD CSH相呈现连续分布,在复合材料性能计算中可以看作水泥净浆的基体相,而含量较少的PP、CH和UHP相以夹杂相形式散布在水泥净浆中.

图6  3个水灰比试样的矿相分布图

Fig.6  Main mineral phases mapping of cement samples with three watercement ratios

表2可见:

(1)0.2水灰比水泥净浆中UHP相较多而PP相较少,其体积分数分别为13%和10%;0.6水灰比水泥净浆中UHP相的体积分数只有约1%,PP相的体积分数较高,为27%.这说明高水灰比有利于水泥颗粒的充分水化,但同时生成了较多的孔隙,导致其强度降低.

(2)CH相的体积分数也与水泥的水化程度对应,0.2、0.4、0.6水灰比的水泥试样,CH相的体积分数分别是12%、13%、17%,随着水化程度的提高而上升.随着水灰比的增大,LD CSH相的体积分数增加,HD CSH相的体积分数减少.在0.2、0.4、0.6水灰比试样中,LD CSH相的体积分数分别为17%、25%、31%,HD CSH相的体积分数按48%、39%、24%依次递减.

2.3 矿相的尺寸分布

得到水泥净浆主要矿相的几何分布图后,可以对不连续夹杂相的尺寸分布进行统计分析,而CSH相由于含量高且呈连续性分布,暂时不进行统计分析.测量每个单独夹杂的面积(S),可以根据式(7)计算出等效直径(Deq),并统计出不同直径的夹杂数量.

Deq=4×Sπ (7)

水泥净浆的CH相在水泥颗粒周围生成,存在与其他相互相包含的情况,因此上述方法并不适用.此时,可采用图像细化方

29,在保留图像的拓扑特征的同时,提取二值图像的中心线.采用二值化和图像细化法,提取CH相骨架,用面积/骨架长度计算CH相的晶粒聚集层的厚度.

将PP相、UHP相的颗粒粒径采用对数正态分布对累积分布函数(CDF)(式(8))进行拟合,结果如图78所示,拟合优度(R2)均在0.98以上.

Fxμ,σ=                1σ2π0x1texp-logt-μ22σ2dt, x>0 (8)

图7  不同水灰比水泥净浆多孔相的尺寸分布

Fig.7  Size distribution of PP phase in cement paste with different watercement ratios

图8  不同水灰比水泥净浆未水化水泥颗粒相的尺寸分布

Fig.8  Size distribution of UHP phase in cement paste with different watercement ratios

由图78可以看出:

(1)随着水灰比的增大,PP相的平均尺寸也在逐渐增大,在0.2、0.4、0.6水灰比的水泥净浆试样中,其平均粒径分别为3.41、3.91、5.33 μm.UHP相代表的未水化水泥颗粒相在水灰比为0.2和0.4时的平均粒径接近,分别为4.57、4.50 μm,水灰比0.6时为2.56 μm.与图1所示水泥颗粒初始粒径集中分布在30 μm相比,未水化水泥颗粒的粒径大幅减小.

(2)在尺寸分布上,0.2水灰比UHP相的最大粒径27.97 μm,0.4、0.6水灰比UHP相的最大粒径分别为18.14、12.94 μm,大尺寸的未水化水泥颗粒的含量随着水灰比的增大而减小.采用图像细化计算时,在水灰比为0.2、0.4、0.6的试样中,CH聚集体的厚度分别为1.49、2.11、3.00 μm,说明CH产物随水灰比的增大而增多.

3 结论

(1)高速纳米压痕3Dmapping技术可以快速获得较高分辨率的水泥基材料力学性能图像.采用Kmedoids聚类分析处理3Dmapping云图,可实现矿相划分并得到矿相的几何分布,为深入分析水泥基材料的矿相性能及分布提供了一种新思路与新方法.

(2)典型水灰比水泥基材料3Dmapping云图的聚类分析结果与传统结论一致,随着水灰比的增加,水泥的水化程度提高,氢氧化钙(CH)相的含量增加;同时多孔相(PP)和低密度CSH相的含量逐渐增加,高密度CSH相和未水化水泥颗粒相的含量逐渐减少.

(3)夹杂相的尺寸分布符合对数正态分布.随着水灰比的增加,PP相、CH相聚集体的尺寸逐渐增大,未水化水泥相的含量逐渐减少.

参考文献

1

CONSTANTINIDES GULM F. The effect of two types of C‑S‑H on the elasticity of cement‑based materials: Results from nanoindentation and micromechanical modeling[J]. Cement and Concrete Research2004341):6780. [百度学术] 

2

MAEKAWA KISHIDA TKISHI T. Multi‑scale modeling of concrete performance:Integrated material and structural mechanics[J]. Journal of Advanced Concrete Technology200312):91126. [百度学术] 

3

TAYLOR H F W. Cement chemistry[M]. LondonThomas Telford1997. [百度学术] 

4

WEI YLIANG SGAO X. Phase quantification in cementitious materials by dynamic modulus mapping[J]. Materials Characterization2017127348356. [百度学术] 

5

赵素晶孙伟. 纳米压痕在水泥基材料中的应用与研究进展[J]. 硅酸盐学报2011391):164176. [百度学术] 

ZHAO SujingSUN Wei. Application and research progress of nanoindentation in the field of cement‑based materials[J]. Journal of the Chinese Ceramic Society2011391):164176. (in Chinese) [百度学术] 

6

李一凡管学茂刘松辉. 压痕点数及解卷积法对水泥纳米压痕试验的影响[J]. 建筑材料学报2021242):291295. [百度学术] 

LI YifanGUAN XuemaoLIU Songhuiet al. Effect of number of indentation points and deconvolution method on cement nanoindentation test[J]. Journal of Building Materials2021242):291295. (in Chinese) [百度学术] 

7

刁益彤麻海燕余红发. 珊瑚混凝土骨料净浆界面区强度[J]. 建筑材料学报2021243):516524. [百度学术] 

DIAO YitongMA HaiyanYU Hongfaet al. Strength of aggregate‑paste interface transition zone of coral concrete[J]. Journal of Building Materials2021243):516524. (in Chinese) [百度学术] 

8

OLIVER W CPHARR G M. An improved technique for determining hardness and elastic modulus using load and displacement sensing indentation experiments[J]. Journal of Materials Research199276):15641583. [百度学术] 

9

OLIVER W CPHARR G M. Measurement of hardness and elastic modulus by instrumented indentation:Advances in understanding and refinements to methodology[J]. Journal of Materials Research2004191):320. [百度学术] 

10

CONSTANTINIDES GULM F J. The nanogranular nature of C‑S‑H[J]. Journal of the Mechanics and Physics of Solids2007551):6490. [百度学术] 

11

ZHU W ZHUGHES J JBICANIC Net al. Nanoindentation mapping of mechanical properties of cement paste and natural rocks[J]. Materials Characterization20075811/12):11891198. [百度学术] 

12

NĚMEČEK JLUKEŠ JNĚMEČEK J. High‑speed mechanical mapping of blended cement pastes and its comparison with standard modes of nanoindentation[J]. Materials Today Communications202023100806. [百度学术] 

13

VRANJES‑WESSELY SMISCH DKIENER Det al. High‑speed nanoindentation mapping of organic matter‑rich rocks:A critical evaluation by correlative imaging and machine learning data analysis[J]. International Journal of Coal Geology2021247103847. [百度学术] 

14

VIGNESH BOLIVER W CKUMAR G Set al. Critical assessment of high speed nanoindentation mapping technique and data deconvolution on thermal barrier coatings[J]. Materials and Design2019181108084. [百度学术] 

15

HINTSALA ERISAN JDIETRICH Ret al. Correlated EBSD and high speed nanoindentation mapping[J]. Microscopy and Microanalysis201723Suppl 1):780781. [百度学术] 

16

JANAKIRAM SPHANI P SUMMETHALA Get al. New insights on recovery and early recrystallization of ferrite‑pearlite banded cold rolled high strength steels by high speed nanoindentation mapping[J]. Scripta Materialia2021194113676. [百度学术] 

17

MUGHAL M ZMOSCATELLI RSEBASTIANI M. Load displacement and high speed nanoindentation data set at different state of charge (SoC) for spinel LixMn2O4 cathodes[J]. Data in Brief20168203206. [百度学术] 

18

ROA J JSUDHARSHAN PHANI POLIVER W Cet al. Mapping of mechanical properties at microstructural length scale in WC‑Co cemented carbides:Assessment of hardness and elastic modulus by means of high speed massive nanoindentation and statistical analysis[J]. International Journal of Refractory Metals and Hard Materials201875211217. [百度学术] 

19

CUY J LMANN A BLIVI K Jet al. Nanoindentation mapping of the mechanical properties of human molar tooth enamel[J]. Archives of Oral Biology2002474):281291. [百度学术] 

20

JAIN A KMURTY M NFLYNN P J. Data clustering:A review[J]. ACM Computing Surveys1999313):264323. [百度学术] 

21

CHEN X WHOU D WHAN Y Det al. Clustering analysis of grid nanoindentation data for cementitious materials[J]. Journal of Materials Science20215621):12238‑12255. [百度学术] 

22

PARK HJUN C. A simple and fast algorithm for K‑medoids clustering[J]. Expert Systems with Applications2009362):33363341. [百度学术] 

23

LOTHENBACH B. Thermodynamic equilibrium calculations in cementitious systems[J]. Materials and Structures20104310):14131433. [百度学术] 

24

ULM F JVANDAMME MBOBKO Cet al. Statistical indentation techniques for hydrated nanocomposites:Concrete, bone, and shale[J]. Journal of the American Ceramic Society2007909):26772692. [百度学术] 

25

HU C L. Nanoindentation as a tool to measure and map mechanical properties of hardened cement pastes[J]. MRS Communications201551):8387. [百度学术] 

26

SUDHARSHAN PHANI POLIVER W C. A critical assessment of the effect of indentation spacing on the measurement of hardness and modulus using instrumented indentation testing[J]. Materials and Design2019164107563. [百度学术] 

27

VANDAMME MULM FFONOLLOSA P. Nanogranular packing of C‑S‑H at substochiometric conditions[J]. Cement and Concrete Research2010401):1426. [百度学术] 

28

魏亚高翔梁思明. 硬化水泥净浆基于纳米压痕的相态识别与水化程度计算[J]. 复合材料学报2017345):11221129. [百度学术] 

WEI YaGAO XiangLIANG Siming. Nanoindentation‑based study of the microstructure and degree of hydration in hardened cement paste[J]. Acta Materiae Compositae Sinica2017345):11221129. (in Chinese) [百度学术] 

29

LEE TKASHYAP R L. Building skeleton models via 3‑D medial surface axis thinning algorithms[J]. CVGIP:Graphical Models and Image Processing1994566):462478. [百度学术]