Structural anti-seismic performance analysis method and device considering cross-fault seismic oscillation, and terminal

Abstract

The invention is suitable for the technical field of seismic analysis, and provides a structure seismic performance analysis method and device considering cross-fault seismic oscillation, and a terminal. The method comprises the following steps: acquiring earthquake magnitude data, and calculating a recurrence period of each earthquake magnitude; establishing an empirical relationship between the magnitude and the recurrence period; calculating geometric dimension parameters of the fault through the empirical relation, and determining a research area crossing the fault and a fault fracture model; according to the fault fracture model, the digital elevation model and the soil layer information model, broadband seismic oscillation simulation is carried out, and an engineering site cross-fault seismic oscillation field is obtained; the method comprises the following steps: establishing a finite element model containing a cross-fault structure and surrounding rock, converting seismic oscillation into an equivalent force load on a boundary node of the finite element model, applying the equivalent force load on the boundary node of the finite element model, and performing cross-fault seismic response analysis on the cross-fault structure. According to the method, a closed-loop process from seismic oscillation selection input to performance index evaluation is realized, and a systematic solution is provided for anti-seismic performance evaluation analysis of a cross-fault structure.

Description

Technical Field

[0001] The present application relates to the technical field of earthquake analysis, and in particular to a method, device and terminal for analyzing the seismic performance of a structure taking into account cross-fault seismic motion. Background Art

[0002] The construction of cross-fault structures (such as cross-fault tunnels and bridges) is a key link in the modern transportation network and plays a vital role in rescue work after earthquake disasters. Therefore, it is particularly important to strengthen their seismic performance. Some areas are in high-seismic zones, and active faults and complex geological structures make the construction of cross-fault tunnels and bridges almost inevitable. However, the current seismic design standards for underground structures usually require the construction of tunnels and bridges to avoid active faults, but in actual projects, due to geographical restrictions, faults cannot be completely avoided, resulting in a large contradiction between the standards and actual engineering needs. Therefore, when constructing cross-fault tunnels and bridges in these high-risk areas, scientifically evaluating the seismic performance of cross-fault tunnels and bridges and improving their seismic toughness and safety are not only the key to ensuring project quality, but also a necessary measure to improve the seismic resistance of the overall transportation infrastructure.

[0003] In existing research, the seismic performance standard of cross-fault structures (such as cross-fault tunnels and bridges) is determined based on the seismic fortification zoning map, which provides the response spectrum parameters of the acceleration peak and characteristic period. The inventors of the present application found in their research that the existing standard response spectrum can only represent general seismic input and cannot meet the seismic input requirements of cross-fault structures. Summary of the invention

[0004] In view of this, the embodiments of the present application provide a method, device and terminal for analyzing the seismic performance of structures taking into account cross-fault seismic motions, realizing a complete method for evaluating the seismic performance of cross-fault structures from cross-fault seismic motion selection to quantitative evaluation of seismic performance, realizing a closed-loop process from seismic motion selection input to performance index evaluation, overcoming the limitations of traditional earthquake fortification zoning, and providing a systematic solution for the evaluation and analysis of the seismic performance of cross-fault structures.

[0005] In order to achieve the above objectives, this application adopts the following technical solutions: In a first aspect, an embodiment of the present application provides a method for analyzing seismic performance of a structure considering cross-fault ground motion, comprising: Obtain earthquake magnitude data caused by specific fault activity and calculate the return period of each magnitude; Fitting the earthquake magnitude data with the corresponding return period to establish an empirical relationship between magnitude and return period; The geometrical dimension parameters of the target fault are calculated by the empirical relationship between magnitude and recurrence period and the empirical relationship between source parameters, and the research area and fault rupture model across the target fault are determined according to the geometrical dimension parameters of the target fault; wherein the empirical relationship between source parameters represents the relationship between magnitude and the length and width of the fault rupture surface; Broadband seismic motion simulation is performed based on the fault rupture model, digital elevation model and soil layer information model of the study area to obtain the cross-fault seismic motion field of the engineering site; Establishing a finite element model including a cross-fault structure and surrounding rock, converting the simulated earthquake motion into an equivalent force load on the boundary nodes of the finite element model, and then applying the equivalent force load on the boundary nodes of the finite element model to perform a cross-fault seismic response analysis on the cross-fault structure; The damage of the cross-fault structure is evaluated based on the results of the cross-fault seismic response analysis, and the seismic performance of the cross-fault structure is determined.

[0006] Based on the first aspect, in some embodiments, obtaining earthquake magnitude data caused by specific fault activity and calculating the recurrence period of each magnitude includes: Select the active faults spanned by the cross-fault structure and count the magnitudes and occurrence times of all earthquakes in the area affected by the fault; Classify earthquake events according to magnitude range and calculate the recurrence period of each magnitude; The magnitude is divided into 0.5 levels. Calculate the return period of each magnitude, T is the return period, T total is the total time span of earthquake records, N is the total number of earthquake events that occurred during this time span.

[0007] Based on the first aspect, in some embodiments, fitting the earthquake magnitude data with the corresponding return period to establish an empirical relationship between magnitude and return period includes: pass lgT ( M )=- a t + B The earthquake magnitude data are fitted with the corresponding return period to obtain the best fitting curve representing the relationship between magnitude and return period. T ( M ) is a certain magnitude M Or the average interval between earthquakes above t is the total length of time for statistics, a t and b is an empirical constant.

[0008] Based on the first aspect, in some embodiments, the geometric size parameters of the target fault are calculated by using the empirical relationship between magnitude and recurrence period and the preset empirical relationship between source parameters, and the research area and fault rupture model across the target fault are determined according to the geometric size parameters of the target fault, including: Determine the magnitude values ​​corresponding to extremely rare, rare, fortified and frequent occurrence of the target fault according to the best fitting curve, and calculate the length and width of the rupture surface of the target fault according to the magnitude value and the preset empirical relationship between the source parameters; The epicenter location, focal depth, burial depth, sliding angle and dip angle are searched through literature data, and the fault rupture model and study area are determined based on the epicenter location, focal depth, burial depth, sliding angle and dip angle, as well as the length and width of the target fault rupture surface.

[0009] Based on the first aspect, in some embodiments, performing broadband seismic simulation according to the fault rupture model, digital elevation model and soil layer information model of the study area to obtain the cross-fault seismic field of the engineering site includes: Combining the elevation dataset and the surrounding rock unified seismic tomography model USTClitho2.0, a digital elevation model and soil layer information model of the study area were constructed; The fault rupture model, the soil layer information model and the digital elevation model are integrated, and based on the integrated model, a broadband seismic motion simulation is performed on the study area using computing software to obtain the cross-fault seismic motion field in the study area.

[0010] Based on the first aspect, in some embodiments, establishing a finite element model including a cross-fault structure and surrounding rock, and converting the simulated ground motion into an equivalent effective load on the boundary nodes of the finite element model includes: A finite element model including the cross-fault structure and its surrounding rock was established in ANSYS, and viscoelastic artificial boundaries were added to the five surfaces of the finite element model except the ground surface; The cross-fault seismic field data extracted from the cross-fault seismic simulation model are converted into equivalent forces on the boundary nodes of the finite element model using the area reduction method.

[0011] Based on the first aspect, in some embodiments, converting the cross-fault seismic field data extracted from the cross-fault seismic simulation model into equivalent forces on the boundary nodes of the finite element model using a region reduction method includes: pass The cross-fault seismic field data is converted into equivalent forces on the boundary nodes of the finite element model. For equal effect, M , C and K are the mass matrix, damping matrix and stiffness matrix respectively,e and b They are respectively + and nodes on region Γ, where region Γ represents the internal and external intersection boundaries of the interaction model of the local area intercepted in the study area, and region Ω + represents the external region of the interaction model of the local region, Represents Ω + The acceleration field in the region, Represents Ω + The velocity field of the region, Represents Ω + The displacement field of the region, represents the acceleration field in region Γ, represents the velocity field in the Γ region, represents the displacement field in the Γ region, Represents area Ω + The mass submatrix stored in the mass matrix of is in the order of node b and node e. Represents area Ω + The mass submatrix stored in the mass matrix according to the order of node e and node b, Represents area Ω + The damping submatrix stored in the damping matrix according to the order of node b and node e, Represents area Ω + The damping submatrix stored in the damping matrix according to the order of node e and node b, Represents area Ω + The stiffness submatrix stored in the stiffness matrix according to the order of node b and node e, Represents area Ω + The stiffness submatrix stored in the stiffness matrix according to the order of node e and node b.

[0012] Based on the first aspect, in some embodiments, the cross-fault structure includes a cross-fault tunnel and a bridge, and the evaluating the damage of the cross-fault structure according to the cross-fault seismic response analysis result and determining the seismic performance of the cross-fault structure includes: Calculating the evaluation index value of the cross-fault tunnel or bridge according to the response of the finite element model of the cross-fault tunnel or bridge under earthquake excitation, the evaluation index includes one or more of the relative diameter deformation rate, the diameter deformation rate and the tunnel inclination angle, the displacement ductility ratio of the bridge support, and the curvature ductility ratio of the bridge pier; The evaluation index value is compared with the evaluation index threshold to determine the seismic performance of the cross-fault tunnel or bridge.

[0013] In a second aspect, an embodiment of the present application provides a structural seismic performance analysis device considering cross-fault ground motion, comprising: A data acquisition module is used to obtain earthquake magnitude data caused by specific fault activities and calculate the recurrence period of each magnitude; A fitting module, used to fit the earthquake magnitude data with the corresponding return period to establish an empirical relationship between the magnitude and the return period; A regional model determination module is used to calculate the geometric size parameters of the target fault through the empirical relationship between the magnitude and the recurrence period and the empirical relationship between the source parameters, and determine the research area and the fault rupture model across the target fault according to the geometric size parameters of the target fault; wherein the empirical relationship of the source parameters represents the relationship between the magnitude and the length and width of the fault rupture surface; The vibration simulation module is used to simulate broadband seismic motion based on the fault rupture model, digital elevation model and soil layer information model of the study area to obtain the cross-fault seismic field of the engineering site; An analysis module is used to establish a finite element model including a cross-fault structure and surrounding rock, convert the simulated earthquake motion into an equivalent force load on the boundary nodes of the finite element model, and then apply the equivalent force load to the boundary nodes of the finite element model to perform a cross-fault seismic response analysis on the cross-fault structure; The seismic performance determination module is used to evaluate the lining damage of the cross-fault structure based on the cross-fault seismic response analysis results and determine the seismic performance of the cross-fault structure.

[0014] In a third aspect, an embodiment of the present application provides a terminal comprising a memory and a processor, wherein the memory stores a computer program that can be executed on the processor, and when the processor executes the computer program, it implements the structural seismic performance analysis method considering cross-fault seismic motion as described in any one of the first aspects.

[0015] Compared with the prior art, the embodiments of the present invention have the following beneficial effects: The embodiments of the present application establish a series of cross-fault structure seismic performance evaluation standards, which can systematically analyze the seismic response of cross-fault structures under multi-level performance targets such as slight damage, moderate damage and severe damage, and realize the seismic performance evaluation of cross-fault structures under different magnitude inputs. In addition, based on the relationship between magnitude and recurrence period and the construction of cross-fault seismic motion input, the present application realizes a complete cross-fault structure seismic performance evaluation method from cross-fault seismic motion selection to quantitative evaluation of seismic performance, and realizes a closed-loop process from seismic motion selection input to performance index evaluation, overcoming the limitations of traditional earthquake fortification zoning, and providing a systematic solution for the evaluation and analysis of seismic performance of cross-fault structures. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0017] Figure 1 It is a flow chart of a method for analyzing seismic performance of a structure considering cross-fault ground motion provided in an embodiment of the present application; Figure 2 is an empirical relationship curve between magnitude and return period provided in the embodiment of the present application; Figure 3 is a schematic diagram of the structure of a structural seismic performance analysis device considering cross-fault ground motion provided in an embodiment of the present application; Figure 4 It is a schematic diagram of a terminal provided in an embodiment of the present application. DETAILED DESCRIPTION

[0018] The present application is described more clearly below in conjunction with specific embodiments. The following embodiments will help those skilled in the art to further understand the effects of the present application, but are not intended to limit the present application in any form. It should be noted that, for those of ordinary skill in the art, several variations and improvements may be made without departing from the concept of the present application. These all fall within the scope of protection of the present application.

[0019] In order to make the purpose, technical solutions and advantages of the present application clearer, specific embodiments will be described below in conjunction with the accompanying drawings.

[0020] Traditional earthquake fortification standards are usually based on earthquake fortification zoning maps, which provide seismic motions (such as peak acceleration and characteristic period) in specific areas. These response spectra are suitable for general seismic motion inputs, but cannot reflect the seismic motion characteristics unique to cross-fault structures (such as cross-fault tunnels and bridges) (cross-fault seismic motions are often not uniform, especially when the fault is displaced or slipped, the direction, frequency and intensity of the seismic motion may be asymmetrically distributed). Different magnitudes, source locations, fault characteristics and other factors will affect the characteristics of seismic motion. Therefore, traditional response spectra cannot meet the seismic motion input requirements of cross-fault structures.

[0021] In order to more accurately simulate the seismic motion experienced by cross-fault structures, this application introduces the relationship between magnitude and return period (i.e., magnitude and frequency). Through historical earthquake data research, the frequency of occurrence of earthquakes of different magnitudes (such as M5.0, M6.0, etc.) (i.e., the distribution of magnitude) and the seismic motion characteristics under different magnitudes can be obtained, providing a more accurate seismic input for the seismic analysis of cross-fault structures than the standard response spectrum in traditional methods. Through the relationship between magnitude and return period, the seismic input that conforms to the characteristics of seismic activity in the area can be selected based on the historical earthquake data at the location of the project. According to the above relationship between magnitude and return period, the seismic input under different magnitudes can be obtained, which does not rely on a certain standard seismic motion, but comprehensively considers earthquakes of different magnitudes, so as to better capture the possible seismic motion characteristics.

[0022] The implementation scheme of the present application is described in detail below with reference to the accompanying drawings.

[0023] See also Figure 1 The structural seismic performance analysis method considering cross-fault ground motion provided in the embodiment of the present application may include the following steps: Step 101, obtaining earthquake magnitude data caused by specific fault activity, and calculating the recurrence period of each magnitude.

[0024] In some embodiments, step 101 may include: selecting an active fault spanned by a cross-fault structure, and counting the magnitudes and occurrence times of all earthquakes in the area affected by the fault; classifying earthquake events according to the magnitude range, and calculating the recurrence period of each magnitude; dividing the magnitude into 0.5 levels, and Calculate the return period of each magnitude, T is the return period, T total is the total time span of earthquake records, N is the total number of earthquake events that occurred during this time span.

[0025] For example, for the Menyuan area (located at the junction of Province A and Province B), earthquake data from a specific area (a fault near the Menyuan area) is collected, as shown in Table 1: Table 1 Information on earthquakes that occurred in Menyuan area from 1927 to 2017

[0026] Divide the magnitude into 0.5 levels and use the formula The recurrence period of each magnitude bin is calculated and the results are shown in Table 2.

[0027] Table 2 Earthquake recurrence period in Menyuan area

[0028] Step 102, fitting the earthquake magnitude data with the corresponding return period to establish an empirical relationship between magnitude and return period.

[0029] In some embodiments, step 102 may include: lgT ( M )=- a t + B + lgt The earthquake magnitude data are fitted with the corresponding return period to obtain the best fitting curve representing the relationship between magnitude and return period. T ( M ) is a certain magnitude M Or the average interval between earthquakes above t is the total length of time for statistics, a t and b is an empirical constant, a t and b It can be estimated by the least squares method.

[0030] For example, taking the data in Table 2 as an example, lgT ( M )=- a t + B Establish an empirical relationship between magnitude and return period, and obtain Figure 2 The empirical relationship curve is shown. Figure 2 In the figure, the horizontal axis is the return period (in years) and the vertical axis is the magnitude M.

[0031] Step 103, calculating the geometrical dimension parameters of the target fault through the empirical relationship between magnitude and recurrence period and the empirical relationship between source parameters, and determining the research area and fault rupture model across the target fault according to the geometrical dimension parameters of the target fault.

[0032] The empirical relationship of the earthquake source parameters represents the relationship between the magnitude and the length and width of the fault rupture surface. Specifically, the empirical relationship of the earthquake source parameters indicates that the length and width of the fault rupture surface are inferred from the magnitude of the earthquake. A formula of the empirical relationship of the earthquake source parameters is log (SRL) = c + d × M, where SRL is the length of the fault surface rupture, M is the magnitude, and c and d are empirical parameters.

[0033] In some embodiments, step 103 may include: determining the magnitude values ​​corresponding to extremely rare, rare, fortified and frequent occurrences of the target fault according to the best fitting curve, and calculating the length and width of the target fault rupture surface according to the magnitude values ​​and the preset empirical relationship between the source parameters; searching the epicenter location, source depth, burial depth, sliding angle and inclination through literature data, and determining the fault rupture model and study area according to the epicenter location, source depth, burial depth, sliding angle and inclination, and the length and width of the target fault rupture surface.

[0034] In the embodiments of the present application, extremely rare refers to an earthquake with a magnitude of 0.5% exceeding probability in 50 years, rare refers to an earthquake with a magnitude of 2% exceeding probability in 50 years, defensive refers to an earthquake with a magnitude of 10% exceeding probability in 50 years, and frequent refers to an earthquake with a magnitude of 63% exceeding probability in 50 years.

[0035] Step 104 , performing broadband seismic motion simulation based on the fault rupture model, digital elevation model and soil layer information model of the study area to obtain the cross-fault seismic motion field of the engineering site.

[0036] In some embodiments, step 104 may include: combining the elevation data set and the surrounding rock unified seismic tomography model USTClitho2.0 to construct a digital elevation model and a soil layer information model of the study area; integrating the fault rupture model, the soil layer information model and the digital elevation model, and performing a broadband seismic motion simulation on the study area using computing software based on the integrated model to obtain the cross-fault seismic motion field of the engineering site (i.e., the study area).

[0037] Step 105, establish a finite element model including the cross-fault structure and surrounding rock, convert the simulated seismic motion into an equivalent force load on the boundary nodes of the finite element model, and then apply the equivalent force load to the boundary nodes of the finite element model to perform a cross-fault seismic response analysis on the cross-fault structure.

[0038] In some embodiments, "establishing a finite element model including a cross-fault structure and surrounding rock, and converting the simulated seismic motion into equivalent force loads on the boundary nodes of the finite element model" in step 105 may include: establishing a finite element model including a cross-fault structure and its surrounding rock in ANSYS, adding viscoelastic artificial boundaries on five surfaces other than the ground surface in the finite element model; converting the cross-fault seismic field data extracted from the cross-fault seismic motion simulation model into equivalent forces on the boundary nodes of the finite element model using the area reduction method. After the equivalent forces on the boundary nodes of the finite element model, the equivalent forces on the boundary nodes of the finite element model are applied to the boundary nodes of the finite element model to perform cross-fault seismic response analysis on the cross-fault structure.

[0039] In this embodiment, the improved region reduction method MDRM can be used to convert the cross-fault seismic field data extracted from the cross-fault seismic simulation model into equivalent forces on the boundary nodes of the finite element model.

[0040] Specifically, the above-mentioned conversion of the cross-fault seismic field data extracted from the cross-fault seismic simulation model into equivalent forces on the boundary nodes of the finite element model using the area reduction method may include: pass The cross-fault seismic field data is converted into equivalent forces on the boundary nodes of the finite element model. For equal effect, M , C and K are the mass matrix, damping matrix and stiffness matrix respectively, e and b They are respectively + and nodes on region Γ, where region Γ represents the internal and external intersection boundaries of the interaction model of the local area intercepted in the study area, and region Ω + represents the external region of the interaction model of the local region, Represents Ω + The acceleration field in the region, Represents Ω + The velocity field of the region, Represents Ω + The displacement field of the region, represents the acceleration field in region Γ, represents the velocity field in the Γ region, represents the displacement field in the Γ region, Represents area Ω + The mass submatrix stored in the mass matrix of is in the order of node b and node e. Represents area Ω + The mass submatrix stored in the mass matrix according to the order of node e and node b, Represents area Ω + The damping submatrix stored in the damping matrix according to the order of node b and node e, Represents area Ω + The damping submatrix stored in the damping matrix according to the order of node e and node b, Represents area Ω + The stiffness submatrix stored in the stiffness matrix according to the order of node b and node e, Represents area Ω + The stiffness submatrix stored in the stiffness matrix according to the order of node e and node b.

[0041] Step 106, evaluating the damage of the lining of the cross-fault structure according to the cross-fault seismic response analysis results, and determining the seismic performance of the cross-fault structure.

[0042] In some embodiments, the cross-fault structure may include a cross-fault tunnel and a bridge, and step 106 may include: calculating the evaluation index value of the cross-fault tunnel or bridge according to the response of the finite element model of the cross-fault tunnel or bridge under seismic excitation, the evaluation index including one or more of the relative diameter deformation rate, the diameter deformation rate and the tunnel inclination angle, the displacement ductility ratio of the bridge support, and the curvature ductility ratio of the pier; comparing the evaluation index value with the evaluation index threshold to determine the seismic performance of the cross-fault tunnel or bridge.

[0043] In this embodiment, the seismic performance of cross-fault tunnels and bridges can be evaluated by setting performance evaluation indicators, and the performance evaluation indicators can include relative diameter deformation rate, diameter deformation rate and tunnel inclination angle, displacement ductility ratio of bridge supports, curvature ductility ratio of bridge piers, etc. Among them, relative diameter deformation rate and diameter deformation rate, displacement ductility ratio of bridge supports, and curvature ductility ratio of bridge piers are all used to evaluate the degree of damage to tunnels and bridges, and are seismic performance indicators of tunnels and bridges. The difference between the two is that the relative diameter deformation rate is more suitable for evaluating the seismic performance of horseshoe-shaped tunnels, the diameter deformation rate is more suitable for evaluating the seismic performance of circular tunnels, and the displacement ductility ratio of bridge supports and curvature ductility ratio of bridge piers are more suitable for evaluating the seismic performance of simply supported beam bridges.

[0044] For example, for horseshoe tunnels, the relative diameter deformation rate is usually used as a performance evaluation index. Specifically, when the relative diameter deformation rate of the horseshoe tunnel reaches 0.6‰, the tunnel structure is considered to be in a state of slight damage; if the relative diameter deformation rate of the horseshoe tunnel reaches 2.5‰, the tunnel structure is considered to be in a state of moderate damage; and when the relative diameter deformation rate of the horseshoe tunnel reaches 4‰, the tunnel structure is considered to be in a state of severe damage. For simply supported beam bridges, the displacement ductility ratio of the bearing and the curvature ductility ratio of the pier are usually used as performance evaluation indicators. Specifically, when the curvature ductility ratio of the pier of a simply supported beam bridge reaches 1.259, the bridge structure is considered to be in a state of slight damage; if the curvature ductility ratio of the pier of a simply supported beam bridge reaches 1.429, the bridge structure is considered to be in a state of moderate damage; and when the curvature ductility of the pier of a simply supported beam bridge reaches 2.444, the bridge structure is considered to be in a state of severe damage.

[0045] In this way, the degree of structural damage can be determined from the perspective of structural deformation rather than structural damage, thereby analyzing the seismic performance of cross-fault tunnels and bridges.

[0046] This application combines the improved regional reduction method MDRM to analyze the seismic response of cross-fault structures, and establishes a series of cross-fault structure seismic performance evaluation standards, which can systematically analyze the seismic response of cross-fault structures under multi-level performance targets such as slight damage, moderate damage and severe damage, and realize the seismic performance evaluation of cross-fault structures under different magnitude inputs. In addition, based on the relationship between magnitude and recurrence period and the construction of cross-fault seismic motion input, this application realizes a complete cross-fault structure seismic performance evaluation method from cross-fault seismic motion selection to quantitative evaluation of seismic performance, and realizes a closed-loop process from seismic motion selection input to performance index evaluation, overcoming the limitations of traditional earthquake fortification zoning, and providing a systematic solution for the evaluation and analysis of seismic performance of cross-fault structures.

[0047] It should be understood that the size of the serial numbers of the steps in the above embodiments does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present application.

[0048] Corresponding to the structural seismic performance analysis method considering cross-fault ground motion described in the above embodiment, Figure 3 A structural block diagram of a structural seismic performance analysis device considering cross-fault seismic motion provided in an embodiment of the present application is shown. For ease of explanation, only the parts related to the embodiment of the present application are shown.

[0049] See also Figure 3 An embodiment of the present application provides a structural seismic performance analysis device taking into account cross-fault ground motions, including a data acquisition module 201, a fitting module 202, a regional model determination module 203, a vibration simulation module 204, an analysis module 205 and a seismic performance determination module 206.

[0050] Specifically, the data acquisition module 201 is used to acquire earthquake magnitude data caused by specific fault activities and calculate the recurrence period of each magnitude.

[0051] The fitting module 202 is used to fit the earthquake magnitude data with the corresponding return period to establish an empirical relationship between the magnitude and the return period.

[0052] The regional model determination module 203 is used to calculate the geometric size parameters of the target fault through the empirical relationship between magnitude and recurrence period and the empirical relationship between source parameters, and determine the research area and fault rupture model across the target fault according to the geometric size parameters of the target fault. The empirical relationship between source parameters represents the relationship between magnitude and the length and width of the fault rupture surface.

[0053] The vibration simulation module 204 is used to perform broadband seismic simulation based on the fault rupture model, digital elevation model and soil layer information model of the study area to obtain the cross-fault seismic field of the engineering site.

[0054] The analysis module 205 is used to establish a finite element model including the cross-fault structure and surrounding rock, convert the simulated seismic motion into an equivalent force load on the boundary nodes of the finite element model, and then apply the equivalent force load to the boundary nodes of the finite element model to perform a cross-fault seismic response analysis on the cross-fault structure.

[0055] The seismic performance determination module 206 is used to evaluate the lining damage of the cross-fault structure according to the cross-fault seismic response analysis results and determine the seismic performance of the cross-fault structure.

[0056] Optionally, the data acquisition module 201 is specifically used to: select active faults spanned by cross-fault structures, and count the magnitudes and occurrence times of all earthquakes in the area affected by the faults; classify earthquake events according to magnitude ranges, and calculate the recurrence period of each magnitude; divide the magnitude into 0.5-level bins, and Calculate the return period of each magnitude, T is the return period, T total is the total time span of earthquake records, N is the total number of earthquake events that occurred during this time span.

[0057] Optionally, the fitting module 202 is specifically used to: lgT ( M )=- a t + B The earthquake magnitude data are fitted with the corresponding return period to obtain the best fitting curve representing the relationship between magnitude and return period. T ( M ) is a certain magnitude M Or the average interval between earthquakes above t is the total length of time for statistics, a t and b is an empirical constant.

[0058] Optionally, the regional model determination module 203 is specifically used to: determine the magnitude values ​​corresponding to extremely rare, rare, fortified and frequent occurrences of the target fault according to the best fitting curve, and calculate the length and width of the target fault rupture surface according to the magnitude value and the preset empirical relationship between the source parameters; search the epicenter position, source depth, burial depth, sliding angle and inclination through literature data, and determine the fault rupture model and study area according to the epicenter position, source depth, burial depth, sliding angle and inclination, and the length and width of the target fault rupture surface.

[0059] Optionally, the vibration simulation module 204 is specifically used to: combine the elevation data set and the unified seismic tomography model of surrounding rock USTClitho2.0 to construct a digital elevation model and a soil layer information model of the study area; integrate the fault rupture model, the soil layer information model and the digital elevation model, and use computing software to perform broadband seismic motion simulation on the study area based on the integrated model to obtain the cross-fault seismic motion field in the study area.

[0060] Optionally, the analysis module 205 is specifically used to: establish a finite element model including a cross-fault structure and its surrounding rock in ANSYS, and add viscoelastic artificial boundaries on five surfaces other than the ground surface in the finite element model; use the area reduction method to convert the cross-fault seismic field data extracted from the cross-fault seismic motion simulation model into equivalent forces on the boundary nodes of the finite element model.

[0061] Exemplarily, the analysis module 205 may specifically be The cross-fault seismic field data is converted into equivalent forces on the boundary nodes of the finite element model. For equal effect, M , C and K are the mass matrix, damping matrix and stiffness matrix respectively, e and b They are respectively + and nodes on region Γ, where region Γ represents the internal and external intersection boundaries of the interaction model of the local area intercepted in the study area, and region Ω + represents the external region of the interaction model of the local region, Represents Ω + The acceleration field in the region, Represents Ω + The velocity field of the region, Represents Ω + The displacement field of the region, represents the acceleration field in region Γ, represents the velocity field in the Γ region, represents the displacement field in the Γ region, Represents area Ω + The mass submatrix stored in the mass matrix of is in the order of node b and node e. Represents area Ω + The mass submatrix stored in the mass matrix according to the order of node e and node b, Represents area Ω + The damping submatrix stored in the damping matrix according to the order of node b and node e, Represents area Ω +The damping submatrix stored in the damping matrix according to the order of node e and node b, Represents area Ω + The stiffness submatrix stored in the stiffness matrix according to the order of node b and node e, Represents area Ω + The stiffness submatrix stored in the stiffness matrix according to the order of node e and node b.

[0062] Optionally, the cross-fault structure includes a cross-fault tunnel and a bridge, and the seismic performance determination module 206 is specifically used to: calculate the evaluation index value of the cross-fault tunnel or bridge according to the response of the finite element model of the cross-fault tunnel or bridge under seismic excitation, the evaluation index including one or more of the relative diameter deformation rate, the diameter deformation rate and the tunnel inclination angle, the support displacement ductility ratio of the bridge, and the pier curvature ductility ratio; compare the evaluation index value with the evaluation index threshold to determine the seismic performance of the cross-fault tunnel or bridge.

[0063] Figure 4 FIG. 1 is a schematic diagram of a terminal provided by an embodiment of the present invention. Figure 4 As shown, the terminal 300 of this embodiment includes a processor 310 and a memory 320, wherein the memory 320 stores a computer program that can be run on the processor 310, such as a structural seismic performance analysis program considering cross-fault ground motion. When the processor 310 executes the computer program, the steps in the above-mentioned structural seismic performance analysis method embodiment considering cross-fault ground motion are implemented, such as Figure 1 Alternatively, when the processor 310 executes the computer program, the functions of the modules in the above-mentioned device embodiments are implemented, for example Figure 3 The functions of the data acquisition module 201 to the seismic performance determination module 206 are shown.

[0064] The terminal 300 may be a computing device such as a desktop computer, a notebook, a PDA, or a cloud server. The terminal may include, but is not limited to, a processor 310 and a memory 320. Those skilled in the art will appreciate that Figure 4 It is only an example of terminal 300 and does not constitute a limitation of terminal 300. It may include more or fewer components than shown in the figure, or a combination of certain components, or different components. For example, the terminal may also include input and output devices, network access devices, buses, etc.

[0065] The processor 310 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.

[0066] The memory 320 may be an internal storage unit of the terminal 300, such as a hard disk or memory of the terminal 300. The memory 320 may also be an external storage device of the terminal 300, such as a plug-in hard disk, a smart media card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the terminal 300. Further, the memory 320 may also include both an internal storage unit of the terminal 300 and an external storage device. The memory 320 is used to store the computer program and other programs and data required by the terminal. The memory 320 may also be used to temporarily store data that has been output or is to be output.

[0067] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0068] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0069] The embodiments described above are only used to illustrate the technical solutions of the present invention, rather than to limit the same. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that the technical solutions described in the aforementioned embodiments may still be modified, or some of the technical features may be replaced by equivalents. Such modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be included in the protection scope of the present invention.

Claims

1. A method for analyzing the seismic performance of a structure considering cross-fault ground motion, characterized in that: include: Obtain earthquake magnitude data caused by specific fault activity and calculate the return period of each magnitude; Fitting the earthquake magnitude data with the corresponding return period to establish an empirical relationship between magnitude and return period; The geometrical dimension parameters of the target fault are calculated by the empirical relationship between magnitude and recurrence period and the empirical relationship between source parameters, and the research area and fault rupture model across the target fault are determined according to the geometrical dimension parameters of the target fault; wherein the empirical relationship between source parameters represents the relationship between magnitude and the length and width of the fault rupture surface; Broadband seismic motion simulation is performed based on the fault rupture model, digital elevation model and soil layer information model of the study area to obtain the cross-fault seismic motion field of the engineering site; Establishing a finite element model including a cross-fault structure and surrounding rock, converting the simulated earthquake motion into an equivalent force load on the boundary nodes of the finite element model, and then applying the equivalent force load on the boundary nodes of the finite element model to perform a cross-fault seismic response analysis on the cross-fault structure; The damage of the cross-fault structure is evaluated based on the results of the cross-fault seismic response analysis, and the seismic performance of the cross-fault structure is determined.

2. The structural seismic performance analysis method considering cross-fault ground motion according to claim 1 is characterized in that: The method of obtaining earthquake magnitude data caused by specific fault activity and calculating the recurrence period of each magnitude includes: Select the active faults spanned by the cross-fault structure and count the magnitudes and occurrence times of all earthquakes in the area affected by the fault; Classify earthquake events according to magnitude range and calculate the recurrence period of each magnitude; The magnitude is divided into 0.5 levels. Calculate the return period of each magnitude, T is the return period, T total is the total time span of earthquake records, N is the total number of earthquake events that occurred during this time span.

3. The structural seismic performance analysis method considering cross-fault ground motion according to claim 1 is characterized in that: The step of fitting the earthquake magnitude data with the corresponding return period to establish an empirical relationship between magnitude and return period includes: pass lgT ( M )=- a t + B The earthquake magnitude data are fitted with the corresponding return period to obtain the best fitting curve representing the relationship between magnitude and return period. T ( M ) is a certain magnitude M Or the average interval between earthquakes above t is the total length of time for statistics, a t and b is an empirical constant.

4. The structural seismic performance analysis method considering cross-fault ground motion according to claim 3 is characterized in that: The method of calculating the geometric size parameters of the target fault by using the empirical relationship between magnitude and recurrence period and the preset empirical relationship between earthquake source parameters, and determining the research area and fault rupture model across the target fault according to the geometric size parameters of the target fault, includes: Determine the magnitude values ​​corresponding to extremely rare, rare, fortified and frequent occurrence of the target fault according to the best fitting curve, and calculate the length and width of the rupture surface of the target fault according to the magnitude value and the preset empirical relationship between the source parameters; The epicenter location, focal depth, burial depth, sliding angle and dip angle are searched through literature data, and the fault rupture model and study area are determined based on the epicenter location, focal depth, burial depth, sliding angle and dip angle, as well as the length and width of the target fault rupture surface.

5. The structural seismic performance analysis method considering cross-fault ground motion according to claim 4 is characterized in that: The broadband seismic motion simulation is performed based on the fault rupture model, digital elevation model and soil layer information model of the study area to obtain the cross-fault seismic motion field of the engineering site, including: Combining the elevation dataset and the surrounding rock unified seismic tomography model USTClitho2.0, a digital elevation model and soil layer information model of the study area were constructed; The fault rupture model, the soil layer information model and the digital elevation model are integrated, and based on the integrated model, a broadband seismic motion simulation is performed on the study area using computing software to obtain the cross-fault seismic motion field in the study area.

6. The structural seismic performance analysis method considering cross-fault ground motion according to claim 1 is characterized in that: The step of establishing a finite element model including the cross-fault structure and surrounding rock, and converting the simulated ground motion into an equivalent effective load on the boundary nodes of the finite element model, includes: A finite element model including the cross-fault structure and its surrounding rock was established in ANSYS, and viscoelastic artificial boundaries were added to the five surfaces of the finite element model except the ground surface; The cross-fault seismic field data extracted from the cross-fault seismic simulation model are converted into equivalent forces on the boundary nodes of the finite element model using the area reduction method.

7. The structural seismic performance analysis method considering cross-fault ground motion according to claim 6, characterized in that: The cross-fault seismic field data extracted from the cross-fault seismic simulation model is converted into equivalent forces on the boundary nodes of the finite element model using a region reduction method, including: pass The cross-fault seismic field data is converted into equivalent forces on the boundary nodes of the finite element model. For equal effect, M , C and K are the mass matrix, damping matrix and stiffness matrix respectively, e and b They are respectively + and nodes on region Γ, where region Γ represents the internal and external intersection boundaries of the interaction model of the local area intercepted in the study area, and region Ω + represents the external region of the interaction model of the local region, Represents Ω + The acceleration field in the region, Represents Ω + The velocity field of the region, Represents Ω + The displacement field of the region, represents the acceleration field in region Γ, represents the velocity field in the Γ region, represents the displacement field in the Γ region, Represents area Ω + The mass submatrix stored in the mass matrix of is in the order of node b and node e. Represents area Ω + The mass submatrix stored in the mass matrix according to the order of node e and node b, Represents area Ω + The damping submatrix stored in the damping matrix according to the order of node b and node e, Represents area Ω + The damping submatrix stored in the damping matrix according to the order of node e and node b, Represents area Ω + The stiffness submatrix stored in the stiffness matrix according to the order of node b and node e, Represents area Ω + The stiffness submatrix stored in the stiffness matrix according to the order of node e and node b.

8. The structural seismic performance analysis method considering cross-fault ground motion according to claim 1 is characterized in that: The cross-fault structure includes a cross-fault tunnel and a bridge. The evaluation of the damage of the cross-fault structure according to the cross-fault seismic response analysis result and the determination of the seismic performance of the cross-fault structure include: Calculating the evaluation index value of the cross-fault tunnel or bridge according to the response of the finite element model of the cross-fault tunnel or bridge under earthquake excitation, the evaluation index includes one or more of the relative diameter deformation rate, the diameter deformation rate and the tunnel inclination angle, the displacement ductility ratio of the bridge support, and the curvature ductility ratio of the bridge pier; The evaluation index value is compared with the evaluation index threshold to determine the seismic performance of the cross-fault tunnel or bridge.

9. A structural seismic performance analysis device considering cross-fault ground motion, characterized in that: include: A data acquisition module is used to obtain earthquake magnitude data caused by specific fault activities and calculate the recurrence period of each magnitude; A fitting module, used to fit the earthquake magnitude data with the corresponding return period to establish an empirical relationship between the magnitude and the return period; A regional model determination module is used to calculate the geometric size parameters of the target fault through the empirical relationship between the magnitude and the recurrence period and the empirical relationship between the source parameters, and determine the research area and the fault rupture model across the target fault according to the geometric size parameters of the target fault; wherein the empirical relationship of the source parameters represents the relationship between the magnitude and the length and width of the fault rupture surface; The vibration simulation module is used to simulate broadband seismic motion based on the fault rupture model, digital elevation model and soil layer information model of the study area to obtain the cross-fault seismic motion field of the engineering site; An analysis module is used to establish a finite element model including a cross-fault structure and surrounding rock, convert the simulated earthquake motion into an equivalent force load on the boundary nodes of the finite element model, and then apply the equivalent force load to the boundary nodes of the finite element model to perform a cross-fault seismic response analysis on the cross-fault structure; The seismic performance determination module is used to evaluate the damage of the cross-fault structure based on the cross-fault seismic response analysis results and determine the seismic performance of the cross-fault structure.

10. A terminal comprising a memory and a processor, wherein the memory stores a computer program that can be run on the processor, characterized in that: When the processor executes the computer program, the steps of the method for analyzing the seismic performance of a structure considering cross-fault ground motion as described in any one of claims 1 to 8 are implemented.

Images