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

Верификационные примеры
Эквивалентные напряжения определяются по формуле:
σэкв(ρ) = (σρ (ρ))2 + (σt (ρ))2 − 2 σρ (ρ) σt (ρ) .
Эквивалентные напряжения на внутренней поверхности сферы равны σэкв(r) =1148 МПа. Выполнив расчет при помощи T-FLEX Анализ, получаем следующие результаты:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 1. |
||
|
|
|
|
|
|
|
|
|
|
|
Параметры конечно-элементной сетки |
|||||||||||||||
|
|
|
|
|
|
|
Число конечных элементов |
|
|
|||||||||||||||||
Тип конечных элементов |
|
Число узлов |
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
5338 |
|
|
|
|
|
|
|
|
|
|||||||||
квадратичный тетраэдр (10 узлов) |
|
1846 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2. |
||
|
|
|
|
|
|
|
|
|
|
|
|
Результат «Перемещение» |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
w − w* |
|
|
|
|
|
|
|
|||||||||
Численное решение |
|
Аналитическое решение |
|
Ошибка δ = |
|
|
|
|
|
|
|
×100% |
|
|
||||||||||||
Перемещение w |
* |
, м |
|
|
Перемещение w , м |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
w |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
0.03 |
|
|
|
|
|
|
|
|
|
|||||||
1.4067E-003 |
|
|
|
|
|
1.4063E-003 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 3. |
||
|
|
|
|
|
|
|
|
|
|
Результат «Эквивалентные напряжения» |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
Численное решение |
|
Аналитическое решение |
|
|
|
|
|
σ −σ* |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
Напряжение σ |
* |
, |
|
Н |
|
|
Напряжение σ , |
Н |
|
|
Ошибка δ = |
|
|
|
|
|
|
|
|
|
×100% |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
σ |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
м2 |
|
м2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
0.17 |
|
|
|
|
|
|
|
|
|
|||||||||
1.1499E+009 |
|
|
|
|
|
1.1480E+009 |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Квадратная пластина под силой в центре
Рассматривается жестко закреплённая квадратная пластина под сосредоточенной силой в центре [4].
Задача статического анализа решается с использованием пластинчатых конечных элементов. В расчёте рассматривается 14 часть пластины с
наложением условий симметрии на соответствующих гранях (ограничение перемещений в направлении оси локальной системы координат, нормальной к плоскости грани; ограничение поворотов).
183

Руководство пользователя T-FLEX Анализ
Конечно-элементная модель пластины с нагрузками и закреплениями
Решение этой задачи для прогиба под силой даёт формула: w = 0.0224 |
Pa |
2 |
Eh3 |
|
|
|
, D = |
|
. |
||
D |
12(1−ν 2 ) |
||||
|
|
|
Примем следующие исходные данные: длина и ширина пластины a = 500мм, толщина пластины h = 3мм, приложенная сосредоточенная сила P = 50кгс (или P = 490.3325Н ).
Характеристики материала примем E = 2.1×1011 Па, ν = 0.28. Аналитическое решение имеет вид: w = 5.3557×10-3 м.
Выполнив расчет при помощи T-FLEX Анализ, получили следующие результаты:
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 1. |
|
|
|
|
Параметры конечно-элементных сеток |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Номер |
Тип конечных элементов |
|
Число узлов |
|
Число конечных |
||||||||
сетки |
|
|
элементов |
||||||||||
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
квадратичный тетраэдр (10 узлов) |
|
8083 |
|
|
24222 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
линейный треугольник (6 узлов) |
|
3818 |
|
|
7426 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
квадратичный треугольник(6 узлов) |
|
3818 |
|
|
7426 |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2. |
|
|
|
|
|
|
Результат «Перемещение» |
|||||||
|
|
|
|
|
|
|
|
|
|||||
Номер |
Численное решение |
Аналитическое решение |
Ошибка δ = |
|
w − w* |
|
×100% |
||||||
сетки |
Перемещение w* , м |
Перемещение w , м |
|
|
|
|
|
|
|||||
|
|
w |
|
|
|
||||||||
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
1 |
5.3700E-003 |
|
5.3557E-003 |
|
|
0.27 |
|
|
|
|
|||
2 |
5.4078E-003 |
|
5.3557E-003 |
|
|
0.97 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|||
3 |
5.3763E-003 |
|
5.3557E-003 |
|
|
0.38 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Цилиндрический резервуар со стенками постоянной толщины
Резервуар подвергается действию давления жидкости, как показано на рисунке. Днище резервуара защемлёно в абсолютно жёсткий фундамент.
184

Верификационные примеры
На практике толщина стенки резервуара h бывает в большинстве случаев мала в сравнении как с радиусом R , так и с глубиной d резервуара. Учитывая это и тот факт, что днище резервуара не испытывает деформаций, в качестве конечно-элементной модели резервуара можно рассматривать цилиндрическую оболочку, нижний край которой защемлён.
Конечно-элементная модель конструкции с нагрузками и закреплениями
Аналитическое решение этой задачи имеет вид [5]:
w = e−β zˆ (C1 cos(βzˆ)+C2 sin(βzˆ))− ρg(d − zˆ)R2 , Eh
185

Руководство пользователя T-FLEX Анализ
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3(1−ν 2 ) |
|
|
|
|
|
|
|
|||||
|
|
|
z = z + d , |
C1 = |
|
ρgR2d |
, |
|
ρgR2 |
|
1 |
|
, |
β |
4 |
= |
, |
|
|
|
|
|
|
||||||||||
|
|
|
|
C2 = |
|
|
d |
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
ˆ |
|
|
Eh |
Eh |
|
β |
|
|
2 |
h |
2 |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
|
|||||||
ρ - плотность жидкости, |
g - ускорение свободного падения ( ≈ 9.8 |
м |
) |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
с2 |
|
|
|
|
|
|
|
|
|
|
||
Примем |
следующие исходные |
данные: |
глубина |
резервуара |
d =1000мм, радиус |
R = 200мм , |
|||||||||||||||||||||||||||
толщина резервуара h = 3мм, плотность жидкости ρ =1000 кг |
м |
3 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Характеристики материала примем: E = 2.1×1011 Па, ν = 0.28. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
Таким образом, w = 6.1366×10-7 |
|
м (при zˆ = 0.056 м). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
Выполнив расчет при помощи T-FLEX Анализ, получили следующие результаты: |
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 1. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Параметры конечно-элементных сеток |
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Номер |
|
|
Тип конечных элементов |
|
|
|
Число узлов |
|
|
|
Число конечных |
||||||||||||||||||||||
сетки |
|
|
|
|
|
|
|
|
элементов |
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
1 |
|
|
квадратичный тетраэдр (10 узлов) |
|
|
|
|
23092 |
|
|
|
|
|
|
|
70280 |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
2 |
|
|
линейный треугольник (6 узлов) |
|
|
|
|
|
11676 |
|
|
|
|
|
|
|
23140 |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
3 |
|
|
квадратичный треугольник(6 узлов) |
|
|
|
|
11676 |
|
|
|
|
|
|
|
23140 |
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 2. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Результат «Перемещение» |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
Номер |
|
|
Численное решение |
|
Аналитическое решение |
|
|
Ошибка δ = |
|
w − w* |
|
×100% |
|||||||||||||||||||||
сетки |
|
|
Перемещение w* , м |
|
Перемещение w , м |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
w |
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
1 |
|
|
6.2345E-007 |
|
|
|
|
|
6.1366E-007 |
|
|
|
|
|
|
|
|
|
|
|
1.60 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
2 |
|
|
6.1938E-007 |
|
|
|
|
|
6.1366E-007 |
|
|
|
|
|
|
|
|
|
|
|
0.93 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
3 |
|
|
6.1743E-007 |
|
|
|
|
|
6.1366E-007 |
|
|
|
|
|
|
|
|
|
|
|
0.61 |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Кручение бруса с круглым поперечным сечением
Рассмотрим брус с круглым поперечным сечением радиуса R . Длина бруса l (см. рис.).
186

Верификационные примеры
Выберем систему координат с осью z вдоль оси бруса, при этом z = 0 на его левом конце.
Брус находится под действием внешнего скручивающего момента Mкр . Момент приложен на правом конце бруса, левый конец бруса жёстко закреплён.
Примем следующие исходные данные: длина l бруса равна 0.6м, радиус поперечного сечения R бруса равен 0.02м, величина приложенного крутящего момента M кр =100 Н м .
Характеристики материала: E = 2.1×1011 Па, ν = 0.28.
Для определения угла закручивания воспользуемся соотношением [6]:
|
|
z |
M |
êð |
|
|
|
|
|
|
|
|
ϕ = |
∫0 |
|
dz +ϕ0 |
, |
|
|
|
|||
|
|
|
|
|
|
||||||
|
|
GJ p |
|
|
|
|
|
||||
где ϕ0 |
– угол поворота сечения z = 0 , G = |
|
E |
|
|
- модуль сдвига, |
J p = |
πR4 |
- полярный момент |
||
2(1+ν) |
2 |
||||||||||
|
|
|
|
|
|
инерции круглого сечения.
Так как, по условию, левый конец бруса жёстко закреплён, то ϕ0 = 0 . Тогда на расстоянии z = 0.5l от жёстко закреплённого конца бруса угол закручивания ϕ определяется формулой:
|
|
ϕ0.5l = |
0.5Mкрl |
|
|
|
|
|
|
GJ p |
|
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, ϕ0.5l =1.4551×10−3 рад. |
|
|
|
|
|
|
|
Выполнив расчет при помощи T-FLEX Анализ, получаем следующие результаты. |
Таблица 1. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Параметры конечно-элементной сетки |
||
|
|
|
|
Число узлов при |
|
|
|
Тип конечных элементов |
|
Число главных |
|
|
Число конечных |
||
|
узлов |
|
расчёте задачи |
|
элементов |
||
|
|
|
|
||||
|
|
|
|
30188 |
|
|
|
квадратичный тетраэдр (10 узлов) |
|
4526 |
|
|
18345 |
187