Method for evaluating seismic resilience of medium and large railway passenger station building

Abstract

The application discloses a kind of middle and large railway passenger station building seismic toughness evaluation method, using the engineering demand parameter matrix expansion algorithm based on FEMA-P58, toughness index calculation algorithm based on Monte Carlo sampling simulation, and formulates the post-earthquake repair strategy of "first structure, then enclosure, then electromechanical", form toughness evaluation grade determination and post-earthquake function recovery simulation based on repair time, finally realizes the quantitative evaluation of middle and large railway passenger station building seismic toughness, improves the feasibility and accuracy of railway station seismic toughness evaluation.

Description

Technical Field

[0001] This invention relates to the field of building structure monitoring and evaluation technology, specifically to a method for evaluating the seismic toughness of medium and large-sized railway passenger station buildings. Background Technology

[0002] Medium and large railway passenger stations, as key nodes in the transportation network, play a vital role in passenger flow, material transport, and train scheduling. Ensuring their post-earthquake functionality is particularly crucial for earthquake relief and disaster recovery. Historical earthquake damage shows that damage to passenger stations can lead to prolonged disruptions in transportation services and trigger severe economic and social consequences. Therefore, to ensure the post-earthquake recovery capabilities of railway passenger stations and mitigate the chain reaction of disasters caused by functional disruptions, the stations must possess sufficient seismic resilience.

[0003] Currently, several standards and specifications are used for seismic toughness assessment of typical civil buildings, including FEMA-P58, the REDi rating system, the USRC building rating system, and the "Standard for Seismic Toughness Evaluation of Buildings" (GB / T 38591-2020). However, it should be noted that these standards and specifications are mainly applicable to typical civil buildings. For buildings with special functions, such as railway passenger stations, hospitals, and substations, these standards cannot be directly applied to the quantitative analysis of their seismic toughness. Therefore, considering the special nature of building functions, in recent years, scholars both domestically and internationally have further conducted research on seismic toughness assessment methods for buildings with specific functions. These research results have laid a solid foundation for the seismic toughness assessment of both conventional building structures and buildings with specific functions.

[0004] However, railway stations possess unique architectural structural characteristics, functional support equipment, and complex functional hierarchy transmission mechanisms and post-earthquake recovery processes, making existing resilience assessment methods unsuitable for assessing the seismic resilience of railway stations. On one hand, the lack of dedicated databases and recovery strategies weakens the assessment foundation: existing assessment systems lack vulnerability databases for special components of railway passenger stations (large-span steel space frames, steel-concrete composite columns, etc.) and specialized equipment (such as signaling systems and power supply systems). They also lack standardized functional recovery strategies and resource scheduling models that align with their operational characteristics, resulting in inaccurate inputs and reliable evidence for the assessment. On the other hand, the core assessment process is complex, and functional dependencies are insufficiently characterized: the generation, expansion, and verification of the engineering requirement parameter matrix are cumbersome, and existing methods struggle to quantify and simulate the complex hierarchical dependencies and transmission relationships between "structure-equipment-function" within railway passenger stations. This leads to a break in the assessment chain from component damage to system functional loss, affecting the accuracy of the overall resilience assessment. Summary of the Invention

[0005] This application provides a method for assessing the seismic toughness of medium and large-sized railway passenger station buildings, in order to solve the problem that existing toughness assessment methods are difficult to apply to the seismic toughness assessment of railway stations.

[0006] To achieve the above objectives, the technical solution adopted by the present invention is as follows:

[0007] A method for assessing the seismic toughness of medium-to-large railway passenger station buildings is provided, the method comprising:

[0008] Based on the finite element analysis of railway passenger station buildings under seismic conditions, the engineering demand parameter matrix of railway passenger station buildings is obtained. By expanding and verifying the engineering demand parameter matrix of railway passenger station buildings, the expanded engineering demand parameter matrix is ​​obtained.

[0009] Based on the Building Information Modeling (BIM) of railway passenger station buildings, obtain railway passenger station building information model data including the number of station building components or equipment, station building room area, and engineering quantity;

[0010] A vulnerability parameter database for railway passenger station building components or equipment is constructed. The vulnerability parameter database contains the expected values ​​and variances of engineering requirement parameters for different components or equipment under different damage states.

[0011] The post-earthquake repair strategy of "structure first, then enclosure, then electromechanical" was formulated, which means repairing the structural components and vertical transportation components, the enclosure structure and water supply and drainage system, the heating, ventilation and air conditioning system and power and main signal equipment, and the secondary signal equipment in four stages in sequence.

[0012] Based on the established post-earthquake repair strategy, and according to the expanded engineering demand parameter matrix, building information model data, and vulnerability parameter database, the Monte Carlo simulation function was used to calculate the resilience assessment indicators of railway passenger station buildings, including repair time and post-earthquake function.

[0013] Based on the obtained resilience assessment index, a post-earthquake functional recovery function model that varies with repair time is fitted and established, and a post-earthquake functional recovery curve is plotted. The seismic resilience index is calculated based on the area of ​​the post-earthquake functional recovery curve obtained by integration, and the seismic resilience assessment level is classified according to the seismic resilience index.

[0014] Furthermore, based on finite element analysis of railway passenger station buildings under seismic conditions, the engineering demand parameter matrix of railway passenger station buildings is obtained, specifically including:

[0015] A finite element model of the railway passenger station structure was established, and seismic conditions were set to perform elastoplastic time history analysis on the railway passenger structure. After the calculation, the inter-story drift angle and peak floor acceleration of each floor under each seismic condition were extracted, and a parameter matrix of engineering requirements for the railway passenger station was constructed. , where m is the row of the matrix, each row corresponds to a seismic condition; n is the column of the matrix, each column corresponds to an engineering requirement parameter, which includes the inter-story drift angle and peak floor acceleration of each floor.

[0016] Furthermore, the engineering demand parameter matrix for railway passenger station buildings is expanded and verified to obtain an expanded engineering demand parameter matrix, which specifically includes:

[0017] Read the original project requirement parameter matrix ;

[0018] Calculate the original engineering requirement parameter matrix covariance matrix and mean matrix ;

[0019] Calculate the covariance matrix rank R, eigenvalue matrix eigenvector matrix ;

[0020] The eigenvalue matrix is ​​divided according to rank R. eigenvector matrix ,get and Further segmentation of the eigenvalue matrix Square root and transform into a diagonal matrix ;

[0021] Determine the number of expansions x, and generate a dimension of Independent standard normal random variable matrix And based on the mean matrix Generate a dimension of Mean expansion matrix ;

[0022] Calculate the expanded engineering requirement parameter matrix using the following formula. :

[0023] .

[0024] Furthermore, based on the Building Information Model (BIM) of railway passenger station buildings, data including the quantity of station building components or equipment, the area of ​​station building rooms, and the amount of work completed are obtained. Specifically, this includes:

[0025] The station building components or equipment include:

[0026] Structural components include: reinforced concrete frame columns, reinforced concrete frame beams, reinforced concrete shear walls, reinforced concrete connecting beams, steel structure beams, steel structure columns, steel bracing components, steel-concrete composite columns, steel-concrete composite beams, steel space frames, and steel trusses.

[0027] Displacement-sensitive non-structural components include: infill walls, glass curtain walls, and staircases;

[0028] Acceleration-sensitive non-structural components include: ceilings, elevators, suspended lighting fixtures, switchgear, distribution boxes, water supply pipes, fire sprinkler pipes, sprinkler head risers, HVAC ducts, air outlets, HVAC fans, and air conditioning system fans;

[0029] Specialized equipment includes: battery cabinets, cabling systems, IFS cabinets (Interface and Frequency Synthesizer), IFSI cabinets (Interface and Frequency Synthesizer Interface), ISFS cabinets (Integrated Signal and Frequency Synthesizer), IFSSI cabinets (Integrated Signal and Frequency Synthesizer Interface), IRC cabinets (Interface and Route Controller), ISRC cabinets (Integrated Signal and Route Controller), train control cabinets, CTC (Centralized Traffic Control) cabinets, interlocking cabinets, RBC (Radio Block Center) cabinets, and TSRS (Temporary Speed ​​Restriction Server) cabinets. Server (temporary rate-limited server) rack.

[0030] Furthermore, a post-earthquake repair strategy of "structure first, then enclosure, then electromechanical" was formulated, which means repairing structural components and vertical transportation components, enclosure structure and water supply and drainage system, HVAC system and power and main signal equipment, and secondary signal equipment in four stages in sequence, specifically including:

[0031] The structural components and vertical transportation components include: reinforced concrete frame columns, reinforced concrete frame beams, reinforced concrete shear walls, reinforced concrete connecting beams, steel structure beams, steel structure columns, steel support components, steel-concrete composite columns, steel-concrete composite beams, steel space frames, steel trusses, stairs, and elevators.

[0032] The enclosure structure and water supply and drainage system include: infill walls, glass curtain walls, suspended ceilings, water supply pipes, fire sprinkler pipes, and sprinkler head risers;

[0033] The HVAC system and power and main signaling equipment include: suspended lighting fixtures, switchgear, distribution boxes, air outlets, HVAC ducts, HVAC fans, air conditioning system fans, battery cabinets, cable systems, ISFS cabinets, ISFSI cabinets, ISRC cabinets, train control cabinets, CTC cabinets, interlocking cabinets, RBC cabinets, and TSRS cabinets.

[0034] The secondary signaling equipment includes: IFS cabinet, IFSI cabinet, and IRC cabinet.

[0035] Furthermore, based on the established post-earthquake repair strategy, and according to the expanded engineering demand parameter matrix, building information model data, and vulnerability parameter database, Monte Carlo simulation functions were used to calculate resilience assessment indicators for railway passenger station buildings, including repair time and post-earthquake functionality. Fitted values ​​for each indicator were calculated at specified confidence levels, specifically including:

[0036] Step A: Define the number of damaged components or equipment. A floor index used to store the number of damaged components or equipment on each floor, initialized to 0, and used to initialize the inter-floor displacement angle and peak floor acceleration. , All are 0;

[0037] Step B: Traverse the expanded engineering requirement parameter matrix row by row. This yields the value in the i-th row and j-th column, i.e. ;

[0038] Step C, if If the inter-floor displacement angle is used, then the number of components or equipment on the current floor is traversed and filtered. The term "structural component" or "displacement-sensitive non-structural component" is used; if To determine the peak floor acceleration, the number of components or equipment on each floor is then iterated through and filtered. The "acceleration-sensitive non-structural components" or "specialized equipment" mentioned above;

[0039] Step D: If the number of components or equipment on the current floor is not zero after filtering, then iterate through the damage states of the corresponding components or equipment, and calculate the log-normal cumulative distribution probability of the p-th type of component or equipment under the d-th damage state using the norm.cdf function of the Python Scipy library. As shown in the following formula:

[0040] ;

[0041] in, Let f be the expected value of the p-th component or equipment in the fragility parameter database F under the d-th damage state; Let be the variance of the p-th component or equipment in the fragility parameter database F under the d-th damage state;

[0042] Step E: Determine the quantity of the p-th type of component or equipment in the current layer. Repeated random sampling Each time, a random number between 0 and 1 is generated. Based on the probability interval of the random number falling into each range, determine the number of components or equipment of type p in each damage state.

[0043] Step F, if If the inter-story drift angle is the value, then update DC. n The Middle Calculate the number of components or equipment of type p in the d-th damage state, and update the floor index. Repeat steps B through E; if For the peak floor acceleration, update The Middle Calculate the number of components or equipment of type p in the d-th damage state, and update the floor index. Repeat steps B through E;

[0044] Step G: Calculate the repair time, repair cost, casualties, and post-earthquake functionality of railway passenger station buildings based on the sampling results, and return the calculation results.

[0045] Furthermore, in step E, the probability interval division is specifically as follows:

[0046] If the probability If it contains 4 elements, it means that the damage state of the component is of type 4. Indicates no damage ; Indicates minor injury ; Indicates moderate damage ; Indicates severe injury ; Indicates complete destruction ;

[0047] if If it contains 3 elements, it indicates that the damage state of the component is of type 3. Indicates no damage ; Indicates minor injury ; Indicates moderate damage ; Indicates severe injury ;

[0048] if If it contains two elements, it indicates that the damage state of the component is of type 2. Indicates no damage ; Indicates minor injury ; Indicates severe injury ;

[0049] if If a component contains one element, it indicates that the damage state is type 1. Indicates no damage ; Indicates severe injury .

[0050] Furthermore, in step G, the specific steps for calculating the repair time are as follows:

[0051] G1, Read the number of damaged components or equipment Determine the repair time Q for components or equipment, the number of workers q required to repair a single component or piece of equipment, the station building area A, and the repair cost reduction factor. ;

[0052] G2. First, calculate the total time to complete the first phase of repair work. The details are as follows:

[0053] First, calculate the number of workers required to repair the vertical transportation components on the k-th floor. and total working hours The repair time is obtained by dividing the total working hours by the number of workers required. Where p represents the p-th type of component or equipment; P represents the total number of component or equipment types; d represents the d-th damage state; D represents the total number of damage state types; k represents the k-th floor; and K represents the total number of floors. The number of workers required to repair component or equipment of type p; This represents the number of the p-th type of component or equipment in the k-th layer under the d-th damage state. Repair time for component or equipment of type p under damage condition d;

[0054] Next, calculate the maximum allowed number of workers on the k-th floor (0.026A(k)), the required number of workers to repair structural components (0.02A(k)), and the total working hours for repairing structural components. Provided that the sum of the number of workers needed to repair structural components and the number of workers needed to repair vertical transportation components does not exceed the maximum number of workers allowed per floor, the time required to repair the structural components is obtained by dividing the total working hours by the number of workers needed. ;

[0055] Compare and extract and The larger of the two, serving as the total time for repairing structural components and vertical transportation components. Therefore, the total time to complete the first phase of repair work is... , recorded as ;in The total time for repairing structural components and vertical transportation components on the kth floor; This is the maximum value in the total time of the first phase of repair work, that is, the longest time required to repair the entire station building structure and vertical transportation components;

[0056] G3. Next, calculate the total time T2 for completing the first and second phases of repair work, as follows:

[0057] Calculate the number of workers required (0.01A(k)) and the total man-hours for repairing the k-th floor enclosure structure and water supply and drainage system. The repair time is obtained by dividing the total working hours by the number of workers required. Therefore, the total time to complete the first and second phases of repair work is... , recorded as ;in, This is the maximum value in the total time of the second phase of repair work, that is, the longest time required to repair the entire station building envelope and water supply and drainage system; It is the sum of the maximum value of the total time of the first phase of repair work and the maximum value of the total time of the second phase of repair work, that is, the longest time required to repair the entire station building structure components, vertical transportation components, enclosure structure and water supply and drainage system.

[0058] G4. Recalculate the total time required to complete the repair work in phases 1, 2, and 3. The details are as follows:

[0059] First, calculate the number of workers required to repair the power and main signal equipment on the k-th floor. and total working hours The repair time is obtained by dividing the total working hours by the number of workers required. ;

[0060] Next, calculate the maximum number of workers allowed to repair the k-th floor (0.026A(k)), the number of workers required to repair the HVAC system (0.01A(k)), and the total man-hours for repairing the HVAC system. Provided that the sum of the number of workers needed to repair the HVAC system and the number of workers needed to repair the electrical and main signaling equipment does not exceed the maximum number of workers allowed on the floor, the repair time for the HVAC system can be obtained by dividing the total working hours by the number of workers needed. ;

[0061] Comparison Extraction and The larger of the two, representing the total time required to repair the HVAC system and electrical and major signaling equipment. Therefore, the total time to complete the first, second, and third phases of repair work is: , recorded as ;in, This is the maximum value in the total time of the third phase of repair work, that is, the longest time required to repair the entire station building's HVAC system, power supply, and main signaling equipment; It is the sum of the maximum value of the total repair time in the first phase, the maximum value of the total repair time in the second phase, and the maximum value of the total repair time in the third phase. That is, the longest time required to repair the entire station building structure and vertical transportation components, enclosure structure and water supply and drainage system, as well as HVAC system and power and main signal equipment.

[0062] G5. Finally, calculate the total time required to complete the repair work in phases 1, 2, 3, and 4. The details are as follows:

[0063] Calculate the number of workers required to repair the signal equipment at level k, 0.01A(k), and the total working hours. The repair time is obtained by dividing the total working hours by the number of workers required. Therefore, the total time to complete the repair work in stages 1, 2, 3, and 4 is... , recorded as ;in, This is the maximum value in the total time of the fourth phase of repair work, that is, the longest time required to repair the secondary signal equipment of the entire station building; It is the sum of the maximum values ​​of the total repair time in the first phase, the second phase, the third phase, and the fourth phase, which is the maximum time required to repair the entire station building structure and vertical transportation components, the enclosure structure and water supply and drainage system, the HVAC system and power and main signal equipment, and the secondary signal equipment.

[0064] G6, Return to the repair time of each stage , , , .

[0065] Furthermore, the specific steps for post-earthquake functional calculation in step G are as follows:

[0066] Read the number of components or equipment on each floor Read the corresponding expected values ​​from the vulnerability parameter library of various components or equipment. Then the functional loss reduction coefficient of the p-th type of component or equipment under the d-th damage state is: , recorded as ;in Let be the expected value of the vulnerability parameter of the p-th type of component or equipment under the D-th damage state, where D is the damage state index corresponding to the failure of the component or equipment, which is related to the type of component or equipment; Let be the expected value of the p-th type of component or equipment under the d-th damage state;

[0067] Read the number of damaged components or equipment Quantity of components or equipment on each floor Calculate the damage probability of components or equipment under each damage state. , recorded as Then calculate the building function at the time of the earthquake. , recorded as ;in, This represents the number of the p-th type of component or equipment in the k-th layer under the d-th damage state. The quantity of the p-th type of component or equipment in the k-th layer; The functional loss reduction coefficient for the p-th type of component or equipment under the d-th damage state; The probability of damage to the p-th type of component or equipment on the k-th floor under the d-th damage state at the time after the earthquake; The number of components or equipment in the k-th layer;

[0068] copy Update the damage probabilities corresponding to structural components and vertical transportation components. If the damage state is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then calculate the building function after the first repair phase is completed. , recorded as ;in, The damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state after repairing structural components and vertical transportation components;

[0069] copy Update the damage probabilities corresponding to the building envelope and water supply and drainage system. If the damage state is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then, calculate the building functions after the completion of the first and second repair phases. , recorded as ;in, For complex structural components and vertical transportation components, enclosure structures and water supply and drainage systems, the damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state is given.

[0070] copy Update the damage probability corresponding to the HVAC system, power supply, and main signal equipment. If the damage status is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as P. r3 Then, the building functions after the completion of the first, second, and third repair phases are calculated. , recorded as ;in, The damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state after repairing structural components and vertical transportation components, enclosure structures and water supply and drainage systems, HVAC systems and power and main signal equipment.

[0071] Copy P r3 Update the damage probability corresponding to the secondary signal device. If the damage state is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then, the building functions after the completion of the first, second, third, and fourth repair phases are calculated. , recorded as ;in, The damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state after repairing structural components and vertical transportation components, enclosure structures and water supply and drainage systems, HVAC systems and power and main signal equipment, and secondary signal equipment.

[0072] Return to post-earthquake function .

[0073] Furthermore, based on the obtained resilience assessment indices, a post-earthquake functional recovery function model that varies with repair time is fitted and established, and a post-earthquake functional recovery curve is plotted. The seismic toughness index is calculated based on the area of ​​the post-earthquake functional recovery curve obtained by integration, and the seismic toughness assessment level is classified according to the seismic toughness index, specifically including:

[0074] Read repair time and post-earthquake functions and determine the work stoppage time. ;

[0075] Define a nonlinear post-earthquake functional recovery function model that considers the impact of downtime. ,as follows:

[0076] ;

[0077] Based on the established functional recovery function model, the post-earthquake functional recovery curve is plotted, the area of ​​the post-earthquake functional recovery curve is obtained by integration, and the seismic toughness index R is obtained by normalization.

[0078] Based on the seismic toughness index R, the seismic toughness assessment level of railway passenger station buildings is determined. If the evaluation level is excellent; If the assessment level is good, then the assessment level is good; if If the assessment level is medium, then the assessment level is medium; if If the assessment level is poor, then the assessment level is poor. If the result is negative, the evaluation level is extremely poor.

[0079] This application provides a method for assessing the seismic toughness of medium and large-sized railway passenger station buildings, which has the following beneficial effects:

[0080] (1) This invention constructs a vulnerability parameter database containing special components of railway passenger stations (steel space frame, steel truss, steel-concrete composite columns, etc.) and professional equipment (such as signaling system, power system, etc.), formulates a post-earthquake repair strategy of "structure first, then enclosure, then electromechanical", constructs a toughness assessment level determination method based on repair time, realizes the quantitative assessment of the seismic toughness level of medium and large railway passenger station buildings, and improves the seismic damage prediction and recovery decision-making ability of the building structure system of newly built or existing large-scale integrated transportation hubs.

[0081] (2) This invention provides a Python-based engineering requirement parameter matrix expansion and verification step, forms a repair time calculation method based on post-earthquake repair strategy, establishes a nonlinear post-earthquake functional recovery model considering the impact of downtime, and a post-earthquake functional calculation method with self-updating damage probability. It clarifies the hierarchical dependence and transmission relationship between post-earthquake station building structural damage and functional loss, accurately characterizes the seismic toughness assessment level of medium and large railway passenger station buildings, and ensures their ability to maintain and quickly restore their use functions after an earthquake. Attached Figure Description

[0082] Figure 1 This is an overall flowchart of a method for assessing the seismic toughness of medium and large-sized railway passenger station buildings, provided in one embodiment of the present invention.

[0083] Figure 2 A flowchart of the expansion and verification of the engineering requirement parameter matrix based on Python in a method for assessing the seismic toughness of medium and large railway passenger station buildings, provided as an embodiment of the present invention;

[0084] Figure 3 A flowchart illustrating the "structure first, then enclosure, then electromechanical" post-earthquake repair strategy in a seismic toughness assessment method for medium and large railway passenger station buildings provided in an embodiment of the present invention.

[0085] Figure 4 The flowchart of the calculation of the Monte Carlo simulation function in a method for assessing the seismic toughness of medium and large railway passenger station buildings provided in an embodiment of the present invention;

[0086] Figure 5 A flowchart illustrating the repair time calculation in a seismic toughness assessment method for medium and large-sized railway passenger station buildings, provided in one embodiment of the present invention;

[0087] Figure 6 A flowchart of post-earthquake functional calculation for a method for assessing the seismic toughness of medium and large railway passenger station buildings, provided as an embodiment of the present invention;

[0088] Figure 7 The diagram shows the functional recovery curve of a railway passenger station in a method for assessing the seismic toughness of medium and large railway passenger station buildings, as provided in one embodiment of the present invention. Detailed Implementation

[0089] The present invention will now be described in further detail with reference to specific embodiments and accompanying drawings. Similar elements in different embodiments are referred to by associated similar element reference numerals. In the following embodiments, many details are described to facilitate a better understanding of this application. However, those skilled in the art will readily recognize that some features may be omitted in different situations, or may be replaced by other elements, materials, or methods. In some cases, certain operations related to this application are not shown or described in the specification. This is to avoid obscuring the core parts of this application with excessive description. For those skilled in the art, detailed description of these related operations is not necessary; they can fully understand the related operations based on the description in the specification and general technical knowledge in the art.

[0090] Furthermore, the features, operations, or characteristics described in the specification can be combined in any suitable manner to form various embodiments. At the same time, the steps or actions in the method description can be rearranged or adjusted in a manner obvious to those skilled in the art. Therefore, the various orders in the specification and drawings are only for the clear description of a particular embodiment and do not imply a necessary order, unless otherwise stated that a particular order must be followed.

[0091] This invention provides a method for assessing the seismic toughness of medium-to-large railway passenger station buildings. It employs an algorithm based on the FEMA-P58 engineering demand parameter matrix expansion algorithm and a toughness index calculation algorithm based on Monte Carlo sampling simulation. Furthermore, it establishes a post-earthquake repair strategy of "structure first, then enclosure, and finally electromechanical systems," forming a toughness level determination and functional recovery simulation based on repair time. Ultimately, this method achieves a quantitative assessment of the seismic toughness of medium-to-large railway passenger station buildings, improving the feasibility and accuracy of seismic toughness assessment for railway passenger stations. The following section combines... Figure 1 Please provide a detailed explanation.

[0092] Step 1: Based on the finite element analysis of railway passenger station buildings under seismic conditions, obtain the engineering demand parameter matrix of railway passenger station buildings, and obtain the expanded engineering demand parameter matrix by expanding and verifying the engineering demand parameter matrix of railway passenger station buildings.

[0093] The engineering requirements parameters include inter-story drift angle and peak floor acceleration. A finite element model of the station structure was established, seismic conditions were set, and elastoplastic time history analysis was performed. After calculation, the inter-story drift angle and peak floor acceleration of each floor under each seismic condition were extracted and used as engineering requirements parameters to assemble an engineering requirements parameter matrix. , where m is the row of the matrix, each row corresponds to a seismic condition; n is the column of the matrix, each column corresponds to an engineering requirement parameter, which includes the inter-story drift angle and peak floor acceleration of each floor.

[0094] Among them, the steps of expanding and verifying the engineering requirement parameter matrix (such as...) Figure 2 The specifics are as follows:

[0095] Step 1.1: Read the original project requirement parameter matrix ;

[0096] Step 1.2: For the matrix Taking the logarithm, we obtain the transformation matrix. The covariance matrix is ​​calculated using the `cov` and `mean` functions of Python NumPy array manipulation. and mean matrix ;

[0097] Step 1.3: Calculate the covariance matrix using the `linalg.matrix_rank` and `linalg.eigh` functions of Python NumPy. rank R, eigenvalues eigenvectors ;

[0098] Step 1.4: Partition the eigenvalue matrix according to rank R eigenvector matrix ,have to and Then, the eigenvalue matrix after segmentation is processed using the sqrt and diag functions of Python NumPy array operations. Square root and transform into a diagonal matrix ;

[0099] Step 1.5: Determine the expansion count x, and generate an array with dimension x using the random.randn function of Python NumPy array operations. Independent standard normal random variable matrix And based on the mean matrix of the original engineering requirement parameter matrix Use the dot and ones functions to generate a dimension of matrix ;

[0100] Step 1.6: Using the dot and exp functions of Python NumPy array manipulation, calculate the expanded engineering requirement parameter matrix according to the following formula. And calculate in step 1.2. covariance matrix and mean matrix ;

[0101] (1)

[0102] Step 1.7: Calculate the absolute error of the covariance before and after expansion using the abs function of Python NumPy array manipulation. absolute error of expectation ;

[0103] Step 1.8: Determine the original covariance matrix through screening. and The mask of non-zero elements;

[0104] Step 1.9: Create the covariance matrix using the zeros_like functions of Python NumPy array manipulation. and Given a zero matrix of the same dimension, calculate the relative error of the non-zero elements based on the non-zero element mask, set the relative error of other elements to None, and update the corresponding elements in the zero matrix to obtain the relative error matrix of covariance and expectation. , ;

[0105] Step 1.10: Return the covariance and expected relative error , The maximum, minimum, and average values ​​are determined and verified (generally, the average error is used to check whether it exceeds a preset threshold). After verification, the expanded engineering requirement parameter matrix is ​​generated. Save to the specified directory.

[0106] Step 2: Based on the Building Information Model (BIM) of railway passenger station buildings, obtain the quantity C of station building components or equipment. n Station building room area A, project volume U.

[0107] The quantity U is the number of components or equipment counted based on the BIM model, for example:

[0108] The volume of the reinforced concrete frame column is 100 cubic meters; the volume of the steel structure column is 100 tons.

[0109] The ceiling installation area is 100 square meters;

[0110] The number of IFS racks is 10;

[0111] And so on, counting each layer.

[0112] Step 3: Construct a vulnerability parameter database for railway passenger station components or equipment. The vulnerability parameter database contains the expected values ​​and variances of engineering requirement parameters for different components or equipment under different damage states.

[0113] Specifically, the types of damage states are determined based on the component. Different components or equipment have different types of damage states. For example, reinforced concrete frame columns generally have four types of damage states; cable systems generally have two types of damage states. The damage states of these components are determined according to specifications and experimental studies. The expected values ​​and standard deviations of the engineering requirement parameters corresponding to different components or equipment under different damage states are also determined according to specifications and experimental studies.

[0114] The components or equipment include the following:

[0115] (1) Structural components: reinforced concrete frame columns, reinforced concrete frame beams, reinforced concrete shear walls, reinforced concrete connecting beams, steel structure beams, steel structure columns, steel bracing components, steel-concrete composite columns, steel-concrete composite beams, steel space frames, and steel trusses;

[0116] (2) Displacement-sensitive non-structural components: infill walls, glass curtain walls, stairs;

[0117] (3) Acceleration-sensitive non-structural components: ceiling, elevator, suspended light fixture, switchgear, distribution box, water supply pipe, fire sprinkler pipe, sprinkler head riser, HVAC duct, air outlet, HVAC fan, air conditioning system fan;

[0118] (4) Specialized equipment: battery cabinet, cable system, IFS cabinet, IFSI cabinet, ISFS cabinet, IFSSI cabinet, IRC cabinet, ISRC cabinet, train control cabinet, CTC cabinet, interlocking cabinet, RBC cabinet, TSRS cabinet.

[0119] Step Four: Develop a post-earthquake repair strategy of "structure first, then enclosure, then electromechanical systems," which involves repairing structural components and vertical transportation components, enclosure structures and water supply and drainage systems, HVAC systems and electrical and main signaling equipment, and secondary signaling equipment (such as...) in four phases. Figure 3 ).

[0120] In step four, the structural components and vertical transportation components include the following: reinforced concrete frame columns, reinforced concrete frame beams, reinforced concrete shear walls, reinforced concrete connecting beams, steel structure beams, steel structure columns, steel support components, steel-concrete composite columns, steel-concrete composite beams, steel space frames, steel trusses, stairs, and elevators.

[0121] In step four, the building envelope and water supply and drainage system include the following: infill walls, glass curtain walls, suspended ceilings, water supply pipes, fire sprinkler pipes, and sprinkler head risers;

[0122] In step four, the HVAC system and power and main signal equipment include the following: suspended lighting fixtures, switchgear, distribution boxes, air outlets, HVAC ducts, HVAC fans, air conditioning system fans, battery cabinets, cable systems, ISFS cabinets, ISFSI cabinets, ISRC cabinets, train control cabinets, CTC cabinets, interlocking cabinets, RBC cabinets, and TSRS cabinets.

[0123] In step four, the secondary signal equipment includes the following: IFS cabinet, IFSI cabinet, and IRC cabinet.

[0124] Step 5: Based on the established post-earthquake repair strategy, and according to the expanded engineering requirement parameter matrix, building information model data, and vulnerability parameter database, call the Monte Carlo simulation function to calculate the resilience assessment indicators of the station, including repair time, repair cost, casualties, and post-earthquake function, and calculate the fitted values ​​of each resilience assessment indicator according to the specified confidence level.

[0125] Among them, the specific steps of the Monte Carlo simulation function (such as...) Figure 4 )as follows:

[0126] Step 5.1: Define the initial number of damaged components or equipment. (Initially set to 0, used to store the number of damaged components or equipment on each floor), and initialize the floor indices for inter-story drift angle and peak floor acceleration. ;

[0127] Step 5.2: Traverse the extended matrix of project requirement parameters row by row. , got the first Line (operating condition), number The value of the column (parameter), i.e. .

[0128] Step 5.3: If If the inter-story drift angle is used, then the number of components or equipment in that story is determined by iterative filtering. The term "structural component" or "displacement-sensitive non-structural component" is used; if To determine the peak floor acceleration, the number of components or equipment on each floor is then iterated through and filtered. The "acceleration-sensitive non-structural components" or "specialized equipment" mentioned above;

[0129] Step 5.4: If the number of components or devices in the selected layer is not zero, then iterate through the damage states of the component or device and calculate the log-normal cumulative distribution probability of the p-th type of component or device in the d-th damage state using the norm.cdf function of the Python Scipy library. As shown in the following formula;

[0130] (2)

[0131] in, Let f be the expected value of the p-th component or equipment in the fragility parameter database F under the d-th damage state; Let be the variance of the p-th component or equipment in the fragility parameter database F under the d-th damage state.

[0132] Step 5.5: Determine the quantity of the p-th type of component or equipment in this layer. Repeated random sampling Each time, a random number between 0 and 1 is generated. Based on the probability interval of the random number falling into each range, determine the number of components or equipment of type p in each damage state.

[0133] The probability intervals are divided as follows:

[0134] (1) If the probability If it contains 4 elements, it means that the damage state of the component is of type 4. Indicates no damage ; Indicates minor injury ; Indicates moderate damage ; Indicates severe injury ; Indicates complete destruction .

[0135] (2) If If it contains 3 elements, it indicates that the damage state of the component is of type 3. Indicates no damage ; Indicates minor injury ; Indicates moderate damage ; Indicates severe injury .

[0136] (3) If If it contains two elements, it indicates that the damage state of the component is of type 2. Indicates no damage ; Indicates minor injury ; Indicates severe injury .

[0137] (4) If If a component contains one element, it indicates that the damage state is type 1. Indicates no damage ; Indicates severe injury .

[0138] Step 5.6: If For inter-story drift angle, update The Middle The number of components or equipment of type p in the d-th damage state is determined, and the floor index is updated. Repeat steps 5.2 to 5.5 above; if Update for peak floor acceleration The Middle The p-th type of component or equipment in the layer Count the number of components under each damage condition and update the floor index. Repeat steps 5.2 to 5.5 above;

[0139] Step 5.7: Calculate the repair time, repair cost, casualties, and post-earthquake functionality of the station building based on the sampling results, and return the calculation results.

[0140] Among them, the specific steps for calculating the repair time (such as...) Figure 5 )as follows:

[0141] Step (1): Count the number of damaged components or equipment. Determine the repair time Q for components or equipment, the number of workers q required to repair a single component or piece of equipment, the station building area A, and the repair cost reduction factor. The station building room area A was obtained from the BIM model, and other parameters were determined according to the "Evaluation Standard for Seismic Toughness of Buildings" (GB / T 38591-2020).

[0142] Step (2): Calculate the number of workers required to repair the vertical transportation components on the k-th floor. and total working hours The repair time can be obtained by dividing the total working hours by the number of workers required. Where p represents the p-th type of component or equipment; P represents the total number of component or equipment types; d represents the d-th damage state; D represents the total number of damage state types; k represents the k-th floor; and K represents the total number of floors. The number of workers required to repair component or equipment of type p; This represents the number of the p-th type of component or equipment in the k-th layer under the d-th damage state. Repair time for component or equipment of type p under damage condition d.

[0143] Step (3): Calculate the maximum number of workers allowed on the k-th floor (0.026A(k)), the number of workers required to repair the structural components (0.02A(k)), and the total working hours for repairing the structural components. Provided that the sum of the number of workers needed to repair structural components and the number of workers needed to repair vertical transportation components does not exceed the maximum number of workers allowed per floor, the time required to repair the structural components can be obtained by dividing the total working hours by the number of workers needed. ;

[0144] Step (4): Extraction and The larger of the two, serving as the total time for repairing structural components and vertical transportation components. Therefore, the total time to complete the first phase of repair work is... Recorded as ;in The total time for repairing structural components and vertical transportation components on the kth floor; The maximum value in the total time for repairing the structural components and vertical transportation components of each floor (Phase 1), that is, the longest time required to repair the entire station building's structural components and vertical transportation components;

[0145] Step (5): Calculate the number of workers required for repairing the k-th floor enclosure structure and water supply and drainage system, 0.01A(k), and the total working hours. By dividing the total working hours by the number of workers required, the repair time T2(k) can be obtained, and thus the total time to complete the repair work in stages 1 and 2 is calculated. , recorded as ;in, The maximum value in the total time for repairing the building envelope and water supply and drainage system of each floor (Phase 2), that is, the longest time required to repair the entire station building envelope and water supply and drainage system; The sum of the maximum total time required to repair the structural components and vertical transportation components of each floor and the maximum total time required to repair the enclosure structure and water supply and drainage system of each floor is the maximum time required to repair the entire station building's structural components, vertical transportation components, enclosure structure, and water supply and drainage system.

[0146] Step (6): Calculate the number of workers required to repair the power and main signal equipment on the k-th floor. and total working hours The repair time can be obtained by dividing the total working hours by the number of workers required. ;

[0147] Step (7): Calculate the maximum number of workers allowed to repair the k-th floor (0.026A(k)), the number of workers required to repair the HVAC system (0.01A(k)), and the total man-hours for repairing the HVAC system. Provided that the sum of the number of workers needed to repair the HVAC system and the number of workers needed to repair the electrical and main signaling equipment does not exceed the maximum number of workers allowed on the floor, the time required to repair the HVAC system can be obtained by dividing the total working hours by the number of workers needed. ;

[0148] Step (8): Extraction and The larger of the two, representing the total time required to repair the HVAC system and electrical and major signaling equipment. Therefore, the total time to complete the first, second, and third phases of repair work is: , recorded as ;in, The maximum value in the total time required to repair the HVAC system, power and main signaling equipment on each floor (Phase 3), i.e. the longest time required to repair the entire station building's HVAC system, power and main signaling equipment; The sum of the maximum total time required to repair the structural components and vertical transportation components of each floor, the maximum total time required to repair the enclosure structure and water supply and drainage system of each floor, and the maximum total time required to repair the heating, ventilation and air conditioning system and power and main signaling equipment of each floor, is the maximum time required to repair the entire station building's structural components and vertical transportation components, enclosure structure and water supply and drainage system, and heating, ventilation and air conditioning system and power and main signaling equipment.

[0149] Step (9): Calculate the number of workers required to repair the signal equipment at level k, 0.01A(k), and the total working hours. The repair time can be obtained by dividing the total working hours by the number of workers required. Therefore, the total time to complete the repair work in stages 1, 2, 3, and 4 is... , recorded as ;in, The maximum value in the total time required to repair the secondary signal equipment on each floor (Phase 4) is the maximum time required to repair the secondary signal equipment in the entire station building. The maximum time required to repair the structural components and vertical transportation components of each floor, the maximum time required to repair the enclosure structure and water supply and drainage system of each floor, the maximum time required to repair the heating, ventilation and air conditioning system, power and main signaling equipment of each floor, and the maximum time required to repair the main signaling equipment of each floor is the sum of the four values. In other words, it is the maximum time required to repair the entire station building's structural components and vertical transportation components, enclosure structure and water supply and drainage system, heating, ventilation and air conditioning system, power and main signaling equipment, and secondary signaling equipment.

[0150] Step (10): Return the repair time for each stage .

[0151] The specific steps for calculating the repair costs are as follows:

[0152] Step (1): Count the number of damaged components or equipment. Quantity of components or equipment on each floor Quantity U, determine the component or equipment loss coefficient. Cost and expenses C, repair coefficient Repair cost reduction factor Floor location influence coefficient ;

[0153] Step (2): Based on the number of damaged components or equipment and the number of components or equipment on each floor The component damage ratio was obtained. Then, calculate the repair cost of the p-th component in the k-th layer under the d-th damage condition. , recorded as ;in, For the quantity of components or equipment of type p; The cost of the p-th component in the k-th layer under the d-th damage condition; The component damage ratio of the p-th type of component or equipment in the k-th layer under the d-th damage state; Let be the loss coefficient of the p-th component or equipment under the d-th damage state;

[0154] Step (3): Calculate the repair cost of the p-th type of component in the k-th layer. , recorded as ;in, Let be the repair coefficient of the p-th component or equipment under the d-th damage state;

[0155] Step (4): Calculate the total repair cost of all components on the k-th floor. , recorded as ;

[0156] Step (5): Calculate the total repair cost , and return; where K is the total number of floors; is the influence coefficient of the floor location on the kth floor.

[0157] The specific steps for calculating casualties are as follows:

[0158] Step (1): Count the number of damaged components or equipment. Quantity of components or equipment on each floor The number of structural components and non-structural components that could cause injury or death (including infill walls, glass curtain walls, and suspended ceilings) should be counted separately. , Determine indoor occupancy density Station building room area A;

[0159] Step (2): Calculate the proportion of structural components and non-structural components that may cause injury or death on each floor under different damage states. This allows for the determination of the damage level of each floor, and the calculation of the nominal casualty rate based on the damage level. and nominal mortality rate ;in, The number of structural components of type p in the k-th layer under type d damage state; Let p be the number of non-structural components that can cause injury or death in the k-th layer;

[0160] Specifically, the proportion of structural components and non-structural components that could cause injury or death on each floor under different damage states was analyzed. By consulting the "Standard for Seismic Toughness Evaluation of Buildings," the floor damage level can be determined. Based on the floor damage level, the nominal casualty rate can then be determined. and nominal mortality rate ;

[0161] Step (3): Calculate the population density of the k-th layer. , recorded as Where i represents the room type, including: waiting room, shop, restaurant, office, equipment room, and restroom; Let i be the area of ​​the i-th type of room on the k-th floor, where i is the room type. The indoor occupancy density of type i housing was determined through a survey; Let be the area of ​​the k-th type of house;

[0162] Step (4): Calculate the number of injured. and death toll , respectively denoted as and ;in, Let K be the population density at the kth layer. The nominal mortality rate for the k-th layer; The nominal casualty rate for the k-th layer;

[0163] Step (5): Calculate the true casualty rate and the actual mortality rate , and return.

[0164] Among them, the specific steps of post-earthquake functional calculation (such as...) Figure 6 )as follows:

[0165] Step (1): Read the quantity of components or equipment on each floor Read the corresponding expected values ​​from the vulnerability parameter library of various components or equipment. Then the functional loss reduction coefficient of the p-th type of component or equipment under the d-th damage state is: Let α be the denoted α; where Let be the expected value of the vulnerability parameter of the p-th type of component or equipment under the D-th damage state, where D is the damage state index corresponding to the failure of the component or equipment, which is related to the type of component or equipment; Let be the expected value of the p-th type of component or equipment under the d-th damage state;

[0166] Step (2): Count the number of damaged components or equipment. Quantity of components or equipment on each floor Calculate the damage probability of components or equipment under each damage state. , recorded as Then calculate the building function at the time of the earthquake. , denoted as F s0 ;in, This represents the number of the p-th type of component or equipment in the k-th layer under the d-th damage state. The quantity of the p-th type of component or equipment in the k-th layer; The functional loss reduction coefficient for the p-th type of component or equipment under the d-th damage state; The probability of damage to the p-th type of component or equipment on the k-th floor under the d-th damage state at the time after the earthquake; The number of components or equipment in the k-th layer;

[0167] Step (3): Copy Update the damage probabilities corresponding to structural components and vertical transportation components. If the damage state is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then calculate the building function after the first repair phase is completed. , denoted as F s1 ;in, The damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state after repairing structural components and vertical transportation components;

[0168] Step (4): Copy P r1 Update the damage probabilities corresponding to the building envelope and water supply and drainage system. If the damage state is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then, calculate the building functions after the completion of the first and second repair phases. , recorded as ;in, For complex structural components and vertical transportation components, enclosure structures and water supply and drainage systems, the damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state is given.

[0169] Step (5): Copy Update the damage probability corresponding to the HVAC system, power supply, and main signal equipment. If the damage status is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then, the building functions after the completion of the first, second, and third repair phases are calculated. , recorded as ;in, The damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state after repairing structural components and vertical transportation components, enclosure structures and water supply and drainage systems, HVAC systems and power and main signal equipment.

[0170] Step (6): Copy Update the damage probability corresponding to the secondary signal device. If the damage state is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then, the building functions after the completion of the first, second, third, and fourth repair stages are calculated. , recorded as ;in, The damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state after repairing structural components and vertical transportation components, enclosure structures and water supply and drainage systems, HVAC systems and power and main signal equipment, and secondary signal equipment.

[0171] Step (7): Return to post-earthquake functionality .

[0172] The specific steps for calculating the fitted value are as follows:

[0173] Step (1): Initialize three lists to store the fitted values. Expected fit value , standard deviation of fitted values And read the repair time, repair cost, casualties, and post-earthquake functional data from the Monte Carlo simulation results;

[0174] Step (2): Define a function that can calculate the fitted value of the given data at a specified confidence level using the stats.lognorm.fit and stats.lognorm.ppf functions in the Python Scipy library;

[0175] Step (3): Based on the data from step (1), call the fitting function from step (2) to calculate the fitting values, expected values, and standard deviations of the fitting values ​​for repair time, repair cost, casualties, post-earthquake function, and save them to the corresponding lists.

[0176] Step 6: Based on the obtained resilience assessment index, fit and establish a post-earthquake functional recovery function model that varies with repair time, and plot the post-earthquake functional recovery curve. Calculate the seismic resilience index based on the area of ​​the post-earthquake functional recovery curve obtained by integration, and classify the seismic resilience assessment level according to the seismic resilience index.

[0177] The specific steps of the toughness grade determination method based on repair time are as follows:

[0178] Step 6.1: Read the repair time and post-earthquake functions Determine the work stoppage time (Customizable based on actual needs);

[0179] Step 6.2: Define a nonlinear function recovery model that considers the impact of downtime. ,as follows:

[0180] (3)

[0181] Step 6.3: Based on the repair time and functional recovery model, generate time series and post-earthquake functional series, plot the post-earthquake functional recovery curve, integrate to obtain the area of ​​the functional recovery curve, and obtain the post-earthquake toughness index R through normalization.

[0182] Step 6.4: Determine the seismic toughness assessment level of the railway passenger station based on the post-earthquake toughness index R. If it is excellent; if If it is good; if If it is, then it is considered medium; if If it is poor, then it is considered poor; if If the range is 0, then it is the extreme range.

[0183] Application examples are as follows:

[0184] Applying the technology of this invention, the original engineering demand parameter matrix (39 rows × 16 columns) of a railway passenger station was read. The matrix was expanded 1000 times, taking approximately 5-10 seconds. The expansion results showed that the maximum relative error of the variance was 2.08%, the minimum relative error was 0.36%, and the average relative error was 1.41%. The expected maximum relative error was 0.79%, the minimum relative error was 0.00%, and the average relative error was 0.21%. Further reading of the expanded engineering demand parameter matrix and building information model data, the Monte Carlo simulation was performed 1000 times, with the city size set to a large city and a fitting confidence level of 0.84. The time to complete 1000 Monte Carlo simulations was approximately 15-20 minutes. The simulation results showed that the fitted value for the first stage repair time was 60.01 days (expected 6.69, variance 2.21), and the fitted value for the second stage repair time was 133.36 days (expected 67.08, variance 0.69). The fitted value for the repair time in Phase 3 was 133.69 days (expected 67.73 days, variance 0.68), and the fitted value for the repair time in Phase 4 was 133.73 days (expected 68.33 days, variance 0.68). The fitted value for repair cost was 11.027 million yuan (expected 7.2804 million yuan, variance 0.42 million yuan). The fitted value for the casualty rate was 0.25% (expected 0.03%, variance 2.31%), and the fitted value for the mortality rate was 0.19% (expected 0.09%, variance 0.77%). Using the resilience index data obtained from the Monte Carlo simulation, and setting the downtime to 0 days, a recovery curve for the railway passenger station's functions was plotted. Figure 7 As shown in the figure. The evaluation results indicate that the seismic toughness index of the station building is 0.753, and the seismic toughness level is medium (three stars). This invention takes into account the unique building components and professional equipment of railway passenger stations, and establishes a post-earthquake repair strategy of "structure first, then enclosure, and then electromechanical". It realizes the seismic toughness assessment process of medium and large railway passenger station buildings in a clear and intuitive way, and improves the convenience and feasibility of seismic toughness assessment of railway passenger stations.

[0185] The above examples illustrate the present invention only to aid in understanding it and are not intended to limit the scope of the invention. Those skilled in the art can make various simple deductions, modifications, or substitutions based on the principles of this invention.

Claims

1. A method for assessing the seismic toughness of medium and large-sized railway passenger station buildings, characterized in that, The method includes: Based on the finite element analysis of railway passenger station buildings under seismic conditions, the engineering demand parameter matrix of railway passenger station buildings is obtained. By expanding and verifying the engineering demand parameter matrix of railway passenger station buildings, the expanded engineering demand parameter matrix is ​​obtained. Based on the Building Information Model (BIM) of railway passenger station buildings, obtain railway passenger station building information model data including the number of station building components or equipment, station building room area, and engineering quantity. A vulnerability parameter database for railway passenger station building components or equipment is constructed. The vulnerability parameter database contains the expected values ​​and variances of engineering requirement parameters for different components or equipment under different damage states. The post-earthquake repair strategy of "structure first, then enclosure, then electromechanical" was formulated, which means repairing the structural components and vertical transportation components, the enclosure structure and water supply and drainage system, the heating, ventilation and air conditioning system and power and main signal equipment, and the secondary signal equipment in four stages in sequence. Based on the established post-earthquake repair strategy, and according to the expanded engineering demand parameter matrix, building information model data, and vulnerability parameter database, the Monte Carlo simulation function was used to calculate the resilience assessment indicators of railway passenger station buildings, including repair time and post-earthquake function. Based on the obtained resilience assessment index, a post-earthquake functional recovery function model that varies with repair time is fitted and established, and a post-earthquake functional recovery curve is plotted. The seismic resilience index is calculated based on the area of ​​the post-earthquake functional recovery curve obtained by integration, and the seismic resilience assessment level is divided according to the seismic resilience index. in: Based on finite element analysis of railway passenger station buildings under seismic conditions, the engineering requirement parameter matrix of railway passenger station buildings is obtained, specifically including: A finite element model of the railway passenger station structure was established, and seismic conditions were set to perform elastoplastic time history analysis on the railway passenger structure. After the calculation, the inter-story drift angle and peak floor acceleration of each floor under each seismic condition were extracted, and a parameter matrix of engineering requirements for the railway passenger station was constructed. Where m is the row of the matrix, each row corresponds to a seismic condition; n is the column of the matrix, each column corresponds to an engineering requirement parameter, which includes the inter-story drift angle and peak floor acceleration of each floor. The engineering demand parameter matrix for railway passenger station buildings was expanded and verified to obtain an expanded engineering demand parameter matrix, which specifically includes: Read the original project requirement parameter matrix ; Calculate the original engineering requirement parameter matrix covariance matrix and mean matrix ; Calculate the covariance matrix rank R, eigenvalue matrix eigenvector matrix ; The eigenvalue matrix is ​​divided according to rank R. eigenvector matrix ,get and Further segmentation of the eigenvalue matrix Square root and transform into a diagonal matrix ; Determine the number of expansions x, and generate a dimension of Independent standard normal random variable matrix And based on the mean matrix Generate a dimension of Mean expansion matrix ; Calculate the expanded engineering requirement parameter matrix using the following formula. : 。 2. The method for assessing the seismic toughness of medium and large-sized railway passenger station buildings as described in claim 1, characterized in that, Based on the Building Information Model (BIM) of railway passenger station buildings, data on the BIM of railway passenger station buildings, including the quantity of station building components or equipment, the area of ​​station building rooms, and the amount of work, are obtained. Specifically, this includes: The station building components or equipment include: Structural components include: reinforced concrete frame columns, reinforced concrete frame beams, reinforced concrete shear walls, reinforced concrete connecting beams, steel structure beams, steel structure columns, steel bracing components, steel-concrete composite columns, steel-concrete composite beams, steel space frames, and steel trusses. Displacement-sensitive non-structural components include: infill walls, glass curtain walls, and staircases; Acceleration-sensitive non-structural components include: ceilings, elevators, suspended lighting fixtures, switchgear, distribution boxes, water supply pipes, fire sprinkler pipes, sprinkler head risers, HVAC ducts, air outlets, HVAC fans, and air conditioning system fans; Specialized equipment, including: battery cabinets, cable systems, IFS cabinets, IFSI cabinets, ISFS cabinets, IFSSI cabinets, IRC cabinets, ISRC cabinets, train control cabinets, CTC cabinets, interlocking cabinets, RBC cabinets, and TSRS cabinets.

3. The method for assessing the seismic toughness of medium and large-sized railway passenger station buildings as described in claim 2, characterized in that, A post-earthquake repair strategy of "structure first, then enclosure, then electromechanical" was adopted, which means repairing structural components and vertical transportation components, enclosure structure and water supply and drainage system, HVAC system and electrical and main signaling equipment, and secondary signaling equipment in four stages in sequence. Specifically, this includes: The structural components and vertical transportation components include: reinforced concrete frame columns, reinforced concrete frame beams, reinforced concrete shear walls, reinforced concrete connecting beams, steel structure beams, steel structure columns, steel support components, steel-concrete composite columns, steel-concrete composite beams, steel space frames, steel trusses, stairs, and elevators. The enclosure structure and water supply and drainage system include: infill walls, glass curtain walls, suspended ceilings, water supply pipes, fire sprinkler pipes, and sprinkler head risers; The HVAC system and power and main signaling equipment include: suspended lighting fixtures, switchgear, distribution boxes, air outlets, HVAC ducts, HVAC fans, air conditioning system fans, battery cabinets, cable systems, ISFS cabinets, ISFSI cabinets, ISRC cabinets, train control cabinets, CTC cabinets, interlocking cabinets, RBC cabinets, and TSRS cabinets. The secondary signaling equipment includes: IFS cabinet, IFSI cabinet, and IRC cabinet.

4. The method for assessing the seismic toughness of medium and large-sized railway passenger station buildings as described in claim 3, characterized in that, Based on the established post-earthquake repair strategy, and using the expanded engineering demand parameter matrix, building information model data, and vulnerability parameter database, Monte Carlo simulation functions were used to calculate resilience assessment indicators for railway passenger station buildings, including repair time and post-earthquake functionality. Fitted values ​​for each indicator were calculated at specified confidence levels, specifically including: Step A: Define the number of damaged components or equipment. A floor index used to store the number of damaged components or equipment on each floor, initialized to 0, and used to initialize the inter-floor displacement angle and peak floor acceleration. , All are 0; Step B: Traverse the expanded engineering requirement parameter matrix row by row. This yields the value in the i-th row and j-th column, i.e. ; Step C, if If the inter-floor displacement angle is used, then the number of components or equipment on the current floor is traversed and filtered. The term "structural member" or "displacement-sensitive non-structural member" is used; if To determine the peak floor acceleration, the number of components or equipment on each floor is then iterated through and filtered. The "acceleration-sensitive non-structural components" or "specialized equipment" mentioned above; Step D: If the number of components or equipment on the current floor is not zero after filtering, then iterate through the damage states of the corresponding components or equipment, and calculate the log-normal cumulative distribution probability of the p-th type of component or equipment under the d-th damage state using the norm.cdf function of the Python Scipy library. As shown in the following formula: ； in, Let f be the expected value of the p-th component or equipment in the fragility parameter database F under the d-th damage state; Let be the variance of the p-th component or equipment in the fragility parameter database F under the d-th damage state; Step E: Determine the quantity of the p-th type of component or equipment in the current layer. Repeated random sampling Each time, a random number r between 0 and 1 is generated, and the number of components or equipment of type p in each damage state is determined according to the probability interval of the random number falling into it. Step F, if If the inter-story drift angle is [value], then update [the algorithm / mechanism]. The Middle Calculate the number of components or equipment of type p in the d-th damage state, and update the floor index. Repeat steps B through E; if For the peak floor acceleration, update The Middle Calculate the number of components or equipment of type p in the d-th damage state, and update the floor index. Repeat steps B through E; Step G: Calculate the repair time, repair cost, casualties, and post-earthquake functionality of railway passenger station buildings based on the sampling results, and return the calculation results.

5. The method for assessing the seismic toughness of medium and large-sized railway passenger station buildings as described in claim 4, characterized in that, In step E, the probability interval is divided as follows: If the probability If it contains 4 elements, it means that the damage state of the component is of type 4. Indicates no damage ; Indicates minor injury ; Indicates moderate damage ; Indicates severe injury ; Indicates complete destruction ; if If it contains 3 elements, it indicates that the damage state of the component is of type 3. Indicates no damage ; Indicates minor injury ; Indicates moderate damage ; Indicates severe injury ; if If it contains two elements, it indicates that the damage state of the component is of type 2. Indicates no damage ; Indicates minor injury ; Indicates severe injury ; if If a component contains one element, it indicates that the damage state is type 1. Indicates no damage ; Indicates severe injury .

6. The method for assessing the seismic toughness of medium and large-sized railway passenger station buildings as described in claim 5, characterized in that, In step G, the specific steps for calculating the repair time are as follows: G1, Read the number of damaged components or equipment Determine the repair time Q for components or equipment, the number of workers q required to repair a single component or piece of equipment, the station building area A, and the repair cost reduction factor. ; G2. First, calculate the total time T1 for completing the first phase of repair work, as follows: First, calculate the number of workers required to repair the vertical transportation components on the k-th floor. and total working hours The repair time is obtained by dividing the total working hours by the number of workers required. Where p represents the p-th type of component or equipment; P represents the total number of component or equipment types; d represents the d-th damage state; D represents the total number of damage state types; k represents the k-th floor; and K represents the total number of floors. The number of workers required to repair component or equipment of type p; This represents the number of the p-th type of component or equipment in the k-th layer under the d-th damage state. Repair time for component or equipment of type p under damage condition d; Next, calculate the maximum allowed number of workers on the k-th floor (0.026A(k)), the required number of workers to repair structural components (0.02A(k)), and the total working hours for repairing structural components. Provided that the sum of the number of workers needed to repair structural components and the number of workers needed to repair vertical transportation components does not exceed the maximum number of workers allowed per floor, the time required to repair the structural components is obtained by dividing the total working hours by the number of workers needed. ; Compare and extract and The larger of the two, serving as the total time for repairing structural components and vertical transportation components. Therefore, the total time to complete the first phase of repair work is... , recorded as ;in The total time for repairing structural components and vertical transportation components on the kth floor; This is the maximum value in the total time of the first phase of repair work, that is, the longest time required to repair the entire station building structure and vertical transportation components; G3. Next, calculate the total time to complete the first and second phases of repair work. The details are as follows: Calculate the number of workers required (0.01A(k)) and the total man-hours for repairing the k-th floor enclosure structure and water supply and drainage system. The repair time is obtained by dividing the total working hours by the number of workers required. Therefore, the total time to complete the first and second phases of repair work is... , recorded as ;in, This is the maximum value in the total time of the second phase of repair work, that is, the longest time required to repair the entire station building envelope and water supply and drainage system; It is the sum of the maximum value of the total time of the first phase of repair work and the maximum value of the total time of the second phase of repair work, that is, the longest time required to repair the entire station building structure components, vertical transportation components, enclosure structure and water supply and drainage system. G4. Recalculate the total time required to complete the repair work in phases 1, 2, and 3. The details are as follows: First, calculate the number of workers required to repair the power and main signal equipment on the k-th floor. and total working hours The repair time T is obtained by dividing the total working hours by the number of workers required. 3_eq (k); Next, calculate the maximum number of workers allowed to repair the k-th floor (0.026A(k)), the number of workers required to repair the HVAC system (0.01A(k)), and the total man-hours for repairing the HVAC system. Provided that the sum of the number of workers needed to repair the HVAC system and the number of workers needed to repair the electrical and main signaling equipment does not exceed the maximum number of workers allowed on the floor, the repair time for the HVAC system can be obtained by dividing the total working hours by the number of workers needed. ; Comparison Extraction and The larger of the two, representing the total time required to repair the HVAC system and electrical and major signaling equipment. Therefore, the total time to complete the first, second, and third phases of repair work is: , recorded as ;in, This is the maximum value in the total time of the third phase of repair work, that is, the longest time required to repair the entire station building's HVAC system, power supply, and main signaling equipment; It is the sum of the maximum value of the total repair time in the first phase, the maximum value of the total repair time in the second phase, and the maximum value of the total repair time in the third phase. That is, the longest time required to repair the entire station building structure and vertical transportation components, enclosure structure and water supply and drainage system, as well as HVAC system and power and main signal equipment. G5. Finally, calculate the total time T4 for completing the repair work in stages 1, 2, 3, and 4, as follows: Calculate the number of workers required to repair the signal equipment at level k, 0.01A(k), and the total working hours. The repair time is obtained by dividing the total working hours by the number of workers required. Therefore, the total time to complete the repair work in stages 1, 2, 3, and 4 is... , recorded as ;in, This is the maximum value in the total time of the fourth phase of repair work, that is, the longest time required to repair the secondary signal equipment of the entire station building; It is the sum of the maximum values ​​of the total repair time in the first phase, the second phase, the third phase, and the fourth phase, which is the maximum time required to repair the entire station building structure and vertical transportation components, the enclosure structure and water supply and drainage system, the HVAC system and power and main signal equipment, and the secondary signal equipment. G6, Return to the repair time of each stage , , , .

7. The method for assessing the seismic toughness of medium and large-sized railway passenger station buildings according to claim 6, characterized in that, In step G, the specific steps for post-earthquake functional calculation are as follows: Read the number of components or equipment on each floor Read the corresponding expected values ​​from the vulnerability parameter library of various components or equipment. Then the functional loss reduction coefficient of the p-th type of component or equipment under the d-th damage state is: Let α be the denoted α; where Let be the expected value of the vulnerability parameter of the p-th type of component or equipment under the D-th damage state, where D is the damage state index corresponding to the failure of the component or equipment, which is related to the type of component or equipment; Let be the expected value of the p-th type of component or equipment under the d-th damage state; Read the number of damaged components or equipment Quantity of components or equipment on each floor Calculate the damage probability of components or equipment under each damage state. , recorded as Then calculate the building function at the time of the earthquake. , recorded as ;in, This represents the number of the p-th type of component or equipment in the k-th layer under the d-th damage state. The quantity of the p-th type of component or equipment in the k-th layer; The functional loss reduction coefficient for the p-th type of component or equipment under the d-th damage state; The probability of damage to the p-th type of component or equipment on the k-th floor under the d-th damage state at the time after the earthquake; The number of components or equipment in the k-th layer; Copy P r0 Update the damage probabilities corresponding to structural components and vertical transportation components. If the damage state is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then calculate the building function after the first repair phase is completed. , denoted as F s1 ;in, The damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state after repairing structural components and vertical transportation components; Copy P r1 Update the damage probabilities corresponding to the building envelope and water supply and drainage system. If the damage state is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then, calculate the building functions after the completion of the first and second repair phases. , recorded as ;in, For complex structural components and vertical transportation components, enclosure structures and water supply and drainage systems, the damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state is given. copy Update the damage probability corresponding to the HVAC system, power supply, and main signal equipment. If the damage status is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then, the building functions after the completion of the first, second, and third repair phases are calculated. , recorded as ;in, The damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state after repairing structural components and vertical transportation components, enclosure structures and water supply and drainage systems, HVAC systems and power and main signal equipment. copy Update the damage probability corresponding to the secondary signal device. If the damage state is... If the probability of updating damage is 1, then the probability of updating damage is 0, and the new list is denoted as . Then, the building functions after the completion of the first, second, third, and fourth repair phases are calculated. , recorded as ;in, The damage probability of the p-th type of component or equipment on the k-th floor under the d-th damage state after repairing structural components and vertical transportation components, enclosure structures and water supply and drainage systems, HVAC systems and power and main signal equipment, and secondary signal equipment. Return to post-earthquake function , , , , .

8. The method for assessing the seismic toughness of medium and large-sized railway passenger station buildings according to claim 7, characterized in that, Based on the obtained resilience assessment indices, a post-earthquake functional recovery function model that varies with repair time is fitted and established, and a post-earthquake functional recovery curve is plotted. The seismic resilience index is calculated based on the area of ​​the post-earthquake functional recovery curve obtained by integration, and the seismic resilience assessment level is classified according to the seismic resilience index, specifically including: Read repair time and post-earthquake functions and determine the work stoppage time. ; Define a nonlinear post-earthquake functional recovery function model that considers the impact of downtime. ,as follows: ； Based on the established functional recovery function model, the post-earthquake functional recovery curve is plotted, the area of ​​the post-earthquake functional recovery curve is obtained by integration, and the seismic toughness index R is obtained by normalization. Based on the seismic toughness index R, the seismic toughness assessment level of railway passenger station buildings is determined. If the evaluation level is excellent; If the assessment level is good, then the assessment level is good; if If so, the evaluation level is medium; If the assessment level is poor, then the assessment level is poor. If the result is negative, the evaluation level is extremely poor.

Images