Key-based near-field dynamics numerical simulation method, system, equipment and medium

Abstract

The invention discloses a numerical simulation method, system and equipment for bond-based near-field dynamics and a medium. The method comprises the following steps: randomly assigning linear elastic ultimate elongation and tensile ultimate elongation of each substance point of a to-be-analyzed material subjected to quantitative discrete processing by adopting a Weibull function; based on a discrete processing result, calculating a relative position vector, a relative displacement vector, bonding stiffness and bond elongation of any two substance points in the ethnic group in sequence, and constructing a bilinear damage function and a bond-based near-field dynamic constitutive equation; based on a constitutive equation, circularly updating the coordinate information of each material point and the elongation of the key step by step; obtaining a bilinear damage function value according to the updated elongation of the key, and updating the ethnic group information; and calculating a local damage value of each substance point according to the updated group information so as to carry out fracture condition analysis and obtain a damage judgment result. According to the method provided by the invention, the fracture prediction accuracy of the cement stabilizing material is improved, and the strain softening characteristics of the material can be reflected.

Description

Technical Field

[0001] The present invention relates to the field of mechanical engineering technology, and in particular to a bond-based peridynamic numerical simulation method, system, equipment and medium. Background Art

[0002] In the field of materials engineering, bond-based peridynamic numerical simulations are used to study key mechanical properties of materials, such as failure mechanism, fracture toughness, and fatigue life. By simulating the deformation and failure process of materials under different load conditions, researchers can gain a deeper understanding of the influence of the material's microstructure on its macroscopic mechanical behavior.

[0003] Bond-based peridynamics theory is an important branch of peridynamics (PD) theory. This theory builds a model based on the idea of ​​nonlocal action and describes the mechanical behavior of materials by solving spatial integral equations. Bond-based peridynamics theory discretizes the material into multiple material points and uses the positions between the material points to describe the interaction between the points. The existing bond-based peridynamics numerical simulation method is based on the elastic-brittle model of force and deformation, that is, the relationship between the distance and force between two material points fully conforms to the elastic relationship when the distance between the two points is less than a certain value. When the distance exceeds a certain value, it is considered that the interaction force between the material points disappears, that is, the bond is completely broken.

[0004] However, the existing bond-based peridynamics theory uses a certain numerical definition of the fracture parameters of all material points of the research object, so the existing bond-based peridynamics theory is applicable to homogeneous elastic-brittle materials. Cement-stabilized materials are a mixture of sand, gravel and cement, which is an inhomogeneous material, so cement-stabilized materials are quasi-brittle inhomogeneous materials. The existing bond-based peridynamics numerical simulation method is not suitable for simulating the fracture behavior analysis of quasi-brittle inhomogeneous materials.

[0005] It can be seen that how to solve the problem that the existing bond-based peridynamic numerical simulation method is not suitable for simulating the fracture behavior of quasi-brittle inhomogeneous materials has become a technical problem that needs to be urgently solved by those skilled in the art. Summary of the invention

[0006] The present invention provides a bond-based peridynamic numerical simulation method, system, device and medium to solve the technical problem that the existing bond-based peridynamic numerical simulation method is not suitable for simulating the fracture behavior of quasi-brittle inhomogeneous materials, and to provide a method suitable for analyzing the fracture behavior of quasi-brittle inhomogeneous materials, which can accurately simulate the fracture behavior of quasi-brittle inhomogeneous materials.

[0007] In a first aspect, the present invention provides a bond-based peridynamic numerical simulation method, which is applied to the simulation of fracture behavior of quasi-brittle inhomogeneous materials, and the method comprises:

[0008] In the process of material damage judgment, quantitative discrete information of the material to be analyzed is obtained, wherein the quantitative discrete information includes coordinate information, volume information and group information of each material point;

[0009] The linear elastic limit elongation and the tensile limit elongation of each material point are randomly assigned using a Weibull function to obtain a fracture parameter matrix;

[0010] Based on the coordinate information, the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points in the group are calculated in sequence, and a bilinear damage function is constructed according to the fracture parameter matrix and the bond elongation;

[0011] constructing a bond-based peridynamic constitutive equation according to the relative position vector, the relative displacement vector, the bond elongation, the bond stiffness and the bilinear damage function;

[0012] Within a preset total calculation time step, based on the bond-based peridynamic constitutive equation, the volume information and the group information, the coordinate information of each of the material points is cyclically updated time step by time step, and the elongation of the corresponding bond is updated based on the cyclically updated coordinate information to obtain the updated bond elongation of each of the material points;

[0013] Substituting the updated bond elongation into the bilinear damage function to obtain a bilinear damage function value, and updating the group information based on the bilinear damage function value to obtain updated group information;

[0014] The local damage value of each material point is calculated according to the updated group information, and the damage and fracture state of the material to be analyzed at the corresponding time step is obtained according to the local damage value. The calculation of all time steps is completed in a cycle to obtain the bond-based peridynamic simulation result of the damage and fracture behavior of the material to be analyzed.

[0015] Preferably, the step of obtaining the quantitative discrete information of the material to be analyzed includes:

[0016] Obtaining values ​​of various mechanical parameters required for numerical simulation of bond-based peridynamics, wherein the various mechanical parameters at least include: elastic modulus, density, Poisson's ratio, geometric dimensions and material point size of the material to be analyzed;

[0017] Quantifying and discretizing the material to be analyzed according to the size of the material points to obtain coordinate information and volume information of each material point;

[0018] For each of the material points, the neighboring material points within a preset area radius are taken as a group of the material points, and the group information of the group is acquired and stored.

[0019] Preferably, the Weibull function is used to randomly assign the linear elastic limit elongation and the tensile limit elongation of each material point to obtain a fracture parameter matrix, including:

[0020] Constructing a first matrix having the same size as the number of the material points, and randomly assigning values ​​to the first matrix using a Weibull function;

[0021] Acquire test values ​​of the linear elastic limit elongation and the tensile limit elongation of the material to be analyzed, and use the product of the first matrix after random assignment and the test value of the linear elastic limit elongation as the linear elastic limit elongation matrix of each material point, and use the product of the first matrix after random assignment and the test value of the tensile limit elongation as the tensile limit elongation matrix of each material point;

[0022] The linear elastic ultimate elongation matrix and the tensile ultimate elongation matrix are used as fracture parameter matrices.

[0023] Preferably, based on the coordinate information, the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points are calculated in sequence, and a bilinear damage function is constructed according to the fracture parameter matrix and the bond elongation, including:

[0024] Calculating the relative position vector, relative displacement vector and bonding stiffness between any two material points according to the initial coordinate information;

[0025] The relative position vector expression is:

[0026] ξ=x j -x i

[0027] Among them, x j represents the coordinates of material point j, x i represents the coordinates of material point i;

[0028] The relative displacement vector expression is:

[0029] η=u(x j ,t)-u(x i ,t)

[0030] Among them, u(x j ,t) represents the displacement vector of material point j at time step t, u(x i ,t) represents the displacement vector of material point i at time step t;

[0031] The bond stiffness expression is:

[0032]

[0033] Where, E represents the elastic modulus; δ is the field radius;

[0034] According to the relative position vector and the relative displacement vector, the elongation of the key is calculated, and the expression of the elongation of the key is:

[0035]

[0036] A bilinear damage function is constructed according to the linear elastic limit elongation and the tensile limit elongation of each material point in the fracture parameter matrix and the elongation of the bond. The bilinear damage function is:

[0037]

[0038] Among them, s p Represents the linear elastic limit elongation of each material point, s r represents the ultimate tensile elongation of each material point, ΔOFB represents the area of ​​the triangle formed by the unloading path after damage, the continuous loading path after damage and the ultimate tensile elongation in the damage fracture process, and ΔOAB represents the area of ​​the triangle formed by the loading and unloading path in the elastic stage, the continuous loading path after damage and the unloading path after damage in the damage fracture process.

[0039] Preferably, the bond-based peridynamic constitutive equation is:

[0040]

[0041] Where D is the bilinear damage function, ξ represents the relative position vector, η represents the relative displacement vector, δ is the domain radius, s represents the bond elongation, and c(0,δ) represents the viscous stiffness.

[0042] Preferably, based on the bond-based peridynamic constitutive equation, the volume information and the group information, the coordinate information of each material point is cyclically updated time by time step, and the elongation of the corresponding bond is updated based on the cyclically updated coordinate information to obtain the updated bond elongation of each material point, including:

[0043] Based on the bond-based peridynamic constitutive equation, calculating the acceleration of each of the material points at each time step, and calculating the velocity and displacement of each of the material points at each time step according to the acceleration;

[0044] The expression of acceleration at each time step is:

[0045]

[0046] where ρ represents the mass density, a(t) represents the acceleration of material point i at time step t, f represents the force function obtained by solving the bond-based peridynamic constitutive equation, and f(u(x j ,t)-u(x j ,t),x j -x i ) represents the force exerted by material point i on material point j at time step t, V j represents the volume of material point j, H i represents the population of material point i, b(x i ,t) represents the external force density of material point i at time step t;

[0047] The expression of the velocity at each time step is:

[0048] v(t)=v t-1 +a(t)dt

[0049] Among them, v t-1 represents the velocity at time step t-1;

[0050] The expression of displacement at each time step is:

[0051] u(t)=u t-1 +v(t)dt

[0052] Among them, u t-1 represents the velocity at the t-1 time step;

[0053] updating the coordinate information of each of the material points according to the velocity at each time step and the displacement at each time step;

[0054] Based on the updated coordinate information, the relative position vector is updated according to the relative position vector expression, and the relative displacement vector is updated according to the relative displacement vector expression;

[0055] Based on the updated relative position vector and the updated relative displacement vector, the elongation of the bond is updated according to the elongation expression of the bond to obtain the updated elongation of the bond of each material point.

[0056] Preferably, the updating of the group information based on the bilinear damage function value includes:

[0057] If the value of the bilinear damage function is not equal to 0, the group information is not updated;

[0058] If the value of the bilinear damage function is equal to 0, the two material points are mutually deleted from the population information.

[0059] In a second aspect, the present invention further provides a bond-based peridynamic numerical simulation system, the system comprising: a discrete parameter acquisition unit, a fracture parameter matrix acquisition unit, a loss function construction unit, a constitutive equation construction unit, a bond elongation update unit, a group information update unit and a damage judgment unit;

[0060] The discrete parameter acquisition unit is used to acquire the quantitative discrete information of the material to be analyzed during the material damage judgment process, wherein the quantitative discrete information includes the coordinate information, volume information and group information of each material point;

[0061] The fracture parameter matrix acquisition unit is used to randomly assign the linear elastic limit elongation and the tensile limit elongation of each material point using a Weibull function to obtain a fracture parameter matrix;

[0062] The loss function construction unit is used to calculate the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points in the group in sequence based on the coordinate information, and construct a bilinear damage function according to the fracture parameter matrix and the bond elongation;

[0063] The constitutive equation construction unit is used to construct a bond-based peridynamic constitutive equation according to the relative position vector, the relative displacement vector, the elongation of the bond, the bonding stiffness and the bilinear damage function;

[0064] The bond elongation updating unit is used to cyclically update the coordinate information of each material point in each time step within a preset total calculation time step based on the bond-based peridynamic constitutive equation, the volume information and the group information, and update the corresponding bond elongation based on the cyclically updated coordinate information to obtain the updated bond elongation of each material point;

[0065] The group information updating unit is used to substitute the updated key elongation into the bilinear damage function to obtain a bilinear damage function value, and update the group information based on the bilinear damage function value to obtain updated group information;

[0066] The damage judgment unit is used to calculate the local damage value of each material point according to the updated group information, obtain the damage and fracture state of the material to be analyzed at the corresponding time step according to the local damage value, and complete the calculation of all time steps in a cycle to obtain the bond-based peridynamic simulation result of the damage and fracture behavior of the material to be analyzed.

[0067] In a third aspect, the present invention also provides a computer device, comprising a memory, a processor and a transceiver, which are connected via a bus; the memory is used to store a set of computer program instructions and data, and transmit the stored data to the processor, and the processor executes the program instructions stored in the memory to execute the above-mentioned bond-based near-field dynamics numerical simulation method.

[0068] In a fourth aspect, the present invention further provides a computer-readable storage medium, wherein a computer program is stored in the computer-readable storage medium, and when the computer program is executed, the above-mentioned bond-based peridynamic numerical simulation method is implemented.

[0069] The present invention provides a method, system, device and medium for numerical simulation of bond-based peridynamics. Compared with the prior art, the beneficial effects of the embodiments of the present invention are at least one of the following:

[0070] (1) According to the characteristics of quasi-brittle heterogeneous materials, a method suitable for analyzing the fracture behavior of quasi-brittle heterogeneous materials is provided to improve the accuracy of fracture position prediction of cement-stabilized materials.

[0071] (2) The load-displacement curve obtained by the bond-based peridynamics numerical simulation method can predict the experimental results of the material to be analyzed and reflect the strain softening characteristics of the material to be analyzed. BRIEF DESCRIPTION OF THE DRAWINGS

[0072] Figure 1 It is a schematic diagram of the steps of a bond-based peridynamic numerical simulation method provided by a preferred embodiment of the present invention;

[0073] Figure 2 It is a schematic diagram of the steps of a quantization discretization method provided by a preferred embodiment of the present invention;

[0074] Figure 3 is a schematic diagram of the interaction between discrete material points provided by a preferred embodiment of the present invention;

[0075] Figure 4 It is a schematic diagram of the steps of a method for obtaining a fracture parameter matrix provided by a preferred embodiment of the present invention;

[0076] Figure 5 A preferred embodiment of the present invention provides a bilinear damage function diagram;

[0077] Figure 6 It is a schematic diagram of loading positions of a cement-stabilized gravel beam specimen in four-point bending provided by a preferred embodiment of the present invention;

[0078] Figure 7It is a comparison diagram of the damage cloud diagram of the cement stabilized gravel beam specimen at different time steps during the fracture process and the actual fracture photos of the material after the test provided by a preferred embodiment of the present invention;

[0079] Figure 8 It is a comparison diagram of the numerical results of the displacement-load curve applied in the test and the test results provided by a preferred embodiment of the present invention;

[0080] Fig. 9 It is a schematic diagram of the structure of a bond-based peridynamic numerical simulation system provided by a preferred embodiment of the present invention;

[0081] Fig.10 It is a schematic diagram of a computer device provided by a preferred embodiment of the present invention. DETAILED DESCRIPTION

[0082] The following is a detailed explanation of the implementation mode of the present invention in conjunction with the accompanying drawings. The embodiments are provided only for illustrative purposes and cannot be understood as limiting the present invention. The accompanying drawings are only for reference and illustration purposes and do not constitute a limitation on the scope of patent protection of the present invention. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention. In the description of the present invention, the terms "first", "second", "third", etc. are only used for descriptive purposes and cannot be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, the features defined as "first", "second", "third", etc. may explicitly or implicitly include one or more of the features. In the description of the present invention, unless otherwise specified, the meaning of "multiple" is two or more.

[0083] In the description of the present invention, it should be noted that, unless otherwise clearly specified and limited, the terms "installed", "connected" and "connected" should be understood in a broad sense. For example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be a connection between the two elements. The terms "vertical", "horizontal", "left", "right", "upper", "lower" and similar expressions used herein are only for illustrative purposes, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present invention. The term "and / or" used herein includes any and all combinations of one or more related listed items. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0084] In the description of the present invention, it should be noted that, unless otherwise defined, all technical and scientific terms used in the present invention have the same meaning as those commonly understood by those skilled in the art. The terms used in the specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood by specific circumstances.

[0085] See also Figure 1 In an embodiment of the present invention, a bond-based peridynamic numerical simulation method is provided, the method is applied to the simulation of fracture behavior of quasi-brittle inhomogeneous materials, the method comprising:

[0086] S1. In the process of material damage judgment, the quantitative discrete information of the material to be analyzed is obtained, and the quantitative discrete information includes the coordinate information, volume information and group information of each material point;.

[0087] S2. Using Weibull function to randomly assign the linear elastic limit elongation and the tensile limit elongation of each material point to obtain a fracture parameter matrix.

[0088] S3. Based on the coordinate information, the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points in the group are calculated in turn, and a bilinear damage function is constructed according to the fracture parameter matrix and the bond elongation.

[0089] S4. Constructing a bond-based peridynamic constitutive equation according to the relative position vector, the relative displacement vector, the elongation of the bond, the bonding stiffness and the bilinear damage function.

[0090] S5. Within a preset total calculation time step, based on the bond-based peridynamic constitutive equation, the volume information and the group information, the coordinate information of each of the material points is cyclically updated time step by time step, and the elongation of the corresponding bond is updated based on the cyclically updated coordinate information to obtain the updated bond elongation of each of the material points.

[0091] S6. Substitute the updated elongation of the key into the bilinear damage function to obtain a bilinear damage function value, and update the group information based on the bilinear damage function value to obtain updated group information.

[0092] S7. Calculate the local damage value of each material point according to the updated group information, obtain the damage and fracture state of the material to be analyzed at the corresponding time step according to the local damage value, and complete the calculation of all time steps in a cycle to obtain the bond-based peridynamic simulation result of the damage and fracture behavior of the material to be analyzed.

[0093] In a preferred embodiment of the present invention, a bond-based peridynamic numerical simulation method for cement-stabilized materials is provided. According to the characteristics of quasi-brittle heterogeneous materials, the material to be analyzed is discretized into a series of material points to obtain quantitative discrete information of the material to be analyzed. The quantitative discrete information includes coordinate information, volume information and group information of each material point, such as Figure 2 As shown, the following steps are included:

[0094] S101, obtaining values ​​of various mechanical parameters required for numerical simulation of bond-based peridynamics, wherein the various mechanical parameters at least include: elastic modulus, density, Poisson's ratio, geometric dimensions and material point size of the material to be analyzed.

[0095] S102, quantizing and discretizing the material to be analyzed according to the size of the material points to obtain coordinate information and volume information of each material point.

[0096] S103 . For each of the material points, take the neighboring material points within a preset area radius as a group of the material point, and obtain and store group information of the group.

[0097] In a preferred embodiment of the present application, the numerical values ​​of various mechanical parameters required for the numerical simulation of bond-based peridynamics are obtained, specifically including the material elastic modulus E, density ρ, Poisson's ratio υ, the geometric dimensions of the research object, the material point size dx, and the total calculation time step nt.

[0098] The material to be analyzed is discretized into a series of material points with a diameter of dx according to the material point size, and the coordinate information and volume information of each material point are obtained. Figure 3 The figure shows the interaction between discrete material points. i represents the neighboring material points of material point i in the research object, that is, the group of material point i, and the neighborhood radius is δ. The neighborhood radius δ is a positive number. In the preferred embodiment of the present application, the value is 3.015dx. There is an interaction force between material point i and its neighboring material points. Under three-dimensional conditions, H i The shape of is assumed to be spherical, whereas in two-dimensional plane stress and strain problems, H i The shape is assumed to be circular.

[0099] Cement stabilized materials are quasi-brittle heterogeneous materials, and the damage and fracture process of their bonds can be described by a bilinear damage function. In view of this characteristic of quasi-brittle heterogeneous materials, the linear elastic limit elongation and tensile limit elongation are used as the fracture parameters of the material to be analyzed. The Weibull function is used to randomly assign the linear elastic limit elongation and tensile limit elongation of each material point to obtain the fracture parameter matrix, as shown in Figure 4 As shown, the following steps are included:

[0100] S201, constructing a first matrix having the same number as the material points, and using a Weibull function to randomly assign values ​​to the first matrix.

[0101] S202. Obtain test values ​​of the linear elastic limit elongation and the tensile limit elongation of the material to be analyzed, and use the product of the randomly assigned first matrix and the test value of the linear elastic limit elongation as the linear elastic limit elongation matrix of each material point, and use the product of the randomly assigned first matrix and the test value of the tensile limit elongation as the tensile limit elongation matrix of each material point.

[0102] S203, using the linear elastic limit elongation matrix and the tensile limit elongation matrix as fracture parameter matrices.

[0103] The traditional bond-based peridynamic numerical simulation method assumes that the fracture parameters of all material points are the same, that is, the material is assumed to be uniform. However, cement-stabilized materials are a mixture of sand, gravel and cement, which is an inhomogeneous material. Therefore, random assignment of the fracture parameters of the material can reflect the inhomogeneous characteristics of quasi-brittle heterogeneous materials. In the preferred embodiment of the present application, the Weibull function is used to randomly assign the linear elastic limit elongation and tensile limit elongation of each material point to achieve the simulation of the quasi-brittle heterogeneous characteristics of cement-stabilized materials.

[0104] The Weibull function expression is:

[0105]

[0106] Where e represents the parameter that satisfies the Weibull function distribution, represents the scaling factor, which is related to the mean E(e) of e, but not exactly equal to it. β represents the shape parameter, which determines the shape of the distribution curve.

[0107] The expression of the mean value E(e) of the Weibull function parameter is:

[0108]

[0109] The mean of the Weibull function parameters is set to 1, and the shape parameter β = 4. Based on the number of material points, a first matrix A with the same size as the number of material points is constructed. The size of the first matrix A is [n×1], where n is the total number of material points formed after the material to be analyzed is discretized. The Weibull function is used to randomly assign values ​​to the first matrix.

[0110] Furthermore, for the linear elastic limit elongation and tensile limit elongation of each material point in the present application, first, the linear elastic limit elongation and tensile limit elongation of the material to be analyzed are measured by an experimental method, and the linear elastic limit elongation and tensile limit elongation of each material point of the material to be analyzed are used as the test value of the linear elastic limit elongation and tensile limit elongation of each material point. The product of the first matrix after random assignment and the test value of the linear elastic limit elongation is used as the linear elastic limit elongation matrix of each material point, that is, B1(i,1)=s p ×A(i,1). The product of the first matrix after random assignment and the test value of the tensile elongation is used as the tensile elongation matrix of each material point, that is, B2(i,1)=s r ×A(i,1).

[0111] Furthermore, boundary conditions and loading modes are set, wherein specific boundary conditions include fixed support and hinge support, etc., which are achieved by limiting the displacement of material points at the support position in different directions. Specific loading modes include displacement loading and external force loading, etc., which are achieved by applying displacement or external force to the material points at the loading position.

[0112] Furthermore, based on the set boundary conditions and loading mode, according to the coordinate information of each material point, the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points in the group are calculated in turn, and a bilinear damage function is constructed according to the fracture parameter matrix and the bond elongation. Specifically, according to the coordinate information, the relative position vector, relative displacement vector and bonding stiffness between any two material points at each time step are calculated. When the material to be analyzed is deformed from the initial state to the deformed state, during the entire deformation process, the two material points i and j will produce a displacement u i and u j , the material point coordinate changes to y i and j The relative position vector expression and relative displacement vector expression are as follows:

[0113] The relative position vector expression is:

[0114] ξ=x j -x i

[0115] Among them, x j represents the coordinates of material point j, x i represents the coordinates of material point i.

[0116] The relative displacement vector expression is:

[0117] η=u(x j ,t)-u(x i ,t)

[0118] Among them, u(xj ,t) represents the displacement vector of material point j at time step t, u(x i ,t) represents the displacement vector of material point i at time step t. The displacement vector u(x i ,t) is affected by the material point i and the neighborhood H i The influence of the interaction forces between other material points within.

[0119] The bond stiffness expression is:

[0120]

[0121] Where E is the elastic modulus and δ is the field radius. The bond stiffness is called the bond constant.

[0122] According to the relative position vector and the relative displacement vector, the elongation of the key is calculated. The expression of the elongation of the key is:

[0123]

[0124] A bilinear damage function is constructed based on the linear elastic limit elongation and tensile limit elongation of each material point in the fracture parameter matrix, as well as the elongation of the bond. The bilinear damage function is:

[0125]

[0126] Among them, s p Represents the linear elastic limit elongation of each material point, s r It represents the tensile limit elongation of each material point. Figure 5 The figure shows a bilinear damage function diagram. ΔOFB represents the area of ​​the triangle formed by the post-damage unloading path, the post-damage continuous loading path, and the tensile ultimate elongation during the bond damage fracture process, that is, Figure 5 The area of ​​the triangle OFB in the middle is ΔOAB, which is the area of ​​the triangle formed by the loading and unloading path in the elastic stage, the continuous loading path after damage, and the unloading path after damage during the damage-fracture process, that is, Figure 5 In the embodiment of the present invention, the bilinear damage function is used as the elongation damage fracture judgment criterion of the bond between material points.

[0127] The bond-based peridynamic constitutive equation is constructed based on the relative position vector, relative displacement vector, bond elongation, bonding stiffness and bilinear damage function. The bond-based peridynamic constitutive equation is:

[0128]

[0129] Where D is the bilinear damage function, ξ represents the relative position vector, η represents the relative displacement vector, δ is the domain radius, s represents the bond elongation, and c(0,δ) represents the viscous stiffness.

[0130] In a preferred embodiment of the present application, based on the bond-based peridynamic constitutive equation, volume information and group information, the coordinate information of each material point is cyclically updated time by time step, and the elongation of the corresponding bond is updated based on the cyclically updated coordinate information to obtain the updated bond elongation of each material point, including the following steps:

[0131] Based on the bond-based peridynamic constitutive equation, the acceleration of each material point at each time step is calculated, and the velocity and displacement of each material point at each time step are calculated based on the acceleration.

[0132] The expression of acceleration at each time step is:

[0133]

[0134] where a(t) represents the acceleration of material point i at time step t, ρ represents the mass density, f represents the force function obtained by solving the bond-based peridynamic constitutive equation, and f(u(x j ,t)-u(x j ,t),x j -x i ) represents the force exerted by material point i on material point j at time step t, V j represents the volume of material point j, H i represents the population of material point i, b(x i ,t) represents the external force density acting on material point i at time step t.

[0135] The expression of the velocity at each time step is:

[0136] v(t)=v t-1 +a(t)dt

[0137] Among them, v t-1 represents the velocity at time step t-1.

[0138] The expression of displacement at each time step is:

[0139] u(t)=u t-1 +v(t)dt

[0140] Among them, u t-1 represents the velocity at time step t-1.

[0141] The coordinate information of each material point is updated according to the velocity and displacement at each time step.

[0142] After substituting the updated coordinate information into the relative position vector expression and the relative displacement vector expression, the relative position vector and the relative displacement vector are updated.

[0143] The updated relative position vector and the updated relative displacement vector are substituted into the bond elongation to update the bond elongation, thereby obtaining the updated bond elongation of each material point.

[0144] Substitute the updated bond elongation into the bilinear damage function to obtain the bilinear damage function value. According to the bilinear damage function value, determine whether the elongation of the bond between each material point and the neighboring material point reaches the critical value. When the elongation of the bond between the material point and the neighboring material point is less than the linear elastic limit elongation s of the material point, the elongation of the bond between the material point and the neighboring material point is less than the linear elastic limit elongation s of the material point. p , the bilinear damage function value is 1. According to the bond-based peridynamic constitutive equation, the force between the two material points is calculated. It can be seen that when the bilinear damage function value is 1, the force between the two material points is not weakened, and no damage is caused. The elongation of the bond between the material point and the neighboring material point exceeds the linear elastic limit elongation s of the material point. p , but still less than the ultimate tensile elongation s r , the bilinear damage function value is between (0, 1), corresponding to the weakening of the force between the two material points, but no bond fracture occurs, and the bond is in a state of bond damage. The elongation of the bond between the material point and the neighboring material point exceeds the tensile limit elongation s of the material point. r , the value of the bilinear damage function is 0, the corresponding force between the two material points is 0, and the bond between the two material points is broken. After the bond is broken, the two material points need to delete each other from the group, update the group information, and obtain the updated group information.

[0145] In each iteration of the time step, the interaction forces between material points are recalculated and the group information of each material point is updated. The continuous accumulation of bond breakage eventually causes macroscopic damage to the material to be analyzed.

[0146] In a preferred embodiment of the present invention, after the group information is updated, the local damage value of each material point is calculated according to the updated group information. The local damage is defined as the weighted ratio of the number of neighboring material points whose interaction forces are eliminated in the material point group to the total number of initial neighboring material points in the material point group. The specific expression of the local damage value is:

[0147]

[0148] Where D(i,t) represents the number of neighboring material points whose interaction forces are eliminated in the population of material point i at time step t.

[0149] The calculated local damage value is a constant with a value between [0,1]. It indicates that all local bonds of the material to be analyzed are completely broken, that is, the material is completely broken.

[0150] The local damage value of each material point at each time step of the material to be analyzed is expressed in the form of a cloud diagram to obtain the overall damage and fracture conditions of the material to be analyzed at different time steps.

[0151] In the preferred embodiment of the present application, the four-point bending process of cement-stabilized crushed stone beam is taken as an example to verify the bond-based peridynamic numerical simulation method of the present invention. The dimensions of the four-point bending cement-stabilized crushed stone beam specimen are 400 mm long, 100 mm wide, and 100 mm high. The loading mode is displacement control, the loading speed is 0.1 mm / min, and the loading position is as follows: Figure 6 As shown. Three-dimensional analysis is used in the numerical simulation process. In order to reduce the number of material points and improve the calculation efficiency of the model, only the cement-stabilized gravel beam specimens subjected to bending moment within the dotted line are modeled. The calculation parameters are: material elastic modulus E = 4.5MPa, density ρ = 2000kg / m 3 , Poisson's ratio υ = 1 / 4, material point size dx = 0.005m, total calculation time step nt = 12000, dt = 1, the test value of the linear elastic limit elongation is 0.00011, and the test value of the tensile limit elongation is 0.00045. Figure 6 The test arrangement in the figure is a 2D picture of the front of the specimen, while the actual numerical simulation process is 3D modeling and calculation. Figure 7 The damage cloud diagram of the cement stabilized gravel beam specimen at different time steps during the fracture process is shown in the figure, and the actual fracture photos of the material after the test are compared. Figure 7 It can be seen that the bond-based peridynamic numerical simulation method of the present invention can accurately predict the fracture position and crack morphology of cement-stabilized gravel beam specimens. The bending and deflection of the cracks are related to the heterogeneity of the material. The use of Weibull function to randomly assign values ​​to the linear elastic limit elongation and tensile limit elongation of the material points makes the numerical simulation calculation results reflect the heterogeneity of the material. Figure 8 The figure shows the comparison between the numerical results and the experimental results of the displacement-load curve applied in the test. Figure 8 It can be seen that the displacement-load curve of the cement-stabilized gravel beam specimen predicted by the bond-based peridynamic numerical simulation method of the present invention fits the experimental data well, and the curve has a strain softening characteristic, which is consistent with the quasi-brittle characteristics of cement-stabilized materials.

[0152] In a preferred embodiment of the present invention, for the simulation of the fracture behavior of quasi-brittle inhomogeneous materials, in the process of material damage judgment, the quantitative discrete information of the material to be analyzed is obtained, and the quantitative discrete information includes the coordinate information, volume information and group information of each material point; the linear elastic limit elongation and the tensile limit elongation of each material point are randomly assigned using the Weibull function to obtain a fracture parameter matrix; based on the coordinate information, the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points in the group are calculated in turn, and a bilinear damage function is constructed according to the fracture parameter matrix and the bond elongation; the bond matrix is ​​constructed according to the relative position vector, relative displacement vector, bond elongation, bonding stiffness and bilinear damage function. Peridynamic constitutive equation; within the preset total calculation time step, based on the bond-based peridynamic constitutive equation, volume information and group information, the coordinate information of each material point is cyclically updated step by step, and the elongation of the corresponding bond is updated based on the cyclically updated coordinate information to obtain the updated bond elongation of each material point; the updated bond elongation is substituted into the bilinear damage function to obtain the bilinear damage function value, and the group information is updated based on the bilinear damage function value to obtain the updated group information; the local damage value of each material point is calculated according to the updated group information, and the bond-based peridynamic numerical simulation is performed on the material to be analyzed according to the local damage value to obtain the damage judgment result of the material to be analyzed. The bond-based peridynamic numerical simulation method provided by the present invention improves the accuracy of the fracture position prediction of cement-stabilized materials for the fracture behavior of quasi-brittle heterogeneous materials, and the load-displacement curve obtained by numerical simulation can predict the experimental results of the material and reflect the strain softening characteristics of the material.

[0153] Accordingly, if Fig. 9 As shown, based on a bond-based peridynamic numerical simulation method, an embodiment of the present invention further provides a bond-based peridynamic numerical simulation system, including: a discrete parameter acquisition unit 1, a fracture parameter matrix acquisition unit 2, a loss function construction unit 3, a constitutive equation construction unit 4, a bond elongation update unit 5, a group information update unit 6 and a damage judgment unit 7;

[0154] The discrete parameter acquisition unit 1 is used to acquire the quantitative discrete information of the material to be analyzed during the material damage judgment process, wherein the quantitative discrete information includes the coordinate information, volume information and group information of each material point;

[0155] The fracture parameter matrix acquisition unit 2 is used to randomly assign the linear elastic limit elongation and the tensile limit elongation of each material point using a Weibull function to obtain a fracture parameter matrix;

[0156] The loss function construction unit 3 is used to calculate the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points in the group in sequence based on the coordinate information, and construct a bilinear damage function according to the fracture parameter matrix and the bond elongation;

[0157] The constitutive equation construction unit 4 is used to construct a bond-based peridynamic constitutive equation according to the relative position vector, the relative displacement vector, the elongation of the bond, the bonding stiffness and the bilinear damage function;

[0158] The bond elongation updating unit 5 is used to cyclically update the coordinate information of each material point in each time step within a preset total calculation time step based on the bond-based peridynamic constitutive equation, the volume information and the group information, and update the corresponding bond elongation based on the cyclically updated coordinate information to obtain the updated bond elongation of each material point;

[0159] The group information updating unit 6 is used to substitute the updated key elongation into the bilinear damage function to obtain a bilinear damage function value, and based on the bilinear damage function value, update the group information to obtain updated group information;

[0160] The damage judgment unit 7 is used to calculate the local damage value of each material point according to the updated group information, obtain the damage and fracture state of the material to be analyzed at the corresponding time step according to the local damage value, and complete the calculation of all time steps in a cycle to obtain the bond-based peridynamic simulation result of the damage and fracture behavior of the material to be analyzed.

[0161] For the specific definition of a bond-based peridynamic numerical simulation system, please refer to the above-mentioned definition of a bond-based peridynamic numerical simulation method, which will not be repeated here. A person of ordinary skill in the art can appreciate that the various modules and steps described in conjunction with the embodiments disclosed in the present invention can be implemented in hardware, software, or a combination of both. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0162] like Fig.10 As shown, an embodiment of the present invention provides a computer device, including a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, the steps in the above-mentioned embodiment of the bond-based peridynamic numerical simulation method are implemented, for example Figure 1 Steps S1 to S7 described in .

[0163] Those skilled in the art will understand that the Fig.10 It is only an example of a computer device and does not constitute a limitation of the computer device. It may include more or fewer components than shown in the figure, or a combination of certain components, or different components. For example, the computer device may also include input and output devices, network access devices, buses, etc.

[0164] The processor may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc. The processor is the control center of the computer device, and uses various interfaces and lines to connect various parts of the entire computer device.

[0165] The memory can be used to store the computer program and / or module, and the processor realizes various functions of the computer device by running or executing the computer program and / or module stored in the memory, and calling the data stored in the memory. The memory can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, an application required for at least one function (such as a sound playback function, an image playback function, etc.), etc.; the data storage area can store data created according to the use of the mobile phone (such as audio data, a phone book, etc.), etc. In addition, the memory can include a high-speed random access memory, and can also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), at least one disk storage device, a flash memory device, or other volatile solid-state storage devices.

[0166] Wherein, if the module integrated in the computer device 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 present invention implements all or part of the processes in the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Wherein, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, U disk, mobile hard disk, disk, optical disk, computer memory, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), electrical carrier signal, telecommunication signal and software distribution medium, etc.

[0167] A person skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program, and the program can be stored in a computer-readable storage medium, and when the program is executed, it can include the processes of the embodiments of the above-mentioned methods. The storage medium can be a disk, an optical disk, a read-only memory (ROM) or a random access memory (RAM), etc.

[0168] Accordingly, an embodiment of the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium includes a stored computer program, wherein when the computer program is executed, the device where the computer-readable storage medium is located is controlled to perform the steps in the bond-based peridynamic numerical simulation method of the above embodiment, for example Figure 1 Steps S1 to S7 described in .

[0169] A bond-based peridynamic numerical simulation method, system, computer device and storage medium provided in this embodiment address the technical problem that existing bond-based peridynamic numerical simulation methods are not suitable for simulating the fracture behavior of quasi-brittle inhomogeneous materials. In the process of material damage judgment, the present invention obtains quantitative discrete information of the material to be analyzed, and the quantitative discrete information includes coordinate information, volume information and group information of each material point; the linear elastic limit elongation and tensile limit elongation of each material point are randomly assigned using Weibull function to obtain a fracture parameter matrix; based on the coordinate information, the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points in the group are calculated in turn, and a bilinear damage function is constructed according to the fracture parameter matrix and the bond elongation; a bond-based peridynamic constitutive equation is constructed according to the relative position vector, relative displacement vector, bond elongation, bonding stiffness and bilinear damage function; in the prediction In the total calculation time step set first, based on the bond-based peridynamic constitutive equation, volume information and group information, the coordinate information of each material point is cyclically updated step by step, and the elongation of the corresponding bond is updated based on the cyclically updated coordinate information to obtain the updated bond elongation of each material point; the updated bond elongation is substituted into the bilinear damage function to obtain the bilinear damage function value, and the group information is updated based on the bilinear damage function value to obtain the updated group information; the local damage value of each material point is calculated according to the updated group information, and the bond-based peridynamic numerical simulation is performed on the material to be analyzed according to the local damage value to obtain the damage judgment result of the material to be analyzed. The bond-based peridynamic numerical simulation method provided by the present invention improves the accuracy of the fracture position prediction of cement-stabilized materials for the fracture behavior of quasi-brittle heterogeneous materials, and the load-displacement curve obtained by numerical simulation can predict the experimental results of the material and reflect the strain softening characteristics of the material. .

[0170] The above-mentioned embodiments only express several preferred implementation modes of the present invention, and the description thereof is relatively specific and detailed, but it cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for ordinary technicians in this technical field, several improvements and substitutions can be made without departing from the technical principle of the present invention, and these improvements and substitutions should also be regarded as the protection scope of the present invention. Therefore, the protection scope of the patent of the present invention shall be based on the protection scope of the claims.

Claims

1. A bond-based peridynamic numerical simulation method, characterized in that: The method is applied to the simulation of fracture behavior of quasi-brittle inhomogeneous materials, and the method comprises: In the process of numerical simulation of material damage and fracture, quantitative discrete information of the material to be analyzed is obtained, wherein the quantitative discrete information includes coordinate information, volume information and group information of each material point; The linear elastic limit elongation and the tensile limit elongation of each material point are randomly assigned using a Weibull function to obtain a fracture parameter matrix; Based on the coordinate information, the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points in the group are calculated in sequence, and a bilinear damage function is constructed according to the fracture parameter matrix and the bond elongation; constructing a bond-based peridynamic constitutive equation according to the relative position vector, the relative displacement vector, the bond elongation, the bond stiffness and the bilinear damage function; Within a preset total calculation time step, boundary constraints and loading modes are set, and based on the bond-based peridynamic constitutive equation, the volume information and the group information, the coordinate information of each of the material points is cyclically updated time step by time step, and the elongation of the corresponding bond is updated based on the cyclically updated coordinate information to obtain the updated bond elongation of each of the material points; Substituting the updated bond elongation into the bilinear damage function to obtain a bilinear damage function value, and updating the group information based on the bilinear damage function value to obtain updated group information; The local damage value of each material point is calculated according to the updated group information, and the damage and fracture state of the material to be analyzed at the corresponding time step is obtained according to the local damage value. The calculation of all time steps is completed in a cycle to obtain the bond-based peridynamic simulation result of the damage and fracture behavior of the material to be analyzed.

2. The bond-based peridynamic numerical simulation method according to claim 1, characterized in that: The step of obtaining the quantitative discrete information of the material to be analyzed includes: Obtaining values ​​of various mechanical parameters required for numerical simulation of bond-based peridynamics, wherein the various mechanical parameters at least include: elastic modulus, density, Poisson's ratio, geometric dimensions and material point size of the material to be analyzed; Quantifying and discretizing the material to be analyzed according to the size of the material points to obtain coordinate information and volume information of each material point; For each of the material points, the neighboring material points within a preset area radius are taken as a group of the material points, and the group information of the group is acquired and stored.

3. The bond-based peridynamic numerical simulation method according to claim 2, characterized in that: The Weibull function is used to randomly assign the linear elastic limit elongation and the tensile limit elongation of each material point to obtain a fracture parameter matrix, including: Constructing a first matrix having the same size as the number of the material points, and randomly assigning values ​​to the first matrix using a Weibull function; Acquire test values ​​of the linear elastic limit elongation and the tensile limit elongation of the material to be analyzed, and use the product of the first matrix after random assignment and the test value of the linear elastic limit elongation as the linear elastic limit elongation matrix of each material point, and use the product of the first matrix after random assignment and the test value of the tensile limit elongation as the tensile limit elongation matrix of each material point; The linear elastic ultimate elongation matrix and the tensile ultimate elongation matrix are used as fracture parameter matrices.

4. The bond-based peridynamic numerical simulation method according to claim 3, characterized in that: Based on the coordinate information, the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points are calculated in sequence, and a bilinear damage function is constructed according to the fracture parameter matrix and the bond elongation, including: Calculating the relative position vector, relative displacement vector and bonding stiffness between any two material points according to the initial coordinate information; The relative position vector expression is: ξ=x j -x i Among them, x j represents the coordinates of material point j, x i represents the coordinates of material point i; The relative displacement vector expression is: η=u(x j ,t)-u(x i ,t) Among them, u(x j ,t) represents the displacement vector of material point j at time step t, u(x i ,t) represents the displacement vector of material point i at time step t; The bond stiffness expression is: Where, E represents the elastic modulus; δ is the field radius; According to the relative position vector and the relative displacement vector, the elongation of the key is calculated, and the expression of the elongation of the key is: A bilinear damage function is constructed according to the linear elastic limit elongation and the tensile limit elongation of each material point in the fracture parameter matrix and the elongation of the bond. The bilinear damage function is: Among them, s p Represents the linear elastic limit elongation of each material point, s r represents the ultimate tensile elongation of each material point, ΔOFB represents the area of ​​the triangle formed by the unloading path after damage, the continuous loading path after damage and the ultimate tensile elongation in the damage fracture process, and ΔOAB represents the area of ​​the triangle formed by the loading and unloading path in the elastic stage, the continuous loading path after damage and the unloading path after damage in the damage fracture process.

5. The bond-based peridynamic numerical simulation method according to claim 4, characterized in that: The bond-based peridynamic constitutive equation is: Where D is the bilinear damage function, ξ represents the relative position vector, η represents the relative displacement vector, δ is the domain radius, s represents the bond elongation, and c(0,δ) represents the viscous stiffness.

6. The bond-based peridynamic numerical simulation method according to claim 5, characterized in that: The step of cyclically updating the coordinate information of each material point in a time-step manner based on the bond-based peridynamic constitutive equation, the volume information, and the group information, and updating the elongation of the corresponding bond based on the cyclically updated coordinate information to obtain the updated bond elongation of each material point includes: Based on the bond-based peridynamic constitutive equation, calculating the acceleration of each of the material points at each time step, and calculating the velocity and displacement of each of the material points at each time step according to the acceleration; The expression of acceleration at each time step is: where a(t) represents the acceleration of material point i at time step t, ρ represents the mass density, f represents the force function obtained by solving the bond-based peridynamic constitutive equation, and f(u(x j ,t)-u(x j ,t),x j -x i ) represents the force exerted by material point i on material point j at time step t, V j represents the volume of material point j, H i represents the population of material point i, b(x i ,t) represents the external force density of material point i at time step t; The expression of the velocity at each time step is: v(t)=v t-1 +a(t)dt Among them, v t-1 represents the velocity at time step t-1; The expression of displacement at each time step is: u(t)=u t-1 +v(t)dt Among them, u t-1 represents the velocity at the t-1 time step; updating the coordinate information of each of the material points according to the velocity at each time step and the displacement at each time step; Based on the updated coordinate information, the relative position vector is updated according to the relative position vector expression, and the relative displacement vector is updated according to the relative displacement vector expression; Based on the updated relative position vector and the updated relative displacement vector, the elongation of the bond is updated according to the elongation expression of the bond to obtain the updated elongation of the bond of each material point.

7. The bond-based peridynamic numerical simulation method according to claim 1, characterized in that: The updating of the group information based on the bilinear damage function value includes: If the value of the bilinear damage function is not equal to 0, the group information is not updated; If the value of the bilinear damage function is equal to 0, the two material points are mutually deleted from the population information.

8. A bond-based peridynamic numerical simulation system, characterized in that: The system comprises: a discrete parameter acquisition unit, a fracture parameter matrix acquisition unit, a loss function construction unit, a constitutive equation construction unit, a bond elongation update unit, a group information update unit and a damage judgment unit; The discrete parameter acquisition unit is used to acquire the quantitative discrete information of the material to be analyzed during the material damage judgment process, wherein the quantitative discrete information includes the coordinate information, volume information and group information of each material point; The fracture parameter matrix acquisition unit is used to randomly assign the linear elastic limit elongation and the tensile limit elongation of each material point using a Weibull function to obtain a fracture parameter matrix; The loss function construction unit is used to calculate the relative position vector, relative displacement vector, bonding stiffness and bond elongation of any two material points in the group in sequence based on the coordinate information, and construct a bilinear damage function according to the fracture parameter matrix and the bond elongation; The constitutive equation construction unit is used to construct a bond-based peridynamic constitutive equation according to the relative position vector, the relative displacement vector, the elongation of the bond, the bonding stiffness and the bilinear damage function; The bond elongation updating unit is used to cyclically update the coordinate information of each material point in each time step within a preset total calculation time step based on the bond-based peridynamic constitutive equation, the volume information and the group information, and update the corresponding bond elongation based on the cyclically updated coordinate information to obtain the updated bond elongation of each material point; The group information updating unit is used to substitute the updated key elongation into the bilinear damage function to obtain a bilinear damage function value, and update the group information based on the bilinear damage function value to obtain updated group information; The damage judgment unit is used to calculate the local damage value of each material point according to the updated group information, obtain the damage and fracture state of the material to be analyzed at the corresponding time step according to the local damage value, and complete the calculation of all time steps in a cycle to obtain the bond-based peridynamic simulation result of the damage and fracture behavior of the material to be analyzed.

9. A computer device, characterized in that: The computer device includes a memory, a processor and a transceiver, which are connected via a bus; the memory is used to store a set of computer program instructions and data, and transmit the stored data to the processor, and the processor executes the program instructions stored in the memory to perform the bond-based peridynamics numerical simulation method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed, the bond-based peridynamics numerical simulation method according to any one of claims 1 to 7 is implemented.

Images