автореферат диссертации по строительству, 05.23.17, диссертация на тему:Экспериментально-теоретическое исследование устойчивости пространственных рамных систем и разработка инженерной методики определения критической силы с учетом нелинейности



На правах рукописи

Сон Марк Петрович

ЭКСПЕРИМЕНТАЛЬНО-ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ УСТОЙЧИВОСТИ ПРОСТРАНСТВЕННЫХ РАМНЫХ СИСТЕМ И РАЗРАБОТКА ИНЖЕНЕРНОЙ МЕТОДИКИ ОПРЕДЕЛЕНИЯ КРИТИЧЕСКОЙ СИЛЫ С УЧЕТОМ НЕЛИНЕЙНОСТИ

Специальность 05.23.17 - Строительная механика

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

- 9 ЛЕН 7П10

Москва-2010

004616984

Работа выполнена в Государственном образовательном учреждении высшего профессионального образования «Пермский государственный технический университет».

Научный руководитель: доктор технических наук, доцент

Кашеварова Галина Геннадьевна

Официальные оппоненты: доктор технических наук, профессор

Гайджуров Петр Павлович

доктор технических наук, профессор Косицын Сергей Борисович

Ведущая организация: Институт ОАО "НИЦ "Строительство" - ЦНИИСК им. В.А. Кучеренко, г. Москва

Защита состоится «21 » декабря 2010 г. в 14 — часов на заседании

диссертационного совета Д 212.138.12 при ГОУ ВПО Московский

государственный строительный университет по адресу: 129337, г. Москва, Ярославское шоссе, д. 26, ауд. 420 УЖ.

С диссертацией можно ознакомиться в библиотеке ГОУ ВПО Московский государственный строительный университет.

Автореферат разослан «1% » ноября 2010 г.

Ученый секретарь диссертационного совета

Анохин Н.Н.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

Расчетное обоснование проектных решений - многоэтапный процесс, в котором, как минимум, следует выделять две основные стадии: инженерный проект и выпуск рабочей документации. На стадии инженерного проекта технические решения целесообразно обосновывать приближенными (инженерными) расчетами, а на второй стадии следует применять численные методы и расчеты с использованием ЭВМ.

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

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

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

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

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

Задачи исследования

1. Проведение численного линейного и нелинейного анализа устойчивости равновесия пространственных рамных каркасов в разных программных комплексах.

2. Экспериментальные исследования образцов рамных каркасов на устойчивость и оценка адекватности (алгоритмической надежности) результатов численного расчета и эксперимента.

3. Разработка и обоснование основных теоретических положений приближенной (инженерной) методики оценки критической силы для многопролетных и многоэтажных рам.

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

5. Исследование границ применимости разработанной приближенной методики расчета устойчивости равновесия стержневых каркасов зданий.

Научную новизну диссертации составляют:

- численное исследование устойчивости пространственных рамных каркасов: влияния конечно-элементной дискретизации, модели деформирования материала и организации итерационного процесса в нелинейном анализе на результаты расчета;

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

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

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

Практическая значимость работы состоит:

> в создании авторских программ-макросов, которые могут стать составной частью при построении вычислительных комплексов;

> в разработке инженерной методики и алгоритма расчета обобщенной критической силы для пространственных каркасов зданий с учетом физической нелинейности материала конструкций, влияния поперечных нагрузок и случайных эксцентриситетов;

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

Основные разделы работы выполнялись в рамках тематического плана госбюджетных НИР по заданиям Министерства образования РФ по направлению «Совершенствование методов проектирования зданий, конструкций, оснований и фундаментов, санитарно-технических систем, дорожно-транспортных комплексов с учетом энергосбережения и экологии». Работа поддержана грантом РФФИ № 08-08-00702-а «Механика закритического деформирования и вопросы прочностного анализа».

Использование работы. Разработанные методики, алгоритмы и программы использовались для решения задач расчета строительных конструкций при выполнении научно-исследовательских работ в ПГТУ, а также при написании учебного пособия по строительной механике для студентов строительных специальностей в рамках приоритетного национального проекта «Образование» по программе Пермского государственного технического университета «Создание инновационной системы формирования профессиональных компетенций кадров и центра инновационного развития региона на базе многопрофильного технического университета».

Методы исследования, достоверность и обоснованность результатов. В

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

На защиту выносятся

1. Математическая модель и численная реализация задачи устойчивости пространственных рамных каркасов с учетом геометрической и физической нелинейности материала с привлечением разных программных комплексов промышленного типа;

2. Результаты экспериментальных исследований устойчивости металлических рамных каркасов;

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

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

Апробация работы. Основные положения и результаты работы обсуждались на: научно-практических конференциях строительного факультета ПГТУ «Строительство, архитектура. Теория и практика» (Пермь, 2004, 2006, 2007); 7-й всероссийской конференции «Информация инновации инвестиции» (Пермь, 2006); 25-й юбилейной межвузовской студенческой научно-технической конференции (Самара, 2006); Всесоюзных и международных зимних школах по механике сплошных сред (Пермь, 2007, 2009); Областной дистанционной научно-практической конференции ученых и студентов «Молодежная наука Прикамья» (Пермь, 2007); Втором Международном симпозиуме «Актуальные проблемы компьютерного моделирования конструкций и сооружений» (Пермь, 2008); Научной сессии «Новое в исследовании и проектировании пространственных конструкций» (НИИЖБ, Москва, 2009); XXIII Международной конференции «Математическое моделирование в механике твердых тел и конструкций. Методы граничных и конечных элементов» (Санкт- Петербург, 2009); Научной сессии «Проблемы нелинейного расчета большепролетных пространственных конструкций» (НИИЖБ, Москва, 2010); III Международном симпозиуме «Актуальные проблемы компьютерного моделирования конструкций и сооружений» (Новочеркасск, 2010).

Материалы диссертационной работы в целом обсуждались на научных семинарах кафедры «Строительная механика и вычислительные технологии», на расширенном заседании кафедр «Строительная механика и вычислительные технологии», «Строительные конструкции» и «Динамика и прочность машин» Пермского государственного технического университета и на научном семинаре Научно-образовательного центра компьютерного моделирования уникальных зданий, сооружений и комплексов МГСУ под руководством профессора A.M. Белостоцкого (Москва, 2010г.).

Публикации

Основное содержание работы отображено в 13 научных работах, из них 2 в научных журналах, рекомендованных ВАК.

Структура и объем работы

Диссертационная работа состоит из введения, четырех глав, заключения списка литературы, включающего 130 наименований, и приложения. Работа изложена на 146 листах машинописного текста, 3 страниц приложения и содержит 87 рисунков, 10 таблиц.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обосновывается актуальность темы, сформулированы цели и задачи исследования, отмечается научная новизна и практическая значимость полученных результатов, приведено краткое содержание работы.

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

Дается краткое описание развития теории устойчивости равновесия и экспериментальных исследований отдельного стержня, которая до сих пор не утратила своего значения. Отмечаются теоретические и экспериментальные исследования по данному вопросу И.Г. Бубнова, А.Н. Динника, Саусвелла, Ф.Энгессера, Ф.С. Ясинского, Т. Кармана, Ю.Н. Работнова, А. Дюло, И.Ходкинса, Э. Ламарля, М. Рэнкина, И. Баушингера, JL Тетмайера, А. Консидера, Ф. Шенли, Я.Г. Пановко, A.C. Вольмира, E.JI. Николаи.

Большой вклад в развитие теории устойчивости сооружений внесли A.B. Александров, H.A. Алфутов, В.В. Болотин, Б.М. Броуде, В.З. Власов, A.C. Вольмир, Г.Ю. Джанелидзе, Н.В. Корноухов, В.В. Новожилов, Я.Л. Нудельман, Я.Г. Пановко, П.Ф. Папкович, А.Р. Ржаницин, А.Ф. Смирнов, С.П. Тимошенко, В.И. Феодосьев.

В настоящее время проблема устойчивости равновесия строительных конструкций не потеряла своей актуальности. К ней обращались и обращаются многие исследователи: П.П. Гайджуров, И.Д. Грудев, В.Д. Клюшников, С.Б. Косицын, JI.C. Ляхович, A.B. Перельмутер, В.И. Сливкер, В.Д. Потапов, В.В. Улитин, В.Д. Райзер, Г. Циглер и др.

Своеобразие и нестандартность задач устойчивости равновесия неоднократно приводило к различного рода ошибкам, как в решениях конкретных прикладных задач, так и при обосновании отдельных положений теории устойчивости равновесия конструкций, о чем весьма поучительно изложено в книге Я.Г. Пановко и И.И. Губановой, и в книгах A.B. Перельмутера и В.И. Сливкера.

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

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

Во второй главе исследуются возможности программного комплекса ANS YS по решению задачи устойчивости пространственных рамных каркасов: • расчет потери устойчивости, связанный с вычислением собственных значений (он же линейный), и

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

В линейном расчете определение критической нагрузки сводится к определению наименьшего положительного собственного числа Xi для следующей системы уравнений: ([к]- Х\Кг J){u} = 0 (1)

где: [К] - обычная матрица жесткости конструкции; [К„] - матрица эффективной (или геометрической) жесткости; X - собственное значение (масштабный фактор), а величина Xi является коэффициентом запаса по устойчивости; {и} - собственный вектор, характеризующий форму потери устойчивости системы.

При решении задачи (1) предполагается, что все внутренние усилия растут пропорционально одному параметру - собственному числу, что связано с таким же пропорциональным ростом всех приложенных к системе нагрузок. Этот метод реализован практически во всех конечно-элементных программных системах. Для решения обобщенной проблемы собственных значений нами использовался блочный метод Ланцоша.

Нелинейное поведение конструкции является результатом множества причин, связанных с геометрической и физической нелинейностью. Если конструкция имеет большие деформации, изменение ее геометрической формы может вызвать нелинейный отклик. Геометрическая нелинейность характеризуется большими перемещениями и (или) поворотами

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

[Kg]£ = [Т„]т (3{F'nte }/5{ue}) + (5[Тп]т/3{ие}) {Finte} = = [Т„]т i[Bv]T(<3{ac}/3{uc}) d(vol) + [Т„]т Jö [Bv]T/3{ue}){ae} d(vol)+ (2)

(I) (И)

+ (ö[Tv]T/Ö{ue}){Fime} (III)

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

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

В программе ANSYS реализовано достаточно много вариантов моделей деформирования материалов. Нами были выбраны и исследованы следующие модели, рекомендуемые для описания поведения металлов: модель билинейного кинематического упрочнения (BKIN), использующая критерий текучести Мизеса, ассоциативный закон течения и кинематическое упрочнение; модель полилинейного кинематического упрочнения (MKIN), модель полилинейного изотропного упрочнения (MISO), которая использует сочетание условия текучести Мизеса с изотропным расширением поверхности пластичности и рекомендуется для проведения анализа при больших деформациях.

Для получения решения нелинейных задач устойчивости используется метод Ньютона - Рафсона. При этом вся нагрузка заменяется серией ее небольших приращений и выполнением на каждом таком шаге по нагрузке последовательности линейных приближений до получения состояния равновесия. Каждое линейное приближение требует выполнения равновесных итераций.

Подход, использующий процедуру последовательных приближений Ньютона-Рафсона, приводит к следующему соотношению, справедливому для некоторой равновесной итерации: [К] ы {Au}¡ = {F} - {Fel }ц, (3)

где: [К] м - матрица жесткости на предыдущей итерации; {Au}¡ - вектор, компонентами которого являются приращения перемещений двух последовательных итераций {u}¡ = {u }я - {Au}¡; {u}¡ - вектор перемещений, относящихся к текущей итерации; {F} - вектор приложенных к системе сил; {Fel }í-i - вектор упругих сил, соответствующих перемещениям предыдущей итерации с номером (i -1).

Для улучшения сходимости итерационного процесса применялась комбинация метода Ныотона-Рафсона с техникой метода ограничивающих дуг {arc-length method). Критической нагрузкой (нагрузкой, соответствующей потере устойчивости) является тот ее уровень, при котором решение начинает расходиться.

Для проведения вычислительных экспериментов нами разработаны программы - макросы на языке APDL. Вычислительные эксперименты по устойчивости пространственных конструкций каркасов проводились на моделях двух видов (рис. 1)

Рис. 1. Модели каркасов

Расчетные модели создавались балочными конечными элементами BEAM 189. Это квадратичный, 3-узловой, пространственный элемент, имеющий шесть или семь степеней свободы в каждом узле, седьмая степень свободы учитывает депланацию элемента. В основу описания этого элемента положена теория балки Тимошенко, в нем учитываются эффекты деформации сдвига и поддерживаются модели упругости, ползучести и пластичности материала для разных типов сечений. Этот элемент подходит для анализа линейных и нелинейных задач с учетом больших деформаций, а также для решения задач устойчивости конструкций с учетом изгибных и крутильных форм потери устойчивости. Кроме того, данный тип конечного элемента позволяет, используя разные варианты сечений, разбивать эти сечения, вводя дополнительные узлы по радиусу и окружности или по сторонам сечения.

В результате проведения линейного анализа устойчивости определялись формы потери устойчивости и критическая нагрузка. На рис. 2 показаны первые формы потери устойчивости испытуемых каркасов при действии вертикальной нагрузки Р= 1 на узлы. Обобщенные (суммарные) критические силы получились равными для 1-го и 2-го каркаса соответственно: Р® =123,7

Для уточнения величины критической силы и формы деформирования каркаса далее проводился нелинейный анализ. При этом для получения величины критической силы и графика зависимости нагрузки от перемещений, к узлам рамы (к тем же, что и единичные силы в линейном расчете) прикладывалось постепенно возрастающее кинематическое воздействие.

В процессе проведения вычислительных экспериментов по устойчивости каркаса 1-й модели исследовалось влияние на результаты расчета: модели деформирования материала; размера поперечного сечения стержней; дискретизации модели.

При использовании разных моделей пластичности, получен практически одинаковый характер зависимости и величины обобщенной критической силы (~ 94 кН). Поэтому далее мы использовали модель билинейного кинематического упрочнения (ВКШ).

кН, P'tf = 1044 кН.

Рис.2. Формы потери устойчивости каркасов

и

Исследование влияния размера поперечного сечения стержней (диаметра d) на устойчивость пространственного каркаса выявило следующую картину: при изменении диаметра стержней в конструкции меняется не только величина критической силы, но и форма потери устойчивости. При стержнях диаметром от 5 до 7 мм конструкция деформируется по крутильной форме, при диаметрах 8-9 - по смешанной, при диаметрах от 10 и выше - по изгибной форме.

На рис. 3 показаны зависимости критической силы от диаметра стержней в конструкции в линейном и нелинейном расчетах.

При d=9 мм значения критической силы

практически совпадают (гибкость }„=89).

Рис. 3. Зависимость величины критической силы от диаметра стержней в конструкции

Исследуя влияние конечно - элементной дискретизации конструкции на результаты расчета выявлено, что форма потери устойчивости зависит от числа конечных элементов N, на которое разбивается стержень по длине, в то время, как величина критической нагрузки при этом практически не меняется {Р[р =94 кН).

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

При расчете на устойчивость каркаса 2-го типа использовались разные варианты приложения нагрузки (узловая Р= 1 и распределенная q=l - на ригели (рис.8)). Кроме того, рассматривались разные варианты граничных условий (без закрепления верхних граней каркаса и с учетом контакта).

В линейном расчете первые формы потери устойчивости и значения минимальных обобщенных критических сил Р\р для разных вариантов приложения нагрузки получились достаточно близкими: 1063 кН и 1044 кН соответственно (разница 1,8%).

Критическая сила в нелинейном расчете получилась 692 кН. (к узлам рамы вертикально прикладывалось постепенно возрастающее кинематическое воздействие). При определении деформаций в зонах пластических шарниров использовалась более мелкая cenca по сечению стержня.

300 -—

1 -♦-Р крСм/л)

/ /

} / y

■TTÍC6B I13&284

■vrsf*" «Ss

■«ИЗГ^ 1—— íl, Lili

4 5 В 7 8 9 10 11 12 1

В третьей главе приводится описание проведенных натурных экспериментов двух моделей пространственных каркасов на устойчивость той же геометрии, что и в вычислительных экспериментах (по 3 каркаса каждого вида).

Для испытаний использовался гидравлический 500т пресс, точность по нагрузке - 16 кг, по деформациям - 0,0036 мм (рис.4). В процессе испытаний образцов регистрировались значения внешней нагрузки Р и перемещения в направлении действия нагрузки. Для снижения трения между образцом и верхней плитой пресса укладывался фторопласт. Нагрузка прикладывалась постепенно со скоростью 2мм/мин, а снималась после того, как была пройдена высшая точка на диаграмме деформирования.

На рис. 5 показаны формы потери устойчивости одного из образцов каждого вида (остальные образцы деформировались аналогично), а на рис. 6 приведены диаграммы изменения критической нагрузки в процессе деформирования каркасов. Два опыта из трех (2 и 3) проводились с разгрузкой. Кроме того, экспериментально определялись физико-механические свойства материала образцов, которые в дальнейшем использовались для уточненных расчетов.

Рис.4. Установка для испытания на устойчивость каркаса

Перемещение, мм

Рис. 5. Форма потери Рис.6. Графики экспериментальной зависимости

устойчивости каркасов суммарной нагрузки от перемещений

Значения критических сил в экспериментах на 1-й модели получились равными 97кН, 91кН, 92кН, а на 2-й 679кН, 730кН, 642 кН, соответственно.

Известно, что расчетная модель работы сооружения должна правильно и полно отражать работу реального объекта, а математические модели и методы расчета должны исходить из форм деформаций и разрушений, поэтому были скорректированы размеры, физико-механические свойства материала, граничные условия (учет контакта с прессом) в соответствии с результатами экспериментов. Формы деформирования каркасов показаны на рис. 7, 9, расчетные зависимости суммарной вертикальной силы Ркр от перемещений - на рис.10, значения критических нагрузок сведены в табл.1.

а)

Рис.7. Формы деформирования каркаса 1 (нелинейный расчет): а) без учета контакта; б) с учетом контакта с прессом а)

Рис. 8. Расчетная схема каркаса

т^шш

11 11111 11 »иМ! I II ч

Рис.9. Формы деформирования каркаса 2 (нелинейный расчет): а) без учета контакта; б) с учетом контакта с прессом

Суммарная силз Р, кН 100

0.15 о.:

Перемещения.!

Суммарная сила Р, кН

/

/

1

1

\

1

Рис. 10. Расчетные зависимости суммарной вертикальной силы Ркр от перемещений:

Вычислительные эксперименты устойчивости рамных каркасов были выполнены и в программном комплексе М1сгоРе: для 1-го каркаса по концепции Шенли при различных гипотезах работы материала (по диаграмме Прандтля, по унифицированной диаграмме строительных сталей) с учетом реальных свойств материала. Для расчета второго образца используется концепция "нулевой отпорности". Результаты расчета обобщенной критической силы, полученной теоретически и экспериментально, приведены в табл. 1.

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

Вводится допущение о том, что в плоской одноэтажной раме при потере устойчивости узловые вертикальные силы совершают работу на одинаковом перемещении 6 (см. рис. 11).

Рис.11. Произвольная многопролетная рама с постоянной высотой стоек

На основе введенного допущения доказывается независимость суммарной критической силы от ее распределения между стойками рамы, т.к. независимо от ее распределения в новом положении не меняется равновесие системы (сумма моментов относительно т.1 вырезанной части (рис.116)) Рг1 + (Р2.(11+1)-Р2 ./,) + „. + (*>,.(/,+/,+.. + /. + 1)-Р„ .(/,+/2 +.. + /,)) +

+ НгИ + Н2-И + ... + Ня-}1-М1-М2-...-М1п_1-Ми=0 ±Мк={Р1 +Р2+... + Р„)-1=Р-1,

»•I

и соблюдается вариационный принцип Лагранжа А'""' = Атуя.

Соответственно при расчете многопролетных рам с постоянной высотой стоек на устойчивость предлагается определять не критические нагрузки,

и

приложенные в узлах, а обобщенную критическую силуР = £ Р,. которую рама

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

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

несправедливым, но такие случаи не представляют большого интереса для практики.

В основе предлагаемой методики расчета лежит положение о том, что любую многопролетную раму можно представить совокупностью Г- образных и Т-образных рам, разрезав ее посередине ригелей. Т-образную раму всегда можно заменить Г-образной с учетом суммирования погонных жесткостей ригелей (при жестком сопряжении ригеля с колонной), а критическая сила для произвольной одноэтажной рамы будет равна сумме критических сил для Г-образных рам входящих в нее (рис.12).

При этом доказывается, что расчет любой Г-образной рамы (рис. 12а) можно заменить расчетом стержня с упругоподатливой на изгиб опорой с жесткостью Ир, свободной при этом от горизонтального перемещения.

(рис.12в)

а) б) в)

Рис. 12. Замена Г-образной рамы стержнем с упругоподатливой на изгиб опорой

Решая краевую задачу об устойчивости стержня с упругоподатливой на

изгиб опорой, получено решение и/ — + tan ni = 0 (4) или приближенно с

3г'Р

точностью 2,4%) в виде п =

1-

0.75

1 + -

(5)

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

формула:

р =Е

КО Lmi

1--

0.75

1 +

■ж

EIt

(6)

где ¿-номер стойки или ригеля; -погонная жесткость стойки; ¿/-погонная жесткость ригеля, п -количество участков разбиения, равное числу стоек рамы.

Чтобы показать универсальность предложенной формулы (6) и методики были рассмотрены примеры расчета многопролетных рам на устойчивость при действии вертикальной нагрузки, значения погонных жесткостей элементов и распределение нагрузки в которых принимались произвольно. Результаты

расчета сравнивались с результатами, полученными методом перемещений. Погрешность не превышала 2,5%.

Рассуждения, применяемые при выводе формулы для одноэтажной рамы можно применить и для многоэтажной рамы. При этом критическая сила будет определяться как сумма критических сил для рам, полученных путем рассечения исходной рамы посередине ригеля (рис.13). Приближенная формула для расчета многоэтажной рамы записывается в виде:

г \

1 _ °-75

! | +

Vf^. (7)

где Xt и хг - параметры, зависящие от количества этажей.

-BT-TL 1--—

v- 1 , 1.5/ А 1 ic х, = > —, х, = 1--- >„----, где m-число этажей.

I т 20 I т ip

Выполнена проверка формулы (7) при разных Рис.13. Часть соотношениях погонных жесткостей (if/ic=0,33-l0), и при многоэтажной рамы разном числе этажей (от 1 до 50). Результаты расчета по этой формуле сравнивались с результатами расчета на устойчивость многоэтажной рамы, полученными в программных комплексах SCAD и Лира. Разница не превышала 5%.

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

В основу полученных выше формул (6, 7), положено условие о том, что сжатый стержень работает согласно закону Гука.

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

1. Выполнить расчет (1-ю итерацию) по формуле (7) в предположении упругой работы стержней.

2. В зависимости от полученной критической силы найти гибкости сжатых

- , я I стержней Л =--

V Р

3. По найденным значениям гибкости Я определить коэффициенты приведенного модуля г, (при упругой работе стержня г = 1, при упруго-

пластичной т- — -(1-р(1- гдер = — ).

Е ж Е ж Е ат

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

Итерации продолжаются до тех пор, пока разница между двумя соседними приближениями не окажется в пределах нужной погрешности. В каждой итерации расчет рамы ведется в предположении, что все стержни работают упруго по формуле (7).

5. Зная гибкости стержней, определить а по формуле

Л2

а =ат -(оу - Oy)-- и соответствующую критическую силу.

и (min Я)2

Описанные выше алгоритмы легко реализуются в Excel. Хорошая сходимость достигается на 2-4 шаге расчета.

Для учета случайных

эксцентриситетов необходимо выполнить расчет рамы или рамного каркаса по описанной методике и умножить вычисленные критические силы

(критические напряжения) на

соответствующие коэффициенты <рг (рис.14), которые имеют статистическую природу, или после нахождения гибкостей на последней итерации расчета учесть коэффициенты <р по СНиП.

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

о 100 ш л

Рис.14. Влияние случайных эксцентриситетов в зависимости от гибкости стеожня

пластических деформаций может происходить перераспределение продольных сил в системе.

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

В предлагаемой методике расчета это можно учесть следующим образом. Поскольку продольные силы остаются неизменными, то расчет на устойчивость рамы или рамного каркаса сначала ведем только на действие вертикальной нагрузки по описанной методике. Произведя расчет и найдя приведенные жесткости стержней, находим распределение момента между стойками рамы (пропорционально изгибным жесткостям стоек). Суммарный момент по краям стоек первого этажа равен произведению суммарной горизонтальной реакции умноженной на высоту этажа. Можно считать, что 66% момента воспримет опорное сечение, а 33% - второе.

Зная расчетные продольные силы, моменты, гибкости и форму поперечного сечения, находим коэффициенты <р, по СНиП «Стальные конструкции для случая внецентренно сжатого стержня». Общая критическая сила получается суммированием критических нагрузок для отдельных стержней.

Таблица 1. Сравнение результатов расчетов и экспериментов обобщенной критической силы

Критическая сила, кН

Каркас Инженерная методика МкгоГе Эксперимент

Лин. расч Н/лин. расч. С учетом случ. зксцентр. Лин. расчет Н/лин расчет Лин. расч Н/лин. расч. 3 опыта

1 120,8 96,5 71,3 123,7 94, 96- (контакт) 95,3 90,3- (разн. модели) 97,91,92

2 1048 702 605,8 1063 692, 752- (контакт) 1062,4 678,6 617,4 - (случ. эксцен) 679,730,642

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

1. Анализ численных результатов устойчивости равновесия пространственных рамных систем показал, что линейный расчет в разных программных комплексах (ANSYS, MicroFE, SCAD, ЛИРА-Windows) дает практически одинаковые, но завышенные значения критической нагрузки для конструкций средней и малой гибкости.

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

3. Экспериментальные исследования образцов рамных каркасов на устойчивость показали, адекватность экспериментальных и теоретических результатов, формы потери устойчивости, полученные в программных комплексах, а также по разработанной приближенной методике оказались одинаковыми. Разброс в значениях критических нагрузок вычисленных в ANSYS, MicroFE и по предложенной инженерной методике не превышает 5%. Разница же между теоретическими и экспериментальными величинами нагрузок составляет 12-15%, что можно объяснить неточностью изготовления, установки на прессе каркаса и наличием сил трения в эксперименте.

4. Сформулированы основные теоретические положения и допущения приближенной (инженерной) методики оценки критической силы для многопролетной рамы; выполнена оценка погрешности (границы применимости) введенных допущений.

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

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

Основные положения и результаты диссертации опубликованы в следующих работах (общее количество -18; ниже перечислено - 13):

1. М.П.Сон. Расчет систем с бесконечно жесткими элементами методом перемещений (статья). Молодежная наука Прикамья:Сб. науч. тр. Вып. 6./ Перм. гос. тех. ун-т-Пермь, 2005.С. 170-173

2. Сон МЛ. Новый подход к расчету многопролетных и многоэтажных рам на устойчивость (статья). Информация, инновация, инвестиция. Материалы 7-й Всерос. Конф., гЛермь 29-30 нояб. 2006г. Пермь/М-во пром-сти и энергетики РФ [и др.].-2006. Пермь, с. 108-113.

3. Сон МЛ. Разработка методики расчета каркаса здания с учетом физической нелинейности (статья) Информация, инновация, инвестиция. Материалы 7-й Всерос.

Конф., г.Пермь 29-30 нояб. 2006г. Пермь/М-во пром-сти и энергетики РФ [и др.].-2006. Пермь.с.113-117.

4. Кашееарова Г.Г., Сон М.П. Новый подход к определению прогибов в стержневой системе с учетом образования пластических шарниров (статья) Зимняя школа по механике сплошных сред (пятнадцатая). Сб. ст. В 3-х частях. Ч.2./ Ин-т механики сплошных сред: УрО РАН [и др.].-Пермь, 2007.С.134-137

5. М.П.Сон. Уравнение идеально устойчивого стержня (статья). Строительство и образование: Сб. науч. тр./ ГОУ ВПО Урал, гос.техн. ун-т-УПИ.-Екатеринбург, 2007. с.42-43

6. Кашееарова Г.Г., Сон М.П. Новый метод расчета рамных систем на устойчивость (статья). Научные исследования и инновации.-2007.-№1 Перм. гос. тех. ун-т.- Пермь. - с.46-54.

7. Сон М.П. Расчет рамных каркасов на устойчивость (статья). International journal for computational civil and structural engineering.-2008.-Vol.4, iss. 2.c.125-127.

8. Сон МЛ., Кашееарова Г.Г. Теоретически и экспериментальные исследования устойчивости многопролетных рам (статья, электронный рессурс). Механика сплощных сред как основа современных технологий: тр. XVI Зимней шк. по механике сплошных сред [Пермь 24-27 февраля 2009г]/Урал. отд-ние Рос. акад. наук, [и др.].-Пермь,2009

9. Сон М.П., Кашееарова Г.Г. Теоретические и экспериментальные исследования устойчивости многопролетных рамных каркасов зданий (тезисы). Математическое моделирование в механике деформируемых тел и конструкций. Методы граничных и конечных элементов.:ХХШ междунар. конф.: тез., г. Санк- Петербург, 28сент.-01окт. 2009г./Ин-т проблем машиноведения Рос. акад. наук . [и др.].-СПб.,2009.с.199-201

10. Кашееарова Г.Г.,Сон М.П. Устойчивость рамных каркасов за пределом упругости. Пространственные конструкции зданий и сооружений (Исследование, расчет, проектирование и применение): Сб. статей. Вып. 12 / МОО «Содействие развитию и применению пространственных конструкций в строительстве», Научный совет РААСН "Пространственные конструкции зданий и сооружений". М.: 2009. - 240 с.(С.52-58)

11. Сон М.П. Инженерная методика расчета устойчивости многопролетных многоэтажных рамных каркасов зданий. Научно-технический журнал «Строительная механика и расчет сооружений». № 2.М.2010. с.41-47

12. Кашееарова Г.Г., Сон М.П. Экспериментально-теоретические исследования рамных каркасов на устойчивость. Материалы III международного симпозиума «Актуальные проблемы компьютерного моделирования конструкций и сооружений. Новочеркасск, ЮРГТУ (НПИ), 21-24 июня 2010. с.90-92

13. Сон М.П. Методики расчета на устойчивость металлических рамных каркасов за пределом упругости. Тезисы докладов научной сессии «Проблемы нелинейного расчета большепролетных пространственных конструкций». / МОО «Содействие развитию и применению пространственных конструкций в строительстве», Научный совет РААСН "Пространственные конструкции зданий и сооружений"/ 20 апреля Москва.2010.с.55-5б

Издательство

Пермского государственного технического университета 614990, г. Пермь, Комсомольский пр., 29, к.113 тел. (342) 219-80-33

ВВЕДЕНИЕ.

ГЛАВА 1. АНАЛИТИЧЕСКИЙ ОБЗОР СОВРЕМЕННЫХ МЕТОДОВ РАСЧЕТА ПРОСТРАНСТВЕННЫХ РАМ И ИХ ЭЛЕМЕНТОВ НА УСТОЙЧИВОСТЬ.

1.1. Обзор основных этапов и положений теории устойчивости равновесия стержневых систем.

1.2. Анализ современных подходов, требований нормативных документов и возможностей оценки устойчивости равновесия пространственных рамных систем.

1.3. Аналитический обзор основных методов расчета на устойчивость стержневых систем

1.4. Обзор экспериментальных исследований на устойчивость стержневых систем.

1.5. Обоснование состава и структуры частных задач исследования.

ВЫВОДЫ по главе.

ГЛАВА 2. ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ПРОБЛЕМЫ УСТОЙЧИВОСТИ РАМНЫХ КАРКАСОВ.

2.1. МКЭ в задачах устойчивости строительных конструкций.

2.1.1. Разрешающие уравнения МКЭ в задачах устойчивости равновесия.

2.1.2. Учет физической нелинейности материала. Основные соотношения теории пластичности.

2.2. Возможности и проблемы расчета устойчивости программного комплекса А^УБ

2.2.1. Методы расчета устойчивости, реализованные в

2.2.2. Математические модели неупругого деформирования материала.

2.3. Вычислительные эксперименты по расчету на устойчивость рамных каркасов в ПК АШУБ.

2.3.1. Расчет на устойчивость модели 1 рамного каркаса.

2.3.2. Расчет на устойчивость .модели 2 рамного каркаса.

ВЫВОДЫ по главе.

ГЛАВА 3. ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ МОДЕЛЕЙ

ПРОСТРАНСТВЕННЫХ КАРКАСОВ НА УСТОЙЧИВОСТЬ.

3.1. Натурные испытания образцов рамных каркасов на устойчивость.

3.1.1. Испытания моделей каркаса 1-го типа.

3.1.2. Определение физико-механических свойств материала каркаса 1-го типа.

3.1.3. Испытания моделей каркаса 2-го типа.

3.1.2. Определение физико-механических свойств материала каркаса 2-го типа.

3.2. Сопоставление результатов расчета с результатами экспериментов.

3.2.1. Оценка адекватности расчетных моделей каркасов в ПК ЛА'ЯУЯ.

3.2.2. Проверка устойчивости модели рамного каркаса в программном комтексе М'югоРе.

ВЫВОДЫ по главе.

ГЛАВА 4. РАЗРАБОТКА ТЕОРЕТИЧЕСКИХ ПОЛОЖЕНИЙ И НОВОГО ИНЖЕНЕРНОГО МЕТОДА РАСЧЕТА МНОГОПРОЛЕТНЫХ МНОГОЭТАЖНЫХ РАМ И РАМНЫХ КАРКАСОВ НА УСТОЙЧИВОСТЬ.

4.1. Разработка основных теоретических положений расчета многопролетной рамы на устойчивость.

4.2. Методика расчета многоэтажных рам на устойчивость.

4.3. Методика расчета пространственных каркасов.

4.4. Учет физической нелинейности материала конструкции при расчете каркаса на устойчивость (2 рода).

4.5. Учет влияния на критическую силу случайных эксцентриситетов и поперечных нагрузок.

4.6. Расчет устойчивости сжато-изогнутых и внецентренно-сжатых стержней и конструкций.

4.7. Расчет на устойчивость моделей натурного эксперимента по предлагаемой методике и верификация результатов.

4.7.1. Расчет каркаса 1.

4.7.2. Расчет каркаса 2.

4.7.3. Сравнение результатов расчетов и экспериментов.

ВЫВОДЫ по главе.

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

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

Расчетное обоснование проектных решений - это многоэтапный процесс, в котором, как минимум, следует выделять две основные стадии: инженерный проект и выпуск рабочей документации. На стадии инженерного проекта технические решения целесообразно обосновывать приближенными (инженерными) расчетами, а на второй стадии при рабочем проектировании следует применять численные методы и расчеты с использованием ЭВМ.

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

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

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

Объект исследования — многопролетные многоэтажные рамные каркасы зданий.

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

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

Задачи исследования

1. Проведение численного линейного и нелинейного анализа устойчивости равновесия пространственных рамных каркасов в разных программных комплексах.

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

3. Разработка и обоснование основных теоретических положений приближенной (инженерной) методики оценки критической силы для многопролетной рамы.

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

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

6. Исследование границ применимости разработанной приближенной методики расчета устойчивости равновесия стержневых каркасов зданий.

Научную новизну диссертации составляют:

- численное исследование устойчивости пространственных рамных каркасов: изучение влияния конечно-элементной дискретизации, модели деформирования материала и организации итерационного процесса нелинейного расчета на результаты расчета;

- разработка авторских программ-макросов для расчета на устойчивость рамных каркасов с учетом геометрической и физической нелинейности на языке параметрического программирования АРБЬ; создание приближенной (инженерной) методика определения обобщенной критической силы при расчете на устойчивость многопролетных, многоэтажных рам и пространственных каркасов зданий с учетом физической 6 нелинейности материала, влияния поперечных нагрузок и случайных эксцентриситетов.

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

Практическая значимость работы состоит: в создании авторских программ-макросов, которые могут стать составной частью при построении вычислительных комплексов; в разработке приближенной методики и алгоритма расчета обобщенной критической силы для пространственных каркасов зданий с учетом физической нелинейности материала конструкций, влияния поперечных нагрузок и случайных эксцентриситетов; в выдаче практических рекомендаций для проектировщиков по реализации разработанных алгоритмов нелинейного расчета многопролетных, многоэтажных рам и рамных каркасов на устойчивость.

Основные разделы работы выполнялись в рамках тематического плана госбюджетных НИР по заданиям Министерства образования РФ по направлению «Совершенствование методов проектирования зданий, конструкций, оснований и фундаментов, санитарно-технических систем, дорожно-транспортных комплексов с учетом энергосбережения и экологии». Работа поддержана грантом РФФИ № 08-08-00702-а «Механика закритического деформирования и вопросы прочностного анализа».

Внедрение работы состоит в использовании разработанных методик, алгоритмов и программ для решения задач расчета строительных конструкций при выполнении научно-исследовательских работ в ПГТУ, а также при написании учебного пособия по строительной механике для студентов строительных специальностей в рамках приоритетного национального проекта «Образование» по программе Пермского государственного технического университета «Создание инновационной системы формирования профессиональных компетенций кадров и центра инновационного развития региона на базе многопрофильного технического университета».

Методы исследования, достоверность и обоснованность результатов. В работе использованы теоретические, экспериментальные и эмпирические методы исследования. Решения задач базируются на известных теоретических положениях и принципах строительной механики, теории сопротивления материалов, теоретической механики, математического моделирования и экспериментальных данных. Достоверность полученных результатов подтверждается корректностью математической постановки задачи, использованием известных численных методов решения краевых задач с привлечением верифицированных программных комплексов промышленного типа, сходимостью полученных теоретических результатов с данными натурных экспериментов.

На защиту выносятся

1. Математическая модель и численная реализация задачи устойчивости пространственных рамных каркасов с учетом геометрической и физической нелинейности материала с привлечением разных программных комплексов промышленного типа;

2. Результаты экспериментальных исследований по изучению устойчивости металлических рамных каркасов;

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

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

Апробация работы. Основные положения и результаты работы обсуждались на: научно-практических конференциях строительного факультета ПГТУ «Строительство, архитектура. Теория и практика» (Пермь, 2004, 2006, 2007); 7-й всероссийской конференции «Информация инновации инвестиции» (Пермь, 2006); 25-й юбилейной межвузовской студенческой научно-технической конференции (Самара, 2006); Всесоюзных и международных зимних ппсолах по механике сплошных сред (Пермь, 2007, 2009); Областной дистанционной научно-практической конференции ученых и студентов «Молодежная наука Прикамья» (Пермь, 2007); Втором Международном симпозиуме «Актуальные проблемы компьютерного моделирования конструкций и сооружений» (Пермь, 2008); Научной сессии «Новое в исследовании и проектировании пространственных конструкций» (НИИЖБ, Москва, 2009); XXIII Международной конференции «Математическое моделирование в механике твердых тел и конструкций. Методы граничных и конечных элементов» (Санкг- Петербург, 2009), Научной сессии «Проблемы нелинейного расчета большепролетных пространственных конструкций» (НИИЖБ, Москва, 2010), III Международном симпозиуме «Актуальные проблемы компьютерного моделирования конструкций и сооружений» (Новочеркасск, 2010).

Материалы диссертационной работы в целом обсуждались на научных семинарах кафедры «Строительная механика и вычислительные технологии», на расширенном заседании кафедр «Строительная механика и вычислительные технологии», «Строительные конструкции» и «Динамика и прочность машин» Пермского государственного технического университета и на научном семинаре Научно-образовательного центра компьютерного моделирования уникальных зданий, сооружений и комплексов МГСУ под руководством профессора А.М. Белостоцкого (Москва, 2010).

Имеются некоторые сомнения относительно возможности чисто экспериментального подхода к решению задач устойчивости о расчете жестких рам из-за большого числа и широких интервалов изменения параметров, характеризующих каждую отдельную задачу. Этими параметрами являются осевые силы, отношение жесткостей, форма поперечного сечения и всевозможные условия поперечного нагружения. Успеха можно добиться только при тесной связи теоретических и экспериментальных исследований, что вызывает необходимость проведения серии опытов, условия которых соответствуют допущениям, принятым в аналитическом решении задачи. Первоочередная цель таких опытов состоит в получении экспериментальных данных для определения тех коэффициентов, теоретически значения которых нельзя вывести или вычисление которых громоздко."

Блейх Ф.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

1. Анализ численных результатов устойчивости равновесия пространственных рамных систем показал, что линейный расчет в разных программных комплексах (ANSYS, MicroFe, SCAD, ЛИР A-Windows) дает практически одинаковые, но завышенные значения критической нагрузки для конструкций средней и малой гибкости.

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

3. Экспериментальные исследования образцов рамных каркасов на устойчивость показали, адекватность экспериментальных и теоретических результатов, формы потери устойчивости, полученных в программных комплексах, а также по разработанной приближенной методике оказались одинаковыми. Разброс в значениях критических нагрузок вычисленных в ANS YS, MicroFE и по предложенной методике не превышает 5%. Однако отличие от экспериментальных нагрузок составляет 12-15%, что объясняется наличием сил трения пресса о ригели каркаса.

4. Сформулированы основные теоретические положения и допущения приближенной (инженерной) методики оценки критической силы для многопролетной рамы; выполнена оценка погрешности (границы применимости) введенных допущений.

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

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

1. Александров A.B., Лащеннков Б.Я. О применении энергетического метода в задачах устойчиости упругих систем.// Строительная механика и расчет сооружений. 1965, №5.-с.28-32.

2. Алфутов H.A. Основы расчета на устойчивость упругих систем -Машиностроение, М., 1991.-334с.

3. Алфутов H.A., Колесников К. С. Устойчивость движения и равновесия : Учеб. для вузов / Под ред. К.С. Колесникова. М.: Изд-во МГТУ, 2003.-253с.

4. Анохин H.H. Строительная механика в примерах и задачах: Учебное пособие для вузов. ч.П. Статически неопределимые системы. 2000.- 353с.

5. Басов К.A. ANSYS. Справочник пользователя. М.: ДМК Пресс, 2005. — 640 с.

6. Безухов Н.И., Лужин О.В. Устойчивость и динамика в примерах и задачах.-М.Стройиздат, 1963.-371.

7. Беленя Е.М. Металлические конструкции, изд.б-е, М., Стройиздат, 1986. -560с.

8. Вельский Г.Е., Одесский П.Д. О едином подходе к использованию диаграмм работы строительных сталей // Промышленное строительство, 1980, №7, с.4-6.

9. Биргер И.А., Мавлютов P.P. Сопротивление материалов: Учебное пособие. -М.-.Наука. 1986. 560с.

10. Биргер И.А., Пановко Я.Г. Прочность, устойчивость, колебания. Справочник. ТомЗ.-М.: Машиностроение, 1968.-567с.

11. Блейх Ф. устойчивость металлических конструкций, перев. с англ., Физматиз, 1959.-544с.

12. Болотин В.В. Нелинейная теория упругости и устойчивости в болыном.//Расчеты на прочность.-М.: Машгиз, 1958, №3.-с.6-27.

13. Болотин B.B. О вариационных принципах теории упругой устойчивости.// Проблемы механики твердого деформируемого тела.-Л.: Судостроение, 1973.-С.83-88.

14. Болотин В.В. О понятии устойчивости в строительной механике.//Проблемы устойчивости в строительной механике.-М.:Стройиздат, 1965.-С.6-27.

15. Болотин В.В. О сведении трехмерных задач теории упругой устойчивости к одномерным и двухмерным задачам./ДТроблемы устойчивости в строительной механике. -М.:Стройиздат, 1965.-С. 186-196.

16. Болотин В.В. Динамическая устойчивость упругих систем, Гостехиздат, М., 1956.

17. Болотин В.В. Неконсервативные задачи теории упругой устойчивости, Физматгиз, М., 1961.-339с.

18. Броуде Б.М. Теория устойчивости и принципы расчета конструкций. -Проблемы устойчивости в строительной механике. -М.:Стройиздат, 1965.-С.28-43.

19. Броуде Б.М. Предельные состояния стальных балок, Госстройиздат, М.-Л., 1953.

20. Бубнов И.Г. Строительная механика корабля. Т.1, 1912.

21. Ведеников Г.С. (под ред.) Металлические конструкции: Общий курс: Учебник для вузов. 7-е изд., перераб. идоп.-М. Стройиздат, 1998.,-760с.

22. Власов В.З. Тонкостенные упругие стержни. Физматгиз, М., 1959.

23. Вольмир A.C. Гибкие пластины и оболочки, Гостехиздат, М., 1956.

24. Вольмир A.C. Устойчивость упругих систем. М., 1963.-879с.

25. Гайджуров П.П. Расчет стержневых систем на устойчивость и колебания: учебное пособие. Юж.-Рос. гос. техун-т. -Новочеркасск: ЮРГТУ, 2009. -195с.

26. Галеркин Б.Г. Стержни и пластинки, «Вестник инженера» №19, 1915.

27. Галеркин Б.Г. Теория продольного изгиба и опыт применения теории продольного изгиба к многоэтажным стержням, стойкам с жесткими соединениями и системам стоек, Собр. Соч., т.1, изд. АН СССР М., 1952.

28. Галеркин Б.Г. Стержни и пластинки. Собр. Соч., т.1. изд-во АН СССР, М., 1952.

29. Гастев В. А. Краткий курс сопротивления материалов. Физматгиз., 1959.

30. Гемерлинг A.B. О несущей способности сжатых стальных конструкций, Труды ЦНИПС, №7, М.1952.

31. Геммерлинг A.B. Несущая способность стержневых стальных конструкций, Госстройиздат, 1958.

32. Грудев И.Д. Устойчивость стержневых элементов в составе стальных конструкций. М: МИК, 2005.-320с.

33. Губанова И.И., Пановко Я.Г. О влиянии докритического обжатия стержня на критическое значение сжимающей силы.// Известия АН СССР, МТТ, 1968, №2.-С. 153-154.

34. Джанелидзе Г.Ю. Устойчивость равновесия нелинейно-деформируемых систем/Труды Ленингр. Политехи, ин-та, 1955, №178.

35. Динник А.Н. Продольный изгиб и кручение, М. 1955-392с.

36. Динник А.Н. Избранные труды, т. 1, Киев, 1952;

37. Дыховичный Ю.А. Конструирование и расчет жилых и общественных зданий повышенной этажности. М., Стройиздат, 1970.

38. Жемочкин Б.Н. Расчет рам. Издательство литературы по строительству. М. 1965.

39. Кашеварова Г.Г., Сон М.П Новый метод расчета рамных систем на устойчивость (статья). Научные исследования и инновации.-2007.-№1 Перм. гос. тех. ун-т — Пермь. — с.46-54.

40. Клейн Г.К., Рекач В.Г., Розенблат Г.И. Руководство к проведению занятий по специальному курсу строительной механики, Высшая школа 1964,297 с.

41. Клюшников В.Д. Устойчивость упруго-пластических систем.-М.:Наука, 1980.-240с.

42. Корноухов Н.В. (под редакцией), Исследование устойчивости пространственного каркаса по типу высотной части Дворца Советов, Изд. АН УССР, 1938.

43. Корноухов Н.В. Прочность и устойчивость стержневых систем. —М.: Госстройиздат, 1949.-376с.

44. Корчак М.Д., Галкин С.В., Картопольцев В.М. Основы неустойчивости в теории катастроф инженерных сооружений.-Томск: Изд-во Томского университета, 1997.-124с.

45. Лагранж Ж. Аналитическая механика. T.I-II, М.-Л.:ГИТТЛ, 1950.

46. Лейтес С.Д., Устойчивость сжатых стальных стержней, Госстройиздат, 1954.

47. Моисеев Н.Д. Очерк развития теории устойчивости.-М.: Гостехиздат, 1949.- 218с.

48. Николаи Е.Л., К вопросу об устойчивости скрученного стержня. Вестн. прикл. матем. мех. 1,1929.

49. Николаи Е.Л. Об устойчивости прямолинейной формы равновесия сжатого и скрученного стержня. Труды по механике, Гостехиздат, 1955.

50. Николаи Е.Л., Труды по механике. -М.: Гостехтеориздат, 1955.- 582с.

51. Нормы и технические условия проектирования стальных конструкций (НиТУ 121-55), Госстройиздат, 1952.

52. Нудельман Я. Л. Монография. Методы определения собственных частот и критических сил для стержневых систем, М.-Л. 1949.

53. Нудельман Я.Л., Кронберг В.А. О методе наложения в задачах устойчивости упругих систем. // Исследования по теории сооружений, Вып.12,-М.Госстройиздат, 1963.- С.89-100.

54. Пановко Я.Т., Губанова И.И. Устойчивость и колебания упругих систем. Современные концепции, парадоксы и ошибки.- М.: Наука, 1987.-352с.

55. Папкович П.Ф. Труды по строительной механике корабля. Т.4. Устойчивость стержней, перекрытой и пластин-Л.: Судостроение, 1963.-552 с.

56. Перельмутер A.B., Сливкер В.И. Расчетные модели сооружений и возможность их анализа. -М.: ДМК пресс, 2007.- 597с.

57. Перельмутер A.B., Сливкер В.И. Устойчивость равновесия конструкций и родственные проблемы. Том 1. М.: Издательство СКАД СОФТ, 2007 - 670 с.

58. Подольский Д.М. Пространственный расчет зданий повышенной этажности.-М.: Стройиздат, 1975.-158с.

59. Пособие по проектированию стальных конструкций (к СНиП П-23-81*. Стальные конструкции) / ЦНИИСК им.Кучеренко. -М.: ЦИТП Госстроя СССР, 1989.-148с.

60. Постон Т., Стюарт И. Теория катастроф и ее приложения. М.: Мир, 1980.- 608с.

61. Прокофьев И.П., Смирнов А.Ф. Теория сооружений, ч. III, Трансжелдориздат, 1948.

62. Рабинович И.М. Курс строительной механики, ч. II, Госстройиздат, 1954.

63. Рабинович И.М. Основы строительной механики стержневых систем, Госстройиздат, 1960.

64. Работнов Ю.Н. О равновесии сжатых стержней за пределом пропорциональности // Инженерный сборник.-1952.-Т.Х1.-С. 123-126.

65. Работнов Ю.Н. Сопротивление материалов. Межгосударственное изд-во физико-математической литературы. 1962. 456 с.

66. Раевский А.Н. Основы расчета сооружений на устойчивость. Высшая школа.-1962.

67. Раевский А.Н. Устойчивость одноярусных однопролетных рам, Изд. ЛИСИ, вып. 36, 1958.

68. Раевский А.Н. Формулы и графики для определения расчетных длин сжатых стоек одноярусных рам, "Строительная механика и расчет сооружений", №5, 1959.

69. Ржаницин А.Р. Устойчивость равновесия упругих систем, Гостеортехиздат, 1955.-475с.

70. Рогицкий С.А. Устойчивость стержневых систем с неподвижными узлами, Тр. УПИ, сб. 54, Госстройиздат,1955.

71. Розин Л.Л. и др. Расчет статически неопределимых стержневых систем. Учебное пособие. 1988.

72. Рэнкин У.Д.М. Руководство для инженеров строителей.- СПб.: 1870.

73. Саргсян А.Е. ( под ред.) Строительная механика: Учебник для вузов/ Саргсян А.Е., Дименко А.Т. и др., 2000.

74. Семенов В.А., Лебедев В.Л., Солдатов А.Ю. К устойчивости стержневых строительных конструкций с учетом физической нелинейности // тезисы докладов научной сессии "проблемы нелинейного расчета большепролетных пространственных конструкций".М.: НИИЖБ, 2010.

75. Сливкер В.И. Строительная механика. Вариационные основы. -М.6 Изд-во Ассоциации строительных вузов, 2005.-708с.

76. Смирнов А Ф. Устойчивость и колебания сооружений. -М.: Трансжелдориздат, 1958.-571с.

77. Смирнов А.Ф. Статическая и динамическая устойчивость сооружения. -М.:, 1947.

78. СНиП 11-23-81*. Стальные конструкции. Госстрой СССР. М.:ЦИТП Госстроя СССР, 1991.-96с.

79. Снитко H.K. Строительная механика, изд-во высшая школа, М. 1968.

80. Снитко Н.К. Устойчивость сжатых и сжато-изогнутых стержневых систем, Госстройиздат, 1956.

81. Снитко Н.К. Устойчивость стержневых систем, Госстройиздат, 1952.

82. Сон М.П. Инженерная методика расчета устойчивости многопролетных многоэтажных рамных каркасов зданий. Научно-технический журнал «Строительная механика и расчет сооружений». № 2.М.2010. с.41-47.

83. Сон М.П. Новый подход к расчету многопролетных и многоэтажных рам на устойчивость (статья). Информация, инновация, инвестиция. Материалы 7-й Всерос. Конф., г.Пермь 29-30 нояб. 2006г. Пермь/М-во пром-сти и энергетики РФ и др..-2006. Пермь, с. 108-113.

84. Сон М.П. Разработка методики расчета каркаса здания с учетом физической нелинейности (статья) Информация, инновация, инвестиция. Материалы 7-й Всерос. Конф., г.Пермь 29-30 нояб. 2006г. Пермь/М-во пром-сти и энергетики РФ и др..-2006. Пермь.с.113-117.

85. Сон М.П. Расчет рамных каркасов на устойчивость (статья). International journal for computational civil and structural engineering.-2008.-Vol.4, iss. 2.C.125-127.

86. Сон М.П. Расчет систем с бесконечно жесткими элементами методом перемещений (статья). Молодежная наука Прикамья:Сб. науч. тр. Вып. 6./ Перм. гос. тех. ун-т-Пермь, 2005.с. 170-173.

87. Сон М.П. Уравнение идеально устойчивого стержня (статья). Строительство и образование: Сб. науч. тр./ ГОУ ВПО Урал, гос.техн. ун-т-УПИ.-Екатеринбург, 2007. с.42-43.

88. СП 53-102-2004. Свод правил по проектированию и строительству. Общие правила проектирования стальных конструкций.

89. СТО 02494680-0049-2005. Стандарт организации. Конструкции стальные строительные. Центральный орден трудового красного знамени научно-исследовательский и проектный институт строительных металлоконструкций им. Н.П.Мельникова.

90. Стрелецкий Н.С. Работа сжатых стоек, Госстройиздат, 1959.

91. Тимошенко С.П. Об устойчивости упругих систем, «Известия Киевского политехнического института», 1910.

92. Тимошенко С.П. Устойчивость упругих систем (с дополнением проф. В.З.Власов), Гостехиздат, 1946-567с.

93. Улитин В.В. Физически нелинейный анализ устойчивости конструкций. -СПб.:ГИОРД, 2007. 96 с.

94. Феодосьев В.И. Десять лекций-бесед по сопротивлению материалов. -М. Наука, 1975.-174с.

95. Ханджи В. В. Расчет многоэтажных зданий со связевым каркасом. М., 1977.-187С.

96. Харт Ф., Хенн В., Зонтаг X. Атлас стальных конструкций./Многоэтажные здания. Пер. с нем. М., 1977.

97. Циглер Г. Основы теории устойчивости конструкций. М.: Мир, 1971.-192с.

98. Чигарев A.B., Кравчюк A.C., Смалюк А.Ф. ANS YS для инженеров. Справочное пособие. М.: Машиностроение-1. 2004. 512с.

99. Чувикин Г.М., Устойчивость рам и стержней, Госстройиздат, 1952.

100. Чудновский В.Г., Методы расчета колебаний и устойчивости упругих систем, изд. АН УССР, Киев, 1952.

101. Штаерман И .Я., Пиковский A.A. Основы теории устойчивости строительных конструкций. Госстройиздат, 1939.

102. Эйлер JL Метод нахождения кривых линий, обладающих свойством максимуму или минимума. Русское издание «Гостехиздат, 1934.

103. Ясинский Ф.С. Избранные работы по устойчивости сжатых стержней. -М.: -JL: Госстройиздат, 1952.-427с.

104. Ясинский Ф.С. О сопротивлении продольному изгибу (диссертация). // Ф.С.Ясинский. Избранные работы по устойчивости сжатых стержней.-M.-JI.: Гостехиздат, 1952.-е. 11-137.

105. ANSYS Basic Analysis Procedures Guide. ANSYS Release 11. ANSYS Inc., 2008.

106. Bauschinger J.// Mittheilungen aus dem mechanish technichen Laboratorium der Technichen Hochschule in Munchen.-1886,1887.-HV.

107. Bryan G.H. Proc. Camb. Phil. Soc., VI (1888).

108. Considéré A. G. Resistance des piences comprimées // Congres International des procédés de construction.-T. 3.-Paris: 1891.-P.371.

109. Duleau A. Essai the orique et experimential sur la resistance du fer forge.-Paris: 1820.

110. Engesser F. Ueber Knickfestigkeit gerader Staebe // Zeitschrift des Architekten und Ingeniuer Verein zu Hannover. 1889. Bd. 35- S.456-468.

111. Engesser F. Ueber Knick Flagen // Schweizerishe Bauzetung.- 1895.-Bd.26 -S. 24.

112. HodgkinsonE. Experimental research on the strength of pillars of castjron and other materials // Philosophical Transactions 1840.- Part II.- P. 385-456.

113. Jasinski F. Zu den Knickfragen // Schweizerische Bauzeitung . 1895.-T. 25. -S. 172.

114. Karman T. Ueber Knickfestigkeit H Forschungsarbeiten. 1910. N81.

115. L.Tetmajer, Die Gesetze der Knickungs — und zusammengesetzten Festigkeit der technisch wichtigsten Baustoffe, Leipzig, 1907.

116. Lamarle E. Memo ire sur la flexion du bois // Annales des travaux publics de Belgique. Brussels.- 1845. T. 3; - 1846. -T. 4.

117. Musschenbrock P. van. Physique expérimentales et geometricae. Lugdum : 1729.

118. Rankin W. J. M. Manual of applied machanics . London : 1858.

119. Shanley F. R. Inelastic column theory // Journal of tne Aeronautical Science .-1946.-V. 13.-N12.

120. Southwell Phil Trans. Roy. Soc., London, ser A, v.213 (1913).

121. Tetmajer L. Mittheiliuigen der Anstalt zur Pruefung von Baumaterialen in Zuerich . 1890.-H. IV.129. www.ans.com.ru.130. www.mysopromat.ru.