Shield tunnel anti-seismic toughness evaluation method and device and electronic equipment

Abstract

The invention discloses a shield tunnel anti-seismic toughness evaluation method and device and electronic equipment. The method comprises the following steps: establishing a grid, defining a material unit and setting a boundary, and establishing a tunnel finite element model; a seismic oscillation-input model is selected, and according to the tunnel finite element model, seismic vulnerability analysis is carried out after response is obtained; determining structure recovery information according to the result of the earthquake vulnerability analysis; determining a tunnel performance function curve and a function recovery curve by combining a function recovery function of the structure according to the structure recovery information; and performing tunnel toughness evaluation according to the function recovery curve to obtain an evaluation result of tunnel anti-seismic toughness. According to the embodiment of the invention, based on the anti-seismic toughness design, the influence of the seismic oscillation intensity and structural factors on the shield tunnel damage is evaluated, the time and cost of function recovery after the earthquake are quantitatively reflected, and the technical gap of an existing toughness evaluation method in the aspects of shield tunnel economic loss and recovery time calculation is solved; the method can be widely applied to the technical field of computers.

Description

Technical Field

[0001] The present invention relates to the field of computer technology, and in particular to a shield tunnel seismic toughness assessment method, device and electronic equipment. Background Art

[0002] In recent years, with the saturation of urban ground space, the number of newly built shield tunnels in China has continued to increase. However, due to the special environment, the repair and reconstruction of tunnels after earthquake damage are difficult, resulting in serious casualties, economic losses and traffic interruptions. Especially in urban underground spaces, the economic losses during the suspension 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 on "pre-earthquake prevention" and "earthquake protection", but pays less attention to the post-earthquake function recovery of tunnels. In recent years, the design method is shifting to "toughness-based", emphasizing the system's ability to recover under abnormal conditions, including indicators such as repair cost and recovery time. This shift not only evaluates the impact of earthquake intensity and structural factors on tunnel damage, but also quantitatively analyzes the time and cost of post-earthquake function recovery. However, there is still a technical gap in the evaluation of shield tunnel toughness, especially in the specific calculation method of toughness 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", but pays less attention to the post-earthquake function recovery of tunnels; in addition, there is a lack of accurate toughness assessment models for shield tunnels, and there is also a lack of effective methods for assessing functional loss and repair time.

[0004] In summary, the related technology has the following disadvantages:

[0005] (1) Existing tunnel assessments are mostly static assessments. Risk assessments of dynamic damage to tunnels when they are hit by earthquakes are relatively rare. In addition, although the relevant risk assessment system is scientific and feasible in theoretical research, it is difficult to apply it to sudden earthquake damage.

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

[0007] The main purpose of the embodiments of the present invention is to provide a shield tunnel seismic toughness assessment method, device and electronic equipment, which can improve the accuracy of seismic toughness evaluation.

[0008] To achieve the above object, an embodiment of the present invention provides a method for evaluating the seismic toughness of a shield tunnel, comprising the following steps:

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

[0010] Selecting a seismic motion-input model, and performing seismic vulnerability analysis after obtaining responses based on the tunnel finite element model;

[0011] Determining structural restoration information according to the results of the seismic vulnerability analysis; the structural restoration information includes structural function loss, tunnel repair costs and repair time;

[0012] According to the structural recovery information, combined with the functional recovery function of the structure, a tunnel functional function curve is determined, and then a functional recovery curve is determined;

[0013] The tunnel toughness is evaluated according to the functional recovery curve to obtain an evaluation result of the tunnel seismic toughness.

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

[0015] The saturated porous medium coupled solid-fluid unit is used to simulate the soil, and the elastoplastic model related to or independent of the confining pressure is used to simulate the liquefied soil and non-liquefied soil respectively.

[0016] The bolts and compression pads in the lining joints are simulated by zero-length elements;

[0017] Combining fiber-cross-section-based nonlinear beam elements and quadrilateral solid elements to simulate segments;

[0018] The segmented thin-layer interface element with shear parameter weakening is used to simulate the tunnel-soil interface.

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

[0020] Displacement equal constraints and extremely heavy soil columns are applied to the left and right boundaries of the model to simulate the free field boundary;

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

[0022] In some embodiments, the seismic fragility analysis is performed after obtaining the response based on the tunnel finite element model by selecting the seismic motion-input model, including the following steps: constructing the seismic motion using the IDA incremental dynamic analysis method to obtain 200 seismic motions;

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

[0024] According to the different opening amounts of the joint outside, the damage state is divided into slight damage, moderate damage and severe damage, and the corresponding damage state limits are 2mm, 4mm and 6mm;

[0025] Through fuzzy comprehensive evaluation, it is determined that the optimal intensity index of the current site is the velocity spectrum intensity;

[0026] The seismic response is obtained according to the tunnel finite element model, and seismic vulnerability analysis is performed.

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

[0028]

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

[0030] In some embodiments, determining structural restoration information according to the results of the seismic vulnerability analysis comprises the following steps:

[0031] Monte Carlo simulation method is used to determine the loss ratio;

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

[0033] The repair cost of shield tunnel is estimated according to the tunnel construction cost, and the repair cost of tunnel under different velocity spectrum intensities is obtained;

[0034] The total maintenance time is calculated by summing the products of the probabilities of various damage states and the repair times corresponding to the damage states.

[0035] In some embodiments, the functional loss of the structure is calculated as:

[0036]

[0037] Among them, S loss represents the direct functional loss of the tunnel under a given specific IM situation; R j represents the repair ratio, the percentage required to repair the jth damage state; P[DS j |IM] represents the probability of the tunnel being located in damage state j and earthquake intensity level IM;

[0038] The total maintenance time is calculated as:

[0039]

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

[0041] In some embodiments, the step of determining a tunnel function function curve based on the structural recovery information and combining the functional recovery function of the structure, and then determining a functional recovery curve, comprises the following steps:

[0042] According to the importance of the structure and the repair resources, a corresponding function recovery function is selected, wherein the function recovery function includes: a linear function, a trigonometric function and an exponential function;

[0043] According to the function recovery function, combined with the function loss and the repair time, a tunnel function function curve is drawn;

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

[0045] In some embodiments, performing tunnel toughness assessment according to the functional recovery curve to obtain an evaluation result of tunnel seismic toughness comprises the following steps:

[0046] By formula The seismic toughness of the shield tunnel under the main aftershocks was evaluated, where R represents the seismic toughness of the shield tunnel under the main aftershocks; T RE It indicates the time required for the structure to recover to its full function after the earthquake damage. 0E It indicates the time when the structure is damaged; Q(t) indicates the tunnel function curve;

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

[0048] According to the calculated shield tunnel seismic toughness index, the tunnel seismic toughness is quantitatively evaluated and divided into three levels according to the seismic toughness index.

[0049] Determine the corresponding toughness assessment level based on the response results obtained.

[0050] Another aspect of the embodiment of the present invention further provides a shield tunnel seismic toughness assessment device, comprising:

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

[0052] The second module is used to select a seismic motion-input model and perform seismic vulnerability analysis after obtaining the response according to the tunnel finite element model;

[0053] A third module is used to determine structural restoration information according to the results of the seismic vulnerability analysis; the structural restoration information includes structural function loss, tunnel repair costs and repair time;

[0054] The fourth module is used to determine the tunnel function function curve according to the structural recovery information and the functional recovery function of the structure, and then determine the functional recovery curve;

[0055] The fifth module is used to evaluate the toughness of the tunnel according to the functional recovery curve to obtain an evaluation result of the seismic toughness of the tunnel.

[0056] To achieve the above objective, another aspect of an embodiment of the present invention provides an electronic device, the electronic device comprising a memory and a processor, the memory storing a computer program, and the processor implementing the above-mentioned method when executing the computer program.

[0057] To achieve the above objective, another aspect of an embodiment of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned method is implemented.

[0058] The embodiment of the present invention also discloses a computer program product or a computer program, which includes a computer instruction stored in a computer-readable storage medium. A processor of a computer device can read the computer instruction from the computer-readable storage medium, and the processor executes the computer instruction, so that the computer device executes the above method.

[0059] The embodiments of the present invention include at least the following beneficial effects: the present invention provides a shield tunnel seismic toughness assessment method, device and electronic equipment, the scheme establishes a tunnel finite element model by creating a grid, defining material units and setting boundaries; selects a seismic motion-input model, and performs seismic vulnerability analysis according to the tunnel finite element model after obtaining the response; determines structural recovery information according to the results of the seismic vulnerability analysis; the structural recovery information includes structural function loss, tunnel repair costs and repair time; determines a tunnel function function curve based on the structural recovery information and the functional recovery function of the structure, and then determines a functional recovery curve; performs tunnel toughness assessment according to the functional recovery curve to obtain an evaluation result of the seismic toughness of the tunnel. The embodiments of the present invention evaluate the influence of seismic motion intensity and structural factors on shield tunnel damage based on seismic toughness design, and quantitatively reflects the time and cost of post-earthquake functional recovery, which solves the technical gap of existing toughness assessment methods in calculating economic losses and recovery time of shield tunnels, and improves the accuracy of seismic toughness assessment. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0061] Figure 2 is a flow chart of the overall steps provided by an embodiment of the present invention;

[0062] Figure 3 It is a flowchart of specific implementation steps provided by an embodiment of the present invention;

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

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

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

[0066] Figure 7 It is a schematic diagram of the hardware structure of an electronic device provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0067] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention. When the following description refers to the 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 the present invention. They are only examples of devices and methods consistent with some aspects of the embodiments of the present invention as detailed in the attached claims.

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

[0069] The terms "at least one", "multiple", "each", "any", etc. used in the present invention, at least one includes one, two or more, multiple includes two or more, 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 those commonly understood by those skilled in the art to which the present invention belongs. The terms used herein are only for the purpose of describing the embodiments of the present invention and are not intended to limit the present invention.

[0071] Before describing the embodiments of the present invention in detail, some related technologies involved in the embodiments of the present invention are first described as follows:

[0072] Seismic resilience refers to the ability to resist damage when hit by an earthquake and the ability to quickly recover its functions 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 restoration time of the tunnel structure is defined as the time from the occurrence of earthquake damage to the completion of building repair.

[0075] Functional recovery function refers to the repair function of a structure that gradually recovers to its full function after being damaged by an earthquake.

[0076] The shield tunnel seismic resilience assessment method, device and electronic device provided in the embodiment of the present invention relate to the field of computer technology. The shield tunnel seismic resilience assessment method provided in the embodiment of the present invention can be applied to a terminal, can be applied to a server, and can also be software running in a terminal or a server. In some embodiments, the terminal can be a smart phone, a tablet computer, a laptop computer, a desktop computer, a smart speaker, a smart watch, and a vehicle-mounted terminal, etc., but is not limited to this; the server side can be configured as an independent physical server, or a server cluster or distributed system composed of multiple physical servers, and can also be configured as a cloud server that provides cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms and other basic cloud computing services. The server can also be a node server in a blockchain network; the software can be an application that implements the shield tunnel seismic resilience assessment method, etc., but is not limited to the above forms.

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

[0078] like Figure 1 FIG. 1 is a schematic diagram of an implementation environment provided by an embodiment of the present invention. 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 to a network wirelessly or wired to complete data transmission and exchange.

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

[0080] In addition, the server 101 can also be a node server in the blockchain network. Blockchain is a new application model of computer technologies such as distributed data storage, peer-to-peer transmission, consensus mechanism, and encryption algorithm.

[0081] The terminal 102 may be a smart phone, a tablet computer, a laptop computer, a desktop computer, a smart speaker, a smart watch, etc. The terminal 102 may also be a vehicle-mounted terminal of various device types described above, but is not limited thereto. The terminal 102 and the server 101 may be directly or indirectly connected via wired or wireless communication, which is not limited in this embodiment of the present invention.

[0082] Based on the example Figure 1 In the implementation environment shown, an embodiment of the present invention provides a shield tunnel seismic toughness assessment method. The following is an example of the shield tunnel seismic toughness assessment method being applied in the server 101. It can be understood that the method can also be applied in the terminal 102.

[0083] Reference Figure 2 , Figure 2 The flowchart of the shield tunnel seismic toughness assessment method applied to the server provided in the embodiment of the present invention, the execution subject of the method can be any of the aforementioned computer devices (including servers or terminals). Figure 2 , the method may include the following steps:

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

[0085] Selecting a seismic motion-input model, and performing seismic vulnerability analysis after obtaining responses based on the tunnel finite element model;

[0086] Determining structural restoration information according to the results of the seismic vulnerability analysis; the structural restoration information includes structural function loss, tunnel repair costs and repair time;

[0087] According to the structural recovery information, combined with the functional recovery function of the structure, a tunnel functional function curve is determined, and then a functional recovery curve is determined;

[0088] The tunnel toughness is evaluated according to the functional recovery curve to obtain an evaluation result of the tunnel seismic toughness.

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

[0090] The saturated porous medium coupled solid-fluid unit is used to simulate the soil, and the elastoplastic model related to or independent of the confining pressure is used to simulate the liquefied soil and non-liquefied soil respectively.

[0091] The bolts and compression pads in the lining joints are simulated by zero-length elements;

[0092] Combining fiber-cross-section-based nonlinear beam elements and quadrilateral solid elements to simulate segments;

[0093] The segmented thin-layer interface element with shear parameter weakening is used to simulate the tunnel-soil interface.

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

[0095] Displacement equal constraints and extremely heavy soil columns are applied to the left and right boundaries of the model to simulate the free field boundary;

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

[0097] In some embodiments, the seismic fragility analysis is performed after obtaining the response based on the tunnel finite element model by selecting the seismic motion-input model, including the following steps: constructing the seismic motion using the IDA incremental dynamic analysis method to obtain 200 seismic motions;

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

[0099] According to the different opening amounts of the joint outside, the damage state is divided into slight damage, moderate damage and severe damage, and the corresponding damage state limits are 2mm, 4mm and 6mm;

[0100] Through fuzzy comprehensive evaluation, it is determined that the optimal intensity index of the current site is the velocity spectrum intensity;

[0101] The seismic response is obtained according to the tunnel finite element model, and seismic vulnerability analysis is performed.

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

[0103]

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

[0105] In some embodiments, determining structural restoration information according to the results of the seismic vulnerability analysis comprises the following steps:

[0106] Monte Carlo simulation method is used to determine the loss ratio;

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

[0108] The repair cost of shield tunnel is estimated according to the tunnel construction cost, and the repair cost of tunnel under different velocity spectrum intensities is obtained;

[0109] The total maintenance time is calculated by summing the products of the probabilities of various damage states and the repair times corresponding to the damage states.

[0110] In some embodiments, the functional loss of the structure is calculated as:

[0111]

[0112] Among them, S loss represents the direct functional loss of the tunnel under a given specific IM situation; R j represents the repair ratio, the percentage required to repair the jth damage state; P[DS j |IM] represents the probability of the tunnel being located in damage state j and earthquake intensity level IM;

[0113] The total maintenance time is calculated as:

[0114]

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

[0116] In some embodiments, the step of determining a tunnel function function curve based on the structural recovery information and combining the functional recovery function of the structure, and then determining a functional recovery curve, comprises the following steps:

[0117] According to the importance of the structure and the repair resources, a corresponding function recovery function is selected, wherein the function recovery function includes: a linear function, a trigonometric function and an exponential function;

[0118] According to the function recovery function, combined with the function loss and the repair time, a tunnel function function curve is drawn;

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

[0120] In some embodiments, performing tunnel toughness assessment according to the functional recovery curve to obtain an evaluation result of tunnel seismic toughness comprises the following steps:

[0121] By formula The seismic toughness of the shield tunnel under the main shock and aftershock is evaluated, where E represents the seismic toughness of the shield tunnel under the main shock and aftershock; T RE It indicates the time required for the structure to recover to its full function after the earthquake damage. 0E It indicates the time when the structure is damaged; Q(t) indicates the tunnel function curve;

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

[0123] According to the calculated shield tunnel seismic toughness index, the tunnel seismic toughness is quantitatively evaluated and divided into three levels according to the seismic toughness index.

[0124] Determine the corresponding toughness assessment level based on the response results obtained.

[0125] The specific implementation process of the embodiment of the present invention is described in detail below by taking a specific application scenario as an example:

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

[0127] refer to Figure 3 The shield tunnel seismic toughness assessment process of the embodiment of the present invention includes the following steps:

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

[0129] The water-soil coupling unit in OpenSees, that is, the saturated porous medium coupled solid-fluid unit, is used to simulate the soil. The elastoplastic models related to and independent of the confining pressure are used to simulate the liquefied soil and non-liquefied soil materials respectively. The bolts and compression pads in the lining joints are simulated in the form of zero-length units. The nonlinear beam unit based on the fiber section and the solid unit with quadrilaterals are combined to simulate the pipe segments. The segmented thin layer interface unit with weakened shear parameters is used to simulate the tunnel-soil contact surface. The degrees of freedom of all nodes at the bottom of the finite element model in the vertical direction are fixed, and the nodes on both sides of the bottom are completely fixed. The same displacement constraints and extremely heavy soil columns are applied to the left and right boundaries of the model to simulate the free field boundary. The groundwater level is located 1m below the top surface of the model. A finite element model is established, and the specific numerical model is as follows: Figure 4 .

[0130] Step 2: Select the earthquake motion-input model and perform earthquake vulnerability analysis after obtaining the response.

[0131] The IDA incremental dynamic analysis method was used to construct the earthquake, and 200 seismic motions were obtained. The joint outer opening δ will be used as the structural damage index (DM) to characterize the structural seismic response and damage degree. According to the different joint outer openings, the damage state is divided into slight damage, moderate damage and severe damage, and the corresponding damage state limits are 2mm, 4mm and 6mm. Ou Entong conducted a fuzzy comprehensive evaluation and determined that the optimal strength index of the site is VSI (velocity spectrum intensity).

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

[0133]

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

[0135] Step 3: Combined with the fragility curve, determine the structural function loss, tunnel repair cost and repair time.

[0136] The Monte Carlo simulation method is used to determine the loss ratio.

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

[0138] Table 1

[0139]

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

[0141]

[0142] Among them, S loss represents the direct functional loss of the tunnel under a given specific IM situation; R j represents the repair ratio, the percentage required to repair the jth damage state; P[DS j |IM] represents the probability of the tunnel being located in damage state j and earthquake intensity level IM.

[0143] The repair cost of shield tunnels was estimated based on the domestic tunnel construction cost. It was assumed that the repair cost per meter of tunnel was 1 million yuan. loss The repair cost required for the tunnel under different VSIs can be obtained by multiplying the repair cost and the repair length (in this embodiment, the length is assumed to be 1 m).

[0144] With the help of the expert survey results in FEMA-P58 issued by FEMA (Federal Emergency Management Agency) in 2012, the repair time of the tunnel 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 maintenance time T RE The probability of various damage states O[DS j |IM] and the repair time T corresponding to the damaged state REj It is calculated by taking the sum of the products of , as shown in formula (3).

[0148]

[0149] T RE represents the required repair time of the tunnel under a given specific IM; TREj represents the required repair time of the tunnel structure under different damage conditions; P[DS j |IM] represents the probability of the tunnel being located in the damage state j and earthquake intensity level (IM).

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

[0151] There are three functional recovery functions that can be used to simply evaluate different situations. Different functional recovery functions can be selected according to the importance of the structure and the repair resources. They mainly include the following three categories: linear function, trigonometric function and exponential function.

[0152] Through the above recovery function, combined with the function loss and repair time, the tunnel function function curve Q(t) can be drawn, as shown in Formula 4:

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

[0154] After determining the functional loss and repair time, the functional recovery curve under a given intensity index (IM) can be drawn according to the functional recovery function, such as Figure 5 shown.

[0155] Step 5: Conduct tunnel toughness assessment.

[0156] The seismic toughness of the shield tunnel under the action of the main aftershock is evaluated, and the formula is shown in formula (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 represents the repair time required for the structure to recover to its full function after the i-th earthquake damage, n represents the number of earthquakes that occur in the structure during its life cycle, and t 0E1 …t 0En Indicates the time from the 1st to the nth time the structure is damaged.

[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 evaluation level can be obtained according to the response results.

[0163] Table 3 Toughness grade classification

[0164] Toughness grade Toughness index R range high 0.9≤R＜1 middle 0.6≤R＜0.9 Low R≤0.6

[0165] In summary, the present invention can quantitatively evaluate the seismic toughness of shield tunnels. The method of the present invention can determine the repair cost and repair time based on the probability of different damage states based on the fragility curve established using the joint opening, and then select a suitable recovery curve to evaluate the toughness of the shield tunnel.

[0166] Compared with the prior art, the present invention fills the gap in the existing toughness evaluation method, which lacks a calculation method and process for the seismic toughness index of the shield tunnel. It solves the technical problem that the existing structural toughness evaluation method lacks a specific calculation method and steps for the economic loss and recovery time of the shield tunnel, and provides a shield tunnel seismic toughness evaluation method.

[0167] Another aspect of the embodiment of the present invention further provides a shield tunnel seismic toughness assessment device, comprising:

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

[0169] The second module is used to select a seismic motion-input model and perform seismic vulnerability analysis after obtaining the response according to the tunnel finite element model;

[0170] A third module is used to determine structural restoration information according to the results of the seismic vulnerability analysis; the structural restoration information includes structural function loss, tunnel repair costs and repair time;

[0171] The fourth module is used to determine the tunnel function function curve according to the structural recovery information and the functional recovery function of the structure, and then determine the functional recovery curve;

[0172] The fifth module is used to evaluate the toughness of the tunnel according to the functional recovery curve to obtain an evaluation result of the seismic toughness of the tunnel.

[0173] It can be understood that the contents of the above method embodiments are all applicable to the present system embodiments, the functions specifically implemented by the present system embodiments 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.

[0174] The embodiment of the present invention further provides an electronic device, the electronic device comprising a memory and a processor, the memory storing a computer program, and the processor implementing the above-mentioned shield tunnel seismic toughness assessment method when executing the computer program. The electronic device can be any intelligent terminal including a tablet computer, a vehicle-mounted computer, etc.

[0175] It can be understood that the contents of the above method embodiments are all applicable to the present device embodiments, the functions specifically implemented by the present device embodiments 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] See also Figure 7 , Figure 7 The hardware structure of an electronic device of another embodiment is illustrated, and the electronic device includes:

[0177] The processor 701 may be implemented by a general-purpose CPU (Central Processing Unit), a microprocessor, an 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 in the form of a read-only memory (ROM), a static storage device, a dynamic storage device, or a random access memory (RAM). The memory 702 can store an operating system and other application programs. When the technical solution provided in the embodiment of this specification is implemented by software or firmware, the relevant program code is stored in the memory 702, and the processor 701 calls and executes the shield tunnel seismic toughness assessment method of the embodiment of the present invention;

[0179] Input / output interface 703, used to implement information input and output;

[0180] Communication interface 704, used to realize communication interaction between the device and other devices, which can be realized through wired mode (such as USB, network cable, etc.) or wireless mode (such as mobile network, WIFI, Bluetooth, etc.);

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

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

[0183] An embodiment of the present invention further provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the above-mentioned shield tunnel seismic toughness assessment method is implemented.

[0184] It can be understood that the contents of the above method embodiments are all applicable to the present storage medium embodiments, the functions specifically implemented by the present storage medium embodiments 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.

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

[0186] It should be noted that in various specific embodiments of the present invention, when it comes to the need to perform relevant processing based on data related to user identity or characteristics such as user information, user behavior data, user historical data, and user location information, the user's permission or consent will be obtained first, and the collection, use, and processing of these data will comply with relevant laws, regulations, and standards. In addition, when the embodiment of the present invention needs to obtain the user's sensitive personal information, it will obtain the user's separate permission or consent through a pop-up window or jump to a confirmation page, and after clearly obtaining the user's separate permission or consent, it will obtain the necessary user-related data for the normal operation of the embodiment of the present invention.

[0187] The embodiments described in the embodiments of the present invention are intended to more clearly illustrate the technical solutions of the embodiments of the present invention, and do not constitute a limitation on the technical solutions provided by the embodiments of the present invention. Those skilled in the art can appreciate that with the evolution of technology and the emergence of new application scenarios, the technical solutions provided by the embodiments of the present invention are also applicable to similar technical problems.

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

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

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

[0191] The terms "first", "second", "third", "fourth", etc. (if any) in the specification of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0192] It should be understood that in the present invention, "at least one (item)" refers to one or more, and "plurality" refers to two or more. "And / or" is used to describe the association relationship of associated objects, indicating that three relationships may exist. For example, "A and / or B" can represent: only A exists, only B exists, and A and B exist at the same time, where A and B can be singular or plural. The character " / " generally indicates that the previous and next associated objects are in an "or" relationship. "At least one of the following items" or similar expressions refers to any combination of these items, including any combination of single items or plural items. For example, at least one 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, c can be single or multiple.

[0193] In the several embodiments provided by the present invention, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the above units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

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

[0195] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0196] If the integrated unit is implemented in the form of 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. The computer software product is stored in a storage medium, including multiple instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods of various embodiments of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM), random access memory (RAM), disk or optical disk, etc., various media that can store programs.

[0197] The preferred embodiments of the present invention are described above with reference to the accompanying drawings, but the scope of the rights of the embodiments of the present invention is not limited thereby. Any modification, equivalent substitution and improvement made by those skilled in the art without departing from the scope and essence of the embodiments of the present invention shall be within the scope of the rights of the embodiments of the present invention.

Claims

1. A shield tunnel seismic toughness assessment method, characterized in that: The following steps are involved: Create mesh, define material elements and set boundaries to build a tunnel finite element model; Selecting a seismic motion-input model, and performing seismic vulnerability analysis after obtaining responses based on the tunnel finite element model; determining structural restoration information based on the results of the seismic vulnerability analysis; The structural restoration information includes structural function loss, tunnel restoration cost and restoration time; According to the structural recovery information, combined with the functional recovery function of the structure, a tunnel functional function curve is determined, and then a functional recovery curve is determined; The tunnel toughness is evaluated according to the functional recovery curve to obtain an evaluation result of the tunnel seismic toughness.

2. A shield tunnel seismic toughness assessment method according to claim 1, characterized in that: The step of creating a mesh, defining material units and setting boundaries to establish a tunnel finite element model includes the following steps: The saturated porous medium coupled solid-fluid unit is used to simulate the soil, and the elastoplastic model related to or independent of the confining pressure is used to simulate the liquefied soil and non-liquefied soil respectively. The bolts and compression pads in the lining joints are simulated by zero-length elements; Combining fiber-cross-section-based nonlinear beam elements and quadrilateral solid elements to simulate segments; The segmented thin-layer interface element with shear parameter weakening is used to simulate the tunnel-soil interface. Fix the degrees of freedom of all nodes at the bottom of the finite element model in the vertical direction, and completely fix the nodes on both sides of the bottom; Displacement equal constraints and extremely heavy soil columns are applied to the left and right boundaries of the model to simulate the free field boundary; The groundwater level was set 1 m below the top surface of the model and a finite element model was established.

3. A shield tunnel seismic toughness assessment method according to claim 1, characterized in that: The selecting of the earthquake motion-input model and performing seismic vulnerability analysis after obtaining the response according to the tunnel finite element model comprises the following steps: The IDA incremental dynamic analysis method was used to construct the earthquake, and 200 ground motions were obtained; The outer opening of the joint is used as a structural damage index to characterize the seismic response and damage degree of the structure. According to the different opening amounts of the joint outside, the damage state is divided into slight damage, moderate damage and severe damage, and the corresponding damage state limits are 2mm, 4mm and 6mm; Through fuzzy comprehensive evaluation, it is determined that the optimal intensity index of the current site is the velocity spectrum intensity; The seismic response is obtained according to the tunnel finite element model, and seismic vulnerability analysis is performed.

4. A shield tunnel seismic toughness assessment method according to claim 3, characterized in that: The expression of the seismic vulnerability analysis is: Among them, P[EDP≥LS|IM] represents the probability that the seismic demand parameter exceeds the limit state under a given strength index; LS represents the limit state of the seismic demand parameter under earthquake action; represents the standard normal distribution function; β EDP|IM represents the logarithmic standard deviation of the seismic demand parameter under a given strength index; S D Represents the median value of the seismic demand parameter under a given strength index.

5. A shield tunnel seismic toughness assessment method according to claim 1, characterized in that: Determining structural restoration information according to the result of the seismic vulnerability analysis comprises the following steps: Monte Carlo simulation method is used to determine the loss ratio; After determining the functional loss ratio of the structure, the functional loss of the structure under different damage states is obtained by combining the fragility curve with the velocity spectrum intensity as the strength index; The repair cost of shield tunnel is estimated according to the tunnel construction cost, and the repair cost of tunnel under different velocity spectrum intensities is obtained; The total maintenance time is calculated by summing the products of the probabilities of various damage states and the repair times corresponding to the damage states.

6. A shield tunnel seismic toughness assessment method according to claim 5, characterized in that: The functional loss of the structure is calculated as: Among them, S loss represents the direct functional loss of the tunnel under a given specific IM situation; R j represents the repair ratio, the percentage required to repair the jth damage state; P[DS j |IM] represents the probability of the tunnel being located in damage state j and earthquake intensity level IM; The total maintenance time is calculated as follows: Among them, T RE represents the required repair time of the tunnel under a given specific IM; T REj represents the required repair time of the tunnel structure under different damage conditions; P[DS j |IM] represents the probability of the tunnel being located in damage state j and earthquake intensity level IM.

7. A shield tunnel seismic toughness assessment method according to claim 1, characterized in that: Determining the tunnel function function curve according to the structural recovery information and in combination with the functional recovery function of the structure, and then determining the functional recovery curve, comprises the following steps: According to the importance of the structure and the repair resources, a corresponding function recovery function is selected, wherein the function recovery function includes: a linear function, a trigonometric function and an exponential function; According to the function recovery function, combined with the function loss and the repair time, a tunnel function function curve is drawn; After determining the functional loss and repair time, the functional recovery curve under a given strength index is drawn based on the functional recovery function.

8. A shield tunnel seismic toughness assessment method according to claim 1, characterized in that: The step of evaluating the tunnel toughness according to the functional recovery curve to obtain an evaluation result of the tunnel seismic toughness comprises the following steps: By formula The seismic toughness of the shield tunnel under the main aftershocks was evaluated, where R represents the seismic toughness of the shield tunnel under the main aftershocks; T RE It indicates the time required for the structure to recover to its full function after the earthquake damage. 0E It indicates the time when the structure is damaged; Q(t) indicates the tunnel function curve; Determine the seismic resilience index of the tunnel when multiple earthquakes occur; According to the calculated shield tunnel seismic toughness index, the tunnel seismic toughness is quantitatively evaluated and divided into three levels according to the seismic toughness index. Determine the corresponding toughness assessment level based on the response results obtained.

9. A shield tunnel seismic toughness assessment device, characterized in that: include: The first module is used to create a mesh, define material elements and set boundaries to build a tunnel finite element model; The second module is used to select a seismic motion-input model and perform seismic vulnerability analysis after obtaining the response according to the tunnel finite element model; A third module is used to determine structural restoration information according to the results of the seismic vulnerability analysis; The structural restoration information includes structural function loss, tunnel restoration cost and restoration time; The fourth module is used to determine the tunnel function function curve according to the structural recovery information and the functional recovery function of the structure, and then determine the functional recovery curve; The fifth module is used to evaluate the toughness of the tunnel according to the functional recovery curve to obtain an evaluation result of the seismic toughness of the tunnel.

10. An electronic device, characterized in that: including a processor and a memory; The memory is used to store programs; The processor executes the program to implement the method according to any one of claims 1 to 8.

Images