一种混凝土群裂纹演化的多重分形量化表征方法与流程

1.本发明涉及混凝土结构健康检测与安全预警技术领域，具体为一种混凝土群裂纹演化的多重分形量化表征方法。


背景技术：

2.混凝土结构受材料内在缺陷以及外部环境的影响，容易产生裂缝。裂缝促使混凝土结构碳化加剧，强度降低，甚至造成灾害性事故。混凝土结构裂缝的产生和扩展常被用于度量混凝土结构的生命周期。因此，研究混凝土结构的裂缝演化规律，对混凝土结构性态监控、性能评价和安全预警具有重要的理论意义和工程应用价值。混凝土结构裂缝演化，特别是多因素耦合作用下混凝土结构群裂缝演化，伴随着与外界物质和能量的交换，其本质是一个开放、耗散、具有混沌性质的非线性行为。现有多重分形技术能够有效刻画混凝土群裂缝演化过程中瞬时裂缝分布的非线性和复杂性，但是难以量化表征混凝土结构裂缝渐进扩展过程。


技术实现要素：

3.本发明的目的在于提供一种混凝土群裂纹演化的多重分形量化表征方法，用于量化表征混凝土群裂缝渐进扩展过程。
4.本发明提供了一种混凝土群裂纹演化的多重分形量化表征方法，包括：
5.采集混凝土结构上的目标区域的混凝土群裂缝分布状态序列图像；
6.构建混凝土群裂缝分布状态序列图像中的每个图像的多重分形奇异谱；
7.根据多重分形奇异谱的曲线特征判断裂缝分布的多重分形特征；
8.根据判断结果将混凝土群裂缝分布状态序列图像的多个多重分形奇异谱集结于同一个坐标系构成混凝土群裂缝状态序列图像的多重分形奇异谱簇；
9.从所述多重分形奇异谱簇中提取混凝土群裂缝分布状态序列图像的多参数的分形特征因子集合。
10.优选的，所述采集混凝土结构上的目标区域的群裂缝分布状态序列图像，包括以下步骤：
11.对混凝土结构连续作用多个荷载步，并分别记录所述混凝土结构上的目标区域在每个荷载步作用下的混凝土群裂缝分布状态图像。
12.优选的，所述构建混凝土群裂缝分布状态序列图像中的每个图像的多重分形奇异谱，包括以下步骤：
13.对每个荷载作用下的混凝土群裂缝分布状态图像分别进行二值化处理，得到每个荷载作用下的混凝土群裂缝分布二值图像；
14.分别计算每个所述混凝土群裂缝分布二值图像的奇异性指数和奇异性维数；
15.依据每个所述混凝土群裂缝分布二值图像的奇异性指数和奇异性维数分别绘制每个所述混凝土群裂缝分布二值图像的多重分形奇异谱。
16.优选的，所述计算混凝土群裂缝分布二值图像的奇异性指数和奇异性维数，包括以下步骤：
17.计算所述混凝土群裂缝分布二值图像的总像素m；
18.分别采用不同尺寸盒子覆盖所述混凝土群裂缝分布二值图像，并分别计算当每个不同尺寸盒子的尺寸为l时，所述混凝土群裂缝分布二值图像中的裂缝的在盒子中的局部裂缝比例pi(l)，具体表达公式为：
[0019][0020]
其中mi(l)表示盒子尺寸为l时，第i个盒子内所含裂缝的像素；
[0021]
对每个不同尺寸盒子的局部裂缝比例pi(l)进行归一化处理，具体计算公式如下：
[0022][0023]
其中q为权重因子，其取值一般为[-10,10]；
[0024]
μi(q,l)表示所述混凝土群裂缝分布二值图像局部像素差异性，其取值范围为[0,1]；
[0025]
计算在不同权重因子下的所述混凝土群裂缝分布二值图像的奇异性指数，具体计算公式如下：
[0026][0027]
其中a(q)为在不同权重因子下的所述混凝土群裂缝分布二值图像的奇异性指数；
[0028]
计算在不同权重因子下的所述混凝土群裂缝分布二值图像的奇异性维数，具体计算如下：
[0029][0030]
其中f(q)为在不同权重因子下的所述混凝土群裂缝分布二值图像的奇异性维数。
[0031]
优选的，以在不同权重因子下的所述混凝土群裂缝分布二值图像的奇异性指数a(q)为横坐标，并以f(q)在不同权重因子下的所述混凝土群裂缝分布二值图像的奇异性维数为纵坐标绘制所述混凝土群裂缝分布二值图像的多重分形奇异谱；
[0032]
若当权重因子q》0时，所述混凝土群裂缝分布二值图像的多重分形奇异谱的曲线特征为f(q)随着a(q)的增大而增大，且当权重因子q《0时，所述混凝土群裂缝分布二值图像的多重分形奇异谱的曲线特征为f(q)随着a(q)的增大而减小，则所述混凝土结构上的群裂缝分布具有多重分形特性。
[0033]
优选的，所述混凝土群裂缝分布状态序列图像的多参数的分形特征因子集合，包括：
[0034]
最小奇异性指数a
min
，其表征所述混凝土结构上的目标区域的最大局部裂缝比例影响力；
[0035]
最大奇异性指数a
max
，其表征所述混凝土结构上的目标区域的最小局部裂缝比例影响力；
[0036]
多重分形奇异谱谱宽
△
a,其表征所述混凝土结构上的目标区域的裂缝分布不均匀程度,其表达式为：
[0037]
△
a＝a
max-a
min
ꢀꢀꢀ
(5)
[0038]
最大奇异性维数f
max
，其表征所述混凝土结构上的目标区域的裂缝分布的整体情况；
[0039]
多重分形奇异谱奇异性维数差
△
f，其表征所述混凝土结构上的目标区域的不同局部裂缝比例种类差异所导致的分布不均匀性,其表达式为：
[0040]
△
f＝f(a
max
)-f(a
min
)
ꢀꢀꢀ
(6)
[0041]
其中f(a
max
)为当所述混凝土群裂缝分布二值图像的多重分形奇异谱中的异性指数为a
max
时的奇异性维数值；
[0042]
f(a
min
)为当所述混凝土群裂缝分布二值图像的多重分形奇异谱中的异性指数为a
min
时的奇异性维数值。
[0043]
与现有技术相比，本发明的有益效果：
[0044]
本发明特别适用于混凝土群裂缝渐进扩展过程的量化表征，突破了传统裂缝统计方法难以刻画荷载作用下混凝土群裂纹动态演化的局限性，超越了现有多重分形用于描述混凝土群裂缝分布的技术水平，提出的多参数分形特征因子能够定量表征荷载作用下混凝土群裂缝动态演化，形成混凝土结构性态监控、性能评价和安全预警的新技术。
附图说明
[0045]
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：
[0046]
图1是本发明提出的一种混凝土群裂纹演化的多重分形量化表征方法的原理流程图；
[0047]
图2是本发明的实施例中提供的静态循环荷载作用下钢筋混凝土剪力墙裂缝分布；
[0048]
图3是本发明的实施例中提供的钢筋混凝土剪力墙在每个荷载步作用下的混凝土群裂缝分布二值图像；
[0049]
图4是本发明的实施例中提供的混凝土群裂缝分布二值图像上的某一群裂缝的多参数的分形特征因子；
[0050]
图5是本发明的实施例中提供的5个连续加载步下的裂缝分布的多重分形奇异谱簇。
具体实施方式
[0051]
下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，但应当理解本发明的保护范围并不受具体实施方式的限制。
[0052]
实施例
[0053]
如图1-5所示，一种混凝土群裂纹演化的多重分形量化表征方法，包括：
[0054]
步骤1：采集混凝土结构上的目标区域的混凝土群裂缝分布状态序列图像，包括以下步骤：对混凝土结构连续作用多个荷载步，并分别记录混凝土结构上的目标区域在每个
荷载步作用下的混凝土群裂缝分布状态图像，构成混凝土群裂缝分布状态序列图像。
[0055]
本实施例中以某钢筋混凝土剪力墙在静态循环荷载作用下的表观裂缝分布为研究对象，实际分析5个连续荷载步下裂缝演化过程中多重分形特征。裂缝分布的原始图像是在作用荷载的过程中，通过在墙体上标记，如图2所示，然后按照一定比例采用cad绘制而成，以jpg的格式导出裂缝分布图像。
[0056]
步骤2：构建混凝土群裂缝分布状态序列图像中的每个图像的多重分形奇异谱，具体包括以下步骤：
[0057]
步骤2.1：去除混凝土群裂缝分布状态序列图像中的非裂缝信息，对每个荷载作用下的混凝土群裂缝分布状态图像分别进行二值化处理，得到每个荷载作用下的混凝土群裂缝分布二值图像。由于本实施例中以jpg的格式呈现的裂缝分布图像，已经不包含非裂缝信息，且裂缝分布清晰；只需采用matlab对裂缝图像进行二值化处理，处理后的图像如图3所示。
[0058]
步骤2.2：分别计算每个混凝土群裂缝分布二值图像的奇异性指数和奇异性维数。
[0059]
计算混凝土群裂缝分布二值图像的奇异性指数和奇异性维数，具体包括以下步骤：
[0060]
步骤2.2.1：计算混凝土群裂缝分布二值图像的总像素m。
[0061]
步骤2.2.2：分别采用不同尺寸盒子覆盖混凝土群裂缝分布二值图像，并分别计算当每个不同尺寸盒子的尺寸为l时，混凝土群裂缝分布二值图像中的裂缝的在盒子中的局部裂缝比例pi(l)，具体表达公式为：
[0062][0063]
其中mi(l)表示盒子尺寸为l时，第i个盒子内所含裂缝的像素。
[0064]
步骤2.2.3：对每个不同尺寸盒子的局部裂缝比例pi(r)进行归一化处理，具体计算公式如下：
[0065][0066]
其中q为权重因子，其取值一般为[-10,10]；
[0067]
μi(q,l)表示混凝土群裂缝分布二值图像局部像素差异性，其取值范围为[0,1]；
[0068]
通过引入权重因子q，用于凸显局部像素的差异性：
[0069]
当q》0时，局部裂缝比例较大占主导作用；
[0070]
当q《0时，局部裂缝比例较小者占主导作用；
[0071]
而当q＝0，局部裂缝比例的大小不加以区别。
[0072]
其中，对于每一个权重因子q，以logl为横坐标，为纵坐标构成的点集，采用最小二乘法拟合该点集所得直线斜率即为a(q)；以为坐标构成的点集，采用最小二乘法拟合点集所得直线斜率即为f(q)。
[0073]
步骤2.2.4：计算在不同权重因子下的所述混凝土群裂缝分布二值图像的奇异性指数，具体计算公式如下：
[0074][0075]
其中a(q)为在不同权重因子下的混凝土群裂缝分布二值图像的奇异性指数。
[0076]
步骤2.2.5：计算在不同权重因子下的所述混凝土群裂缝分布二值图像的奇异性维数，具体计算如下：
[0077][0078]
其中f(q)为在不同权重因子下的混凝土群裂缝分布二值图像的奇异性维数。
[0079]
步骤2.3：依据每个混凝土群裂缝分布二值图像的奇异性指数和奇异性维数分别绘制每个混凝土群裂缝分布二值图像的多重分形奇异谱。
[0080]
以在不同权重因子下的混凝土群裂缝分布二值图像的奇异性指数a(q)为横坐标，并以f(q)在不同权重因子下的混凝土群裂缝分布二值图像的奇异性维数为纵坐标绘制混凝土结构上的群裂缝分布二值图像的多重分形奇异谱。
[0081]
步骤3：根据多重分形奇异谱的曲线特征判断裂缝分布的多重分形特征。
[0082]
根据混凝土群裂缝分布二值图像的多重分形奇异谱的曲线特征判断混凝土结构上的裂缝分布的多重分形特征，具体如下：
[0083]
若当权重因子q》0时，混凝土群裂缝分布二值图像的多重分形奇异谱的曲线特征为f(q)随着a(q)的增大而增大，且当权重因子q《0时，混凝土群裂缝分布二值图像的多重分形奇异谱的曲线特征为f(q)随着a(q)的增大而减小，则混凝土结构上的群裂缝分布具有多重分形特性。每一个混凝土群裂缝演化状态均对应着一个多重分形奇异谱。
[0084]
如本实施例中某一群裂缝的多参数的分形特征因子的曲线呈现为开口向下的凸型曲线，且开口具有一定宽度，如图4所示，f(q)与a(q)的单调性满足：q》0部分单调递增，q《0部分单调递减，则认为混凝土结构上的裂缝分布具有多重分形特性。
[0085]
步骤4：根据判断结果将混凝土群裂缝分布状态序列图像的多个多重分形奇异谱集结于同一个坐标系构成混凝土群裂缝状态序列图像的多重分形奇异谱簇。
[0086]
将每个混凝土群裂缝分布二值图像的多重分形奇异谱集结在同一个坐标系中。本实施例中将5个连续荷载步下的混凝土结构上的裂缝分布的多重奇异谱集结于同一个坐标系构成混凝土群裂缝状态序列的多重分形奇异谱簇图如图5所示。
[0087]
步骤5：从多重分形奇异谱簇中提取混凝土群裂缝分布状态序列图像的多参数的分形特征因子集合。
[0088]
当混凝土结构上的裂缝分布具有多重分形特征的裂缝分布，方可继续提取多参数的分形特征因子集合。提取混凝土群裂缝分布状态序列图像的多参数的分形特征因子集合，该分形特征因子集合表征了在荷载作用下混凝土结构上的群裂纹的演化过程。
[0089]
混凝土群裂缝分布状态序列图像的多参数的分形特征因子集合包括：
[0090]
最小奇异性指数a
min
，其表征混凝土结构上的目标区域的最大局部裂缝比例影响力；
[0091]
最大奇异性指数a
max
，其表征混凝土结构上的目标区域的最小局部裂缝比例影响力；
[0092]
多重分形奇异谱谱宽
△
a,其表征混凝土结构上的目标区域的裂缝分布不均匀程度,其表达式为：
[0093]
△
a＝a
max-a
min
ꢀꢀꢀ
(5)
[0094]
最大奇异性维数f
max
，其表征混凝土结构上的目标区域的裂缝分布的整体情况；
[0095]
多重分形奇异谱奇异性维数差
△
f，其表征混凝土结构上的目标区域的不同局部裂缝比例种类差异所导致的分布不均匀性,其表达式为：
[0096]
△
f＝f(a
max
)-f(a
min
)
ꢀꢀꢀ
(6)
[0097]
其中f(a
max
)为当混凝土群裂缝分布二值图像的多重分形奇异谱中的异性指数为a
max
时的奇异性维数值；
[0098]
f(a
min
)为当混凝土群裂缝分布二值图像的多重分形奇异谱中的异性指数为a
min
时的奇异性维数值。
[0099]
表1 5个连续荷载步下群裂缝分布的多参数分形特征因子集合
[0100][0101]
在本实施例中，分别提取5个连续荷载步下混凝土群裂缝状态序列多参数的分形特征因子集合，如上表1所示。
[0102]
最后说明的是：以上公开的仅为本发明的一个具体实施例，但是，本发明实施例并非局限于此，任何本领域的技术人员能思之的变化都应落入本发明的保护范围。