
- •Общие положения
- •Универсальный стержень (КЭ 10)
- •Специальные конечные элементы (КЭ 51, 53,54,55)
- •РЕШЕНИЕ СИСТЕМЫ КАНОНИЧЕСКИХ УРАВНЕНИЙ
- •РАСЧЕТ НА ДИНАМИЧЕСКИЕ ВОЗДЕЙСТВИЯ
- •СУПЕРЭЛЕМЕНТНОЕ МОДЕЛИРОВАНИЕ
- •Стержни
- •Плоское напряженное состояние
- •Плиты
- •Оболочки
- •Объемные элементы
- •Загружения
- •РАСЧЕТ НА УСТОЙЧИВОСТЬ
- •РЕШЕНИЕ НЕЛИНЕЙНЫХ ЗАДАЧ
- •Общие положения
- •Расчет физически нелинейных задач
- •Библиотека законов деформирования материалов
- •Типы дробления сечений стержней
- •Типы арматурных включений
- •Библиотека конечных элементов для физически нелинейных задач
- •Стержневые конечные элементы (КЭ 210 и 205)
- •Универсальный стержневой элемент (КЭ - 310)
- •Конечный элемент предварительного натяжения (КЭ 308)
- •Специальные конечные элементы односторонних связей
- •Одноузловой элемент односторонней связи (тип КЭ-261)
- •Двухузловой элемент одностоpонней связи (тип КЭ - 262)
- •Рациональная разбивка на конечные элементы
- •Объединение перемещений
- •Абсолютно жесткие вставки
- •Угол чистого вращения
- •Моделирование податливости узлов сопряжения элементов
- •Моделирование шарниров в стержневых и плоскостных элементах
- •Расчет на заданные перемещения
- •Введение связей конечной жесткости
- •Расчет на температурные воздействия
- •Моделирование предварительного напряжения
- •Учёт прямой и косой симметрии
- •Вычисление коэффициентов постели упругого основания
- •Учет работы конструкций совместно с упругим основанием
- •Расчет оболочек и плит, подкреплённых рёбрами
- •Задание весов масс и динамических воздействий
- •Сбор нагрузок на фундаменты
- •Расчетные сочетания нагрузок
- •Согласованная система координат для пластин
- •ПРИНЦИПЫ АНАЛИЗА РЕЗУЛЬТАТОВ РАСЧЕТА
- •Правила знаков при чтении результатов расчета.
- •Результаты расчета на динамические воздействия
- •Суммарные усилия от динамических воздействий
- •Документирование
- •ЖЕСТКОСТНЫЕ ХАРАКТЕРИСТИКИ ЭЛЕМЕНТОВ
- •ПРОВЕРКА ПРОЧНОСТИ ПО РАЗЛИЧНЫМ ТЕОРИЯМ
- •Главные напряжения
- •КЭ плоской задачи теории упругости
- •КЭ плиты
- •КЭ объемного НДС
- •КЭ оболочки
- •Стержневые КЭ
- •Вычисление эквивалентных напряжений
- •Результаты расчета
- •РАСЧЕТ И ПРОЕКТИРОВАНИЕ СТАЛЬНЫХ КОНСТРУКЦИЙ
- •Назначение и возможности
- •Проектируемые сечения
- •Задание дополнительных данных для расчета
- •Конструктивные и унифицированные элементы
- •Проверки несущей способности элементов
- •Описание алгоритмов
- •Сквозной расчет
- •Локальный расчет
- •Представление результатов расчета
- •ПОДБОР И ПРОВЕРКА АРМИРОВАНИЯ В ЖЕЛЕЗОБЕТОННЫХ ЭЛЕМЕНТАХ
- •Армирование стержневых элементов
- •Проверка заданного армирования
- •Армирование пластинчатых элементов
6. РАСЧЕТ НА УСТОЙЧИВОСТЬ
Реализованный вариант расчета на устойчивость предполагает, что распределение сил Nо известно из линейного расчета. Требуется найти значение числового параметра λо такое, чтобы при силах (λо * Nо) произошла потеря устойчивости.
Задача определения критических сил и соответствующих им форм потери устойчивости допускает следующую вариационную формулировку: найти перемещение и ≠ 0 и число λо такие, что при всех допустимых
перемещениях v справедливо равенство: |
|
a(u,v)= λ 0d(u,v) |
(6.1) |
где d(u,v ) - возможная работа сил при заданном их распределении No. |
|
Пользуясь выражением (1.3) и обозначив D матрицу с |
элементами |
di, j = d(µi ,µ j ), получим из (6.1) задачу на собственные значения для матриц |
|
KX = λ 0 DX |
(6.2) |
Погрешность МКЭ в задаче устойчивости для критических сил и соответствующих им форм потери устойчивости пропорциональна hτ .
Решение производится методом половинного деления. Этот метод основан на том, что матрицаK(λ) = K − λ D положительно определена лишь
при λ<λ0. Отсутствие положительной определенности матрицы соответствует наличию отрицательных чисел на главной диагонали после исключений по методу Гаусса.
После определения с заданной точностью параметра λ0 форма потери устойчивости находится как собственный вектор матрицы K(λ0) методом
итерации подпространств, изложенным в п.3.
Расчет реализуется в упругой стадии. Значения сжимающих сил и напряжений в элементах схемы уже вычислены с помощью линейного процессора. При выполнении расчета на устойчивость предполагается, что эти значения выражены через критический параметр нагрузки:
Pкpi = λi Pi Nкpij = λi Nij
где
i- номер загружения;
j- номер элемента в схеме;
Pi- суммарная нагрузка в i-том загружении; Pкpi - критическая нагрузка в i-том загружении;
Nij - продольное усилие или напряжение в j-том элементе в i-том
загружении;
Nкpij - критическое продольное усилие в j-том элементе в i-том
загружении;
λi - параметр нагрузки (коэффициент запаса устойчивости).
В процессе счета для каждого загружения определяются первая форма потери устойчивости и соответствующий ей коэффициент запаса.
Допускается также производить проверку устойчивости по линейным комбинациям загружений (РСН).
Если в расчете схемы присутствуют динамические загружения, то
40

проверка устойчивости схемы для них может быть произведена только через комбинации загружений (РСН). Это связано с тем, что разложенные по формам колебаний результаты расчета на динамическое воздействие необходимо преобразовать в суммарные.
Впроцессе расчета общей устойчивости итерационным методом определяется значение λ такое, при котором хотя бы один элемент диагонали матрицы жесткости обращается в ноль. Если λi ≥1, то считается, что схема устойчива в данном загружении или при данной комбинации загружений.
Вкачестве исходных данных задаются ς начальный масштабный
множитель к продольным силам Ni ( по умолчанию ς=2) , а также точность вычислений (по умолчанию равна 0.01). Предполагается, что при λi>ς система абсолютно устойчива.
Врезультате вычисляются коэффициенты запаса устойчивости λi, первая форма потери устойчивости и коэффициенты свободной длины для стержневых элементов, исходя из общей устойчивости, по следующим формулам:
µyij |
= |
π |
EJyj |
(6.3) |
|
|
lj |
Nкрij |
|
µzij |
= |
π |
EJzj , где: |
|
|
|
lj |
Nкрij |
|
µyij, µzij – коэффициенты свободной длины j-того стержня |
|
|||
|
|
соответственно в плоскостях X1oZ1, X1oY1 |
|
|
|
|
для i-того загружения; |
|
|
EJyj, EJzj– изгибные жесткости j-того стержня соответственно в |
|
|||
|
|
плоскостях соответственно X1oZ1, X1oY1; |
|
|
Nкрij = λi*Nij |
- критическое продольное усилие сжатия в j-том стержне |
для i-того загружения;
λi – коэффициент запаса устойчивости для i-того загружения; lj - длина j-того стержня.
41