CC BY
10
ИИ
УДК 624.012
Ф.Х. Ахметзянов - кандидат технических наук, доцент
О.В. Радайкин - аспирант
Казанский государственный архитектурно-строительный университет (КазГАСУ)
К ОПРЕДЕЛЕНИЮ ЗАВИСИМОСТЕЙ ПАРАМЕТРОВ ПОВЕРХНОСТНЫХ ТРЕЩИН В БЕТОННЫХ И ЖЕЛЕЗОБЕТОННЫХ ОБЪЁМНЫХ ЭЛЕМЕНТАХ
АННОТАЦИЯ
Для определения зависимостей параметров поверхностных несквозных трещин рассмотрена задача о растяжении объёмных бетонных элементов (в виде призм). Одна группа элементов содержала нормальную к растягивающим усилиям трещину, а другая - наклонную ( задача моделировалась методом конечных элементов в программном комплексе Ansys). Определены коэффициенты интенсивности напряжений трёх типов, для чего написан макрос (подпрограмма) на языке APDL, вычисленные коэффициенты 2-го и 3-го типа (Kp Кп) в рассматриваемых задачах использованы впервые. На основе этих коэффициентов установлены соотношения между длиной, шириной раскрытия, заглублением трещин, найдены границы зоны предразрушения, дана оценка остаточной несущей способности элементов с трещиной. Показано, что область разгрузки напряжений вокруг трещины имеет форму ромба. Полученные зависимости и процедура получения численных значений K и К составляют научную новизну.
F.H. Ahmetzjanov - candidate of technical sciences, associate professor
O.V. Radaikin - post-graduate student
Kazan State University of Architecture and Engineering (KSUAE)
TO DEFINITION OF DEPENDENCE BETWEEN SURFACE CRACKS’ PARAMETERS IN SOLID CONCRETE AND REINFORCED-CONCRETE ELEMENTS
ABSTRACT
To determine the dependencies of surface non-cutting cracks the problem of the tension of concrete solid elements (in the form of prisms) was considered. One group of elements contained the crack normal to tensile stress, while another contained the inclined crack (the task was simulated by finite elements method in Ansys software). The coefficients of three types of stress intensity factors were defined; for that the macro (Sub) in APDL language was written; the calculated coefficients of the 2nd and 3d type factors (Kt and Kn) in these tasks were used for the first time. Based on these factors, the dependencies between the length, width disclosure and depth of cracks were established, the borders of process zone were found, and the residual bearing capacity of the elements with cracks was evaluated. It was shown that the discharge stress area around the cracks was in the form of diamond. The received dependencies and the procedure for obtaining of the numerical values KI and KII constitute the scientific novelty.
Ранее нами получены зависимости между длиной, заглублением и шириной раскрытия нормальной трещины в стеновой панели. Цель этой работы - найти те же соотношения для объёмной задачи более общего вида (на примере призматических стержней с одиночной несквозной трещиной).
Для достижения цели решаются следующие задачи:
- проводится моделирование напряженно-деформированного состояния объёмных бетонных элементов с одиночной полуэллиптической трещиной (как с нормальной, так и с наклонной) методом конечных элементов;
- вычисляются коэффициенты интенсивности напряжений для трещин с различным соотношением длины и заглубления, определяется ширина раскрытия трещины;
- определяются границы зоны предразрушения в прифронтовой области макротрещины и зоны релаксации напряжений.
Таким образом, рассматривается растяжение объёмных бетонных элементов ( в виде призм 1x1x1 м), имеющих по одной поверхностной полуэллиптической трещине. Одна группа элементов содержит нормальные к растягивающим усилиям трещины, а другая - наклонные (рис. 1).
Методика исследований строится на тех же допущениях, что и для ранее рассмотренной в наших исследованиях стеновой панели, т.е. справедливой остается линейная механика квазихрупкого разрушения для бетона как однородного упругого материала.
ТЕОРИЯ ИНЖЕНЕРНЫХ СООРУЖЕНИЙ И СТРОИТЕЛЬНЫЕ КОНСТРУКЦИИ
Рис.1. Расчетная модель объёмного бетонного элемента с боковой трещиной: а) с трещиной нормального отрыва, б) с наклонной трещиной с углом наклона 300, в) локальная система координат в центре полуэллиптической трещины, г) локальная цилиндрическая система координат в произвольной точке фронта трещины
Для нахождения зависимостей между параметрами несквозной трещины и определения коэффициентов интенсивности напряжений используется метод конечных элементов (МКЭ), реализованный в программном комплексе Ansys.
В работе рассматриваются только макротрещины (сопоставимые с сечением образцов). Так, длина их составляет lcrc = 2а = 0,3м - const, заглубление hcrc = b = 0,1...0,15м (а,b - полуоси эллипса). Отметим, что длина трещины и её заглубление (соответственно большая и малая полуоси эллипса) связаны между собой
^ ~2
2
,2
х~ У 4х у -л. т-г
уравнением — + -— = 1 или —- +------------ = 1, где ху - координаты произвольной точки фронта трещины. Причем
2 Ь2 ¡к
а ь І
Сгс с/с
одним из недостатков предложенной геометрической модели трещины является ограниченность возможного
2а І І
соотношения — = —— = 1,33...3,33 (в работе принимается = 1,5...3). Поэтому МКЭ комплекса Апєує
Ь ^сгс ^сгс
применим только для призматических или массивных элементов (колонны, плотины) и плохо подходит для расчета
2а л
объёмных балок-стенок и плит. Кроме того, необходимо соблюсти условие <1, где Ж - ширина образца.
Пространство в прифронтовой зоне трещины, где, как известно, возникает особенность в полях перемещений и напряжений, разбивается специальными изопараметрическими 20-узловыми конечными элементами Solid186, что предложено в руководстве пользования Ansys (рис. 2).
В качестве примера материал образцов принят из тяжелого бетон В27,5 с характеристиками Rbt ser = 1,02МПа,
Eb = 30200МПа, KIC = 0,405МПа -л[м . Модуль упругости материала при моделировании не варьируется, так как зона концентрации напряжений (псевдопластических деформаций) невелика и нелинейностью модуля в ней можно пренебречь, а в остальной части сечения Eb — const.
Рис.2. Этап твердотельного моделирования полуэллиптической трещины в объёме (слева - разбиение бетонного элемента на участки, справа - конечно-элементное разбиение объема вокруг фронта трещины)
Для вышеприведенных двух типов бетонных элементов (рис. 1 а, б) численно определяются КИН нормального отрыва Кр поперечного Кп и продольного Кш сдвига (для нормальной трещины только К:) при различном значении
полярного угла 00 £ ф £ 1800 (см. рис. 3, 4 а-в) и для различного заглубления трещины 0,1 £ Исгс £ 0,15м. Для
этого используется тот же аппроксимирующий алгоритм, что и в плоской задаче [1], только в более общем виде для формул [2]:
К =
К„ =
К ш =
G Av
л-—- lim^,
1 + к г®о д/ r
I— G Au
12л -—- lim —¡=-,
1 + к г®о д/ r
I— G Awl
7 2 л -—- lim —<=-
1 + к г®о д/r
G
Kj = л[2ж -—- lim ((/ + Bjr);
1 + к г ®о
G
/2л —к lim (Ап + Впг);
1 + к г®о
Кн =
KIII=
(1)
G
12л , , lim ІАш + Вшг )-
1 + к г®о
£ _______________ £ _________________ £
КІ = ^2л-— Аі; Кп =,^к-— Ап; Кш = ^2л-— Аш (2)
1 + к 1 + к 1 + к
где Аы , Ау , Ам - смещение одной поверхности трещины относительно другой в соответствующих
координатных осях (рис. 3г), причем Ау - отвечает расклиниванию берегов трещины, Аы - их относительному
поперечному сдвигу, а Ам - их антиплоскому продольному сдвигу; АІ; АІІ; АІІІ, ВІ; Ви; Вш -постоянные коэффициенты для линейных функций, аппроксимирующих перемещения в узлах конечных
Е
элементов на берегах трещины; О = —,-------г - модуль сдвига; коэффициент к = 3 — 4т для плоско-
2(1+т)
деформированного состояния, т - коэффициент Пуассона.
Отметим, что К1 нормальной трещины определяется по существующей аналитической зависимости [3]:
„ сыРЬ
К, =--------------р
I 7-1 5
Е,,
(3)
где Р8 - аппроксимирующая функция; а - номинальные растягивающие напряжения; Ек - полный эллиптический интеграл второго рода; 1 - толщина образца.
На рис. 3 приведено сопоставление численного и аналитического решения К: в образце первого типа при
напряжении с = 1МПа и соотношениях = 2;2,4;2,73;3 . Расхождение результатов не превышает 5%, что
И
сгс
является удовлетворительным.
По рисунку 3 можно сделать следующие выводы:
- при небольшом заглублении (Ь < а) К: максимален, когда ф = 900 - трещина с ростом нагрузки будет развиваться за счет заглубления центральной части фронта;
1
1,2
- при Ь » а К11 00 = К19=900 - трещина может как заглубляться, так и удлиняться за счет продвижения её
кончиков по поверхности бетона;
- при Ь > а К1 максимален при ф = 00 - трещина будет удлиняться;
- причем при Ь > — а К
1
_1,2
Очевидно, что нижние две кривые графика, построенные по формуле (3), имеют более одного экстремума из-за погрешности аппроксимации функции Б В то же время экспериментальные точки по МКЭ описываются плавной кривой с одним экстремумом.
Рис.3. График сопоставления аналитического [3] и численного определения К1 (кривые линии на графике соответствуют
2а I 2а = 2 2а = 3
формуле (3)); □ - численное решение при — = _й^ = 1,5, ° - то Же при = 2 , А - при = 3
Ь Ь Ь
Поскольку в процессе своего развития трещина может удлиняться и заглубляться только одновременно, то
1 1
b » 12a — const (lcrc » 2,4Acrc). Варианты соотношений b Ф a возможны только при искусственно
созданной трещине, либо когда вдоль фронта трещины имеется некоторая неоднородность жесткостных свойств материала (за счет включений крупного заполнителя или арматуры).
Таким образом, при s = \МПа и b = 1/1,2a Kj
j=900,b=
0,46ffV м . Тогда критические напряжения
KI
= 1,309МПа (ОНС). Полученная ширина раскрытия нормальной трещины равна
1 j=900-b=
асгс = 0,0423сМПа• мм, тогда предельное раскрытие асгс = асгсссг = 0,0554мм лежит в интервале опытных значений
0,001...0,06 мм [4].
Вышеописанный численный метод применяется в работе и для наклонной трещины. Для чего, помимо К (есть аналитическая формула [3]), численно определяются Кп и Кш (рис. 4 а-в) - данные по их вычислению в литературе не приводятся. Этот этап моделирования и процедура нахождения Кп и Кш составляют научную новизну работы и практическую её ценность.
С целью автоматизации построения топологической схемы модели, разбиения на конечные элементы области вокруг фронта трещины, вычисления КИН написан макрос на встроенном в АпБуБ языке АРБЬ.
в)
Кш[<р\МПа
гУ~ N
// \
4\ V \
/ />
// \v \
А Ч- К Ь>а
Ь=1 2*~ Ь<а
0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180
Рис.4. Параметры механики разрушения для наклонной трещины
a
a
Чтобы определить соотношение —, при котором трещина будет удлиняться и заглубляться одновременно,
ь
учитываются все три коэффициента через величину интенсивности высвобождения упругой энергии деформирования (рис. 4 г):
G = G, + G,, + G,,, =
(4)
Рассуждая так же, как и относительно трещины нормального отрыва, установлено, что для наклонной трещины
2a 1
соотношение — непостоянно: для верхней половинки трещины (0 £ j £ 900) b » — a — const (lcrc » 2,4hcrc ),
b 1,2
а для нижней (90 < j £ 180°) - b » a — const (lcrc » 2hcrc). В этих рассуждениях вместо К рассматривается обобщённая характеристика G. Для одновременного роста трещины по всему фронту обеспечивается условие
G (j)» G
j-9<f ъ-а — const (пунктир на рис. 6 г), что возможно только при непостоянстве соотношения '22a .
Таким образом, при с = 1МПа и ь = ^ 2 а при 0 < р < 900, Ь = а при 90< р < 1801, получены следующие
1
значения КИН и интенсивности высвобождения энергии (табл.).
Значения КИН и С в зависимости от полярного угла ф
Таблица
j KI / s, л/м к,II / s , л/м Кш/s, л[м G / s, м
0С С,352 С,18С С,С6С
9 О с С,316 С,СС4 С,16С 4,2 • 10-6
180С 0,353 С, 182 0,061
По рисункам 4 а-в и данным таблицы можно видеть, что в области кончиков фронта трещины в процесс её распространения большой вклад вносит механизм поперечного сдвига (Кп - максимален). Напротив, в середине фронта трещины Кп - минимален, а Кш - максимален, т.е. здесь большую роль играет механизм продольного сдвига (из плоскости трещины).
Интенсивность высвобождения энергии в момент перед разрушением равна Осс = G
Kr
= 4,83 10-6 МПа • м
I j=0c
Критические напряжения (ОНС) s cr = s
кг\
= 1,15МПа, что меньше, чем для нормальной трещины такой
же длины. Ширина раскрытия наклонной трещины равна асгс = 0,0Ъ50сМПа мм, тогда предельное раскрытие а*гс = асгсСсг = 0,0403мм, что также меньше, чем это для нормальной трещины - ассс = 0,0554 мм.
Таким образом, получены следующие зависимости: = 300 = 5415 - нормальная трещина,
ассс 0,0554
=-------= 7444 - наклонная трещина. Эти соотношения согласуются с экспериментальными данными
а 0,0403
работы [5]: acrc =(1600...10200) llcrc и в среднем acrc = 319С 4crc.
Рис. 5. Изополя главных растягивающих напряжений (а - трещины нормального отрыва, б - для наклонной трещины с углом наклона 300).
Из рассмотрения изополей напряжений (рис. 5) вытекает, что вдоль фронта трещины возникает объемное растяжение. Зона предразрушения вдоль фронта по форме сходна с подобной областью для плоской задачи [1], и её размер намного меньше длины трещины и размеров образца. Это является основным условием применимости линейной механики разрушения [6] для рассматриваемой задачи.
Размер зоны предразрушения в вершине нормальной трещины для тяжелого бетона В27,5 равен
К2
й = —т----- ---ут = 12,4мм < / = 300мм, й << Ж = 1000мм . Для наклонной трещины эта зона в образце из
зрЦдб^ТГ
бетона одного класса прочности примерно на 18% меньше, чем для нормальной, и определяется формулой [6]
О Е
й = —^ = 10,6мм . Также отметим, что область разгрузки напряжений (рис. 5) имеет форму не
3р (и6^,.)
окружности, как это принято в классических работах Гриффитца, а форму ромба с большей диагональю длиной 2/
Выводы
Рассмотрена задача о растяжении объёмных бетонных элементов (в виде призм 1x1x1 м), имеющих по одной поверхностной полуэллиптической трещине. Одна группа элементов содержала нормальные к растягивающим усилиям трещины, а другая - наклонные (задача моделировалась методом конечных элементов в программном комплексе Апбуб). Определены коэффициенты интенсивности напряжений трёх типов, для чего написан макрос (подпрограмма) на языке АРБЬ, вычисленные коэффициенты 2-го и 3-го типа в рассматриваемых задачах использованы впервые. Установлены соотношения между длиной, шириной раскрытия, заглублением трещин:
/сгс » 2,4Исгс , ~~ = 5415 (ас,с = 0,0554) - для нормальной трещины; 1сгс ®(2...2,4)йсге , = 7444
с,с ас,с
(асгс = 0,0403) - для наклонной трещины. Установлено, что зона предразрушения в вершине наклонной трещины
й = 10,6мм на 18% меньше, чем для нормальной й = 12,4мм. Показано, что область разгрузки напряжений вокруг трещины имеет форму ромба.
Полученные соотношения, размер зоны предразрушения и ширина раскрытия трещин подтверждается опытными данными [4,5] и имеют научную новизну. Также научную новизну представляет собой процедура численного определения коэффициентов К , К для наклонной трещины.
Литература
1. Радайкин О.В. К оценке остаточной несущей способности стеновых железобетонных стеновых панелей методами
механики разрушения // Материалы 58-ой Республ. науч. конф. - Казань: КГАСУ, 2006. - С. 86-92.
2. Седов Л.И. Механика сплошной среды. Т. 1, Т. 2. - М.: Наука, 1970. - 547 с., 586 с.
3. Мураками Ю. Справочник по коэффициентам интенсивности напряжений. Т. 1, Т. 2. - М.: Изд-во «МИР», 1990.
- 448 с., 566 с.
4. Почтовик Г.Я., Липник В.Г., Филонидов А.М. Дефектоскопия бетона ультразвуком в энергетическом строительстве. - М.: Энергия, 1977. - 121 с.
5. Ахметзянов Ф.Х., Арсентьев Е.З., Мухаметзянов И. А. К закономерностям развития макротрещин в наружных стенах из монолитного керамзитобетона в жилых домах при эксплуатации // Тезисы докладов 3-й Всеросс. конференции «Новое в архитектуре, проектировании строительных конструкций и реконструкции НАСКР-2001», Ч. I. - Чебоксары: Изд-во Чуваш. ун-та, 2001. - С. 160-165.
6. Зайцев Ю.В. Механика разрушения для строителей. - М.: Высшая школа, 1991. - 288 с.
