Сейсмология является важной составной частью наук о Земле. С одной стороны это наука о землетрясениях, с другой - это главный источник информации о глубинном строении земных недр, недоступных для прямого изучения. Сейсмические волны возбуждаются ветром, морскими волнами, падениями метеоритов, запусками ракет, вулканическими извержениями, атмосферными взрывами и даже людьми, проходящими близко от сейсмометров. Сейсмические наблюдения приоткрыли тайну внутреннего строения Луны и Марса. Для широкого внедрения в сейсмологию количественных подходов потребовалось существенно развить ее теоретическую базу, и прежде всего теорию распространения сейсмических волн. Этой цели послужили многочисленные исследования сейсмологов-теоретиков по решению так называемых прямых задач сейсмологии - расчету волновых полей в заданной модели среды при заданном способе возбуждения колебаний. Решим задачу Стокса для изучения распространения сейсмических волн.
Создайте геометрическую модель
Создайте параллелепипед. На панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb731822e841ef8f42a_img_1.png)
Переместите модель. На панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Преобразовать. Из выплывающего списка выберите Переместить Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb7b86f3b93abf97855_img_2.png)
Разделите одно ребро на две равные части. На панели команд выберите Режим - Геометрия, Объект - Кривая, Действие - Изменить. Из выплывающего списка выберите Разбиение. Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb791da53363c72a41c_img_3.png)
Построение сетки
Постройте гексаэдральную сетку. На панели команд выберите Режим - Сетка, Объект - Объемная, Действие - Интервалы. Из выплывающего списка выберите Задать размер. Задайте необходимые параметры. Нажмите Применить, Построить сетку.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb79bceec55c52a343f_img_4.png)
Задание материала и свойств блока
Создайте материал. На панели команд выберите Режим - Материал, Объект - Управление материалами. В окне Материал укажите имя "Материал 1". Перетащите нужные свойства из левого окна свойств. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb7d87d485433c106d5_img_5.png)
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb752e9f86ff3327342_img_6.png)
Создайте блок для объемной модели. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb868416f3f4e2e8560_img_7.png)
Задайте свойства блоку. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb8f92642128af0fe01_img_8.png)
Задание граничных условий
Закрепите сторону по направлению OY. На панели команд выберите Режим - Граничные условия, Объект - Перемещения, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb893a74d153e2927e3_img_9.png)
Закрепите ребро модели по направлению OX. На панели команд выберите Режим - Граничные условия, Объект - Перемещения, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb826b57240a3bb7163_img_10.png)
Закрепите сторону по направлению OZ. На панели команд выберите Режим - Граничные условия, Объект - Перемещения, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb931822e7cd0f8f43d_img_11.png)
Приложите точечную силу к вершине модели. На панели команд выберите Режим - Граничные условия, Объект - Точечная сила, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fb8bd721e29a380dc03_img_12.png)
Задайте зависимость для точечной силы по закону Берлаге. На панели команд выберите Режим - Зависимость ГУ. Задайте необходимые параметры для силы. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fba72023b0051b6c255_img_13.png)
Приложите к четырем сторонам куба неотражающие граничные условия. На панели команд выберите Режим - Граничные условия, Объект - Неотражающее условие, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fba68416f7d1f2e856f_img_14.png)
Запуск расчета
Для запуска на расчет в CAE Fidesys выберите на панели команд Режим - Настройки расчета, Настройки расчета - Временной анализ, Временной анализ - Общие. Задайте требуемые параметры и нажмите Применить.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fba48031456b2980a25_img_15.png)
Далее выберите на панели команд Режим - Настройки расчета, Настройки расчета - Временной анализ, Временной анализ - Поля вывода. Задайте требуемые параметры и нажмите Применить, Начать расчет.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fba313db858d9354a64_img_16.png)
Анализ результатов
Откройте файл с результатами. Это можно сделать тремя способами:
- Нажмите Ctrl+E;
- В главном меню выберите Расчёт - Результаты. Нажмите Открыть последний результат;
- На панели команд выберите Результаты (Режим - Результаты, Результаты - Открыть Результаты).
Появится окно Fidesys Viewer, в котором вы сможете ознакомиться с результатами расчёта.
На панели инструментов выберите необходимое время.
Для построения временной зависимости выберите любую точку на поверхности куба при помощи инструмента Выбрать точки вкл.. Далее на панели инструментов выберите Фильтры - Алфавитный указатель - Построить выделенное в зависимости от времени. Нажмите Применить.
Ниже на рисунках показаны распределение напряжений по Мизесу.
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fba46a5a50e326ff600_img_17.png)
![Fidesys: Распространение сейсмических волн. Задача Стокса](https://uploads-ssl.webflow.com/577f3315340bfadb0a80534f/627d0fbaab294a7b3c8180be_img_18.png)
Использование консольного интерфейса
Построение геометрии, генерацию сетки, задание граничных условий и материалов можно выполнить с использованием консольного интерфейса. Ниже приведён код программы, позволяющий выполнить шаги описанного выше руководства, необходимо только самостоятельно указать полный путь и название сохраняемого файла.
reset
brick x 100
move Volume 1 location 0 50 50 include_merged
partition create curve 6 position 0 50 50
volume all size 10
volume all size 10
mesh volume all
create material 1
modify material 1 name 'материал 1'
modify material 1 set property 'MODULUS' value 2e+08
modify material 1 set property 'POISSON' value 0.3
modify material 1 set property 'DENSITY' value 1900
set duplicate block elements off
block 1 add volume 1
block 1 material 1 cs 1 element solid order 3
create displacement on surface 3 dof 2 fix
create displacement on curve 2 dof 1 fix
create displacement on surface 2 dof 3 fix
create force on vertex 9 force value 1 direction 1 0 0
bcdep force 1 value 'berlage(2.5e+07, 10, time)'
create absorption on surface 1 5 6 4
analysis type dynamic elasticity dim3 preload off
dynamic method full_solution scheme explicit maxtime 0.4 maxsteps 50000
output nodalforce off energy off record3d on material off results everystep 10 fullperiodic off
calculation start path 'D:/Lavrova/result.pvd'