Crack calculation method and system based on parallel computing and matrix near-field dynamics

Abstract

The application belongs to the technical field of numerical simulation, and provides a crack calculation method and system based on parallel calculation and matrix near-field dynamics. The application firstly constructs a neighborhood matrix, a coordinate matrix, a displacement matrix, a load density matrix and a marking matrix; then, the force of a material point at each time step is calculated in combination with the coordinate matrix and the marking matrix; the marking matrix is used to filter the calculation result of the load density matrix, the resultant force of each material point is determined by summing each row of the load density matrix, and the deformation and displacement of the material point are solved by the resultant force of the material point; by converting all variables used in the near-field dynamics calculation process into a matrix format, based on the basic control equation of near-field dynamics, the CPU / GPU parallel matrix operation is used to realize efficient solution of near-field dynamics, and the problem of high time cost in the iterative calculation process by using the traditional method of successive search is avoided.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of numerical simulation, and particularly relates to a crack calculation method and system based on parallel computing and matrix peridynamics. BACKGROUND

[0002] Peridynamics is a method for establishing a model based on non-local action and describing physical and mechanical behavior of matter by solving a spatial integral equation, which avoids the singularity of traditional methods for modeling and solving spatial differential equations based on continuity assumption when facing discontinuous problems, and has significant advantages in material failure and crack propagation, and has become one of the most advanced numerical methods for solving discontinuous problems.

[0003] The inventor finds that, since peridynamics solves integral equations by considering non-local action between material points and solves physical and mechanical behavior of matter through a large number of micro time step cyclic iterations, compared with traditional methods for modeling and solving based on local action idea such as finite element method and discrete element method, the amount of calculation is greatly increased, and the calculation efficiency is greatly reduced with the increase of the number of material points, which limits the further application of peridynamics; at present, using CPU (central processing unit) or GPU (graphics processing unit) parallel computing is an effective way to improve the efficiency of numerical calculation, but the traditional peridynamic equation involves dynamic retrieval of material point fracture bonds and cyclic iteration of force, which is difficult to directly apply CPU or GPU parallel computing, and the time cost of iterative calculation process is high. SUMMARY

[0004] In order to solve the above problems, the application provides a crack calculation method and system based on parallel computing and matrix peridynamics, which stores all physical and mechanical information of material points in matrix format by establishing a specific material point neighborhood matrix, and uses CPU or GPU parallel computing to carry out pure matrix operation, so as to realize the acceleration of peridynamic solving process, thereby achieving the purposes of reducing time cost and improving calculation efficiency.

[0005] In order to achieve the above purposes, the application is implemented by the following technical scheme:

[0006] In the first aspect, the application provides a crack calculation method based on parallel computing and matrix peridynamics, which comprises:

[0007] determining a calculation space of an object to be solved, and discretizing the calculation space into a finite number of material points;

[0008] Number all the material points; build a neighborhood matrix according to the numbering sequence of the material points and the numbering of the adjacent material points that may be associated with each material point; build a coordinate matrix, a displacement matrix and a load density matrix of the near-field dynamics according to the neighborhood matrix; build a marker matrix representing the integrity of the material bond of the near-field dynamics;

[0009] According to the basic equation of the near-field dynamics, parallel matrix operations are used; in combination with the coordinate matrix and the marker matrix, the force on each material point at each time step is calculated; based on the displacement matrix, the scalar values in the marker matrix are judged, the marker matrix is used to filter the calculation results of the load density matrix, the resultant force on each material point is determined by summing each row of the load density matrix, and the deformation and displacement of the material point are solved by the resultant force on the material point.

[0010] Further, based on the updated displacement matrix, the scalar values in the marker matrix are judged, and in combination with the updated coordinate matrix and the marker matrix, the force on each material point at the next time step is solved; the near-field dynamics simulation is realized through cyclic iteration, and the calculation results at each time step are saved to complete the whole process of near-field dynamics numerical calculation.

[0011] Further, the material point refers to a node representing a certain volume and mass in the calculation region; each material point is the average of the physical and mechanical properties of the material in a certain space range.

[0012] Further, according to the numbering sequence of the material points, the numbering of the adjacent material points that may be associated with each material point is stored as a neighborhood matrix; the row number of the neighborhood matrix represents the center material point number, and the value of each row represents the adjacent material point number associated with the center material point; the maximum number of adjacent material points is taken as the number of matrix columns, and 0 is added to the positions without adjacent material point numbers;

[0013] The coordinate matrix refers to the relative coordinate values of the center material point and its adjacent material points, the row number of the coordinate matrix represents the center material point number, and the value of each position is the relative coordinate value of the material point represented by the corresponding position of the neighborhood matrix and the center material point;

[0014] The displacement matrix refers to the relative displacement values of the center material point and its adjacent material points, the row number of the displacement matrix represents the center material point number, and the value of each position is the relative displacement value of the material point represented by the corresponding position of the neighborhood matrix and the center material point;

[0015] The load density matrix refers to the interaction force between the center material point and its adjacent material points, the row number of the load density matrix represents the center material point number, and the value of each position is the interaction force between the material point represented by the corresponding position of the neighborhood matrix and the center material point.

[0016] Further, the row number of the marker matrix represents the center material point number, and the value of each row of the marker matrix represents the integrity state of the material bond between the center material point and the adjacent material point; the marker matrix has the same size as the neighborhood matrix, has the same corresponding material point number sequence, and has no adjacent material point number position, and 0 is supplemented.

[0017] The material bond refers to the interaction relationship between the center material point and other material points within the neighborhood range of the center material point.

[0018] Further, the basic equation of near-field dynamics is:

[0019]

[0020] Wherein, [p] is a material point mass density matrix; is a material point acceleration matrix; [F] is a material point load density matrix; [V i ] is a volume matrix of other adjacent material points in the neighborhood range of the center material point; [b] is a material point body force density matrix; n represents the maximum value of the number of adjacent material points.

[0021] Further, the load density matrix calculation refers to solving the force of the near-field dynamics material point according to the material point relative coordinate matrix and the marker matrix.

[0022]

[0023] Wherein, [F] is a material point load density matrix; [B] is a near-field dynamics marker matrix; [K] is a bulk modulus matrix; [G] is a shear modulus matrix; [m] is a weighted volume scalar coefficient matrix; [w] is an influence function matrix; [X] is a material point coordinate matrix; [theta] is a material point volume expansion rate matrix; [e d ] is a bias tensor matrix of the elongation state of the bond.

[0024] In a second aspect, the present application further provides a crack calculation system based on parallel computing and matrix near-field dynamics, comprising:

[0025] The discrete module is configured to determine the calculation space of the object to be solved, and discretize the calculation space into a finite number of material points.

[0026] The matrix construction module is configured to number all the material points, construct a neighborhood matrix according to the material point number sequence and the adjacent material point numbers that may exist in connection with each material point, construct a coordinate matrix, a displacement matrix and a load density matrix of near-field dynamics according to the neighborhood matrix, and construct a marker matrix representing the integrity of the material bond of near-field dynamics.

[0027] The computing module is configured to: according to a near-field dynamics basic equation, utilize parallel matrix operation; in combination with a coordinate matrix and a marker matrix, calculate the force of each material point at each time step; based on a displacement matrix, judge the scalar value in the marker matrix, filter the calculation result of the load density matrix by using the marker matrix, determine the resultant force of each material point by summing each row of the load density matrix, and solve the deformation and displacement of the material point by the resultant force of the material point.

[0028] In a third aspect, the present application further provides a computer readable storage medium, which stores a computer program, and the program is executed by a processor to implement the steps of the crack calculation method based on parallel computing and matrix near-field dynamics in the first aspect.

[0029] In a fourth aspect, the present application further provides an electronic device, which comprises a memory, a processor and a computer program stored in the memory and executable on the processor, and the processor implements the steps of the crack calculation method based on parallel computing and matrix near-field dynamics in the first aspect when executing the program.

[0030] Compared with the prior art, the present application has the following beneficial effects:

[0031] 1. Firstly, the present application constructs a neighborhood matrix, a coordinate matrix, a displacement matrix, a load density matrix and a marker matrix; then according to a near-field dynamics basic equation, utilizes parallel matrix operation; in combination with the coordinate matrix and the marker matrix, calculates the force of each material point at each time step; based on the displacement matrix, judges the scalar value in the marker matrix, filters the calculation result of the load density matrix by using the marker matrix, determines the resultant force of each material point by summing each row of the load density matrix, and solves the deformation and displacement of the material point by the resultant force of the material point; by converting all variables used in the near-field dynamics calculation process into a matrix format, utilizing CPU / GPU parallel matrix operation based on the near-field dynamics basic control equation, the efficient solution of the near-field dynamics is realized, and the problem of high time cost in the iterative calculation process by using the traditional method of successive search is avoided.

[0032] 2. The present application is not only suitable for CPU parallel, but also suitable for GPU parallel, which significantly improves the calculation efficiency and can realize efficient simulation of millions of material points on a single computer. BRIEF DESCRIPTION OF DRAWINGS

[0033] The drawings constituting a part of this embodiment are used to provide further understanding of the embodiment, and the schematic embodiment and its description are used to explain the embodiment, and do not constitute improper limitation on the embodiment.

[0034] Figure 1 It is a flowchart of the embodiment 1 of the present application;

[0035] Figure 2 Fig. 1 is a schematic diagram of a coordinate matrix of an original object in an embodiment of the present application;

[0036] Figure 3 Fig. 2 is a schematic diagram of a neighborhood matrix of the original object in the embodiment of the present application;

[0037] Figure 4 Fig. 3 is a schematic diagram of a relative coordinate component in an x direction of the original object in the embodiment of the present application;

[0038] Figure 5 Fig. 4 is a schematic diagram of a relative displacement component in the x direction of the original object in the embodiment of the present application;

[0039] Figure 6 Fig. 5 is a schematic diagram of a load density matrix in the x direction of the original object in the embodiment of the present application and a resultant force thereof;

[0040] Figure 7 Fig. 6 is a schematic diagram of a marker matrix of the original object in the embodiment of the present application. DETAILED DESCRIPTION

[0041] The present application will be further described below in conjunction with the accompanying drawings and embodiments.

[0042] It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the present application. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the present application belongs.

[0043] Embodiment 1

[0044] The present embodiment provides a crack calculation method based on parallel computing and matrix near-field dynamics, comprising:

[0045] determining a calculation space of an object to be solved, and discretizing the calculation space into a finite number of material points; optionally, the calculation space is a region or space on a material to be calculated, that is, the object to be solved can be concrete, steel and other materials that can have cracks or models corresponding to the materials;

[0046] numbering all the material points; constructing a neighborhood matrix according to the numbering order of the material points and the numbering of the neighboring material points that each material point can be connected to; constructing a coordinate matrix, a displacement matrix and a load density matrix of near-field dynamics according to the neighborhood matrix; constructing a marker matrix representing the integrity of the material bond of near-field dynamics; wherein the coordinate matrix can also be referred to as a near-field dynamics coordinate matrix or a material point coordinate matrix, the displacement matrix can also be referred to as a near-field dynamics displacement matrix or a material point displacement matrix, the load density matrix can also be referred to as a near-field dynamics load density matrix or a material point load density matrix, and the marker matrix can also be referred to as a near-field dynamics marker matrix or a material point marker matrix;

[0047] According to the basic equation of near-field dynamics, parallel matrix operation is used; in combination with the coordinate matrix and the marker matrix, the force of the material point at each time step is calculated; based on the displacement matrix, the scalar value in the marker matrix is determined, the marker matrix is used to filter the calculation results of the load density matrix, the sum of each row of the load density matrix is determined to determine the resultant force of each material point, and the deformation and displacement of the material point are solved through the resultant force of the material point.

[0048] Specifically, the calculation space is discretized into a finite number of material points, and all the material points are numbered from 1 to m. According to the numbering order of the material points, the coordinate information of all the material points is stored in the original coordinate matrix [X I ] in sequence; the near-field dynamics neighborhood matrix is constructed by material point retrieval, and the numbering of the adjacent material points that may be connected with each material point is stored in the near-field dynamics neighborhood matrix [N] according to the numbering order of the material points; the matrix form of the near-field dynamics coordinate matrix [X], the near-field dynamics displacement matrix [U], the near-field dynamics load density matrix [F] and other physical and mechanical parameter matrices are constructed according to the near-field dynamics neighborhood matrix; the near-field dynamics marker matrix [B] representing the integrity of the near-field dynamics bond is constructed, wherein the row number of the near-field dynamics marker matrix represents the center material point number, and the value stored in each row of the matrix represents the integrity state of the bond between the center material point and the adjacent material point, wherein 1 represents that the bond between the material points is complete, and 0 represents that the bond between the material points is broken; according to the basic equation of near-field dynamics, parallel matrix operation is used to calculate the force of the material point at each time step, the marker matrix is used to filter the calculation results of the load density matrix, the sum of each row of the load density matrix is determined to determine the resultant force of each center material point, the deformation and displacement of the material point are solved through the resultant force of the material point, the scalar value in the near-field dynamics marker matrix is determined based on the updated material point displacement matrix, the force of the material point at the next time step is solved in combination with the updated material point coordinate matrix and the marker matrix; the near-field dynamics simulation is realized through such a cycle iteration, and the calculation results at each time step are saved, the calculation results are output and viewed to complete the whole process of near-field dynamics numerical calculation. The specific steps are as follows:

[0049] S1, model discretization:

[0050] Specifically, the calculation space is discretized into a finite number of material points, each material point contains certain physical and mechanical information such as volume and mass, and all the material points are numbered from 1 to m; according to the numbering order of the material points, the coordinate information of all the material points is stored in the original coordinate matrix [X I ] in sequence, and the original coordinate matrix has a size of m×3, such as Figure 2As shown in the formula (1), wherein the row number of the matrix represents the particle number, and the column number of the matrix represents the particle x coordinate, the particle y coordinate and the particle z coordinate respectively.

[0051] Optionally, the particle refers to a node representing a certain volume, mass and other physical and mechanical information in a calculation region; the particle itself does not occupy a spatial volume, and only the material volume and mass in the space range represented by the particle are recorded in the material parameter matrix, so each particle can be regarded as the average of the physical and mechanical properties of the material in a certain space range.

[0052] Optionally, the original coordinate matrix [X I ] refers to a two-dimensional matrix stored from top to bottom according to the particle number sequence, and the matrix size is m*3, wherein m represents the total number of particles, the row number of the matrix corresponds to the particle number, and the particle x coordinate, the particle y coordinate and the particle z coordinate are stored from left to right in each row of the matrix. For a two-dimensional plane problem, the particle z coordinate can be replaced by 0, and the following is similar.

[0053] S2, constructing a near-field dynamics neighborhood matrix:

[0054] Specifically, the near-field dynamics neighborhood matrix is constructed by particle searching, wherein each particle is in contact with all other particles within a certain distance or neighborhood range; according to the particle number sequence, the neighborhood numbers of the adjacent particles that may be in contact with each particle are stored as a neighborhood matrix [N], and the neighborhood matrix size is m*n, as shown in the formula (2). Figure 3 As shown in the formula (2), wherein the row number of the matrix represents the center particle number, and the value stored in each row of the matrix represents the neighborhood number of the adjacent particle in contact with the center particle; due to the different coordinate positions of the particles, not all particles have the same number of adjacent particles, but in order to ensure the integrity of the matrix, the maximum number of adjacent particles n is taken as the column number of the matrix, and 0 is supplemented to the position without the neighborhood number of the particle; after the neighborhood matrix is determined, it will not change with subsequent calculation.

[0055] Optionally, the neighborhood range refers to the maximum distance between the center particle and other particles that produce interaction force with the center particle, which can be expressed as: ||x i -x j ||≤δ, wherein x i and x j represent the coordinates of the particles numbered i and j respectively, and δ represents the neighborhood range.

[0056] Optionally, the neighborhood matrix [N] refers to a two-dimensional matrix storing the number of each material point and its adjacent material points, with a size of m*n, wherein the row number of the matrix represents the number of the center material point, and the value stored in each column corresponds to the number of the adjacent material point in contact with the center material point. For each material point, there is a certain number of adjacent material points in contact with it. However, due to the different positions of the material points, the number of adjacent material points of the material points is different, so the number of matrix rows is set to the maximum value, and the positions without the number of adjacent material points of the material points are filled with 0, so as to ensure that the matrix can normally participate in the calculation.

[0057] S3, initializing the near-field dynamics parameter matrix:

[0058] Specifically, according to the near-field dynamics neighborhood matrix, the near-field dynamics coordinate matrix [X], the near-field dynamics displacement matrix [U], the near-field dynamics load density matrix [F], and other physical and mechanical parameter matrices are constructed in the form of a matrix, wherein the size of each matrix is m*n, the row number of the matrix represents the number of the center material point, and the value stored in each row of the matrix represents the information of the number of the adjacent material point, and the positions without the information of the material point are filled with 0.

[0059] Optionally, the near-field dynamics coordinate matrix, as shown in Figure 4 , refers to the relative coordinate values of the center material point and its adjacent material points, including an x-direction coordinate component matrix, a y-direction coordinate component matrix, and a z-direction coordinate component matrix, and the size of each matrix is the same as that of the neighborhood matrix, i.e., m*n. Among them, the row number of the matrix represents the number of the center material point, and the value stored in each position is the relative coordinate value of the material point represented by the corresponding position of the neighborhood matrix and the center material point.

[0060] Optionally, the near-field dynamics displacement matrix, as shown in Figure 5 , refers to the relative displacement values of the center material point and its adjacent material points, including an x-direction displacement component matrix, a y-direction displacement component matrix, and a z-direction displacement component matrix, and the size of each matrix is the same as that of the neighborhood matrix, i.e., m*n. Among them, the row number of the matrix represents the number of the center material point, and the value stored in each position is the relative displacement value of the material point represented by the corresponding position of the neighborhood matrix and the center material point.

[0061] Optionally, the near-field dynamics load density matrix [F], as shown in Figure 6 , refers to the interaction force between the center material point and its adjacent material points, including an x-direction load density component matrix, a y-direction load density component matrix, and a z-direction load density component matrix, and the size of each matrix is the same as that of the neighborhood matrix, i.e., m*n. Among them, the row number of the matrix represents the number of the center material point, and the value stored in each position is the interaction force between the material point represented by the corresponding position of the neighborhood matrix and the center material point.

[0062] Optionally, other physical and mechanical parameter matrices mainly refer to various parameters required in the near-field dynamics calculation process, including but not limited to surface correction coefficients, volume correction coefficients, elastic modulus, Poisson's ratio, density, volume weight and volume strain, etc. The size of each matrix is the same as that of the neighborhood matrix, i.e. m x n. Among them, the row number of the matrix represents the center material point number, and the value stored at each position is the physical and mechanical parameter of the material point represented by the corresponding position of the neighborhood matrix.

[0063] S4, constructing a near-field dynamics marker matrix:

[0064] Specifically, a marker matrix representing the integrity of the near-field dynamics bond is constructed according to the neighborhood matrix. The near-field dynamics represents the integrity of the bond between material points by introducing a scalar function μ with historical dependence, μ = 1 or 0, where μ = 1 represents that the bond between material points is complete, and μ = 0 represents that the bond between material points is broken; force transmission can occur between complete bonds, while force transmission no longer occurs between broken bonds; the near-field dynamics marker matrix [B] is constructed in sequence according to the material point number, and the size of the near-field dynamics marker matrix is m x n, as shown in the following formula: Figure 7 wherein the row number of the matrix represents the center material point number, and the value stored in each row of the matrix represents the integrity state of the bond between the center material point and the adjacent material point; the marker matrix and the neighborhood matrix are of the same size, and the order of the corresponding material point numbers is the same, wherein the positions without adjacent material point numbers are supplemented with 0; the near-field dynamics marker matrix needs to be dynamically updated according to the calculation process.

[0065] Optionally, the near-field dynamics material bond refers to the interaction relationship between the center material point and other material points within its neighborhood range, which is called bond, and there is no interaction between the center material point and other material points outside its neighborhood range, i.e. no bond.

[0066] Optionally, the near-field dynamics marker matrix [B] refers to a scalar function representing the integrity of the material point bond, and its value has only two states, 1 or 0, wherein 1 represents that the interaction between the current material point and the center material point exists, and 0 represents that the interaction between the current material point and the center material point no longer exists.

[0067] S5, carrying out parallel calculation:

[0068] Specifically, according to the basic equation of near-field dynamics, the force of the material point at each time step is calculated by using CPU / GPU parallel matrix operation, that is, the load density matrix is calculated, the load density matrix is filtered by using the marker matrix, the resultant force of each central material point is determined by summing each row of the load density matrix, and the actual force of the material point is obtained; the deformation and displacement of the material point are solved by the resultant force of the material point, the scalar value in the near-field dynamics marker matrix [B] is determined based on the updated material point displacement matrix, and the force of the material point at the next time step is solved in combination with the updated material point coordinate matrix and the marker matrix.

[0069] Optionally, the basic equation of near-field dynamics refers to a control equation describing the relationship between the displacement and force of the material point, which is usually expressed in the form of Newton's second law, and has the form of:

[0070]

[0071] Where [p] is the mass density matrix of the material point; is the acceleration matrix of the material point, which is the second derivative of the displacement matrix [U] with respect to time t and the first derivative of the velocity matrix with respect to time t; [F] is the load density matrix of the material point; [V i ] is the volume matrix of other adjacent material points in the neighborhood range of the central material point; [b] is the body force density matrix of the material point; and n represents the maximum value of the number of adjacent material points.

[0072] Optionally, the parallel calculation refers to parallel calculation of the matrix operation process by using the CPU (central processing unit) or GPU (graphics processing unit) of the computer to accelerate the matrix operation process.

[0073] Optionally, the load density matrix calculation refers to solving the force of the near-field dynamics material point according to the relative coordinate matrix and the marker matrix of the material point:

[0074]

[0075] Where [F] is the load density matrix of the material point; [B] is the near-field dynamics marker matrix; [K] is the bulk modulus matrix; [G] is the shear modulus matrix; [m] is the weighted volume scalar coefficient matrix; [w] is the influence function matrix; [X] is the material point coordinate matrix; [theta] is the material point volume expansion rate matrix; and [e d is the bias tensor matrix of the elongation state of the bond.

[0076] Optionally, the bias tensor matrix [e d of the elongation state of the bond has the form of:

[0077]

[0078] where [e] is the elongation state matrix of the material point bond, and its expression is

[0079] [e] = ||[X] + [U]|| - ||X||

[0080] where [X] is the coordinate matrix of the material point, and [U] is the displacement matrix of the material point.

[0081] Optionally, the weighted volume scalar coefficient matrix [m] is shown as follows:

[0082] [m] = [ω][X] · [X]

[0083] where · represents convolution.

[0084] Optionally, the influence function matrix [ω] refers to the correction of the interaction force between adjacent material points on the central material point, which usually includes a linear influence function (ω = 1), an inverse influence function (ω = 1 / r), etc., where r represents the relative distance between the material points.

[0085] Optionally, the volume expansion rate matrix [θ] refers to the volume change rate of the material point under the current stress state, also known as the volume strain, and its expression is:

[0086]

[0087] where · represents convolution.

[0088] Optionally, the load resultant force refers to filtering the calculation results of the load density matrix [F] using the marker matrix, and determining the resultant force [F s ] on each central material point by summing each row of the load density matrix.

[0089] Optionally, the material point deformation and displacement solution refers to first solving the acceleration matrix of the material point using the basic equation of near-field dynamics, and then solving the velocity matrix of the material point using time difference, and its expression is as follows:

[0090]

[0091] Further, the displacement matrix [U] of the material point is solved using time difference, and its expression is as follows:

[0092]

[0093] where Δt is the time step; and l is the lth time step.

[0094] Optionally, the marker matrix update refers to judging whether the fracture condition is met according to the elongation of the near-field dynamic material point bond. In general, the fracture condition of the material point bond is expressed as follows:

[0095]

[0096] Wherein, s0 is the critical elongation of the material point bond, i.e. the maximum allowable elongation, which is determined by the material properties, and s is the elongation of the material point bond, which is expressed as:

[0097]

[0098] That is, when the elongation of the material point bond exceeds the critical elongation, the bond will be permanently fractured, at which time μ = 0. If the critical elongation is not exceeded, the bond remains intact, at which time μ = 1.

[0099] S6, time step loop iteration calculation:

[0100] Specifically, through loop iteration, the material point load density matrix-displacement matrix-load density matrix is solved in turn, the near-field dynamic simulation is realized, and the calculation results of each time step are saved. Through output and viewing of the calculation results, the whole process of near-field dynamic numerical calculation is completed.

[0101] Optionally, the loop iteration refers to dividing the calculation time into a finite number of time steps, calculating the material point load density matrix of the current time step, solving the material point displacement matrix, and then solving the material point load density matrix of the next time step according to the current material point deformation, i.e. through such loop iteration to realize the near-field dynamic solution.

[0102] Optionally, the result output and viewing refers to saving the calculation results of each time step, including material point coordinates, displacement and damage parameters, to a specified file path, and viewing the calculation results by extracting relevant data.

[0103] Embodiment 2:

[0104] The embodiment provides a crack calculation system based on parallel computing and matrix near-field dynamics, comprising:

[0105] The discrete module is configured to determine the calculation space of the object to be solved, and discretize the calculation space into a finite number of material points;

[0106] The matrix construction module is configured to number all the material points, construct a neighborhood matrix according to the material point numbering sequence and the numbering of the adjacent material points that may be associated with each material point, construct a coordinate matrix, a displacement matrix and a load density matrix of the near-field dynamics according to the neighborhood matrix, and construct a marker matrix representing the integrity of the near-field dynamic material bond;

[0107] The computing module is configured to: according to the near-field dynamics basic equation, utilize parallel matrix operation; combine the coordinate matrix and the marker matrix to calculate the force of the material point at each time step; based on the displacement matrix, judge the scalar value in the marker matrix, filter the load density matrix calculation result by using the marker matrix, determine the resultant force of each material point by summing each row of the load density matrix, and solve the deformation and displacement of the material point by the resultant force of the material point.

[0108] The working method of the system is the same as the crack calculation method based on parallel computing and matrix near-field dynamics of embodiment 1, and details are not repeated here.

[0109] Embodiment 3:

[0110] The embodiment provides a computer readable storage medium, which stores a computer program, and the program is executed by a processor to realize the steps of the crack calculation method based on parallel computing and matrix near-field dynamics of embodiment 1.

[0111] Embodiment 4:

[0112] The embodiment provides an electronic device, which comprises a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and the processor realizes the steps of the crack calculation method based on parallel computing and matrix near-field dynamics of embodiment 1 when executing the program.

[0113] The above only describes the preferred embodiments of the present embodiment and is not used to limit the present embodiment. The present embodiment can have various changes and variations for those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present embodiment shall be included in the protection scope of the present embodiment.

Claims

1. A crack calculation method based on parallel computation and matrix near-field dynamics, characterized by, The application relates to a method for simulating the dynamic process of a structure by using near-field dynamics. The method comprises the following steps: determining a calculation space of an object to be solved, and discretizing the calculation space into a finite number of material points; the object to be solved is concrete and steel; numbering all the material points, and constructing a neighborhood matrix according to the numbering sequence of the material points and the numbering of the neighboring material points associated with each material point; constructing a coordinate matrix, a displacement matrix and a load density matrix of near-field dynamics according to the neighborhood matrix; constructing a marker matrix representing the integrity of the material bond of near-field dynamics; calculating the force of each material point at each time step by using parallel matrix operation based on the basic equation of near-field dynamics, combining the coordinate matrix and the marker matrix; 2. The parallel computing and matrix-based near-field dynamics based crack calculation method of claim 1, wherein, judging the scalar value in the marker matrix based on the displacement matrix, filtering the calculation result of the load density matrix by using the marker matrix, determining the resultant force of each material point by summing each row of the load density matrix, and solving the deformation and displacement of the material point by using the resultant force of the material point.

3. The parallel computing and matrix-based near-field dynamics based crack calculation method of claim 1, wherein, judging the scalar value in the marker matrix based on the updated displacement matrix, solving the force of the material point at the next time step by combining the updated coordinate matrix and the marker matrix, realizing the near-field dynamics simulation by cyclic iteration, saving the calculation result of each time step, and completing the whole process of the near-field dynamics numerical calculation.

4. The parallel computing and matrix-based near-field dynamics based crack calculation method of claim 1, wherein, The material point refers to a node representing a certain volume and mass in a calculation region; each material point is the average of the physical and mechanical properties of the material in a certain space range. The numbering of the neighboring material points associated with each material point is stored as a neighborhood matrix according to the numbering sequence of the material points; the row number of the neighborhood matrix represents the numbering of the center material point, and the value of each row represents the numbering of the neighboring material points associated with the center material point; the maximum number of the neighboring material points is taken as the number of the matrix columns, and 0 is supplemented to the positions without the numbering of the neighboring material points. The coordinate matrix refers to the relative coordinate values of the center material point and its neighboring material points; the row number of the coordinate matrix represents the numbering of the center material point, and the value of each position is the relative coordinate value of the material point represented by the corresponding position of the neighborhood matrix and the center material point. The displacement matrix refers to the relative displacement values of the center material point and its neighboring material points; the row number of the displacement matrix represents the numbering of the center material point, and the value of each position is the relative displacement value of the material point represented by the corresponding position of the neighborhood matrix and the center material point.

5. The parallel computing and matrix-based near-field dynamics based crack calculation method of claim 1, wherein, The load density matrix refers to the interaction force of the center material point and its neighboring material points; the row number of the load density matrix represents the numbering of the center material point, and the value of each position is the interaction force of the material point represented by the corresponding position of the neighborhood matrix and the center material point. The row number of the marker matrix represents the numbering of the center material point, and the value of each row of the marker matrix represents the integrity state of the material bond between the center material point and the neighboring material points; the marker matrix and the neighborhood matrix are of the same size, and the numbering sequence of the corresponding material points is the same; 0 is supplemented to the positions without the numbering of the neighboring material points.

6. The parallel computing and matrix-based near-field dynamics based crack calculation method of claim 1, wherein, The material bond refers to the interaction relationship between the center material point and other material points within the neighborhood range of the center material point. The basic equation of near-field dynamics is as follows: where, is the mass density matrix of the material points; is the acceleration matrix of the material points; is the load density matrix of the material points; is the volume matrix of other neighboring material points within the neighborhood range of the central material point; is the body force density matrix of the material points; n represents the maximum value of the number of neighboring material points.

7. The parallel computing and matrix-based near-field dynamics based crack calculation method of claim 1, wherein, The load density matrix calculation refers to solving forces of the material points in the near-field dynamics according to the relative coordinate matrix and the marker matrix of the material points: wherein, is a material point load density matrix; is a near-field dynamics marker matrix; is a bulk modulus matrix; is a shear modulus matrix; is a weighted bulk scalar coefficient matrix; is an influence function matrix; is a material point coordinate matrix; is a material point volumetric dilatation matrix; is a bias tensor matrix of the elongation state of the bond.

8. A crack calculation system based on parallel computation and matrix near field dynamics, characterized by, The method comprises the steps of: a discretization module configured to determine a calculation space of a to-be-solved object, and discretize the calculation space into a finite number of material points; the to-be-solved object is concrete and steel; a matrix construction module configured to number all the material points, and construct a neighborhood matrix according to a numbering sequence of the material points and numbering of adjacent material points associated with each material point; a coordinate matrix, a displacement matrix and a load density matrix of the near-field dynamics are constructed according to the neighborhood matrix; a marker matrix representing integrity of a material bond in the near-field dynamics is constructed; a calculation module configured to utilize parallel matrix operation according to a basic equation of the near-field dynamics, and calculate forces of the material points at each time step in combination with the coordinate matrix and the marker matrix; scalar values in the marker matrix are judged based on the displacement matrix, the marker matrix is used to filter a calculation result of the load density matrix, a resultant force of each material point is determined by summing each row of the load density matrix, and deformation and displacement of the material points are solved based on the resultant force of the material points.

9. A computer readable storage medium having stored thereon a computer program, characterized in that, The program, when executed by the processor, implements the steps of the crack calculation method based on parallel calculation and matrix near-field dynamics according to any one of claims 1-7.

10. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, The processor, when executing the program, implements the steps of the crack calculation method based on parallel calculation and matrix near-field dynamics according to any one of claims 1-7.

Images