
- •Содержание
- •Введение
- •Технические требования
- •Требования к компьютеру
- •Рекомендуемые параметры вычислительной техники для эффективной (профессиональной) работы с T-FLEX Анализом
- •Установка системы T-FLEX Анализ
- •Структурная организация приложения T-FLEX Анализ
- •Этапы анализа конструкций
- •Быстрое начало
- •Шаг 1. Подготовка объёмной твердотельной модели изделия
- •Шаг 2. Создание «Задачи»
- •Шаг 3. Назначение материала
- •Шаг 4.1 Наложение граничных условий. Задание закреплений
- •Шаг 4.2 Наложение граничных условий. Задание нагружений
- •Шаг 5. Выполнение расчёта
- •Шаг 6. Анализ результатов расчёта
- •Подготовка конечно-элементной модели для Анализа (Препроцессор)
- •Виды конечно-элементных моделей
- •Назначение и роль сеток
- •Виды и роль граничных условий
- •Управление «Задачами», команды управления задачами
- •Общие свойства задач
- •Задание материала
- •Построение сетки
- •Параметры сетки
- •Задание ограничений
- •Полное закрепление
- •Частичное закрепление
- •Контакт
- •Упругое основание
- •Задание нагрузок
- •Механические нагрузки
- •Сила
- •Давление
- •Вращение
- •Ускорение
- •Цилиндрическая нагрузка
- •Крутящий момент
- •Осциллятор
- •Дополнительная масса
- •Тепловые нагрузки
- •Температура
- •Тепловой поток
- •Тепловая мощность
- •Конвективный теплообмен
- •Излучение
- •Сводная таблица нагрузок
- •Редактирование нагрузок и закреплений
- •Настройки и сервисные команды
- •Работа с 3D окном при подготовке элементов задач
- •Особенности работы с параметрической моделью
- •Экспорт
- •Обработка результатов (Постпроцессор)
- •Общие принципы работы с результатами
- •Настройки и сервисные команды окна результатов расчёта
- •Настройка окна результатов расчёта
- •Настройка цветовой шкалы
- •Использование датчиков для анализа результатов
- •Использование графиков для анализа результатов
- •Интегральное значение
- •Построение сечений
- •Генерация отчётов
- •Пример интерпретации результата
- •Статический анализ
- •Особенности этапов статического анализа
- •Алгоритм оценки статической прочности по результатам моделирования
- •Настройки процессора линейной и нелинейной статики
- •Задача оптимизации
- •Задача об оптимизации толщины балки
- •Приложение (справочные материалы)
- •Характеристики конструкционных материалов
- •Объёмное напряжённо-деформированное состояние в точке
- •Оценка статической прочности конструкций. Теории прочности
- •Анализ устойчивости
- •Особенности этапов анализа на устойчивость
- •Алгоритм оценки устойчивости по результатам моделирования
- •Настройки Процессора анализа устойчивости
- •Частотный анализ
- •Особенности этапов частотного анализа
- •Настройки Процессора частотного анализа
- •Вынужденные колебания
- •Вводные сведения
- •Особенности этапов анализа вынужденных колебаний
- •Настройки препроцессора анализа вынужденных колебаний
- •Настройки процессора анализа вынужденных колебаний
- •Настройки постпроцессора и анализ результатов вынужденных колебаний
- •Анализ Усталости
- •Цикл напряжений. Основные характеристики
- •Кривая усталости
- •Методы коррекции напряжений
- •Оценка характеристик сопротивления усталости при сложном напряженном состоянии
- •Этапы анализа усталости
- •Результаты усталостного расчёта
- •Примеры расчётов деталей на усталостную прочность
- •Однособытийный усталостный расчет
- •Многособытийный усталостный расчет
- •Примеры результатов однособытийного усталостного расчёта
- •Пример результатов многособытийного усталостного расчёта
- •Тепловой анализ
- •Особенности этапов теплового анализа
- •Настройки Процессора теплового анализа
- •Примеры тепловых расчётов
- •Тепловой расчёт радиатора охлаждения. Установившийся режим
- •Расчёт времени нагревания радиатора охлаждения. Нестационарный режим
- •Расчёт времени остывания радиатора охлаждения. Нестационарный режим
- •Верификационные примеры
- •Примеры расчётов задач статики
- •Изгиб консольно-защемлённой балки под действием сосредоточенной нагрузки
- •Статический расчет круглой пластины, защемленной по контуру
- •Расчет сферического сосуда давления
- •Квадратная пластина под силой в центре
- •Цилиндрический резервуар со стенками постоянной толщины
- •Кручение бруса с круглым поперечным сечением
- •Стержень под действием собственного веса
- •Расчёт вращающегося сплошного диска постоянной толщины
- •Свободно опертая прямоугольная пластинка под синусоидальной нагрузкой
- •Температурные напряжения биметаллического элемента
- •Примеры расчётов задач устойчивости
- •Расчет устойчивости сжатого прямого стержня
- •Устойчивость квадратной пластины
- •Устойчивость прямоугольной пластины
- •Примеры задач частотного анализа
- •Определение собственных частот колебаний балки
- •Определение первой собственной частоты колебаний круглой пластинки
- •Свободные колебания сферического купола
- •Примеры задач теплового анализа
- •Установившаяся температура
- •Поток тепла в шаре
- •Теплопроводность цилиндрической стенки
- •Литература

Руководство пользователя T-FLEX Анализ
Модуль перемещения (при z = 0.5l ) u = 2.9139E-005 м.
Таблица 2. Результат «Угол закручивания»
Численное решение |
|
|
|
ϕ |
−ψ |
|
|
|
|
|||
Угол закручивания |
Аналитическое решение |
Ошибка δ = |
|
|
|
|
×100% |
|||||
|
|
|||||||||||
|
u |
|
Угол закручивания ϕ , рад |
|
|
|
|
|
|
|
||
рад |
|
|
ϕ |
|
|
|
|
|||||
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|||||
ψ = arcsin |
, |
|
|
|
|
|
|
|
||||
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
1.4570E-003 |
|
1.4551E-003 |
1.2608E-001 |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Стержень под действием собственного веса
Рассмотрим стержень радиуса R и длиной l , подвешенный за верхний конец и растягиваемый собственным весом (см. рис.).
Примем следующие исходные данные: длина стержня l равна 1м, радиус поперечного сечения стержняR равен 0.02м.
Характеристики материала: E = 2.1×1011 Па, ν = 0.28, ρ = 7800 мкг3 .
Полное удлинение стержня под действием собственного веса можно определить по формуле [7]:
l = γ l 2 ,
2E
где γ – удельный вес материала стержня, то есть γ = ρ g , g = 9.80665 см2 .
Напряжение в сечении стержня, находящемся на расстоянии x от нижнего (незакрёпленного) конца, можно вычислить по формуле [7]:
σ = γ x .
Таким образом, l =1.8212×10−7 м; σ = 3.8246×104 |
Н |
|
при |
x = |
l |
. |
|
|
|
|
|
|
|
|||||
м2 |
2 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Выполнив расчет при помощи T-FLEX Анализ, получаем следующие результаты: |
|
|
|
|
Таблица 1. |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Параметры конечно-элементной сетки |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
Тип конечных элементов |
Число главных |
|
Число узлов при |
|
Число конечных |
|||||||||||||
узлов |
|
|
расчёте задачи |
|
|
элементов |
||||||||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
квадратичный тетраэдр (10 узлов) |
2482 |
|
|
|
|
|
16089 |
|
|
9518 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2. |
|
|
|
|
|
|
|
|
|
|
|
Результат «Перемещение» |
|||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
Численное решение |
|
Аналитическое решение |
|
Ошибка δ = |
w − w* |
|
×100% |
|||||||||||
Перемещение w |
* |
, м |
|
Перемещение w , м |
|
|
|
|
|
|
||||||||
|
|
|
|
w |
|
|
||||||||||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
1.8177E-007 |
|
|
|
1.8212E-007 |
|
|
|
|
|
1.9410E-001 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
188

Верификационные примеры
Таблица 3. Результат «Напряжение»
Численное решение |
|
Аналитическое решение |
|
|
|
σ −σ* |
|
|
|
|||||||
|
|
|
|
|
||||||||||||
Напряжение σ |
* |
, |
Н |
|
Напряжение σ , |
Н |
|
Ошибка δ = |
|
|
|
|
|
|
|
×100% |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
σ |
|
|
|
|||||||||
|
|
|
|
|
|
|||||||||||
|
м2 |
|
м2 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
8.0139E-003 |
|
||||||||||
3.8249E+004 |
|
|
3.8246E+004 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Расчёт вращающегося сплошного диска постоянной толщины
Рассмотрим диск радиуса R и толщиной h (см. рис.).
В расчёте рассматривается 14 часть диска с наложением условий симметрии на соответствующих
гранях (ограничение перемещений в направлении оси локальной системы координат, нормальной к плоскости грани).
Диск находится под действием центробежной силы F = ρ ω2 R , где ρ - масса единицы объёма материала диска, ω - угловая скорость вращения.b
Примем следующие исходные данные: радиус диска R равен 0.457м , толщина диска h равна
0.01м, величина угловой скорости вращения ω равна 300 радс .
Характеристики материала: E = 2.1×1011 Па, ν = 0.28, ρ = 7800 мкг3 .
189

Руководство пользователя T-FLEX Анализ
Перемещения u в таком случае можно определить по формуле [8]:
|
|
1 |
|
|
|
1 ( |
1−ν |
2 |
) |
2 |
3 |
|
|
|
|
|
|
|
(1−ν )C1r −(1−ν )C2 r − |
|
|
|
|
, |
|||
|
|
|
|
|
|
|
|
|
|||||
|
|
u = E |
8 |
|
|
ρω r |
|
||||||
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где постоянные C = |
3 +ν |
ρω2 R2 |
, C |
2 |
= 0 определяются из граничных условий. |
||||||||
|
|||||||||||||
1 |
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Максимальное перемещение umax |
ожидается при r = R , то есть umax = |
(1−ν )ρω2 R3 . |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
4E |
Компоненты напряжений σr , σθ находятся в виде [8]:
σr = 3 +8ν ρω2 (R2 − r 2 ),
σθ = 3 +8ν ρω2 R2 −1+83ν ρω2r2 .
Эти напряжения принимают максимальное значение в центре диска, где
|
|
|
|
|
|
|
σr |
|
=σθ |
= |
3 +ν |
ρω2 R2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, umax = 5.7430×10−5 м, σ |
|
|
= σθ |
=σ = 6.0111×107 |
|
Н |
. |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
м2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
Выполнив расчет при помощи T-FLEX Анализ, получаем следующие результаты: |
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Параметры конечно-элементной сетки |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Тип конечных элементов |
|
|
Число главных |
|
|
Число узлов при |
|
Число конечных |
|||||||||||||||||||||||
|
|
|
узлов |
|
|
|
расчёте задачи |
|
|
|
|
|
элементов |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
квадратичный тетраэдр (10 узлов) |
|
|
|
1553 |
|
|
|
|
9010 |
|
|
|
|
4372 |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Результат «Перемещение» |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
Численное решение |
|
Аналитическое решение |
|
|
Ошибка δ = |
|
|
|
w − w* |
|
|
×100% |
|||||||||||||||||||
Перемещение w |
* |
, м |
|
|
|
|
Перемещение w , м |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
w |
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
5.8033E-005 |
|
|
|
|
|
|
|
5.7430E-005 |
|
|
|
|
|
|
1.0390E+000 |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 3. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Результат «Напряжение» |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
Численное решение |
|
Аналитическое решение |
|
|
|
|
|
|
|
|
σ −σ* |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
Напряжение σ |
* |
, |
|
Н |
|
|
|
|
|
Напряжение |
σ , |
Н |
|
|
|
Ошибка δ = |
|
|
|
|
|
|
|
|
×100% |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
σ |
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
м2 |
|
|
|
|
м2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
6.4740E+007 |
|
|
|
|
|
|
|
6.0111E+007 |
|
|
|
|
|
|
7.7009E+00 |
|
190

Верификационные примеры
Свободно опертая прямоугольная пластинка под синусоидальной нагрузкой
Рассмотрим прямоугольную пластинку со сторонами a , b и толщиной h (см. рис.).
Толщина пластины h существенно меньше длин её сторон a , b .
Пластина находится под действием нагрузки, распределённой по поверхности пластины согласно
закону: q = q |
sin |
|
π x sin |
|
π y |
, где q0 представляет собой интенсивность нагрузки в центре |
||
0 |
|
|
a |
|
|
b |
|
|
|
|
|
|
|
|
|
пластинки.
Рассмотрим случай, когда края пластины свободно оперты.
Примем следующие исходные данные: длина стороны a пластины равна 0.5м, длина стороны b пластины равна 0.4м, толщина пластины h = 0.003м , интенсивность нагрузки в центре пластинки
q0 =100 мН2 .
Характеристики материала: E = 2.1×1011 Па, ν = 0.28. Аналитическое решение этой задачи имеет вид [5]:
|
|
|
w = |
|
|
|
|
q0 |
|
|
|
|
|
|
|
|
|
π x |
|
|
|
π y |
, |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sin |
|
|
|
|
|
sin |
|
|
|
|
|
||||||
|
|
|
|
|
1 |
|
1 |
|
2 |
a |
|
|
|
|
|||||||||||||||||||
|
|
|
|
π |
4 |
D |
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
2 |
b |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
где D = |
Eh3 |
|
– цилиндрическая жесткость пластины. |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
12(1−ν 2 ) |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
Максимальный прогиб получается в центре пластинки. |
|
Подставив |
x = |
a |
и y = |
b |
в уравнение, |
||||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||
получаем: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
2 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
q0 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
wmax |
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
4 |
|
1 |
|
|
|
1 2 |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
π |
|
D |
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
b |
2 |
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
||||
Таким образом, |
w =1.9059×10−5 |
м. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
max |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Врасчёте рассматривается 14 часть пластины с наложением условий симметрии на
соответствующих гранях (ограничение перемещений в направлении оси локальной системы координат, нормальной к плоскости грани; ограничение поворотов).
191

Руководство пользователя T-FLEX Анализ
|
|
|
Модель пластины с нагрузками и закреплениями |
|
|
|
|
|
|
||||||||
Выполнив расчет при помощи T-FLEX Анализ, получаем следующие результаты. |
|
|
|
Таблица 1. |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
Параметры конечно-элементных сеток |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Номер |
|
|
|
|
|
Число узлов при |
|
|
|
Число |
||||||
|
|
Тип конечных элементов |
|
|
Число узлов |
|
|
|
конечных |
|
|||||||
|
сетки |
|
|
|
|
расчёте задачи |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
элементов |
|
||||
|
1 |
квадратичный тетраэдр (10 узлов) |
|
5836 |
|
|
34698 |
|
|
|
|
17318 |
|
||||
|
2 |
|
линейный треугольник (6 узлов) |
|
2747 |
|
|
2747 |
|
|
|
|
5304 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
квадратичный треугольник(6 узлов) |
2747 |
|
|
10797 |
|
|
|
|
5304 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2. |
||
|
|
|
|
|
|
|
|
|
Результат «Перемещение» |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
Номер |
|
Численное решение |
Аналитическое решение |
|
|
|
w − w* |
|
|
|||||||
|
сетки |
|
Перемещение w* , м |
|
Перемещение w , м |
|
Ошибка δ = |
|
|
|
|
|
×100% |
|
|||
|
|
|
|
|
w |
|
|
|
|||||||||
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||
|
1 |
|
1.8741E-005 |
|
|
1.9059E-005 |
|
|
1.6660E+000 |
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||
|
2 |
|
1.8798E-005 |
|
|
1.9059E-005 |
|
|
1.3693E+000 |
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||
|
3 |
|
1.8752E-005 |
|
|
1.9059E-005 |
|
|
1.6098E+000 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
192