fidesys
В данном примере демонстрируется возможность моделирования горного массива под действием давления и собственного веса, состоящего из трёх слоев различных материалов. При этом средний слой ослаблен горной выработкой (объемной искусственной полостью, сделанной в недрах земли или на поверхности.). Граничные условия задачи определяются условиями нагружения и условиями симметрии.
Построение геометрии
Создайте вершины, для этого на панели команд выберите Режим - Геометрия, Объект - Вершина, Действие - Создать. Из списка геометрических примитивов выберите Координаты и задайте необходимые координаты вершин. Нажмите Применить.
Проделайте тоже самое для координат (-901.50, 1.29, 0), (-902.00, 1.58, 0), (-902.50, 1.71, 0), (-904.00, 1.9, 0), (-904.50, 4, 0), (-905.00, 3.96, 0), (-905.50, 4.02, 0), (-906.00, 3.89, 0), (-906.50, 4.14, 0), (-907.50, 4.05, 0), (-908.10, 4.12, 0), (-908.50, 4.91, 0), (-909.00, 5.24, 0), (-909.90, 5.08, 0), (-911.00, 5.95, 0), (-912.10, 5.8, 0), (-913.10, 5.86, 0), (-914.10, 6.05, 0), (-915.10, 6.24, 0), (-916.10, 6.24, 0), (-917.10, 6.37, 0), (-918.10 ,6.53, 0), (-919.10, 11.81, 0), (-920.10, 13.65, 0), (-921.10, 15.21, 0), (-922.00, 16.19, 0), (-923.00, 17.07, 0), (-924.20, 17.71, 0), (-925.00, 18.32, 0), (-926, 18.96, 0), (-927.2, 19.41, 0), (-928.2, 19.91, 0), (-929, 20.2, 0), (-930, 20.93, 0), (-931, 21.5, 0), (-932.1, 22.41, 0), (-933.1, 23.35, 0), (-934.1, 24.13, 0), (-935.1, 24.92, 0), (-936.1, 25.65, 0), (-937.1, 26.22, 0), (-938.1, 26.93, 0), (-939.1, 27.45, 0), (-940.1, 27.92, 0), (-941.1, 28.56, 0), (-942.3, 29.07, 0), (-943.3, 29.63, 0), (-944.3, 30.16, 0), (-945.3, 30.67, 0), (-946.1, 31.16, 0), (-947.3, 31.76, 0), (-948.3, 32.13, 0), (-949.3, 32.58, 0), (-950.1, 32.86, 0), (-951.1, 33.15, 0), (-952.1, 32.96 ,0), (-953.1, 33.45, 0), (-954.1, 33.96, 0), (-955.1, 34.38, 0), (-956.1, 34.9, 0), (-957.1, 35.2, 0), (-958.1, 35.52, 0), (-959.1, 35.67, 0), (-960.1, 35.76, 0), (-961.1, 35.36, 0), (-962.4, 34.13, 0), (-963.4, 33.96, 0), (-964, 34.13, 0), (-965, 34.33, 0), (-966, 34.63 ,0), (-967, 33.5, 0), (-968.1, 34.08, 0), (-969.1, 34.77, 0), (-970.1, 35.28, 0), (-971.1, 6.2, 0), (-972.1, 37.03, 0), (-973.1, 37.72, 0), (-974.1, 38.32, 0), (-975.1, 38.83, 0), (-976.1, 39.38, 0), (-977.1, 39.65, 0), (-978.1, 39.97, 0), (-979.1, 40.25, 0), (-980.1, 40.43, 0), (-981.1, 40.56, 0), (-982.5, 40.51, 0), (-983.1, 40.48, 0), (-984.1, 40.21, 0), (-985.1, 39.97, 0), (-986.1, 39.78, 0), (-987.1, 39.56, 0), (-988.1, 39.46, 0), (-989.1, 39.26, 0), (-990.1, 39.18, 0), (-991.1, 39.17, 0), (-992.1, 39.05, 0), (-993.1, 39.05, 0), (-994.1, 39.06, 0), (-995.1, 38.97, 0), (-996.1, 38.81, 0), (-997.1, 38.9, 0), (-998, 39.06, 0), (-999.2, 39.3, 0), (-1000, 39.3, 0), (-1001, 39.35, 0), (-1004.1, 38.67, 0), (-1005.1, 38.35, 0), (-1006.1 ,37.96, 0), (-1007.1, 37.66 ,0), (-1008.3, 37.2, 0), (-1009.1, 36.86, 0), (-1010.1, 36.67, 0), (-1011.1, 34.46, 0), (-1012, 33.51, 0), (-1013, 32.98 ,0), (-1014, 32.91, 0), (-1015 ,32.26 ,0), (-1016, 31.92, 0), (-1017, 31.47, 0), (-1018, 31.27, 0), (-1019, 31.13, 0), (-1020, 30.1, 0), (-1021, 29.48, 0), (-1022.1, 28.2, 0), (-1023.3, 27.83, 0), (-1024.1, 25.85, 0), (-1025.1, 26.12, 0), (-1026.1, 26.41, 0), (-1027.1, 26.41, 0), (-1028.1, 26.3, 0), (-1029.1, 25.9, 0), (-1030.1, 25.61, 0), (-1031.1, 24.88, 0), (-1032.1, 24.31, 0), (-1033.1, 23.8, 0), (-1034.1, 23.45, 0), (-1035.1, 22.25, 0), (-1036.1, 22.03, 0), (-1037.1, 20.95, 0), (-1038.1, 19.46, 0), (-1039, 17.62, 0), (-1039, 0, 0).
Создайте сплайн. На панели команд выберите Режим - Геометрия, Объект - Кривая, Действие - Создать. Из выплывающего списка выберите Сплайн и укажите ID вершин: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142. Нажмите Применить.
Создайте кривую, из выпадающего списка выберите Линия. Введите необходимые значения и нажмите Применить.
Проделайте тоже самое для создания второй кривой. Введите ID вершин и нажмите Применить.
Создайте поверхность при помощи ограничивающих кривых. На панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Создать. Из выпадающего списка выберите Ограничивающие кривые и задайте необходимые параметры. Нажмите Применить.
Создайте прямоугольник. На панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Создать. Из выпадающего списка выберите Прямоугольник и задайте необходимые параметры. Нажмите Применить.
Переместите созданный прямоугольник, для этого на панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Преобразовать. Из выпадающего списка выберите Переместить и задайте необходимые параметры. Нажмите Применить.
Вычтите из второй поверхности первую. На панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Логические операции. Из выпадающего списка выберите Вычесть и задайте требуемые параметры. Нажмите Применить.
Разрежьте модель. На панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Разрез. Из списка возможных видов разреза выберите пункт Координатная плоскость. Задайте необходимые параметры. Нажмите Применить.
Еще раз разрежьте модель. В параметрах задайте ID тела 3 и значение смещения -1050. Нажмите Применить.
Создайте кривую. На панели команд выберите Режим - Геометрия, Объект - Кривая, Действие - Создать. Из выплывающего списка выберите Линия и задайте требуемые параметры. Нажмите Применить.
Создайте поверхность при помощи ограничивающих кривых. На панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Создать. Из выпадающего списка выберите Ограничивающие кривые и задайте необходимые параметры. Нажмите Применить.
Создайте объем при помощи вытягивания поверхности. Для этого на панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Создать. Из списка геометрических примитивов выберите Протяжка и задайте требуемые параметры. Нажмите Применить.
Разрежьте модель. На панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Разрез. Из списка возможных видов разреза выберите пункт Координатная плоскость. Задайте необходимые параметры. Нажмите Применить.
Еще раз разрежьте модель. В параметрах задайте ID тел 6 7 и значение смещения -918.1. Нажмите Применить.
Еще раз разрежьте модель. В параметрах задайте ID тела 8 9 и значение смещения -920.1. Нажмите Применить.
Срастите геометрическую модель. Введите в командную строку команду: merge all.
Построение сетки
Постройте сетку на кривых 48 77 106 37. На панели команд выберите Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выплывающего списка выберите Равномерно и укажите необходимые параметры. Нажмите Построить сетку.
Постройте сетку на кривых 41 34 104 79 50 28, задав при этом интервал равный 7.
Постройте сетку на кривых 116 118, задав при этом интервал равный 3.
Постройте сетку на кривых 14 27 22 42, задав при этом интервал равный 4.
Постройте сетку на кривых 54 55, задав при этом интервал равный 3.
Постройте сетку на поверхности 63. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Построение сетки. Из выплывающего списка выберите Карта и укажите необходимые параметры. Нажмите Построить сетку.
Постройте сетку на кривых 82 80. На панели команд выберите Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выплывающего списка выберите Равномерно и укажите необходимые параметры. Нажмите Построить сетку.
Постройте сетку на поверхности 46. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Построение сетки. Из выплывающего списка выберите Карта и укажите необходимые параметры. Нажмите Построить сетку.
Постройте сетку на кривых 117 119. На панели команд выберите Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выплывающего списка выберите Равномерно и укажите необходимые параметры. Нажмите Построить сетку.
Постройте сетку на поверхности 71. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Построение сетки. Из выплывающего списка выберите Карта и укажите необходимые параметры. Нажмите Построить сетку.
Постройте сетку на поверхности 16, из выпадающего списка выберите Автоматическая и укажите необходимые параметры. Нажмите Построить сетку.
Постройте сетку на кривых 18 39. На панели команд выберите Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выплывающих списков выберите Сгущение и Интервалы и сгущения укажите необходимые параметры. Нажмите Построить сетку.
Постройте сетку на кривой 105, задав при этом интервал равный 10, коэффициент сгущения - 1.3 и стартовая вершина - 25.
Постройте сетку на кривой 107, задав при этом интервал равный 15, коэффициент сгущения - 1.3 и стартовая вершина - 206.
Постройте сетку на кривой 76, задав при этом интервал равный 20, коэффициент сгущения - 1.3 и стартовая вершина - 23.
Постройте сетку на кривой 78, задав при этом интервал равный 20, коэффициент сгущения - 1.3 и стартовая вершина - 191.
Постройте тетраэдральную сетку на объеме 10. На панели команд выберите Режим - Сетка, Объект - Объемная, Действие - Построение сетки. Из выплывающего списка выберите Тетраэдральная и укажите необходимые параметры. Нажмите Построить сетку.
Постройте сетку на объеме 6.
Постройте сетку на объеме 8.
Постройте сетку на объеме 3.
Постройте сетку на кривой 13. На панели команд выберите Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выплывающих списков выберите Сгущение и Интервалы и сгущения укажите необходимые параметры. Нажмите Построить сетку
Постройте сетку на кривой 23, задав при этом интервал равный 7, коэффициент сгущения - 1.3 и стартовая вершина - 159.
Постройте сетку на кривых 30 7 40 5. На панели команд выберите Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выплывающего списка выберите Равномерно и укажите необходимые параметры. Нажмите Построить сетку.
Постройте тетраэдральную сетку на объемах 2 4. На панели команд выберите Режим - Сетка, Объект - Объемная, Действие - Построение сетки. Из выплывающего списка выберите Тетраэдральная и укажите необходимые параметры. Нажмите Построить сетку.
Усовершенствование модели
Удалите лишние объемы. На панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Удалить. Укажите ID объемов, которые нужно удалить. Нажмите Применить.
Удалите лишние вершины. На панели команд выберите Режим - Геометрия, Объект - Вершина, Действие - Удалить. Укажите ID вершин: 10 100 101 102 103 104 105 106 107 108 109 11 110 111 112 113 114 115 116 117 118 119 12 120 121 122 123 124 125 126 127 128 129 13 130 131 132 133 134 135 136 137 138 139 14 140 141 15 16 17 18 19 2 20 21 22 24 26 27 28 29 3 30 31 32 33 34 35 36 37 38 39 4 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 6 60 61 62 63 64 65 66 67 68 69 7 70 71 72 73 74 75 76 77 78 79 8 80 81 82 83 84 85 86 87 88 89 9 90 91 92 93 94 95 96 97 98 99. Нажмите Применить.
Поверните геометрическую модель вокруг оси Z. На панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Преобразовать. Из выплывающего списка выберите Повернуть и укажите необходимые параметры. Нажмите Применить.
Поверните геометрическую модель вокруг оси Y. На панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Преобразовать. Из выплывающего списка выберите Повернуть и укажите необходимые параметры. Нажмите Применить.
Задание материала и свойств блоков
Создайте три вида материала. На панели команд выберите Режим - Материал, Объект - Управление материалами. В открывшемся виджете Управление материалами в средней колонке укажите имя первого материала. В колонке свойств откройте список Упругость и перетащите Материал Гука в колонку Свойства материала, затем в левой колонке перейдите в раздел и перетащите Общие → Плотность. Задайте необходимые значения.
Создайте второй материал. В колонке свойств из раздела Упругость перетащите Материал Гука в колонку Свойства материала, затем в левой колонке перейдите в раздел и перетащите Общие → Плотность, Пластичность → Критерий Друкера-Прагера (Ползучесть). Задайте необходимые значения.
Создайте третий материал. В колонке свойств из раздела Упругость перетащите Материал Гука в колонку Свойства материала, затем в левой колонке перейдите в раздел и перетащите Общие → Плотность. Задайте необходимые параметры. Нажмите Применить и закройте окно Управление материалами.
Создайте первый блок. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте необходимые параметры. Нажмите Применить.
Создайте второй блок. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте необходимые параметры. Нажмите Применить.
Создайте третий блок. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте необходимые параметры. Нажмите Применить.
Задайте свойства первому блоку. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Задайте необходимые параметры. Нажмите Применить.
Задайте свойства второму блоку. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Задайте необходимые параметры. Нажмите Применить.
Задайте свойства третьему блоку. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Задайте необходимые параметры. Нажмите Применить.
Задание граничных условий
Закрепите модель в перемещениях по X. На панели команд выберите Режим - Граничные условия, Объект - Перемещение, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
Закрепите модель в перемещениях по X и Z. На панели команд выберите Режим - Граничные условия, Объект - Перемещение, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
Закрепите модель в перемещениях по Z. На панели команд выберите Режим - Граничные условия, Объект - Перемещение, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
Приложите давление на верхнюю грань. На панели команд выберите Режим - Граничные условия, Объект - Давление, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
Приложите давление на поверхность выработки. На панели команд выберите Режим - Граничные условия, Объект - Давление, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
Задайте гравитацию на весь объем. На панели команд выберите Режим - Граничные условия, Объект - Гравитация, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
Задание зависимостей ГУ
Задайте зависимость для давления. На панели команд выберите Режим - Зависимость ГУ. Нажмите на вкладку Давление 1 и укажите, что табличная зависимость будет от времени. Нажмите на 
чтобы получилось 3 строки. Введите необходимые значения. Нажмите Применить.
Затем перейдите на вкладку Давление 2 и укажите, что табличная зависимость будет от времени. Нажмите на чтобы получилось 3 строки. Введите необходимые значения. Нажмите Применить.
Перейдите на вкладку Гравитация 1, степень свободы Y и укажите, что табличная зависимость будет от времени. Нажмите на чтобы получилось 3 строки. Введите необходимые значения. Нажмите Применить. Закройте окно Зависимость ГУ.
Закрепите модель в перемещениях по Y. На панели команд выберите Режим - Граничные условия, Объект - Перемещение, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.
Запуск расчета
Для запуска на расчет в CAE Fidesys выберите на панели команд Режим - Настройки расчета, Настройки расчета - Статический, Действие - Общие. Задайте требуемые параметры и нажмите Применить, Начать расчет.
Анализ результатов
Откройте файл с результатами. Это можно сделать тремя способами:
Нажмите Ctrl+E;
В главном меню выберите Расчёт - Результаты. Нажмите Открыть последний результат;
На панели команд выберите Результаты (Режим - Результаты, Результаты - Открыть Результаты).
Появится окно Fidesys Viewer, в котором вы сможете ознакомиться с результатами расчёта.
На верхней панели выберите данные результата расчета для отображения. Из первого выпадающего списка выберите Перемещения, из второго – X. Затем выберите точки на модели, используя кнопку Выбрать точки и примените фильтр Построить выделенное в зависимости от времени.
В качестве поля для отображения на графике выберите Перемещения_Y (точка: 2066), Перемещения_X (точка: 3252). В результате в правой части экрана отобразятся необходимые графики зависимости интенсивности перемещений от времени:
Далее в поле отображения выберите Деформации_Мизес (точка: 2066), Деформации_Мизес (точка: 3252). В результате отобразятся графики зависимостей интенсивности деформаций от времени:
Использование консольного интерфейса
Построение геометрии, генерацию сетки, задание граничных условий и материалов можно выполнить с использованием консольного интерфейса. Ниже приведён код программы, позволяющий выполнить шаги описанного выше руководства, необходимо только самостоятельно указать полный путь и название сохраняемого файла.
reset create vertex -900.00 0 0 create vertex -901.50 1.29 0 create vertex -902.00 1.58 0 create vertex -902.50 1.71 0 create vertex -904.00 1.9 0 create vertex -904.50 4 0 create vertex -905.00 3.96 0 create vertex -905.50 4.02 0 create vertex -906.00 3.89 0 create vertex -906.50 4.14 0 create vertex -907.50 4.05 0 create vertex -908.10 4.12 0 create vertex -908.50 4.91 0 create vertex -909.00 5.24 0 create vertex -909.90 5.08 0 create vertex -911.00 5.95 0 create vertex -912.10 5.8 0 create vertex -913.10 5.86 0 create vertex -914.10 6.05 0 create vertex -915.10 6.24 0 create vertex -916.10 6.24 0 create vertex -917.10 6.37 0 create vertex -918.10 6.53 0 create vertex -919.10 11.81 0 create vertex -920.10 13.65 0 create vertex -921.10 15.21 0 create vertex -922.00 16.19 0 create vertex -923.00 17.07 0 create vertex -924.20 17.71 0 create vertex -925.00 18.32 0 create vertex -926 18.96 0 create vertex -927.2 19.41 0 create vertex -928.2 19.91 0 create vertex -929 20.2 0 create vertex -930 20.93 0 create vertex -931 21.5 0 create vertex -932.1 22.41 0 create vertex -933.1 23.35 0 create vertex -934.1 24.13 0 create vertex -935.1 24.92 0 create vertex -936.1 25.65 0 create vertex -937.1 26.22 0 create vertex -938.1 26.93 0 create vertex -939.1 27.45 0 create vertex -940.1 27.92 0 create vertex -941.1 28.56 0 create vertex -942.3 29.07 0 create vertex -943.3 29.63 0 create vertex -944.3 30.16 0 create vertex -945.3 30.67 0 create vertex -946.1 31.16 0 create vertex -947.3 31.76 0 create vertex -948.3 32.13 0 create vertex -949.3 32.58 0 create vertex -950.1 32.86 0 create vertex -951.1 33.15 0 create vertex -952.1 32.96 0 create vertex -953.1 33.45 0 create vertex -954.1 33.96 0 create vertex -955.1 34.38 0 create vertex -956.1 34.9 0 create vertex -957.1 35.2 0 create vertex -958.1 35.52 0 create vertex -959.1 35.67 0 create vertex -960.1 35.76 0 create vertex -961.1 35.36 0 create vertex -962.4 34.13 0 create vertex -963.4 33.96 0 create vertex -964 34.13 0 create vertex -965 34.33 0 create vertex -966 34.63 0 create vertex -967 33.5 0 create vertex -968.1 34.08 0 create vertex -969.1 34.77 0 create vertex -970.1 35.28 0 create vertex -971.1 36.2 0 create vertex -972.1 37.03 0 create vertex -973.1 37.72 0 create vertex -974.1 38.32 0 create vertex -975.1 38.83 0 create vertex -976.1 39.38 0 create vertex -977.1 39.65 0 create vertex -978.1 39.97 0 create vertex -979.1 40.25 0 create vertex -980.1 40.43 0 create vertex -981.1 40.56 0 create vertex -982.5 40.51 0 create vertex -983.1 40.48 0 create vertex -984.1 40.21 0 create vertex -985.1 39.97 0 create vertex -986.1 39.78 0 create vertex -987.1 39.56 0 create vertex -988.1 39.46 0 create vertex -989.1 39.26 0 create vertex -990.1 39.18 0 create vertex -991.1 39.17 0 create vertex -992.1 39.05 0 create vertex -993.1 39.05 0 create vertex -994.1 39.06 0 create vertex -995.1 38.97 0 create vertex -996.1 38.81 0 create vertex -997.1 38.9 0 create vertex -998 39.06 0 create vertex -999.2 39.3 0 create vertex -1000 39.3 0 create vertex -1001 39.35 0 create vertex -1004.1 38.67 0 create vertex -1005.1 38.35 0 create vertex -1006.1 37.96 0 create vertex -1007.1 37.66 0 create vertex -1008.3 37.2 0 create vertex -1009.1 36.86 0 create vertex -1010.1 36.67 0 create vertex -1011.1 34.46 0 create vertex -1012 33.51 0 create vertex -1013 32.98 0 create vertex -1014 32.91 0 create vertex -1015 32.26 0 create vertex -1016 31.92 0 create vertex -1017 31.47 0 create vertex -1018 31.27 0 create vertex -1019 31.13 0 create vertex -1020 30.1 0 create vertex -1021 29.48 0 create vertex -1022.1 28.2 0 create vertex -1023.3 27.83 0 create vertex -1024.1 25.85 0 create vertex -1025.1 26.12 0 create vertex -1026.1 26.41 0 create vertex -1027.1 26.41 0 create vertex -1028.1 26.3 0 create vertex -1029.1 25.9 0 create vertex -1030.1 25.61 0 create vertex -1031.1 24.88 0 create vertex -1032.1 24.31 0 create vertex -1033.1 23.8 0 create vertex -1034.1 23.45 0 create vertex -1035.1 22.25 0 create vertex -1036.1 22.03 0 create vertex -1037.1 20.95 0 create vertex -1038.1 19.46 0 create vertex -1039 17.62 0 create vertex -1039 0 0 create curve spline location vertex 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 create curve vertex 142 143 create curve vertex 1 143 create surface curve all create surface rectangle width 546 height 243 zplane move Surface 2 x -970 y 121.5 include_merged subtract surface 1 from surface 2 webcut body 2 with plane xplane offset -880 webcut body 3 with plane xplane offset -1050 create curve vertex 153 151 create surface curve 24 9 8 sweep surface all xaxis angle 90 merge webcut body 3 5 with plane xplane offset -904 webcut body 6 7 with plane xplane offset -918.1 webcut body 8 9 with plane xplane offset -920.1 merge all curve 48 77 106 37 interval 5 curve 48 77 106 37 scheme equal mesh curve 48 77 106 37 curve 41 34 104 79 50 28 interval 7 curve 41 34 104 79 50 28 scheme equal mesh curve 41 34 104 79 50 28 curve 116 118 interval 3 curve 116 118 scheme equal mesh curve 116 118 curve 14 27 22 42 interval 4 curve 14 27 22 42 scheme equal mesh curve 14 27 22 42 curve 54 55 interval 3 curve 54 55 scheme equal mesh curve 54 55 surface 29 scheme tripave mesh surface 29 curve 108 109 interval 7 curve 108 109 scheme equal mesh curve 108 109 surface 63 scheme map mesh surface 63 curve 82 80 interval 5 curve 82 80 scheme equal mesh curve 82 80 surface 46 scheme map mesh surface 46 curve 117 119 interval 20 curve 117 119 scheme equal mesh curve 117 119 surface 71 scheme map mesh surface 71 surface 16 scheme auto mesh surface 16 curve 18 39 interval 10 curve 18 39 scheme bias factor 1.2 start vertex 159 mesh curve 18 39 curve 105 107 interval 15 curve 105 scheme bias factor 1.3 start vertex 25 curve 107 scheme bias factor 1.3 start vertex 206 mesh curve 105 107 curve 76 78 interval 20 curve 76 scheme bias factor 1.3 start vertex 23 curve 78 scheme bias factor 1.3 start vertex 191 volume 10 scheme tetmesh mesh volume 10 volume 6 scheme tetmesh mesh volume 6 volume 8 scheme tetmesh mesh volume 8 volume 3 scheme tetmesh mesh volume 3 curve 13 23 interval 7 curve 13 scheme bias factor 1.3 start vertex 155 curve 23 scheme bias factor 1.3 start vertex 159 mesh curve 13 23 curve 30 7 40 5 interval 5 curve 30 7 40 5 scheme equal mesh curve 30 7 40 5 volume 2 4 scheme tetmesh mesh volume 2 4 delete volume 5 7 9 11 delete vertex 10 100 101 102 103 104 105 106 107 108 109 11 110 111 112 113 114 115 116 117 118 119 12 120 121 122 123 124 125 126 127 128 129 13 130 131 132 133 134 135 136 137 138 139 14 140 141 15 16 17 18 19 2 20 21 22 24 26 27 28 29 3 30 31 32 33 34 35 36 37 38 39 4 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 6 60 61 62 63 64 65 66 67 68 69 7 70 71 72 73 74 75 76 77 78 79 8 80 81 82 83 84 85 86 87 88 89 9 90 91 92 93 94 95 96 97 98 99 rotate Volume all angle 90 about Z include_merged rotate Volume all angle 90 about Y include_merged create material 1 modify material 1 name 'Izvestnyak' modify material 1 set property 'MODULUS' value 18.444e3 modify material 1 set property 'POISSON' value 0.17 modify material 1 set property 'DENSITY' value 0.0022 create material 2 modify material 2 name 'Salt' modify material 2 set property 'MODULUS' value 941 modify material 2 set property 'POISSON' value 0.3 modify material 2 set property 'DENSITY' value 0.00217 modify material 2 set property 'DPC_YIELD_STRENGTH' value 4e+20 modify material 2 set property 'DPC_INTERNAL_FRICTION_ANGLE' value 0 modify material 2 set property 'DPC_A' value 2.28e-08 modify material 2 set property 'DPC_N' value 3.34 modify material 2 set property 'DPC_M' value -0.82 create material 3 modify material 3 name 'Angidrit' modify material 3 set property 'MODULUS' value 46.511e3 modify material 3 set property 'POISSON' value 0.22 modify material 3 set property 'DENSITY' value 0.00291 block 1 add volume 2 set duplicate block elements off block 2 add volume 3 6 8 10 set duplicate block elements off block 3 add volume 4 block 1 material 1 cs 1 element solid order 1 block 2 material 2 cs 1 element solid order 1 block 3 material 3 cs 1 element solid order 1 create displacement on surface 12 30 47 64 70 22 dof 1 fix 0 create displacement on surface 20 69 48 27 9 dof 1 dof 3 fix 0 create displacement on surface 28 4 45 66 68 7 dof 3 fix 0 create pressure on surface 10 magnitude 1 create pressure on surface 71 63 46 29 16 magnitude 1 create gravity global modify gravity 1 dof 2 value 1 create table 1 modify table 1 dependency time modify table 1 insert row 1 modify table 1 insert row 1 modify table 1 insert row 1 modify table 1 cell 2 1 value 1 modify table 1 cell 2 2 value 15.042 modify table 1 cell 3 1 value 12 modify table 1 cell 3 2 value 15.042 bcdep pressure 1 table 1 create table 2 modify table 2 dependency time modify table 2 insert row 1 modify table 2 insert row 1 modify table 2 insert row 1 modify table 2 insert row 1 modify table 2 insert row 1 modify table 2 cell 2 1 value 1 modify table 2 cell 2 2 value 19 modify table 2 cell 3 1 value 2 modify table 2 cell 3 2 value 18 modify table 2 cell 4 1 value 9 modify table 2 cell 4 2 value 18 modify table 2 cell 5 1 value 12 modify table 2 cell 5 2 value 4 bcdep pressure 2 table 2 create table 3 modify table 3 dependency time modify table 3 insert row 1 modify table 3 insert row 1 modify table 3 insert row 1 modify table 3 cell 2 1 value 1 modify table 3 cell 2 2 value -9.81 modify table 3 cell 3 1 value 12 modify table 3 cell 3 2 value -9.81 modify gravity 1 dof 2 table 3 create displacement on surface 19 dof 2 fix 0 analysis type static plasticity elasticity dim3 static steps 12 nonlinearopts maxiters 50 minloadsteps 1 maxloadsteps 30 tolerance 1e-3 targetiter 5
