- •Руководство пользователя
- •10.4. Оценка погрешности вычисления напряжений 134
- •1.2. Описание конструкции
- •1.2.1. Системы координат
- •1.2.2. Геометрия модели
- •1.2.3. Конечные элементы
- •1.2.4. Нагрузки
- •1.2.5. Граничные условия
- •1.2.6. Свойства материалов
- •1.3. Структура входного файла msc/nastran
- •1.3.1. Установки nastraNa
- •1.3.2. Секция управления файлами
- •1.3.3. Секция управления выполнением задания
- •1.3.4. Секция управления расчетными случаями
- •1.3.5. Секция исходных данных
- •1.4. Пример модели в msc/nastran
- •1.5. Файлы, создаваемые в процессе решения
- •2. Положения, принятые в msc/nastran
- •2.1. Единицы измерений
- •2.2. Запись символов, целых и вещественных чисел
- •2.3. Свободный, малый и большой форматы полей
- •2.3.1. Малый формат
- •2.3.2. Свободный формат
- •2.3.3. Большой формат
- •2.4. Сетки разбиений и переходные сетки
- •2.5. Создание модели
- •2.6. Использование тестовых моделей
- •3. Ввод координат
- •3.1. Узловые точки
- •3.2. Скалярные точки
- •3.3. Системы координат
- •3.3.1. Ввод ортогональной системы координат
- •3.3.2. Ввод цилиндрической системы координат
- •3.3.3. Ввод сферической системы координат
- •3.3.4. Системы координат элемента и материала
- •4. Элементы в msc/nastran
- •4.1. Одномерные элементы
- •4.2. Двухмерные элементы
- •4.2.4. Сдвиговой плоский элемент cshear
- •4.2.5. Двухмерный элемент с трещиной crac2d
- •4.3. Трёхмерные элементы
- •4.3.1. Снеха, cpentAи ctetra
- •4.3.3. Трехмерный элемент с трещиной crac3d
- •4.4. Скалярные элементы
- •4.5. Ввод данных с помощью элемента genel
- •5. Ввод свойств материалов
- •5.1. Изотропный материал мат1
- •5.2. Двунаправленный анизотропный материал мат2
- •5.3. Осесимметричный ортотропный материал матз
- •5.4. Двунаправленный ортотропный материал mat8
- •5.5. Материал с пространственной анизотропией мат9
- •5.6. Карта ввода свойств оболочки pshell
- •5.7. Элемент из композиционного материала рсомр
- •6. Статические нагрузки
- •6.1. Задание нагрузок в узлах
- •6.2. Нагрузки, распределенные на одномерных элементах
- •6.3. Нагрузки, распределенные на поверхностях
- •6.4. Гравитационная и центробежная сила (grav, rforce)
- •6.4.1. Определение массовых характеристик модели
- •6.5. Предварительный натяг
- •6.6. Комбинирование нагрузок
- •6.7. Температурные нагрузки
- •6.7.1. Использование оператора subcom
- •7. Граничные условия
- •7.1. Закрепления перемещений в отдельных узлах
- •7.2. Автоматическое закрепление в узлах (autospc)
- •7.3. Задание вынужденных перемещений в узловых точках (spcd, spc)
- •7.4. Связь перемещений в нескольких узлах
- •8. Жёсткие элементы
- •8.1. Описание жёстких элементов
- •8.2. Элемент rbar
- •8.3. Элемент rbe2
- •8.4. Элемент rbe3
- •9. Руководство по моделированию
- •9.1. Правильный выбор элемента
- •9.1.1. Общие положения
- •9.1.2. Точечные элементы
- •9.1.3. Одномерные элементы
- •9.1.4. Двухмерные элементы
- •9.1.5. Трёхмерные элементы
- •9.1.6. Жёсткие элементы
- •9.2. Частота сетки разбиения
- •9.3. Сетки в переходных зонах
- •9.3.1. Переход от крупной к мелкой сетке
- •9.3.2. Переходные зоны между элементами разных типов
- •9.4. Напряжения в узловых точках
- •9.4.1. Вычисление и вывод напряжений в узловых точках
- •9.4.2. Топологический метод
- •9.4.3. Геометрический метод
- •9.4.4. Напряжения в особых точках
- •Интерфейс пользователя
- •9.5. Правильно заданое нагружение
- •9.6. Симметрия
- •Интерфейс пользователя
- •10. Верификация модели
- •10.1. Использование возможностей графического препроцессора
- •10.1.1. Масштабирование элементов
- •10.1.7. Распределение свойств материалов и элементов
- •10.1.8. Изображение нормалей плоских элементов
- •10.2. Вывод энергии деформаций
- •10.3. Инструменты диагностики в nastraNe
- •10.3.1. Проверка элементов
- •10.3.2. Проверка геометрии элементов cquad4
- •10.3.3. Проверка геометрии элементов снеха
- •10.3.4. Определение центра тяжести и моментов инерции
- •10.3.5. Выявление механизмов и вырожденных степеней свободы в модели
- •10.3.6. Проверка заданных нагрузок.
- •10.3.7. Проверка сил реакций
- •10.3.8. Проверка конструкции нагружением "lg"
- •10.3.9. Проверка равновесием незакрепленной модели
- •10.3.10. Проверка равновесия от температурного нагружения
- •10.3.11. Равновесие сил в узловой точке
- •10.4. Оценка погрешности вычисления напряжений
- •10.4.1. "Скачки" напряжений в узловых точках
- •10.4.2. "Скачки" напряжений вдоль элементов
- •10.4.3. Обсуждение способов измерений погрешностей
- •Интерфейс пользователя
- •10.5. Проверки с использованием постпроцессора
- •11. Инерционное уравновешивание.
- •11.1. Описание применения инерционного уравновешивания
- •11.2. Задание инерционного уравновешивания вMsc/nastran
- •12. Матричные операции.
- •12.1 Определение набора
- •12.1.1 Глобальный набор перемещений
- •12.1.2. Поднаборы глобального набора перемещений.
- •12.2 Статическая конденсация (редуцирование по Гайяну)
- •Интерфейс пользователя
- •12.3 Прямой ввод матриц
- •Интерфейс Bulk Data для dmig
- •Интерфейс для раздела Case Control
- •13. Задача устойчивости в линейной постановке
- •13.1. Подход метода конечных элементов
- •13.2. Методы определения собственных значений
- •13.2.1. Метод обратных итераций (inverse power)
- •13.2.2. Усовершенствованный метод обратных итераций
- •13.2.3. Метод Ланцоша
- •13.2.4. Сравнение методов
- •Интерфейс пользователя
- •13.3. Допущения и ограничения расчета устойчивости при линейном подходе
- •14. Повторные запуски расчетов - restarts
- •14.1. Типы повторных запусков
- •Старые рестарты
- •Новые рестарты
- •14.2. Структура входного файла msc/nastran
- •Интерфейс пользователя
- •14.3.1. "Холодный" запуск расчёта
- •14.3.2. Повторный запуск расчёта
- •14.3 Определение версии рестарта
- •14.4 Прочая информация по рестартам
- •15. Обслуживание баз данных
- •15.1 Определения
- •15.2 База данных msc/nastran
- •15.3. Секция File Managemernt Statements
- •15.4 Рекомендации для решения задач большой размерности
- •Приложение Оценка ресурсов
- •Оценка размеров базы данных
4.3.3. Трехмерный элемент с трещиной crac3d
Трехмерный элемент с трещиной используется для моделирования объемных тел, имеющих нарушения сплошности в виде трещины. Как и CRAC2D элемент CRAC3D реализован как "dummy''-элемент, поэтому в Bulk Data необходимо задавать карту ADUM9. Топология элемента задается в карте CRAC3D. Форматы карт CRAC3D и ADUM9 следующие.
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
CRACD2 |
EID |
PID |
G1 |
G2 |
G3 |
G4 |
G5 |
G6 |
|
|
G7 |
G8 |
G9 |
G10 |
G11 |
G12 |
G13 |
G14 |
|
|
G15 |
G16 |
G17 |
G18 |
G19 |
G20 |
G21 |
G22 |
|
|
G23 |
G24 |
G25 |
G26 |
G27 |
G28 |
G29 |
G30 |
|
|
G31 |
G32 |
G33 |
G34 |
G35 |
G36 |
G37 |
G38 |
|
|
G39 |
G40 |
G41 |
G42 |
G43 |
G44 |
G45 |
G46 |
|
|
G47 |
G48 |
G49 |
G50 |
G51 |
G52 |
G53 |
G54 |
|
|
G55 |
G56 |
G57 |
G58 |
G59 |
G60 |
G61 |
G62 |
|
|
G63 |
G64 |
|
|
|
|
|
|
|
ADUM9 |
18 |
0 |
5 |
0 |
CRAC3D |
|
|
|
|
EID - идентификационный номер элемента (целое, > 0)
PID - идентификационный номер свойств - номер карты PSHEAR
Gi - идентификационные номера узловых угловых точек.
Задание узловых точек с 1 по 10 и с 19 по 28 обязательно, остальные точки могут не задаваться.
Свойства элемента задаются картой PRAC3D, имеющей следующий формат
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
PRAC3D |
PID |
MID |
GAMMA |
PH1 |
|
|
|
|
|
Figure 4-32. CRAC3D Solid Crack Tip Element.
Здесь
PID - идентификационный номер свойств
MID - идентификационный номер материала МАТ1
GAMMA - величина экспоненты, определяющая поле перемещений
PHI - угол (в градусах) относительно оси х элемента, вдоль которой определяется коэффициент интенсивности напряжений.
Можно использовать 2 опции - вариант бруска (brick) и симметричную опцию. В первом случае программа автоматически разбивает элемент на 8 пятигранников (wedges), во втором - на 4. Для опции "брусок" напряжения вычисляются путем осреднения напряжений в wedges 4 и 5, коэффициенты интенсивности напряжений Кi и Кii вычисляются по wedges 1 и 8. Для симметричной опции напряжения вычисляются в wedge 4, коэффициенты интенсивности - в wedge 1.
Элемент выполнен в трехмерной формулировке. Как стороны (с узловыми точками с 1 по 18 и с 19 по 36) так и срединная поверхность (точки с 37 по 46) должны быть плоскими, значительные отклонения могут вызвать печать сообщения об ошибке.
4.4. Скалярные элементы
Скалярный элемент есть элемент, связывающий две степени свободы, либо одну степень свободы и основание (жесткую плиту). Степень свободы может быть любой из 6 степеней свободы узловой точки или скалярной точки, имеющей 1 компонент. В отличие от одно- двух- и трехмерных элементов скалярный элемент не требует геометрического описания. Поэтому этот элемент не имеет и элементной координатной системы.
Скалярные элементы обычно используются совместно с обычными конечными элементами для моделирования тех участков, где физический объект неизвестен, либо его задание не обязательно. Типовые примеры таких объектов - это демпфер, упругая связь в соединении, изолирующая прокладка и т.п. Если используются скалярные элементы для связи узловых точек, то весьма рекомендуется, чтобы данные узловые точки совпадали. Если точки не совпадают, то силы, возникающие в узловой точке от действия скалярных элементов могут приводить к появлению значительных моментов, искажающих реальное состояние.
Для статических расчетов используется элементы типа линейной пружины (CELASi, I=1-4) и сосредоточенной массы (CMASSi, I=1-4). Существуют 4 типа карт определения пружины и массы как скалярных элементов. Форматы карт CELASi (упругие пружины) следующие
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
CELAS1 |
EID |
PID |
G1 |
C1 |
G2 |
C2 |
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
CELAS2 |
EID |
K |
G1 |
C1 |
G2 |
C2 |
CE |
S |
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
CELAS3 |
EID |
PID |
S1 |
S2 |
|
|
|
|
|
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
CELAS4 |
EID |
K |
S1 |
S2 |
|
|
|
|
|
где
EID - уникальный идентификационный номер материала
PID - идентификационный номер карты свойств - PELAS (для CELAS1 и CELAS3
G1,G2 - идентификационные номера узловых или скалярных точек (для CELAS1 и CELAS2)
С1 ,С2 - номера компонент (для CELAS1 и CELAS2)
S1,S2 - идентификационные номера скалярных точек (для CELAS3 и CELAS4)
К - жесткость скалярной пружины
S - коэффициент напряжений (CELAS2)
GE - коэффициент демпфирования (CELAS2)
Первый тип скалярных элементов (CELAS1, CMASS1) может ссылаться как на узловые, так и на скалярные точки. При описании топологии указываются идентификационные номера соответствующих элементов - аналогично заданию CROD. Во втором типе скалярных элементов (CELAS2, CMASS2) тоже может использовать узловые и скалярные точки, однако свойства элемента задаются непосредственно в карте CELAS. Третий и четвертый типы подобны соответственно первому и второму, за исключением того, что используют только скалярные точки. Если модель имеет много скалярных элементов, соединенных только со скалярными точками, то предпочтительно использовать эти последние типы элементов. Формат карты PELAS следующий
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
PELAS |
PID1 |
K1 |
CE1 |
S1PID2 |
K2 |
S2 |
|
|
|
PIDi - идентификационный номер свойств
Ki - величина, характеризующая упругость
Gei - коэффициент демпфирования
Si - коэффициент напряжений