Рябенький В.С. Різницеві схеми. М .: Наука, 1973.
.Огнева В.В. Метод прогонки для вирішення різницевих рівнянь.// ЖВМіМФ.- 1967. -№4.
.Потемкін В.Г. Система інженерних і наукових розрахунків MATLAB 5.х.- М .: Діалог-МІФІ, 1 999.
. Іванов А. В., Мастрюков Б. С. Прогнозування поширення аварійних викидів аміаку в атмосфері міста.// Изв. Вищих навчальних закладів. Кольорова металургія - 1998. - №6.
. Методичні рекомендації щодо оцінки соціально-економічного збитку від порушення здоров'я населення, зумовленої забрудненням атмосферного повітря./А. А. Биков, Л. Соленова, Г. Земляна та ін.//Управління ризиком.- 1999. - № 3.
. Інструкція щодо складання планів ліквідації (локалізації) аварій в металургійних і коксохімічних виробництвах./Держнаглядохоронпраці.- М .: тисяча дев'ятсот дев'яносто чотири.
. Основи інженерної психології./Под ред. Ломова Б.Ф.- М .: Вища школа, 1986.
. Довідник практичного лікаря./Ю. Є. Вельтищев, Ф. І. Комаров, С. М. Навашин та ін. - М .: Баян, 1992.
53. Rosser R., Kund P.//Int. J. Of Epidemiology.- 1978, № 7, p. 347-358.
54. Биков А. А., Мурзін Н.В. Проблеми аналізу безпеки людини, суспільства, природи.- СПб .: Наука, +1997.
. Основи інженерної психології./Под ред. Ломова Б.Ф.- М .: Вища школа, 1986.
. Мастрюков Б.С. Комп'ютер - друг чи ворог?// Сталь.- №15-16.- Тисяча дев'ятсот дев'яносто вісім.
.Сідоров А.П. Секрети здоров'я. М .: Фізкультура і спорт, 1990.
.Воробьев В.І. Оздоровлення. М .: Центр, 1992.
ДОДАТОК 1
ПРОГРАМА РОЗРАХУНКУ ПОЛІВ ШВИДКОСТЕЙ
Errortrap ononshort e on
% Тротиловий еквівалент вибуху посудини з токсикантом
ctnt=26.9
% Розмір елементарної комірки
x0=1
y0=1
% Розміри розрахункової області
dl=100
shir=100
% Інтервал розбиття розрахункової області
dx=x0
% Розрахунок надлишкового тиску у фронті ударної хвилі
for x=1:dx:dly=1:dy:shir(x,y)=101300+(113.4*(ctnt)^0.33)/sqrt(x^2+y^2)+185.9*(ctnt^0.67)/(x^2+y^2)+9.02*ctnt/(x^2+y^2)^1.5)*1000(x,y)=(P0(x,y)- 101300)/101300 (x, y)=1.29 * (1 + DP0 (x, y) * 0.85)/(1 + DP0 (x, y) * 0.147) (x, y)=(P0 (x, y ) - 101300) * DP0 (x, y)/(0.4 * DP0 (x, y) +2.8) (x, y)=sqrt (2 * PSK (x, y)/R0 (x, y))
End
% Побудова графіка поля швидкостей у фронті хвилі
[X, Y]=meshgrid ([1: dx: dl]) (1) tall (1,1,1) gray fill
contour (V0)
title (Поле швидкостей в прямій хвилі, zet=0)
xlabel (Координата X, дам)
ylabel (Координата Y, дам)
% Тривалість фази стиснення
for x=1: dx: dl
for y=1: dy: shir (x, y)=ctnt ^ 0.167 * (x ^ 2 + y ^ 2) ^ 0.25 * 0.001
End
% Зміна швидкостей за фронтом ударної хвилі
for zet=0.1: 0.1: 0.9x=1: dx: dly=1: dy: shirDP0 (x, y) lt;=3 (x, y)=(- 1.33) * DP0 (x, y) DP0 (x, y) gt; 10 (x, y)=0 (x, y)=0.33 * DP0 (x, y) - 5.6DP0 (x, y) lt;=1 (x, y)=3.2 * DP0 (x, y) - 0.75 (x, y)=de (x, y) + 6.4 * DP0 (x, y)/(1 + 0.725 * DP0 (x, y) * zet) ( x, y)=PSK (x, y) * (1-zet) ^ 2 * exp ((- 1) * be (x, y) * zet) (x, y)=sqrt (2 * PSKNOV (x, y)/R0 (x, y))
End
% Побудова графіка поля швидкостей за фронтом хвилі
[X, Y]=meshgrid ([1: dx: dl]) (zet * 10 + 1) tall (1,1,1) grayfill (VNOV) (Поле швидкостей за фронтом прямої хвилі , zet =)
xlabel (Координата X, дам)
ylabel (Координата Y, дам)
end
% Відстань дії відбитої від перепони хвилі
xwol=0.7 * dl
% Розрахунок поля швидкостей в зоні відбитої хвилі
for x=xwol: dx: dly=1: dy: shir (x, y)=(P0 (x, y) * 101 300 * 3.2-0.4 * P0 (x, y))/ (0.4 * P0 (x, y) + 1.4 * 101300) (x, y)=(P0 (x, y) - 101300)/101300 (x, y)=1.29 * (1 + DP0 (x, y) * 0.85)/(1 + DP0 (x, y) * 0.147) (x, y)=(P0 (x, y) - 101300) * DP0 (x, y)/(0.4 * DP0 (x, y) +2.8 ) (x, y)=sqrt (2 * PSK (x, y)/R0 (x, y))
End
% Побудова графіка поля швидкостей в відбитої хвилі
[X, Y]=meshgrid ([1: dx: dl]) (1) tall (1,1,1) gray fill
contour (V0)
title (Поле швидкостей в відбитої хвилі, zet=0)
