Fidesys: Распространение сейсмических волн. Задача Стокса

05/2022

Сейсмология является важной составной частью наук о Земле. С одной стороны это наука о землетрясениях, с другой - это главный источник информации о глубинном строении земных недр, недоступных для прямого изучения. Сейсмические волны возбуждаются ветром, морскими волнами, падениями метеоритов, запусками ракет, вулканическими извержениями, атмосферными взрывами и даже людьми, проходящими близко от сейсмометров. Сейсмические наблюдения приоткрыли тайну внутреннего строения Луны и Марса. Для широкого внедрения в сейсмологию количественных подходов потребовалось существенно развить ее теоретическую базу, и прежде всего теорию распространения сейсмических волн. Этой цели послужили многочисленные исследования сейсмологов-теоретиков по решению так называемых прямых задач сейсмологии - расчету волновых полей в заданной модели среды при заданном способе возбуждения колебаний. Решим задачу Стокса для изучения распространения сейсмических волн.

Создайте геометрическую модель

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

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

Fidesys: Распространение сейсмических волн. Задача Стокса
Fidesys: Распространение сейсмических волн. Задача Стокса

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

Fidesys: Распространение сейсмических волн. Задача Стокса

Задайте зависимость для точечной силы по закону Берлаге. На панели команд выберите Режим - Зависимость ГУ. Задайте необходимые параметры для силы. Нажмите Применить.

Fidesys: Распространение сейсмических волн. Задача Стокса

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

Fidesys: Распространение сейсмических волн. Задача Стокса

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

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

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

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

На панели инструментов выберите необходимое время.

Для построения временной зависимости выберите любую точку на поверхности куба при помощи инструмента Выбрать точки вкл.. Далее на панели инструментов выберите Фильтры - Алфавитный указатель - Построить выделенное в зависимости от времени. Нажмите Применить.

Ниже на рисунках показаны распределение напряжений по Мизесу.

Fidesys: Распространение сейсмических волн. Задача Стокса
Fidesys: Распространение сейсмических волн. Задача Стокса

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

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

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'