Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ANSYS Mechanical

.pdf
Скачиваний:
914
Добавлен:
31.08.2019
Размер:
61.81 Mб
Скачать

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 (“исследовательские” задачи)

нагрузка, кН

-5

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

ANSYS Mechanical. Верификационный отчет. Том 3 (“исследовательские

” задачи)

 

 

vk.com/club152685050 | vk.com/id446425943

 

 

 

1

 

 

 

 

 

 

 

NODAL SOLUTION

 

 

 

 

MAY

6 2009

STEP=2

 

 

 

 

 

 

 

 

 

 

 

20:01:40

SUB =3

 

 

 

 

 

 

 

TIME=1.087

 

 

 

 

 

 

SY

(AVG)

 

 

 

 

 

 

RSYS=0

 

 

 

 

 

 

 

DMX =.002209

 

 

 

 

 

 

SMN =-.246E+07

 

 

 

 

 

 

SMX =609099

 

 

 

 

 

 

 

 

MX

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

MN Y

 

 

 

 

 

 

 

X

 

 

 

 

 

 

-.246E+07

-.178E+07

-.110E+07

 

-414022

268058

609099

 

 

-.212E+07

-.144E+07

-755062

-72982

 

Рис.5.6 Изополя напряжений SY до появления первых трещин

 

1

 

 

 

 

 

 

 

NODAL SOLUTION

 

 

 

 

MAY

6 2009

STEP=2

 

 

 

 

 

 

 

 

 

 

 

20:02:09

SUB =3

 

 

 

 

 

 

 

TIME=1.087

 

 

 

 

 

 

SZ

(AVG)

 

 

 

 

 

 

RSYS=0

 

 

 

 

 

 

 

DMX =.002209

 

 

 

 

 

 

SMN =-.835E+07

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

MX Y

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

MN

 

 

 

 

-.835E+07

-.649E+07

-.557E+07

-.464E+07

-.278E+07

-927769

0

 

 

-.742E+07

-.371E+07

 

-.186E+07

 

Рис.5.7 Изополя напряжений SZ до появления первых трещин

 

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

V-6

ANSYS Mechanical. Верификационный отчет. Том 3 (“исследовательские

” задачи)

 

 

vk.com/club152685050 | vk.com/id446425943

 

 

 

1

 

 

 

 

 

 

 

NODAL SOLUTION

 

 

 

 

MAY

6 2009

STEP=2

 

 

 

 

 

20:02:41

SUB =3

 

 

 

 

 

 

 

TIME=1.087

MX

 

 

 

 

 

USUM

(AVG)

 

 

 

 

 

RSYS=0

 

 

 

 

 

 

 

DMX =.002209

 

 

 

 

 

 

SMX =.002209

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

Y

 

 

 

 

 

 

 

X

MN

 

 

 

 

 

 

 

 

 

 

 

0

 

.491E-03

.982E-03

 

.001472

.001963

.002209

 

 

.245E-03

.736E-03

.001227

.001718

Рис.5.8 Изополя перемещений USUM (м) до появления первых трещин

 

1

 

 

 

 

 

 

NODAL SOLUTION

 

 

 

MAY

6 2009

STEP=2

 

 

 

 

 

19:57:44

SUB =4

 

 

 

 

 

 

TIME=1.144

 

 

 

 

 

SX

(AVG)

 

 

 

 

 

RSYS=0

 

 

 

 

 

 

DMX =.003629

 

 

 

 

 

SMN =-.408E+07

 

MX

 

 

 

SMX =550074

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

Y

 

 

 

 

 

 

X

 

 

 

 

 

 

 

MN

 

 

 

 

-.408E+07

-.305E+07

-.202E+07

-992304

35948

550074

 

 

-.356E+07

-.253E+07

-.151E+07

-478178

 

Рис. 5.9 Изополя напряжений SX при появлении первых трещин

 

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

V-7

ANSYS Mechanical. Верификационный отчет. Том 3 (“исследовательские

” задачи)

 

 

vk.com/club152685050 | vk.com/id446425943

 

 

 

1

 

 

 

 

 

 

 

NODAL SOLUTION

 

 

 

 

MAY

6 2009

STEP=2

 

 

 

 

 

 

 

 

 

 

 

19:56:33

SUB =4

 

 

 

 

 

 

 

TIME=1.144

 

 

 

 

 

 

SZ

(AVG)

 

 

 

 

 

 

RSYS=0

 

 

 

 

 

 

 

DMX =.003629

 

 

 

 

 

 

SMN =-.138E+08

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

MX Y

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

MN

 

 

 

 

 

-.138E+08

-.107E+08

-.765E+07

 

-.459E+07

-.153E+07

 

 

-.122E+08

-.918E+07

-.612E+07

-.306E+07

0

 

Рис.5.10 Изополя напряжений SZ при появлении первых трещин

 

1

 

 

 

 

 

 

NODAL SOLUTION

 

 

 

MAY

6 2009

STEP=2

 

 

 

 

19:55:47

SUB =4

 

 

 

 

 

 

TIME=1.144

MX

 

 

 

 

USUM

(AVG)

 

 

 

 

RSYS=0

 

 

 

 

 

 

DMX =.003629

 

 

 

 

 

SMX =.003629

 

 

 

 

 

 

Z

 

 

 

 

 

 

Y

 

 

 

 

 

X

MN

 

 

 

 

 

 

 

 

 

 

0

.806E-03

.001613

.002419

.003225

.003629

 

.403E-03

.00121

.002016

.002822

 

Рис.5.11 Изополя перемещений USUM (м) при появлении первых трещин

 

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

V-8