10.2.1. Исходные данные:

1. Вид взрыва: а) воздушный, б) наземный;

2. Мощность взрыва;

3. Расстояние до ОЭ;

4. Раcстояние до района рассредоточения;

5. Коэффицент ослабления атмосферы;

6. Скорость ветра;

7. Угол между осью следа радиоактивного облака и линией, проведенной через ОЭ и эпицентр взрыа;

10.2.2. Выходные данные:

1. Избыточное давление во фронте ударной волны;

2. Импульс светового излучения;

3. Суммарная доза гамма-излучения;

4. Мощность дозы Г-излучения;

5. Поток нейтронов;

6. Вертикальная составляющая эл.поля ЭМИ;

7. Уровень радиации радиоактивного заражения;

10.3. Текст программы

program voina;

var

Pvzr,vid,rr,vv,bet:real;

R_onx,ko,dPf,qy,R,rs,U,Fn,Pg,Dz,Dosk,Dg,E,P0,alf:real;

begin

writeln('Введите вид взрыва,если взрыв воздушный -> нажми 1');

writeln(' если взрыв наземный -> нажми 2');

read(vid);

writeln('Введите мощность взрыва,Kт');

read(Pvzr);

writeln('Введите расстояние до ОЭ,км');

read(R_onx);

writeln('Введите раcстояние до района рассредоточения,км');

read(rr);

writeln('Введите коэфицент ослабления');

read(ko);

writeln('Введите скорость ветра,км/ч');

read(vv);

writeln('Введите угол,град');

read(bet);

qy:=0.5*Pvzr*1000000;

R:=R_onx*1000;

dPf:=105/R*exp(1/3*ln(qy))+410/R/R*exp(2/3*ln(qy))+1370/r/r/r*qy;

if vid=1 then

rs:=0.052*exp(0.4*ln(Pvzr))

else

rs:=0.068*exp(0.4*ln(Pvzr));

R:=R_onx;

U:=111*Pvzr/R/R*exp(-ko*(R-rs));

R:=R*1000;

Fn:=7.5*exp(22*ln(10))/R/R*Pvzr*exp(-R/190);

Pg:=exp(13*ln(10))/R/R*Pvzr*exp(-R/200);

Dz:=5*exp(8*ln(10))/R/R*Pvzr*exp(-R/410);

Dosk:=1.4*exp(9*ln(10))*Pvzr*(1+0.2*exp(0.65*ln(Pvzr)))/R/R*exp(-R/300);

Dg:=Dz+Dosk;

R:=R_onx;

alf:=Pi/4-2*bet*Pi/180;

E:=5*exp(3*ln(10))*(1+2*R)/R/R/R*ln(14.5*Pvzr)/ln(10);

P0:=10*Pvzr/(exp(1.5*ln(rr/22))*exp(sqrt(rr/vv)))*sqr(sqr(sin(alf)/cos(alf)));

writeln('Избыточное давление во фронте ударной волны: ',dPf:1:3,' кПа');

writeln('Импульс светового излучения: ',U:1:3,' кДж/м¤');

writeln('Суммарная доза гамма-излучения: ',Dg:1:3);

writeln('Мощность дозы Г-излучения: ',Pg:1:3);

writeln('Поток нейтронов: ',Fn:3,' н/м¤');

if vid=2 then

writeln('Вертикальная составляющая эл.поля ЭМИ: ',E:1:3,' В/м');

writeln('Уровень радиации радиоактивного заражения: ',P0:1:3);

writeln('');

end.

10.4. Проврка работоспособности

Избыточное давление во фронте

ударной волны P_ф,кПа

R/q	100	200	300	400	500
3	21.614	31.134	39.068	46.196	52.813
4	14.219	19.846	24.393	28.396	32.057
5	10.51	14.377	17.432	20.082	22.478
6	8.31	11.211	13.465	15.398	17.13
7	6.861	9.164	10.931	12.432	13.769
8	5.839	7.740	9.184	10.402	11.479
9	5.079	6.694	7.91	8.93	9.828
10	4.493	5.894	6.942	7.817	8.584

где R - расстояние до центра взрыва, км;

q - мощность взрыва, кт.

10.5. Выводы:

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