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

В данном примере рассматривается статическое нагружение сборки из двух отсеков ракеты - носителя. Нижний отсек выполнен в виде трехслойной оболочки с заполнителем, верхний переходный отсек - в виде фермы, на верхнем торце которой размещен Г- образный металлический шпангоут. На шпангоуте закреплена масса 500 кг, имитирующая полезную нагрузку. Трехслойная оболочка имеет обшивки из углепластика со схемой армирования [20°2мм/-20°2мм]. Толщина заполнителя 80 мм. Стержни ферменной конструкции выполнены также из углепластика с такой же схемой армирования. Стержни объединены с трехслойной оболочкой и массой жесткой связью по всем степеням свободы.

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

Импортирование геометрической модели

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

Геометрическая модель ракеты-носителя

Задайте для задачи необходимые настройки. Для этого в командной строке пропишите вручную:

  • merge surface all. Нажмите Enter;

  • merge curve all. Нажмите Enter.

  • merge vertex all. Нажмите Enter

Данные команды выполняют слияние указанных объектов.

Создайте новую систему координат. На панели команд выберите Режим - Системы координат, Системы координат - Создать. Задайте необходимые параметры и нажмите Применить.

Создайте вершину. На панели команд выберите Режим - Геометрия, Объект - Вершина, Действие - Создать. Из выплывающего списка выберите Расположение и задайте необходимые параметры. Нажмите Применить.

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

Создайте новые системы координат. На панели команд выберите Режим - Системы координат, Системы координат - Создать. Задайте необходимые параметры и нажмите Применить.

Проделайте тоже самое для вершин: (91, 92, 89); (155, 156, 153); (98, 97, 100); (163, 164, 161); (42, 41, 43); (50, 108, 105); (50, 49, 51); (115, 116, 114); (58, 57, 59); (123, 124, 122); (66, 65, 67); (131, 132, 129); (74, 73, 75); (82, 140, 137); (82, 81, 84); (147, 148, 146).

Построение сетки

Создайте сетку на поверхностях. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Интервалы. Из выплывающего списка выберите Автоматический размер и задайте необходимые параметры. Нажмите Задать размер, Построить сетку.

Постройте сетку на других поверхностях 16 19 18 21 20 23 22 25 24 27 26 29 28 31 30 33 32 35 34 37 11 36 39 38 13 41 12 15 40 14 42 17.

Постройте сетку на кривых. На панели команд выберите Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выплывающего списка выберите Равномерно и задайте необходимые параметры. Нажмите Применить, Построить сетку.

Постройте сетку на вершине. На панели команд выберите Режим - Сетка, Объект - Узловая, Действие - Построение сетки. Из выплывающего списка выберите Построение сетки и задайте необходимые параметры. Нажмите Построить сетку.

Задание материала

Создайте материал. На панели команд выберите модуль задания свойств материала Режим - Материал, Объект - Управление материалами. Перетащите из правой колонки надпись Кевлар, а также из раздела Общие надпись Плотность в колонку Свойства материала. Задайте необходимые параметры.

Создайте второй материал Заполнитель. Перетащите из левой колонки надпись Материал Гука, а так же в разделе Общие надпись Плотность в колонку Свойства материала. Задайте необходимые параметры.

Перетащите из правой колонки надпись Углеродистая сталь.

Нажмите Применить и закройте окно Управление материалами.

Создание блоков и задание свойств

Создайте первый блок. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте нужные параметры. Нажмите Применить.

Задайте параметры для созданного блока. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Задайте свойства для блока. Затем необходимо задать свойства оболочке, для этого нажмите на кнопку с троеточием. Введите необходимые параметры. Нажмите Применить.

Закройте окно Задать свойства оболочки и нажмите Применить.

Создайте второй блок. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте нужные параметры. Нажмите Применить.

Создайте еще несколько блоков, для этого проделайте тоже самое для ID объекта(ов): 39 40, 25 26, 41 42, 11 12, 27 28, 13 14, 29 30, 16 15, 32 31, 17 18, 33 34, 19 20, 35 36, 21 22, 37 38.

Задайте свойства оболочки. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Нажмите на кнопку с троеточием. Введите необходимые параметр и нажмите Применить. Закройте окно Задать свойства оболочки.

Для второго блока задайте необходимые параметры и нажмите Применить.

Далее для блоков 3-17 задайте необходимые параметры. Порядок, Категория и ID свойств блока остаются прежними, а система координат (СИ) меняется. Блок 3 - СИ 5, блок 4 - СИ 6, блок 5 - СИ 7, блок 6 - СИ 8, блок 7 - СИ 9, блок 8 - СИ 10, блок 9 - СИ 11, блок 10 - СИ 12, блок 11 - СИ 13, блок 12 - СИ 14, блок 13 - СИ 15, блок 14 - СИ 16, блок 15 - СИ 17, блок 16 - СИ 18, блок 17 - СИ 19.

Создайте еще один блок. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте нужные параметры. Нажмите Применить.

Задайте свойства балки. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Задайте свойства для блока. Затем необходимо задать свойства балки, для этого нажмите на кнопку с троеточием. Введите необходимые параметры. Нажмите Применить.

Закройте окно Задать свойства балки и нажмите Применить.

Создайте последний блок. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте нужные параметры. Нажмите Применить.

Задайте свойства точечной массы. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Задайте свойства для блока. Затем необходимо задать свойства точечной массы, для этого нажмите на кнопку с троеточием. Введите необходимые параметры. Нажмите Применить.

Закройте окно Задать свойства точечной массы и нажмите Применить.

Задание граничных условий

Добавьте связи. На панели команд выберите Режим - Граничные условия, Объект - Связи, Действие - Создать. Задайте нужные параметры. Нажмите Применить.

Проделайте тоже самое для таких параметров: 175 и 62 124 126 60; 177 и 68 70 132 134; 179 и 140 142 76 78; 181 и 84 86 148 150; 170 и 156 158 92 94; 169 и 166 100 102 164; 171 и 46 108 110 44; 22 и 114 120 64 58 ;13 и 72 128 122 66;16 и 136 74 130 80; 28 и 82 138 88 144; 17 и 152 90 146 96; 9 и 98 154 104 160; 3 и 168 48 42 162; 8 и 112 50 106 56.

Задайте нужные параметры. Нажмите Применить.

Закрепите часть ракеты по всем осям. На панели команд выберите Режим - Граничные условия, Объект - Перемещение, Действие - Создать. Задайте нужные параметры. Нажмите Применить.

Задайте точечную силу. На панели команд выберите Режим- Граничные условия, Объект - Точечная сила, Действие - Создать. Задайте требуемые параметры. Нажмите Применить.

Задайте силу тяжести. На панели команд выберите Режим - Граничные условия, Объект - Гравитация, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Запуск расчета

Для запуска на расчет в CAE Fidesys выберите на панели команд Режим - Настройки расчета, Настройки расчета - Потеря устойчивости, Потеря устойчивости - Общие. Задайте требуемые параметры и нажмите Применить, Начать расчет.

Анализ результатов

Откройте файл с результатами. Это можно сделать тремя способами:

  • Нажмите Ctrl+E;

  • В главном меню выберите Расчёт - Результаты. Нажмите Открыть последний результат;

  • На панели команд выберите Результаты (Режим - Результаты, Результаты - Открыть Результаты).

Появится окно Fidesys Viewer, в котором вы сможете ознакомиться с результатами расчёта.

Можно посмотреть, как деформируется тело под приложенной силой. Для этого выберите Фильтр → Алфавитный указатель → Деформировать по вектору. Во вкладке Свойства установите необходимые параметры и нажмите Применить. В результате отобразится деформированное тело.

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

Выберите Фильтр → Алфавитный указатель → Порог. Во вкладке Свойства установите необходимые параметры и нажмите Применить.

Затем примените еще один фильтр. Выберите Фильтр → Алфавитный указатель → Системы координат. Во вкладке Свойства установите необходимые параметры и нажмите Применить. В результате отобразятся нормальные напряжения вдоль волокон в слое 1 верхней обшивки трехслойной оболочки.

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

Перейдите на фильтр Порог, поменяйте параметры во вкладке Свойства и нажмите Применить.

Вернитесь на фильтр Системы координат и во вкладке Свойства установите необходимые параметры и нажмите Применить. В результате отобразятся нормальные напряжения вдоль волокон в слое 1 стержня фермы.

Во вкладке Свойства установите другие параметры и нажмите Применить. В результате отобразятся нормальные напряжения вдоль волокон в слое 2 стержня фермы.

Удалите все ранее созданные фильтры, выберите Фильтр → Алфавитный указатель → Деформировать по вектору. Задайте необходимые параметры и нажмите Применить. В результате отобразится деформированный вид при статическом нагружении.

Использование консольного интерфейса

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

reset
import iges 'C:/Users/Rocket_otsek2.igs' nofreesurfaces heal
merge surface all
merge curve all
merge vertex all
create cs type cartesian origin vertex 4 dir1 vertex 3 dir2 vertex 1
create vertex location 0 0 0 
create vertex center curve 169 
create cs type cartesian origin vertex 185 dir1 vertex 186 dir2 vertex 49
create cs type cartesian origin vertex 91 dir1 vertex 92 dir2 vertex 89
create cs type cartesian origin vertex 155 dir1 vertex 156 dir2 vertex 153
create cs type cartesian origin vertex 98 dir1 vertex 97 dir2 vertex 100
create cs type cartesian origin vertex 163 dir1 vertex 164 dir2 vertex 161
create cs type cartesian origin vertex 42 dir1 vertex 41 dir2 vertex 43
create cs type cartesian origin vertex 50 dir1 vertex 108 dir2 vertex 105
create cs type cartesian origin vertex 50 dir1 vertex 49 dir2 vertex 51
create cs type cartesian origin vertex 115 dir1 vertex 116 dir2 vertex 114
create cs type cartesian origin vertex 58 dir1 vertex 57 dir2 vertex 59
create cs type cartesian origin vertex 123 dir1 vertex 124 dir2 vertex 122
create cs type cartesian origin vertex 66 dir1 vertex 65 dir2 vertex 67
create cs type cartesian origin vertex 131 dir1 vertex 132 dir2 vertex 129
create cs type cartesian origin vertex 74 dir1 vertex 73 dir2 vertex 75
create cs type cartesian origin vertex 82 dir1 vertex 140 dir2 vertex 137
create cs type cartesian origin vertex 82 dir1 vertex 81 dir2 vertex 84
create cs type cartesian origin vertex 147 dir1 vertex 148 dir2 vertex 146
surface 7 8 9 6 3 4 10 2 5 1 size auto factor 5
undo group begin
surface 7 8 9 6 3 4 10 2 5 1 size auto factor 5
mesh surface 7 8 9 6 3 4 10 2 5 1
undo group end
surface 16 19 18 21 20 23 22 25 24 27 26 29 28 31 30 33 32 35 34 37 11 36 39 38 13 41 12 15 40 14 42 17 size auto factor 5
undo group begin
surface 16 19 18 21 20 23 22 25 24 27 26 29 28 31 30 33 32 35 34 37 11 36 39 38 13 41 12 15 40 14 42 17 size auto factor 5
mesh surface 16 19 18 21 20 23 22 25 24 27 26 29 28 31 30 33 32 35 34 37 11 36 39 38 13 41 12 15 40 14 42 17
undo group end
curve 169 171 170 173 172 175 174 176 interval 10
curve 169 171 170 173 172 175 174 176 scheme equal
curve 169 171 170 173 172 175 174 176 interval 10
curve 169 171 170 173 172 175 174 176 scheme equal
mesh curve 169 171 170 173 172 175 174 176
mesh vertex 185
create material 1 from 'Кевлар'
modify material 1 set property 'ORTHOTROPIC_E_X' value 2e+11
modify material 1 set property 'DENSITY' value 1500
create material 2
modify material 2 name 'Заполнитель'
modify material 2 set property 'DENSITY' value 100
modify material 2 set property 'POISSON' value 0.3
modify material 2 set property 'MODULUS' value 100000
create material 3 from 'Углеродистая сталь'
set duplicate block elements off
block 1 add surface 7 8 9 6 3 4 10 2 5 1
create shell properties 1
modify shell properties 1 layer count 5
modify shell properties 1 layer 1 thickness 0.002
modify shell properties 1 layer 1 material 1
modify shell properties 1 layer 1 angle -20
modify shell properties 1 layer 2 thickness 0.002
modify shell properties 1 layer 2 material 1
modify shell properties 1 layer 2 angle 20
modify shell properties 1 layer 3 thickness 0.08
modify shell properties 1 layer 3 material 2
modify shell properties 1 layer 3 angle 0
modify shell properties 1 layer 4 thickness 0.002
modify shell properties 1 layer 4 material 1
modify shell properties 1 layer 4 angle -20
modify shell properties 1 layer 5 thickness 0.002
modify shell properties 1 layer 5 material 1
modify shell properties 1 layer 5 angle 20
modify shell properties 1 eccentricity 0.5
modify shell properties 1 layer direction normal
modify shell properties 1 thickness_change on
block 1 cs 2 element shell order 1
block 1 shell properties 1
set duplicate block elements off
block 2 add surface 23 24
set duplicate block elements off
block 3 add surface 39 40
set duplicate block elements off
block 4 add surface 25 26
set duplicate block elements off
block 5 add surface 41 42
set duplicate block elements off
block 6 add surface 11 12
set duplicate block elements off
block 7 add surface 27 28
set duplicate block elements off
block 8 add surface 13 14
set duplicate block elements off
block 9 add surface 29 30
set duplicate block elements off
block 10 add surface 16 15
set duplicate block elements off
block 11 add surface 32 31
set duplicate block elements off
block 12 add surface 17 18
set duplicate block elements off
block 13 add surface 33 34
set duplicate block elements off
block 14 add surface 19 20
set duplicate block elements off
block 15 add surface 35 36
set duplicate block elements off
block 16 add surface 21 22
set duplicate block elements off
block 17 add surface 37 38
create shell properties 2
modify shell properties 2 layer count 2
modify shell properties 2 layer 1 thickness 0.002
modify shell properties 2 layer 1 material 1
modify shell properties 2 layer 1 angle 0
modify shell properties 2 layer 2 thickness 0.002
modify shell properties 2 layer 2 material 1
modify shell properties 2 layer 2 angle 0
modify shell properties 2 eccentricity 0.5
modify shell properties 2 layer direction normal
modify shell properties 2 thickness_change on
block 2 cs 4 element shell order 1
block 2 shell properties 2
block 3 cs 5 element shell order 1
block 3 shell properties 2
block 4 cs 6 element shell order 1
block 4 shell properties 2
block 5 cs 7 element shell order 1
block 5 shell properties 2
block 6 cs 8 element shell order 1
block 6 shell properties 2
block 7 cs 9 element shell order 1
block 7 shell properties 2
block 8 cs 10 element shell order 1
block 8 shell properties 2
block 9 cs 11 element shell order 1
block 9 shell properties 2
block 10 cs 12 element shell order 1
block 10 shell properties 2
block 11 cs 13 element shell order 1
block 11 shell properties 2
block 12 cs 14 element shell order 1
block 12 shell properties 2
block 13 cs 15 element shell order 1
block 13 shell properties 2
block 14 cs 16 element shell order 1
block 14 shell properties 2
block 15 cs 17 element shell order 1
block 15 shell properties 2
block 16 cs 18 element shell order 1
block 16 shell properties 2
block 17 cs 19 element shell order 1
block 17 shell properties 2
set duplicate block elements off
block 18 add curve 169 171 170 173 172 175 174 176
create beam properties 1
modify beam properties 1 type 'L-Beam'
modify beam properties 1 angle 0.0
modify beam properties 1 ey 0.0
modify beam properties 1 ez 0.0
modify beam properties 1 geom_B 0.055
modify beam properties 1 geom_H 0.1
modify beam properties 1 geom_c1 0.0072
modify beam properties 1 geom_d 0.0072
modify beam properties 1 mesh_quality 2
modify beam properties 1 warping_dof off
block 18 material 3 cs 1 element beam order 1
block 18 beam properties 1
set duplicate block elements off
block 19 add vertex 185
create lumpmass properties 1
modify lumpmass properties 1 mass 500
block 19 cs 1 element lumpmass
block 19 lumpmass properties 1
create coupling master vertex 173 slave curve 52 54 116 118 distance
create coupling master vertex 175 slave curve 62 124 126 60 distance
create coupling master vertex 177 slave curve 68 70 132 134 distance
create coupling master vertex 179 slave curve 140 142 76 78 distance
create coupling master vertex 181 slave curve 84 86 148 150 distance
create coupling master vertex 170 slave curve 156 158 92 94 distance
create coupling master vertex 169 slave curve 166 100 102 164 distance
create coupling master vertex 171 slave curve 46 108 110 44 distance
create coupling master vertex 22 slave curve 114 120 64 58 distance
create coupling master vertex 13 slave curve 72 128 122 66 distance
create coupling master vertex 16 slave curve 136 74 130 80 distance
create coupling master vertex 28 slave curve 82 138 88 144 distance
create coupling master vertex 17 slave curve 152 90 146 96 distance
create coupling master vertex 9 slave curve  98 154 104 160 distance
create coupling master vertex 3 slave curve 168 48 42 162, distance
create coupling master vertex 8 slave curve 112 50 106 56 distance
create coupling master vertex 173 175 177 179 181 170 169 171 slave vertex 185 distance
create displacement  on curve 38 32 36 24 26 18 14 6 4 10 dof all fix  
create force  on vertex 185 vector 1e4 -1e5 0 0 0 0 
create gravity global
modify gravity 1 dof 1 value 20
modify gravity 1 dof 2 value -20
analysis type buckling elasticity dim3
eigenvalue find 1 smallest