vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)
Решение задачи S12 в объёмной постановке, с применением элементов SOLID45
Рис. 4.13 КЭ-модель рамы (элементы SOLID45) с указанием закреплений и нагрузок. Жёсткие стойки показаны красным. 3D-визуализация
Рис. 4.14 1-ая форма потери устойчивости для КЭ-модели. Критическая сила Pкр = 1071,000 т
ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009 |
IV-9 |
vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)
Рис. 4.15 2-ая форма потери устойчивости для КЭ-модели. Критическая сила Pкр = 137806 т
ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009 |
IV-10 |
vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)
Задача 5. Разрушение образца кирпичной кладки (сравнение с экспериментом)
Источник |
1. Кашеварова Г.Г., Акулова А.Н. Проведение натурных |
|
|
экспериментов для исследования влияния физико- |
|
|
механических свойств материалов кирпичной кладки на |
|
|
процесс разрушения. Сборник трудов Вестник УГТУ- |
|
|
УПИ. Строительство и образование. №12 (83), – |
|
|
Екатеринбург, 2006, с.с.110-112 |
Тип задачи: |
Статический расчет с учетом физической нелинейности |
|
|
(ANTYPE = 0) |
Тип |
верифицируемых |
SOLID65 (3-D 8-узловой объемный элемент) |
КЭ: |
|
|
Входные файлы: |
kashevarova_1.mac |
Постановка задачи
Рассматривается процесс разрушения образца представительного объема многорядной кирпичной кладки (рис. 5.1) при сжатии. Вертикально по направлению оси Z к верхней поверхности образца прикладывается кинематическое воздействие в виде перемещения U0. Учитывается собственный вес.
Рис. 5.1 Образец кирпичной кладки
Определяется напряженно-деформированное состояние образца и исследуется процесс развития трещин при нагружении с учетом сгущения сетки. Результаты расчета сравниваются с результатами натурных экспериментов.
ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009 |
V-1 |
vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)
Свойства материалов |
Геометрические размеры |
Кирпич |
Длина кирпича – 0,25 м; |
Е = 7×108 Па |
ширина кирпича – 0,12 м; |
ν = 0,25 |
высота кирпича 0,065 м. |
ρ = 1900 кг/м3 |
|
σсж = 2×107 Па – предел прочности на |
|
одноосное сжатие; |
|
σр = 1,2×106 Па – предел прочности на |
|
одноосное растяжение. |
Толщина растворного шва – 0,012 м. |
Раствор |
Е = 3,5×108 Па |
|
ν = 0,2 |
|
ρ = 2000 кг/м3 |
|
σсж = 1,2×107 Па – предел прочности на |
|
одноосное сжатие; |
|
σр = 7×105 Па – предел прочности на |
|
одноосное растяжение. |
|
Натурные эксперименты
Работы выполнены сотрудниками строительного факультета ПГТУ в лаборатории механических испытаний Горного института РАН. Для испытаний использовано два пресса грузоподъемностью 200 т (УП-200) и 25 т (ZVICK).
В результате проведенных натурных экспериментов по разрушению представительного объема кирпичной кладки (3 образца) и отдельных кирпичей получены полные диаграммы деформирования. Процесс разрушения образца кирпичной кладки показан на рис. 5.2, диаграмма деформирования одного из образцов приведена на рис.5.3.
ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009 |
V-2 |
vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)
Рис.5.2 Процесс разрушения образца кирпичной кладки
ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009 |
V-3 |
vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)
500 |
|
|
|
|
|
450 |
|
|
|
|
|
400 |
|
|
|
|
|
350 |
|
|
|
|
|
300 |
|
|
|
|
|
250 |
|
|
|
|
|
200 |
|
|
|
|
|
150 |
|
|
|
|
|
100 |
|
|
|
|
|
50 |
|
|
|
|
|
0 |
|
|
|
|
|
-50 0 |
5 |
10 |
15 |
20 |
25 |
перемещения, мм
Рис.5.3 Диаграмма деформирования образца
Ниспадающая ветвь диаграммы соответствует реальному поведению кирпичных строений, способных воспринимать внешнюю нагрузку при наличии систем трещин и разрывов. Наличие отдельных трещин – это еще не разрушение конструкции. Накапливаясь, трещины могут привести к потере несущей способности здания. Однако этот момент не является внезапным, а является результатом накопления повреждений на разных структурных уровнях.
Методика расчетных исследований
Разбиение представительного объема кладки на конечные элементы (SOLID65) производилось так, чтобы границы элементов попадали на границы разделов раствор – кирпич; в каждом конечном элементе упругие свойства однородны и изотропны и соответствуют свойствам кирпича или раствора.
Число степеней свободы – 6588 |
Число степеней свободы – 96540 |
Рис.5.4 Конечно-элементные модели
ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009 |
V-4 |
vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)
Экспериментально показано, что первые трещины в кирпичной кладке появляются при нагрузке, составляющей 40-60% от разрушающей, что приводит к нелинейному характеру деформирования материала. Нелинейный расчет выполнялся методом Ньютона-Рафсона с использованием полной процедуры NROPT,FULL, в которой матрица жесткости модифицируется на каждом шаге и процедуры NROPT, AUTO, в которой программа сама выбирает необходимую опцию метода. Кроме того, деление шагов решения на подшаги производилось разными способами (использовались команды
NSUBST и DELTIM).
Картина развития трещин для разных сеток и разных процедур решения несколько различалась, но моменты появления первой трещины и разрушения образца (разрушения 50% объема – подшаг 9999) оказались достаточно близкими.
Результаты расчетов
Ниже приведены результаты нелинейного расчета образца кирпичной кладки при кинематическом воздействии U0 для отдельных подшагов решения. Для верификации фиксировались перемещения USUM в характерных точках диаграммы деформирования:
1– до появления первых трещин (рис.5.8);
2– в момент появления первых трещин (рис.5.11);
3– в момент, соответствующий критической точке на диаграмме деформирования
(рис.5.14);
4– в момент разрушения образца (прекращение счета при разрушении 50% объема) – рис.5.17.
На рис. 5–7, 9–10 приведены также компоненты тензора напряжений до появления трещин и после появления первых трещин. Показаны картины трещин и разрушенных конечных элементов для разных подшагов решения.
1 |
|
|
|
|
|
|
NODAL SOLUTION |
|
|
|
MAY |
6 2009 |
STEP=2 |
|
|
|
|
|
|
|
|
|
20:01:10 |
SUB =3 |
|
|
|
|
|
|
TIME=1.087 |
|
|
|
|
|
SX |
(AVG) |
|
|
|
|
|
RSYS=0 |
|
|
|
|
|
|
DMX =.002209 |
|
|
|
|
|
SMN =-.246E+07 |
|
|
|
|
|
SMX =611461 |
|
|
|
|
|
|
|
MX |
|
|
|
|
|
Z |
|
|
|
|
|
|
Y |
|
|
|
|
|
X |
|
|
|
|
|
|
|
MN |
|
|
|
|
-.246E+07 |
-.178E+07 |
-.110E+07 |
-412586 |
270112 |
611461 |
|
-.212E+07 |
-.144E+07 |
|
-753935 |
-71237 |
|
Рис.5.5 Изополя напряжений SX до появления первых трещин |
|
ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009 |
V-5 |