All-direction isobaric simulation method and device, electronic equipment and storage medium
The flexible constraint model was constructed by the discrete unit method and the synchronous pressure processing was performed using the command flow code, which solved the problem of discoherence of axial and lateral pressures in the three-axis test, and achieved more accurate isotropic isopressurization simulation.
Patent Information
- Application Number
- CN202510217753.6
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-02-26
- Publication Date
- 2025-06-27
AI Technical Summary
In the prior art, the axial pressure and lateral pressure loading mechanisms of the triaxial test are not coordinated, resulting in the sample being synchronized from the synchronized isobaric pressure, and it is impossible to accurately reflect the meticulous mechanical characteristics of the sample during the anisobaric isobaric loading and unloading process.
The flexible constraint model is constructed through the discrete unit method, including the top loading plate, the bottom loading plate and the lateral constraint structure, and the command flow code is used to perform axial and lateral synchronous pressure loading and unloading processing to ensure the accuracy of the isotropic isopressurization simulation results.
The accuracy of the simulation is improved, ensuring the coordinated loading or unloading of axial and lateral pressures is ensured, and the problems of slow loading and different stresses of traditional servo mechanisms can be avoided, and the meticulous mechanical properties of the test materials can be more accurately simulated during the loading and unloading process.
Smart Images
Figure CN120217485A_ABST
Abstract
Description
Technical Field
[0001] This application relates to the technical field of geotechnical engineering, and particularly to an isotropic pressure simulation method, device, electronic device, and storage medium. Background Art
[0002] With the rapid development of geotechnical engineering, the research on the mechanical properties of geotechnical materials has become increasingly refined. In engineering practice, geotechnical materials are often subjected to complex stress conditions, and these stress paths may not be limited to uniaxial or triaxial loading. Especially in frozen soil areas and high-stress environments, the stress response and deformation law of materials are directly related to the safety and durability of engineering structures. Therefore, the isotropic pressure loading and unloading test, as an important mechanical test method, has been gradually widely used to analyze the deformation characteristics and failure mechanisms of geotechnical materials under complex stress paths.
[0003] In related technologies, there is a triaxial test for rigid constraints, but the axial pressure and lateral pressure loading mechanisms of this method are not coordinated, resulting in asynchronous isotropic pressure of the specimen, which does not conform to the loading process of the specimen in the laboratory test and cannot accurately reflect the meso-mechanical properties of the specimen during the isotropic pressure loading and unloading process. In summary, the technical problems existing in related technologies need to be improved. Summary of the Invention
[0004] The main purpose of the embodiments of this application is to propose an isotropic pressure simulation method, device, electronic device, and storage medium, which can improve the accuracy of simulation.
[0005] To achieve the above object, on the one hand, an isotropic pressure simulation method is proposed in the embodiments of this application. The method includes:
[0006] Construct a flexible constraint model according to the discrete element method. The flexible constraint model includes a top loading plate, a bottom loading plate, and a lateral constraint structure;
[0007] Perform synchronous pressure loading and unloading on the axial and lateral directions of the flexible constraint model according to the command stream code to obtain an isotropic pressure simulation result. The command stream code is used to synchronously control the top loading plate, the bottom loading plate, and the lateral constraint structure.
[0008] In some embodiments, the constructing a flexible constraint model according to the discrete element method includes the following steps:
[0009] Determine a wall generation command and a structure generation command according to the discrete element method;
[0010] Construct the top loading plate and the bottom loading plate according to the wall generation command;
[0011] The lateral constraint structure is constructed according to the command generated based on the said structure;
[0012] The flexible constraint model is constructed based on the said top loading plate, the said bottom loading plate and the said lateral constraint structure.
[0013] In some embodiments, constructing the said top loading plate and the said bottom loading plate according to the command generated based on the wall body includes the following steps:
[0014] Performing area simulation processing according to the command generated based on the wall body to generate two parallel planar walls;
[0015] Performing servo control system establishment processing on the said planar walls to obtain the said top loading plate and the said bottom loading plate.
[0016] In some embodiments, constructing the said lateral constraint structure according to the command generated based on the structure includes the following steps:
[0017] Performing finite element simulation processing according to the command generated based on the structure to obtain shell elements;
[0018] Performing material definition and condition constraint processing on the said shell elements to obtain the said lateral constraint structure.
[0019] In some embodiments, the command stream code includes a loading function and an unloading function. Processing the flexible constraint model with synchronous axial and lateral pressure loading and unloading according to the command stream code to obtain the isotropic pressure simulation result includes the following steps:
[0020] Performing synchronous axial and lateral pressure loading processing on the flexible constraint model according to the loading function to obtain a pressure loading result;
[0021] Performing synchronous axial and lateral pressure unloading processing on the flexible constraint model according to the unloading function to obtain a pressure unloading result;
[0022] Obtaining the isotropic pressure simulation result based on the pressure loading result and the pressure unloading result.
[0023] In some embodiments, performing synchronous axial and lateral pressure loading processing on the flexible constraint model according to the loading function to obtain a pressure loading result includes the following steps:
[0024] Obtaining the current stress, the target stress and the stress increment according to the loading function;
[0025] Judge the current stress and the target stress according to the first cycle condition. When the current stress is less than the target stress, update the current stress according to the stress increment, and cyclically control the top loading plate, the bottom loading plate and the lateral constraint structure of the flexible constraint model to apply pressure until the current stress is equal to the target stress, and obtain the pressure loading result.
[0026] In some embodiments, the cyclically controlling the top loading plate, the bottom loading plate and the lateral constraint structure of the flexible constraint model to apply pressure includes the following steps:
[0027] Apply pressure to the top loading plate and the bottom loading plate according to the wall servo command;
[0028] At the same time, apply confining pressure to the lateral constraint structure according to the shell element command.
[0029] To achieve the above object, another aspect of the embodiments of the present application provides an isotropic pressure simulation device, the device includes:
[0030] A first module, configured to construct a flexible constraint model according to the discrete element method, the flexible constraint model includes a top loading plate, a bottom loading plate and a lateral constraint structure;
[0031] A second module, configured to perform synchronous pressure loading and unloading on the flexible constraint model axially and laterally according to the command stream code to obtain an isotropic pressure simulation result, and the command stream code is used to synchronously control the top loading plate, the bottom loading plate and the lateral constraint structure.
[0032] To achieve the above object, another aspect of the embodiments of the present application provides an electronic device, the electronic device includes a memory and a processor, the memory stores a computer program, and when the processor executes the computer program, the method described above is implemented.
[0033] To achieve the above object, another aspect of the embodiments of the present application provides a computer-readable storage medium, the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the method described above is implemented.
[0034] The embodiments of the present application at least include the following beneficial effects: The present application provides an isotropic pressure simulation method, device, electronic device and storage medium. According to the discrete element method, a flexible constraint model is constructed. The flexible constraint model includes a top loading plate, a bottom loading plate and a lateral constraint structure. By constructing the flexible constraint model, the efficiency of applying pressure to the lateral constraint structure is increased, and the calculation efficiency is improved. In addition, according to the command stream code, axial and lateral synchronous pressure loading and unloading processing are performed on the flexible constraint model to obtain the isotropic pressure simulation result. The command stream code is used to synchronously control the top loading plate, the bottom loading plate and the lateral constraint structure, and can synergistically increase the axial pressure and the lateral pressure in the discrete element numerical simulation study of the isotropic pressure loading and unloading test, avoiding the problems of slow loading and different stresses during the loading process of the traditional servo mechanism, and can more accurately simulate the meso-mechanical properties of the test material during the loading and unloading process. Description of the Drawings
[0035] Figure 1 is the flowchart of the isotropic pressure simulation method provided by the embodiments of the present application;
[0036] Figure 2 is the structural schematic diagram of the flexible constraint model provided by the embodiments of the present application;
[0037] Figure 3 is the stress-strain curve relationship diagram of loading and unloading provided by the embodiments of the present application;
[0038] Figure 4 is the structural schematic diagram of the isotropic pressure simulation device provided by the embodiments of the present application;
[0039] Figure 5 is the hardware structural schematic diagram of the electronic device provided by the embodiments of the present application. Detailed Embodiments
[0040] In order to make the objectives, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application. When the following description refers to the drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with the embodiments of the present application. They are only examples of devices and methods consistent with some aspects of the embodiments of the present application detailed in the appended claims.
[0041] It can be understood that the terms "first", "second", etc. used in this application may be used herein to describe various concepts, but unless otherwise specified, these concepts are not limited by these terms. These terms are only used to distinguish one concept from another. For example, without departing from the scope of the embodiments of this application, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Depending on the context, the words "if", "when" as used herein may be interpreted as "when...", "while...", or "in response to determining".
[0042] The terms "at least one", "a plurality of", "each", "any one", etc. used in this application, at least one includes one, two or more than two, a plurality of includes two or more than two, each refers to each one of the corresponding plurality, and any one refers to any one of the plurality.
[0043] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by those of ordinary skill in the technical field to which this application belongs. The terms used herein are only for the purpose of describing the embodiments of this application and are not intended to limit this application.
[0044] Before elaborating on the embodiments of this application in detail, some nouns and terms involved in the embodiments of this application are first explained, and the nouns and terms involved in the embodiments of this application are applicable to the following explanations.
[0045] 1) The Discrete Element Method (DEM) is a numerical simulation method used to simulate the movement and interaction of a large number of discrete particles under different conditions. This method is usually used to study the mechanical behavior of particulate materials, such as particle flow, particle packing, particle dispersion, etc. It regards the rock mass as a discrete medium, allows the rock blocks to move and deform, so as to simulate the characteristics of nonlinear large deformation.
[0046] 2) In geotechnical engineering, rigid constraints usually manifest as strict restrictions on the movement of geotechnical bodies or structures, such as displacement, deformation, and rotation. Such restrictions can be provided by rigid materials such as rocks, soils, concrete, steel, etc., or can be provided by artificial structures such as retaining walls and support structures.
[0047] 3) In geotechnical engineering, flexible constraints mainly refer to the use of soft and deformable materials or structures to limit the deformation of geotechnical bodies or structures, while allowing them to have a certain amount of displacement and adaptability.
[0048] With the rapid development of geotechnical engineering, the research on the mechanical properties of rock and soil masses has become increasingly refined. In engineering practice, rock and soil masses are often subjected to complex stress conditions, and these stress paths may not be limited to uniaxial or triaxial loading. Especially in permafrost regions and high-stress environments, the stress response and deformation law of materials are directly related to the safety and durability of engineering structures. Therefore, the isotropic compression loading and unloading test, as an important mechanical test method, has been gradually widely used to analyze the deformation characteristics and failure mechanisms of rock and soil materials under complex stress paths.
[0049] However, the current research on the isotropic compression loading and unloading test focuses on experimental research, lacking the study of the meso-mechanical properties of soil during the loading and unloading process. Existing technologies usually focus on rigid confinement triaxial tests. The isotropic compression loading and unloading test is applied through the upper and lower loading plates of the specimen and the rigid lateral confinement wall. The servo mechanism of the upper and lower loading plates and the rigid lateral confinement wall is realized by writing command flow codes. The traditional method of isotropic compression loading and unloading does not have a coordinated loading mechanism for axial pressure and lateral pressure. They are loaded separately, resulting in the non-synchronization of isotropic pressure on the specimen, which does not conform to the loading process of the specimen in the laboratory test and cannot accurately reflect the meso-mechanical properties of the specimen during the isotropic compression loading and unloading process.
[0050] In view of this, in the embodiments of the present application, an isotropic compression simulation method, device, electronic device, and storage medium are provided. The scheme constructs a flexible confinement model based on the discrete element method. The flexible confinement model includes a top loading plate, a bottom loading plate, and a lateral confinement structure. By constructing the flexible confinement model, the efficiency of applying pressure to the lateral confinement structure is increased, and the calculation efficiency is improved. In addition, the scheme performs synchronous pressure loading and unloading on the flexible confinement model axially and laterally according to the command flow code to obtain the isotropic compression simulation result. The command flow code is used to synchronously control the top loading plate, the bottom loading plate, and the lateral confinement structure, and can coordinate the increase of axial pressure and lateral pressure in the discrete element numerical simulation study of the isotropic compression loading and unloading test, avoiding the problems of slow loading and different stresses during the loading process in the traditional servo mechanism, and being able to more accurately simulate the meso-mechanical properties of the test material during the loading and unloading process.
[0051] An isotropic pressure simulation method provided by an embodiment of the present application relates to the field of information technology. The isotropic pressure simulation method provided by an embodiment of the present application can be applied to a terminal, a server, or software running on a terminal or a server. In some embodiments, the terminal can be a smart phone, a tablet computer, a notebook computer, a desktop computer, a smart speaker, a smart watch, a vehicle-mounted terminal, etc., but is not limited thereto; the server side can be configured as an independent physical server, or can be configured as a server cluster or a distributed system composed of multiple physical servers, or can be configured as a cloud server providing basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communications, middleware services, domain name services, security services, CDN, and big data and artificial intelligence platforms. The server can also be a node server in a blockchain network; the software can be an application implementing an isotropic pressure simulation method, etc., but is not limited to the above forms.
[0052] The present application can be used in many general or special computer system environments or configurations. For example: personal computers, server computers, handheld or portable devices, tablet-type devices, multi-processor systems, microprocessor-based systems, set-top boxes, programmable consumer electronic devices, network PCs, minicomputers, mainframe computers, distributed computing environments including any of the above systems or devices, and so on. The present application can be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. The present application can also be practiced in a distributed computing environment where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules can be located in local and remote computer storage media including storage devices.
[0053] Figure 1 is an optional flowchart of an isotropic pressure simulation method provided by an embodiment of the present application, Figure 1 The method in can include but is not limited to steps S101 to S102.
[0054] Step S101, constructing a flexible constraint model according to the discrete element method, where the flexible constraint model includes a top loading plate, a bottom loading plate, and a lateral constraint structure;
[0055] Step S102, performing synchronous pressure loading and unloading on the flexible constraint model axially and laterally according to the command stream code to obtain an isotropic pressure simulation result, where the command stream code is used to synchronously control the top loading plate, the bottom loading plate, and the lateral constraint structure.
[0056] Steps S101 to S102 shown in the embodiments of the present application construct a flexible constraint model through the discrete element method. Among them, the flexible constraint model includes a top loading plate, a bottom loading plate, and a lateral constraint structure. The lateral constraint structure is constructed by using shell elements in finite elements, and stress can be directly applied to the shell elements. According to the target stress input by the target, the stress of the lateral shell elements is applied step by step according to the increment size, avoiding the defects of the traditional servo loading mechanism. Then, according to the command flow code, synchronous pressure loading and unloading processes in the axial and lateral directions are performed on the flexible constraint model, which can unify the axial pressurization and lateral pressurization processes, enabling the multi-directional pressures to be loaded or unloaded cooperatively. In the embodiments of the present application, the rigid constraint of the specimen is changed to a flexible constraint. The top and bottom loading plates of the specimen are constructed through the wall generation command in the discrete element software. Using the structure generation command, a shell element body is established as the lateral constraint wall, and the command flow code is written to unify the axial pressurization and lateral pressurization processes, enabling the multi-directional pressures to be loaded or unloaded cooperatively. To make the multi-directional pressure loading and unloading more accurate, the corresponding loading and unloading command flow code is written. By setting the initial stress and the target stress, and specifying the pressure increment at the same time, as the calculation time increases, the multi-directional pressures are loaded step by step to the target stress state according to the size of the increment, which is closer to the process of uniformly applying pressure in the indoor test.
[0057] In step S101 of some embodiments, the flexible constraint model constructed according to the discrete element method may include the following steps:
[0058] Determine the wall generation command and the structure generation command according to the discrete element method;
[0059] Construct the top loading plate and the bottom loading plate according to the wall generation command;
[0060] Construct the lateral constraint structure according to the structure generation command;
[0061] Construct the flexible constraint model according to the top loading plate, the bottom loading plate, and the lateral constraint structure.
[0062] In the embodiments of the present application, according to the discrete element method, a corresponding flexible constraint model can be constructed for the test material using the corresponding discrete element software. Among them, the discrete element software is mainly used to simulate and analyze the mechanical behavior and motion characteristics of granular materials, and software such as EDEM, YADE, and PFC3D can be used. The test material refers to the target that needs to be subjected to isotropic compression simulation, such as geotechnical materials. By using a flexible constraint structure to perform an isotropic compression model on the test material, the stress response and deformation law of the material can be analyzed. In the embodiments of the present application, by changing the rigid constraint of the specimen material to a flexible constraint, specifically, the top and bottom loading plates of the specimen are constructed using the wall generation command in the discrete element software, and the lateral constraint structure is constructed using the structure generation command. The lateral constraint structure is a lateral constraint wall obtained by establishing a shell element body. Specifically, the top loading plate and the bottom loading plate are constructed according to the wall generation command. The wall usually serves as a boundary condition or an obstacle in the discrete element simulation, and can be fixed or movable. The wall generation command needs to specify parameters such as the type, position, size, and attributes of the wall to generate the simulated wall, and obtain the top loading plate and the bottom loading plate. Then, the lateral constraint structure is constructed according to the structure generation command. The structure generation command is usually used to create complex geometric shapes or structures, such as spheres, blocks, polyhedrons, etc. Since the test material may be of an irregular shape, when performing lateral pressure loading and unloading, by establishing a flexible constraint lateral constraint structure based on the shape of the test material, it is possible to better wrap the test material for isotropic compression simulation. As Figure 2 shown, a flexible constraint model is constructed according to the generated top loading plate, bottom loading plate, and lateral constraint structure.
[0063] One of the technical solutions in the above technical solutions has the following advantages or beneficial effects: The flexible constraint model constructed by the discrete element method in the embodiments of the present application can better monitor the meso-mechanical properties of the specimen material and can accurately reflect the meso-mechanical properties of the specimen during the isotropic compression loading and unloading process.
[0064] In some embodiments, the step of constructing the top loading plate and the bottom loading plate according to the wall generation command includes the following steps:
[0065] Perform regional simulation processing according to the wall generation command to generate two parallel plane walls;
[0066] Perform servo control system establishment processing on the plane wall to obtain the top loading plate and the bottom loading plate.
[0067] In the embodiments of the present application, first, the dimensions of the specimen material, as well as the dimensions and positions of the loading plates, are determined according to the requirements of the isotropic pressure simulation. Then, the corresponding top loading plate and bottom loading plate are constructed according to the wall generation command. Specifically, regional simulation processing is carried out according to the wall generation command to generate a calculation region slightly larger than the dimensions of the specimen material to generate a wall. For the top and bottom loading plates, generally, two parallel planar walls need to be generated, and the dimensions and positions of the walls are set to ensure that they can completely cover the top and bottom of the specimen. At the same time, the properties of the walls, such as the friction coefficient and stiffness, are adjusted according to the needs. Then, servo control system establishment processing is carried out on the planar walls. Through the servo control system, the moving speeds of the top loading plate and the bottom loading plate can be controlled to maintain the constancy of the stress on the specimen material.
[0068] In some embodiments, constructing the lateral constraint structure according to the structure generation command includes the following steps:
[0069] Carrying out finite element simulation processing according to the structure generation command to obtain shell elements;
[0070] Carrying out material definition and condition constraint processing on the shell elements to obtain the lateral constraint structure.
[0071] In the embodiments of the present application, finite element analysis is carried out according to the structure generation command. Shell elements are established based on the geometric dimensions, material properties, etc. of the test material. The geometric shape of the shell elements can be created using the geometric modeling tools of the software. The shell elements can simulate thin plate structures and can consider bending and membrane effects. Then, the material properties and section properties of the shell elements are defined, and boundary conditions are applied at appropriate positions of the shell elements to simulate the actual constraint situation of the wall. For example, the shell elements are fixed to ensure that the wall has the required constraint conditions in the lateral direction. For example, displacement constraints or force loads are applied at the edges or specific nodes of the wall, and finally, the lateral constraint structure is obtained.
[0072] In some embodiments, the command stream code includes a loading function and an unloading function. Carrying out axial and lateral synchronous pressure loading and unloading processing on the flexible constraint model according to the command stream code to obtain the isotropic pressure simulation result includes the following steps:
[0073] Carrying out axial and lateral synchronous pressure loading processing on the flexible constraint model according to the loading function to obtain the pressure loading result;
[0074] Carrying out axial and lateral synchronous pressure unloading processing on the flexible constraint model according to the unloading function to obtain the pressure unloading result;
[0075] Obtaining the isotropic pressure simulation result according to the pressure loading result and the pressure unloading result.
[0076] In the embodiments of the present application, in order to make the loading and unloading of the pressure in all directions more accurate, the corresponding loading and unloading command stream codes are written. By setting the command stream codes, the top loading plate, the bottom loading plate and the lateral constraint structure of the flexible constraint model are synchronously controlled, so as to perform isotropic pressure simulation on the test material. Among them, the command stream code includes a loading function and an unloading function. The loading function sets the initial stress and the target stress, and at the same time specifies the pressure increment. As the calculation time increases, the pressure in all directions is loaded step by step to the target stress state according to the size of the increment, which is closer to the process of applying pressure uniformly in the indoor test. Similarly, the unloading function sets the initial stress and the target stress, and at the same time specifies the pressure increment. As the calculation time increases, the pressure in all directions is unloaded step by step to the target stress state according to the size of the increment, which is closer to the process of reducing the pressure uniformly in the indoor test. The embodiments of the present application perform synchronous pressure loading on the axial and lateral directions of the flexible constraint model according to the loading function, so as to simulate the process of applying pressure uniformly in the test. It can be imagined that the embodiments of the present application can perform multiple isotropic pressure simulation processes on the test material according to the loading function and the unloading function to obtain the corresponding isotropic pressure simulation results.
[0077] One of the technical solutions in the above technical solutions has the following advantages or beneficial effects: In the embodiments of the present application, by writing the command stream code, the axial pressurization process and the lateral pressurization process are unified, so that the pressures in all directions can be loaded or unloaded cooperatively, and the loading and unloading of the pressures in all directions are more accurate.
[0078] In some embodiments, the performing synchronous pressure loading on the axial and lateral directions of the flexible constraint model according to the loading function to obtain a pressure loading result includes the following steps:
[0079] Obtain the current stress, the target stress and the stress increment according to the loading function;
[0080] Judge the current stress and the target stress according to the first loop condition. When the current stress is less than the target stress, update the current stress according to the stress increment, and cyclically control the top loading plate, the bottom loading plate and the lateral constraint structure of the flexible constraint model to perform pressurization processing until the current stress is equal to the target stress to obtain the pressure loading result.
[0081] In the embodiment of the present application, a loading function is defined, and input variables are determined: current stress, target stress, stress increment. A loop statement is used, and the judgment condition is that the current stress is less than the target stress. When the stress meets the judgment condition, the current stress is increased step by step according to the magnitude of the stress increment. According to the current stress value in each cycle, confining pressure is directly applied to the shell element. According to the contact area between the top and bottom loading plates and the specimen and the target stress value, a command of wall servo force value is used for the top and bottom walls, and the target contact force is directly applied. The cycle continues until the target stress value is reached, and the pressure loading result is obtained. Similarly, by defining an unloading function, input variables are determined: current stress, target stress, stress increment. A loop statement is used, and the judgment condition is that the current stress is greater than the target stress. When the stress meets the judgment condition, the current stress is decreased step by step according to the magnitude of the stress increment. According to the current stress value in each cycle, confining pressure is directly applied to the shell element. According to the contact area between the top and bottom loading plates and the specimen and the target stress value, a command of wall servo force value is used for the top and bottom walls, and the target contact force is directly applied. The cycle continues until the target stress value is unloaded. By performing pressure loading and unloading on the test material in the embodiment of the present application, the stress-strain curve relationship of loading and unloading can be obtained, as Figure 3 shown. By writing the command flow code for loading and unloading in the embodiment of the present application, when the stress meets the judgment condition, the current stress is increased or decreased step by step according to the magnitude of the stress increment. According to the current stress value in each cycle, confining pressure is directly applied to the shell element. According to the contact area between the top and bottom loading plates and the specimen and the target stress value, a command of wall servo force value is used for the top and bottom walls, and the target contact force is directly applied. The cycle continues until the target stress value is loaded or unloaded. It can realize the coordinated increase of axial pressure and lateral pressure in the discrete element numerical simulation study of isotropic equal-pressure loading and unloading tests, ensure the synchronous increase of all-directional pressures, and more accurately simulate the action of geotechnical materials under hydrostatic pressure.
[0082] In some embodiments, the cycle controls the top loading plate, the bottom loading plate, and the lateral constraint structure of the flexible constraint model to perform pressure application, including the following steps:
[0083] Perform pressure application on the top loading plate and the bottom loading plate according to the wall servo command;
[0084] At the same time, perform confining pressure application on the lateral constraint structure according to the shell element command.
[0085] In the embodiments of the present application, for the top loading plate and the bottom loading plate, the wall servo command is directly used, and the top loading plate and the bottom loading plate are moved through the servo control system to apply pressure to the test material. At the same time, according to the shell element command, the lateral restraint structure is subjected to confining pressure application treatment. By using the shell element to load the confining pressure of the specimen, the stress magnitude on the surface of the shell element can be directly applied by the Command command. The Command command is used to call and execute the specified instruction, which greatly simplifies the complexity of the code. At the same time, it can effectively avoid the disadvantages of the traditional servo loading mechanism, increase the efficiency of applying the lateral restraint wall pressure, and improve the calculation efficiency.
[0086] Next, in combination with specific application examples, the solutions of the embodiments of the present application will be introduced and described in detail:
[0087] The embodiments of the present application can be applied to the field of geotechnical engineering technology and are applicable to the simulation test of test materials, such as geotechnical materials. In the embodiments of the present application, the rigid constraint of the specimen is changed to a flexible constraint. The top and bottom loading plates of the specimen are constructed through the wall generation command in the discrete element software. The structure generation command is used to establish a Shell element body as the lateral constraint wall. The axial pressurization and lateral pressurization processes are unified by writing command flow codes, so that the stresses in all directions can be loaded or unloaded synergistically. Specifically, by defining a loading function and an unloading function, the loading process of the specimen is realized through the loading function, and the unloading process of the specimen is realized through the unloading function. According to the numerical simulation requirements, the current stresses, target stresses, and increments in the axial and lateral directions of the specimen are input. Usually, the current stress value can be directly obtained. It is 0 before the start of loading and becomes the target value after loading to the target stress. The current stress value can be directly determined during the simulation process. The input data is judged. If the current stress is less than the target stress, by controlling the top and bottom loading plates and the lateral constraint wall, the top loading plate moves downward, the bottom loading plate moves upward, and the lateral constraint wall shrinks inward, so that the stresses in all directions increase. The unloading process is the opposite. The current stress is gradually increased according to the increment, ensuring that within the specified time step, the current stress can increase by the corresponding stress increment, and the process is repeated continuously until the current stress of the specimen reaches the target stress, then the loading stops. During the loading process, to ensure the uniform application of the axial stress and the lateral stress, the command of the shell element is used for the outer flexible membrane, that is, the lateral wall adopts the shell element in the finite element. Utilizing the characteristic that stress can be directly applied to the shell element, according to the target stress input by the target, the stress of the lateral shell element is applied step by step according to the increment size. This avoids the defects of the traditional servo loading mechanism. The traditional servo loading mechanism uses FISH language to write command flow codes to control the movement of the loading plates in all directions. Since the FISH code has a lower operation efficiency and slower calculation compared to the Command command, and it is prone to particle overflow during the loading process. When using the shell element to load the confining pressure of the specimen, the Command command can be used. The specific command is the "structure shell apply[begin]" command. Among them, structure represents the context or field of the command, that is, structural analysis. Structural analysis is a branch of engineering mechanics used to determine the response of a structure under external forces, such as stress, strain, and displacement. Shell is used to indicate that the next operation is for the shell structure. apply is used to indicate that the next operation is to apply a certain condition or load to the structure. In structural analysis, this usually involves applying forces, pressures, temperatures, displacements, etc. [begin] represents an optional parameter or flag of the command, used to indicate the start of a series of operations or settings.The stress magnitude directly applied to the surface of the shell element in the embodiments of the present application greatly simplifies the complexity of the code. At the same time, it can effectively avoid the disadvantages of the traditional servo loading mechanism, improve the efficiency of applying the pressure on the lateral constraint wall and enhance the calculation efficiency. For the top and bottom loading plates, the wall servo Command command is directly used. First, calculate the contact area between the top and bottom loading plates and the specimen. Multiply the two according to the target stress magnitude to calculate the magnitude of the target axial force, and directly apply the same force value to the top and bottom walls. The wall servo command is "wall servo force(0,0,[begin*area])activate true", where wall refers to a boundary or wall in the simulation domain, that is, the top loading plate and the bottom loading plate. Servo refers to a closed-loop control system used to precisely control the output of a certain process or system, which can dynamically adjust the force applied to the wall according to preset conditions. Force refers to the force applied to the wall. (0,0,[begin*area]) represents the definition of the direction and magnitude of the force, where the first two components of the force (the components in the x and y directions) are set to 0, meaning the force is applied along the third direction (assumed to be the z direction). [begin*area] is an expression, where begin is a variable or parameter, area is the area of the wall, and * represents multiplication. This means that the magnitude of the force applied to the wall is proportional to the value of begin and the area of the wall. Through continuous cyclic loading until the specimen reaches the target stress, the unloading process of the specimen is the same. The embodiments of the present application use the shell element as the lateral constraint of the specimen and use the servo command for the upper and lower loading plates, avoiding the deficiencies of the traditional loading method using the command flow servo mechanism, writing the corresponding command flow code, enabling the upper and lower loading plates and the lateral constraint wall to control the loading rate in an incremental loading manner, achieving the purpose of loading or unloading, and at the same time ensuring that the loading stress values in all directions are the same, avoiding the situation of unequal pressures in all directions.
[0088] Please refer to Figure 4 , the embodiments of the present application further provide an isotropic pressure simulation device, which can implement the above isotropic pressure simulation method. The device includes:
[0089] The first module 401 is used to construct a flexible constraint model according to the discrete element method. The flexible constraint model includes a top loading plate, a bottom loading plate, and a lateral constraint structure;
[0090] The second module 402 is used to perform synchronous pressure loading and unloading processing on the axial and lateral directions of the flexible constraint model according to the command flow code to obtain an isotropic pressure simulation result. The command flow code is used to synchronously control the top loading plate, the bottom loading plate, and the lateral constraint structure.
[0091] It can be understood that the content in the above method embodiments is applicable to the device embodiments of the present application. The functions specifically implemented by the device embodiments of the present application are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those of the above method embodiments.
[0092] An embodiment of the present application also provides an electronic device, which includes a memory and a processor. The memory stores a computer program, and when the processor executes the computer program, the above-mentioned isotropic pressure simulation method is implemented. The electronic device can be any intelligent terminal including a tablet computer, a vehicle-mounted computer, etc.
[0093] It can be understood that the content in the above method embodiments is applicable to the device embodiments of the present application. The functions specifically implemented by the device embodiments of the present application are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those of the above method embodiments.
[0094] Please refer to Figure 5 , Figure 5 which shows the hardware structure of an electronic device in another embodiment. The electronic device includes:
[0095] A processor 501, which can be implemented in a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (ASIC), or one or more integrated circuits, etc., and is used to execute relevant programs to implement the technical solutions provided by the embodiments of the present application;
[0096] A memory 502, which can be implemented in the form of a read-only memory (ROM), a static storage device, a dynamic storage device, or a random access memory (RAM), etc. The memory 502 can store an operating system and other application programs. When implementing the technical solutions provided by the embodiments of this specification through software or firmware, the relevant program codes are stored in the memory 502 and are called by the processor 501 to execute the isotropic pressure simulation method of the embodiments of the present application;
[0097] An input / output interface 503, which is used to implement information input and output;
[0098] A communication interface 504, which is used to implement communication interaction between this device and other devices. Communication can be achieved through wired means (such as USB, network cable, etc.) or through wireless means (such as mobile network, WIFI, Bluetooth, etc.);
[0099] A bus 505 that transmits information between the various components of the device (e.g., the processor 501, the memory 502, the input / output interface 503, and the communication interface 504);
[0100] The processor 501 , the memory 502 , the input / output interface 503 and the communication interface 504 are connected to each other in communication within the device via the bus 505 .
[0101] An embodiment of the present application further provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the above-mentioned isotropic pressure simulation method is implemented.
[0102] It can be understood that the contents of the above method embodiments are all applicable to the present storage medium embodiments, the functions specifically implemented by the present storage medium embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.
[0103] The memory, as a non-transient computer-readable storage medium, can be used to store non-transient software programs and non-transient computer executable programs. In addition, the memory may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some embodiments, the memory may optionally include a memory remotely disposed relative to the processor, and these remote memories may be connected to the processor via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.
[0104] The embodiment of the present application provides an isotropic pressure simulation method, device, electronic device and storage medium. The scheme constructs a flexible constraint model according to the discrete element method. The flexible constraint model includes a top loading plate, a bottom loading plate and a lateral constraint structure. By constructing the flexible constraint model, the efficiency of applying pressure to the lateral constraint structure is increased, and the calculation efficiency is improved. In addition, the scheme performs axial and lateral synchronous pressure loading and unloading processing on the flexible constraint model according to the command stream code to obtain an isotropic pressure simulation result. The command stream code is used to synchronously control the top loading plate, the bottom loading plate and the lateral constraint structure, and can coordinately increase the axial pressure and lateral pressure in the discrete element numerical simulation study of the isotropic pressure loading and unloading test, avoiding the problems of slow loading of the traditional servo mechanism and different stresses during the loading process, and can more accurately simulate the microscopic mechanical properties of the test material during the loading and unloading process.
[0105] The embodiments described in the embodiments of the present application are to more clearly illustrate the technical solutions of the embodiments of the present application, and do not constitute a limitation on the technical solutions provided by the embodiments of the present application. As can be known to those skilled in the art, with the evolution of technology and the emergence of new application scenarios, the technical solutions provided by the embodiments of the present application are equally applicable to similar technical problems.
[0106] Those skilled in the art can understand that the technical solutions shown in the figures do not constitute a limitation on the embodiments of the present application, and may include more or fewer steps than shown in the figures, or combine certain steps, or different steps.
[0107] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated, that is, they may be located in one place, or may be distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment.
[0108] Those of ordinary skill in the art can understand that all or some of the steps in the methods disclosed above, and the functional modules / units in the systems and devices, can be implemented as software, firmware, hardware, and their appropriate combinations.
[0109] The terms "first", "second", "third", "fourth", etc. (if any) in the specification of the present application and the above-mentioned drawings are used to distinguish similar objects, and do not necessarily need to be used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged under appropriate circumstances, so that the embodiments of the present application described here can be implemented in an order other than those illustrated or described here. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but may include other steps or units not clearly listed or inherent to these processes, methods, products, or devices.
[0110] It should be understood that in this application, "at least one (item)" means one or more, and "a plurality" means two or more. "And / or" is used to describe the association relationship of associated objects and indicates that there can be three relationships. For example, "A and / or B" can mean: only A exists, only B exists, and both A and B exist at the same time. Here, A and B can be singular or plural. The character " / " generally indicates that the associated objects before and after are in an "or" relationship. "At least one (item) of the following" or a similar expression means any combination of these items, including any combination of single item (item) or plural items (items). For example, at least one (item) of a, b, or c can mean: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", where a, b, c can be single or multiple.
[0111] In several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the above division of units is only a logical function division. In actual implementation, there can be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling or direct coupling or communication connection between each other can be through some interfaces. The indirect coupling or communication connection of devices or units can be in electrical, mechanical or other forms.
[0112] The units described above as separate components may or may not be physically separated. The components displayed as units may or may not be physical units, that is, they can be located in one place or distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.
[0113] In addition, in each embodiment of this application, each functional unit can be integrated in a processing unit, or each unit can exist physically alone, or two or more units can be integrated in one unit. The above integrated units can be implemented in the form of hardware or in the form of software functional units.
[0114] When the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or all or part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes multiple instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods of various embodiments of this application. The aforementioned storage medium includes: various media that can store programs, such as USB flash drives, mobile hard disks, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical discs.
[0115] The preferred embodiments of the embodiments of this application have been described above with reference to the accompanying drawings, which does not limit the scope of the rights of the embodiments of this application. Any modifications, equivalent replacements, and improvements made by those skilled in the art without departing from the scope and essence of the embodiments of this application shall be within the scope of the rights of the embodiments of this application.
Claims
1. An isotropic pressure simulation method, characterized in that: The method comprises the following steps: A flexible constraint model is constructed according to a discrete element method, wherein the flexible constraint model includes a top loading plate, a bottom loading plate and a lateral constraint structure; The flexible constraint model is subjected to axial and lateral synchronous pressure loading and unloading processing according to a command flow code to obtain an isotropic pressure simulation result. The command flow code is used to synchronously control the top loading plate, the bottom loading plate and the lateral constraint structure.
2. The method according to claim 1, characterized in that The flexible constraint model is constructed according to the discrete element method, comprising the following steps: Determine a wall generation command and a structure generation command according to the discrete element method; Constructing the top loading plate and the bottom loading plate according to the wall generation command; Constructing the lateral constraint structure according to the structure generation command; The flexible constraint model is constructed according to the top loading plate, the bottom loading plate and the lateral constraint structure.
3. The method according to claim 2, characterized in that The step of constructing the top loading plate and the bottom loading plate according to the wall generation command comprises the following steps: Performing regional simulation processing according to the wall generation command to generate two parallel plane walls; The servo control system is established for the plane wall to obtain the top loading plate and the bottom loading plate.
4. The method according to claim 2, characterized in that: The step of constructing the lateral constraint structure according to the structure generation command comprises the following steps: Perform finite element simulation processing according to the structure generation command to obtain shell elements; The shell element is subjected to material definition and condition constraint processing to obtain the lateral constraint structure.
5. The method according to claim 1, characterized in that The command flow code includes a loading function and an unloading function. The command flow code is used to perform axial and lateral synchronous pressure loading and unloading processing on the flexible constraint model to obtain an isotropic pressure simulation result, including the following steps: Performing synchronous axial and lateral pressure loading processing on the flexible constraint model according to the loading function to obtain a pressure loading result; According to the unloading function, the flexible constraint model is subjected to synchronous axial and lateral pressure unloading processing to obtain a pressure unloading result; The isotropic pressure simulation result is obtained according to the pressure loading result and the pressure unloading result.
6. The method according to claim 5, characterized in that The step of performing axial and lateral synchronous pressure loading processing on the flexible constraint model according to the loading function to obtain a pressure loading result comprises the following steps: Obtaining current stress, target stress and stress increment according to the loading function; The current stress and the target stress are judged according to the first cycle condition; when the current stress is less than the target stress, the current stress is updated according to the stress increment; the top loading plate, the bottom loading plate and the lateral constraint structure of the flexible constraint model are cyclically controlled to perform pressurization until the current stress is equal to the target stress, thereby obtaining the pressure loading result.
7. The method according to claim 6, characterized in that The cyclic control of the top loading plate, the bottom loading plate and the lateral constraint structure of the flexible constraint model to perform a pressurization process comprises the following steps: Pressurizing the top loading plate and the bottom loading plate according to a wall servo command; At the same time, the lateral restraint structure is subjected to confining pressure according to the shell unit command.
8. An isotropic pressure simulation device, characterized in that: The device comprises: A first module is used to construct a flexible constraint model according to a discrete element method, wherein the flexible constraint model includes a top loading plate, a bottom loading plate and a lateral constraint structure; The second module is used to perform axial and lateral synchronous pressure loading and unloading processing on the flexible constraint model according to the command flow code to obtain isotropic pressure simulation results. The command flow code is used to synchronously control the top loading plate, the bottom loading plate and the lateral constraint structure.
9. An electronic device, characterized in that: The electronic device comprises a memory and a processor, the memory stores a computer program, and the processor implements the method according to any one of claims 1 to 7 when executing the computer program.
10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.