橡胶基混凝土铺装材料细观分析与数值模拟.docVIP

橡胶基混凝土铺装材料细观分析与数值模拟 刘泽佳，汤立群，黄小清，邱燕红 （，广州 51064）[1-7]。圆形与骨料实际形状相差较大，而随机多边形也都没有考虑边数和角度等骨料形状的影响。 近年来图像分析技术在混合料内部微细观结构研究方面得到了广泛的应用[8-11]。本文的研究对象新型钢桥面铺装材料是一种橡胶混凝土，从视觉角度来看，可利用色彩描述它的混合料体积组分。本文将基于图像分析技术，研究骨料颗粒的几何特性等的统计规律，并据此生成随机骨料模型对试件剖面内的粗骨料及水泥砂浆基体直接进行有限元网格剖分，然后根据骨料在网格中的位置判定单元类型，并依单元类型赋予相应的材料特性，对这种材料的性质进行数值模拟。目前对边界层的研究还不够透彻，因此只将混凝土及其它铺装材料看作一级二相材料。 2新型桥面铺装材料的结构几何描述 钢筋混凝土以及沥青混合料等铺装材料图像的灰度直方图通常具有明显的双峰性。运用数字图像处理技术从灰度图像中提取骨料颗粒，将灰度图像转化为二值图像，利用图像处理专业软件对粒子进行轮廓跟踪、边界搜索、边缘检测、颗粒识别及描述，从而获得颗粒质心、面积、周长等。 图1 原始图像（灰度图） 图2去除小颗粒后的边界曲线 图1为橡胶混凝土断面的原始图像。图2是运用专业软件对原始图像进行开运算处理后，成功去除图像的颗粒噪声，且减小了骨料边缘细小的“突刺” 等影响获取的单像素宽的骨料边界轮廓。 3 不规则二维图形的多边形拟合 常用的边界表达方法有链码法、边界段法、多边形法和标记法。由于混凝土多采用碎石骨料，其截面的二维图形为不规则多边形，多边形法是用一系列线段的封闭集合来逼近真实曲线，以直代曲，以多边形来拟合目标图像不规则的边界曲线，与其他三种方法比较，多边形法在保证骨料形状失真很小的前提下，更准确、易于理解而且数据量较少。得到这个拟合多边形后，再统计多边形的边数，边长，面积等，进而统计描述目标图像形状的特征值。目前常用的多边形拟合技术是迭代端点拟合法，这种方法在编程调试时比较困难，而且运算量大，处理速度慢。 本文针对混凝土骨料形状随机性和不确定性的特点，提出一种简单高效的多边形拟合快速算法。 通过运用包络矩形寻找不规则图形四条相切边的方法，再对四条相切边之间的弧段进行拟合，算法步骤及原理如下： (1)确定每个骨料的包络矩形MNOP； (2)确定骨料边界与其包络矩形相切的轮廓线AB,CD,EF,GH； (3)分别对四条边之间的弧段BC,DE,FG, HA进行拟合。将每个弧段均分为n份， n的取值具体根据不规则图形的形状而定，通过BI和BJ的夹角α判断该边界点是否为所求多边形的其中一个拐点，若αβ，表示B、I、J近似在一条线，继续将BI和BK的夹角与β进行比较；若αβ，表示I是所求多边形的一个拐点，后面以此类推。 步骤(1)-(3)循环可以得到第一次拟合的多边形。 (4)对第一次搜到的边界点进行第二次多边形拟合。同样通过角度判断拐点。在进行第二次多边形拟合的角度判断时，可以适当增大标准值β。 图3 多边形拟合快速算法示意图 图4 多边形拟合后的橡胶混凝土图形 表1 橡胶混凝土的骨料边数、面积对照表 Table 1 the correlation table of the aggregate side number and area of rubber concrete 编号 边数 第一次拟合 边数 第二次拟合 原始骨料面积（pixel） 拟合多边形骨料面积（pixel） 误差 （％） 1 9 8 1535.00 1543.500000 0.55 2 11 6 2860.00 2750.500000 3.83 3 11 7 1588.00 1676.500000 5.57 4 8 5 6179.00 6650.500000 7.62 5 10 7 4141.00 4127.000000 0.34 6 7 4 2274.00 2330.000000 2.46 7 8 5 1654.00 1710.000000 3.39 8 9 5 1967.00 1807.000000 8.13 9 7 6 1342.00 1419.000000 5.74 10 10 7 1678.00 1727.500000 2.95 对骨料进行多边形拟合后，能获取断面所有骨料的形态学指标，图4是对图2中骨料的边界二值图进行两次多边形拟合后的结果图。表1是用本文的多边形拟合方法得到的结果与专业软件的处理结果对比。二次拟合后的多边形骨料边数比第一次拟合结果减少了，从而减少了数据量。对原始图像的骨料进行多边形拟合后，拟合多边形的骨料面积与真实面积相差不超过10%