A shield tunnel anti-seismic toughness evaluation method, device and electronic equipment

Abstract

The application discloses a kind of shield tunnel seismic toughness evaluation method, device and electronic equipment, the method is established by creating grid, defining material unit and setting boundary, establishes tunnel finite element model;Selecting seismic motion-input model, according to tunnel finite element model, after obtaining response, seismic vulnerability analysis is carried out;According to the result of seismic vulnerability analysis, determine structure recovery information;According to structure recovery information, in combination with the function recovery function of structure, determine tunnel function function curve and function recovery curve;According to function recovery curve, tunnel toughness is evaluated, and the evaluation result of tunnel seismic toughness is obtained.The embodiment of the application is based on seismic toughness design, evaluates the influence of seismic motion intensity and structure factor on shield tunnel damage, and quantitatively reflects the time and cost of function recovery after earthquake, solves the technical gap of existing toughness evaluation method in the calculation of economic loss and recovery time of shield tunnel, improves the accuracy of seismic toughness evaluation, and can be widely applied in computer technology field.

Description

Technical Field

[0001] This invention relates to the field of computer technology, and in particular to a method, apparatus, and electronic device for assessing the seismic toughness of shield tunnels. Background Technology

[0002] In recent years, with the saturation of urban above-ground space, the number of newly built shield tunnels in China has been increasing. However, due to the unique environmental characteristics, the repair and reconstruction of tunnels after earthquake damage faces difficulties, leading to serious casualties, economic losses, and traffic disruptions. Especially in urban underground spaces, the economic losses during the shutdown period often exceed the repair costs. Therefore, it is urgent to study the workload and recovery time required for post-earthquake repair. Traditional seismic design focuses more on "pre-earthquake prevention" and "earthquake protection," while paying less attention to the functional recovery of tunnels after earthquakes. In recent years, design methods have been shifting towards "resilience-based" approaches, emphasizing the system's resilience under anomalous loads, including indicators such as repair costs and recovery time. This shift not only assesses the impact of earthquake intensity and structural factors on tunnel damage but also quantitatively analyzes the time and cost of post-earthquake functional recovery. However, there are still technical gaps in the evaluation of shield tunnel resilience, especially in the specific calculation methods for resilience indicators.

[0003] Current research on the seismic vulnerability of shield tunnels mainly relies on damage indicators related to use and internal forces, but does not consider the impact of joint opening on the lining structure. Traditional seismic design focuses more on "pre-earthquake prevention" and "earthquake protection," while paying less attention to the post-earthquake functional recovery of tunnels; in addition, there is a lack of accurate toughness assessment models for shield tunnels, as well as effective methods for assessing functional loss and repair time.

[0004] In summary, the relevant technologies have the following drawbacks:

[0005] (1) Existing tunnel assessments are mostly static assessments, and risk assessments of dynamic damage to tunnels when subjected to earthquakes are rare. Although the relevant risk assessment systems are scientific and feasible in terms of theoretical research, they are difficult to apply in the face of sudden earthquake damage.

[0006] (2) There is a lack of accurate toughness assessment models for shield tunnels, as well as an effective method for assessing functional loss and repair time. Summary of the Invention

[0007] The main objective of this invention is to provide a method, device, and electronic equipment for assessing the seismic toughness of shield tunnels, which can improve the accuracy of seismic toughness evaluation.

[0008] To achieve the above objectives, one aspect of this invention proposes a method for assessing the seismic toughness of shield tunnels, comprising the following steps:

[0009] Create a mesh, define material elements, and set boundaries to build a finite element model of the tunnel;

[0010] Select the ground motion-input model, and perform seismic vulnerability analysis after obtaining the response based on the tunnel finite element model;

[0011] Based on the results of the earthquake vulnerability analysis, structural recovery information is determined; the structural recovery information includes structural functional loss, tunnel repair costs, and repair time.

[0012] Based on the structural recovery information and the structural functional recovery function, the tunnel functional function curve is determined, and then the functional recovery curve is determined.

[0013] The tunnel's seismic toughness is evaluated based on the functional recovery curve, yielding the evaluation results.

[0014] In some embodiments, creating a mesh, defining material elements, and setting boundaries to establish a tunnel finite element model includes the following steps:

[0015] Soil simulation was performed using saturated porous media coupled with solid-fluid elements. Liquefied and non-liquefied soils were simulated using elastoplastic models that were related to and unrelated to confining pressure, respectively.

[0016] The bolts and pressure-resistant gaskets in the lining joint are simulated using zero-length elements;

[0017] The tunnel lining segments are simulated by combining nonlinear beam elements based on fiber cross sections and quadrilateral solid elements;

[0018] A segmented thin-layer interface element with weakened shear parameters is used to simulate the tunnel-soil contact surface.

[0019] Fix the vertical degrees of freedom of all nodes at the bottom of the finite element model, and completely fix the nodes on both sides of the bottom.

[0020] The model is subjected to the same displacement constraint and extremely heavy soil columns on the left and right boundaries to simulate the free field boundary.

[0021] The groundwater level was set 1m below the top surface of the model to establish a finite element model.

[0022] In some embodiments, the selection of the ground motion-input model and the seismic vulnerability analysis after obtaining the response based on the tunnel finite element model include the following steps: constructing the earthquake using the IDA incremental dynamic analysis method to obtain 200 ground motions;

[0023] The opening amount on the outer side of the joint is used as a structural damage index to characterize the seismic response and degree of damage of the structure.

[0024] Based on the different opening amounts on the outer side of the joint, the damage state is divided into slight damage, moderate damage and severe damage, with corresponding damage state limits of 2mm, 4mm and 6mm.

[0025] Through fuzzy comprehensive evaluation, the optimal strength index for the current site is determined to be the velocity spectrum intensity.

[0026] The seismic response was obtained based on the finite element model of the tunnel, and seismic vulnerability analysis was performed.

[0027] In some embodiments, the expression for the seismic vulnerability analysis is:

[0028]

[0029] Where P[EDP≥LS|IM] represents the probability that the seismic demand parameter exceeds the limit state under a given intensity index; LS represents the limit state of the seismic demand parameter under seismic action. Represents the standard normal distribution function; β EDP|IM S represents the logarithmic standard deviation of the seismic demand parameter under a given strength index; D This represents the median value of the seismic demand parameters under a given strength index.

[0030] In some embodiments, determining structural recovery information based on the results of the seismic vulnerability analysis includes the following steps:

[0031] The loss ratio was determined using Monte Carlo simulation.

[0032] After determining the functional loss ratio of the structure, the functional loss of the structure under different failure states is obtained by combining the vulnerability curve with the velocity spectrum intensity as the strength index.

[0033] The repair costs of shield tunnels are estimated based on the tunnel construction costs, and the required repair costs for tunnels under different velocity spectrum intensities are obtained.

[0034] The total maintenance time is calculated by summing the products of the probabilities of various failure states and the corresponding repair times for each failure state.

[0035] In some embodiments, the formula for calculating the functional loss of the structure is:

[0036]

[0037] Among them, S loss R represents the direct functional loss of the tunnel under a given specific IM condition; j P[DS] represents the repair ratio, the percentage required to repair the j-th damage state; j [IM] represents the probability that the tunnel is in a state of failure j and an earthquake intensity level IM;

[0038] The formula for calculating the total repair time is:

[0039]

[0040] Among them, T RE T represents the required tunnel repair time under a given specific IM; REj P[DS] represents the repair time required for the tunnel structure under different damage conditions; j [IM] represents the probability that the tunnel is in a damaged state j and the earthquake intensity level IM.

[0041] In some embodiments, determining the tunnel function function curve based on the structural recovery information and the structural function recovery function, and then determining the function recovery curve, includes the following steps:

[0042] Based on the importance of the structure and the repair resources, the corresponding functional recovery function is selected. The functional recovery function includes: linear function, trigonometric function and exponential function.

[0043] Based on the aforementioned function recovery function, and considering the function loss and repair time, a tunnel function function curve is plotted.

[0044] After determining the functional loss and repair time, the functional recovery curve under the given intensity index is plotted based on the functional recovery function.

[0045] In some embodiments, the process of evaluating tunnel toughness based on the functional recovery curve to obtain the evaluation result of tunnel seismic toughness includes the following steps:

[0046] Through formula The seismic toughness of shield tunnels under mainshock and aftershock is evaluated, where R represents the seismic toughness of shield tunnels under mainshock and aftershock; T RE Indicates the repair time required for a structure to return to full functionality after an earthquake; T 0E The time at which the structure fails is represented; Q(t) represents the tunnel function curve.

[0047] Determine the seismic toughness index of the tunnel when multiple earthquakes occur;

[0048] Based on the calculated seismic toughness index of the shield tunnel, the seismic toughness of the tunnel is quantitatively evaluated, and the seismic toughness of the tunnel is divided into three levels through the seismic toughness index.

[0049] Based on the obtained response results, the corresponding resilience assessment level is determined.

[0050] Another aspect of this invention provides a device for assessing the seismic toughness of a shield tunnel, comprising:

[0051] The first module is used to create the mesh, define material elements and set boundaries to build the tunnel finite element model.

[0052] The second module is used to select the seismic motion input model and perform seismic vulnerability analysis based on the tunnel finite element model after obtaining the response.

[0053] The third module is used to determine structural recovery information based on the results of the seismic vulnerability analysis; the structural recovery information includes structural functional loss, tunnel repair costs, and repair time.

[0054] The fourth module is used to determine the tunnel function function curve based on the structural recovery information and the structural function recovery function, and then to determine the function recovery curve.

[0055] The fifth module is used to evaluate the tunnel's toughness based on the functional recovery curve, and to obtain the evaluation results of the tunnel's seismic toughness.

[0056] To achieve the above objectives, another aspect of the present invention provides an electronic device, which includes a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the method described above.

[0057] To achieve the above objectives, another aspect of the present invention provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the methods described above.

[0058] This invention also discloses a computer program product or computer program, which includes computer instructions stored in a computer-readable storage medium. A processor of a computer device can read the computer instructions from the computer-readable storage medium and execute the computer instructions, causing the computer device to perform the aforementioned method.

[0059] The embodiments of this invention include at least the following beneficial effects: This invention provides a method, device, and electronic equipment for assessing the seismic toughness of shield tunnels. This scheme establishes a tunnel finite element model by creating a mesh, defining material elements, and setting boundaries; selecting a seismic motion-input model; and performing seismic vulnerability analysis after acquiring the response based on the tunnel finite element model; determining structural recovery information based on the results of the seismic vulnerability analysis; the structural recovery information includes structural functional loss, tunnel repair costs, and repair time; determining the tunnel functional function curve based on the structural recovery information and the structural functional recovery function, and thus determining the functional recovery curve; and assessing the tunnel toughness based on the functional recovery curve to obtain the evaluation result of the tunnel's seismic toughness. The embodiments of this invention, based on seismic toughness design, assess the impact of seismic motion intensity and structural factors on shield tunnel damage, and quantitatively reflect the time and cost of post-earthquake functional recovery. This addresses the technical gap in existing toughness evaluation methods regarding the calculation of economic losses and recovery time for shield tunnels, and improves the accuracy of seismic toughness evaluation. Attached Figure Description

[0060] Figure 1 This is a schematic diagram of an implementation environment provided by an embodiment of the present invention;

[0061] Figure 2 This is a flowchart of the overall steps provided in the embodiments of the present invention;

[0062] Figure 3 This is a flowchart illustrating the specific implementation steps provided in the embodiments of the present invention;

[0063] Figure 4 This is a schematic diagram of a tunnel finite element model provided in an embodiment of the present invention;

[0064] Figure 5 This is a schematic diagram of the recovery model provided in an embodiment of the present invention;

[0065] Figure 6 This is a schematic diagram of the seismic toughness provided in an embodiment of the present invention;

[0066] Figure 7 This is a schematic diagram of the hardware structure of the electronic device provided in an embodiment of the present invention. Detailed Implementation

[0067] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In the following description, when referring to the accompanying drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the embodiments of this invention; they are merely examples of apparatuses and methods consistent with some aspects of the embodiments of this invention as detailed in the appended claims.

[0068] It is understood that the terms “first,” “second,” etc., used in this invention may be used herein to describe various concepts, but unless specifically stated otherwise, these concepts are not limited by these terms. These terms are used only to distinguish one concept from another. For example, first information may also be referred to as second information without departing from the scope of embodiments of the invention, and similarly, second information may also be referred to as first information. Depending on the context, the words “if,” “when,” or “in response to determination” as used herein may be interpreted as “when…” or “when…” or “in response to determination.”

[0069] The terms “at least one,” “multiple,” “each,” “any,” etc., used in this invention, “at least one” includes one, two, or more than two; “multiple” includes two or more than two; “each” refers to each of the corresponding multiple; and “any” refers to any one of the multiple.

[0070] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains. The terminology used herein is for the purpose of describing embodiments of the invention only and is not intended to limit the invention.

[0071] Before providing a detailed description of the embodiments of the present invention, some related technologies involved in the embodiments of the present invention will be described first, as follows:

[0072] Seismic toughness refers to the ability to resist damage when subjected to earthquakes and the ability to quickly restore function after an earthquake.

[0073] The loss ratio is defined as the ratio of the function lost by the structure to the initial function.

[0074] The total recovery time of a tunnel structure is defined as the time between the occurrence of earthquake damage and the completion of building repairs.

[0075] The functional recovery function refers to the repair function that allows a structure to gradually recover to its full function after being damaged by an earthquake.

[0076] The shield tunnel seismic toughness assessment method, apparatus, and electronic equipment provided in this invention relate to the field of computer technology. The shield tunnel seismic toughness assessment method provided in this invention can be applied to a terminal, a server, or software running on a terminal or server. In some embodiments, the terminal can be a smartphone, tablet, laptop, desktop computer, smart speaker, smartwatch, or vehicle-mounted terminal, but is not limited to these. The server can be configured as an independent physical server, a server cluster or distributed system composed of multiple physical servers, or a cloud server providing basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms. The server can also be a node server in a blockchain network. The software can be an application implementing the shield tunnel seismic toughness assessment method, but is not limited to the above forms.

[0077] This invention can be used in a wide variety of general-purpose or special-purpose computer system environments or configurations. Examples include: personal computers, server computers, handheld or portable devices, tablet devices, multiprocessor systems, microprocessor-based systems, set-top boxes, programmable consumer electronics, network PCs, minicomputers, mainframe computers, and distributed computing environments including any of the above systems or devices. This invention can be described in the general context of computer-executable instructions, such as program modules, that are executed by a computer. Generally, program modules include routines, programs, objects, components, data structures, etc., that perform specific tasks or implement specific abstract data types. This invention can also be practiced in distributed computing environments where tasks are performed by remote processing devices connected via a communication network. In distributed computing environments, program modules can reside in local and remote computer storage media, including storage devices.

[0078] like Figure 1 The diagram shown is a schematic representation of an implementation environment provided by an embodiment of the present invention. (Refer to...) Figure 1 The implementation environment includes at least one terminal 102 and a server 101. The terminal 102 and the server 101 can be connected via a network, either wirelessly or via a wired connection, to complete data transmission and exchange.

[0079] Server 101 can be a standalone physical server, a server cluster or distributed system composed of multiple physical servers, or a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, CDN (Content Delivery Network), and big data and artificial intelligence platforms.

[0080] Additionally, server 101 can also be a node server in a blockchain network. Blockchain is a novel application model of computer technologies such as distributed data storage, peer-to-peer transmission, consensus mechanisms, and encryption algorithms.

[0081] Terminal 102 can be a smartphone, tablet, laptop, desktop computer, smart speaker, smartwatch, etc. It can also be a vehicle-mounted terminal of the various device types described above, but is not limited to these. Terminal 102 and server 101 can be directly or indirectly connected via wired or wireless communication, and this embodiment of the invention does not impose any limitations.

[0082] Exemplary based on Figure 1 The implementation environment shown in this embodiment of the invention provides a method for assessing the seismic toughness of a shield tunnel. The following description uses the application of this method to server 101 as an example. It can be understood that this method can also be applied to terminal 102.

[0083] Reference Figure 2 , Figure 2 This is a flowchart illustrating a method for assessing the seismic toughness of a shield tunnel applied to a server, provided in an embodiment of the present invention. The execution subject of this method can be any of the aforementioned computer devices (including a server or a terminal). (Refer to...) Figure 2 The method may include the following steps:

[0084] Create a mesh, define material elements, and set boundaries to build a finite element model of the tunnel;

[0085] Select the ground motion-input model, and perform seismic vulnerability analysis after obtaining the response based on the tunnel finite element model;

[0086] Based on the results of the earthquake vulnerability analysis, structural recovery information is determined; the structural recovery information includes structural functional loss, tunnel repair costs, and repair time.

[0087] Based on the structural recovery information and the structural functional recovery function, the tunnel functional function curve is determined, and then the functional recovery curve is determined.

[0088] The tunnel's seismic toughness is evaluated based on the functional recovery curve, yielding the evaluation results.

[0089] In some embodiments, creating a mesh, defining material elements, and setting boundaries to establish a tunnel finite element model includes the following steps:

[0090] Soil simulation was performed using saturated porous media coupled with solid-fluid elements. Liquefied and non-liquefied soils were simulated using elastoplastic models that were related to and unrelated to confining pressure, respectively.

[0091] The bolts and pressure-resistant gaskets in the lining joint are simulated using zero-length elements;

[0092] The tunnel lining segments are simulated by combining nonlinear beam elements based on fiber cross sections and quadrilateral solid elements;

[0093] A segmented thin-layer interface element with weakened shear parameters is used to simulate the tunnel-soil contact surface.

[0094] Fix the vertical degrees of freedom of all nodes at the bottom of the finite element model, and completely fix the nodes on both sides of the bottom.

[0095] The model is subjected to the same displacement constraint and extremely heavy soil columns on the left and right boundaries to simulate the free field boundary.

[0096] The groundwater level was set 1m below the top surface of the model to establish a finite element model.

[0097] In some embodiments, the selection of the ground motion-input model and the seismic vulnerability analysis after obtaining the response based on the tunnel finite element model include the following steps: constructing the earthquake using the IDA incremental dynamic analysis method to obtain 200 ground motions;

[0098] The opening amount on the outer side of the joint is used as a structural damage index to characterize the seismic response and degree of damage of the structure.

[0099] Based on the different opening amounts on the outer side of the joint, the damage state is divided into slight damage, moderate damage and severe damage, with corresponding damage state limits of 2mm, 4mm and 6mm.

[0100] Through fuzzy comprehensive evaluation, the optimal strength index for the current site is determined to be the velocity spectrum intensity.

[0101] The seismic response was obtained based on the finite element model of the tunnel, and seismic vulnerability analysis was performed.

[0102] In some embodiments, the expression for the seismic vulnerability analysis is:

[0103]

[0104] Where P[DEP≥LS|IM] represents the probability that the seismic demand parameter exceeds the limit state under a given intensity index; LS represents the limit state of the seismic demand parameter under seismic action. Represents the standard normal distribution function; β EDP|IM S represents the logarithmic standard deviation of the seismic demand parameter under a given strength index; D This represents the median value of the seismic demand parameters under a given strength index.

[0105] In some embodiments, determining structural recovery information based on the results of the seismic vulnerability analysis includes the following steps:

[0106] The loss ratio was determined using Monte Carlo simulation.

[0107] After determining the functional loss ratio of the structure, the functional loss of the structure under different failure states is obtained by combining the vulnerability curve with the velocity spectrum intensity as the strength index.

[0108] The repair costs of shield tunnels are estimated based on the tunnel construction costs, and the required repair costs for tunnels under different velocity spectrum intensities are obtained.

[0109] The total maintenance time is calculated by summing the products of the probabilities of various failure states and the corresponding repair times for each failure state.

[0110] In some embodiments, the formula for calculating the functional loss of the structure is:

[0111]

[0112] Among them, S loss R represents the direct functional loss of the tunnel under a given specific IM condition; j P[DS] represents the repair ratio, the percentage required to repair the j-th damage state; j [IM] represents the probability that the tunnel is in a state of failure j and an earthquake intensity level IM;

[0113] The formula for calculating the total repair time is:

[0114]

[0115] Among them, T RE T represents the required tunnel repair time under a given specific IM; REj P[DS] represents the repair time required for the tunnel structure under different damage conditions; j [IM] represents the probability that the tunnel is in a damaged state j and the earthquake intensity level IM.

[0116] In some embodiments, determining the tunnel function function curve based on the structural recovery information and the structural function recovery function, and then determining the function recovery curve, includes the following steps:

[0117] Based on the importance of the structure and the repair resources, the corresponding functional recovery function is selected. The functional recovery function includes: linear function, trigonometric function and exponential function.

[0118] Based on the aforementioned function recovery function, and considering the function loss and repair time, a tunnel function function curve is plotted.

[0119] After determining the functional loss and repair time, the functional recovery curve under the given intensity index is plotted based on the functional recovery function.

[0120] In some embodiments, the process of evaluating tunnel toughness based on the functional recovery curve to obtain the evaluation result of tunnel seismic toughness includes the following steps:

[0121] Through formula The seismic toughness of the shield tunnel under the action of the main shock and aftershock is evaluated, where E represents the seismic toughness of the shield tunnel under the action of the main shock and aftershock; T RE Indicates the repair time required for a structure to return to full functionality after an earthquake; T 0E The time at which the structure fails is represented; Q(t) represents the tunnel function curve.

[0122] Determine the seismic toughness index of the tunnel when multiple earthquakes occur;

[0123] Based on the calculated seismic toughness index of the shield tunnel, the seismic toughness of the tunnel is quantitatively evaluated, and the seismic toughness of the tunnel is divided into three levels through the seismic toughness index.

[0124] Based on the obtained response results, the corresponding resilience assessment level is determined.

[0125] The following describes the specific implementation process of the embodiments of the present invention in detail using a specific application scenario as an example:

[0126] This invention provides a specific calculation method and steps for toughness indicators such as economic loss and recovery time of shield tunnels. Based on seismic toughness design, this method assesses the impact of seismic intensity and structural factors on shield tunnel damage, and quantitatively reflects the time and cost of post-earthquake functional recovery. This innovation addresses the technical gap in existing toughness evaluation methods for calculating economic loss and recovery time of shield tunnels, providing an effective method for seismic toughness evaluation.

[0127] refer to Figure 3 The seismic toughness assessment process for shield tunnels according to embodiments of the present invention includes the following steps:

[0128] Step 1: Based on existing technology, create a mesh, define material elements and set boundaries to establish a finite element model of the tunnel.

[0129] The soil was simulated using water-soil coupled elements (i.e., saturated porous media coupled solid-fluid elements) in OpenSees. Liquefied and non-liquefied soil materials were simulated using elasto-plastic models related to and independent of confining pressure, respectively. Bolts and pressure-resistant gaskets in the lining joints were simulated using zero-length elements. Nonlinear beam elements based on fiber sections and quadrilateral solid elements were combined to simulate tunnel segments. Shear parameter-weakened segmented thin-layer interface elements were used to simulate the tunnel-soil contact surface. All nodes at the bottom of the finite element model had their vertical degrees of freedom fixed, while nodes on the bottom two sides were completely fixed. Equal displacement constraints and extremely heavy soil columns were applied to the left and right boundaries of the model to simulate free-field boundaries. The groundwater level was located 1m below the top surface of the model. A finite element model was established; the specific numerical model is as follows: Figure 4 .

[0130] Step 2: Select ground motion - input model, obtain the response and then perform earthquake vulnerability analysis.

[0131] The IDA incremental dynamic analysis method was used to construct 200 ground motions. The joint lateral opening δ was used as the structural damage index (DM) to characterize the seismic response and damage degree of the structure. Based on different joint lateral opening values, the damage state was divided into minor, moderate, and severe damage, with corresponding damage state thresholds of 2 mm, 4 mm, and 6 mm, respectively. Ou Entong used fuzzy comprehensive evaluation to determine that the optimal strength index for this site was VSI (velocity spectrum intensity).

[0132] Under the condition of seismic intensity (IM=x), the specific expression is shown in formula (1).

[0133]

[0134] Where P[EDP≥LS|IM] represents the probability that the seismic demand parameter exceeds the limit state under a given intensity index; LS (limit states) represents the limit state of the seismic demand parameter under seismic action. Represents the standard normal distribution function; β EDP|IM S represents the logarithmic standard deviation of the seismic demand parameter under a given strength index; D This represents the median value of the seismic demand parameters under a given strength index.

[0135] Step 3: Combine the fragility curve to determine the structural functional loss, tunnel repair cost, and repair time.

[0136] Monte Carlo simulation was used to determine the loss ratio.

[0137] Table 1 shows the ratio of functional losses of tunnel lining under different conditions in the FEMA-P58 specification issued by the FEMA (Federal Emergency Management Agency) in 2012.

[0138] Table 1

[0139]

[0140] After determining the functional loss ratio of the structure, and combining the vulnerability curve with VSI as the strength index, the functional loss of the structure under different failure states can be obtained according to formula (2).

[0141]

[0142] Among them, S loss R represents the direct functional loss of the tunnel under a given specific IM condition; j P[DS] represents the repair ratio, the percentage required to repair the j-th damage state; j [IM] represents the probability that the tunnel is in a damaged state j and the earthquake intensity level IM.

[0143] The repair costs of shield tunnels were estimated based on domestic tunnel construction costs, assuming a repair cost of 1 million yuan per meter of tunnel. According to S... loss The repair cost required for the tunnel under different VSI conditions can be obtained by multiplying the repair cost and the repair length (in this embodiment, the length is set to 1m).

[0144] Based on the findings of expert surveys in FEMA-P58, issued by the FEMA (Federal Emergency Management Agency) in 2012, the repair time for tunnels is divided into minor, moderate, and severe damage, as shown in Table 2.

[0145] Table 2 Repair time under different damage states

[0146]

[0147] Total repair time T RE The probability O[DS] of various destructive states can be used to determine the destructive state. j |IM] and the repair time T corresponding to the damaged state. REj The calculation is performed by summing the products of the products, as shown in equation (3).

[0148]

[0149] T RE T represents the required tunnel repair time under a given specific IM;REj P[DS] represents the repair time required for the tunnel structure under different damage conditions; j [IM] represents the probability that the tunnel is in a state of failure j and the earthquake intensity level (IM).

[0150] Step 4: Determine the functional recovery curve of the structure.

[0151] Three functional recovery functions can be used to easily evaluate different situations. Depending on the importance of the structure and the resources available for repair, different functional recovery functions can be selected. These mainly include the following three categories: linear functions, trigonometric functions, and exponential functions.

[0152] Using the recovery function described above, combined with the functional loss and repair time, the tunnel function curve Q(t) can be plotted, as shown in Equation 4:

[0153] Q(t) = 1 - S loss ×f rec (t) (4)

[0154] After determining the functional loss and recovery time, a functional recovery curve under a given intensity index (IM) can be plotted based on the functional recovery function, such as... Figure 5 As shown.

[0155] Step 5: Conduct a tunnel toughness assessment.

[0156] The seismic toughness of shield tunnels under main shock and aftershock is evaluated, and the formula is shown in equation (5):

[0157]

[0158] The seismic toughness index R can be used Figure 6 As shown:

[0159] When multiple earthquakes occur, the seismic toughness of the tunnel is defined as shown in formula (6):

[0160]

[0161] Among them, T RE,i t represents the repair time required for the structure to return to full functionality after the i-th earthquake, n represents the number of earthquakes that occur during the structure's entire lifespan, and t represents the total time required for the structure to recover. 0E1 …t 0En This indicates the time when the structure is destroyed from the 1st to the nth time.

[0162] The seismic toughness index R of the shield tunnel can be calculated by formula (5). The seismic toughness of the tunnel can be quantitatively evaluated by the toughness index, and the seismic toughness of the tunnel can be divided into three levels by the toughness index, as shown in Table 3. The corresponding toughness assessment level can be obtained according to the response results.

[0163] Table 3 Toughness Grade Classification

[0164] Toughness rating Toughness index R range High 0.9≤R＜1 Medium 0.6≤R＜0.9 Low R≤0.6

[0165] In summary, this invention can quantitatively assess the seismic toughness of shield tunnels. The method of this invention can determine repair costs and repair time based on the vulnerability curve established using joint opening amount, by calculating the probability of different failure states, and then select an appropriate recovery curve to assess the toughness of the shield tunnel.

[0166] Compared to existing technologies, this invention fills the gap in existing toughness evaluation methods by lacking calculation methods and procedures for seismic toughness indices of shield tunnels. It addresses the technical problem of existing structural toughness evaluation methods lacking specific calculation methods and steps for economic losses and recovery time of shield tunnels, and provides a seismic toughness evaluation method for shield tunnels.

[0167] Another aspect of this invention provides a device for assessing the seismic toughness of a shield tunnel, comprising:

[0168] The first module is used to create the mesh, define material elements and set boundaries to build the tunnel finite element model.

[0169] The second module is used to select the seismic motion input model and perform seismic vulnerability analysis based on the tunnel finite element model after obtaining the response.

[0170] The third module is used to determine structural recovery information based on the results of the seismic vulnerability analysis; the structural recovery information includes structural functional loss, tunnel repair costs, and repair time.

[0171] The fourth module is used to determine the tunnel function function curve based on the structural recovery information and the structural function recovery function, and then to determine the function recovery curve.

[0172] The fifth module is used to evaluate the tunnel's toughness based on the functional recovery curve, and to obtain the evaluation results of the tunnel's seismic toughness.

[0173] It is understood that the content of the above method embodiments is applicable to this system embodiment. The specific functions implemented in this system embodiment are the same as those in the above method embodiments, and the beneficial effects achieved are also the same as those achieved in the above method embodiments.

[0174] This invention also provides an electronic device, which includes a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the aforementioned method for assessing the seismic toughness of shield tunnels. This electronic device can be any smart terminal, including tablet computers, in-vehicle computers, etc.

[0175] It is understood that the content of the above method embodiments is applicable to this device embodiment. The specific functions implemented by this device embodiment are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0176] Please see Figure 7 , Figure 7 The hardware structure of an electronic device according to another embodiment is illustrated. The electronic device includes:

[0177] The processor 701 can be implemented using a general-purpose CPU (Central Processing Unit), microprocessor, application-specific integrated circuit (ASIC), or one or more integrated circuits, and is used to execute relevant programs to implement the technical solutions provided in the embodiments of the present invention.

[0178] The memory 702 can be implemented as a read-only memory (ROM), a static storage device, a dynamic storage device, or a random access memory (RAM). The memory 702 can store the operating system and other application programs. When the technical solutions provided in the embodiments of this specification are implemented through software or firmware, the relevant program code is stored in the memory 702 and is called and executed by the processor 701 to execute the shield tunnel seismic toughness assessment method of the embodiments of this invention.

[0179] The input / output interface 703 is used to implement information input and output;

[0180] The communication interface 704 is used to enable communication and interaction between this device and other devices. Communication can be achieved through wired means (such as USB, Ethernet cable, etc.) or wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0181] Bus 705 transmits information between various components of the device (e.g., processor 701, memory 702, input / output interface 703, and communication interface 704);

[0182] The processor 701, memory 702, input / output interface 703, and communication interface 704 are connected to each other within the device via bus 705.

[0183] This invention also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for assessing the seismic toughness of shield tunnels.

[0184] It is understood that the content of the above method embodiments is applicable to this storage medium embodiment. The specific functions implemented in this storage medium embodiment are the same as those in the above method embodiments, and the beneficial effects achieved are also the same as those achieved in the above method embodiments.

[0185] Memory, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs and non-transitory computer-executable programs. Furthermore, memory may include high-speed random access memory, and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some embodiments, memory may optionally include memory remotely located relative to the processor, and these remote memories can be connected to the processor via a network. Examples of such networks include, but are not limited to, the Internet, intranets, local area networks, mobile communication networks, and combinations thereof.

[0186] It should be noted that in various specific embodiments of the present invention, when processing data related to user identity or characteristics, such as user information, user behavior data, user historical data, and user location information, user permission or consent is obtained first. Furthermore, the collection, use, and processing of this data comply with relevant laws, regulations, and standards. In addition, when embodiments of the present invention require access to sensitive personal information of users, separate permission or consent from the user is obtained through pop-ups or redirection to a confirmation page. Only after obtaining the user's separate permission or consent is the necessary user-related data for the normal operation of the embodiments of the present invention acquired.

[0187] The embodiments described in this invention are for the purpose of more clearly illustrating the technical solutions of the embodiments of this invention, and do not constitute a limitation on the technical solutions provided by the embodiments of this invention. As those skilled in the art will know, with the evolution of technology and the emergence of new application scenarios, the technical solutions provided by the embodiments of this invention are also applicable to similar technical problems.

[0188] Those skilled in the art will understand that the technical solutions shown in the figures do not constitute a limitation on the embodiments of the present invention, and may include more or fewer steps than shown, or combine certain steps, or different steps.

[0189] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs.

[0190] Those skilled in the art will understand that all or some of the steps in the methods disclosed above, as well as the functional modules / units in the systems and devices, can be implemented as software, firmware, hardware, or suitable combinations thereof.

[0191] The terms "first," "second," "third," "fourth," etc. (if present) in the specification and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0192] It should be understood that in this invention, "at least one (item)" refers to one or more, and "more than one" refers to two or more. "And / or" describes the relationship between related objects, indicating that three relationships can exist. For example, "A and / or B" can represent three cases: only A exists, only B exists, and both A and B exist simultaneously, where A and B can be singular or plural. The character " / " generally indicates that the preceding and following related objects are in an "or" relationship. "At least one (item) of the following" or similar expressions refer to any combination of these items, including any combination of single or plural items. For example, at least one (item) of a, b, or c can represent: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", where a, b, and c can be single or multiple.

[0193] In the several embodiments provided by this invention, it should be understood that the disclosed apparatus and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative; for instance, the division of the units described above is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be through some interfaces; the indirect coupling or communication connection between devices or units may be electrical, mechanical, or other forms.

[0194] The units described above as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0195] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0196] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes multiple instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing programs, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0197] The preferred embodiments of the present invention have been described above with reference to the accompanying drawings, but this does not limit the scope of the claims of the present invention. Any modifications, equivalent substitutions, and improvements made by those skilled in the art without departing from the scope and spirit of the present invention should be within the scope of the claims of the present invention.

Claims

1. A method for assessing the seismic toughness of shield tunnels, characterized in that, Includes the following steps: Create a mesh, define material elements, and set boundaries to build a finite element model of the tunnel; Select the ground motion-input model, and perform seismic vulnerability analysis after obtaining the response based on the tunnel finite element model; Based on the results of the earthquake vulnerability analysis, structural recovery information is determined; The structural restoration information includes structural functional loss, tunnel repair costs, and repair time; Based on the structural recovery information and the structural functional recovery function, the tunnel functional function curve is determined, and then the functional recovery curve is determined. Based on the functional recovery curve, the tunnel toughness is assessed to obtain the evaluation results of the tunnel's seismic toughness. The process of creating a mesh, defining material elements, and setting boundaries to establish a finite element model of the tunnel includes the following steps: Soil simulation was performed using saturated porous media coupled with solid-fluid elements. Liquefied and non-liquefied soils were simulated using elastoplastic models that were related to and unrelated to confining pressure, respectively. The bolts and pressure-resistant gaskets in the lining joint are simulated using zero-length elements; The tunnel lining segments are simulated by combining nonlinear beam elements based on fiber cross sections and quadrilateral solid elements; A segmented thin-layer interface element with weakened shear parameters is used to simulate the tunnel-soil contact surface. Fix the vertical degrees of freedom of all nodes at the bottom of the finite element model, and completely fix the nodes on both sides of the bottom. The model is subjected to the same displacement constraint and extremely heavy soil columns on the left and right boundaries to simulate the free field boundary. The groundwater level was set 1m below the top surface of the model, and a finite element model was established. The determination of structural recovery information based on the results of the seismic vulnerability analysis includes the following steps: The loss ratio was determined using Monte Carlo simulation. After determining the functional loss ratio of the structure, the functional loss of the structure under different failure states is obtained by combining the vulnerability curve with the velocity spectrum intensity as the strength index. The repair costs of shield tunnels are estimated based on the tunnel construction costs, and the required repair costs for tunnels under different velocity spectrum intensities are obtained. The total maintenance time is calculated by summing the products of the probability of each failure state and the repair time corresponding to each failure state. The formula for calculating the functional loss of the structure is as follows: in, This represents the direct functional loss of the tunnel under a given specific IM condition; Represents the repair ratio, the percentage required to repair the j-th damage state; This indicates that the tunnel is in a state of damage. and earthquake intensity level The probability of; The formula for calculating the total repair time is: in, This represents the required repair time for a tunnel under a given specific IM; This represents the repair time required for the tunnel structure under different damage conditions. This indicates that the tunnel is in a state of damage. and earthquake intensity level The probability of.

2. The method for evaluating the seismic toughness of a shield tunnel according to claim 1, characterized in that, The selected ground motion-input model, based on the tunnel finite element model, performs seismic vulnerability analysis after obtaining the response, including the following steps: constructing the earthquake using the IDA incremental dynamic analysis method to obtain 200 ground motions; The opening amount on the outer side of the joint is used as a structural damage index to characterize the seismic response and degree of damage of the structure. Based on the different opening amounts on the outer side of the joint, the damage state is divided into slight damage, moderate damage and severe damage, with corresponding damage state limits of 2mm, 4mm and 6mm. Through fuzzy comprehensive evaluation, the optimal strength index for the current site is determined to be the velocity spectrum intensity. The seismic response was obtained based on the finite element model of the tunnel, and seismic vulnerability analysis was performed.

3. The method for evaluating the seismic toughness of a shield tunnel according to claim 2, characterized in that, The expression for the seismic vulnerability analysis is: in, This represents the probability that the seismic demand parameters exceed the limit state under a given strength index. This represents the limit state of the seismic demand parameters under seismic loading; Represents the standard normal distribution function; This represents the logarithmic standard deviation of the seismic demand parameter under a given strength index. This represents the median value of the seismic demand parameters under a given strength index.

4. The method for evaluating the seismic toughness of a shield tunnel according to claim 1, characterized in that, The process of determining the tunnel function function curve based on the structural recovery information and the structural function recovery function, and then determining the function recovery curve, includes the following steps: Based on the importance of the structure and the repair resources, the corresponding functional recovery function is selected. The functional recovery function includes: linear function, trigonometric function and exponential function. Based on the aforementioned function recovery function, and considering the function loss and repair time, a tunnel function function curve is plotted. After determining the functional loss and repair time, the functional recovery curve under the given intensity index is plotted based on the functional recovery function.

5. The method for evaluating the seismic toughness of a shield tunnel according to claim 1, characterized in that, The process of evaluating tunnel toughness based on the functional recovery curve to obtain the evaluation result of tunnel seismic toughness includes the following steps: Through formula An assessment of the seismic toughness of shield tunnels under main shock and aftershock was conducted, including... This indicates the seismic toughness of the shield tunnel under the action of the main shock and aftershocks; This indicates the repair time required for a structure to return to full functionality after earthquake damage. Indicates the time when the structure failed; Represents the tunnel function curve; Determine the seismic toughness index of the tunnel when multiple earthquakes occur; Based on the calculated seismic toughness index of the shield tunnel, the seismic toughness of the tunnel is quantitatively evaluated, and the seismic toughness of the tunnel is divided into three levels through the seismic toughness index. Based on the obtained response results, the corresponding resilience assessment level is determined.

6. An apparatus for implementing the method for assessing the seismic toughness of shield tunnels as described in any one of claims 1-5, characterized in that, include: The first module is used to create the mesh, define material elements and set boundaries to build the tunnel finite element model. The second module is used to select the seismic motion input model and perform seismic vulnerability analysis based on the tunnel finite element model after obtaining the response. The third module is used to determine structural recovery information based on the results of the earthquake vulnerability analysis. The structural restoration information includes structural functional loss, tunnel repair costs, and repair time; The fourth module is used to determine the tunnel function function curve based on the structural recovery information and the structural function recovery function, and then to determine the function recovery curve. The fifth module is used to evaluate the tunnel's toughness based on the functional recovery curve, and to obtain the evaluation results of the tunnel's seismic toughness.

7. An electronic device, characterized in that, Including the processor and memory; The memory is used to store programs; The processor executes the program to implement the method as described in any one of claims 1 to 5.

Images