книги / Мониторинг состояния цементобетонных дорожных конструкций
..pdfния целесообразно положить многопараметрическую оценку, считая на первом этапе все свойства одинаково значимыми. Так как целью работы является изучение воздействия климата на бетон, то возрастные изменения свойств, присущие бетонам в нормальных условиях, целесообразно исключить из рассмотрения.
Из возможных параметров, характеризующих отдельные свойства или состояние структуры бетона в количественном или качественном отношении, оставлены только количественные, методика определения которых достаточно отработана, что дает возможность применить упрощенные методы оценки и последующего выделения значимых параметров методами математической статистики и теории вероятности.
В соответствии с нормативными документами для производства строительных материалов и строительства в качестве таких параметров приняты: кубиковая (R) и призменная прочности (Rb), модуль упругости
(Eb) и коэффициент (γ) упругости, нижний (Rcrc0 ) и верхний (Rcrcγ ) преде-
лы микротрещинообразования, водопоглощение (W). Таким образом, считается достаточным математически свойства бетона задавать семимерным вектором Z = (R, Rb, Eb, γ, Rcrc0 , Rcrcγ , W)Т.
Этот вектор может быть представлен в виде суммы:
Z = Z0 + ∆Zt (t)+ Z*,
где Z0 – вектор начального состояния; ∆Zt(t) – вектор, характеризующий изменения от времени; Z* – вектор, определяемый влиянием воздействием условий среды. Для выделения части Z* изменения названных параметров, которые определяются воздействием условий среды, возрастные изменения значений этих параметров в нормальных условиях исключены из рассмотрения. Исследование изменений большого числа параметров требует больших объемов работ по испытаниям. Значительное их сокращение возможно, если учесть взаимосвязи этих изменений. Например, возможны довольно тесные взаимосвязи пары параметров, и тогда можно ограничиться исследованием изменений только одного из них, а изменения другого определять расчетно.
Нужно иметь в виду, что взаимосвязи не являются жесткими, детерминированными, они проявляются при статистической обработке. Принято рассматривать парные взаимосвязи, математически они проявляются в их корреляции и количественно характеризуются коэффициентами корреляции. Близость этих коэффициентов к +1 или к –1 свидетельствует
151
о тесной связи: чем дальше от этих значений, тем меньше взаимная зависимость. Дальнейшие исследования направлены на получение и анализ матрицы коэффициентов корреляции.
Далее, вместо абсолютных значений характеристик свойств бетона, испытывающего воздействие климата (R, Rb, Eb, γ, Rcrc0 , Rcrcγ , W), будем
рассматривать относительные отклонения номинальных значений этих параметров от соответствующих характеристик бетонов нормального твердения. Их совокупность представляется в виде вектора Х = (x1, x2, x3, x4, x5, x6, x7)Т. Это позволит исключить возможные систематические ошибки, которые могли иметь место из-за несовершенства определения значений названных характеристик свойств. Расчетные данные о взаимосвязях изменений отдельных параметров бетона естественного твердения под воздействием климата в виде матрицы коэффициентов корреляции приведены в табл. 4.1–4.3.
Таблица 4.1
Парные коэффициенты корреляции параметров бетона естественного твердения
|
х1 |
х2 |
х3 |
х4 |
х5 |
х6 |
х7 |
х1 |
1,0 |
0,895 |
0,653 |
0,439 |
0,740 |
0,790 |
–0,289 |
х2 |
0,895 |
1,0 |
0,772 |
0,442 |
0,874 |
0,742 |
–0,421 |
х3 |
0,653 |
0,772 |
1,0 |
0,826 |
0,618 |
0,421 |
–0,762 |
х4 |
0,439 |
0,442 |
0,826 |
1,0 |
0,112 |
0,029 |
–0,622 |
х5 |
0,740 |
0,874 |
0,618 |
0,112 |
1,0 |
0,899 |
–0,375 |
х6 |
0,790 |
0,742 |
0,421 |
0,029 |
0,899 |
1,0 |
–0,041 |
х7 |
–0,289 |
0,421 |
–0,762 |
–0,622 |
–0,375 |
–0,041 |
1,0 |
Предварительно для вычисления коэффициентов корреляции значения параметров были нормированы и центрированы, так что значения коэффициентов корреляции определялись по формуле
|
|
|
|
|
|
n |
|
|
|
n |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n∑xi x j −∑xi ∑x j |
|
|
|
|
|
|
|
||||||||
rx x |
= |
|
|
|
i |
=1 |
|
|
|
i=1 |
j=1 |
|
|
|
|
|
|
|
. |
(4.1) |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2 |
||||
i |
j |
n |
|
|
n |
|
|
|
n |
|
|
|
n |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
n |
x2 |
x |
|
n |
|
− |
x |
|
|
|
||||||||
|
|
∑ |
− |
∑ |
|
|
x |
2 |
|
|
|
|
||||||||
|
|
|
j |
|
|
j |
|
|
∑ i |
|
|
∑ i |
|
|
|
|||||
|
|
|
j=1 |
|
|
j=1 |
|
|
|
|
i=1 |
|
|
i=1 |
|
|
|
|
|
152
Таблица 4.2 Парные коэффициенты корреляции параметров пропаренных бетонов
|
х1 |
х2 |
х3 |
х4 |
х5 |
х6 |
х7 |
х1 |
1,0 |
0,912 |
0,866 |
0,046 |
0,826 |
0,647 |
–0,938 |
х2 |
0,912 |
1,0 |
0,834 |
0,292 |
0,911 |
0,733 |
–0,746 |
х3 |
0,866 |
0,834 |
1,0 |
0,135 |
0,663 |
0,434 |
–0,866 |
х4 |
0,046 |
0,292 |
0,135 |
1,0 |
0,370 |
0,442 |
–0,359 |
х5 |
0,826 |
0,911 |
0,663 |
0,370 |
1,0 |
0,913 |
–0,577 |
х6 |
0,647 |
0,733 |
0,434 |
0,442 |
0,913 |
1,0 |
–0,040 |
х7 |
–0,938 |
–0,746 |
–0,866 |
–0,359 |
–0,577 |
–0,040 |
1,0 |
Таблица 4.3
Парные коэффициенты корреляции параметров пропаренных бетонов с добавкой
|
х1 |
х2 |
х3 |
х4 |
х5 |
х6 |
х7 |
х1 |
1,0 |
0,945 |
0,272 |
0,487 |
0,868 |
0,806 |
–0,894 |
х2 |
0,945 |
1,0 |
0,093 |
0,576 |
0,945 |
0,885 |
–0,849 |
х3 |
0,272 |
0, 093 |
1,0 |
0,027 |
0,074 |
0,053 |
–0,199 |
х4 |
0,487 |
0,576 |
0,027 |
1,0 |
0,365 |
0,216 |
–0,588 |
х5 |
0,868 |
0, 945 |
0,074 |
0,365 |
1,0 |
0,939 |
–0,754 |
х6 |
0,806 |
0,885 |
0,053 |
0,216 |
0,939 |
1,0 |
–0,623 |
х7 |
–0,894 |
–0,849 |
–0,199 |
–0,588 |
–0,754 |
–0,623 |
1,0 |
Была проведена компьютерная обработка массивов данных, полученных при испытаниях. Функциональные связи между характеристиками свойств бетонов, испытывающих воздействие климата, могут быть представлены в виде полученных нами зависимостей:
|
|
0,17 −0,176W |
* |
|
|
Eb* = |
|
Eb , |
|||
|
|
|
W |
|
|
R* |
|
0,041 |
−0,214W |
* |
|
= |
R, |
||||
|
|
|
W |
|
153
R* = 0,848
b
0* |
|
|
|
Rcrc |
−0,034 R . |
||
Rcrc0 |
|||
|
b |
Остальные аналитические зависимости между характеристиками свойств цементных бетонов, испытывающих воздействие природной климатической среды, представлены ниже:
а) для бетонов естественного твердения
R* =(0,688Rb* |
Rb −0,057)R |
Eb* =(1,43ν* |
ν+0,005)Eb |
|
|
|||||||||||||
R* =(0,75E* E−0,230)R |
Eb* =(−0,164W * |
W +0,212)Eb |
|
|||||||||||||||
R* =(1,86Rcrc0* |
Rcrc0 +0,132)R |
Rb =(3,87Rcrc0 |
Rcrc28 −2,54)Rb28 |
|
||||||||||||||
R* =(2,81Rcrcν* |
Rcrcν +0,023)R |
Rcrc0 =(0,198Rb |
Rb28 +0,748)Rcrc0 |
28 |
||||||||||||||
R* |
1,18E* E |
− |
0,232 R |
W = −1,16R |
R |
|
|
+2,01 W |
|
|||||||||
b |
=( |
b |
b |
|
) b |
( |
b |
|
b28 |
|
|
|
) 28 |
|
||||
Rb* =(2,84Rcrc0* |
Rcrc0 +0,325)Rb |
Rb =(−0,86W W28 +1,73)Rb28 |
|
|||||||||||||||
Rb* =(3,44Rcrcν* |
|
Rcrcν +0,085)Rb |
|
|
|
|
|
|
|
|
|
|
|
|
||||
б) для бетонов с добавкой пластификатора П-20 |
|
|
|
|
|
|
|
|
|
|
||||||||
Rb* =(3,5Rcrc0* Rcrc0 +0,132)Rb |
Rcrc0* =(0,939Rcrcν* |
|
|
Rcrcν −0,048)Rcrc0 |
||||||||||||||
Rb* =(3,29Rcrcν* |
Rcrcν −0,032)Rb |
Rcrc0* =(0,032−0,126W * W )Rcrc0 |
|
|||||||||||||||
Rb* =(0,054−0,530W * |
W )Rb |
Rb =(4,6Rcrc0 |
|
Rcrc0 |
28 −3,59)Rb28 |
|
||||||||||||
R* =(0,905Rb* |
Rb −0,136)R |
Rcrc0 =(0,163Rb |
Rb28 +0,818)Rcrc0 |
28 |
||||||||||||||
* |
|
0* |
|
0 |
|
|
W = −0,696R |
|
R |
28 |
+1,61 W |
|
||||||
R |
=(3,08Rcrc |
Rcrc −0,025)R |
( |
|
b |
|
|
b |
|
) 28 |
|
|||||||
R* =(2,86Rcrcν* |
Rcrcν −0,17)R |
W =(−4,56Rcrc0 |
|
Rcrc0 |
28 +5,49)W28 |
|||||||||||||
R* =(0,032−0,533W * |
W )R |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
в) для пропаренных бетонов без добавки |
|
|
|
|
|
|
|
|
|
|
|
|||||||
R* =(1,095Rb* |
Rb −0,053)R |
Rb* =(0,848Rcrc0* |
|
|
Rcrc0 −0,034)Rb |
|||||||||||||
R* =(0,97Eb* Eb −0,168)R |
Rb* =(1,98Rcrcν* |
|
Rcrcν −0,032)Rb |
|||||||||||||||
R* =(2,98Rcrc0* |
Rcrc0 +0,101)R |
Rb* =(0,031−0,142W * W )Rb |
||||||||||||||||
R* =(−0,214W * |
W +0,041)R |
Eb* =(0,17−0,176W * |
W )Eb |
|
||||||||||||||
Rb* =(0,775Eb* |
Eb −0,106)Rb |
Rcrc0* =(0,82Rcrcν* |
|
|
Rcrcν −0,06)Rcrc0 |
154
Предварительно группировка параметров произведена методом корреляционных плеяд. Оценка параметров при этом осуществлялась по сравнению с пороговым значением r = 0,72, что больше критических значений r = 0,707, определённых при надёжности P = 0,95 из выражения
rn n −1 ≥1,87 [141].
Корреляционные плеяды для нашего случая изображены на рис. 4.1.
а
б
в
Рис. 4.1. Корреляционные плеяды: а – для бетонов естественного твердения; б – для пропаренных бетонов; в – для пропаренных бетонов с добавкой ПЯ-01
На основании корреляционных плеяд была осуществлена разбивка параметров на три группы. Задача разбиения параметров на группы неявно предполагает справедливость гипотезы о том, что в основе сильной связи между параметрами лежит некоторый неизвестный внутренний
155
фактор, который и обуславливает эту связь. Из этого следует, что группировку целесообразно проводить таким образом, чтобы в одну группу попадали наиболее связанные параметры и, наоборот, параметры далёкие, непохожие попадали в разные группы.
Мера близости между параметрами принята в виде K =1+1αR , где
|
|
|
|
|
|
|
|
1 |
|
|
|
n |
|
( p) |
(q) |
|
2 |
2 |
|
R = |
|
|
|
− xi |
|
|
представлена в табл. 4.4. |
||
|
∑i=1 xi |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
Из данных, представленных в табл. 4.4, следует, что в одну группу должны попасть (х1, х2, х5); (х3, х4), прослеживается автономность параметров (х7, х6). Следовательно, нужно применять более строгий способ группировки параметров.
|
|
|
|
|
|
Таблица 4.4 |
||
Мера близости между параметрами пропаренного бетона |
|
|||||||
|
|
|
|
|
|
|
|
|
Мера |
|
|
С параметрами |
|
|
|
||
близости |
1 |
2 |
3 |
4 |
5 |
6 |
|
7 |
параметров |
|
|||||||
|
|
|
|
|
|
|
|
|
1 |
1,0 |
0,932 |
0,321 |
0,334 |
0,792 |
0,389 |
|
0,09 |
2 |
0,932 |
1,0 |
0,352 |
0,328 |
0,822 |
0,341 |
|
0,083 |
|
|
|
|
|
|
|
|
|
3 |
0,321 |
0,352 |
1,0 |
0,608 |
0,269 |
0,156 |
|
0,078 |
4 |
0,334 |
0,328 |
0,608 |
1,0 |
0,262 |
0,171 |
|
0,097 |
|
|
|
|
|
|
|
|
|
5 |
0,792 |
0,822 |
0,269 |
0,262 |
1,0 |
0,279 |
|
0,092 |
6 |
0,389 |
0,341 |
0,156 |
0,171 |
0,279 |
1,0 |
|
0,129 |
|
|
|
|
|
|
|
|
|
7 |
0,09 |
0,083 |
0,078 |
0,097 |
0,092 |
0,129 |
|
1,0 |
Такая более строгая группировка параметров произведена с помощью алгоритма агрегирования, предложенного в [87]. При таком алгоритме объединяются такие две группы параметров Ap, Aq, для которых
K(Ap, Aq) = max(Al, At) l, t = 1, 2, … n – k + 1 l ≠ t. |
(4.2) |
При этом мера близости вычислялась по формуле
K(Ap, Aq) = |
1 |
∑x |
|
|
∑x |
|
K (xi , x j ), |
(4.3) |
|
nA |
nA |
A |
|
A |
|||||
|
i |
|
p |
j |
|
p |
|
||
|
p |
q |
|
|
|
|
|
|
|
где K (xi , x j ) – мера связи между двумя параметрами.
156
K (x |
, x |
|
) = |
|
|
1 |
, |
(4.4) |
|
1 |
+αR2 |
||||||
i |
|
j |
|
|
|
|||
|
|
|
R(xi , x j ) = |
n |
||||
где α – свободный член алгоритма; |
∑(xip x jp )2 – расстояние |
i=1
между параметрами.
Реализация алгоритма агрегирования позволила предусмотреть разбивку параметров на три группы, а именно: (х1, х2); (х3, х4); (х5, х6, х7).
Так как признаки, принятые для оценки состояния бетона, зависимы по своему физическому смыслу (мерой их статистической зависимости служат коэффициенты парной корреляции), то целесообразно из первоначальной системы исходных признаков сформировать систему статистически независимых признаков.
Известно, что одним из методов такого формирования признаков является метод главных компонент. Каждый из независимых признаков является комбинацией исходных, а коэффициенты при исходных признаках представляют собой компоненты собственных векторов матрицы взаимной корреляции признаков.
Собственные числа матриц связи для всех исследуемых бетонов представлены в табл. 4.5.
Таблица 4.5 Собственные числа матриц коэффициентов парной корреляции признаков
Бетон |
|
Значения собственных чисел |
|
|||||
1 |
2 |
3 |
4 |
5 |
6 |
7 |
||
|
||||||||
|
|
|
|
|
|
|
|
|
Естественного твердения |
4,13 |
1,84 |
0,76 |
0,43 |
0,29 |
0,08 |
–0,53 |
|
Пропаренный |
4,71 |
1,76 |
0,42 |
0,13 |
0,08 |
0,02 |
–0,12 |
|
С добавкой ПЯ-01 пропаренный |
4,70 |
1,075 |
0,84 |
0,24 |
0,055 |
0,05 |
0,0 |
Мерой информативности главных компонентов принято считать отношение вида
|
n |
|
|
|
η= |
∑λi |
|
|
|
i=1 |
, |
(4.5) |
||
S p A |
||||
|
|
|
где λi – собственные числа матрицы корреляции; SpA – сумма её диагональных элементов.
157
В нашем случае при учёте двух собственных чисел
η = 4,13 +1,84 |
= 0,853, |
η |
2 |
= 4,71+1,76 |
= 0,924, |
|
1 |
7 |
|
|
7 |
|
|
|
|
|
|
|
η = 4,7 +1,075 = 3 0,825.
7
Таким образом, можно считать, что первые два главных компонента содержат от 82,5 до 92,4 % всей информации об исходной совокупности, а система из семи первичных признаков может быть заменена двумя обобщёнными. При этом информативность главных компонентов для пропаренного бетона (0,924) выше информативности компонентов бетона естественного твердения (0,853) и бетона с добавкой ПЯ-01 (0,825).
Метод главных компонентов использован для сопоставления результатов испытаний лабораторных образцов бетона, твердевших в естественных условиях, в моделированных условиях, и образцов кернов, отобранных из конструкций. Графическая интерпретация результатов представлена на рис. 4.2 в координатах U1, U2.
Рис. 4.2. Геометрическая интерпретация метода главных компонентов
Для сравнения были взяты образцы следующих бетонов: естественного твердения, пропаренного, пропаренного с добавкой пластификатора ПЯ-01, бетон кернов, отобранных из плит ПДН, образцов, претерпевших 70 циклов воздействия, смоделированных в лабораторных условиях. Идентификационные признаки названных бетонов представлены в табл. 4.6.
158
|
|
|
|
|
|
|
|
Таблица 4.6 |
|
Идентификационные признаки бетонов |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Бетоны |
|
|
|
|
|
кернов, |
естественного |
|
пропаренные из состава |
||||
|
|
|
|
|
1 и помещённые |
|||
|
отобранных |
|
|
|
|
|||
xi |
твердения |
|
1 |
|
2 |
|
в лабораторные |
|
из плит ПДН |
|
|
|
|||||
|
|
|
|
|
|
|
|
условия |
|
|
|
Номер системы |
|
||||
|
1 |
2 |
|
3 |
|
4 |
|
5 |
x1 |
–0,24 |
–0,175 |
|
–0,208 |
|
–0,259 |
|
–0,396 |
x2 |
–0,225 |
–0,284 |
|
–0,21 |
|
–0,176 |
|
–0,28 |
x3 |
–0,215 |
–0,156 |
|
–0,133 |
|
–0,111 |
|
–0,11 |
x4 |
–0,085 |
–0,048 |
|
–0,007 |
|
0,01 |
|
–,015 |
x5 |
–0,08 |
–0,222 |
|
–0,133 |
|
–0,111 |
|
–0,23 |
x6 |
–0,25 |
–0,065 |
|
–0,006 |
|
–0,074 |
|
–0,11 |
x7 |
1,9 |
2,348 |
|
1,0 |
|
0,146 |
|
2,52 |
За меру близости между системами может быть принято: среднеквадратическое расстояние
S= 1 ∑n d 2 (xi( p) , xi(q) ),
n i=1
евклидово расстояние
n 2 d = ∑(xi( p) − xi(q) )
i=1
(4.6)
(4.7)
или косинус угла между векторами в пространстве признаков
|
n |
|
|
|
|
|
cosϕ = |
∑ xi( p) xi(q) |
, |
(4.8) |
|||
i=1 |
|
|
||||
|
|
|
|
|
|
|
|
∑(xi( p) )2 |
1 |
∑(xi(q) )2 |
1 |
|
|
|
2 |
2 |
|
|||
|
n |
|
n |
|
|
|
|
i=1 |
|
i=1 |
|
|
где хi – координата состояния в n-мерном пространстве признаков; p, q – соответственно индексы принадлежности к состоянию.
Расстояние между системами в пространстве признаков, оцениваемое по показателю «сos φ» по (4.8) и для сравнения – в виде евклидового расстояния «d» по (4.7) при просчитывании по всем параметрам, представлено в табл. 4.7.
159
|
|
|
|
|
|
|
Таблица 4.7 |
|
|
|
Расстояние между системами |
|
|
|
|||
|
|
|
|
|
|
|
|
|
Расстояние |
1–2 |
|
1–3 |
1–5 |
2–3 |
2–5 |
|
3–5 |
Cos φ |
0,986 |
|
0,986 |
0,992 |
0,985 |
0,995 |
|
0,99 |
d |
0,515 |
|
0,925 |
0,684 |
1,353 |
0,306 |
|
1,552 |
|
|
|
|
|
|
|
|
|
Из анализа данных, представленных в табл. 4.6 и 4.7, можно сделать вывод о том, что характер зависимостей изменений свойств бетона образцов или конструкций в реальной климатической среде от свойств бетона носит «похожий» характер, и при дальнейшем исследовании бетон изделий может быть заменён бетоном образцов в реальных или смоделированных климатических условиях.
При этом изменение под воздействием климата свойств бетона естественного твердения и пластифицированного добавкой ПЯ-01 существенно отличается от соответствующих изменений свойств пропаренного бетона. Это пример того, как состояние структуры бетона до начала воздействия климата оказывает влияние на характер изменения состояния под воздействием климата.
Вышеприведенный способ способствует группировке исходных параметров, но не исключает необходимости определения всей их совокупности. При использовании метода главных компонентов количество исходных параметров не уменьшается, так как они полностью входят в состав главных компонентов.
В то же время группировка параметров позволила разбить совокупность исходных параметров на группы, несущие информацию одного вида и характеризующие бетон с одной из сторон (прочностные, деформационные и группа параметров, характеризующих структуру). Сильная корреляционная связь между параметрами одной группы может привести в случае учёта всех параметров при оценке состояния бетона к неоднократному учёту информации одного вида, при этом значимость отдельных факторов будет неоправданно завышаться.
При оценивании качества бетона оперирование большим числом параметров увеличивает риск ошибок и объём вычислительной работы и, как правило, не обеспечивает точность оценки состояния. Среди исходных параметров, очевидно, имеются малоинформативные, малозначимые в решении конкретной задачи, которые целесообразно исключить из рассмотрения. Но для такого исключения требуются серьезные обоснования.
160