>> %--Интенсивность землетрясения
clear all;
close all;
m=8;
h=12;
h2=h^2;
l=33;
l1=0.7*l;
l2=7*l;
del=1.5;
i1=1.5*m-3.5*log10(h)+3
q=exp(2.3026*(1.5*m+4.8))
x=l1:1:l2;
y=i1-6*log10(sqrt(x.^2+h2)/h)+del;
plot(x,y,'r-','LineWidth',2);
grid on;
xlabel('L, км');
ylabel('J, баллы')
i1 =
11.2229
q =
6.3112e+16
>> %--Подъем уровня воды в реке
clear all;
close all;
j=15;
f=250;
n=0.032;
h0=4.7;
b0=160;
b1=140;
r=140;
k=0.7;
i=1.2e-4;
v0=h0^(2/3)*i^0.333/n
m=2*b0/(b1+b0);
s0=h0*(b0+b1)/2;
q0=s0*v0;
qm=q0+k*j*f/3.6;
a=3/(3*m+2);
hm=h0*(((qm/q0)^a)-1)
a1=(1+hm/h0)^(2/3);
vm=v0*a1
vg=1.5*vm
t=r/vg/3.6
a2=k*f/3.6/q0;
x=10:1:30;
y=((((x.*a2)+1).^a)-1)*h0;
plot(x,y,'m-','LineWidth',2);
grid on;
xlabel('Интенсивность осадков, мм/час');
ylabel('Подъем уровня воды, м');
v0 =
4.3379
hm =
0.6171
vm =
4.7098
vg =
7.0647
t =
5.5047
>> % --Затопление прибрежной территории
clear all;
close all;
h0=4.7;
f1=0.026;
f2=0.32;
b0=160;
b1=140;
hm=0.62;
h=2.2;
l=84;
tga=(b0-b1)/2/h0;
b2=(h0+h)*tga;
b3=(h0+hm)*tga;
s=(b1+b2+b3)*(hm-h);
tgf1=tan(f1);
tgf2=tan(f2);
a=0.5*tgf1+0.5*tgf1*tgf1/tgf2;
c=(b1+2*b2)*tgf1;
ca=c/2/a;
ca2=ca*ca;
p=[a c -s];
lz1=roots(p)
x=2.5:0.1:4.5;
y=sqrt(ca2+(b1+b2+(h0+x).*tga).*(x-h)/a)-ca;
plot(x,y,'m-','LineWidth',2);
grid on;
xlabel('Максимальный уровень подъема воды hm, м');
ylabel('Глубина затопления территории Lz, м')
lz1 =
-234.2253
-79.8509
>> %--Интенсивность землетрясения
clear all;
close all;
m=8;
h=18;
h2=h^2;
l=30;
l1=0.7*l;
l2=7*l;
del=1.5;
i1=1.5*m-3.5*log10(h)+3
q=exp(2.3026*(1.5*m+4.8))
x=l1:1:l2;
y=i1-6*log10(sqrt(x.^2+h2)/h)+del;
plot(x,y,'r-','LineWidth',2);
grid on;
xlabel('L, км');
ylabel('J, баллы')
i1 =
10.6065
q =
6.3112e+16
>> %--Интенсивность землетрясения
clear all;
close all;
m=8;
h=12;
h2=h^2;
l=33;
l1=0.7*l;
l2=7*l;
del=1.5;
i1=1.5*m-3.5*log10(h)+3
q=exp(2.3026*(1.5*m+4.8))
x=l1:1:l2;
y=i1-6*log10(sqrt(x.^2+h2)/h)+del;
plot(x,y,'r-','LineWidth',2);
grid on;
xlabel('L, км');
ylabel('J, баллы')
i1 =
11.2229
q =
6.3112e+16
>> %--Подъем уровня воды в реке
clear all;
close all;
j=15;
f=250;
n=0.032;
h0=4.7;
b0=160;
b1=140;
r=140;
k=0.7;
i=1.2e-4;
v0=h0^(2/3)*i^0.333/n
m=2*b0/(b1+b0);
s0=h0*(b0+b1)/2;
q0=s0*v0;
qm=q0+k*j*f/3.6;
a=3/(3*m+2);
hm=h0*(((qm/q0)^a)-1)
a1=(1+hm/h0)^(2/3);
vm=v0*a1
vg=1.5*vm
t=r/vg/3.6
a2=k*f/3.6/q0;
x=10:1:30;
y=((((x.*a2)+1).^a)-1)*h0;
plot(x,y,'m-','LineWidth',2);
grid on;
xlabel('Интенсивность осадков, мм/час');
ylabel('Подъем уровня воды, м');
v0 =
4.3379
hm =
0.6171
vm =
4.7098
vg =
7.0647
t =
5.5047
>> %--Подъем уровня воды в реке
clear all;
close all;
j=15;
f=250;
n=0.032;
h0=4.7;
b0=160;
b1=140;
r=140;
k=0.7;
i=1.2e-4;
v0=h0^(2/3)*i^0.333/n
m=2*b0/(b1+b0);
s0=h0*(b0+b1)/2;
q0=s0*v0;
qm=q0+k*j*f/3.6;
a=3/(3*m+2);
hm=h0*(((qm/q0)^a)-1)
a1=(1+hm/h0)^(2/3);
vm=v0*a1
vg=1.5*vm
t=r/vg/3.6
a2=k*f/3.6/q0;
x=10:1:30;
y=((((x.*a2)+1).^a)-1)*h0;
plot(x,y,'m-','LineWidth',2);
grid on;
xlabel('Интенсивность осадков, мм/час');
ylabel('Подъем уровня воды, м');
v0 =
4.3379
hm =
0.6171
vm =
4.7098
vg =
7.0647
t =
5.5047
>> % --Затопление прибрежной территории
clear all;
close all;
h0=4.7;
f1=0.026;
f2=0.32;
b0=160;
b1=140;
hm=0.62;
h=2.2;
l=84;
tga=(b0-b1)/2/h0;
b2=(h0+h)*tga;
b3=(h0+hm)*tga;
s=(b1+b2+b3)*(hm-h);
tgf1=tan(f1);
tgf2=tan(f2);
a=0.5*tgf1+0.5*tgf1*tgf1/tgf2;
c=(b1+2*b2)*tgf1;
ca=c/2/a;
ca2=ca*ca;
p=[a c -s];
lz1=roots(p)
x=2.5:0.1:4.5;
y=sqrt(ca2+(b1+b2+(h0+x).*tga).*(x-h)/a)-ca;
plot(x,y,'m-','LineWidth',2);
grid on;
xlabel('Максимальный уровень подъема воды hm, м');
ylabel('Глубина затопления территории Lz, м')
lz1 =
-234.2253
-79.8509
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=130; % --Левая точка интерполяции из таблицы
d1=40;
w2=120; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.08;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=90:1:160; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(150.3-r1)/5;
pr=log(t*wi1^1.333)*2.99-11.4; %--Пробит для ожога 2 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
50.8937
p =
0.8792
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.08;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=70:1:150; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(150.3-r1)/5;
pr=log(t*wi1^1.333)*2.99-11.4; %--Пробит для ожога 2 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
50.8937
p =
0.8792
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.08;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=40:1:150; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(150.3-r1)/5;
pr=log(t*wi1^1.333)*2.99-11.4; %--Пробит для ожога 2 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
50.8937
p =
0.8792
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.08;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=30:1:150; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(150.3-r1)/5;
pr=log(t*wi1^1.333)*2.99-11.4; %--Пробит для ожога 2 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
50.8937
p =
0.8792
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1) %--Плотность теплового излучения
vm=0.08;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=40:1:150; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(40-r1)/5;
pr=log(t*wi1^1.333)*2.99-11.4; %--Пробит для ожога 2 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
w =
11.7574
hp =
50.8937
Error using erf
Input must be real and full.
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.04;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=42:1:80; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(150.3-r1)/5;
pr=log(t*wi1^1.333)*2.99-11.4; %--Пробит для ожога 2 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
33.3454
p =
0.8792
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.04;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=42:1:80; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(42-r1)/5;
pr=log(t*wi1^1.333)*2.99-9.16; %--Пробит для ожога 1 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
33.3454
Error using erf
Input must be real and full.
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.04;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=42:1:80; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(42-r1)/5;
pr=log(t*wi1^1.333)*2.99-9.16; %--Пробит для ожога 1 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
33.3454
Error using erf
Input must be real and full.
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.04;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=42:1:80; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(42-r1)/5;
pr=log(t*wi1^1.333)*2.99-9.16; %--Пробит для ожога 1 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
33.3454
Error using erf
Input must be real and full.
>> %--Пожар пролива
clear all;
close all;
m=6e4; % --Масса ГЖ 36 т
r1=80; %--Расстояние до человека 100 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.04;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=42:1:80; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(60-r1)/5;
pr=log(t*wi1^1.333)*2.99-9.16; %--Пробит для ожога 1 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
33.3454
p =
3.2694e-08
>> %--Пожар пролива
clear all;
close all;
m=60000; % --Масса ГЖ 60 т
r1=80; %--Расстояние до человека 80 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.04;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=40:1:80; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(60-r1)/5;
pr=log(t*wi1^1.333)*2.99-9.16; %--Пробит для ожога 1 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
33.3454
p =
3.2694e-08
>> %--Пожар пролива
clear all;
close all;
m=60000; % --Масса ГЖ 60 т
r1=80; %--Расстояние до человека 80 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.04;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=40:1:80; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(41-r1)/5;
pr=log(t*wi1^1.333)*2.99-9.16; %--Пробит для ожога 1 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
33.3454
Error using erf
Input must be real and full.
>> %--Пожар пролива
clear all;
close all;
m=60000; % --Масса ГЖ 60 т
r1=80; %--Расстояние до человека 80 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.04;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=40:1:80; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(42-r1)/5;
pr=log(t*wi1^1.333)*2.99-9.16; %--Пробит для ожога 1 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
33.3454
Error using erf
Input must be real and full.
>> %--Пожар пролива
clear all;
close all;
m=60000; % --Масса ГЖ 60 т
r1=80; %--Расстояние до человека 80 м
ro=900;
ob=m/ro; %-- Объем жидкости
pl=ob/0.05;
d=sqrt(4*pl/3.14) %--Диаметр пятна, м
w1=12; % --Левая точка интерполяции из таблицы
d1=40;
w2=10; % --Правая точка интерполяции из таблицы
d2=50;
w=w1+(w2-w1)*(d-d1)/(d2-d1); %--Плотность теплового излучения
vm=0.04;
hp=42*d*(vm/1.29/sqrt(9.81*d))^0.61 % --Высота пламени
h=2*hp/d;
r=40:1:80; % --Расстояние до людей
s=2*r/d;
a=0.5*(h*h+s.*s+1)./s;
b=0.5*(1+s.*s)./s;
art=atan(sqrt((a+1).*(s-1)./(a-1)./(s+1)));
fv=(atan(h./sqrt(s.*s-1))-h.*(atan(sqrt((s-1)./(s+1)))...
-a./sqrt(a.*a-1).*art))./s/3.14;
fh=((b-1./s)./sqrt(b.*b-1).*atan(sqrt((b+1).*(s-1)./(b-1)./(s+1)))...
-(a-1./s)./sqrt(a.*a-1).*art)/3.14;
fi=sqrt(fv.*fv+fh.*fh);
k=exp(-(r-0.5*d)*7e-4);
wi=w*fi.*k; % Интенсивность излучения, кВт/м2
plot(r,wi,'m-','LineWidth',2);
grid on;
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2'); %--Определение интенсивности излучения на r1 - 9кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
%--По точному графику это расстояние 150,3 м
wi1=9;
t=5+(43-r1)/5;
pr=log(wi1^1.333)*2.99-9.16; %--Пробит для ожога 1 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии 100 м
d =
41.2130
hp =
33.3454
p =
3.2694e-08
>> %--Пожар "Огненный шар"
clear all;
close all;
m=85; %--Масса газа, кг
r1=60; %--Расстояние, м
d=5.33*m^0.327; %--Диаметр шара
h=d/2;
hd=h/d+0.5;
t=0.92*m^0.303; %--Время горения
r=50:1:80; %--Расстояние до человека
rd=r./d;
fi=hd/4./(hd*hd+rd.*rd).^1.5;
k=exp(-7e-4*(sqrt(r.*r+h*h)-d/2));
wi=450*fi.*k; %--Интенсивность облучения
plot(r,wi,'m-','LineWidth',2);
grid on;
title('Пожар "Огненный шар"');
xlabel('Расстояние до человека, м');
ylabel('Интенсивность теплового излучения, кВт/м2');
%--Определение интенсивности излучения на r1 - 4.9 кВт/м2
%--Определение расстояния, где wi=4 кВт/м2
wi1=4.9;
t=5+(65-r1)/5;
pr=log(t*wi1^1.333)*2.99-11.4; %--Пробит для ожога 2 степени
p=0.5*(1+erf((pr-5)/1.414)) %--Вероятность ожога на расстоянии r1
p =
1.2436e-06
Дата добавления: 2015-08-29; просмотров: 33 | Нарушение авторских прав
| <== предыдущая лекция | | | следующая лекция ==> |
| Московская ложа Святого Мафусаила 1 | | | Военная доктрина Российской Федерации |