Method, electronic device and storage medium for establishing high volume ratio target polyhedral aggregate model

Abstract

The invention discloses a method, an electronic device and a storage medium for establishing a target polyhedral aggregate model with a high volume rate, and belongs to the technical field of establishing a finite element numerical model of concrete. The invention provides a polyhedral random aggregate model with a high volume rate. First, according to the size and the aggregate volume fraction of the target polyhedral aggregate model, the basic polyhedral aggregate model parameters with a larger size and a lower aggregate volume fraction are calculated, and a transition aggregate with a slightly larger particle size is generated and placed; then, the material parameters, mesh parameters and load parameters of the basic polyhedral aggregate model components are set, and the load calculation is performed to move all aggregates within the target size range; finally, the final coordinates of the transition aggregate vertices are extracted, and based on the extension distance of the set aggregate vertices, the transition aggregate is scaled and processed to restore the aggregate, and the final target polyhedral aggregate model is obtained. The volume rate of the model aggregate of the invention is high and determined.

Description

Technical Field

[0001] The invention belongs to the technical field of establishing a concrete finite element numerical model, and in particular relates to a method, an electronic device and a storage medium for establishing a target polyhedral aggregate model with a high volume ratio. Background Art

[0002] The practicality and reliability of the results of concrete numerical simulation first depend on the established numerical model. In the past research on concrete three-dimensional numerical models, researchers have been able to generate aggregates in the form of spheres and ellipsoids, and the study found that the strength results obtained by model simulation using spheres and ellipsoid aggregates are relatively safe. At present, it is possible to use polyhedral aggregates that reflect the diversity and angularity of aggregates to construct concrete numerical models. Some studies use the "take-and-place" method to place polyhedral aggregates into the model area, but since each aggregate placed needs to be judged for intrusion with the aggregates that have been successfully placed, the volume rate of polyhedral aggregates in the concrete model constructed by this method is usually low due to the limitation of computer computing power. In order to obtain a model with a higher volume rate, some scholars have calculated the distance between aggregates and used measures such as translating and rotating aggregates to improve the density of aggregates. However, due to the difficulty of calculating the distance between polyhedral aggregates, this method is mostly applied to regular aggregates. Some scholars have carried out random aggregate modeling based on software such as MATLAB, PYTHON, and FORTRAN, and then imported the established models into finite element software such as ABAQUS, ANSYS, and LS-DYNA for pre-processing and post-processing, and then simulated calculations to encrypt the aggregates. However, the aggregate volume ratio of the final model established by this method cannot be determined, and the gaps between aggregates required for subsequent simulation calculations cannot be guaranteed. In addition, the model needs to be converted between different software, and the operation is relatively cumbersome. Summary of the invention

[0003] The problem to be solved by the present invention is to optimize the volume ratio of a target polyhedral aggregate model, and a method, electronic equipment and storage medium for establishing a target polyhedral aggregate model with a high volume ratio are proposed.

[0004] To achieve the above object, the present invention is implemented through the following technical solutions:

[0005] A method for establishing a high volume ratio target polyhedral aggregate model comprises the following steps:

[0006] S1. Determine the parameters required to establish the target polyhedron aggregate model, including length, width, height, volume ratio of aggregate, proportion of aggregate gradation segment to aggregate in the gradation segment, and extension distance of aggregate vertices, and then determine the parameters of the basic polyhedron aggregate model based on the target polyhedron aggregate model parameters to generate the basic polyhedron aggregate model space;

[0007] S2. According to the target polyhedral aggregate model parameters and the basic polyhedral aggregate model parameters determined in step S1, the coordinate data of the vertices and placement positions of the initial polyhedral aggregate that meet the volume fraction requirements are generated in the basic polyhedral aggregate model space;

[0008] S3. Based on the extension distance of the aggregate vertices set in step S1, the coordinate data of the aggregate vertices after extension is calculated based on the initial polyhedral aggregate coordinate data obtained in step S2, and then the transition aggregate is generated and the delivery is completed to generate a basic polyhedral aggregate model, and the components of the basic polyhedral aggregate model include transition aggregate, outer frame and cover plate;

[0009] S4. Pre-processing the basic polyhedral aggregate model obtained in step S3, setting the material parameters, mesh parameters, load parameters of the basic polyhedral aggregate model components, and performing load calculation to obtain the calculation results;

[0010] S5. Based on the calculation results obtained in step S4, the initial coordinate data and displacement data of the transition aggregate vertex in the last frame are extracted to calculate the final coordinates of the transition aggregate vertex;

[0011] S6. Based on the extension distance of the aggregate vertices set in step S1 and the final coordinates of the transition aggregate vertices obtained in step S5, scaling processing is performed to restore the aggregate to obtain the final target polyhedral aggregate model.

[0012] Furthermore, the specific implementation method of step S1 includes the following steps:

[0013] S1.1. Set the parameters of the target polyhedral aggregate model, including the length a, width b, height c, aggregate volume ratio P, aggregate grading section [d i ,d i+1 ] and the proportion of aggregate in this grading section v i , the extension distance t of the aggregate vertex;

[0014] S1.2. Determine the parameters of the basic polyhedron aggregate model based on the target polyhedron aggregate model parameters, including the length a, width b, and height h of the basic polyhedron aggregate model, and calculate the volume ratio P' of the basic polyhedron aggregate, which is expressed as:

[0015]

[0016] S1.3. Generate a basic polyhedral aggregate model space based on the parameters of the basic polyhedral aggregate model in step S1.2.

[0017] Furthermore, the specific implementation method of step S2 includes the following steps:

[0018] S2.1. Based on the target polyhedral aggregate model parameters determined in step S1, define the initial aggregate particle size function, and in the aggregate grading section [d i ,d i+1 The initial aggregate size d is randomly generated within the range of i0 ;

[0019] S2.2. Determine the spatial range of aggregate to be placed, the expression is:

[0020] x i ∈(t+r i0 ,atr i0 )

[0021] y i ∈(t+r i0 ,btr i0 )

[0022] z i ∈(t+r i0 ,htr i0 )

[0023] Among them, (x i ,y i ,z i ) is the randomly generated coordinate of the center point of the i-th aggregate to be placed (x i ,y i ,z i ), r i0 is the initial radius of the i-th aggregate to be placed;

[0024] S2.3. Make an intersection judgment on the i-th aggregate to be placed and the j-th aggregate that has been generated and placed. The expression is:

[0025]

[0026] Among them, (x j ,y j ,z j ) is the coordinate of the center point of the jth generated aggregate, r j0 is the initial radius of the jth generated aggregate;

[0027] If the inequality is satisfied, the coordinates of the center point are stored; if the inequality is not satisfied, the process returns to S2.2 to regenerate the coordinates of the center point of the i-th aggregate to be placed;

[0028] S2.4. randomly select points on a sphere with the initial particle size of the i-th aggregate to be placed as the diameter, and obtain the vertex coordinates of the aggregate to be placed;

[0029] S2.5. Determine whether the volume of the aggregate generated for the grading is smaller than the volume of the aggregate to be placed for the grading. If so, return to S2.1 to continue generating the aggregate to be placed until the volume fraction requirement of the aggregate for the grading is met. Traverse all gradings until the volume fraction requirements of all graded aggregates are met to obtain the initial polyhedral aggregate.

[0030] Furthermore, the specific implementation method of step S3 includes the following steps:

[0031] S3.1. Calculate the centroid coordinates (x iQ ,y iQ ,z iQ ), the expression is:

[0032]

[0033] Where m is the total number of vertices of the i-th initial polyhedron aggregate, (x ik0 ,y ik0 ,z ik0 ) is the coordinate of the kth vertex of the ith initial polyhedron aggregate;

[0034] S3.2. Calculate the unit vectors of the lines connecting the vertices and the center of gravity of the i-th initial polyhedral aggregate:

[0035]

[0036] in, is the unit vector in the direction of the line connecting the kth vertex of the i-th initial polyhedron aggregate and the center of gravity, k = 1, 2, …, m.

[0037] S3.3. Calculate the coordinates of each vertex of the i-th initial polyhedral aggregate after extension, and obtain the coordinates of each vertex of the i-th transition aggregate:

[0038] Among them, (x ik ,y ik ,z ik ) is the coordinate of the kth vertex of the ith transition aggregate, k = 1, 2, …, m;

[0039] S3.4. Generate transition aggregate based on the vertex coordinates of the transition aggregate obtained in step S3.3, and then place the transition aggregate based on the center point coordinates obtained in step S2.3;

[0040] S3.5. Based on the parameters of the basic polyhedron aggregate model obtained in step S1, generate an outer frame of the basic polyhedron aggregate model with a thickness of d and a cover plate of the basic polyhedron aggregate model with a size of a×b×d to obtain a basic polyhedron aggregate model.

[0041] Furthermore, the specific implementation method of step S4 includes the following steps:

[0042] S4.1. Set the material properties of the components of the basic polyhedral aggregate model, and keep the density and elastic modulus of the transition aggregate, outer frame and cover plate consistent;

[0043] S4.2. Set the analysis step: Use the dynamic explicit analysis step and determine the time length t;

[0044] S4.3. Set up the interaction between the components of the basic polyhedral aggregate model: the interaction between the transition aggregate, the outer frame and the cover plate adopts explicit general contact, and the tangential behavior considers friction;

[0045] S4.4. Apply load: Apply load to all transition aggregates in layers and to the cover plate to ensure that all transition aggregates eventually move to the range of a×b×c under the load. The load includes body force and gravity.

[0046] S4.5. Discretize the components of the basic polyhedral aggregate model into finite element meshes and submit the calculation to obtain the result file.

[0047] Furthermore, the specific implementation method of step S5 includes the following steps:

[0048] S5.1. Extract the initial coordinate data and displacement data of the transition aggregate vertex in the last frame from the result file obtained in step S4, and then calculate the final coordinates of the transition aggregate vertex The expression is:

[0049]

[0050] in, is the coordinate of the final position of the kth vertex of the i-th transition aggregate after loading, is the distance moved by the kth vertex of the ith transition aggregate under load, k = 1, 2, …, m.

[0051] Furthermore, the specific implementation method of step S6 includes the following steps:

[0052] S6.1. Based on the extension distance of the aggregate vertex set in step S1 and the final coordinates of the transition aggregate vertex obtained in step S5, the final coordinates of the unextended aggregate are calculated by scaling. The expression is:

[0053]

[0054] in, is the final coordinate of the kth vertex of the ith unextended aggregate, k = 1, 2, …, m;

[0055] S6.2. Based on the final coordinates of the unextended aggregate obtained in step S6.1, generate the final target polyhedral aggregate model.

[0056] An electronic device comprises a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of a method for establishing a target polyhedral aggregate model with a high volume ratio when executing the computer program.

[0057] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements a method for establishing a target polyhedral aggregate model with a high volume ratio.

[0058] Beneficial effects of the present invention:

[0059] The method for establishing a target polyhedral aggregate model with a high volume ratio described in the present invention can obtain a polyhedral aggregate model with a relatively high volume ratio by encrypting the aggregate, such as a polyhedral aggregate model with an aggregate volume ratio of 50%.

[0060] The method for establishing a target polyhedral aggregate model with a high volume ratio described in the present invention can obtain a polyhedral aggregate model with a certain expected aggregate volume ratio by moving all transition aggregates to the area of ​​the model to be built through load; by first extending to transition aggregates for encryption and then restoring to unextended aggregates, this method can ensure the gaps between aggregates required for subsequent simulation calculations;

[0061] The method for establishing a target polyhedral aggregate model with a high volume ratio described in the present invention does not require linkage between multiple softwares and is easy to operate. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] Figure 1 A flow chart of a method for establishing a high volume ratio target polyhedral aggregate model according to the present invention;

[0063] Figure 2 Schematic diagram of the basic polyhedral aggregate model, with an aggregate volume fraction of 17%;

[0064] Figure 3 It is a schematic diagram of the internal aggregate distribution of the basic polyhedral aggregate model;

[0065] Figure 4 It is a schematic diagram that all aggregates are moved into the range of the target polyhedral aggregate model;

[0066] Figure 5 Schematic diagram of the random aggregate model with an aggregate volume fraction of 50%. DETAILED DESCRIPTION

[0067] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention, that is, the specific embodiments described are only part of the embodiments of the present invention, rather than all of the specific embodiments. The components of the specific embodiments of the present invention described and shown in the drawings herein can be arranged and designed in various different configurations, and the present invention can also have other embodiments.

[0068] Therefore, the following detailed description of the specific embodiments of the present invention provided in the accompanying drawings is not intended to limit the scope of the invention claimed for protection, but merely represents the selected specific embodiments of the present invention. Based on the specific embodiments of the present invention, all other specific embodiments obtained by those skilled in the art without making creative work are within the scope of protection of the present invention.

[0069] In order to further understand the content, features and effects of the present invention, the following specific implementation methods are given as examples, and the attached Figure 1 -Attached Figure 5 The detailed instructions are as follows:

[0070] Embodiment 1:

[0071] A method for establishing a high volume ratio target polyhedral aggregate model comprises the following steps:

[0072] S1. Determine the parameters required to establish the target polyhedron aggregate model, including length, width, height, volume ratio of aggregate, proportion of aggregate gradation segment to aggregate in the gradation segment, and extension distance of aggregate vertices, and then determine the parameters of the basic polyhedron aggregate model based on the target polyhedron aggregate model parameters to generate the basic polyhedron aggregate model space;

[0073] Furthermore, the specific implementation method of step S1 includes the following steps:

[0074] S1.1. Set the parameters of the target polyhedral aggregate model, including the length a, width b, height c, aggregate volume ratio P, aggregate grading section [d i ,d i+1 ] and the proportion of aggregate in this grading section v i , the extension distance t of the aggregate vertex;

[0075] S1.2. Determine the parameters of the basic polyhedron aggregate model based on the target polyhedron aggregate model parameters, including the length a, width b, and height h of the basic polyhedron aggregate model, and calculate the volume ratio P' of the basic polyhedron aggregate, which is expressed as:

[0076]

[0077] S1.3. Generate a basic polyhedral aggregate model space based on the parameters of the basic polyhedral aggregate model in step S1.2;

[0078] S2. According to the target polyhedral aggregate model parameters and the basic polyhedral aggregate model parameters determined in step S1, the coordinate data of the vertices and placement positions of the initial polyhedral aggregate that meet the volume fraction requirements are generated in the basic polyhedral aggregate model space;

[0079] Furthermore, the specific implementation method of step S2 includes the following steps:

[0080] S2.1. Based on the target polyhedral aggregate model parameters determined in step S1, define the initial aggregate particle size function, and in the aggregate grading section [d i ,d i+1 The initial aggregate size d is randomly generated within the range of i0 ;

[0081] S2.2. Determine the spatial range of aggregate to be placed, the expression is:

[0082] x i ∈(t+r i0 ,atr i0 )

[0083] y i ∈(t+r i0 ,btr i0 )

[0084] z i ∈(t+r i0 ,htr i0 )

[0085] Among them, (x i ,y i ,z i ) is the randomly generated coordinate of the center point of the i-th aggregate to be placed (x i ,y i ,z i ), r i0 is the initial radius of the i-th aggregate to be placed;

[0086] S2.3. Make an intersection judgment on the i-th aggregate to be placed and the j-th aggregate that has been generated and placed. The expression is:

[0087]

[0088] Among them, (x j ,y j ,zj ) is the coordinate of the center point of the jth generated aggregate, r j0 is the initial radius of the jth generated aggregate;

[0089] If the inequality is satisfied, the coordinates of the center point are stored; if the inequality is not satisfied, the process returns to S2.2 to regenerate the coordinates of the center point of the i-th aggregate to be placed;

[0090] S2.4. randomly select points on a sphere with the initial particle size of the i-th aggregate to be placed as the diameter, and obtain the vertex coordinates of the aggregate to be placed;

[0091] S2.5. Determine whether the volume of the aggregate generated for the grading is less than the volume of the aggregate to be placed for the grading. If so, return to S2.1 to continue generating the aggregate to be placed until the volume fraction requirement of the grading aggregate is met. Traverse all gradings until the volume fraction requirements of all grading aggregates are met to obtain the initial polyhedral aggregate.

[0092] S3. Based on the extension distance of the aggregate vertices set in step S1, the coordinate data of the aggregate vertices after extension is calculated based on the initial polyhedral aggregate coordinate data obtained in step S2, and then the transition aggregate is generated and the delivery is completed to generate a basic polyhedral aggregate model, and the components of the basic polyhedral aggregate model include transition aggregate, outer frame and cover plate;

[0093] Furthermore, the specific implementation method of step S3 includes the following steps:

[0094] S3.1. Calculate the centroid coordinates (x iQ ,y iQ ,z iQ ), the expression is:

[0095]

[0096] Where m is the total number of vertices of the i-th initial polyhedron aggregate, (x ik0 ,y ik0 ,z ik0 ) is the coordinate of the kth vertex of the ith initial polyhedron aggregate;

[0097] S3.2. Calculate the unit vectors of the lines connecting the vertices and the center of gravity of the i-th initial polyhedral aggregate:

[0098]

[0099] in, is the unit vector in the direction of the line connecting the kth vertex of the i-th initial polyhedron aggregate and the center of gravity, k = 1, 2, …, m.

[0100] S3.3. Calculate the coordinates of each vertex of the i-th initial polyhedral aggregate after extension, and obtain the coordinates of each vertex of the i-th transition aggregate:

[0101] Among them, (x ik ,y ik ,z ik ) is the coordinate of the kth vertex of the ith transition aggregate, k = 1, 2, …, m;

[0102] S3.4. Generate transition aggregate based on the vertex coordinates of the transition aggregate obtained in step S3.3, and then place the transition aggregate based on the center point coordinates obtained in step S2.3;

[0103] S3.5. Based on the parameters of the basic polyhedron aggregate model obtained in step S1, an outer frame of the basic polyhedron aggregate model with a thickness of d and a cover plate of the basic polyhedron aggregate model with a size of a×b×d are generated to obtain a basic polyhedron aggregate model;

[0104] S4. Pre-processing the basic polyhedral aggregate model obtained in step S3, setting the material parameters, mesh parameters, load parameters of the basic polyhedral aggregate model components, and performing load calculation to obtain the calculation results;

[0105] Furthermore, the specific implementation method of step S4 includes the following steps:

[0106] S4.1. Set the material properties of the components of the basic polyhedral aggregate model, and keep the density and elastic modulus of the transition aggregate, outer frame and cover plate consistent;

[0107] S4.2. Set the analysis step: Use the dynamic explicit analysis step and determine the time length t;

[0108] S4.3. Set up the interaction between the components of the basic polyhedral aggregate model: the interaction between the transition aggregate, the outer frame and the cover plate adopts explicit general contact, and the tangential behavior considers friction;

[0109] S4.4. Apply load: Apply load to all transition aggregates in layers and to the cover plate to ensure that all transition aggregates eventually move to the range of a×b×c under the load. The load includes body force and gravity.

[0110] S4.5. Discretize the components of the basic polyhedral aggregate model into finite element meshes and submit the calculation to obtain the result file.

[0111] S5. Based on the calculation results obtained in step S4, the initial coordinate data and displacement data of the transition aggregate vertex in the last frame are extracted to calculate the final coordinates of the transition aggregate vertex;

[0112] Furthermore, the specific implementation method of step S5 includes the following steps:

[0113] S5.1. Extract the initial coordinate data and displacement data of the transition aggregate vertex in the last frame from the result file obtained in step S4, and then calculate the final coordinates of the transition aggregate vertex The expression is:

[0114]

[0115] in, is the coordinate of the final position of the kth vertex of the i-th transition aggregate after loading, is the distance moved by the kth vertex of the i-th transition aggregate under load, k = 1, 2, …, m;

[0116] S6. Based on the extension distance of the aggregate vertices set in step S1 and the final coordinates of the transition aggregate vertices obtained in step S5, scaling processing is performed to restore the aggregate to obtain the final target polyhedral aggregate model.

[0117] Furthermore, the specific implementation method of step S6 includes the following steps:

[0118] S6.1. Based on the extension distance of the aggregate vertex set in step S1 and the final coordinates of the transition aggregate vertex obtained in step S5, the final coordinates of the unextended aggregate are calculated by scaling. The expression is:

[0119]

[0120] in, is the final coordinate of the kth vertex of the ith unextended aggregate, k = 1, 2, …, m;

[0121] S6.2. Based on the final coordinates of the unextended aggregate obtained in step S6.1, generate the final target polyhedral aggregate model.

[0122] Embodiment 2:

[0123] The difference between this embodiment and embodiment 1 is that in the method for establishing a target polyhedral aggregate model with a high volume rate, step S2.4 randomly selects points on a spherical surface with an initial radius of the i-th aggregate to be placed as the diameter, obtains the vertex coordinates of the aggregate to be placed, and then randomly fluctuates the vertex coordinates of the aggregate to be placed along the coordinates of the aggregate and the center point to obtain a concave-convex polyhedral aggregate, and calculates the volume V of the concave-convex polyhedral aggregate. i .

[0124] Embodiment 3:

[0125] An electronic device includes a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of a method for establishing a target polyhedral aggregate model with a high volume ratio described in Example 1 or Example 2 are implemented.

[0126] The computer device of the present invention may be a device including a processor and a memory, such as a single chip microcomputer including a central processing unit, etc. Furthermore, the processor is used to implement the steps of the above-mentioned method for establishing a target polyhedral aggregate model with a high volume ratio when executing the computer program stored in the memory.

[0127] The processor may be a central processing unit (CPU), 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.

[0128] The memory may mainly include a program storage area and a data storage area, wherein the program storage area may 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 may store data created according to the use of the mobile phone (such as audio data, a phone book, etc.), etc. In addition, the memory may include a high-speed random access memory, and may 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.

[0129] Embodiment 4:

[0130] A computer-readable storage medium having a computer program stored thereon, characterized in that when the computer program is executed by a processor, the method for establishing a target polyhedral aggregate model with a high volume ratio as described in any one of claims 1 to 7 is implemented.

[0131] The computer-readable storage medium of the present invention can be any form of storage medium that can be read by a processor of a computer device, including but not limited to non-volatile memory, volatile memory, ferroelectric memory, etc. The computer-readable storage medium stores a computer program. When the processor of the computer device reads and executes the computer program stored in the memory, the steps of the above-mentioned method for establishing a target polyhedral aggregate model with a high volume ratio can be implemented.

[0132] The computer program includes computer program code, which may 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, USB flash drive, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electric carrier signal, telecommunication signal and software distribution medium, etc.

[0133] It should be noted that relational terms such as "first" and "second" are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the existence of other identical elements in the process, method, article or device including the elements.

[0134] Although the present application has been described above with reference to specific embodiments, various modifications may be made thereto and parts thereof may be replaced with equivalents without departing from the scope of the present application. In particular, as long as there is no structural conflict, the various features in the specific embodiments disclosed in the present application may be used in combination with each other in any manner, and the fact that these combinations are not exhaustively described in this specification is only for the sake of omitting space and saving resources. Therefore, the present application is not limited to the specific embodiments disclosed herein, but includes all technical solutions falling within the scope of the claims.

Claims

1. A method for establishing a high volume ratio polyhedral aggregate model, characterized in that: The steps include: S1. Determine the parameters required to establish the target polyhedron aggregate model, including length, width, height, volume ratio of aggregate, proportion of aggregate gradation segment to aggregate in the gradation segment, and extension distance of aggregate vertices, and then determine the parameters of the basic polyhedron aggregate model based on the target polyhedron aggregate model parameters to generate the basic polyhedron aggregate model space; S2. According to the target polyhedral aggregate model parameters and the basic polyhedral aggregate model parameters determined in step S1, the coordinate data of the vertices and placement positions of the initial polyhedral aggregate that meet the volume fraction requirements are generated in the basic polyhedral aggregate model space; S3. Based on the extension distance of the aggregate vertices set in step S1, the coordinate data of the aggregate vertices after extension is calculated based on the initial polyhedral aggregate coordinate data obtained in step S2, and then the transition aggregate is generated and the delivery is completed to generate a basic polyhedral aggregate model, and the components of the basic polyhedral aggregate model include transition aggregate, outer frame and cover plate; S4. Pre-processing the basic polyhedral aggregate model obtained in step S3, setting the material parameters, mesh parameters, load parameters of the basic polyhedral aggregate model components, and performing load calculation to obtain the calculation results; S5. Based on the calculation results obtained in step S4, the initial coordinate data and displacement data of the transition aggregate vertex in the last frame are extracted to calculate the final coordinates of the transition aggregate vertex; The specific implementation method of step S5 includes the following steps: S5.

1. Extract the initial coordinate data and displacement data of the transition aggregate vertex in the last frame from the result file obtained in step S4, and then calculate the final coordinates of the transition aggregate vertex The expression is: in, is the coordinate of the final position of the kth vertex of the i-th transition aggregate after loading, is the distance moved by the kth vertex of the i-th transition aggregate under load, k = 1, 2, ..., m, (x ik ,y ik ,z ik ) is the coordinate of the kth vertex of the ith transition aggregate; S6. Based on the extension distance of the aggregate vertices set in step S1 and the final coordinates of the transition aggregate vertices obtained in step S5, scaling is performed to restore the aggregate to obtain the final target polyhedral aggregate model; The specific implementation method of step S6 includes the following steps: S6.

1. Based on the extension distance of the aggregate vertex set in step S1 and the final coordinates of the transition aggregate vertex obtained in step S5, the final coordinates of the unextended aggregate are calculated by scaling. The expression is: in, is the final coordinate of the kth vertex of the ith unextended aggregate, k = 1, 2, …, m; S6.

2. Based on the final coordinates of the unextended aggregate obtained in step S6.1, generate the final target polyhedral aggregate model.

2. The method for establishing a high volume ratio polyhedral aggregate model according to claim 1, characterized in that: The specific implementation method of step S1 includes the following steps: S1.

1. Set the parameters of the target polyhedral aggregate model, including the length a, width b, height c, aggregate volume ratio P, aggregate grading section [d i ,d i+1 ] and the proportion of aggregate in this grading section v i , the extension distance t of the aggregate vertex; S1.

2. Determine the parameters of the basic polyhedron aggregate model based on the target polyhedron aggregate model parameters, including the length a, width b, and height h of the basic polyhedron aggregate model, and calculate the volume ratio P' of the basic polyhedron aggregate, which is expressed as: S1.

3. Generate a basic polyhedral aggregate model space based on the parameters of the basic polyhedral aggregate model in step S1.

2.

3. A method for establishing a high volume ratio polyhedral aggregate model according to claim 2, characterized in that: The specific implementation method of step S2 includes the following steps: S2.

1. Based on the target polyhedral aggregate model parameters determined in step S1, define the initial aggregate particle size function, and in the aggregate grading segment [d i ,d i+1 The initial aggregate size d is randomly generated within the range of i0 ; S2.

2. Determine the spatial range of aggregate to be placed, the expression is: x i ∈(t+r i0 ,a-t-r i0 ) y i ∈(t+r i0 ,b-t-r i0 ) z i ∈(t+r i0 ,h-t-r i0 ) Among them, (x i ,y i ,z i ) is the randomly generated coordinate of the center point of the i-th aggregate to be placed (x i ,y i ,z i ), r i0 is the initial radius of the i-th aggregate to be placed; S2.

3. Make an intersection judgment on the i-th aggregate to be placed and the j-th aggregate that has been generated and placed. The expression is: Among them, (x j ,y j ,z j ) is the coordinate of the center point of the jth generated aggregate, r j0 is the initial radius of the jth generated aggregate; If the expression is satisfied, the coordinates of the center point are stored; if the expression is not satisfied, the process returns to S2.2 to regenerate the coordinates of the center point of the i-th aggregate to be placed; S2.

4. randomly select points on a sphere with the initial particle size of the i-th aggregate to be placed as the diameter, and obtain the vertex coordinates of the aggregate to be placed; S2.

5. Determine whether the volume of the aggregate generated for the grading section is less than the volume of the aggregate to be placed in the grading section. If so, return to S2.1 to continue generating the aggregate to be placed until the volume fraction requirement of the aggregate in the grading section is met. Traverse all grading sections until the volume fraction requirements of the aggregate in all grading sections are met to obtain the initial polyhedral aggregate.

4. The method for establishing a high volume ratio polyhedral aggregate model according to claim 3, characterized in that: The specific implementation method of step S3 includes the following steps: S3.

1. Calculate the centroid coordinates (x iQ ,y iQ ,z iQ ), the expression is: Where m is the total number of vertices of the i-th initial polyhedron aggregate, (x ik0 ,y ik0 ,z ik0 ) is the coordinate of the kth vertex of the i-th initial polyhedron aggregate; S3.

2. Calculate the unit vectors of the lines connecting the vertices and the center of gravity of the i-th initial polyhedral aggregate: in, is the unit vector of the line connecting the kth vertex of the i-th initial polyhedron aggregate and the center of gravity, k = 1, 2, ..., m, S3.

3. Calculate the coordinates of each vertex of the i-th initial polyhedral aggregate after extension, and obtain the coordinates of each vertex of the i-th transition aggregate: Among them, (x ik ,y ik ,z ik ) is the coordinate of the kth vertex of the i-th transition aggregate, k = 1, 2, …, m; S3.

4. Generate transition aggregate based on the vertex coordinates of the transition aggregate obtained in step S3.3, and then place the transition aggregate based on the center point coordinates obtained in step S2.3; S3.

5. Based on the parameters of the basic polyhedron aggregate model obtained in step S1, generate an outer frame of the basic polyhedron aggregate model with a thickness of d and a cover plate of the basic polyhedron aggregate model with a size of a×b×d to obtain a basic polyhedron aggregate model.

5. The method for establishing a high volume ratio polyhedral aggregate model according to claim 4, characterized in that: The specific implementation method of step S4 includes the following steps: S4.

1. Set the material properties of the components of the basic polyhedral aggregate model, and keep the density and elastic modulus of the transition aggregate, outer frame and cover plate consistent; S4.

2. Set the analysis step: Use the dynamic explicit analysis step and determine the time length t; S4.

3. Set up the interaction between the components of the basic polyhedral aggregate model: the interaction between the transition aggregate, the outer frame and the cover plate adopts explicit general contact, and the tangential behavior considers friction; S4.

4. Apply load: Apply load to all transition aggregates in layers and to the cover plate to ensure that all transition aggregates eventually move to the range of a×b×c under the load. The load includes body force and gravity. S4.

5. Discretize the components of the basic polyhedral aggregate model into finite element meshes and submit the calculation to obtain the result file.

6. An electronic device, characterized in that: The invention comprises a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of a method for establishing a high volume ratio polyhedral aggregate model as described in any one of claims 1 to 5 when executing the computer program.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for establishing a high volume ratio polyhedral aggregate model as described in any one of claims 1 to 5 is implemented.

Images