Method, system, device and medium for evaluating dynamic stability of high-cold rock landslide

Abstract

The application provides a high-cold rock landslide dynamic stability evaluation method, system, device and medium, comprising: obtaining a global large-scale temperature field and stress field boundary efficiently by using a macroscopic slope model, then accurately mapping it to a mesoscopic discrete fracture network model in a potential damage area, adopting a "global rough calculation + local precise calculation" nesting strategy, which avoids the huge calculation overhead brought by full model discretization, and ensures the high precision of fracture evolution calculation in the key disaster-causing area. The high-cold rock landslide dynamic stability evaluation method can realize dynamic space-time tracking of the rock mass deterioration process under freeze-thaw cycles, and further realize dynamic prediction of the slope life and long-term stability.

Description

Technical Field

[0001] This application belongs to the interdisciplinary field of geological disaster prevention and computational geotechnical mechanics, and in particular relates to a method, system, equipment and medium for assessing the dynamic stability of cold-resistant rock landslides. Background Technology

[0002] Traditional thermo-hydraulic-mechanical (THM) coupled models treat rock masses as continuous media, failing to reflect the local stress concentrations and discontinuous displacements caused by ice crystal expansion within joint / fracture networks. Furthermore, current safety factor calculations do not consider the non-uniform distribution of ice pressure on potential slip surfaces and its weakening effect on anti-slip forces, resulting in a significant lag in predicting instability of high-altitude, cold-climate landslides. Summary of the Invention

[0003] In view of this, this application aims to propose a method, system, equipment and medium for assessing the dynamic stability of cold-resistant rock landslides, in order to solve at least one of the above-mentioned problems.

[0004] To achieve the above objectives, the technical solution of this application is implemented as follows:

[0005] Firstly, this application provides a method for assessing the dynamic stability of high-altitude, cold-weather rock landslides, including:

[0006] Acquire slope geometry data, rock mass property data and multi-period climate data of the target cold region slope, establish a macro slope model based on the slope data, and perform transient heat conduction analysis on the macro slope model to obtain the macro temperature field distribution and thermal stress evolution process of the slope under freeze-thaw cycle.

[0007] Based on the obtained macroscopic temperature field and stress field, local key areas inside the slope that are severely affected by freeze-thaw action and have potential instability are identified, and dynamic boundary conditions at the boundaries of the local key areas are extracted.

[0008] Within the specified key local area, a microscopic discrete fracture network model is constructed based on the actual joint statistical distribution characteristics of the target rock mass; wherein, the microscopic discrete fracture network model consists of complete discrete rock blocks and contact surfaces that cut between the blocks;

[0009] The dynamic boundary conditions are mapped to the microscopic discrete fracture network model to perform microscopic freeze-thaw fluid-structure interaction analysis, the generated dynamic ice shearing force is applied to the boundary of adjacent discrete rock blocks, and the slope instability of the local rock mass is judged.

[0010] In response to local rock mass slope instability, the microscopic failure surface characteristics formed by the convergence and expansion of microcracks are extracted and recorded. Combined with the macroscopic slope model, the overall slope instability analysis and failure surface identification are carried out to output the dynamic stability assessment results of the high-altitude cold rock slope.

[0011] Secondly, based on the same inventive concept, this application also provides a dynamic stability assessment system for high-altitude cold-resistant rock landslides, including:

[0012] The macroscopic model construction module is configured to acquire slope geometry data, rock mass property data and multi-period climate data of the target cold region slope, establish a macroscopic slope model based on the slope data, and perform transient heat conduction analysis on the macroscopic slope model to obtain the macroscopic temperature field distribution and thermal stress evolution process of the slope under freeze-thaw cycle.

[0013] The boundary condition extraction module is configured to identify local key areas inside the slope that are severely affected by freeze-thaw action and have potential instability based on the obtained macroscopic temperature field and stress field, and extract the dynamic boundary conditions at the boundary of the local key area.

[0014] The microscopic model construction module is configured to construct a microscopic discrete fracture network model based on the actual joint statistical distribution characteristics of the target rock mass within the local key area; wherein, the microscopic discrete fracture network model is composed of complete discrete rock blocks and contact surfaces that cut each other between the blocks;

[0015] The fluid-structure interaction analysis module is configured to map the dynamic boundary conditions to the microscopic discrete fracture network model to perform microscopic freeze-thaw fluid-structure interaction analysis, apply the generated dynamic ice shearing force to the boundary of adjacent discrete rock blocks, and determine the slope instability of the local rock mass.

[0016] The instability analysis module is configured to extract and record the microscopic failure surface features formed by the convergence and expansion of microcracks in response to local rock mass slope instability, and combine it with the macroscopic slope model to perform overall slope instability analysis and failure surface identification, so as to output the dynamic stability assessment results of the high-altitude cold rock slope.

[0017] Thirdly, based on the same inventive concept, this application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the method described in the first aspect.

[0018] Fourthly, based on the same inventive concept, this application also provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer instructions for causing the computer to perform the method as described in the first aspect.

[0019] Compared with existing technologies, the dynamic stability assessment method, system, equipment, and medium for high-altitude cold-weather rock landslides described in this application have the following advantages:

[0020] The dynamic stability assessment method for high-altitude cold-resistant rock landslides described in this application utilizes a macroscopic slope model to efficiently obtain the global large-scale temperature and stress field boundaries, and then accurately maps them to a microscopic discrete fracture network model of the potential failure zone. This nested strategy of "global coarse calculation + local fine calculation" avoids the huge computational overhead caused by the discretization of the entire model, and ensures extremely high accuracy in the calculation of fracture evolution in key disaster-causing areas. At the same time, it can realize the dynamic spatiotemporal tracking of the rock mass deterioration process under freeze-thaw cycles, thereby realizing the dynamic prediction of slope life and long-term stability. Attached Figure Description

[0021] The accompanying drawings, which form part of this application, are used to provide a further understanding of this application. The illustrative embodiments and descriptions of this application are used to explain this application and do not constitute an undue limitation of this application. In the drawings:

[0022] Figure 1 This is a flowchart of a dynamic stability assessment method for high-altitude rocky landslides as described in an embodiment of this application.

[0023] Figure 2 This is a schematic diagram illustrating the nested relationship between the macroscopic slope model and the microscopic discrete fracture network model described in the embodiments of this application;

[0024] Figure 3 This is a schematic diagram illustrating the mechanism of water-ice phase transition and dynamic ice shearing force within microscopic discrete fissures as described in the embodiments of this application.

[0025] Figure 4 This is a schematic diagram of the structure of a dynamic stability assessment system for high-altitude rock landslides as described in an embodiment of this application;

[0026] Figure 5 This is a schematic diagram of the hardware structure of the electronic device described in an embodiment of this application. Detailed Implementation

[0027] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with specific embodiments and the accompanying drawings.

[0028] It should be noted that, unless otherwise defined, the technical or scientific terms used in the embodiments of this application should have the ordinary meaning understood by one of ordinary skill in the art to which this application pertains. The terms "first," "second," and similar terms used in the embodiments of this application do not indicate any order, quantity, or importance, but are merely used to distinguish different components. Terms such as "comprising" or "including" mean that the element or object preceding the word encompasses the elements or objects listed after the word and their equivalents, without excluding other elements or objects. Terms such as "connected" or "linked" are not limited to physical or mechanical connections, but can include electrical connections, whether direct or indirect. Terms such as "upper," "lower," "left," and "right" are only used to indicate relative positional relationships; when the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0029] The method described in this embodiment can overcome the limitations of traditional continuous medium models and realistically depict the microscopic ice shearing disaster mechanism. Existing technologies typically simplify freeze-thaw cycles to equivalent temperature stress or parameter reduction in a macroscopic continuous medium, failing to reflect the true mechanical effects of water-ice phase transition on rock fractures. This application establishes a microscopic discrete fracture network model to directly simulate the water-ice phase transition process and the volume expansion of ice crystals at the microfracture level, deriving and applying realistic dynamic ice shearing forces. This method fundamentally reveals the disaster-causing mechanical origin of "fracture water freezing - volume expansion - ice shearing expansion" within high-altitude cold rock masses, significantly improving the physical realism of the mechanical mechanism.

[0030] The embodiments of this application are described in detail below with reference to the accompanying drawings.

[0031] Please see Figure 1 As shown in the figure, this embodiment provides a method for assessing the dynamic stability of high-altitude rock landslides, which specifically includes the following steps:

[0032] Step S101: Obtain slope geometry data, rock mass property data, and multi-period climate data of the target cold region slope. Establish a macro slope model based on the slope data and perform transient heat conduction analysis on the macro slope model to obtain the macro temperature field distribution and thermal stress evolution process of the slope under freeze-thaw cycle.

[0033] Specifically, in this embodiment, the main purpose of this step is to establish a macroscopic slope model and perform heat conduction analysis. First, basic input data for the target cold-region slope is acquired, including slope geometry data, rock mass property data (including rock mechanics and thermal material properties), and multi-period climate data (such as atmospheric temperature fluctuation sequences over time). A macroscopic slope model (i.e., a macroscopic geometric model) is established based on finite element software and the geometry data. Material property data is assigned to different material partitions of the geometric model. Then, multi-field coupled finite element calculations are performed on the model based on the multi-period climate data. Transient heat conduction analysis is performed on the macroscopic slope model, using finite element software to simulate the macroscopic temperature and stress field distribution of the slope under freeze-thaw cycles over a certain period. The next step proceeds after the finite element software simulation converges. If it does not converge, the calculation parameters are adjusted and the solution is recalculated.

[0034] Step S102: Based on the obtained macroscopic temperature field and stress field, identify the local key areas inside the slope that are severely affected by freeze-thaw action and have potential instability, and extract the dynamic boundary conditions at the boundary of the local key areas.

[0035] Specifically, in this embodiment, the main purpose of this step is to extract the boundary conditions and cross-scale data transfer of local key regions. Based on the temperature and stress fields calculated in step S101, local key regions (i.e., sub-regions) that are severely affected by freeze-thaw cycles and have potential instability risks are identified, and dynamic temperature history data at the boundaries of these sub-regions are extracted. And stress boundary conditions, serving as a "data bridge" connecting macroscopic and mesoscopic models, are illustrated in the diagram below, showing the nested relationship between the macroscopic slope model and the mesoscopic discrete fracture network model. Figure 2 As shown.

[0036] Step S103: In the local key area, construct a micro-discrete fracture network model based on the actual joint statistical distribution characteristics of the target rock mass; wherein, the micro-discrete fracture network model is composed of complete discrete rock blocks and contact surfaces that cut each other between the blocks.

[0037] Specifically, in this embodiment, the main purpose of this step is to establish a micro-discrete fracture network model. For the local key areas identified in step S102, within these areas, a micro-discrete fracture network (DFN) model is constructed, which includes complete discrete blocks and a joint network, using discrete element methods (such as UDEC) and based on the actual joint statistical distribution characteristics of the target rock mass. This model consists of complete discrete rock blocks (such as Voronoi polygon blocks) and the contact surfaces (joints) that cut between the blocks.

[0038] This embodiment innovatively proposes a "macro-micro" cross-scale nested mechanism, balancing computational efficiency and local accuracy. Addressing the challenge of computational overload caused by directly using discrete element method (DEM) calculations for large-scale slopes, this embodiment proposes a cross-scale coupled computational architecture. It efficiently obtains the global large-scale temperature and stress field boundaries using a macroscopic continuous medium model (such as the finite element method), and then accurately maps them to a microscopic discrete fracture network model (such as the DEM) of the potential damage zone. This "global coarse calculation + local fine calculation" nested strategy avoids the enormous computational overhead of full model discretization while ensuring extremely high accuracy in calculating fracture evolution in key disaster-causing areas.

[0039] Step S104: Map the dynamic boundary conditions to the microscopic discrete fracture network model to perform microscopic freeze-thaw fluid-structure interaction analysis, apply the generated dynamic ice shearing force to the boundary of adjacent discrete rock blocks, and judge the slope instability of the local rock mass.

[0040] Specifically, in this embodiment, the dynamic boundary conditions extracted in step S102 are mapped to a mesoscopic discrete fracture network model (a schematic diagram of the water-ice phase transition and dynamic ice wedging force mechanism within the mesoscopic discrete fracture network model is shown below). Figure 3 As shown), identify the water-filled fracture network in the model and calculate the ice-water phase transition and ice wedging force based on the following physical mechanisms:

[0041] (1) Dynamic evolution of unfrozen water content: The volume of unfrozen water inside the fissure is controlled by negative temperature, and the calculation model is as follows:

[0042] ;

[0043] In the formula, For temperature The unfrozen water content below saturated moisture content This is the attenuation parameter.

[0044] (2) Calculation of frozen water volume increment: Since water expands by about 9% when it freezes into ice, the frozen water volume increment is calculated as follows: The calculation formula is:

[0045] ;

[0046] In the formula, The initial water-bearing volume of the fracture. , The densities of liquid water and solid ice are respectively. This represents the increase in the volume of frozen water.

[0047] (3) Generate and apply dynamic ice wedging force (i.e., effective ice wedging force): when the temperature inside the fissure... At temperatures below 0°C, a phase transition is triggered, generating effective ice-weeping force. The following formula is used to calculate and apply the ice shearing force to the boundary of adjacent discrete rock blocks, where ice shearing force is... The formula is:

[0048] ;

[0049] In the formula, This represents the initial hydraulic aperture of the fracture. Joint normal stiffness; ∈[0.6,0.9] represents the ice pressure transmission efficiency coefficient (which can be calibrated through indoor single-crack frost heave tests); This is the upper limit of the tensile strength of the rock mass, used to prevent the occurrence of non-physical infinite expansion forces in the program calculation.

[0050] (4) Joint parameter update: Along with the generation of ice weft force, the joint stiffness (normal stiffness) of the fracture is updated synchronously. Tangential stiffness and hydraulic opening .

[0051] Step S105: In response to the local rock mass slope instability, the microscopic failure surface characteristics formed by the convergence and expansion of microcracks are extracted and recorded, and the overall slope instability analysis and failure surface identification are carried out in combination with the macroscopic slope model to output the dynamic stability assessment results of the high-altitude cold rock slope.

[0052] Specifically, in this embodiment, the main purpose of this step is to perform instability analysis and failure surface identification, and to output stability assessment results. Under the combined action of ice shearing force and external load, the geotechnical engineering numerical analysis software (UDEC, Universal Distinct Element Code, mainly used to simulate the mechanical behavior of discontinuous media) is used to perform mechanical solutions to determine whether the local rock mass (here referring to discrete rock blocks) has undergone relative slippage, and whether the joints have produced normal tensile yielding or tangential shear failure (i.e., slope instability judgment).

[0053] If no instability occurs (No), time advances to the next freeze-thaw cycle, returns to perform the heat conduction and fluid-structure interaction analysis for the next freeze-thaw cycle, updates the boundary temperature, and accumulates the damage.

[0054] If local instability occurs and intersects with a crack (Yes), the propagation path of the microcrack is recorded, and the microscopic damage surface features are identified and extracted.

[0055] This embodiment enables dynamic spatiotemporal tracking of rock mass degradation processes under freeze-thaw cycles. Traditional stability assessments are mostly static or single-failure analyses, which are difficult to reflect the unique periodic degradation characteristics of cold regions. This embodiment binds the time step to the depth of the freeze-thaw cycle, dynamically updating not only the boundary temperature in each cycle but also simultaneously updating micro-parameters such as joint stiffness (normal stiffness and tangential stiffness) and hydraulic aperture. This mechanism allows this application to continuously track the cumulative damage process of slopes under multi-year freeze-thaw cycles, achieving dynamic prediction of slope life and long-term stability.

[0056] Based on the comprehensive macroscopic and microscopic stress states and the geometric characteristics of the failure surface, the dynamic safety factor of the slope considering the ice shearing force effect is calculated. :

[0057] ;

[0058] In the formula, and The first The cohesion and internal friction angle of each sliding surface segment; The area of ​​the base; This represents the normal resistance of the bottom surface of the slider. It is a tangential sliding force; This is the dynamic ice shearing force acting on the slip surface. This formula quantitatively reveals the ice shearing force ( As an unloading effect, it effectively offsets the positive resistance of the slip surface. This leads to a disastrous mechanical mechanism that causes a sharp decline in shear strength. The final output determines the stability of the slope after multiple freeze-thaw cycles.

[0059] This embodiment can accurately reproduce the progressive failure and true sliding morphology of complex rock mass instability. Existing macroscopic limit equilibrium methods or strength reduction methods often only assume a single circular or polygonal sliding surface. The microscopic discrete fracture network model of this application allows rock blocks to undergo realistic sliding, separation, and rotation, and can intuitively simulate the entire process of microfracture initiation, joint expansion and convergence, and eventual connection to form a macroscopic failure surface. This method breaks free from the constraints of a pre-set sliding surface and can accurately capture the progressive failure morphology of high-altitude cold rock slopes, from local ice-fracturing to overall collapse and instability.

[0060] It should be noted that the above description describes some embodiments of this application. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in a different order than that shown in the above embodiments and still achieve the desired result. Furthermore, the processes depicted in the drawings do not necessarily require a specific or sequential order to achieve the desired result. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0061] Based on the same inventive concept, and corresponding to any of the above embodiments, the embodiments of this application also provide a dynamic stability assessment system for high-altitude cold rock landslides.

[0062] like Figure 4 As shown, the dynamic stability assessment system for high-altitude cold-weather rock landslides includes:

[0063] The macro model construction module 11 is configured to acquire slope geometry data, rock mass property data and multi-period climate data of the target cold region slope, establish a macro slope model based on the slope data, and perform transient heat conduction analysis on the macro slope model to obtain the macro temperature field distribution and thermal stress evolution process of the slope under freeze-thaw cycle.

[0064] The boundary condition extraction module 12 is configured to identify local key areas inside the slope that are severely affected by freeze-thaw action and have potential instability based on the obtained macroscopic temperature field and stress field, and extract the dynamic boundary conditions at the boundary of the local key area.

[0065] The microscopic model construction module 13 is configured to construct a microscopic discrete fracture network model based on the statistical distribution characteristics of the actual joints of the target rock mass within a local key region; wherein, the microscopic discrete fracture network model consists of complete discrete rock blocks and contact surfaces that cut each other between the blocks;

[0066] The fluid-structure interaction analysis module 14 is configured to map dynamic boundary conditions to a micro-discrete fracture network model to perform micro-freeze-thaw fluid-structure interaction analysis, apply the generated dynamic ice shearing force to the boundary of adjacent discrete rock blocks, and determine the slope instability of local rock masses.

[0067] The instability analysis module 15 is configured to extract and record the microscopic failure surface features formed by the convergence and expansion of microcracks in response to local rock mass slope instability, and combine the macroscopic slope model to perform overall slope instability analysis and failure surface identification, so as to output the dynamic stability assessment results of the high-altitude cold rock slope.

[0068] For ease of description, the above system is described by dividing it into various modules based on their functions. Of course, in implementing the embodiments of this application, the functions of each module can be implemented in one or more software and / or hardware.

[0069] The system described in the above embodiments is used to implement the corresponding method in any of the foregoing embodiments and has the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0070] Based on the same inventive concept, corresponding to the methods of any of the above embodiments, embodiments of this application also provide an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the methods described in any of the above embodiments.

[0071] Figure 5 This embodiment illustrates a more specific hardware structure of an electronic device, which may include a processor 1010, a memory 1020, an input / output interface 1030, a communication interface 1040, and a bus 1050. The processor 1010, memory 1020, input / output interface 1030, and communication interface 1040 are interconnected internally via the bus 1050.

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

[0073] The memory 1020 can be implemented in the form of ROM (Read Only Memory), RAM (Random Access Memory), static storage device, dynamic storage device, etc. The memory 1020 can store the operating system and other applications. When the technical solutions provided in the embodiments of this specification are implemented by software or firmware, the relevant program code is stored in the memory 1020 and is called and executed by the processor 1010.

[0074] The input / output interface 1030 is used to connect input / output modules to realize information input and output. The input / output modules can be configured as components in the device (not shown in the figure) or externally connected to the device to provide corresponding functions. Input devices may include keyboards, mice, touch screens, microphones, various sensors, etc., and output devices may include displays, speakers, vibrators, indicator lights, etc.

[0075] The communication interface 1040 is used to connect a communication module (not shown in the figure) to enable communication between this device and other devices. The communication module can communicate via wired means (such as USB, Ethernet cable, etc.) or wireless means (such as mobile network, WIFI, Bluetooth, etc.).

[0076] Bus 1050 includes a pathway for transmitting information between various components of the device, such as processor 1010, memory 1020, input / output interface 1030, and communication interface 1040.

[0077] It should be noted that although the above-described device only shows the processor 1010, memory 1020, input / output interface 1030, communication interface 1040, and bus 1050, in specific implementations, the device may also include other components necessary for normal operation. Furthermore, those skilled in the art will understand that the above-described device may only include the components necessary for implementing the embodiments of this specification, and not necessarily all the components shown in the figures.

[0078] The electronic devices described above are used to implement the corresponding methods in any of the foregoing embodiments and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0079] Based on the same inventive concept, corresponding to the methods of any of the above embodiments, this application also provides a non-transitory computer-readable storage medium that stores computer instructions for causing the computer to perform the methods described in any of the above embodiments.

[0080] The computer-readable medium of this embodiment includes permanent and non-permanent, removable and non-removable media, and information storage can be implemented by any method or technology. Information can be computer-readable instructions, data structures, program modules, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transfer medium that can be used to store information accessible by a computing device.

[0081] The computer instructions stored in the storage medium of the above embodiments are used to cause the computer to perform the methods described in any of the above embodiments, and have the beneficial effects of the corresponding method embodiments, which will not be repeated here.

[0082] Those skilled in the art should understand that the discussion of any of the above embodiments is merely exemplary and is not intended to imply that the scope of this application (including the claims) is limited to these examples; within the framework of this application, the technical features of the above embodiments or different embodiments can also be combined, the steps can be implemented in any order, and there are many other variations of different aspects of the embodiments of this application as described above, which are not provided in the details for the sake of brevity.

[0083] Although this application has been described in conjunction with specific embodiments thereof, many substitutions, modifications, and variations of these embodiments will be apparent to those skilled in the art from the foregoing description. For example, other memory architectures (e.g., dynamic RAM (DRAM)) may be used with the embodiments discussed.

[0084] The embodiments of this application are intended to cover all such substitutions, modifications, and variations that fall within the broad scope of the appended claims. Therefore, any omissions, modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the embodiments of this application should be included within the protection scope of this application.

Claims

1. A method for evaluating dynamic stability of high-cold rock landslide, characterized in that, include: Acquire slope geometry data, rock mass property data and multi-period climate data of the target cold region slope, establish a macro slope model based on the slope data, and perform transient heat conduction analysis on the macro slope model to obtain the macro temperature field distribution and thermal stress evolution process of the slope under freeze-thaw cycle. Based on the obtained macroscopic temperature field and stress field, local key areas inside the slope that are severely affected by freeze-thaw action and have potential instability are identified, and dynamic boundary conditions at the boundaries of the local key areas are extracted. Within the specified key local area, a microscopic discrete fracture network model is constructed based on the actual joint statistical distribution characteristics of the target rock mass; wherein, the microscopic discrete fracture network model consists of complete discrete rock blocks and contact surfaces that cut between the blocks; The dynamic boundary conditions are mapped to the microscopic discrete fracture network model to perform microscopic freeze-thaw fluid-structure interaction analysis, the generated dynamic ice shearing force is applied to the boundary of adjacent discrete rock blocks, and the slope instability of the local rock mass is judged. In response to local rock mass slope instability, the microscopic failure surface characteristics formed by the convergence and expansion of microcracks are extracted and recorded. Combined with the macroscopic slope model, the overall slope instability analysis and failure surface identification are carried out to output the dynamic stability assessment results of the high-altitude cold rock slope.

2. The method according to claim 1, characterized in that: The dynamic boundary conditions include time-temperature history data and stress boundary conditions, serving as a data bridge between the macroscopic slope model and the microscopic discrete fracture network model.

3. The method of claim 1, wherein, The micro-scale freeze-thaw fluid-structure interaction analysis includes: The model identifies a network of fissures filled with water, determines the water-ice phase transition process and the volume increment of frozen water within the fissures based on the temperature field, and generates dynamic ice-cutting force based on the volume increment of frozen water.

4. The method according to claim 3, characterized in that: The volume of unfrozen water inside the fissure is controlled by negative temperature. The dynamic evolution model of the unfrozen water content is as follows: ； wherein is the temperature is the unfrozen water content at the temperature is the saturation moisture content is the decay parameter The formula for the volume increase of frozen water is: ； wherein is the initial water volume in the fracture, , are the densities of liquid water and solid ice, respectively, is the incremental volume of frozen water. When the temperature inside the fissure At temperatures below 0°C, a phase transition is triggered to generate dynamic ice sheathing force, wherein the generated dynamic ice sheathing force... It is expressed as follows: ； In the formula, This represents the initial hydraulic aperture of the fracture. Joint normal stiffness; The ice pressure transmission efficiency coefficient; This represents the upper limit of the tensile strength of the rock mass.

5. The method according to claim 3, characterized in that: The generated dynamic ice shearing force is applied to the boundary of adjacent discrete rock blocks, and the joint stiffness and hydraulic aperture parameters of the fracture are updated simultaneously.

6. The method according to claim 1, characterized in that: Based on a comprehensive analysis of macroscopic and microscopic stress states and the geometric characteristics of the failure surface, a dynamic safety factor for the slope considering the effect of ice shearing force is determined. The dynamic safety factor... The formula is: ； In the formula, and The first The cohesion and internal friction angle of each sliding surface segment; The area of ​​the base; This represents the normal resistance of the bottom surface of the slider. It is a tangential sliding force; This refers to the dynamic ice-cutting force acting on the slip surface.

7. The method according to claim 1, characterized in that, Also includes: If no slope instability occurs in the local rock mass, the time step is advanced, and the heat conduction and fluid-structure interaction analysis for the next freeze-thaw cycle is performed to update the boundary temperature and accumulate the damage.

8. A dynamic stability assessment system for high-altitude cold-weather rock landslides, characterized in that, include: The macroscopic model construction module is configured to acquire slope geometry data, rock mass property data and multi-period climate data of the target cold region slope, establish a macroscopic slope model based on the slope data, and perform transient heat conduction analysis on the macroscopic slope model to obtain the macroscopic temperature field distribution and thermal stress evolution process of the slope under freeze-thaw cycle. The boundary condition extraction module is configured to identify local key areas inside the slope that are severely affected by freeze-thaw action and have potential instability based on the obtained macroscopic temperature field and stress field, and extract the dynamic boundary conditions at the boundary of the local key area. The microscopic model construction module is configured to construct a microscopic discrete fracture network model based on the actual joint statistical distribution characteristics of the target rock mass within the local key area; wherein, the microscopic discrete fracture network model is composed of complete discrete rock blocks and contact surfaces that cut each other between the blocks; The fluid-structure interaction analysis module is configured to map the dynamic boundary conditions to the microscopic discrete fracture network model to perform microscopic freeze-thaw fluid-structure interaction analysis, apply the generated dynamic ice shearing force to the boundary of adjacent discrete rock blocks, and determine the slope instability of the local rock mass. The instability analysis module is configured to extract and record the microscopic failure surface features formed by the convergence and expansion of microcracks in response to local rock mass slope instability, and combine it with the macroscopic slope model to perform overall slope instability analysis and failure surface identification, so as to output the dynamic stability assessment results of the high-altitude cold rock slope.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the method as described in any one of claims 1-7.

10. A non-transitory computer-readable storage medium, characterized in that, in, The non-transitory computer-readable storage medium stores computer instructions for causing a computer to perform the method described in any one of claims 1-7.

Images