State-dependent toughness calculation method and system for earthquake-damaged structures considering hysteresis time

Abstract

This invention provides a method and system for calculating the state-related toughness of earthquake-damaged structures considering hysteresis time. The method includes: selecting suitable mainshock and aftershock earthquake records and finite element modeling of the engineering structure; sequentially applying mainshock and aftershock records to the finite element model of the engineering structure to obtain damage data; calculating the hysteresis time caused by earthquake damage based on the structure's damage state; calculating the aftershock vulnerability surface based on the mainshock damage state of the structure; incorporating the hysteresis time into the toughness assessment formula, proposing a formula for calculating the state-related toughness of earthquake-damaged structures considering hysteresis time, and calculating the corresponding toughness surface. This method for calculating the state-related toughness of earthquake-damaged structures considering hysteresis time can be used to establish calculation models for the post-earthquake hysteresis time and corresponding toughness indices of engineering structures, evaluate the dynamic toughness level of target structures in real time, and provide a scientific reference for urban managers to formulate post-disaster relief plans.

Description

Technical Field

[0001] The present invention belongs to the technical field of seismic toughness assessment, and in particular relates to a method and system for calculating toughness associated with a state of earthquake-damaged structures taking into account hysteresis time. Background Art

[0002] A large number of earthquake events have shown that a strong earthquake is accompanied by multiple aftershocks. Earthquake records in recent years show that the frequency of strong earthquakes has increased year by year. Engineering structures can be severely damaged by strong earthquakes, and even directly cause the collapse of the structure, which poses a great hidden danger to my country's economic development and the safety of people's lives and property.

[0003] In a real earthquake, estimating the delay time caused by structural damage is a complex task that requires comprehensive consideration of various influencing factors. Since it is difficult to accurately assess the danger level of the structure after the earthquake, it is difficult to estimate the delay time and the impact of the delay time on the structural toughness. Most of the current toughness assessment methods simply estimate the repair time of structural components and regard the repair time of the structure as the delay time. This approach is undoubtedly very incomplete and inaccurate. It is impossible to calculate the delay time caused by post-disaster mobilization activities after the earthquake. This defect will undoubtedly significantly underestimate the delay time caused by the earthquake, which will lead to deviations in the results of the seismic toughness assessment of the structure. The inaccuracy of the calculation results of the seismic toughness of engineering structures will greatly hinder the formulation of rescue plans and the implementation of rescue tasks. This is undoubtedly an urgent problem facing urban managers and post-disaster rescue plan makers.

[0004] In view of this, it is very necessary to comprehensively calculate the post-earthquake hysteresis time and incorporate the hysteresis time into the seismic toughness evaluation formula to consider the impact of the hysteresis time on the seismic toughness level of the structure. Therefore, a method for calculating the state-related toughness of earthquake-damaged structures that takes into account the hysteresis time is urgently needed. Summary of the invention

[0005] The purpose of the present invention is to solve the problem that the existing earthquake resilience assessment cannot fully consider the hysteresis time caused by various post-earthquake rescue work, and propose a method for calculating the state-related resilience of earthquake-damaged structures that considers the hysteresis time. Based on the damage data of the structure, the nominal hysteresis time and the non-nominal hysteresis time are calculated respectively, and the key factor of the hysteresis time is introduced into the framework of seismic resilience analysis, providing a more comprehensive method for estimating the post-earthquake disaster relief time in actual engineering and evaluating the seismic resilience considering the hysteresis time, and the method is applied to the dynamic evaluation of the seismic performance of earthquake-damaged structures.

[0006] On the one hand, the present invention provides a method for calculating the state-related toughness of a damaged structure taking into account a hysteresis time, comprising the following steps:

[0007] S1: Selection of appropriate mainshock and aftershock earthquake records and finite element modeling of engineering structures;

[0008] S2: Apply the main shock and aftershock records to the finite element model of the engineering structure in S1 in sequence to obtain the damage data of the structure;

[0009] S3: Based on the main shock damage data and main aftershock cumulative damage data of S2, according to the damage state of the structure, calculate the hysteresis time of the structure caused by earthquake damage;

[0010] S4: Classify the damage data of the structure in S2 under the main shock, and calculate the aftershock vulnerability surface based on the main and aftershock cumulative damage data under different main shock damage states;

[0011] S5: Based on the calculation results of S3 and S4, the hysteresis time is incorporated into the toughness evaluation formula, and a calculation formula for the toughness related to the state of the damaged structure considering the hysteresis time is proposed, and the corresponding toughness surface is calculated.

[0012] Preferably, the S1 specifically includes:

[0013] Based on the public seismic database, appropriate main shock and aftershock records are selected, and the finite element model of the target engineering structure is established.

[0014] Preferably, the S2 specifically includes:

[0015] Based on the mainshock and aftershock records obtained by S1 and the finite element model of the target engineering structure, the mainshock and aftershock records are first input into the structural finite element model for nonlinear time history analysis. Then, the mainshock damage data and main-aftershock cumulative damage data of the structure are obtained according to the finite element analysis results.

[0016] Preferably, the S3 specifically includes:

[0017] Based on the mainshock and aftershock records obtained by S1 and the finite element model of the target engineering structure, the mainshock and aftershock records are first input into the structural finite element model for nonlinear time history analysis. Then, the mainshock damage data and main-aftershock cumulative damage data of the structure are obtained according to the finite element analysis results.

[0018] Preferably, the delay time (R Downtime ), as shown below:

[0019]

[0020] In the formula, R r and R ir Represents the rated delay time and non-rated delay time respectively; R t,Rebuild and R t,Repair R represents the structural reconstruction time and the earthquake-damaged structure repair time respectively; t,M Represents the post-disaster mobilization time; RCOT Indicates the material allocation time, which can be 14 days.

[0021] Based on the total probability theorem, the rated delay time can be calculated by the following formula:

[0022]

[0023] In the formula, E t,Rebuild and E t,Repair Represent the mathematical expectation of the structural reconstruction time and the damaged structure repair time respectively; P(Rebuild|IM AS ,t,DS MS ) and P(Repair|IM AS ,t,DS MS ) represent the structure in service time t and main shock damage state DS MS When the intensity is IM AS The probability of needing reconstruction and repair under aftershocks. A log-linear relationship is used to fit the mathematical relationship between the reconstruction time and the aftershock intensity and time, as shown in the following formula:

[0024] ln(E t,Repair )=α1ln(IM AS )+α2t+α3

[0025] Where α1, α2 and α3 are fitting coefficients.

[0026] Based on the structural damage status obtained by S2, it is divided into three levels from low to high in terms of danger, represented by colors, namely green, yellow and red. According to the classification level, the post-disaster mobilization time is calculated as shown in the following formula:

[0027] R t,M =E(R t,M |TAG G )P(TAG G |IM AS ,t)

[0028] +E(R t,M |TAG Y )P(TAG Y |IM AS ,t)+E(R t,M |TAG R )P(TAG R |IM AS ,t)

[0029] In the formula, E(R t,M |TAG G )、E(R t,M |TAG Y ) and E(Rt,M |TAG R ) represent green (TAG G )、Yellow(TAG Y ) and Red Level (TAG R ) level; P(TAG|·) represents the probability of reaching the corresponding danger level.

[0030] Preferably, the S4 specifically includes:

[0031] Based on the structural damage state obtained by S2, the earthquake-damaged structures are classified according to the main shock damage results of the structure, and the probabilistic earthquake demand model is established for different main shock damage conditions, as shown in the following formula:

[0032] ln[D(IM AS ,t,Θ|DS MS )]=θ1ln(IM AS )+h(t,θ2)+θ3

[0033]

[0034] h(t,θ2)=θ 2,0 +θ 2,1 t+θ 2,2 t 2 +θ 2,3 t 3 +…

[0035] Where D(·) represents the cumulative damage of the structure; θ1, θ2, and θ3 are model parameters, and θ = (θ1, θ2, θ3); h(·) is a function of time t; θ2 represents a parameter vector, and θ2 = (θ 2,0 ,θ 2,1 ,θ 2,2 …); EDP MA (·) represents the cumulative structural damage parameter; N represents the number of aftershock records.

[0036] Based on the obtained probabilistic earthquake demand model, the aftershock vulnerability surface is calculated using the following formula:

[0037]

[0038] In the formula, represents the seismic capacity parameter of the structure; β represents the logarithmic standard deviation.

[0039] Preferably, the calculation formula for the toughness related to the state of the damaged structure considering the hysteresis time in step 6 is as follows:

[0040]

[0041] Q(IMAS ,t|DS MS )=

[0042] 1-L(IM AS ,t|DS MS )·[H(t-t0)-H(t-t0+R Downtime )]·f(t,t0,R Downtime )

[0043] Where R(·) is the toughness index; IM AS Indicates aftershock intensity; DS MS The damage state of the main shock; t is the service life of the structure; t0 represents the time when the earthquake occurred; T LC To control time; Q(IM AS ,T|DS MS ) represents the performance function; L(IMAS,t|DSMS) represents the loss function, which can be expressed by P(LS j |IM AS ,t,DS MS ) is calculated; R Downtime is the hysteresis time, H(·) and f(·) are the Heaviside step function and recovery function respectively.

[0044] The second aspect of the present invention provides a system for calculating toughness related to the state of earthquake-damaged structures taking into account hysteresis time, comprising:

[0045] Selection module, used to select appropriate main shock and aftershock earthquake records and finite element modeling of engineering structures;

[0046] An acquisition module is used to apply the main shock and aftershock records to the finite element model of the engineering structure in sequence to obtain the damage data of the structure;

[0047] The first calculation module calculates the hysteresis time of the structure caused by earthquake damage based on the main shock damage data and the main aftershock cumulative damage data according to the damage state of the structure;

[0048] The second calculation module is used to classify the damage data of the structure under the main shock, and calculate the aftershock vulnerability surface based on the main and aftershock cumulative damage data under different main shock damage states;

[0049] The third calculation module, based on the calculation results of the first calculation module and the second calculation module, incorporates the hysteresis time into the toughness evaluation formula, proposes a toughness calculation formula related to the state of the damaged structure considering the hysteresis time, and calculates the corresponding toughness surface.

[0050] A third aspect of the present invention provides a non-transitory computer-readable storage medium, which stores computer instructions. When the computer instructions are executed by a processor, they implement the method for calculating the toughness related to the state of earthquake-damaged structures taking into account the hysteresis time as described in the first aspect.

[0051] A fourth aspect of the present invention provides an electronic device, comprising: a memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the method for calculating toughness related to the state of earthquake-damaged structure taking into account the hysteresis time as described in the first aspect by executing the computer instructions.

[0052] The method, system, electronic device and computer-readable storage medium provided by the present invention have the following beneficial effects:

[0053] The present invention first divides the hysteresis time into two categories: quota hysteresis time and non-quota hysteresis time. For the post-disaster mobilization time which is the most difficult to calculate among the non-quota hysteresis times, the present invention adopts the total probability theorem to comprehensively consider the three structural danger levels of green, yellow and red to calculate the post-disaster mobilization time, and fully considers the contribution of each danger level to the post-disaster mobilization time, so as to comprehensively evaluate the various time hysteresis phenomena caused by earthquakes. Finally, the present invention incorporates the key concept of hysteresis time into the existing resilience assessment method, solving the problem that the existing resilience assessment method cannot comprehensively assess the impact of time hysteresis on structural seismic resilience. It can more comprehensively assess the hysteresis time and toughness level of engineering structures under different damage states. The earthquake-damaged structural state-related toughness measurement method considering the hysteresis time of the present invention can be used to establish a calculation model for the post-earthquake hysteresis time and corresponding toughness index of engineering structures, evaluate the dynamic toughness level of the target structure in real time, and provide a scientific reference for urban managers to formulate post-disaster rescue plans. BRIEF DESCRIPTION OF THE DRAWINGS

[0054] Figure 1 It is a flow chart of the method for calculating toughness related to the state of earthquake-damaged structure considering the hysteresis time of the present invention;

[0055] Figure 2 is the mathematical expectation of the reconstruction time and the repair time of the damaged structure in one embodiment of the present invention;

[0056] In the figure, Figure 2 (a) The main shock damage condition is DS0, Figure 2 (b) The main shock damage condition is DS1; Figure 2 (c) The main shock damage condition is DS2, Figure 2 (d) The main shock damage condition is DS3 (corresponding to the post-earthquake reconstruction condition);

[0057] Figure 3is the calculation result of the total hysteresis time in one embodiment of the present invention;

[0058] In the figure, Figure 3 (a) The main shock damage condition is DS0, Figure 3 (b) The main shock damage condition is DS1; Figure 3 (c) the main shock damage condition is DS2;

[0059] Figure 4 is the calculation result of the aftershock vulnerability surface in one embodiment of the present invention;

[0060] In the figure, Figure 4 (a) The main shock damage condition is DS0, Figure 4 (b) The main shock damage condition is DS1; Figure 4 (c) The main shock damage condition is DS2, Figure 4 (d) the main shock damage condition is DS3;

[0061] Figure 5 is a toughness surface related to the state of the damaged structure considering the hysteresis time in one embodiment of the present invention;

[0062] In the figure, Figure 5 (a) The main shock damage condition is DS0, Figure 5 (b) The main shock damage condition is DS1; Figure 5 (c) the main shock damage condition is DS2;

[0063] Figure 6 is an architecture diagram of a system for calculating toughness associated with earthquake-damaged structural states taking into account hysteresis time according to a preferred embodiment of the present invention;

[0064] Figure 7 A schematic structural diagram of an embodiment of an electronic device provided by the present invention. DETAILED DESCRIPTION

[0065] The specific implementation of the present invention is further described in detail below in conjunction with the accompanying drawings and examples. The following examples are used to illustrate the present invention, but are not intended to limit the scope of the present invention.

[0066] Embodiment 1

[0067] like Figure 1 As shown, the method for calculating the state-related toughness of earthquake-damaged structures taking into account the hysteresis time provided in this embodiment incorporates the hysteresis time into the seismic toughness evaluation formula, and considers the influence of the hysteresis time on the seismic toughness level of the structure, so as to more comprehensively evaluate the real-time toughness level of the post-earthquake engineering structure.

[0068] The method for calculating the state-related toughness of earthquake-damaged structures considering the hysteresis time is as follows:

[0069] like Figure 1 As shown, the method for calculating toughness related to the state of earthquake-damaged structure considering the hysteresis time in this embodiment includes:

[0070] Step S1: Based on the public earthquake motion database, select appropriate main shock and aftershock records, and establish a finite element model of the target engineering structure.

[0071] Step S2: applying the main shock and aftershock records to the finite element model of the engineering structure in S1 in sequence to obtain damage data of the structure.

[0072] Step S3, based on the main shock damage data and the main aftershock cumulative damage data of S2, according to the damage state of the structure, calculate the hysteresis time of the structure caused by the earthquake damage.

[0073] Specifically in this embodiment, the above step S3 also includes the following steps:

[0074] Hysteresis time (R Downtime ), the calculation formula is as follows:

[0075]

[0076] In the formula, R r and R ir Represents the rated delay time and non-rated delay time respectively; R t,Rebuild and R t,Repair R represents the structural reconstruction time and the earthquake-damaged structure repair time respectively; t,M Represents the post-disaster mobilization time; R COT Indicates the material allocation time, which can be 14 days.

[0077] Based on the total probability theorem, the rated delay time can be calculated by the following formula:

[0078]

[0079] In the formula, E t,Rebuild and E t,Repair Represent the mathematical expectation of the structural reconstruction time and the damaged structure repair time respectively; P(Rebuild|IM AS ,t,DS MS ) and P(Repair|IM AS ,t,DS MS ) represent the structure in service time t and main shock damage state DS MS When the intensity is IM AS The probability of needing reconstruction and repair under aftershocks can be obtained from aftershock vulnerability analysis. In this embodiment, the aftershock peak acceleration (PGA) is used. AS ) to indicate IM ASThe mathematical expectation of the reconstruction time and the repair time of the damaged structure can be calculated using the current "Building Seismic Resilience Evaluation Standard" in my country. A logarithmic linear relationship is used to fit the mathematical relationship between the reconstruction time and the aftershock intensity and time, as shown in the following formula:

[0080] ln(E t,Repair )=α1ln(PGA AS )+α2t+α3

[0081] Where α1, α2 and α3 are fitting coefficients. The calculation results of reconstruction time and damage structure repair time are as follows: Figure 2 shown.

[0082] Based on the structural damage status obtained by S2, it is divided into three levels from low to high in terms of danger, represented by colors, namely green, yellow and red. According to the classification level, the post-disaster mobilization time is calculated as shown in the following formula:

[0083] R t,M =E(R t,M |TAG G )P(TAG G |IM AS ,t)

[0084] +E(R t,M |TAG Y )P(TAG Y |IM AS ,t)+E(R t,M |TAG R )P(TAG R |IM AS ,t)

[0085] In the formula, E(R t,M |TAG G )、E(R t,M |TAG Y ) and E(R t,M |TAG R ) represent green (TAG G )、Yellow(TAG Y ) and Red Level (TAG R ) level; P(TAG|·) represents the probability of reaching the corresponding danger level. The total delay time calculation result is as follows: Figure 3 shown.

[0086] Step S4: classify the damage data of the structure in S2 under the main shock, and calculate the aftershock vulnerability surface based on the main shock and aftershock cumulative damage data under different main shock damage states.

[0087] Specifically in this embodiment, the above step S4 also includes the following steps:

[0088] The probabilistic earthquake demand model is established for different main shock damage conditions, as shown in the following formula:

[0089] ln[D(IM AS ,t,Θ|DS MS )]=θ1ln(IM AS )+h(t,θ2)+θ3

[0090]

[0091] h(t,θ2)=θ 2,0 +θ 2,1 t+θ 2,2 t 2 +θ 2,3 t 3 +…

[0092] Where D(·) represents the cumulative damage of the structure; θ1, θ2, and θ3 are model parameters, and θ = (θ1, θ2, θ3); h(·) is a function of time t; θ2 represents a parameter vector, and θ2 = (θ 2,0 ,θ 2,1 ,θ 2,2 …); EDP MA (·) represents the cumulative damage parameter of the structure; N represents the number of aftershock records. In this embodiment, the Park-Ang damage index is used to represent EDP MA The calculation results of aftershock vulnerability surface are as follows: Figure 4 shown.

[0093] Step S5: Based on the calculation results of S3 and S4, the hysteresis time is incorporated into the toughness evaluation formula, a calculation formula for the toughness related to the state of the damaged structure considering the hysteresis time is proposed, and the corresponding toughness surface is calculated.

[0094] Specifically in this embodiment, the above step S5 also includes the following steps:

[0095] Based on the calculation results of S3 and S4, the calculated hysteresis time is substituted into the toughness evaluation formula, and a calculation formula for the state-related toughness of the damaged structure considering the hysteresis time is proposed, as shown in the following formula:

[0096]

[0097] Q(IM AS ,t|DS MS )=

[0098] 1-L(IM AS ,t|DS MS)·[H(t-t0)-H(t-t0+R Downtime )]·f(t,t0,R Downtime )

[0099] Where R(·) is the toughness index; IM AS Indicates aftershock intensity; DS MS The damage state of the main shock; t is the service life of the structure; t0 represents the time when the earthquake occurred; T LC To control time; Q(IM AS ,T|DS MS ) represents the functional function; L(IM AS ,t|DS MS ) represents the loss function, which can be expressed by P(LS j |IM AS ,t,DS MS ) is calculated; R Downtime H(·) is the hysteresis time and f(·) is the Heaviside step function and recovery function respectively. The calculation results of the toughness surface related to the state of the damaged structure considering the hysteresis time are shown in Figure 5 shown.

[0100] The method for calculating the state-related toughness of earthquake-damaged structures taking into account the hysteresis time provided in this embodiment can be used to evaluate the hysteresis time for the repair or reconstruction of earthquake-damaged structures, and to establish a corresponding aftershock toughness dynamic analysis model, thereby providing a theoretical basis for the toughness assessment of engineering structures.

[0101] Embodiment 2

[0102] like Figure 6 As shown in the figure, the toughness calculation system of earthquake-damaged structure state considering hysteresis time includes:

[0103] A selection module 101 is used to select appropriate main shock and aftershock earthquake records and finite element modeling of engineering structures;

[0104] An acquisition module 102 is used to apply the main shock and aftershock records to the finite element model of the engineering structure in sequence to obtain damage data of the structure;

[0105] The first calculation module 103 calculates the delay time of the structure caused by earthquake damage based on the main shock damage data and the main aftershock cumulative damage data according to the damage state of the structure;

[0106] The second calculation module 104 is used to classify the damage data of the structure under the main shock, and calculate the aftershock vulnerability surface according to the main shock and aftershock cumulative damage data under different main shock damage states;

[0107] The third calculation module 105, based on the calculation results of the first calculation module and the second calculation module, incorporates the hysteresis time into the toughness evaluation formula, proposes a toughness calculation formula related to the damaged structure state considering the hysteresis time, and calculates the corresponding toughness surface.

[0108] This embodiment incorporates the key concept of hysteresis time into the existing resilience assessment method, solving the problem that the existing resilience assessment method cannot comprehensively assess the various time hysteresis phenomena caused by earthquakes. It can more comprehensively assess the hysteresis time and resilience level caused by repair or reconstruction of engineering structures in different damage states. It can evaluate the dynamic resilience level of the target structure in real time, providing a scientific reference for city managers to formulate post-disaster rescue plans.

[0109] An embodiment of the present invention also provides a non-transitory computer-readable storage medium, which stores computer instructions. When the computer instructions are executed by a processor, a method for calculating toughness related to the state of earthquake-damaged structures taking into account the hysteresis time as described in Example 1 is implemented.

[0110] like Figure 7 As shown, the present invention also provides an electronic device, including: a memory 302 and a processor 301, the memory 302 and the processor 301 are connected to communicate with each other, the memory 302 stores computer instructions, and the processor 301 executes the method for calculating toughness related to the state of earthquake-damaged structure considering the hysteresis time as in Example 1 by executing the computer instructions.

[0111] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. It should be pointed out that a person skilled in the art can make several improvements and modifications without departing from the technical principles of the present invention, and these improvements and modifications should also be regarded as within the scope of protection of the present invention.

Claims

1. A method for calculating the state-related toughness of earthquake-damaged structures considering hysteresis time, characterized in that: The steps include: S1: Selection of main shock and aftershock earthquake records and finite element modeling of engineering structures; S2: Apply the main shock and aftershock records to the finite element model of the engineering structure in S1 in sequence to obtain the damage data of the structure; S3: Based on the main shock damage data and main aftershock cumulative damage data of S2, according to the damage state of the structure, calculate the hysteresis time of the structure caused by earthquake damage; The S3 specifically includes: Based on the structural damage state obtained by S2, calculate the hysteresis time R caused by structural damage Downtime , the calculation formula is as follows In the formula, R r and R ir Represents the rated delay time and non-rated delay time respectively; R t,Rebuild and R t,Repair R represents the structural reconstruction time and the earthquake-damaged structure repair time respectively; t,M Represents the post-disaster mobilization time; R COT Indicates the material allocation time, which is 14 days; Based on the total probability theorem, the rated delay time is calculated by the following formula: In the formula, E t,Rebuild and E t,Repair Represent the mathematical expectation of the structural reconstruction time and the damaged structure repair time respectively; P(Rebuild|IM AS ,t,DS MS ) and P(Repair|IM AS ,t,DS MS ) represent the structure in service time t and main shock damage state DS MS When the intensity is IM AS The probability of needing reconstruction and repair under aftershocks; a log-linear relationship is used to fit the mathematical relationship between the reconstruction time and the aftershock intensity and time, as shown in the following formula: ln(E t,Repair )=α1ln(IM AS )+α2t+α3 Where α1, α2 and α3 are fitting coefficients; The structural damage status obtained based on S2 is divided into three levels from low to high in terms of danger, represented by colors, namely green, yellow and red. According to the classification level, the post-disaster mobilization time is calculated as shown in the following formula: R t,M =E(R t,M |DAY G )P(DAY G |IM AS ,t) +E(R t,M |DAY Y )P(DAY Y |IM AS ,t)+E(R t,M |DAY R )P(DAY R |IM AS ,t) In the formula, E(R t,M |TAG G )、E(R t,M |TAG Y ) and E(R t,M |TAG R ) represent green tags respectively G , yellow tag Y and red tag R The mathematical expectation of post-disaster mobilization time under the level condition; P(TAG|·) represents the probability of reaching the corresponding danger level; S4: Classify the damage data of the structure in S2 under the main shock, and calculate the aftershock vulnerability surface based on the main and aftershock cumulative damage data under different main shock damage states; S5: Based on the calculation results of S3 and S4, the hysteresis time is incorporated into the toughness evaluation formula, and a calculation formula for the toughness related to the state of the damaged structure considering the hysteresis time is proposed, and the corresponding toughness surface is calculated; The S5 specifically includes: Based on the calculation results of S3 and S4, the calculated hysteresis time is substituted into the toughness evaluation formula, and a calculation formula for the state-related toughness of the damaged structure considering the hysteresis time is proposed, as shown in the following formula: Q(IM AS ,t|DS MS )= 1-L(IM AS ,t|DS MS )·[H(t-t0)-H(t-t0+R Downtime )]·f(t,t0,R Downtime ) Where R(·) is the toughness index; IM AS Indicates aftershock intensity; DS MS The damage state of the main shock; t is the service life of the structure; t0 represents the time when the earthquake occurred; T LC To control time; Q(IM AS ,T|DS MS ) represents the functional function; L(IM AS ,t|DS MS ) represents the loss function, P(LS j |IM AS ,t,DS MS ) is calculated; R Downtime is the hysteresis time, H(·) and f(·) are the Heaviside step function and recovery function respectively.

2. The method for calculating toughness related to the state of earthquake-damaged structure considering hysteresis time according to claim 1 is characterized in that: The S1 specifically includes: Based on the public seismic motion database, main shock and aftershock records are selected; and a finite element model of the target engineering structure is established.

3. The method for calculating the state-related toughness of earthquake-damaged structures considering hysteresis time according to claim 1 is characterized in that: The S2 specifically includes: Based on the mainshock and aftershock records obtained by S1 and the finite element model of the target engineering structure, the mainshock and aftershock records are first input into the structural finite element model for nonlinear time history analysis. Then, the mainshock damage data and main-aftershock cumulative damage data of the structure are obtained according to the finite element analysis results.

4. The method for calculating toughness related to the state of earthquake-damaged structure considering hysteresis time according to claim 1 is characterized in that: The S4 specifically includes: Based on the structural damage state obtained by S2, the earthquake-damaged structures are classified according to the main shock damage results of the structure, and the probabilistic earthquake demand model is established for different main shock damage conditions, as shown in the following formula: ln[D(IM AS ,t,Θ|DS MS )]=θ1ln(IM AS )+h(t,θ2)+θ3 h(t,θ2)=θ 2,0 +θ 2,1 t+θ 2,2 t 2 +θ 2,3 t 3 +… Where D(·) represents the cumulative damage of the structure; θ1, θ2, and θ3 are model parameters, and θ = (θ1, θ2, θ3); h(·) is a function of time t; θ2 represents the parameter vector, and θ2 = (θ 2,0 ,θ 2,1 ,θ 2,2 …); EDP MA (·) represents the cumulative structural damage parameter; N represents the number of aftershock records; Based on the obtained probabilistic earthquake demand model, the aftershock vulnerability surface is calculated using the following formula: In the formula, represents the seismic capacity parameter of the structure; β represents the logarithmic standard deviation.

5. A system for calculating toughness related to the state of earthquake-damaged structures taking into account hysteresis time, applied to the method according to any one of claims 1 to 4, characterized in that: include: Selection module, used to select main shock and aftershock seismic records and finite element modeling of engineering structures; An acquisition module is used to apply the main shock and aftershock records to the finite element model of the engineering structure in sequence to obtain the damage data of the structure; The first calculation module calculates the hysteresis time of the structure caused by earthquake damage based on the main shock damage data and the main aftershock cumulative damage data according to the damage state of the structure; The second calculation module is used to classify the damage data of the structure under the main shock, and calculate the aftershock vulnerability surface based on the main and aftershock cumulative damage data under different main shock damage states; The third calculation module, based on the calculation results of the first calculation module and the second calculation module, incorporates the hysteresis time into the toughness evaluation formula, proposes a toughness calculation formula related to the state of the damaged structure considering the hysteresis time, and calculates the corresponding toughness surface.

6. A non-transitory computer-readable storage medium, characterized in that: The non-transitory computer-readable storage medium stores computer instructions, which, when executed by a processor, implement the method for calculating toughness related to the state of a damaged structure taking into account the hysteresis time as described in any one of claims 1 to 4.

7. An electronic device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the method for calculating toughness related to the state of earthquake-damaged structure taking into account the hysteresis time as described in any one of claims 1 to 4 by executing the computer instructions.

Images