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

ANSYS Mechanical

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

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 4 (опыт применения в проектной и экспертной практике)

Рис. 9. Эквивалентные напряжения (МПа) для М=9.44 тсм

Полученные результаты численного моделирования подтверждают и конкретизируют несовершенство проектного и реализованного вариантов верхнего узла колонн, “рабочим” состоянием которого является либо “докритическое” состояние, с восприятием и передачей на колонны со связями значимого момента до 10-15 тсм), либо «закритическое», с разрушением по сварке или по «телу» блюмса, сопровождаемым «мгновенным» падением момента и динамическим нагружением системы.

1.6 Нижний узел колонн

Была составлена также объемная КЭ-модель нижнего узла колонны, состоящая из 6000 узлов (для ½ части с учетом симметрии задачи). Использовались КЭ второго порядка с промежуточными узлами, а также контактные элементы. В качестве основных нагрузок были заданы сдвигающие усилия, которые по результатам расчетов достигли 400кН и выше для эксплуатационного состояния, и угол поворота.

1.7 Устойчивость колонны со связями

Разработаны “оболочечные” модели колонны с двумя коробами с использованием восьмиузловых изопараметрических КЭ второго порядка, поддерживающих сдвиг. Верхняя и нижняя зоны колонны (блюмс с ребрами) моделировались жесткими стержневыми КЭ. Внизу были заданы условия шарнирного опирания, вверху и по границе коробов–расчетные нагрузки и (или) перемещения. Задача решалась в физически (пластичность металла колонн) и геометрически (значимые перемещения) нелинейной постановке. Результаты расчетных исследований, представленные на рис. 10, свидетельствуют, что потеря устойчивости сжатой колонны с растянутыми связями происходит не по “традиционной” стержневой форме при высокой критической нагрузке, а по сугубо оболочечной форме в зоне коробов (с изменением кольцевой формы сечения на сложную эллипсообразную c большим радиальным прогибом со стороны внешней стороны колонны) при весьма высоком уровне пластических деформаций при нагрузках, близких к заданным (от 0,96 до 1,14 для различных подвариантов).

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

I-9

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 4 (опыт применения в проектной и экспертной практике)

Рис. 10. “Оболочечное” деформирование и потеря устойчивости колонны со связями СОК:

(a) – перемещения наружной, внутренней и верхней точек колонны, (б) – эквивалентные пластические деформации, εmax =19,8%

1.8 Опорный ж/б контур оболочки

Выполненное моделирование ребристой оболочки толщиной 70-100 мм пластинчатооболочечными и стержневыми КЭ, по нашему мнению, является обоснованным. Тем не менее, в условиях достаточно высоких сжимающих напряжений на большей части площади необходим анализ прочности по растрескиванию в критических зонах стыка ребер и скорлупы на объемных моделях. Как показали проведенные тестовые расчеты, применение объемных элементов бетона даже при одном слое элементов по толщине возможно без потери точности.

При 10%-й нагрузке (от полной нагрузки распалубки) выявлено резкое изменение жесткости системы, вызванное появлением большой зоны трещинообразования. Как показал анализ результатов расчетов, НДС бетонной части в нелинейной постановке существенно зависит от жесткости колонны со связями. Важным для общей жесткости системы является изгиб относительно слабо армированной зоны переменной толщины (от 200 мм к 500 мм). При 12-15%-й нагрузке произошло резкое уменьшение жесткости системы, вызванное появлением большой зоны трещинообразования в опорном контуре (рис. 12).

Рис.12. Опорный контур: (a) – эквивалентные напряжения в арматуре и закладной пластине со стержнями (σmax=69 МПа). (б) – Растягивающие напряжения в оголовке (МПа)

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

I-10

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 4 (опыт применения в проектной и экспертной практике)

1.9 Выводы

Результатом проведенных экспертных исследований является подтвержденная возможность осуществления двух групп сценариев разрушения сооружения СОК “Трансваальпарк”, связанных с ошибками проектирования:

1.В соответствии с первым вариантом “спусковым крючком” явилось подготовленное всей статической и циклической историей нагружения разрушение одного из тяжело нагруженных узлов опорной системы или потеря устойчивости одной из колонн (вследствие больших пластических деформаций в зоне соединения с распорками), которая через цепь событий привела к образованию, распространению и раскрытие субрадиальной трещины в оболочке. В завершающей стадии падающая оболочка «потянула» за собой и всю опорную систему, не рассчитанную на восприятие таких нагрузок.

2.По второму возможному сценарию вначале произошла потеря устойчивости гибкой ж/б ребристой оболочки (вследствие больших перемещений, вызванных, в том числе, проявлением эффектов ползучести бетона при действии весовой и снеговой нагрузок), сопровождаемая образованием локальных и магистральной трещин. После раскрытия указанной и, возможно, кольцевой магистральной трещины падающая оболочка “потянула” за собой и всю опорную систему.

Архитектурное, инженерно-техническое и административно-бюрократическое сообщества должны усвоить уроки подобных аварий-обрушений. Один из них: эффективное, не сдерживающее архитектурное творчество решение задач обеспечения надежности и безопасности уникальных сооружений возможно силами квалифицированных специалистоврасчетчиков на основе математического моделирования с привлечением современных численных методов и программных комплексов. Эти математические модели должны сопровождать объекты на всех этапах их жизненного цикла (проектирования, строительства, эксплуатации, ремонта и реконструкции), обеспечивая адекватный анализ и прогноз состояния в составе информационно-диагностических систем мониторинга.

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

I-11

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 4 (опыт применения в проектной и экспертной практике)

2. Конечноэлементное моделировании поведения и сценариев обрушения системы “грунтовое основание–фундаментные конструкции–стальные колонны со связями

вантовая оболочка покрытия” Басманного рынка (2006-2007 гг.)

Врамках строительно-технической экспертизы был выполнен анализ причин обрушения большепролетного здания Басманного рынка (23 февраля 2006г.).

Здесь представлен визуальный “дайджест” результатов численного моделирования поведения и сценариев обрушения системы “грунтовое основание–фундаментные конструк- ции–стальные колонны со связями – вантовая оболочка покрытия” и критических узлов, выполненного с применением ПК ANSYS, Лира и СТАДИО.

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

II-1

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 4 (опыт применения в проектной и экспертной практике)

Этапы возведения и эксплуатации (СТАДИО, Лира)

Трехмерные модели

 

фрагментов-узлов (АNSYS)

Сценарий локального разрушения и прогрессирующего обрушения (СТАДИО)

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

II-2

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 4 (опыт применения в проектной и экспертной практике)

3. Анализ причин обрушения промздания. Трехмерное численное моделирование поведения системы “колонны–фермы–связи” и критических зон

3.1. При проведении экспертизы причин “мгновенного” обрушения всех стропильных ферм со связями и кровлей (один из цехов завода), произошедшего 22.12.2005г. на заключительном этапе строительства, использовался методологический подход, ранее выработанный и систематизированный в ходе экспертиз причин обрушения зданий СОК “Трансвааль-парк” [1–3] и Басманного рынка. Этапы экспертизы:

сбор, анализ и обобщение научной, проектной и исполнительной документации по системе “здание–основание”, данных после разрушения;

выработка гипотез и возможных сценариев локальных разрушений и прогрессирующего обрушения конструкций покрытия;

разработка и верификация альтернативных расчетных многопараметрических моделей системы “здание–основание” и критических узлов;

проверка расчетов НДС, прочности и устойчивости конструкций цеха, выполненных Проектантом, на соответствие проектным предположениям, реализованным конструктивным решениям и требованиям действовавших СНиП’ов [4, 5];

выполнение комплексных расчетных исследований НДС, прочности и устойчивости несущих конструкций цеха при проектных, фактических, гипотетических эксплуатационных

и“катастрофических” сочетаниях нагрузок с учетом данных натурных исследований;

итоговый анализ и обобщение в экспертном заключении результатов математического моделирования и натурных исследований.

3.2. Расчетные исследования НДС и прочности выполнены для 3-х выделенных групп нагрузок на покрытие, каждая из которых несет свою смысловую экспертную “нагрузку” (таблица 1,4):

проектные (проверка проектных решений),

фактические, с учетом выявленных отступлений от проекта, на момент обрушения (поиск истинных причин обрушения),

гипотетические при так и не “случившейся” эксплуатации (прогноз последствий,

риски).

Рассмотрены 2 варианта проектных нагрузок: на момент обрушения и при эксплуатации (принятые в проекте).

Фактические нагрузки на момент обрушения оценивались по собранным и проанализированным проектным, натурным (замеры) и экспертным данным. С учетом выявленного существенного “разброса” в качестве расчетных принимались минимальные, средние и максимально возможные значения.

Гипотетические нагрузки при эксплуатации также принимались в диапазоне, границы которого характеризуют минимально и максимально возможные значения с учетом установленных отступлений от проекта.

Снеговые и ветровые нагрузки за весь период строительства цеха принимались по данным Росгидромета. Момент обрушения характеризовался максимальными снеговыми отложениями (до 32 см без учета “мешков”) и незначительным ветром.

Были инициированы натурные испытания, в частности, измерения вибраций от проходящего внутризаводского железнодорожного транспорта, не показавшие скольконибудь значимых амплитуд виброперемещений.

3.3. Моделирование поведения конструкций покрытия цеха при проектных, фактических и гипотетических эксплуатационных нагрузках выполнялось по ПК Robot Millennium [6]. Построены и верифицированы подробные пространственно-стержневые конечноэлементные модели, реалистично описывающие геометрико-жесткостные и нагрузочные характеристики и НДС системы “колонны–фермы–связи” (рис.1).

Некоторые результаты многопараметрических расчетных исследований представлены в таблица 2 и на рис.2.

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

III-1

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 4 (опыт применения в проектной и экспертной практике)

При “средних” и, тем более, “максимальных” фактических (22.12.2005г.) нагрузках выявлен опасный уровень усилий в нижнем растянутом поясе, в раскосах Р2, Р3 и стойках С2 стропильных ферм ФС1 (таблица 2). Моделирование поведение системы при локальном разрушении наиболее напряженных раскосов Р2 (таблица 3, рис.3,4) подтверждает прогрессирующий характер дальнейшего обрушения конструкций покрытия.

3.4. Большое и, как показал опыт экспертизы СОК “Трансвааль-парк” и здания Басманного рынка, оправданное внимание уделено моделированию реального поведения наиболее напряженных и/или “нетрадиционно” решенных конструктивных узлов. Фокус нашего интереса – сварные соединения “пояс П1 – накладка Н1 – раскосы Р” наиболее нагруженной стропильной фермы ФС1.

Построены и исследованы последовательно уточняющиеся объемные физически (пластичность металла), геометрически (большие перемещения) и структурно (контакты размыканием-смыканием и трением) нелинейные конечноэлементные модели (ПК ANSYS [7]). Результаты расчетов подтвердили опасения относительно прочности самого узла и его влияния на работу системы в целом при фактических нагрузках

(рис.4-7).

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

III-2

vk.com/club152685050 | vk.com/id446425943

ANSYS Mechanical. Верификационный отчет. Том 4 (опыт применения в проектной и экспертной практике)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1

 

 

 

 

 

 

 

 

Фактическая нагрузка на покрытие цеха в момент обрушения (22.12.2005г.)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

№№

 

 

Наименование элемента покрытия

Нормативная

Коэффициент

Фактическая нагрузка на покрытие, кг/м2

п/п

 

 

нагрузка

надёжности

“Минимальная”

 

“Средняя”

“Максимальная”

 

 

 

 

 

 

 

1

 

Собственный вес металлоконструкций

80

1,05

84

 

84

84

2

 

Профлист Н78-750-0,8

 

12

1,05

13

 

13

13

3*

 

Бетон

в

профнастиле

Нприв=8,92см (при

 

 

 

 

 

 

 

Н=75+60=135мм)

 

0,0892×2186=195кг/м2

205

1,1

225,5

 

 

 

 

 

 

 

арматура в теле бетона 14281кг:1440=9,9кг/м2

 

 

 

 

 

 

3**

 

Бетон

в

профнастиле

Нприв=9,375см (при

205

1,1

 

 

225,5

225,5

 

Н=75+60=135мм)

0,09375×2186=205кг/м2

 

 

 

 

 

 

 

 

 

 

арматура в теле бетона 14281кг:1440=9,9кг/м2

9,9

1,1

 

 

10,89

10,89

4

 

Пароизоляция – плёнка

 

0,4

1,3

0,52

 

0,52

 

5

**

 

Утеплитель

– пенополистерол ПСБс-35 –

2,655

1,3

3,45

 

3,45

3,45

 

 

 

150мм. 0,15×17,7кг/м3=2,655 кг/м2

 

6**

 

1 слой гидроизоляционного ковра

3,5

1,3

4,55

 

4,55

4.55

7*

 

Песок разуклонки

0,091×1570=143кг/м2

143

1,3

186

 

 

 

7**

 

0,119×1800=214кг/м2

214,2

1,3

 

 

278,5

 

7***

 

 

 

 

0,119×1980=235,6кг/м2

235,6

1,3

 

 

 

306,3

8

*

 

Асфальтовая стяжка

– 40мм

82

1,3

106,6

 

 

 

 

 

0,04×2050 кг/м3=82кг/м2

 

 

 

 

8

*

 

Асфальтовая стяжка – 43,4мм

88,97

1,3

 

 

115,66

115,66

 

 

0,0434×2050 кг/м3=89кг/м2

 

 

 

 

 

 

 

Итого, постоянная нагрузка

 

 

622,8

 

736,07

736,07

 

 

 

 

 

Временная нагрузка

 

 

 

 

 

 

Снег толщиной 32см (отсутствия снегового мешка)

 

 

80,6

 

80,6

 

Снег толщиной 32см (со снеговым мешком)

 

 

 

 

 

120

 

 

 

 

 

 

Всего

 

 

 

704,1

 

816,67

884,07

 

 

 

Отношение к проектной (с технологической

 

 

0,998

 

1,161

1,254

 

 

 

нагрузкой) постоянная+снег=525+180=705,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Примечания:

 

 

 

 

 

 

 

 

 

1.При подсчёте нагрузки по п.п. 3**, 5**, 6** - приняты данные по результатам натурных замеров;

2.При подсчёте нагрузки по п. 7** приняты данные Проектировщика.

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

III-3

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 4 (опыт применения в проектной и экспертной практике)

Рис. 1. Пространственно-стержневая КЭ-Robot-модель несущих металлоконструкций цеха. Общий вид

Рис. 2. КЭ-Robot-модель несущих металлоконструкций цеха. Перемещения общие от нормативной фактической нагрузки

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

III-4

vk.com/club152685050ANSYS Mechanical. Верификационный| vk.com/id446425943отчет. Том 4 (опыт применения в проектной и экспертной практике)

Таблица 2

Отношения запаса прочности сечений фермы ФС1 (конструкция по проекту, нагрузки “средние” фактические, на 16,1% больше проектных)

Стержень

Cечение

Материал

Lay

Laz

Отношение

Нагружение

680

Tг 200x12

С255

19.78

19.78

0.03

7 Расчетная

684

Tг 200x12

С255

19.78

19.78

0.03

7 Расчетная

661

Tг 200x12

С255

19.78

19.78

0.04

7 Расчетная

662

Tг 200x12

С255

19.78

19.78

0.04

7 Расчетная

679

Tг 200x12

С255

19.78

19.78

0.04

7 Расчетная

902

Tг 120x6

С255

68.22

68.22

0.73

7 Расчетная

923

Tг 100x6

С255

81.82

81.82

0.73

7 Расчетная

929

Tг 100x6

С255

81.82

81.82

0.73

7 Расчетная

951

Tг 100x6

С255

83.56

83.56

0.73

7 Расчетная

799

Tг 180x9

С255

43.13

43.13

0.74

7 Расчетная

850

Tг 180x9

С255

43.13

43.13

0.74

7 Расчетная

882

Tг 120x6

С255

66.80

66.80

0.74

7 Расчетная

908

Tг 120x6

С255

68.22

68.22

0.74

7 Расчетная

921

Tг 100x6

С255

81.82

81.82

0.74

7 Расчетная

879

Tг 120x6

С255

66.80

66.80

0.75

7 Расчетная

883

Tг 120x6

С255

66.80

66.80

0.75

7 Расчетная

890

Tг 120x6

С255

66.80

66.80

0.75

7 Расчетная

903

Tг 120x6

С255

68.22

68.22

0.75

7 Расчетная

910

Tг 120x6

С255

68.22

68.22

0.75

7 Расчетная

788

Tг 180x9

С255

43.13

43.13

0.76

7 Расчетная

811

Tг 180x9

С255

43.13

43.13

0.76

7 Расчетная

886

Tг 120x6

С255

66.80

66.80

0.76

7 Расчетная

907

Tг 120x6

С255

68.22

68.22

0.76

7 Расчетная

909

Tг 120x6

С255

68.22

68.22

0.76

7 Расчетная

917

Tг 120x6

С255

68.22

68.22

0.76

7 Расчетная

914

Tг 120x6

С255

68.22

68.22

0.77

7 Расчетная

915

Tг 120x6

С255

68.22

68.22

0.77

7 Расчетная

826

Tг 180x9

С255

43.13

43.13

0.78

7 Расчетная

911

Tг 120x6

С255

68.22

68.22

0.78

7 Расчетная

816

Tг 180x9 НП1

С255

43.13

43.13

0.80

7 Расчетная

817

Tг 180x9

С255

43.13

43.13

0.80

7 Расчетная

820

Tг 180x9

С255

43.13

43.13

0.80

7 Расчетная

878

Tг 120x6

С255

66.80

66.80

0.80

7 Расчетная

887

Tг 120x6

С255

66.80

66.80

0.80

7 Расчетная

896

Tг 120x6

С255

66.80

66.80

0.80

7 Расчетная

892

Tг 120x6

С255

66.80

66.80

0.81

7 Расчетная

898

Tг 120x6

С255

68.22

68.22

0.81

7 Расчетная

792

Tг 180x9

С255

43.13

43.13

0.83

7 Расчетная

880

Tг 120x6

С255

66.80

66.80

0.83

7 Расчетная

888

Tг 120x6

С255

66.80

66.80

0.83

7 Расчетная

901

Tг 120x6

С255

68.22

68.22

0.83

7 Расчетная

905

Tг 120x6

С255

68.22

68.22

0.83

7 Расчетная

916

Tг 120x6

С255

68.22

68.22

0.83

7 Расчетная

793

Tг 180x9

С255

43.13

43.13

0.84

7 Расчетная

827

Tг 180x9

С255

43.13

43.13

0.84

7 Расчетная

829

Tг 180x9

С255

43.13

43.13

0.84

7 Расчетная

840

Tг 180x9

С255

43.13

43.13

0.84

7 Расчетная

841

Tг 180x9

С255

43.13

43.13

0.84

7 Расчетная

853

Tг 180x9

С255

43.13

43.13

0.84

7 Расчетная

854

Tг 180x9

С255

43.13

43.13

0.85

7 Расчетная

897

Tг 120x6

С255

66.80

66.80

0.85

7 Расчетная

900

Tг 120x6

С255

68.22

68.22

0.85

7 Расчетная

889

Tг 120x6 Р2

С255

66.80

66.80

0.86

7 Расчетная

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

III-5