Timestep Control Parameters

SUBGRID DEFINITION

Subgrid Name: название сетки

- ввести название сетки

далее распахнется таблица 2.02:

Processor type - тип процессора

Типы процессоров:

1. Lagrange

2. Euler

3. ALE

4. Shell

5. Euler - FCT

6. SPH

просмотр типов процессора осуществляется клавишей 'пробел' по кругу. Выбрать тип процессора и нажать Enter.

далее распахнется таблица 2.03:

IMAX - количество ячеек по оси I (OX) от 1 до IMAX (например, 101, т.к.1 ячейка фиктивная)

JMAX - количество ячеек по оси J (OY) от 1 до JMAX (например, 101, т.к.1 ячейка фиктивная)

- ввести количество ячеек по осям.

Далее распахивается меню 3:

Меню 3

Zoning

Fill

Boundary

Edit/Targ

Options

Delete

Subgrid

Zoning - нанесение сетки

Fill - заливка материалами расчетной сетки

Boundary - геометрическое расположение граничных условий

Edit/Targ - расположение датчиков на расчетной сетке

Options - дополнительные функции (свойства)

Delete - удаление сетки в которой находимся

Subgrid - переход к другой сетке

- выбрать Zoning. Далее распахивается меню 4:

Меню 4

Generate

Transform

Copy

Delete

IJ-Range

Save

Restore

View

- выбрать Generate (генерирование сетки). Далее распахивается меню 5:

Меню 5

Predefs

Node

Edge

Blоck

Transitio

IJ-Range

View/Mod

- выбрать Predefs (готовые заготовки). Далее распахивается меню 6, характеризующее виды сеток:

Меню 6

Box

Circle

Ogive

Wedge

User

IJ-Range

View

Box - прямоугольная сетка

Circle - круглая сетка

Ogive - сетка в виде оживала

Wedge - грани

User - пользовательская сетка (сделанная самим)

IJ-Range - в каком диапазоне будем работать (настройка диапазона)

View - просмотр сетки

- выбрать Box. Далее распахивается таблица 6.01:

Box data:

XMIN -

XMAX -

YMIN -

YMAX -

физические координаты расчетной сетки

INC - инкримент при неравномерной сетке

Times -?

Pos -?

(выбрать по умолчанию)

Далее высвечиваются характеристики сгенерированной сетки.

В меню 6 при нажатии пункта View сгенерированная сетка отрисовывается. При этом возникает меню 6.1, с помощью которого можно детально рассмотреть сгенерированную сетку:

Меню 6.1

Zoom

Pan

Reset

Axes

Nodes

I-lines

J-lines

Examine

 Zoom - увеличение с помощью мыши

Pan - уменьшение всей картинки на 30%

Reset - возврат к общему виду всей сетки

Axes - отрисовка осей

Nodes - перемещать узлы -?

I -lines - переход к конкретной I-ой линии сетки

J-lines - переход к конкретной J-ой линии сетки

Examine - определение с помощью мыши характеристик точки на расчетной области (i,j и т.д.) -?

При нажатии несколько раз кнопки ESC возвращаемся в меню 7.2:

Меню 7.2

Global

Subgrid

Jain

Interact

Library

View

 

3. Задание материалов

- в меню 7.2 выберем пункт Global. Далее распахивается меню 8:

Меню 8

Symmetry

Timestep

Wrapup

Material

Boundary

Edit

Options

Heading

Symmetry - тип симметрии

Timestep - шаг по времени

Wrapup - количество циклов

Material - материалы

Boundary - тип граничных условий

Edit - задание протокола результатов расчетов на экране и в файлах

Options - дополнительные параметры (в том числе точка инициирования!)

Heading - редактирование заголовка файла

- в меню 8 выберем пункт Material. Далее распахивается меню 9:

 

Меню 9

Add/Mod

Delete

ChngeName

Review

Library

Options

 

- в меню 9 выберем пункт Library. Далее распахивается таблица 9.01:

Таблица 9.01

Current material library

Library Name:

в окошке последовательно появляются названия библиотек материалов:

ALLOYS, APPLIC, COMPDS, ELEMENT, EXAMPL, EXPLOS, GRANMT, HULOC, JNCOOK, JNHOLM, LEETAR, PUFF01, PHTCON, STEIN, SYNTHS, TILLOT,ZERARM,SESAME

 Клавишей 'пробел' можно найти из списка библиотек нужную и клавишей Enter выбрать ее (Выберем библиотеку EXAMPL). При этом распахивается меню 10:

Меню 10

Store

Retrieve

Delete

Print

Library

 

- выберем Retrieve (просмотр материалов, содержащихся в библиотеке). При этом распахивается таблица 10, содержащая названия материалов. Выберем нужный нам материал и клавишами 'Пробел' или Backspase установим YES напротив нужного пункта. После просмотра всех материалов библиотеки в результате на экране дисплея высвечиваются характеристики выбранного материала из данной библиотеки. На экране дисплея также появляется надпись:

Retrieve this material Yes No

Выберем нужное.

Повторным нажатием кнопки ESC вернемся в меню 7.2 или меню 8.

Пример:

В библиотеке EXAMPLE выберем материалы: воздух - GAIR и ВВ - GEXPLOSIVE. Модифицируем свойства ВВ. для этого выйдем в меню 9 и выберем пункт Add/Mod. На экране распахнется таблица 9.02:

MATERIAL

Material name:

в окошке последовательно появляются названия выбранных материалов (в нашем случае GAIR и GEXPLOSIVE). Клавишей 'пробел' найдем из этого списка материал GEXPLOSIVE и клавишей Enter выберем его. При этом распахивается таблица с характеристиками материала GEXPLOSIVE (таблица имеет продолжение на следующей странице экрана дисплея). Следует найти пункт:

  Auto. convert to Ideal Gas

и кнопкой 'пробел' поставить напротив этого пункта: Yes.

Конец примера.

4. Заливка областей расчетной сетки соответствующими материалами.

С помощью кнопки ESC вернемся в меню 7.2 и выберем пункт Subgrid. Распахивается таблица:

subgrid definition

subgrid name

в окошке последовательно появляются названия расчетных сеток (если их несколько). Клавишей 'пробел' найдем нужную и клавишей Enter выберем ее. При этом распахивается меню 3. Выбираем пункт FILL (заливка). После этого возникает меню 11:

 

 

Меню 11

Block

Rectangle

Ellipse

Parabola

Halfspace

Quad

Datafile

Sub. Fil

Unused

View

Block - заливка прямоугольника (используются номера ячеек)

Rectangle - заливка прямоугольника (используются физические координаты)

Ellipse - заливка эллипса (круга, части круга)

Parabola - заливка параболы

Halfspase - заливка полупространства

Quad - произвольный четырехугольник (идти против часовой стрелки)

Datafile - загрузка из внешнего файла

Sub. Fil -

Unused - исключение из расчета части расчетной области

View - просмотр полученной картинки

- выберем пункт меню 11, например Block. Возникает таблица:

Block region definition

 Заполним ее. В пункте "material name:" укажем название материала 'воздух' (GAIR) с помощью кнопки 'пробел', далее распахивается таблица "material conditions:". В пункте "initial int. energy:" надо обязательно указать значение начальной внутренней энергии воздуха (Для нормальной атмосферы оно равно E=2,068*10⁵mДж/mгр.!!!). Далее аналогично заливаем область ВВ. Выбрав пункт View из меню 11, просмотрим сформированную расчетную область уже заполненную материалами. Возникает меню 12:

Меню 12

Materials

Density

Energy

Velocity

-выберем пункт Materials, возникает меню 13:

Меню 13

Location

Status

- выберем пункт Location, на экране появляется картинка расчетной области с залитыми материалами, а также меню 6.2:

Меню 6.2

Ggid

Zoom

Pan

Examine

Reset

Axes

с помощью меню 6.2 можно детально просмотреть область расчета.

5. Точка инициирования

С помощью кнопки ESC вернемся в меню 7.2 и выберем пункт Global. При этом распахивается меню 8. Выбираем пункт Options. При этом распахивается меню 14:

Меню 14

Explode

Gravity

Damping

Cutoffs

Transport

Zone Side

UserVar.

- выбираем пункт Explode (взрыв). При этом распахивается меню 15:

Меню 15

Node

Line

Circle

Path

Times

Delete

Review

View

Выбираем пункт Node (узел, точка). При этом распахивается таблица с координатами точки инициирования, заполняем ее. После этого на экране возникает надпись "Limit range of application:", выбираем "No". Далее в этом же меню 15 выбираем Path (путь детонационной волны). Распахивается таблица:

Path used compute det. times

Path will be

имеются две возможности: Direct и Indirect. Кнопкой 'пробел', а затем Enter выбираем Indirect (огибание детонационной волной инертных препятствий, например, линз). Для визуализации точки инициирования на области расчета в меню 15 выбираем пункт View. При этом распахивается меню 6.3:

Меню 6.3

Next

Zoom

Pan

Reset

Examine

Axes

ShowAll

При надобности используем это меню.

6. Постановка граничных условий

6.1. Выбор типов граничных условий

Выберем нужные нам типы ГУ. Для этого в меню 7.2 выберем пункт Global, при этом откроется меню 8, выберем пункт Boundary, при этом откроется меню 9.1:

Меню 9.1

Add/Mod

Delete

Review

выберем пункт Add/Mod. При этом распахивается таблица:

Boundary condition:

Впишем в свободное поле рабочее название типа ГУ, например, "FLOW OUT", при этом распахнется новая таблица:

Boundary conditions: FLOW

Boundary type:

В рабочем поле таблицы с помощью клавиши 'пробел' устанавливаем нужный тип граничных условий. Имеются следующие типы ГУ:

Stress

X- velocity

Y- velocity

R- velocity

A- velocity

Angle

Flow in

Flow out

Transmit

- выберем, например, "Flow out". При этом распахивается таблица:

Boundary condition: FLOW- Flow

Allow reverse flow?:   (втекание потока назад в расчетную область)

Имеются две возможности: Yes и No. Выберем No. При этом распахнется таблица:

Flow out: Flow out (Euler)

Flow material:

В рабочем поле таблицы содержатся все объявленные при построении задачи материалы (GAIR, GEXPLOSIVE,VOID) и пункт ALL EQUEL (все материалы). Выберем ALL EQUEL.

6.2. Расстановка ГУ на расчетной области

В меню 7.2 выбираем пункт Subgrid. Возникает меню 3, выбираем пункт Boundary. При этом возникает таблица с названием расчетной области (возможно их будет несколько). Выбраем нужное название. При этом возникает меню 16:

Меню 16

Node

I-Line

J-Line

Block

Clear

View

Расставляем ГУ на границах расчетной области: сперва вдоль оси I, потом вдоль оси J. Расставим ГУ, например, вдоль оси I. Для этого открываем пункт I-Line в меню 16. При этом распахивается таблица: 

I-Line definition

I index:     

From J index: to J index:

Boundary condition:

Заполняем ее. Далее если нужно выбираем в меню 16 пункт J-Line и проводим аналогичные операции. Для просмотра на экране полученного результата выбираем пункт View в меню 16. При этом распахивается сервисное меню 6.4:

 

Меню 6.4

Next

Zoom

Pan

Reset

Axes

Замечание:

1. При отсутствии явно заданного граничного условия на каком то краю расчетной сетки, по умолчанию предполагается, что здесь используется условие непротекания.

2. По умолчанию при аксиальной симметрии горизонтальная нижняя ось ОХ (OI) является осью симметрии задачи.

7. Установка параметров расчета

Установим служебные параметры для выполнения расчета построенной задачи.

7.1 Тип симметрии.

В меню 7.2 выберем пункт Global, при этом откроется меню 8, выберем пункт Symmetry. При этом распахивается таблица:

Type of symmetry

Symmetry type:

При этом в рабочей области содержатся параметры Axial (аксиальная симметрия) и Planar (плоская). Устанавливаем с помощью кнопки 'пробел' параметр Axial. По умолчанию используется параметр Axial, поэтому при осевой симметрии задачи задание типа симметрии (пункт Symmetry) не обязательно.

7.2 Задание шага интегрирования по времени.

 В меню 7.2 выберем пункт Global, при этом откроется меню 8, выберем пункт Timestep. При этом распахивается таблица:

Timestep Control Parameters

Start time:

Minimum time step:

Maximum time step:

Initial time step:

Safety factor: (число Куранта, по умолчанию задано 0.6666)

В рабочих полях таблицы по умолчанию заданы все эти параметры, поэтому их можно не задавать.

7.3 Время окончания расчета

В меню 7.2 выберем пункт Global, при этом откроется меню 8, выберем пункт Wrapup (закончить). При этом распахивается таблица: