Explosion-resistant analysis method for building single structure, computer equipment and storage medium

Abstract

The present invention belongs to the technical field of explosion-resistant analysis of building structures, and specifically relates to an explosion-resistant analysis method, computer equipment, and storage medium for a single building structure. The method calculates the equivalent action time of shock wave overpressure and the equivalent explosion load peak value based on building unit information and explosion-resistant information. The single building structure is simplified as a concentrated mass point at each floor, a simplified structural calculation model and corresponding motion equations are established, the dynamic response of the vertex is calculated, and an explosion spectrum is obtained. The equivalent earthquake spectrum is obtained by adjusting earthquake parameters, and an elastic-plastic static analysis method is used to obtain a capacity spectrum. The explosion spectrum, equivalent earthquake spectrum, and capacity spectrum are plotted in the same coordinate system for comparative analysis, resulting in the explosion-resistant performance of the single building structure and the earthquake action equivalent to the explosion. The method simplifies the process of analyzing the explosion-resistant performance of a single building structure and the method for evaluating the explosion-resistant performance, improves the efficiency of the explosion-resistant performance analysis, and enables the parallel analysis of the explosion-resistant performance and the earthquake-resistant performance.

Description

Technical Field

[0001] The present invention belongs to the technical field of explosion-resistant analysis of building structures, and in particular relates to an explosion-resistant analysis method for a single building structure, a computer device and a storage medium. Background Art

[0002] Petrochemical companies face numerous explosion hazards during production and product storage. In actual projects, explosion-resistant analysis reports within chemical parks clearly stipulate that buildings must undergo explosion-resistant analysis. Buildings in office and living areas are subject to shock waves that have less energy than those at the center of the explosion, and therefore often require structural response analysis under blast loads. However, in actual design, explosion-resistant analysis is often too complex.

[0003] my country's petrochemical specifications propose simplified methods such as closed-form solutions, graphical methods, and numerical integration methods to solve single-degree-of-freedom systems. However, these methods can only verify single components such as walls, beams, and columns, and lack methods for analyzing the explosion-resistant performance of the entire building structure. Due to the lack of research on the evaluation of the overall explosion-resistant performance of building structures, many single-story or multi-story buildings such as comprehensive buildings, dormitories, and guard rooms in chemical parks lack explosion-resistant performance analysis, resulting in serious loss of life and property after the explosion. In response to this, Chinese patent application number CN202311295875.4 discloses an explosion response analysis method based on the mode decomposition response spectrum. This method only uses the mode decomposition response spectrum method to perform dynamic analysis on the structure, and cannot achieve parallel analysis of explosion performance and structural seismic performance, and therefore cannot obtain the relationship between the explosion resistance and seismic performance of the structure.

[0004] Therefore, how to simplify the explosion-resistant analysis method of building single structures, realize more rapid and effective explosion-resistant performance analysis of building single structures, and obtain the relationship between the explosion-resistant performance and seismic performance of the structure is a technical problem that needs to be solved urgently by technical personnel in this field. Summary of the Invention

[0005] The present invention provides a method for analyzing the explosion resistance of a single building structure. The method simplifies the process of analyzing the explosion resistance performance of a single building structure and the method for evaluating the explosion resistance performance of a single building structure, thereby improving the efficiency of analyzing the explosion resistance performance of a single building structure and realizing the parallel analysis of the explosion resistance performance and the seismic resistance performance of a single building structure, thereby obtaining the relationship between the explosion resistance and seismic resistance performance of a single building structure.

[0006] The method comprises the following steps:

[0007] S1. Obtain building unit information and explosion-proof information;

[0008] S2. Calculating the equivalent overpressure action time of the shock wave and the equivalent explosion load peak value based on the building unit information and explosion resistance information;

[0009] S3. Based on the building unit information, simplify the building unit structure in the form of concentrated mass points at each floor, and establish a simplified structural calculation model;

[0010] S4. Determine the stress condition of the simplified structural calculation model according to the equivalent action time of the shock wave overpressure and the equivalent explosion load peak value, and establish the corresponding motion equation;

[0011] S5. Calculate the dynamic response of the vertices of the simplified structural calculation model based on the motion equation to obtain an explosion spectrum; adjust the earthquake parameters so that the vertices of the explosion spectrum are located on the earthquake general curve to obtain an equivalent earthquake spectrum; and obtain a capacity spectrum using an elastic-plastic static analysis method based on the building unit information;

[0012] S6. Plot the explosion spectrum, equivalent earthquake spectrum and capacity spectrum into the same coordinate system for comparative analysis, and obtain the explosion resistance performance of the building structure and the earthquake effect equivalent to the explosion.

[0013] In a specific embodiment, in step S1, by PKPM The structural design software obtains the building unit information, which includes: representative value of particle gravity load at the floor, floor height, plane length of the building, width of the building and floor stiffness; the explosion-proof information is obtained through the project explosion safety assessment report, which includes: building name, location of the explosion, peak incident overpressure of the explosion shock wave and time of positive pressure action of the explosion shock wave.

[0014] In a specific embodiment, in step S2, the explosion is preset to occur directly in front of the building monomer structure, and the explosion loads borne by both sides of the building monomer structure are equal in magnitude and opposite in direction. Only the explosion forces of the front wall and the rear wall of the building monomer structure are considered, and the calculated shock wave overpressure equivalent action time is the shock wave overpressure equivalent action time of the front wall of the building monomer structure; the explosion impulse superposition of the front wall and the rear wall of the building monomer structure is also calculated, and based on the explosion impulse superposition of the front wall and the rear wall of the building monomer structure and the shock wave overpressure equivalent action time of the front wall of the building monomer structure, the equivalent explosion load peak is calculated in a constant impulse manner.

[0015] In a specific embodiment, in step S3, the simplified structural calculation model includes a single-degree-of-freedom structural calculation model and a multi-degree-of-freedom structural calculation model, and in step S4, the motion equations correspondingly include a single-degree-of-freedom system motion equation and a multi-degree-of-freedom system motion equation.

[0016] In a specific embodiment, in step S5, the motion equation of the single-degree-of-freedom system is calculated using a general solution plus a special solution to obtain the dynamic response of the vertex of the single-degree-of-freedom structural calculation model.

[0017] In a specific embodiment, in step S5, the motion equation of the multi-degree-of-freedom system is Pushover Analysis technology is used, and the modal analysis method is adopted to first simplify the multi-degree-of-freedom structural calculation model into a single-degree-of-freedom structural calculation model, and then the dynamic response of the vertex of the multi-degree-of-freedom structural calculation model is calculated.

[0018] In a specific embodiment, in step S5, after calculating the dynamic response of the vertex of the simplified structural calculation model, the dynamic response of the vertex is converted into the relationship between the spectral displacement and the spectral acceleration, and plotted as ADRS Coordinate curve diagram to obtain the structural explosion spectrum.

[0019] In a specific embodiment, in step S5, first, according to the earthquake influence coefficient response spectrum, ADRS The earthquake spectrum is obtained by transformation, and the earthquake parameters are adjusted so that the vertex of the explosion spectrum is located on the earthquake general curve to obtain the equivalent earthquake spectrum.

[0020] The present invention further provides a computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the above method when executing the computer program.

[0021] The present invention further provides a computer-readable storage medium storing a computer program, wherein the computer program implements the steps of the above method when executed by a processor.

[0022] The present invention simplifies the building monomer structure in the form of concentrated mass points at floors, establishes a simplified structural calculation model, determines the stress condition of the simplified structural calculation model according to the equivalent explosion load peak and the equivalent action time of the shock wave overpressure, and establishes a corresponding motion equation. According to the motion equation, the dynamic response of the vertex of the simplified structural calculation model is calculated, and the explosion spectrum is obtained according to the dynamic response of the vertex. By adjusting the seismic parameters, the vertex of the explosion spectrum is located on the seismic general curve, and the equivalent seismic spectrum is obtained. The explosion spectrum is combined with the equivalent seismic spectrum and the structural capacity spectrum in the same coordinate system for comparative analysis, and the explosion resistance of the building monomer structure and the earthquake action equivalent to the explosion are obtained. The process of analyzing the explosion resistance of the building monomer structure and the evaluation method of the explosion resistance of the building monomer structure are simplified, the efficiency of the explosion resistance analysis of the building monomer structure is improved, and the explosion resistance analysis of the building monomer structure and the seismic resistance analysis of the building monomer structure are carried out in parallel, so as to obtain the relationship between the explosion resistance and seismic resistance of the building monomer structure. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] Figure 1 Flowchart of an embodiment of the present invention. DETAILED DESCRIPTION

[0024] In order to enable those skilled in the art to better understand the technical solutions in this application, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of this application.

[0025] See also Figure 1 The present invention provides a method for analyzing explosion resistance of a building monomer structure, which specifically includes the following steps:

[0026] First, through PKPM Structural design software captures building unit information, including representative values ​​of point-of-mass gravity loads at each floor, floor height, building length, building width, and floor stiffness. Project explosion safety assessment reports also capture explosion resistance information, including building name, explosion location, peak incident overpressure of the blast wave, and duration of the blast wave's positive pressure.

[0027] Then, assuming the explosion occurs directly in front of the building's monolithic structure, since the explosion loads on both sides of the building's monolithic structure are equal in magnitude and opposite in direction, only the explosion forces on the front and rear walls of the building's monolithic structure are considered. Based on the explosion location, the peak incident overpressure of the explosion shock wave, the duration of the positive pressure of the explosion shock wave, the floor height, and the planar length and width of the building, the equivalent duration of the shock wave overpressure on the front wall of the building's monolithic structure is calculated, as well as the superposition of the impulses on the front and back of each floor under the action of the explosion. Based on the constant impulse, the peak equivalent explosion load of the particle at the floor is calculated. The calculation formula for the equivalent explosion load is:

[0028] P u =BH ( P r t e +P b ( t rb +t d )) / t e ;

[0029] Where, Pu represents the equivalent explosion load acting on the mass point at the floor, B Indicates the width of the building corresponding to the explosion source, H Indicates the floor height, P r represents the peak incident overpressure of the front wall shock wave, t e represents the equivalent action time of the front wall shock wave overpressure, P b represents the peak incident overpressure of the back wall shock wave, t rb represents the equivalent action time of the back wall shock wave overpressure, t d represents the time required for the shock wave to propagate from the front wall to the back wall, and / represents a division sign.

[0030] Based on the calculated equivalent explosion load peak value and the calculated equivalent overpressure action time of the front wall shock wave of the building structure, the equivalent explosion load time history function expression is established: P ( t ) =P u (1- t / t e );

[0031] Where, P ( t ) represents the equivalent explosion load that varies with time, P u represents the equivalent explosion load acting on the mass point at the floor, t represents the time variable, t e represents the equivalent duration of the front wall shock wave overpressure, and / represents the division sign. This time history function can more realistically simulate the time history of the equivalent blast load, capture the dynamic response of a building's individual structures, and thus more accurately predict the dynamic response of a building's individual structures.

[0032] Then, based on the building unit information, the building unit structure is simplified in the form of concentrated mass points at each floor. The mass point gravity is the representative value of the gravity load of the floor. The number of simplified mass points at each floor determines whether the building unit structure is a single-degree-of-freedom system or a multi-degree-of-freedom system. If there is only one mass point, it is a single-degree-of-freedom system; otherwise, it is a multi-degree-of-freedom system. The blast loads on both sides of the building unit structure are equal in magnitude and opposite in direction, and can offset each other. Simply superimpose the blast loads of the front and rear walls and load them onto the building unit structure to establish a simplified structural calculation model under the equivalent blast load. This simplified structural calculation model corresponds to both a single-degree-of-freedom structure calculation model and a multi-degree-of-freedom structure calculation model.

[0033] Then, based on the equivalent explosion load peak and the equivalent action time of the front wall shock wave overpressure, the stress conditions of the simplified structural calculation model are determined, and the corresponding motion equations are established. The motion equations correspond to the single-degree-of-freedom system motion equations and the multi-degree-of-freedom system motion equations.

[0034] For a single degree of freedom system, the motion equation of the single degree of freedom system is directly established as the following undamped motion equation:

[0035] ;

[0036] Where, m Indicates quality, represents acceleration, k represents the floor stiffness, u ( t ) represents displacement, P ( t ) represents the time-varying equivalent explosion load function.

[0037] The general solution plus special solution of the motion equation of the single degree of freedom system is used to calculate the dynamic response of the vertex of the single degree of freedom structural calculation model. After calculating the dynamic response of the vertex of the simplified structural calculation model, the dynamic response of the vertex is converted into the relationship between spectral displacement and spectral acceleration and plotted as ADRS The explosion spectrum of the structure is obtained by plotting the coordinate curve. The explosion spectrum is the relationship between the spectral displacement and spectral acceleration of the top mass point of the structure when the equivalent explosion load acts on the simplified structural calculation model. The spectral coordinate expression of the single degree of freedom system is:

[0038] The spectral shift is:

[0039] ;

[0040] Where, D ( t ) represents the single degree of freedom spectral displacement, P u represents the equivalent explosion load acting on the mass point at the floor, k represents the floor stiffness, ω represents the angular frequency, t u Indicates the total time that the blast load acts on the building, t Indicates the total time of explosion.

[0041] The spectral acceleration is:

[0042] ;

[0043] Where A(t) represents the single degree of freedom spectral acceleration, ω represents the angular frequency, D (t ) represents the single degree of freedom spectral displacement.

[0044] For multi-degree-of-freedom systems, the Pushover The analysis technology converts the multi-degree-of-freedom system into a single-degree-of-freedom system for explosion-resistant response evaluation. The vibration of the building's single structure is controlled by the first vibration mode. Regardless of the magnitude of the earthquake and the deformation of the building's single structure, the shape vector of the vibration mode remains unchanged. The motion equation of the multi-degree-of-freedom system is established:

[0045] ;

[0046] Where, M represents the mass matrix, C represents the damping matrix, K represents the stiffness matrix, represents the acceleration vector, represents the velocity vector, u ( t ) represents the displacement vector, P ( t ) represents the time-varying equivalent explosion load function.

[0047] According to the basic principle of vibration mode, the modal motion equation can be obtained:

[0048] ;

[0049] Where, M n Indicates the n The mass of the first mode, C n Indicates the n The damping of the first mode, K n Indicates the n The stiffness matrix of the first mode, Indicates the n The modal displacement vector of the order state, Indicates the n The modal velocity vector of the order, Indicates the n Modal acceleration vector of the order state.

[0050] According to the modal motion equation, the dynamic response of the vertex of the multi-degree-of-freedom structural calculation model is calculated. After the dynamic response of the vertex of the simplified structural calculation model is calculated, the dynamic response of the vertex is converted into the relationship between the spectral displacement and the spectral acceleration, and plotted as ADRS Coordinate curve diagram to obtain the structural explosion spectrum.

[0051] Specifically, the spectral coordinate expression of the multi-degree-of-freedom system is obtained according to the modal transformation and Duhamel integral formula, where the spectral displacement is:

[0052] ;

[0053] Where, Indicates the n The spectral displacement of the top particle of the hierarchical structure, Indicates the n The modal displacement vector of the order state, represents the modal mass of the top mass point of the n-th order structure, represents the angular frequency corresponding to the first vibration mode of the structure, represents the transpose of the first mode shape vector of the structure, P ( t ) represents the time-varying equivalent explosion load function, t Indicates the total time of explosion action, Represents the integral variable in the Duhamel integral.

[0054] The calculation formula of spectral acceleration is:

[0055] ;

[0056] Where, represents the spectral acceleration of the particle at the top of the multi-degree-of-freedom structure, represents the angular frequency corresponding to the first vibration mode of the structure, Indicates the n Spectral displacement of the top particle of the hierarchical structure.

[0057] Then, according to the seismic parameters and seismic code of the building's single structure, the seismic influence coefficient curve is obtained, and then according to the seismic influence coefficient curve, ADRS Convert to get the seismic spectrum, where

[0058] The spectral acceleration is: ;

[0059] Where, A (t) represents the earthquake demand spectrum acceleration, represents the earthquake influence coefficient, g Represents the acceleration due to gravity.

[0060] Spectral shift: ;

[0061] Where, D ( t ) represents the demand spectrum shift, T n represents the natural vibration period of the structure, represents pi, A (t) represents the earthquake demand spectrum acceleration.

[0062] Then, by adjusting the seismic parameters so that the apex of the explosion spectrum lies on the seismic curve, an equivalent seismic spectrum is obtained. This is a special seismic spectrum that is created by adjusting the seismic parameters so that the apex of the explosion spectrum lies on the seismic spectrum curve. The response of the equivalent seismic spectrum can represent the maximum displacement response of the structure to the blast load.

[0063] Then, according to the building unit information, use PKPM The elastoplastic static analysis module analyzes the capacity spectrum of the building's single structure. The capacity spectrum is the gradual overturning of the building's single structure by rectangular load loading. PKPM The capacity curve of the structure to resist lateral force obtained by elastoplastic static analysis is converted into the relationship between spectral displacement and spectral acceleration through the relationship between the top floor displacement and the base shear force of the building.

[0064] Finally, the explosion spectrum, equivalent earthquake spectrum, and capacity spectrum are plotted in the same coordinate system for comparative analysis, resulting in the explosion resistance of the building's individual structures and the equivalent earthquake effect of the explosion. For example, based on the comparative analysis of the explosion spectrum and capacity spectrum, when the explosion spectrum is below the capacity spectrum, it indicates that the structure's ability to resist lateral forces is greater than the lateral forces generated by the explosion load, and the structure is safe. When the explosion spectrum is above the capacity spectrum, it indicates that the structure's ability to resist lateral forces is less than the lateral forces generated by the explosion load, and the structure is unsafe, requiring structural adjustments. Based on the comparative analysis of the equivalent earthquake spectrum and capacity spectrum, the explosion resistance of the structure is determined by performance points.

[0065] Take the central control room of a project in Jingzhou, Hubei as an example. The central control room is a concrete frame structure with a construction area of ​​210.64 m 2 , 1 floor above ground, building length 18.5 m , width 11 m , height 7 m The representative value of the roof gravity load is 3148 kN , Y Directional stiffness 35860.47 kN / m According to the simulation analysis report of the overpressure consequences of steam cloud explosion, the occurrence of 1×10 -4 The project uses the Norwegian Classification Society ( DNV ) developed by Safeti 8.9" steam cloud explosion wave effect model, calculated the accident load and simulation results of explosion-resistant buildings under different meteorological conditions. The explosion overpressure of the central control room is 1 kPa , duration 12.7 ms , the extent of damage to the building can be assessed as the typical pressure at which glass breaks.

[0066] The present invention also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the above method when executing the computer program, thereby realizing visualization of explosion-resistant analysis results of a single structure.

[0067] The present invention further provides a computer-readable storage medium storing a computer program, which implements the steps of the above method when executed by a processor.

[0068] The above content is a further detailed description of the present invention in conjunction with specific preferred embodiments, and the specific implementation of the present invention should not be considered to be limited to these descriptions. For those skilled in the art of the present invention, several simple deductions and substitutions can be made without departing from the concept of the present invention, and all of these should be considered to fall within the scope of protection of the present invention.

Claims

1. A method for explosion-resistant analysis of a building structure, characterized in that: The steps include: S1. Obtain building unit information and explosion-proof information; S2. Calculating the equivalent overpressure action time of the shock wave and the equivalent explosion load peak value based on the building unit information and explosion resistance information; S3. Based on the building unit information, simplify the building unit structure in the form of concentrated mass points at each floor, and establish a simplified structural calculation model; S4. Determine the stress condition of the simplified structural calculation model according to the equivalent action time of the shock wave overpressure and the equivalent explosion load peak value, and establish the corresponding motion equation; S5. Calculate the dynamic response of the vertices of the simplified structural calculation model based on the motion equation to obtain an explosion spectrum; adjust the earthquake parameters so that the vertices of the explosion spectrum are located on the earthquake spectrum curve to obtain an equivalent earthquake spectrum; and obtain a capacity spectrum using an elastic-plastic static analysis method based on the building unit information; First, according to the earthquake influence coefficient response spectrum, the earthquake spectrum is obtained through ADRS transformation, and then the earthquake parameters are adjusted so that the peak of the explosion spectrum is located on the earthquake spectrum curve to obtain the equivalent earthquake spectrum; The explosion spectrum is the relationship between the spectral displacement and spectral acceleration of the top mass point of the structure when the equivalent explosion load acts on the simplified structural calculation model; According to the seismic parameters and seismic code of the building's single structure, the seismic influence coefficient curve is obtained, and then the seismic spectrum is obtained through ADRS conversion based on the seismic influence coefficient curve; By adjusting the seismic parameters so that the apex of the explosion spectrum is located on the seismic curve, an equivalent seismic spectrum is obtained. The equivalent seismic spectrum is a special seismic spectrum in which the apex of the explosion spectrum is located on the seismic spectrum curve after adjusting the seismic parameters. The response of the equivalent seismic spectrum can represent the maximum displacement response of the explosion load on the structure. The capacity spectrum is the capacity curve of the structure to resist lateral force obtained by the PKPM elastoplastic static analysis when the building single structure is gradually pushed over by rectangular load loading. The relationship between the top displacement and the base shear force of the building is converted into the relationship between the spectral displacement and the spectral acceleration. S6. Plot the explosion spectrum, equivalent earthquake spectrum, and capacity spectrum into the same coordinate system for comparative analysis to obtain the explosion resistance of the building's single structure and the earthquake effect generated by the explosion.

2. The explosion-resistant analysis method for a building monomer structure according to claim 1, characterized in that: In step S1, the building unit information is obtained through the PKPM structural design software, and the building unit information includes: the representative value of the particle gravity load at the floor, the floor height, the plane length of the building, the width of the building and the floor stiffness; the explosion-proof information is obtained through the project explosion safety assessment report, and the explosion-proof information includes: the building name, the location of the explosion, the peak incident overpressure of the explosion shock wave and the positive pressure action time of the explosion shock wave.

3. The explosion-resistant analysis method for a building monomer structure according to claim 2, characterized in that: In step S2, it is assumed that the explosion occurs directly in front of the building monomer structure, and the explosion loads borne by both sides of the building monomer structure are equal in magnitude and opposite in direction. Only the explosion forces of the front wall and the rear wall of the building monomer structure are considered, and the calculated shock wave overpressure equivalent action time is the shock wave overpressure equivalent action time of the front wall of the building monomer structure; the superposition of the explosion impulses of the front wall and the rear wall of the building monomer structure is also calculated, and based on the superposition of the explosion impulses of the front wall and the rear wall of the building monomer structure and the shock wave overpressure equivalent action time of the front wall of the building monomer structure, the equivalent explosion load peak is calculated in a constant impulse manner.

4. The explosion-resistant analysis method for a building monomer structure according to claim 3, characterized in that: In step S3, the simplified structural calculation model includes a single-degree-of-freedom structural calculation model and a multi-degree-of-freedom structural calculation model. In step S4, the motion equations include a single-degree-of-freedom system motion equation and a multi-degree-of-freedom system motion equation.

5. The explosion-resistant analysis method for a building monomer structure according to claim 4, characterized in that: In step S5, the single-degree-of-freedom system motion equation is solved by using the general solution plus the special solution to calculate the dynamic response of the vertex of the single-degree-of-freedom structural calculation model.

6. The explosion-resistant analysis method for a building monomer structure according to claim 4, characterized in that: In step S5, the motion equation of the multi-degree-of-freedom system adopts Pushover analysis technology and modal analysis method to first simplify the multi-degree-of-freedom structural calculation model into a single-degree-of-freedom structural calculation model, and then calculate the dynamic response of the vertex of the multi-degree-of-freedom structural calculation model.

7. The explosion-resistant analysis method for a building monomer structure according to claim 4, characterized in that: In step S5, after the dynamic response of the vertex of the simplified structural calculation model is calculated, the dynamic response of the vertex is converted into the relationship between spectral displacement and spectral acceleration, and plotted into an ADRS coordinate curve to obtain the structural explosion spectrum.

8. A computer device, characterized in that: The method comprises a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

9. A computer-readable storage medium, characterized in that A computer program is stored, and when the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images