金属学报, 2021, 57(12): 1595-1606 DOI: 10.11900/0412.1961.2020.00432

研究论文

高碳钢连铸坯大区域C元素分布不均匀度

郭中傲1,2, 彭治强1,2, 柳前1,2, 侯自兵,1,2

1.重庆大学 材料科学与工程学院 重庆 400044

2.重庆大学 钒钛冶金及新材料重庆市重点实验室 重庆 400044

Nonuniformity of Carbon Element Distribution of Large Area in High Carbon Steel Continuous Casting Billet

GUO Zhongao1,2, PENG Zhiqiang1,2, LIU Qian1,2, HOU Zibing,1,2

1.College of Materials Science and Engineering, Chongqing University, Chongqing 400044, China

2.Chongqing Key Laboratory of Vanadium-Titanium Metallurgy and New Materials, Chongqing University, Chongqing 400044, China

通讯作者: 侯自兵,houzibing@cqu.edu.cn,主要从事凝固组织与偏析特征及精细化控制研究

收稿日期: 2020-10-30   修回日期: 2020-12-02   网络出版日期: 2021-11-18

基金资助: 国家自然科学基金委员会-中国宝武钢铁集团有限公司钢铁联合研究基金项目.  U1860101

Corresponding authors: HOU Zibing, associate professor, Tel:(023)65105202, E-mail:houzibing@cqu.edu.cn

Received: 2020-10-30   Revised: 2020-12-02   Online: 2021-11-18

作者简介 About authors

郭中傲,男,1996年生,硕士生

摘要

基于典型高碳钢铸坯低倍组织灰度图,引入标准差、差分盒维数和惯性矩,探讨一种大区域主要偏析元素(C)分布不均匀度定量表征的简易方法。结果表明:标准差、差分盒维数和惯性矩可以有效地综合表征C元素分布不均匀度;且等轴晶区C元素分布不均匀度高于柱状晶区,平均相差20.85%。对比发现,标准差主要基于灰度值统计特性,差分盒维数和惯性矩结合了灰度值统计信息和空间分布信息;并且差分盒维数具有尺度独立性,受灰度图尺寸和分辨率的影响较小,而惯性矩对微区C元素分布不均匀度变化会更加敏感。此外,标准差主要受到大型偏析点(> 1 mm2)影响,而差分盒维数和惯性矩主要受中型偏析点(0.1~1 mm2)影响。本工作可为高碳钢铸坯大区域C元素分布不均匀度的全面衡量及精细化质量评判提供新的参考方法。

关键词: 偏析 ; 不均匀度 ; 低倍组织 ; 连铸 ; 高碳钢

Abstract

High carbon steel is prone to macroscopic and semimacroscopic segregation. Segregation is a phenomenon that results in the nonuniformity of solute element distribution during solidification; segregation severely impacts the quality of the casting billet. Based on qualitative analysis (macroscopic quality rating) and relatively simple quantitative analysis (segregation index, mean square error), it is found that the accuracy of existing technologies is limited for determining the nonuniformity of carbon element distribution in large areas of cast high carbon steel billets; this causes difficulties while applying such existing technologies for the evaluation of large samples under actual steelworks usage conditions. Therefore, it is essential to study the nonuniformity of carbon element distribution to precisely evaluate and optimize the quality of high carbon steel. Based on the grayscale image of the casting blank macrostructure of typical high carbon steels (mass fraction of carbon = 0.7%), the standard deviation, differential box-counting, and moment of inertia are introduced to discuss simple methods for the quantitative characterization of the nonuniformity of the major segregation element (C) in large areas. Then, the validity of the characterization results was verified using statistical homogeneity and the near-equilibrium solidification model, and the similarities and differences of the three characterization methods were discussed. The results showed that the standard deviation, differential box-counting, and moment of inertia can effectively reflect the nonuniformity of carbon element distribution. Furthermore, the mean value of the nonuniformity of carbon element distribution in the equiaxed area was 20.85% higher than that in the columnar area. The standard deviation is mainly based on the statistical characteristics of the grayscale value, while differential box-counting and moment of inertia combine the grayscale value statistical and spatial distribution information on a grayscale image of the casting blank macrostructure. In addition, differential box-counting has scale independence, which is less affected by the size and resolution of the grayscale image; the moment of inertia is sensitive to the variation of the nonuniformity of carbon element distribution in the microregion. Finally, the standard deviation is mainly influenced by large segregation points (e.g., area larger than 1 mm2), while differential box-counting and the moment of inertia are mainly influenced by medium segregation points (e.g., area of 0.1-1 mm2 ). Therefore, the three characterization methods can be combined to comprehensively evaluate the nonuniformity of carbon element distribution in large areas of cast high carbon steel billets. This research can provide a new theoretical reference for comprehensively evaluating the nonuniformity of carbon element distribution in a large area and for the fine quality evaluation of cast high carbon steel billets.

Keywords: segregation ; nonuniformity ; macrostructure ; continuous casting ; high carbon steel

PDF (2952KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

郭中傲, 彭治强, 柳前, 侯自兵. 高碳钢连铸坯大区域C元素分布不均匀度. 金属学报[J], 2021, 57(12): 1595-1606 DOI:10.11900/0412.1961.2020.00432

GUO Zhongao, PENG Zhiqiang, LIU Qian, HOU Zibing. Nonuniformity of Carbon Element Distribution of Large Area in High Carbon Steel Continuous Casting Billet. Acta Metallurgica Sinica[J], 2021, 57(12): 1595-1606 DOI:10.11900/0412.1961.2020.00432

高碳钢具有高的强度和硬度,高的弹性极限和疲劳极限,被广泛应用于船舶、化工、机械等领域[1,2]。由于其C含量高、溶质分配系数低[3],导致高碳钢的连铸过程凝固区间大,易产生偏析和缩孔等内部质量缺陷,因此C元素是高碳钢中的主要偏析元素。偏析是合金凝固过程中溶质元素在固相和液相中再分配的结果,表现为从铸坯表面到中心或沿铸坯轴向大区域范围(如厘米范围以上)元素分布的不均匀性[4]。偏析元素分布不均匀现象易使连铸坯产生裂纹或失效[5~7]等质量问题,从而降低最终产品的力学性能。因此,为了更精细地评价与优化高碳钢产品质量,需要对其大区域C元素分布不均匀度展开专门研究。

以往,研究者[8,9]对C元素分布不均匀度的研究主要从形貌和成分含量2个方面进行。前者一般是通过热酸洗获取铸坯的偏析形貌,然后根据偏析区域的尺寸、数目的分布特征以及与宏观评级图的对比来定性判断偏析元素分布不均匀度;对于后者,主要基于钻孔取样化学分析法[10]、电子探针[11]和金属原位分析仪[12]等方法获取元素含量,然后通过偏析指数[13]、均方差[14]等对铸坯的元素分布不均匀度进行定量评判。但是,以上的定性方法易造成较大偏差;而定量方法中钻孔取样化学分析法往往是通过局部某几点元素含量反映整体元素分布不均匀度,且电子探针等电子显微方法测量区域太小,而金属原位分析仪试样准备相对繁琐,对C含量测量准确度有限,难以应用于实际钢厂大试样的评判。因此,针对不同组织区域内元素分布不均匀度的判断可能会存在不同的结果。例如,有学者[15]认为等轴晶区内部质量更好,成分分布更均匀,希望增加等轴晶比例。然而,另一些研究[16]指出柱状晶的偏析比等轴晶要少,结构更为致密,希望得到定向凝固组织。可见,柱状晶区与等轴晶区元素分布不均匀度高低的问题尚存在争议,这使得人们无法精确评价铸坯不同组织的质量。目前定量表征高碳钢铸坯大区域C元素分布不均匀度仍是一个难题,其一是难以获取大区域C元素含量分布;其二是没有具体方法衡量大区域C元素分布不均匀度。另外,二维大区域元素分布不均匀度不是一个单一的概念,其包括整体平均不均匀度、局部分布不均匀度等方面;同时元素分布不均匀度会因数据尺度与分析范围的不同而改变。因此,有必要寻找一套有效且能综合表征铸坯大区域范围内偏析元素分布不均匀度的实用方法。本工作提出的大区域偏析元素分布不均匀度和宏观偏析虽均指宏观区域化学成分不均匀性,但不同点在于宏观偏析一般侧重于描述较大尺度范围内某一区域(如元素富集形成的偏析带)与原始钢水中的成分差异(如中心宏观偏析);而大区域偏析元素分布不均匀度更倾向于描述较大尺度范围内所有位置的化学成分分布的不均匀性。

本文作者前期工作[17]针对高碳钢铸坯低倍组织灰度图进行研究,发现低倍组织灰度值与C含量存在较明显的负相关关系,即灰度值越大,C含量越低。同时,利用标准热盐酸酸蚀实验所得到的低倍组织灰度图,易近似获得大区域范围内C含量分布,从而为分析不同组织区域内C元素分布不均匀度提供基础。虽然低倍组织不同位置灰度值可能并不能完全准确表示高碳钢中主要偏析元素(C元素)的含量,但基于实际低倍组织灰度图的分析可一定程度上客观地反映C含量的不均匀分布特征;并且可以根据低倍组织灰度图尺寸方便且灵活地实现不同范围大区域的C元素分布不均匀度表征。另外,基于灰度图进行研究,目前已被应用于材料研究的相关领域,例如估计金属的腐蚀程度[18]、定量金相分析[19]和金属表面缺陷检测[20,21]等。但是,将高碳钢铸坯低倍组织灰度图应用于C元素分布不均匀度的研究还鲜见报道。

本工作以含C量为0.7% (质量分数)的高碳钢连铸方坯为例,尝试通过铸坯低倍组织灰度图探讨一种简易且系统表征高碳钢铸坯大区域C元素分布不均匀度的有效方法。首先,基于低倍组织灰度图,采用标准差(standard deviation)、差分盒维数(differential box-counting)以及惯性矩(moment of inertia)来综合定量表征铸坯大区域C元素分布的不均匀度;然后,结合统计均匀度和近平衡凝固模型对C元素分布不均匀度表征结果的有效性进行验证;最后,探讨3种表征方法的异同。

1 实验方法

本实验采用尺寸为170 mm × 170 mm的70钢连铸方坯,其主要成分(质量分数,%)为:C 0.7,Si 0.2,Mn 0.65,P 0.011,S 0.0023,Fe余量。为了获取铸坯低倍组织灰度图,从铸坯中心纵断面取样并进行标准热酸洗实验。热酸洗实验中酸液为容积比1∶1的工业盐酸水溶液,酸蚀温度60~80℃,酸蚀时间25 min。热酸洗实验后得到整体中心纵断面低倍组织,然后根据图1所示位置从铸坯表面到中心依次选取6个尺寸为10 mm × 10 mm的小方格(1#~6#),在同样光源条件下使用NVC4030型二维影像仪获取每个小方格的低倍组织图像。在得到原始的低倍组织图像后,利用Matlab2018a将原始图像转换为灰度图,灰度图尺寸为256 px × 256 px。

图1

图1   铸坯中心纵断面取样位置示意图

Fig.1   Schematic of sampling location in the central plane of billet (unit: mm)


对于高碳钢铸坯大区域C元素分布不均匀度不同角度的定量表征,本工作从标准差、差分盒维数和惯性矩等方法展开。首先,灰度值标准差通过计算与平均灰度值的波动差值反映铸坯C元素分布的整体平均不均匀度。其计算公式为[22]

St=u=0L-1zu-m2P(zu)

式中,St为标准差;zu为图像灰度值的随机变量;P(zu)为对应区域中灰度值的像素分布概率(u = 0, 1, ..., L - 1);m为灰度图的灰度均值;L为可区分的灰度级数。

分形理论是研究大自然不规则复杂对象(例如三维形貌)的一种有效工具[23]。由于灰度图像的灰度值分布是表面真实三维形貌的直接反应,因此基于表面灰度图,分形理论已被应用于各种材料表面三维形貌的定量表征,如材料磨损表面[24]、材料断裂表面[25]和材料腐蚀形貌[26]等。分形维数作为分形的重要特征和度量,将图像的空间信息和灰度信息简单而有机地结合起来。分形理论的核心问题是如何准确地估计分形维数,前人已经提出了一些计算灰度图分形维数的方法[26],不同的计算方法从不同的角度对分形对象进行描述,可以得到不同的分形维数值。其中由Sarkar和Chaudhuri[27]提出的差分盒维数具有动态取值稳定和计算效率高等优点,在计算灰度图的分形维数被广泛应用,本工作也使用差分盒维数计算铸坯组织灰度图的分形维数。其基本原理为:对于给定区域为M × M的图像,假设它己被分解成S × S的小块,S是介于1~M/2之间的一个整数。此时尺度因子r = S / M。可把M × M的灰度图像看成是一个三维空间,(x, y)表示各点的平面位置,第三维用来表示像素点的灰度值。当把图像平面划分成若干S × S的网格,把第三维图像的灰度值划分成长度为S'的小段,则在每一个网格上存在一系列尺寸为S × S × S'的小盒子。设图像的灰度级为L,取S'满足[L / S'] = [M / S] = l / r。例如,如图2所示,S = S' = 3,组织灰度图平面的第(i, j)个网格内灰度值最高和最低的像素点分别落在第5和第2个盒子,则覆盖这个网格需要3个盒子。假设在图像的第(i, j)个网格内灰度值最高和最低的像素点分别落在第h和第k个盒子内,则尺度因子为r时的总盒子数N(r)在第(i, j)个网格内的分布nr (i, j)由式(2)给出:

nr(i,  j)=h-k+1

图2

图2   差分盒维数算法示意图

Fig.2   Schematic of differential box-counting algorithm (The size of a box is 3 × 3 × 3, which is used to cover gray value)


对所有的网格计算其分布,得到:

N(r)=inr(i,  j)

对应于不同的r,可以得到不同的N(r)值。由分形维数的计算方法,采用最小二乘法拟合出lnN(r)-ln(1/r)的斜率,即为对应的铸坯组织灰度图的差分盒维数(D),计算公式如下:

D=lnN(r)ln(1 / r)

根据差分盒维数的计算原理可知,差分盒维数通过不同尺度下(局部到整体)的灰度最大值与最小值之差反映C元素分布不均匀程度。同时分形维数作为具有自相似性和无标度性特征的表征参数,具有尺度独立性[26]

Haralick等[28]在20世纪70年代初期提出的灰度共生矩阵(GLCM)是用来分析图像纹理特征的经典二阶统计法,其反映图像纹理关于方向、相邻间隔、变化幅度的综合信息,目前已被应用于材料研究的许多相关领域,例如金属形貌特征提取[29,30]和断口形貌分析[25]等。灰度共生矩阵的基本原理是:取M × M的图像灰度值为i的像素点(x, y)及偏离它的另一点灰度值为j的像素点(x + a, y + b),其中,ab分别为两像素点间的水平距离和垂直距离。若(x, y)与(x + a, y + b)间距为d,2者与坐标横轴的的夹角为θ,统计(i, j)同时出现的频率,记为Pδ (i, j)。其中δij的某个固定位置关系(d, θ),ij = 0, 1, 2, ..., n-1, n表示灰度值,称由Pδ (i, j)组成的矩阵为灰度共生矩阵。通过灰度共生矩阵可以分析图像灰度值分布的局部模式和排列规则,但一般不直接应用得到的共生矩阵,而是在此基础上获取特征参数,灰度共生矩阵特征参数惯性矩度量图像的局部灰度值变化总量,可以作为衡量灰度值局部分布不均匀程度的指标。惯性矩(CON)的计算公式为:

CON=i=0L-1j=0L-1(i-j)2Pδ(i, j)

式中,(i, j)为图像上灰度值为ij的组合。

2 实验结果与讨论

2.1 低倍组织灰度图

根据热盐酸酸蚀原理,酸液与高碳钢铸坯接触发生反应时,偏析区域C含量高其Gibbs自由能大,导致与酸液反应剧烈且颜色变暗,即可通过铸坯低倍组织灰度图中的亮暗区域区分偏析与凝固组织[31]。根据实际生产现场广泛使用的标准热盐酸酸蚀实验,易得到铸坯大区域范围的低倍组织灰度图,如图3所示。从铸坯表面向中心凝固组织呈现不同形貌,即从柱状晶向等轴晶逐渐转变,其中1#~3# (图3a~c)位置主要为柱状晶形貌,4#~6# (图3d~f)位置主要为等轴晶形貌。根据GB/T 226-1991《钢的低倍组织及缺陷酸蚀检验法》与前期研究[17,31],图中灰度值较大的部分(亮度较高)为凝固组织,灰度值较低的部分(亮度较暗)为偏析区域;且灰度值与C含量近似呈现较好的负相关关系。故接下来将利用低倍组织灰度图来分析铸坯大区域C元素分布不均匀度变化。

图3

图3   铸坯不同位置的低倍组织灰度图

Fig.3   Grayscale images of macrostructures of samples 1#-6# at different locations of billet (a-f)


使用Matlab2018a编程绘制铸坯不同位置低倍组织灰度图的灰度曲面,即把灰度图放置在XOY平面上,像素的灰度值作为Z轴上的值,由此构成的曲面叫做图像的灰度曲面,灰度曲面即近似为铸坯C元素含量的空间分布。图4为铸坯内部不同位置的低倍组织灰度曲面。可以看出,每个位置的低倍组织灰度曲面无论从整体或局部进行观察,均可以看出灰度值分布非常不均匀。由于灰度值分布不均匀是铸坯C元素含量分布不均匀的一种表现,故可以利用高碳钢低倍组织灰度图这种简易方法来定性比较不同位置不均匀度的差别,但分布不均匀现象非常复杂,仅凭借定性观察很难对比不同位置灰度值分布的不均匀程度。因此很有必要借助于上述方法准确定量表征铸坯不同位置灰度值分布不均匀度的不同方面,并以此来研究高碳钢铸坯大区域C元素分布不均匀度。

图4

图4   铸坯不同位置的低倍组织灰度曲面

Fig.4   Grayscale surfaces of macrostructures of samples 1#-6# at different locations of billet (a-f) (The unit of X and Y axes is pixel)


2.2 C元素分布不均匀度表征结果

2.2.1 标准差

图5为铸坯不同位置的灰度值标准差变化。可见,从铸坯表面到中心,灰度值标准差逐渐增大;同时灰度值标准差在柱状晶区(1#~3#)和等轴晶区(4#~6#)变化趋势一致,但柱状晶区灰度值标准差增长幅度更大,等轴晶区增长幅度较为平缓。图3中从铸坯表面到中心凝固组织形貌依次为柱状晶和等轴晶。灰度值标准差通过衡量选定区域内灰度值与平均灰度值的波动差值表征C元素分布整体平均不均匀度。因此,等轴晶区C元素分布整体平均不均匀度相较于柱状晶区更高。同时,通过计算等轴晶区和柱状晶区标准差均值之差与柱状晶区标准差均值的比值可知,等轴晶区C元素分布整体平均不均匀度比柱状晶区高25.38%。

图5

图5   铸坯不同位置的灰度值标准差

Fig.5   Standard deviation of grayscale value at different locations of billet


2.2.2 差分盒维数

根据差分盒维数的计算原理,使用Matlab2018a编程对铸坯不同位置的低倍组织灰度图的差分盒维数进行计算。图6所示为铸坯低倍组织1#试样的lnN(r)-ln(1 / r)拟合关系。图中分布点拟合直线的斜率为2.6177,拟合系数R2为0.9946。根据差分盒维数的计算原理可知,1#试样的低倍组织灰度图的差分盒维数是其拟合直线的斜率,即2.6177。根据同样的方法可以得到铸坯不同位置的低倍组织灰度图差分盒维数及对应的拟合系数,具体数值见图7表1。可以看出,铸坯不同位置的拟合系数均接近1,这表明低倍组织灰度图具有很好的分形特征,引入差分盒维数这种方法来定量描述铸坯不同位置灰度值分布特征是有效的。差分盒维数通过不同尺度下的灰度最大值与最小值之差反映不同尺度下的C元素分布不均匀程度。对于图像的灰度分布曲面,差分盒维数介于2~3之间[27]。差分盒维数越接近于2,灰度分布曲面越光滑;差分盒维数越接近3,灰度分布曲面越粗糙。因此,当差分盒维数越大时,代表此位置不同尺度下的C元素分布不均匀程度越大,即此位置从局部到整体C元素分布不均匀程度越大。

图6

图6   铸坯低倍组织1#试样lnN(r)-ln(1 / r)关系

Fig.6   Relationship between lnN(r) and ln(1 / r) of 1# sample (R2—fitting coefficients, r—size coefficient, N(r)—number of total boxes)


图7

图7   铸坯不同位置的差分盒维数(D)

Fig.7   Differential box-counting (D) at different locations of billet


表1   不同位置的低倍组织灰度图的D及对应的拟合系数(R2)

Table 1  D and corresponding R2 at different locations

Sample No.DR2
1#2.61770.9946
2#2.66220.9950
3#2.67860.9946
4#2.68200.9943
5#2.71400.9946
6#2.67910.9946

新窗口打开| 下载CSV


图7为铸坯不同位置的组织灰度图的差分盒维数变化。可见,从铸坯表面到中心组织灰度图差分盒维数先逐渐增大,而后在中心等轴晶区发生波动。图3中从铸坯表面到中心凝固组织形貌依次为柱状晶和等轴晶。差分盒维数在柱状晶区(1#~3#)逐渐增大,在等轴晶区(4#~6#)存在波动。差分盒维数可以定量表征不同尺度下的C元素分布不均匀程度。因此,等轴晶区不同尺度下C元素分布不均匀程度相较于柱状晶区更高。同时,通过计算等轴晶区和柱状晶区差分盒维数均值之差与柱状晶区差分盒维数均值的比值可知,等轴晶区不同尺度下C元素分布不均匀程度比柱状晶区高1.57%。

2.2.3 惯性矩

本工作主要通过惯性矩定量表征微区C元素分布不均匀度,其中,间距d取1,夹角θ分别为0°、45°、90°、135°构建灰度共生矩阵,并且以惯性矩均值作为衡量微区C元素分布不均匀度。图8为惯性矩表征局部微区C元素分布不均匀度示意图。可见,所构建微区对应的组织灰度值分布,即为3 × 3的灰度值矩阵。根据惯性矩的计算公式(5)可知,通过累加组织灰度图上所有3 × 3的灰度值矩阵的中心值i与其相邻值j的差值,反映微区C元素分布不均匀度。惯性矩越大,即累加得到的组织灰度图上所有3 × 3的灰度值矩阵的中心值i与其相邻值j的差值越大,微区C元素分布不均匀度越大;惯性矩越小,即累加得到的组织灰度图上所有3 × 3的灰度值矩阵的中心值i与其相邻值j的差值越小,微区C元素分布不均匀度越小。

图8

图8   惯性矩表征微区C元素分布不均匀度示意图

Fig.8   Schematic of the nonuniformity of carbon micro distribution characterized by moment of inertia (d, θ—distance and angle between i and j, respectively )


图9为铸坯不同位置的组织灰度图惯性矩变化。可见,从铸坯表面到中心,惯性矩先逐渐增大,而后在中心等轴晶区发生波动。图3中从铸坯表面到中心凝固组织形貌依次为柱状晶和等轴晶。因此,惯性矩在柱状晶区(1#~3#)逐渐增大,在等轴晶区(4#~6#)存在波动。惯性矩定量表征微区C元素分布不均匀度。因此,等轴晶区C元素分布不均匀程度相较于柱状晶区更高。同时,通过计算等轴晶区和柱状晶区惯性矩均值之差与柱状晶区惯性矩均值的比值可知,等轴晶区微区C元素分布不均匀度相较于柱状晶区高35.61%。

图9

图9   铸坯不同位置的惯性矩

Fig.9   Moment of inertia at different locations of billet


2.3 C元素分布不均匀度表征结果有效性讨论

灰度直方图是将低倍组织灰度图的所有像素,按照灰度值的大小,统计其出现的频数。低倍组织灰度图可以近似代表区域内C元素空间分布趋势,因此灰度直方图可以认为是关于C元素分布的函数,是对区域内C元素分布的统计。

图10为铸坯不同位置的灰度直方图。可见,1#、2#位置的灰度直方图,其峰型较窄,说明C元素含量分布较为集中,表明这2个位置的C元素分布比较均匀;而对于3#~6#位置处的灰度直方图,峰型较宽,说明C元素含量分布较为分散,则表示C元素分布比较不均匀。通过不同位置的灰度直方图可以较为直观地衡量不同位置C元素分布不均匀度。根据王海舟等[32]提出的连铸坯原位统计分布的分析研究,通过对各原位置元素成分的权重比率统计分布进行解析,可得到某元素在特定含量区间(如允许差范围内)所占的权重比率(或称所占的频率),用以定量表征其统计均匀度(H),反映大区域内元素定量统计分布规律,其数值越大,表明在此材料中该元素分布的均匀性越好。H即为在元素平均含量的国家标准规定允许差范围内所占频率。根据GB/T 222-2006《钢的成品化学成分允许偏差》可知,70钢C含量允许差范围为0.7% ± 0.04%。灰度直方图可以代表铸坯各位置C元素含量的权重比率统计分布,因此可以以铸坯不同位置的灰度值均值g0作为C元素平均含量,根据70钢C含量允许差范围计算出对应灰度值允许差范围为(1 ± 5.714%)g0,然后计算出H表2为铸坯不同位置的灰度值允许差范围和H。可见,1#~3#位置的统计均匀度较明显大于4#~6#。图3中从铸坯表面到中心凝固组织形貌依次为柱状晶和等轴晶,因此,等轴晶的C元素分布不均匀度大于柱状晶区。同时,图11分别为统计均匀度随标准差、差分盒维数以及惯性矩的变化。可见,随着标准差、差分盒维数、惯性矩的增加,统计均匀度均明显减小。

图10

图10   铸坯不同位置的灰度直方图

Fig.10   Gray histograms of samples 1#-6# at different locations of billet (a-f)


表2   灰度值在含量允许差范围内所占的权重比率(统计均匀度(H))

Table 2  The weight ratios of grayscale values in the range of content tolerance

Sample No.H / %

Grayscale value in the range

of content tolerance

1#11.98[97,110]
2#9.76[103,116]
3#8.39[115,130]
4#7.35[121,136]
5#6.02[121,136]
6#7.33[118,133]

Note:H—degree of statistical homogeneity

新窗口打开| 下载CSV


图11

图11   统计均匀度与标准差、差分盒维数、惯性矩的关系

Fig.11   Relationship between H and standard deviation (a), differential box-counting (b), and moment of inertia (c)


以上表明,标准差、差分盒维数和惯性矩的表征结果与统计均匀度计算方法的结果是基本一致的。但与此同时,图11并没有呈现出完全负相关的关系,这表明本工作所提出的3种表征方法对不均匀度从不同角度进行了更为精细化的表达。

连铸坯在凝固过程中溶质元素分布除了受到温度梯度、浓度梯度、凝固组织、液相流动等方面的影响,更主要的是受到凝固过程固/液界面溶质再分配的影响。因此,为了有效描述高碳钢铸坯从表面到中心不同位置溶质元素再分布,可通过偏析元素固相无扩散且液相有限扩散的近平衡凝固模型[33~35]对其凝固过程进行分析。Burton等[35]提出近平衡凝固过程中固/液界面处的溶质有效分配系数(kE)可用式(6)表示,而此时固/液界面液相的溶质浓度可用式(7)表示:

kE=CS*C0=k0k0+(1-k0)exp-vDLδN

式中,DL是溶质扩散系数;C0是合金的原始成分;k0是溶质平衡分配系数;CS*是固/液界面前沿固相中的溶质浓度;v为枝晶生长速率;δN为扩散层厚度。

CL=C0(1-fs)(kE-1)

式中,fs是固相体积分数;CL是固/液界面液相中C的质量分数。

式(6)可以看出,溶质有效分配系数主要受到枝晶生长速率和扩散层厚度2方面因素影响。高碳钢连铸坯由表层柱状晶向中心等轴晶区凝固过程中,固/液界面前沿的枝晶生长速率逐渐减小,同时中心等轴晶区由于凝固后期中心线的抽吸收缩作用导致枝晶间液相流动速率增加,从而导致δN减小;以上2个因素均会导致kE减小。kE的减小将导致铸坯从表面向中心凝固过程中不同位置(1#~6#)液相溶质浓度分布发生较大程度的变化,从而影响C元素分布不均匀度。

图12为不同kECLfs的关系。可见,随着fs的增加,CL呈指数递增;同时,kE越小,CL由小到大变化的剧烈程度(切线斜率)越大,也即CL波动越大,导致C元素分布不均匀程度增加;另外,kE越小,CL越大,导致偏析程度增加。当凝固由外往内逐渐发生时,kE逐渐减小,从而C元素分布不均匀程度和偏析程度整体均呈现递增趋势;且C元素分布不均匀度越大,对应区域内偏析程度越严重。图3中从铸坯表面到中心凝固组织形貌依次为柱状晶区和等轴晶区,因此,本工作70高碳钢连铸坯中等轴晶区的C元素分布不均匀度相对于柱状晶区更高;且基于标准差、惯性矩和差分盒维数的计算结果,其不均匀度相差比例的平均值为20.85%。

图12

图12   液相中C的质量分数(CL)与固相体积分数(fs)的关系

Fig.12   Relationship between mass fraction of carbon in the liquid phase (CL) and solid fraction (fs) (kE—effective distribution coefficient of solute)


2.4 标准差、差分盒维数和惯性矩的对比讨论

上述研究表明,标准差、惯性矩和差分盒维数都能有效地反映大区域C元素分布不均匀度,以及铸坯从表面向中心凝固过程中不均匀度变化的整体趋势。根据3种方法的原理和计算方法可知,标准差是基于低倍组织灰度图上所有灰度值与灰度均值的统计特性计算得出的结果;差分盒维数和惯性矩结合了低倍组织灰度图上灰度值统计信息和空间分布信息,但前者只计算不同尺度下的灰度最大值和最小值,而舍弃了其他灰度值数据,导致损失了许多灰度值信息,可能使不同位置的差分盒维数计算结果不敏感;而后者的计算则充分利用了灰度图上所有灰度值信息,对计算结果会更加敏感。另外,相对于标准差和惯性矩,差分盒维数具有尺度独立性,受灰度图尺寸和分辨率的影响较小。从3种方法的整体变化趋势来看,铸坯不同位置的惯性矩和差分盒维数虽然变化趋势相同,但变化程度存在很大差别。可通过计算表征参数最大值与最小值之差与平均值之比作为定量比较差分盒维数与惯性矩的整体变化程度。相较于平均值来说,差分盒维数的整体变化程度很小,仅有3.6%,而惯性矩的变化程度达到了68.7%。因此,惯性矩的变化相对于差分盒维数更加敏感。标准差、差分盒维数和惯性矩均表明等轴晶区C元素分布不均匀度相较于柱状晶区更高,但相差程度存在较大差别,其中等轴晶区标准差比柱状晶区高25.38%;差分盒维数比柱状晶区高1.57%;惯性矩比柱状晶区高35.61%。

同时,通过图57和9可以看出,惯性矩与差分盒维数在5#到6#位置处的变化趋势与标准差存在差别。从5#到6#位置,标准差增大,差分盒维数和惯性矩降低。在分析过程中也发现,铸坯同一位置的碳偏析点尺寸不同,因此本工作基于偏析点的尺寸将偏析点分为I、II、III 3种类型。3种类型对应的偏析点面积分别为:> 1 mm2、0.1~1 mm2、< 0.1 mm2。根据实际铸坯的偏析程度,本研究把 I 类归为大型偏析点,II类归为中型偏析点,III 类归为小型偏析点。根据偏析点的尺寸,可由此计算不同类型偏析点面积比(不同类型偏析点面积与偏析点总面积之比)。图13为5#和6#位置处的I、II、III 3类偏析点的面积比。可以看出,I和II类偏析点面积比变化剧烈,III类偏析点面积比变化程度接近。其中,6#相对5#位置,I类碳偏析点面积比增加25.82%,II类偏析点面积比减小23.03%。5#到6#位置,标准差增大,表明标准差可以作为衡量I类偏析点面积比的指标,即标准差越大,对应区域内I类偏析点面积比越大。相反,5#到6#位置,差分盒维数和惯性矩降低。表明差分盒维数和惯性矩可以作为衡量II类偏析点面积比的指标,即差分盒维数和惯性矩越大,对应区域内II类偏析点面积比越大。另外,I类偏析点面积比越大或II类偏析点面积比越大均可能造成铸坯后续热处理均匀化的难度增加。

图13

图13   5#、6#位置I、II、III 3类偏析点的面积比

Fig.13   Area ratios of I, II, and III three kinds of segregation points (5#-6#)


综上所述,标准差、差分盒维数和惯性矩从不同角度定量表征C元素分布不均匀度,3种方法的变化趋势、变化程度以及等轴晶区与柱状晶区C元素分布不均匀度之差均存在不同差别;同时标准差可以作为衡量I类偏析点面积比的指标,差分盒维数和惯性矩可以作为衡量II类偏析点面积比的指标。因此,对于高碳钢铸坯大区域C元素分布不均匀度的全面衡量可考虑3种表征方法结合使用,从而实现对C元素分布不均匀度的综合评价。

3 结论

(1) 针对定量表征高碳钢铸坯大区域C元素分布不均匀度的不足,基于铸坯低倍组织灰度图,可利用标准差、差分盒维数和惯性矩分别从整体平均不均匀度、不同尺度下的不均匀度和局部不均匀度等角度对C元素分布不均匀度进行综合有效表达。

(2) 由连铸坯外部柱状晶区至内部等轴晶区,C元素分布不均匀度整体呈逐渐增加的趋势,且等轴晶区C元素分布不均匀度均值相对柱状晶区高20.85%。

(3) 标准差主要基于灰度值统计特性,差分盒维数和惯性矩结合了低倍组织灰度图上灰度值统计信息和空间分布信息;并且差分盒维数具有尺度独立性,受灰度图尺寸和分辨率的影响较小,而惯性矩对微区C元素分布不均匀度变化会更加敏感。此外,标准差主要受大型偏析点(> 1 mm2)影响,而差分盒维数和惯性矩主要受中型偏析点(0.1~1 mm2)影响。对于高碳钢铸坯大区域C元素分布不均匀度的全面衡量可考虑3种表征方法结合使用。

参考文献

Wang G D.

Technology innovation and development direction of iron and steel industry

[J]. Iron Steel, 2015, 50(9): 1

[本文引用: 1]

王国栋.

钢铁行业技术创新和发展方向

[J]. 钢铁, 2015, 50(9): 1

[本文引用: 1]

Song X F, He A R, Qiu Z S.

Influence of initial temperature difference on high strength strip buckling in laminar cooling

[J]. Iron Steel, 2016, 51(5): 52

[本文引用: 1]

宋小芳, 何安瑞, 邱增帅.

基于初始温差的高强钢板层冷过程前屈曲分析

[J]. 钢铁, 2016, 51(5): 52

[本文引用: 1]

Li J, Wu M H, Ludwig A, et al.

Simulation of macrosegregation in a 2.45-ton steel ingot using a three-phase mixed columnar-equiaxed model

[J]. Int. J. Heat Mass Transf., 2014, 72: 668

[本文引用: 1]

Ge H H, Ren F L, Li J, et al.

Four-phase dendritic model for the prediction of macrosegregation, shrinkage cavity, and porosity in a 55-ton ingot

[J]. Metall. Mater. Trans., 2017, 48A: 1139

[本文引用: 1]

Pickering E J.

Macrosegregation in steel ingots: The applicability of modelling and characterisation techniques

[J]. ISIJ Int., 2013, 53: 935

[本文引用: 1]

Choudhary S K, Ganguly S.

Morphology and segregation in continuously cast high carbon steel billets

[J]. ISIJ Int., 2007, 47: 1759

Mayer F, Wu M, Ludwig A.

On the formation of centreline segregation in continuous slab casting of steel due to bulging and/or feeding

[J]. Steel Res. Int., 2010, 81: 660

[本文引用: 1]

Miyamura K, Kitamura S Y, Sakaguchi S, et al.

Development of segregation etch print method and its application to investigation of CC slab segregation

[J]. Trans. Iron Steel Inst. Japan, 1984, 24: 718

[本文引用: 1]

Sang B G, Kang X H, Li D Z.

A novel technique for reducing macrosegregation in heavy steel ingots

[J]. J. Mater. Process. Technol., 2010, 210: 703

[本文引用: 1]

Lesoult G.

Macrosegregation in steel strands and ingots: Characterisation, formation and consequences

[J]. Mater. Sci. Eng., 2005, A413-414: 19

[本文引用: 1]

Yoshida N, Umezawa O, Nagai K.

Influence of phosphorus on solidification structure in continuously cast 0.1 mass% carbon steel

[J]. ISIJ Int., 2003, 43: 348

[本文引用: 1]

Wang H Z.

Original position statistic distribution analysis (original position analysis)—A new analytical method in research and quality evaluation of materials

[J]. Sci. China, 2003, 46B: 119

[本文引用: 1]

Chen Y, Yang S B, Zhu M Y.

Key technologies of internal quality control for continuously cast high-speed rail steel bloom

[J]. Iron Steel, 2006, 41(12): 36

[本文引用: 1]

陈 永, 杨素波, 朱苗勇.

高速轨用钢连铸坯内部质量控制的关键技术

[J]. 钢铁, 2006, 41(12): 36

[本文引用: 1]

Chang Y, Hou Z B, Wang W, et al.

Characteristics and mechanism of carbon element and manganese element distribution along casting direction in high-carbon billet

[J]. Iron Steel, 2016, 51(11): 43

[本文引用: 1]

常 毅, 侯自兵, 王 伟.

高碳钢方坯碳锰元素纵向分布及机理

[J]. 钢铁, 2016, 51(11): 43

[本文引用: 1]

Harada H, Miyazawa K I, Matsumiya T, et al.

Numerical prediction of solidification structure taking account of molten steel flow

[J]. Tetsu Hagané, 2003, 89: 265

[本文引用: 1]

原田寛, 宮沢憲一, 松宮徹.

流動の影響を考慮した凝固組織推定モデル

[J]. 鐵と鋼. 2003, 89: 265

[本文引用: 1]

Zhang N, Yang P, Mao W M.

Influence of columnar grains on the cold rolling texture evolution in Fe-3%Si electrical steel

[J]. Acta Metall. Sin., 2012, 48: 782

[本文引用: 1]

张 宁, 杨 平, 毛卫民.

柱状晶对Fe-3%Si电工钢冷轧织构演变规律的影响

[J]. 金属学报, 2012, 48: 782

[本文引用: 1]

Hou Z B, Guo Z A, Guo D W, et al.

A new method for carbon content distribution based on grayscale image of casting blank macrostructure in carbon steel

[J]. J. Iron Steel Res., 2019, 31: 620

[本文引用: 2]

侯自兵, 郭中傲, 郭东伟.

利用碳钢铸坯组织灰度图获取C含量分布的方法

[J]. 钢铁研究学报, 2019, 31: 620

[本文引用: 2]

Lai Y L, Chen W F, Hu Z W, et al.

Corrosion of magnesium alloy and its corrosion morphology fractal

[J]. Mater. Prot., 2018, 51(6): 28

[本文引用: 1]

赖勇来, 陈文峰, 胡志文.

镁合金腐蚀及其腐蚀形貌分形

[J]. 材料保护, 2018, 51(6): 28

[本文引用: 1]

Wang G T, Deng Y H, Wu L M.

Metallographic quantitative analysis using digital image process

[J]. Heat Treat. Met., 2006, 31(2):66

[本文引用: 1]

王桂棠, 邓耀华, 吴黎明.

基于数字图像处理技术的金相组织定量分析

[J]. 金属热处理, 2006, 31(2): 66

[本文引用: 1]

Liao K W, Lee Y T.

Detection of rust defects on steel bridge coatings via digital image recognition

[J]. Automat. Constr., 2016, 71: 294

[本文引用: 1]

Huang Z L.

The high-temperature billet crack detection research based on digital image processing

[D] Chengdu: University of Electronic Science and Technology of China, 2013

[本文引用: 1]

黄志良.

基于数字图像处理的高温钢坯裂纹检测研究

[D]. 成都: 电子科技大学, 2013

[本文引用: 1]

Gong H J.

Estimating paddy yield based on fractal and image texture analysis

[D]. Nanjing Agricultural University, 2008

[本文引用: 1]

龚红菊.

基于分形理论及图像纹理分析的水稻产量预测方法研究

[D]. 南京: 南京农业大学, 2008

[本文引用: 1]

Mandelbrot B B. The Fractal Geometry of Nature [M]. San Francisco: Freeman, 1982: 206

[本文引用: 1]

Wang Y Q, Huang L P, Li J.

The application of fractal theory in the surface analyzing of worn materials

[J]. Mater. Rep., 1998, 12(2): 8

[本文引用: 1]

王乙潜, 黄立萍, 郦 剑.

分形理论在材料磨损表面分析中的应用

[J]. 材料导报, 1998, 12(2): 8

[本文引用: 1]

Mandelbrot B B, Passoja D E, Paullay A J.

Fractal character of fracture surfaces of metals

[J]. Nature, 1984, 308: 721

[本文引用: 2]

Wang S Y, Song S Z.

Corrosion morphology diagnosing syste of metallic materials in seawater based on fractal

[J]. Acta Metall. Sin., 2004, 40: 94

[本文引用: 3]

王守琰, 宋诗哲.

基于分形的金属材料海水腐蚀形貌图像分析诊断系统

[J]. 金属学报, 2004, 40: 94

[本文引用: 3]

Sarkar N, Chaudhuri B B.

An efficient differential box-counting approach to compute fractal dimension of image

[J]. IEEE Trans. Syst. Man Cybern., 1994, 24: 115

[本文引用: 2]

Haralick R M.

Statistical and structural approaches to texture

[J]. Proc. IEEE, 1957, 67: 786

[本文引用: 1]

Renzetti F R, Zortea L.

Use of a gray level co-occurrence matrix to characterize duplex stainless steel phases microstructure

[J]. Frattura ed Integrità Strutturale, 2013, 5: 43

[本文引用: 1]

Tong X Y, Li H X, Yao L J, et al.

Feature extraction and analysis of surface microscopic image of pure copper subjecting low cycle fatigue

[J]. Mech. Sci. Technol. Aerosp. Eng., 2015, 34: 1446

[本文引用: 1]

童小燕, 李洪旭, 姚磊江.

纯铜低周疲劳表面显微形貌的特征提取与分析

[J]. 机械科学与技术, 2015, 34: 1446

[本文引用: 1]

Hou Z B, Cao J H, Chang Y, et al.

Morphology characteristics of carbon segregation in die steel billet by ESR based on fractal dimension

[J]. Acta Metall. Sin., 2017, 53: 769

[本文引用: 2]

侯自兵, 曹江海, 常 毅.

基于分形维数的模具钢电渣重熔铸坯碳偏析形貌特征研究

[J]. 金属学报, 2017, 53: 769

[本文引用: 2]

Wang H Z, Zhao P, Chen J W, et al.

In situ statistical distribution analysis of low alloy steel continuous casting slab

[J]. Sci. China, 2005, 35E: 260

[本文引用: 1]

王海舟, 赵 沛, 陈吉文.

低合金钢连铸坯的原位统计分布分析研究

[J]. 中国科学, 2005, 35E: 260

[本文引用: 1]

Bowe T F, Brody H D, Flemings M C.

Measurements of solute redistribution in dendritic solidification

[J]. Trans. Metall. Soc. AIME, 1996, 236: 624

[本文引用: 1]

Hou Z B, Guo D W, Cao J H, et al.

A method based on the centroid of segregation points: A Voronoi polygon application to solidification of alloys

[J]. J. Alloys Compd., 2018, 762: 508

Burton J A, Prim R C, Slichter W P.

The distribution of solute in crystals grown from the melt. Part I. Theoretical

[J]. J. Chem. Phys., 1953, 21: 1987

[本文引用: 2]

/