Расчет воздействия взрывной волны на препятствие

[Airat_M 0]

Доброго времени суток!

Есть необходимость просчитать воздействие взрывной волны на препятствие.

никогда не занимался подобными расчетами, а в справке не нашел подобных примеров.

подскажите пожалуйста где можно почитать как это сделать в ФВ.

буду признателен за любую помощь в данном вопросе!

[a_schelyaev 366]

Тут ничего сакрального - все как у всех.

Детонацию в чистом виде FV не ситает, поэтому придется придумывать такие начальные условия, чтобы хотя бы повторить давление во фронет ударной волны.

ВО всех методиках пляшут обычно от радиуса сферы с поверхности которой стартует ударная волна. В зависимости от тротилового эквивалента или массы облака и данного радиуса, высчитывают начальную температуру и давление.

Мы считали просто, брали фильтр (на шаблоне бокса или конуса), через наго задавали в ячейке температуру и давление, ставили Курант = 1 с начальным ограничением 1e-7 и вперед с песней.

И не забудь неотражающие граничные условия расставить на внешних стенках.

[Airat_M 0]

спасибо за оперативный ответ! будем пробовать ))

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

пока поставлена задача выбора формы защитного препятствия.

есть ли смысл создавать модель в масштабе 1к 1? или можно создать масштабную модель?

[a_schelyaev 366]

1:1

Масштабную то зачем? Все равно решающим фактором является не габарит, а количество расчетных ячеек.

:)

[Airat_M 0]

Начал пробовать ))

по ходу выполнения возникают вопросы

1. в физ параметрах-начальные значения-давление нужно выставлять атмосферное? т.е. 100кПа или оставить 0? я пока поставил на расчет с Р=0

2. "И не забудь неотражающие граничные условия расставить на внешних стенках."

это есть свободный выход или симметрия?

3. "Мы считали просто, брали фильтр (на шаблоне бокса или конуса), через наго задавали в ячейке температуру и давление"

создаю фильтр с установкой переменной в параллелепипед-а там только одна переменная задается...вопрос-нужно создать 2 параллелепипеда в одной и той же области, один с температурой, второй с давлением, так? я сделал так и задал давление 100 кПа и температуру пока наобум 700 град.

Изменено 15 июля 2010 пользователем Airat_M

[a_schelyaev 366]

1. Ставь какое хошь, но имей в виду, что в ГУ и НУ и Постпроцессоре переменные задаются и выводятся относительно них как приращения.

2. Нет, это неотражающие. Ну если у тебя версия нормальная (неломанная или не старая).

3. Нужно задавать два и незабывать режим Одноразовый включать.

[Airat_M 0]

a_schelyaev

Большое вам спасибо за помощь!

буду дальше разбираться )

[Airat_M 0]

Доброго времени суток!

У меня вопрос по построению графиков.

имеется следующая модель

параметры:

модель - полностью сжимаемая жидкость (галки скорость и турбулентность)

среда - воздух

куранта ставил=1 но потом установил фикс. шаг=1e-006 (1e-008 нет возможности столько ждать)

установил в шаблон конуса переменную давление Р=5 МПа с одноразовым выбросом.

В представленном расчете я создал на стенках здания и преграды супергруппы на которых установил слои «заливка давление» и «характеристики из давления». Характеристики записываются в файл *.GLO из которых планировалось построить графики «время - давление»

Однако открыв эти файлы(GLO) через эксель (данные -> получение внешних данных -> из текста) я получил значительное количество данных.

Обработку в экселе задавал следующую

1. формат данных - с разделителями

2. символом-разделителем является - знак табуляции

3. Указал формат данных «текстовый» для всех данных

4. Получил данные со следующими столбцами.

Время

Итерация

Поток массы

ИнтегралX

ИнтегралY

ИнтегралZ

СилаСТрениемX

СилаСТрениемY

СилаСТрениемZ

ОбщаяПлощадь

ТочкаПриложенияМоментаX

ТочкаПриложенияМоментаY

ТочкаПриложенияМоментаZ

МоментX

МоментY

МоментZ

МоментБезТренияX

МоментБезТренияY

МоментБезТренияZ

Площадь

<f> по площади

Интеграл

Поток частиц черезвыход, kg/s

Теперь главный вопрос: как из этих данных получить давление на каждой итерации для построения графиков?

[a_schelyaev 366]

Давление ГДЕ?

[Airat_M 0]

мне нужно получить график изменения давления на стенке здания по времени

Изменено 21 мая 2011 пользователем Airat_M

[Airat_M 0]

Уважаемый a_schelyaev! Не могли бы вы пояснить что я сделал не так?

Все что я делал я описал выше, но не пойму что за величины в гло файле.

Либо подскажите пожалуйста, в каком столбце данные по давлению на стенке?

[Fedor 1 597]

"неотражающие граничные условия расставить на внешних стенках" - в теории заколачивания гвоздя логика такая. Идет ударный импульс до конца, а потом отражается и в этом месте происходит удвоение давления на препятствие и оно и воздействует на доску...

Если логически перевернуть и считать, что препятствие внедряется в податливую воздушную среду, то получается что-то подобное...

[a_schelyaev 366]

мне нужно получить график изменения давления на стенке здания по времени

Если речь идет о стенке, то там получается куча точек, где есть давление.

Поэтому речь может идти об среднем давлении.

Если нужен max и min, то нужно построить параллепипед, в который будет попадать нужное ГУ и выводить по объему характеристику.

Все эти вещи обычно выводятся в слое Характеристиа по величине Давление по Объекту Супергруппа на интересующем вас ГУ.

Если Характеристика кладется в GLO файл, то на каждой итерации файл будет дополняться и вы получите всю разблюдовку в течении процесса.

Уважаемый a_schelyaev! Не могли бы вы пояснить что я сделал не так?

Все что я делал я описал выше, но не пойму что за величины в гло файле.

Либо подскажите пожалуйста, в каком столбце данные по давлению на стенке?

В GLO файл выводятся значаения характеристик из Характеристики.

Про Характеристику все подробно разжевано в документации в слоях Постпроцессора.

[alann 0]

Ну как всегда, сначала делаем, потом мануалы читаем :)

Вообще-то, FV слабовато ударные волны считает. Размазывание сильное, сходимость (и по сетке и по времени) никакая. Но если нужна оценка "по-быстрому", то сойдёт.

[a_schelyaev 366]

Очень очень очень удивлен такому мнению.

Этот вопрос исселдовали не раз - сходимость по сетке и шагу интегрирования нормальная.

В стационарной структуре точно. В бегущей волне подробно не смотрели, но там точно зависимость от сетки и шага присутствует.

[alann 0]

... исселдовали не раз - сходимость по сетке и шагу интегрирования нормальная.

В стационарной структуре точно. В бегущей волне подробно не смотрели...

В стационарной - да. И профиль и параметры скачка неплохие. Но речь шла про нестационарный случай.

Как пример, продемонстрируйте пожалуйста решение задачи Сода.

[a_schelyaev 366]

А можно для дураков наводку?

[alann 0]

На водку не дам а ссылочку пожалуйста: <noindex>http://en.wikipedia.org/wiki/Sod_shock_tube</noindex>

Очень интересно, на сколько ячеек у вас размажутся ступеньки в областях 3 и 4? И какая будет сходимость по сетке?

P.S. задача одномерная, т.е. по осям OY и OZ вроде как достаточно одной ячейки. Однако мой опыт эксплуатации FV2 говорит, что у вас с шаблоном не все так просто и по оси OY нужно 5 (ПЯТЬ!) ячеек. Хотя, не исключено, что это недостаток старой версии.