Anti-explosion analysis method of building monomer structure, computer equipment and storage medium

Abstract

The invention belongs to the technical field of anti-explosion analysis of building structures, and particularly relates to an anti-explosion analysis method of a building monomer structure, computer equipment and a storage medium. According to the method, shock wave overpressure equivalent action time and an equivalent explosion load peak value are calculated through building monomer information and anti-explosion information; the method comprises the following steps: simplifying a building monomer structure in the form of concentrated mass points at floors, establishing a simplified structure calculation model and a corresponding motion equation, calculating to obtain dynamic response of a vertex, obtaining an explosion spectrum, adjusting seismic parameters to obtain an equivalent seismic spectrum, and obtaining a capability spectrum by adopting an elastic-plastic static analysis method. And drawing the explosion spectrum, the equivalent seismic spectrum and the capability spectrum into the same coordinate system for comparative analysis to obtain the anti-explosion performance of the building monomer structure and the seismic action equivalent to the explosion. The process of analyzing the anti-knock performance of the single building structure and the evaluation method of the anti-knock performance are simplified, the analysis efficiency of the anti-knock performance is improved, and parallel analysis of the anti-knock performance and the anti-knock performance is achieved.

Description

Technical Field

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

[0002] There are a large number of explosion hazards during the production and product storage of petrochemical enterprises. In actual projects, the explosion-proof analysis report in the chemical park clearly stipulates that buildings need to undergo explosion-proof analysis. The shock waves borne by buildings in office and living areas are less energy than those at the center of the explosion, and structural response analysis under explosion loads is often required. However, in the actual design process, explosion-proof analysis is too complicated.

[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, but 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 this regard, a Chinese patent application with application number CN202311295875.4 discloses an explosion response analysis method based on the modal decomposition response spectrum. This method only uses the modal decomposition response spectrum method to perform dynamic analysis on the structure, and cannot achieve parallel analysis of explosion performance and structural seismic performance analysis, 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 monomer structures, realize more rapid and effective explosion-resistant performance analysis of building monomer structures, and obtain the relationship between the explosion-resistant performance and seismic performance of the structure is a technical problem that needs to be urgently solved by technical personnel in this field. Summary of the invention

[0005] The present invention provides an explosion-proof analysis method for a building monomer structure, which simplifies the process of analyzing the explosion-proof performance of a building monomer structure and the evaluation method of the explosion-proof performance of a building monomer structure, improves the efficiency of analyzing the explosion-proof performance of a building monomer structure, and realizes the parallel analysis of the explosion-proof performance of a building monomer structure and the seismic performance analysis of a building monomer structure, thereby obtaining the relationship between the explosion-proof performance and the seismic performance of a building monomer structure.

[0006] The method comprises the following steps: S1. Obtain building monomer information and explosion-proof information; S2. Calculate the equivalent overpressure action time of the shock wave and the equivalent explosion load peak value according to the building monomer information and explosion resistance information; S3. According to the building monomer information, simplify the building monomer structure in the form of concentrated mass points at the floors, 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. According to the motion equation, the dynamic response of the vertex of the simplified structural calculation model is calculated to obtain an explosion spectrum; by adjusting the earthquake parameters, the vertex of the explosion spectrum is located on the earthquake general curve to obtain an equivalent earthquake spectrum; according to the building monomer information, the capacity spectrum is obtained by using the elastic-plastic static analysis method; S6. Plot the explosion spectrum, equivalent earthquake spectrum and capacity spectrum into the same coordinate system for comparative analysis, and obtain the explosion-resistant performance of the building's single structure and the earthquake effect equivalent to the explosion.

[0007] 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 positive pressure action time of the explosion shock wave.

[0008] In a specific embodiment, in step S2, it is preset that the explosion occurs directly in front of the building monomer structure, and the explosion loads borne on both sides of the building monomer structure are equal in magnitude and opposite in direction. Only the explosion forces on 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.

[0009] 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.

[0010] In a specific embodiment, in step S5, the single-degree-of-freedom system motion equation uses a general solution plus a special solution to calculate the dynamic response of the vertex of the single-degree-of-freedom structural calculation model.

[0011] In a specific embodiment, in step S5, the motion equation of the multi-degree-of-freedom system is adopted 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 vertices of the multi-degree-of-freedom structural calculation model is calculated.

[0012] 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.

[0013] 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.

[0014] The present invention also 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.

[0015] The present invention also 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.

[0016] The present invention simplifies the building monomer structure in the form of concentrated particles 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 value and the equivalent action time of the shock wave overpressure, and establishes a corresponding motion equation, calculates the dynamic response of the vertex of the simplified structural calculation model according to the motion equation, obtains the explosion spectrum according to the dynamic response of the vertex, and adjusts the earthquake parameters so that the vertex of the explosion spectrum is located on the earthquake general curve to obtain the equivalent earthquake spectrum, combines the explosion spectrum with the equivalent earthquake spectrum and the structural capacity spectrum in the same coordinate system for comparative analysis, obtains the explosion resistance of the building monomer structure and the earthquake action equivalent to the explosion, simplifies 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, improves the efficiency of analyzing the explosion resistance of the building monomer structure, realizes the parallel analysis of the explosion resistance of the building monomer structure and the seismic resistance of the building monomer structure, and obtains the relationship between the explosion resistance and the seismic resistance of the building monomer structure. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 is a flow chart of an embodiment of the present invention. DETAILED DESCRIPTION

[0018] 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 them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of this application.

[0019] See also Figure 1 The present invention provides a method for analyzing explosion resistance of a building monomer structure, which specifically comprises the following steps: First, through PKPM The structural design software obtains the building unit information of the building unit structure, including: the representative value of the mass point 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 of the building unit structure is obtained through the project explosion safety assessment report, including: the name of the building, 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.

[0020] Then, it is assumed that the explosion occurs in front of the building's single structure. Since the explosion loads on both sides of the building's single structure are equal in magnitude and opposite in direction, only the explosion forces on the front and back walls of the building's single structure are considered. According to the location of the explosion, the peak incident overpressure of the explosion shock wave, the positive pressure action time of the explosion shock wave, the floor height, the plane length and width of the building, the equivalent overpressure action time of the front wall of the building's single structure is calculated, as well as the impulse superposition of the front and back of each floor under the explosion, and the equivalent explosion load peak of the particle at the floor is calculated based on the constant impulse method. The calculation formula for the equivalent explosion load is: P u =BH ( P r t e +P b ( t rb +t d )) / t e ; In the formula, P u 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 overpressure of the front wall shock wave, P b represents the peak incident overpressure of the back wall shock wave, t rb represents the equivalent action time of the overpressure of the rear wall shock wave, t d represents the time required for the shock wave to propagate from the front wall to the rear wall, and / represents a division sign.

[0021] According to the calculated equivalent explosion load peak value and the calculated equivalent overpressure action time of the front wall shock wave of the building monomer structure, the equivalent explosion load time history function expression is established: P ( t ) =P u (1- t / t e ); In the formula, 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 action time of the front wall shock wave overpressure, and / represents the division sign. Through this time history function, the time history of the equivalent explosion load can be simulated more realistically, the dynamic response of the building monomer structure can be captured, and the dynamic response of the building monomer structure can be predicted more accurately.

[0022] Then, according to the building unit information, the building unit structure is simplified in the form of concentrated mass points at the floor, and the mass point gravity is the representative value of the gravity load of the floor. The number of simplified mass points at the floor is used to determine 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 explosion loads on both sides of the building unit structure are equal in magnitude and opposite in direction, and can offset each other. It is only necessary to superimpose the explosion loads of the front wall and the back wall and load them to the building unit structure to establish a simplified structural calculation model under the equivalent explosion load. The simplified structural calculation model corresponds to a single-degree-of-freedom structural calculation model and a multi-degree-of-freedom structural calculation model.

[0023] Then, according to the equivalent explosion load peak and the equivalent action time of the front wall shock wave overpressure, the stress condition of the simplified structural calculation model is determined, and the corresponding motion equation is established. The motion equation corresponds to the single-degree-of-freedom system motion equation and the multi-degree-of-freedom system motion equation.

[0024] 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: ; In the formula, m Indicates quality, represents acceleration, k represents the floor stiffness, u ( t ) represents displacement, P ( t ) represents the equivalent explosion load function that varies with time.

[0025] 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 the spectral displacement and the 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: The spectral shift is: ; In the formula, 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 represents the total time that the blast load acts on the building, t Indicates the total time of explosion.

[0026] The spectral acceleration is: ; 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.

[0027] For multi-degree-of-freedom systems, use Pushover The analysis technology converts the multi-degree-of-freedom system into a single-degree-of-freedom system for explosion-proof 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: ; In the formula,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 equivalent explosion load function that varies with time.

[0028] According to the basic principle of vibration mode, the modal motion equation can be obtained: ; In the formula, M n Indicates n The mass of the first mode, C n Indicates n The damping of the first mode, K n Indicates n The stiffness matrix of the order mode, Indicates n The modal displacement vector of the order, Indicates n The modal velocity vector of the order, Indicates n Modal acceleration vector of the order state.

[0029] 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.

[0030] 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: ; In the formula, Indicates n The spectral displacement of the top particle of the hierarchical structure, Indicates n The modal displacement vector of the order, 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 vibration mode vector of the structure, P ( t ) represents the equivalent explosion load function that varies with time, trepresents the total time of the explosion, Represents the integration variable in the Duhamel integral.

[0031] The calculation formula of spectral acceleration is: ; In the formula, 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 n Spectral displacement of the top particle of the hierarchical structure.

[0032] 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 The spectral acceleration is: ; In the formula, A (t) represents the earthquake demand spectrum acceleration, represents the earthquake influence coefficient, g Represents the acceleration due to gravity.

[0033] Spectral shift: ; In the formula, 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.

[0034] Then, by adjusting the earthquake parameters, the apex of the explosion spectrum is located on the earthquake general curve, and the equivalent earthquake spectrum is obtained. The equivalent earthquake spectrum is a special earthquake spectrum that adjusts the earthquake parameters so that the apex of the explosion spectrum is located on the earthquake spectrum curve. The response of the equivalent earthquake spectrum can represent the maximum displacement response of the explosion load on the structure.

[0035] 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 in the form of 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.

[0036] Finally, the explosion spectrum, equivalent earthquake spectrum and capacity spectrum are plotted in the same coordinate system for comparative analysis, and the explosion resistance of the building's single structure and the earthquake action equivalent to the explosion are obtained. For example, according to the comparative analysis of the explosion spectrum and the capacity spectrum, when the explosion spectrum is below the capacity spectrum, it means that the ability of the structure to resist lateral force is greater than the lateral force produced by the explosion load, and the structure is safe. When the explosion spectrum is above the capacity spectrum, it means that the ability of the structure to resist lateral force is less than the lateral force produced by the explosion load, and the structure is unsafe and needs to adjust the structural layout. According to the comparative analysis of the equivalent earthquake spectrum and the capacity spectrum, the explosion resistance of the structure is judged by the performance points.

[0037] 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 on the overpressure consequences of steam cloud explosion, a 1×10 -4 The project uses the Norwegian Classification Society (DNV) to identify hazard sources, calculate leakage, and simulate the consequences of vapor cloud explosions. DNV ) developed by Safeti 8.9” vapor 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 of glass breakage.

[0038] 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.

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

[0040] The above contents are further detailed descriptions of the present invention in combination with specific preferred embodiments, and it cannot be determined that the specific implementation of the present invention is limited to these descriptions. For ordinary technicians in the technical field to which the present invention belongs, several simple deductions and substitutions can be made without departing from the concept of the present invention, which should be regarded as belonging to the protection scope of the present invention.

Claims

1. A method for explosion-proof analysis of a building monomer structure, characterized in that: The steps include: S1. Obtain building monomer information and explosion-proof information; S2. Calculate the equivalent overpressure action time of the shock wave and the equivalent explosion load peak value according to the building monomer information and explosion resistance information; S3. According to the building unit information, the building unit structure is simplified in the form of concentrated mass points at the floor level, and a simplified structural calculation model is established; 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. According to the motion equation, the dynamic response of the vertex of the simplified structural calculation model is calculated to obtain an explosion spectrum; by adjusting the earthquake parameters, the vertex of the explosion spectrum is located on the earthquake general curve to obtain an equivalent earthquake spectrum; according to the building monomer information, the capacity spectrum is obtained by using the elastic-plastic static analysis method; S6. Plot the explosion spectrum, equivalent earthquake spectrum and capacity spectrum into the same coordinate system for comparative analysis, and obtain the explosion-resistant performance of the building's single structure and the earthquake effect equivalent to the explosion.

2. The explosion-resistant analysis method for a building monomer structure according to claim 1, characterized in that: 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 positive pressure action time of the explosion shock wave.

3. The explosion-resistant analysis method for a building monomer structure according to claim 2 is 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 on both sides of the building monomer structure are equal in magnitude and opposite in direction. Only the explosion forces on 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.

4. The explosion-resistant analysis method for a building monomer structure according to claim 3 is 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 is characterized in that: In step S5, the single-degree-of-freedom system motion equation uses a general solution plus a 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 multi-degree-of-freedom system motion equation is adopted 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 vertices of the multi-degree-of-freedom structural calculation model is calculated.

7. The explosion-resistant analysis method for a building monomer structure according to claim 4, characterized in that: 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.

8. The explosion-resistant analysis method for a building monomer structure according to claim 4, characterized in that: In step S5, firstly 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.

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

10. 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 8 are implemented.

Images