A simulation method of non-decaying creep discrete element contact model for frozen soil

Abstract

The present invention discloses a simulation method of a non-decaying creep discrete element contact model of frozen soil, including the steps of establishing a first viscoelastic model and a damage element of frozen soil, establishing a non-decaying creep discrete element contact model of frozen soil according to the first viscoelastic model and the damage element, and performing frozen soil creep simulation. The non-decaying creep discrete element contact model of frozen soil used in the present invention includes a first viscoelastic model, so the non-decaying creep discrete element contact model of frozen soil can simulate frozen soil that has not entered a damaged state; because the non-decaying creep discrete element contact model of frozen soil includes a damage element, and the damage element stiffness of the damage element changes with time, the non-decaying creep discrete element contact model of frozen soil can simulate frozen soil that has entered a damaged state, i.e., the third stage; through the simulation method of the non-decaying creep discrete element contact model of frozen soil of the present invention, a more sophisticated discrete element numerical simulation of frozen soil creep can be performed. The present invention is widely used in the field of frozen soil creep technology.

Description

Technical Field

[0001] The invention relates to the technical field of frozen soil creep, in particular to a simulation method of a frozen soil non-attenuation creep discrete element contact model. Background Art

[0002] Creep phenomenon is widely present in nature and engineering applications, especially in projects under high temperature, high pressure, or low temperature frozen soil environment, such as geological landslides, underground tunnels, earth-rock dams, cold region road foundations, etc. Traditionally, continuum mechanics (such as finite element method) is used to simulate the creep characteristics of geotechnical materials, but the macroscopic continuity assumption is not applicable to multiphase medium materials such as frozen soil. Therefore, the discrete element method is more widely used in the simulation of geotechnical materials. Since the creep of frozen soil has more complex mechanical properties, the current discrete element method simulation technology is not refined enough for the simulation of frozen soil creep. For example, it can only simulate some stages of frozen soil creep, but cannot accurately simulate the non-attenuation creep stage of frozen soil. Summary of the invention

[0003] In view of the technical problems that the current discrete element simulation of frozen soil creep is not refined enough, the purpose of the present invention is to provide a simulation method of a discrete element contact model of frozen soil non-attenuation creep.

[0004] The embodiment of the present invention includes a simulation method of a frozen soil non-attenuation creep discrete element contact model, and the simulation method of the frozen soil non-attenuation creep discrete element contact model includes the following steps:

[0005] Establish the first viscoelastic model of frozen soil;

[0006] Establishing a damaged element; the stiffness of the damaged element is the damaged element stiffness, and the damaged element stiffness changes with time during the initial damage time period;

[0007] A discrete element contact model of frozen soil non-attenuated creep is established according to the first viscoelastic model and the damage element;

[0008] The discrete element numerical simulation of frozen soil creep is carried out according to the frozen soil non-attenuated creep discrete element contact model.

[0009] Furthermore, the establishing of the first viscoelastic model of frozen soil includes:

[0010] Establish the Burgers model of frozen soil;

[0011] The Burgers model is used as the first viscoelastic model.

[0012] Furthermore, the establishment of a frozen soil non-attenuated creep discrete element contact model based on the first viscoelastic model and the damage element includes:

[0013] The first viscoelastic model and the damage element are connected in series to obtain the frozen soil non-attenuated creep discrete element contact model.

[0014] Furthermore, the stiffness of the damaged element is constant in a first time period, and changes with time in a second time period; the first time period is a time period before the second time period.

[0015] Furthermore, the stiffness of the damaged element decreases as time increases during the second time period.

[0016] Further, the first time period is a time period indicating that the damaged element has not entered a damaged state, and the second time period is a time period indicating that the damaged element has entered a damaged state.

[0017] Furthermore, a discrete element numerical simulation of frozen soil creep is performed according to the frozen soil non-attenuated creep discrete element contact model, including:

[0018] The non-attenuated creep discrete element contact model of frozen soil is derived using the central difference method to obtain the contact force-displacement relationship information between frozen soil particles;

[0019] Obtaining a first source code corresponding to the first viscoelastic model;

[0020] According to the contact force-displacement relationship information, secondary development is performed on the first source code to obtain a second source code;

[0021] Compiling the second source code to obtain a dynamic link library file;

[0022] The dynamic link library file is imported into discrete element simulation software for simulation.

[0023] Furthermore, the contact force-displacement relationship information corresponding to the first time period includes:

[0024]

[0025] in, The first time period The contact force of the damaged component at the moment, The first time period The contact force of the damaged component at time +1, The first time period The total displacement of the frozen soil non-attenuated creep discrete element contact model at time, The first time period The total displacement of the frozen soil non-decayed creep discrete element contact model at time +1, The first time period The displacement of the viscous element of the Kelvin model in the Burgers model at time t, , , , , as a parameter.

[0026] Furthermore, the contact force-displacement relationship information corresponding to the second time period includes:

[0027]

[0028] in, The second time period The contact force of the damaged component at the moment, The second time period The contact force of the damaged component at time +1, The second time period The total displacement of the frozen soil non-attenuated creep discrete element contact model at time, The second time period The total displacement of the frozen soil non-decayed creep discrete element contact model at time +1, The second time period The displacement of the viscous element of the Kelvin model in the Burgers model at time t, The second time period The displacement of the damaged element at time , , , , as a parameter.

[0029] Furthermore, the frozen soil creep simulation according to the frozen soil non-attenuated creep discrete element contact model further includes:

[0030] The dynamic link library file is double-ball verified using the simulation software.

[0031] The beneficial effects of the present invention are: the simulation method of the non-decaying creep discrete element contact model of frozen soil in the embodiment, since the used frozen soil non-decaying creep discrete element contact model includes a first viscoelastic model, the non-decaying creep discrete element contact model of frozen soil can simulate frozen soil that has not entered the damaged state, i.e., the first stage and the second stage; since the non-decaying creep discrete element contact model of frozen soil includes damaged elements, and the damaged element stiffness of the damaged elements changes with time, the non-decaying creep discrete element contact model of frozen soil can simulate frozen soil that has entered the damaged state, i.e., the third stage; therefore, through the simulation method of the non-decaying creep discrete element contact model of frozen soil in the embodiment, a more sophisticated discrete element numerical simulation of frozen soil creep can be performed. BRIEF DESCRIPTION OF THE DRAWINGS

[0032] Figure 1 Schematic diagram of the steps of the simulation method of the frozen soil non-attenuated creep discrete element contact model in the embodiment;

[0033] Figure 2 The first time period in the embodiment Schematic diagram of the corresponding frozen soil non-attenuated creep discrete element contact model;

[0034] Figure 3 The second time period in the embodiment Schematic diagram of the corresponding frozen soil non-attenuated creep discrete element contact model;

[0035] Figure 4 and Figure 5 A schematic diagram of debugging and compiling the second source code in an embodiment;

[0036] Figure 6 A schematic diagram of importing a dll file into a simulation software for trial calculation in an embodiment;

[0037] Figure 7 Schematic diagram of double-ball verification in the embodiment. DETAILED DESCRIPTION

[0038] Continuum mechanics such as the finite element method can be used to simulate the creep characteristics of geotechnical materials such as frozen soil. However, when dealing with granular materials, the continuity assumption is not applicable, and the discrete element method (DEM) becomes an effective tool. DEM is based on the interaction between soil particles in geotechnical materials. By calculating the force and displacement between particles, the movement between particles and their overall mechanical behavior are obtained. In order to describe the creep of granular materials, time-dependent viscoelastic elements are usually added to DEM to simulate the creep effect.

[0039] For frozen soil, the Burgers contact model can be used for simulation. The Burgers model is a viscoelastic model composed of multiple springs and viscous elements. Specifically, the Burgers model connects elastic and viscous elements in series and in parallel to form a comprehensive response of the material. It is particularly suitable for simulating material deformation under long-term loads and can well describe the mechanical properties of frozen soil materials in the initial creep and stable creep stages.

[0040] However, in addition to the initial creep and stable creep stages, the creep of frozen soil materials also has a third stage, the non-attenuated creep stage. In the third stage of the frozen soil creep test, the creep rate of the frozen soil shows a nonlinear increasing trend, while the strain rate used in the Burgers model will eventually tend to a constant value. Therefore, the simple Burgers creep contact model cannot reflect the third stage of frozen soil creep. The simulation effect of the simple Burgers creep contact model does not match the actual rheological properties of frozen soil, and is quite different from the mechanical properties of frozen soil, which has obvious limitations.

[0041] Based on the above principle, in this embodiment, a simulation method of a frozen soil non-attenuated creep discrete element contact model is provided. Figure 1 , the simulation method of the discrete element contact model of frozen soil non-decaying creep includes the following steps:

[0042] S1. Establish the first viscoelastic model of frozen soil;

[0043] S2. Establishing damage elements;

[0044] S3. Based on the first viscoelastic model and the damage element, a non-attenuated creep discrete element contact model of frozen soil is established;

[0045] S4. Frozen soil creep simulation is performed based on the frozen soil non-attenuated creep discrete element contact model.

[0046] In this embodiment, each step in the simulation method of the frozen soil non-decaying creep discrete element contact model including S1-S4 can be executed by a computer, wherein the first viscoelastic model, the damage element and the frozen soil non-decaying creep discrete element contact model can be in data form.

[0047] In step S1, the first viscoelastic model established is a model that can simulate the viscosity and elasticity of frozen soil particles. In this embodiment, the Burgers model can be established as the first viscoelastic model.

[0048] In this embodiment, refer to Figure 2 or Figure 3, the first viscoelastic model, namely the Burgers model, is composed of a Maxwell model and a Kelvin model in series. Among them, a Maxwell model is composed of a viscous element (such as a viscoelastic pot) and an elastic element (such as a spring) in series, and a Kelvin model is composed of a viscous element and an elastic element in parallel. Among them, the viscous element has a damping effect, and its parameter can be expressed by viscosity; the elastic element has an elastic effect, and its parameter can be expressed by stiffness. In this embodiment, unless otherwise specified, the viscosity of a viscous element or the stiffness of an elastic element is a constant, that is, it does not change with time.

[0049] In step S2, a damage element is established. In this embodiment, the damage element has the same elastic effect as the elastic element, and its parameters can be expressed by the damage element stiffness. In this embodiment, the stiffness of the damaged element Changes with time during the damage state period, that is, the stiffness of the damaged element It does not maintain a constant at all times, but only maintains a constant for a part of the time period at most, while damaging the stiffness of the component in other time periods. The value of changes with time, and even damages the stiffness of the component over the entire time period The values ​​of change with time, that is, the stiffness of the damaged element The value of has been changing with time.

[0050] In step S3, refer to Figure 2 or Figure 3 , the first viscoelastic model and the damage element can be connected in series to obtain the non-attenuated creep discrete element contact model of frozen soil. Specifically, Figure 2 and Figure 3 The left side of represents a mass The frozen soil particles and a mass The interaction between the frozen soil particles (discrete elements) in the normal direction, and the right side represents a mass The frozen soil particles and a mass The interaction between frozen soil particles (discrete elements) in the shear direction.

[0051] In step S4, use Figure 2 and Figure 3The non-decaying creep discrete element contact model of frozen soil shown can simulate the interaction between particles in the frozen soil material, thereby simulating the non-decaying creep of frozen soil. Specifically, since the non-decaying creep discrete element contact model of frozen soil includes a first viscoelastic model, and the first viscoelastic model can be a conventional model such as the Burgers model, the non-decaying creep discrete element contact model of frozen soil in step S4 can simulate frozen soil that has not entered a damaged state (which can be further subdivided into an initial creep state belonging to the first stage and a stable creep state belonging to the second stage); since the non-decaying creep discrete element contact model of frozen soil includes a damaged element, and the damaged element stiffness of the damaged element changes with time, the non-decaying creep discrete element contact model of frozen soil in step S4 can simulate frozen soil that has entered a damaged state (equivalent to the third stage); therefore, by executing steps S1-S4, a more detailed simulation of frozen soil creep can be performed.

[0052] In this embodiment, the frozen soil non-decaying creep discrete element contact model can be used to simulate the frozen soil creep behavior in the following process: the frozen soil sample is continuously loaded and the timing is started at the same time; in the first time period, In the first time period, the frozen soil samples were not damaged and did not enter the damaged state. The frozen soil samples were in the initial creep and stable creep. The duration is the destruction time of the frozen soil sample. After the end, the loading time of the frozen soil sample exceeds the destruction time of the frozen soil sample, the frozen soil sample is damaged and enters the damaged state. The frozen soil sample is in accelerated creep and the second time period begins. . That is, the first time period It is the time period when the frozen soil sample is in initial creep and stable creep, corresponding to the first and second stages of frozen soil creep; it is the time period when the frozen soil sample is in accelerated creep, corresponding to the third stage of frozen soil creep.

[0053] In order to adapt to the creep behavior of frozen soil in the above process and achieve a more realistic simulation effect, in this embodiment, the damaged element stiffness of the damaged element is Set as the first time period The constant within and the second time period Specifically, the stiffness of the damaged element In the first period is a constant, which enables the frozen soil non-attenuated creep discrete element contact model to simulate the first and second stages of the frozen soil specimen, namely the initial creep and stable creep; the stiffness of the damaged element In the second period The inside is a variable, which can be set to the stiffness of the damaged element In the second period The stiffness of the damaged element decreases with time. It can decrease linearly with time, which enables the frozen soil non-attenuated creep discrete element contact model to simulate the third stage of frozen soil samples, namely accelerated creep.

[0054] In this embodiment, when executing step S4, that is, performing frozen soil creep simulation according to the frozen soil non-attenuated creep discrete element contact model, the following steps may be specifically performed:

[0055] S401. Use the central difference method to derive the non-attenuated creep discrete element contact model of frozen soil and obtain the contact force-displacement relationship information between frozen soil particles;

[0056] S402. Obtaining a first source code corresponding to a first viscoelastic model;

[0057] S403. Perform secondary development on the first source code according to the contact force-displacement relationship information to obtain a second source code;

[0058] S404. Compile the second source code to obtain a dynamic link library file;

[0059] S405. Import the dynamic link library file into the simulation software for simulation.

[0060] As mentioned above, due to the first period is the time period when the frozen soil is not damaged, and the stiffness of the damaged element is is a constant, and the second time period is the time period when frozen soil is damaged, and the stiffness of the damaged element is decreases with time, while the frozen soil non-decaying creep discrete element contact model uses a damaged element stiffness The damage element is established, that is, the frozen soil non-attenuated creep discrete element contact model is established in the first time period and in the second period is different, so when executing step S401, the first time period can be The corresponding frozen soil non-decayed creep discrete element contact model and the second time period The corresponding frozen soil non-attenuated creep discrete element contact model is modeled separately.

[0061] 1. For the first time period The corresponding frozen soil non-attenuated creep discrete element contact model is used for modeling:

[0062] First time period The corresponding frozen soil non-attenuated creep discrete element contact model is as follows: Figure 2 shown.

[0063] for Figure 2As far as the damage element D between particles is concerned, it is considered that the damage element D between particles has not started to be damaged. K D It does not change with time. Therefore, according to the mechanical properties of the spring and combined with the central difference method, the contact force formula of the damaged element at this time can be obtained, and the displacement of the damaged element at time t+1 can be obtained by further calculation. The specific derivation process is as follows:

[0064] (1)

[0065] Where: f is the contact force, K D is the stiffness of the damaged element, u mD is the displacement of the damaged component, u mD t+1 is the displacement of the damaged element at time t+1, u mD t is the displacement of the damaged element at time t, f t+1 is the contact force of the damaged component at time t+1, f t is the contact force of the damaged component at time t.

[0066] Then the total displacement of the overall element model is further calculated:

[0067] (2)

[0068] Further organized form:

[0069] (3)

[0070] For the convenience of abbreviation, let ,

[0071] The contact force expression at time t+1 when damage has not started is obtained:

[0072] (4)

[0073] Formula (4) represents the first time period In the first time period, that is, when the frozen soil is not damaged, the relationship between the contact force and displacement between the frozen soil particles is given by formula (4). Corresponding contact force displacement information.

[0074] 2. For the second time period The corresponding frozen soil non-attenuated creep discrete element contact model is used for modeling:

[0075] Second time period The corresponding frozen soil non-attenuated creep discrete element contact model is as follows: Figure 3 shown.

[0076] At this time, the expressions of component displacement, velocity and acceleration are:

[0077] (5)

[0078] Forces are expressed in Kelvin elements:

[0079] (6)

[0081] For series elements express forces:

[0082] K m spring:

[0083] (7)

[0084] C m Sticky pot:

[0085] (8)

[0086] K D spring:

[0087] (9)

[0088]

[0089] (10)

[0091] The damaged component is calculated by central difference, and the displacement formula of the damaged component is obtained:

[0092] (11)

[0093] make bIna / 2=G , we can get:

[0094] (12)

[0095] make K D =E , bIna△t / 2=F , and the displacement of the damaged component at time t+1 is obtained:

[0096] (13)

[0097] Calculate the total displacement of the component and solve the contact force calculation formula at time t+1:

[0098] (14)

[0099] make ,and

[0100] Simplified:

[0101] (15)

[0102] For the convenience of abbreviation, let , .

[0103] (16)

[0104] (17)

[0105] (18)

[0106] (19)

[0107] The contact force calculation formula at time t+1 is obtained:

[0108] (20)

[0109] To further simplify its form:

[0110] (twenty one)

[0111] The simplified form of the contact force between particles under the damaged state is obtained:

[0112] (twenty two)

[0113] To simplify the displacement calculation at time t+1, we can introduce L :

[0114] (twenty three)

[0115] Formula (22) represents the second time period In the second time period, that is, when the frozen soil is damaged, the relationship between the contact force and displacement between the frozen soil particles is given by formula (22). Corresponding contact force displacement information.

[0116] In this embodiment, although the first time period The intermediate and final results of the modeling of the corresponding frozen soil non-attenuated creep discrete element contact model, namely formulas (1)-(4), are consistent with the second time period. The intermediate and final results of the modeling of the corresponding frozen soil non-attenuated creep discrete element contact model, namely formulas (5)-(23), are used , , , , , , However, we can refer to the principle of local variables in computer programming, where the same symbols can have different meanings. For example, in formulas (1)-(4), , The symbol is for the first time period The corresponding local variables of the frozen soil non-attenuated creep discrete element contact model are effective, among which Indicates the first time period In moment; in formulas (5)-(23), , The symbol is for the second time period The corresponding local variables of the frozen soil non-attenuated creep discrete element contact model are effective, among which Indicates the second time period In time.

[0117] After executing step S401, the contact force displacement information (corresponding to the first time period) shown in formula (4) is obtained. ) and the contact force displacement information shown in formula (22) (corresponding to the second time period ) After that, steps S402-S405 may be executed for secondary development.

[0118] Before executing steps S402-S405, the PFC contact model secondary development compilation environment can be built using PFC software, Visual Studio 2017 and Qt5.12.0 to ensure that the header file and source file of the contact model can be correctly loaded.

[0119] In this embodiment, the first viscoelastic model can be a conventional model such as the Burgers model, so when executing step S402, the source code corresponding to the first viscoelastic model, i.e., the first source code, can be obtained from the technical data. Specifically, the first source code of the first viscoelastic model, i.e., the Burgers model, is found from the help document of PFC, imported into the Visual Studio software, and a header file (.h file) and a source file (.cpp file) are newly created.

[0120] In step S403, in the header file (.h file) of the first source code, the newly introduced variables of the contact force-displacement relationship information obtained in step S401 are declared, and the calculation method of the constants is explained; the new variables are initialized in the source file (.cpp) of the first source code, and the force-displacement law therein is modified to the contact force-displacement relationship information obtained in step S401, so as to realize the secondary development of the first source code and obtain the second source code.

[0121] In step S404, refer to Figure 4 and Figure 5 , debug the second source code obtained in step S403 to obtain a compiled dynamic link library file (.dll).

[0122] In step S405, by writing command flow code, the dll file obtained by secondary development in step S403 is imported into the simulation software. Figure 6 As shown, a trial calculation is performed to verify the simulation effect. In this embodiment, the simulation software used in step S405 may be PFC.

[0123] When executing step S405, double ball verification is performed through PFC, and the result is as follows: Figure 7 As shown. Figure 7 It can be seen that the frozen soil creep contact model after secondary development can better simulate the third stage of frozen soil creep. Figure 7 The third stage of creep, i.e., the accelerated creep failure stage, is obvious in the simulation, which is more in line with the creep mechanical properties of special soils such as frozen soil, and can make the discrete element method more accurate in simulating the creep mechanical properties of special soils such as frozen soil. It can be seen that the simulation method of the frozen soil non-attenuated creep discrete element contact model in this embodiment has the following advantages:

[0124] 1. On the basis of the original Burgers model, the actual characteristics of damage inside the frozen soil caused by the effect of loading time are taken into account. By adding damage elements, the Burgers creep model considering damage is realized, making its mechanical properties more accurate in describing the frozen soil.

[0125] 2. By using the central difference method, the contact force and displacement formulas in the states before and after the onset of damage are obtained, which can be further used for secondary development of the PFC contact model in Visual studio 2017 and Qt 5.12.0.

[0126] 3. The existing contact model cannot simulate the third stage of the creep curve of frozen soil, and the simulation results are not accurate. The contact model obtained through secondary development can better simulate the third stage of creep of special soils such as frozen soil, which is of great significance for the application expansion of the discrete element method.

[0127] A computer program for executing the simulation method of the non-attenuation creep discrete element contact model of frozen soil in this embodiment can be written and written into a computer device or storage medium. When the computer program is read out and executed, the simulation method of the non-attenuation creep discrete element contact model of frozen soil in this embodiment is executed, thereby achieving the same technical effect as the simulation method of the non-attenuation creep discrete element contact model of frozen soil in the embodiment.

[0128] It should be noted that, unless otherwise specified, when a feature is referred to as "fixed" or "connected" to another feature, it may be directly fixed or connected to the other feature, or it may be indirectly fixed or connected to the other feature. In addition, the descriptions of up, down, left, right, etc. used in the present disclosure are only relative to the relative positional relationship of the components of the present disclosure in the accompanying drawings. The singular forms of "a" and "the" used in the present disclosure are also intended to include the plural forms, unless the context clearly indicates other meanings. In addition, unless otherwise defined, all technical and scientific terms used in this embodiment have the same meaning as those generally understood by technicians in this technical field. The terms used in the specification of this embodiment are only for describing specific embodiments, not for limiting the present invention. The term "and / or" used in this embodiment includes any combination of one or more related listed items.

[0129] It should be understood that, although the term first, second, third etc. may be adopted to describe various elements in the present disclosure, these elements should not be limited to these terms. These terms are only used to distinguish the same type of elements from each other. For example, without departing from the scope of the present disclosure, the first element may also be referred to as the second element, and similarly, the second element may also be referred to as the first element. The use of any and all examples or exemplary language ("for example", "such as" etc.) provided by the present embodiment is only intended to better illustrate embodiments of the present invention, and unless otherwise required, will not impose limitations on the scope of the present invention.

[0130] It should be appreciated that embodiments of the present invention may be implemented or enforced by computer hardware, a combination of hardware and software, or by computer instructions stored in a non-transitory computer-readable memory. The method may be implemented in a computer program using standard programming techniques - including a non-transitory computer-readable storage medium configured with a computer program, wherein the storage medium so configured causes the computer to operate in a specific and predefined manner - according to the methods and drawings described in the specific embodiments. Each program may be implemented in a high-level procedural or object-oriented programming language to communicate with a computer system. However, if desired, the program may be implemented in assembly or machine language. In any case, the language may be a compiled or interpreted language. In addition, the program may be run on a programmed dedicated integrated circuit for this purpose.

[0131] In addition, the operations of the process described in this embodiment may be performed in any suitable order, unless otherwise indicated in this embodiment or otherwise clearly contradicted by the context. The process described in this embodiment (or variations and / or combinations thereof) may be performed under the control of one or more computer systems configured with executable instructions, and may be implemented as codes (e.g., executable instructions, one or more computer programs, or one or more applications) executed jointly on one or more processors, by hardware, or a combination thereof. A computer program includes multiple instructions that may be executed by one or more processors.

[0132] Further, the method can be implemented in any type of computing platform that is operably connected to a suitable computer, including but not limited to a personal computer, a minicomputer, a mainframe, a workstation, a network or distributed computing environment, a separate or integrated computer platform, or in communication with a charged particle tool or other imaging device, etc. Various aspects of the present invention can be implemented in machine-readable code stored on a non-transitory storage medium or device, whether removable or integrated into a computing platform, such as a hard disk, an optical read and / or write storage medium, a RAM, a ROM, etc., so that it can be read by a programmable computer, and when the storage medium or device is read by the computer, it can be used to configure and operate the computer to perform the process described herein. In addition, the machine-readable code, or part thereof, can be transmitted via a wired or wireless network. When such media includes instructions or programs that implement the above steps in conjunction with a microprocessor or other data processor, the invention of this embodiment includes these and other different types of non-transitory computer-readable storage media. When programmed according to the methods and techniques of the present invention, the present invention also includes the computer itself.

[0133] The computer program can be applied to input data to perform the functions of the present embodiment, thereby converting the input data to generate output data stored in a non-volatile memory. The output information can also be applied to one or more output devices such as a display. In a preferred embodiment of the present invention, the converted data represents a physical and tangible object, including a specific visual depiction of the physical and tangible object produced on the display.

[0134] The above are only preferred embodiments of the present invention. The present invention is not limited to the above embodiments. As long as the technical effects of the present invention are achieved by the same means, any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included in the scope of protection of the present invention. Within the scope of protection of the present invention, its technical solutions and / or implementation methods may have various modifications and changes.

Claims

1. A simulation method for a non-attenuated creep discrete element contact model of frozen soil, characterized in that: The simulation method of the frozen soil non-attenuated creep discrete element contact model includes: Establish the first viscoelastic model of frozen soil; Establishing a damaged element; the stiffness of the damaged element is the damaged element stiffness, and the damaged element stiffness changes with time in a period of time after the damage begins; A discrete element contact model of frozen soil non-attenuated creep is established according to the first viscoelastic model and the damage element; Performing frozen soil creep simulation according to the frozen soil non-attenuated creep discrete element contact model; The stiffness of the damaged element is constant in a first time period and changes with time in a second time period; the first time period is a time period before the second time period; The stiffness of the damaged element decreases as time increases during the second time period; The first time period is a time period indicating that the damaged element has not entered a damaged state, and the second time period is a time period indicating that the damaged element has entered a damaged state; The frozen soil creep simulation according to the frozen soil non-attenuated creep discrete element contact model comprises: The non-attenuated creep discrete element contact model of frozen soil is derived using the central difference method to obtain the contact force-displacement relationship information between frozen soil particles; Obtaining a first source code corresponding to the first viscoelastic model; According to the contact force-displacement relationship information, secondary development is performed on the first source code to obtain a second source code; Compiling the second source code to obtain a dynamic link library file; Importing the dynamic link library file into simulation software for simulation; The contact force-displacement relationship information corresponding to the first time period includes: ； in, The first time period The contact force of the damaged component at the moment, The first time period The contact force of the damaged component at time +1, The first time period The total displacement of the frozen soil non-attenuated creep discrete element contact model at time, The first time period The total displacement of the frozen soil non-decayed creep discrete element contact model at time +1, The first time period The displacement of the viscous element of the Kelvin model in the Burgers model at the moment; , , , satisfy ， ； and Determined by the Kelvin model in the Burgers model, and Determined by the Maxwell model in the Burgers model; The contact force-displacement relationship information corresponding to the second time period includes: ； in, The second time period The contact force of the damaged component at the moment, The second time period The contact force of the damaged component at time +1, The second time period The total displacement of the frozen soil non-attenuated creep discrete element contact model at time, The second time period The total displacement of the frozen soil non-decayed creep discrete element contact model at time +1, The second time period The displacement of the viscous element of the Kelvin model in the Burgers model at time t, The second time period The displacement of the damaged element at the time , , satisfy , satisfy ， ， , satisfy K D =E ， ， , satisfy 。 2. The simulation method of the frozen soil non-attenuated creep discrete element contact model according to claim 1 is characterized in that: The step of establishing a first viscoelastic model of frozen soil comprises: Establish the Burgers model of frozen soil; The Burgers model is used as the first viscoelastic model.

3. The simulation method of the frozen soil non-attenuated creep discrete element contact model according to claim 2 is characterized in that: The step of establishing a frozen soil non-attenuated creep discrete element contact model based on the first viscoelastic model and the damage element comprises: The first viscoelastic model and the damage element are connected in series to obtain the frozen soil non-attenuated creep discrete element contact model.

4. The simulation method of the frozen soil non-attenuation creep discrete element contact model according to any one of claims 1 to 3 is characterized in that: The frozen soil creep simulation according to the frozen soil non-attenuated creep discrete element contact model also includes: The dynamic link library file is double-ball verified using the simulation software.

Images