【爆破载荷参数】基于 UFC 3-340-02 / TM 5-855-02 的爆炸压力效应研究（Matlab代码实现）

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码及数据

💥1 概述

本文涵盖了以下内容：
- 爆炸载荷参数：包括爆炸药剂的特性、产生的爆炸压力、冲击波、破片及热辐射等信息。
- 建筑物和结构物的响应：介绍了建筑物和结构物受到爆炸作用时的响应及破坏机制，如振动、应力、位移等。
- 安全设计准则：提供了设计工程和基础设施时需要遵循的安全准则，以减轻爆炸压力效应对结构物、人员和设备的影响。
- 防爆保护措施：包括爆炸抗冲击设计、隔离区域规划、材料选用和防爆装置的使用等，以提高建筑物和设备的抗爆性能。

📚2 运行结果

FAB = matfile('FreeAirBurstParametersPositivePhaseMetric.mat');
figure; hold on; box on; grid on;
title({'Positive phase shock wave parameters for a spherical';...
'TNT explosion in free air at sea level'})
xlabel('Scaled Distance Z = R/W^{1/3} [m/kg^{1/3}]')
axis([0.05 50 0.005 1e6])
plot(FAB.ScaledDistance1 ,FAB.ReflectedPressure)
plot(FAB.ScaledDistance2 ,FAB.PeakIncidentOverpressure)
plot(FAB.ScaledDistance3 ,FAB.ScaledReflectedImpuls)
plot(FAB.ScaledDistance4 ,FAB.ScaledIncidentImpuls)
plot(FAB.ScaledDistance5 ,FAB.ScaledArivalTime)
plot(FAB.ScaledDistance6 ,FAB.ScaledPositivePhaseDuration)
plot(FAB.ScaledDistance7 ,FAB.WaveFrontSpeed)
plot(FAB.ScaledDistance8 ,FAB.ScaledWaveLength)
set(gca,'yscale','log');
set(gca,'xscale','log');
legend('P_r, kPa','P_{so}, kPa','i_r/W^{1/3}, kPa-ms/kg^{1/3}','i_s/W^{1/3}, kPa-ms/kg^{1/3}',...
       't_a/W^{1/3}, ms/kg^{1/3}','t_0/W^{1/3}, ms/kg^{1/3}','U, m/ms','L_w/W^{1/3}, m/kg^{1/3}')

HSB = matfile('HemisphericalBurstParametersPositivePhaseMetric.mat');
figure; hold on; box on; grid on;
title({'Positive phase shock wave parameters for a hemispherical';...
'TNT surface explosion at sea level'})
xlabel('Scaled Distance Z = R/W^{1/3} [m/kg^{1/3}]')
axis([0.05 50 0.005 1e6])
plot(HSB.ScaledDistance1 ,HSB.ReflectedPressure)
plot(HSB.ScaledDistance2 ,HSB.PeakIncidentOverpressure)
plot(HSB.ScaledDistance3 ,HSB.ScaledReflectedImpuls)
plot(HSB.ScaledDistance4 ,HSB.ScaledIncidentImpuls)
plot(HSB.ScaledDistance5 ,HSB.ScaledArivalTime)
plot(HSB.ScaledDistance6 ,HSB.ScaledPositivePhaseDuration)
plot(HSB.ScaledDistance7 ,HSB.WaveFrontSpeed)
plot(HSB.ScaledDistance8 ,HSB.ScaledWaveLength)
set(gca,'yscale','log');
set(gca,'xscale','log');
legend('P_r, kPa','P_{so}, kPa','i_r/W^{1/3}, kPa-ms/kg^{1/3}','i_s/W^{1/3}, kPa-ms/kg^{1/3}',...
       't_a/W^{1/3}, ms/kg^{1/3}','t_0/W^{1/3}, ms/kg^{1/3}','U, m/ms','L_w/W^{1/3}, m/kg^{1/3}')

RPC = matfile('ReflectedPressureCoefficientVsAngleOfIncidenceMetric.mat');
figure; hold on; box on; grid on;
title('Reflected Pressure Coefficient')
view(3)
surf(RPC.AngleIncrements, RPC.PressureIncrements,RPC.AngleOfInclanationCoeficientMatrix)
xlabel('angle [deg]')
ylabel('Peak incident overpressure [kPa]')
zlabel('C_{r \alpha} = {P_r \alpha} / P_{SO}')
set(gca,'yscale','log');

% or as a line plot
figure; hold on; box on; grid on;
plot(RPC.AngleIncrements,RPC.AngleOfInclanationCoeficientMatrix)
title('Reflected Pressure Coefficient')
xlabel('angle [deg]')
ylabel('C_{r \alpha} = {P_r \alpha} / P_{SO}')

RSI = matfile('ReflectedScaledImpulseVsAngleOfIncidenceMetric.mat');
figure; hold on; box on; grid on;
title('Reflected Scaled Impulse')
view(3)
surf(RSI.AngleIncrements, RSI.PressureIncrements,RSI.ReflectedScaledImpulseMatrix)
xlabel('angle [deg]')
ylabel('Peak incident overpressure [kPa]')
zlabel('i_{r \alpha} / W^{1/3} [kPa-msec/kg^{1/3}]')
set(gca,'yscale','log');
set(gca,'zscale','log');

% or as a line plot
figure; hold on; box on; grid on;
plot(RSI.AngleIncrements, RSI.ReflectedScaledImpulseMatrix)
title('Reflected Scaled Impulse')
xlabel('angle [deg]')
ylabel('i_{r \alpha} / W^{1/3} [kPa-msec/kg^{1/3}]')
set(gca,'yscale','log');

MSH = matfile('MachStemHightMetric.mat'); 
figure; hold on; box on; grid on; grid minor
title('Scaled height of triple point')
xlabel('Scaled Horizontal Distance from Charge, m/kg^{1/3}')
ylabel('Scaled Height of Triple Point, H_{T}/W{^1/3} (m/kg^{1/3})')
axis([0 9 0 5])
plot(MSH.ScaledHorizontalDistanceMatrix', MSH.ScaledTriplePointHightMatrix','k')
text(MSH.ScaledHorizontalDistanceMatrix(:,end),MSH.ScaledTriplePointHightMatrix(:,end),...
     num2str(round(MSH.ScaledChargeHight',3)))
str1 = {'Number adjacent to curves','indicate scaled charge','height, H_c/W^{1/3}'};
text(8.8, 4.8, str1, 'HorizontalAlignment', 'right', 'VerticalAlignment', 'top',...
    'EdgeColor', 'k', 'Margin', 1.5, 'LineWidth', 1,'BackgroundColor','w')

figure; hold on; box on; grid on;
title('Angular reflected pressure')
xlabel('Horisonal distance form charge [m]')
ylabel('Pressure [kPa]')
plot(X,Pra_line)

 

figure; hold on; box on; axis equal;
colormap jet
ylabel(colorbar,'Pressure [kPa]')
xlabel('Distance from center [m]')
ylabel('Distance from center [m]')
title('Angular reflected pressure')
contourf(X,flipud(Y'),Pra_contour)

 

W = 0.9; % [kg] Charge
R = 2.5; % [m] Distance
Type = 'Surface'; % Detonation type
[Pr, Pso, ir, is, ta, to, U, Lw, Z] = BlastParameters(R, W, Type);
b = DecayCoefficient(Pso,is,to); % Decay coefficient for the Friedlander's equation
t = linspace(0,to,25); % [ms]
Ps = FriedlandersEquation(Pso, to, b, t);

figure; hold on; box on; grid on;
title('Incident overpressure history')
xlabel('Time, ta = 0, [ms]')
ylabel('Pressure [kPa]')
plot(t,Ps)

部分代码：

function Ps = FriedlandersEquationReal(Pso, ta, to, b, t)
%FriedlandersEquationReal Determines the pressure at time t following the
%Friedlander's equation for blast waveform.
%Ps = FriedlandersEquationReal(Pso, ta, to, b, t)
% INPUT
% Pso [kPa] Peak overpressure
% to  [ms]  Positive phase duration
% b   [-]   Decay coefficient for the waveform
% t   [ms]  Time 0 < t < ta+to (zero is time of detonation)
% 
% OUTPUT
% Ps [kPa] Overpressure at time t
% 
% The Friedlander's equation is given in the form of:
%           Ps(t) = Pso*(1-(t-ta)/to)*exp(-b*(t-ta)/to)
% 
% The units are only recommendations and may be changed at own risk.
% The start time, ta, is NOT set to zero.

% Determine the pressure at time t
Ps = Pso*(1-(t-ta)./to).*exp(-b*(t-ta)./to);

% Correct time for special cases
Ps(t<ta) = 0;
Ps(t>(ta+to)) = 0;

end

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1] Lawver D , Weeks J , Forman S ,et al.Cost Savings Using Stirrup Reinforcement Instead of Laced Wall Reinforcement Per UFC 3-340-02 (TM 5-1300) December, 2008[J].  2010.

[2] Jr W H Z , Acosta P F .SUMMARY OF THE NEW REINFORCED CONCRETE BLAST DESIGN PROVISIONS IN UFC 3-340-02, 'STRUCTURES TO RESIST THE EFFECTS OF ACCIDENTAL EXPLOSIONS'[C]//ACI Fall convention.DoD Explosives Safety Board in Alexandria, Virginia; US Army Engineering and Support Center, Huntsville, 2010.

[3]彭琦,吴昊,方秦,等.长持时平面爆炸波作用下RC梁动力响应研究[J].建筑结构学报, 2023, 44(3):15.

🌈4 Matlab代码及数据