Method for predicting the effect of convex polyhedral reinforcement shape on the elastic deformation behavior of composites

Abstract

The application discloses a kind of convex polyhedron reinforcing phase shape prediction method to the elastic deformation behavior of composite material, belong to novel composite material model design field.The method is as follows:(1) establish original convex polyhedron model;(2) select cutting vertex, set cutting parameter;(3) calculate cutting point and cutting plane;(4) cut convex polyhedron;(5) construct reinforcing phase model;(6) establish composite material model and carry out mesh division;(7) boundary condition is applied to composite material model;(8) based on finite element simulation technology, the mechanical behavior of composite material is simulated numerically, obtains the elastic coefficient of composite material at room temperature.The method has the characteristics such as simple operation, reinforcing phase shape is rich, simulation precision is high.

Description

Technical Field

[0001] This invention belongs to the field of novel composite material model design, specifically relating to a method for predicting the elastic deformation behavior of composite materials based on the shape of a convex polyhedral reinforcing phase. Background Technology

[0002] In composite materials, reinforcing phase particles can take various shapes. The shape of these particles affects the mechanical properties of the composite. Sharp reinforcing phases can easily lead to stress concentration, causing premature local damage and reducing the material's ductility and toughness. Rounded reinforcing phases, on the other hand, can significantly alleviate stress concentration and greatly reduce the stress gradient in their vicinity. However, it is difficult to systematically measure the impact of different reinforcing phase particle shapes on the stress-strain distribution and elastic modulus of the composite material through experimental trial and error.

[0003] In recent years, finite element method (FEM) technology has been increasingly applied to the simulation of the mechanical behavior of composite materials, achieving relatively accurate performance prediction results. A novel modeling method for controlling the shape of reinforcing phase particles, combined with FEM simulation technology, is expected to achieve high-throughput calculations of the influence of different reinforcing phase particle shapes on the elastic behavior of composite materials, revealing the relationship between particle shape and the elastic properties of composite materials, meeting the needs of scientific research, and the calculation results can guide the development of new materials in engineering applications. Summary of the Invention

[0004] The purpose of this invention is to address the difficulty in modeling the shape of reinforcing phase particles in the existing technology by providing a method for predicting the elastic deformation behavior of composite materials based on the shape of convex polyhedral reinforcing phases.

[0005] The specific technical solution adopted in this invention is as follows:

[0006] This invention provides a method for predicting the elastic deformation behavior of composite materials based on the shape of the convex polyhedral reinforcing phase. The specific method is as follows:

[0007] S1: Create the original convex polyhedron in 3D modeling software;

[0008] S2: Select the center point of the original convex polyhedron created in step S1;

[0009] S3: Select all vertices of the original convex polyhedron created in step S1 that need to be cut, determine the area weight of the cut, and calculate the cutting planes; create a set of planes for all cutting planes, denoted as Set-Slice-Plane;

[0010] S4: Based on the calculated cutting plane, cut all vertices of the original convex polyhedron that need to be cut to obtain the set of cut polyhedra Set-Poly; if the cut polyhedron in the set Set-Poly is a convex polyhedron and its center point is within the polyhedron, then the cut polyhedron is recorded as a new convex polyhedron POLYHEDRON, and the remaining cutting results are deleted; if the number of cyclic cutting T = 0, then continue to step S5; if T > 0, then use the new convex polyhedron POLYHEDRON and execute steps S2 to S4 T times;

[0011] S5: Based on the new convex polyhedron POLYHEDRON after cutting in step S4, the geometric model of the enhanced phase is obtained;

[0012] S6: Based on the reinforcing phase geometric model, establish the composite material model and perform mesh generation;

[0013] S7: Input elastic properties and apply boundary conditions to the composite material model;

[0014] S8: Based on the composite material model processed in step S7, numerical simulation of the mechanical behavior of the composite material is performed to obtain the elastic coefficient of the composite material at room temperature.

[0015] Preferably, the 3D modeling software mentioned in step S1 includes AutoCAD, C4D, SolidWorks, NX, Rhinoceros, MAYA, 3DMAX, and Blender.

[0016] Preferably, the area weight of the cut is 0, 1 or 2.

[0017] Preferably, in step S3, all vertices of the original convex polyhedron are selected and denoted as the set Set-Vertex.

[0018] Furthermore, the specific method for calculating the cutting plane in step S3 is as follows:

[0019] Select a vertex p0(x0, y0, z0) within the set Set-Vertex, where vertex p0 is the intersection of the three edges and the three surfaces;

[0020] Obtain the normal vectors and areas of the three surfaces containing vertex p0: the normal vector of the first surface is n1(x1, y1, z1), and the area of ​​the first surface is A1; the normal vector of the second surface is n2(x2, y2, z2), and the area of ​​the second surface is A2; the normal vector of the third surface is n3(x3, y3, z3), and the area of ​​the third surface is A3; where n1, n2, and n3 are all directions away from the center point;

[0021] The cutting vector N(xn, yn, zn) is calculated based on the area weight of the cut; combined with the cutting depth S, where S>0, the coordinates of the cutting point ps are calculated as follows:

[0022] ps = N*S + p0

[0023] Based on the cutting point ps and the cutting vector N, the equation of the cutting plane is: A*x + B*y + C*z + D = 0; where A, B, C, and D are respectively:

[0024] A = nx

[0025] B = ny

[0026] C = nz

[0027] D = -A*x0 - B*y0 - C*z0.

[0028] Furthermore, the calculation process of the cutting vector N(xn, yn, zn) is as follows:

[0029] The vectors are summed based on the area weights of the cuts to obtain the summation vector NT(xnt, ynt, znt):

[0030] 1) If the area weight of the cut is 0, the summation vector NT is:

[0031] NT = n1 + n2 + n3

[0032] 2) If the area weight of the cut is 1, calculate the summation vector NT by combining the first surface area A1, the second surface area A2, and the third surface area A3:

[0033]

[0034] 3) If the area weight of the cut is 2, calculate the summation vector NT by combining the first surface area A1, the second surface area A2, and the third surface area A3:

[0035] NT = n1*A1 + n2*A2 + n3*A3

[0036] Find the length LNT of the summing vector NT:

[0037]

[0038] Then the XYZ components of the cutting vector N(xn, yn, zn) are:

[0039] xn = xnt / LNT

[0040] yn = ynt / LNT

[0041] zn = znt / LNT.

[0042] As a preferred embodiment, the specific process for establishing the composite material model in step S6 is as follows:

[0043] S61: Repeat steps S1 to S5 several times to cut and obtain several convex polyhedra POLYHEDRON, denoted as Set-Reins; create a cube CUBE with side length L and volume L×L×L; set the volume fraction VL;

[0044] S62: Sum the volumes of all convex polyhedra POLYHEDRON in the set-Reins, denoted as VOLUME_ALL; Calculate the target volume V = VL*L*L*L based on the side length L of the cube CUBE and the set volume fraction VL;

[0045] S63: Calculate the scaling factor The volume of all convex polyhedra POLYHEDRON in the set Set-Reins is scaled using the proportional scaling factor SCALE.

[0046] S64: Create an empty Set-Model;

[0047] S65: Select the center point of all convex polyhedra POLYHEDRON in the set Set-Reins, and denote it as CENTER; move the convex polyhedra POLYHEDRON according to the vector -CENTER so that its center point is the origin O(0,0,0);

[0048] S66: Take three positive integers tx, ty, and tz in the range of 0 to 360 as the angle values ​​of rotation; rotate the convex polyhedron POLYHEDRON by tx along the X-axis, rotate the convex polyhedron POLYHEDRON by ty along the Y-axis, and rotate the convex polyhedron POLYHEDRON by tz along the Z-axis.

[0049] S67: Take three positive integers px, py, and pz as coordinate values ​​within the range of 0 to L; the convex polyhedron POLYHEDRON moves according to the vector (px, py, pz);

[0050] S68: If the convex polyhedron POLYHEDRON does not intersect with the set Set-Model, and all vertices of the convex polyhedron POLYHEDRON are inside the cube CUBE, then put the convex polyhedron POLYHEDRON into the set Set-Model; otherwise, repeat steps S65 to S68.

[0051] S69: Assemble the Set-Model and CUBE sets into a single entity to obtain the composite material model.

[0052] Furthermore, the composite material model is meshed with 4-node tetrahedral elements, and the total number of mesh elements is not less than 600,000.

[0053] Furthermore, the elastic properties of the reinforcing phase and the matrix are input into the Set-Model and the CUBE respectively; the elastic properties include Young's modulus and Poisson's ratio.

[0054] As a preferred embodiment, the boundary conditions applied to the composite material model in step S7 are as follows:

[0055] Apply principal strain ε along the X, Y, Z directions xx ε yy ε zz A shear strain ε is applied to the X=L surface of the model along the Y direction and the Y=L surface along the X direction. xy A shear strain ε is applied to the surface of model Y=L along the Z direction and the surface of model Z=L along the Y direction. yz A shear strain ε is applied to the Z=L surface of the model along the X direction and the X=L surface along the Z direction. zx ; principal strain ε xx ε yy ε zz and shear strain ε xy ε yz ε zx The value is 0 or 2%.

[0056] Compared with the prior art, the present invention has the following advantages:

[0057] (1) The method provided by the present invention can be applied to convex polyhedral particle models of arbitrary shapes, which is convenient for modeling, simple in operation steps and has few control parameters;

[0058] (2) This invention uses finite element simulation technology to predict the elastic mechanical properties of composite material models. The convex polyhedral particles used in the method provided by this invention produce a variety of shapes for the reinforcing phase, and the geometric model of the composite material can be obtained by adjusting the cutting process to closely resemble the actual shape of the reinforcing phase particles produced in industrial production. Attached Figure Description

[0059] Figure 1 This is a schematic diagram of the 20 convex polyhedra created by the 3D software in Example 1;

[0060] Figure 2 This is a schematic diagram showing the center point positions of 20 convex polyhedra selected in Example 1;

[0061] Figure 3 These are the 20 new convex polyhedra POLYHEDRONs after cutting in Example 1;

[0062] Figure 4The images show a comparison of the 20 convex polyhedra before and after cutting in Example 1.

[0063] Figure 5 This is a comparison chart of the predicted Young's modulus for Example 1. Detailed Implementation

[0064] The present invention will be further described and illustrated below with reference to the accompanying drawings and specific embodiments. The technical features of each embodiment of the present invention can be combined accordingly, provided that there is no mutual conflict.

[0065] This invention provides a method for predicting the elastic deformation behavior of composite materials based on the shape of a convex polyhedral reinforcing phase. The specific operation method is as follows:

[0066] 1. Set the number of particles (COUNT) and the volume of the polyhedron (VOL).

[0067] 2. Create primitive convex polyhedra of arbitrary shape in 3D modeling software.

[0068] In a preferred embodiment of the present invention, the step is as follows:

[0069] In 3D modeling software (AutoCAD, C4D, SolidWorks, NX, Rhinoceros, MAYA, 3DMAX, or Blender, etc.), create a convex polyhedron that simultaneously satisfies the following three conditions: (1) The volume of the convex polyhedron is VOL. (2) Each vertex of the convex polyhedron is formed by the intersection of three edges. (3) Each vertex of the convex polyhedron is the intersection of three surfaces.

[0070] 3. Select the center point of the original convex polyhedron created in step 2, and denote it as CENTER.

[0071] 4. Select all vertices of the original convex polyhedron created in step 2 that need to be cut, determine the area weights of the cuts, and calculate the cutting planes. Create a set of all cutting planes, denoted as Set-Slice-Plane.

[0072] In a preferred embodiment of the present invention, the step is as follows:

[0073] (1) Select all vertices of the original convex polyhedron and denote them as the set Set-Vertex.

[0074] (2) The area weight value of the cut is 0, 1 or 2.

[0075] (3) Calculate the cutting plane, the specific method is as follows:

[0076] Select a vertex p0(x0, y0, z0) within the set Set-Vertex. Vertex p0 is the intersection of the three edges and the three surfaces.

[0077] Obtain the normal vectors and areas of the three surfaces containing vertex p0: The first surface has a normal vector n1(x1, y1, z1) and an area of ​​A1. The second surface has a normal vector n2(x2, y2, z2) and an area of ​​A2. The third surface has a normal vector n3(x3, y3, z3) and an area of ​​A3. Here, n1, n2, and n3 represent directions away from the center point.

[0078] The cutting vector N(xn, yn, zn) is calculated based on the area weights of the cuts, as follows:

[0079] The vectors are summed based on the area weights of the cuts to obtain the summation vector NT(xnt, ynt, znt):

[0080] 1) If the area weight of the cut is 0, the summation vector NT is:

[0081] NT = n1 + n2 + n3

[0082] 2) If the area weight of the cut is 1, calculate the summation vector NT by combining the first surface area A1, the second surface area A2, and the third surface area A3:

[0083]

[0084] 3) If the area weight of the cut is 2, calculate the summation vector NT by combining the first surface area A1, the second surface area A2, and the third surface area A3:

[0085] NT = n1*A1 + n2*A2 + n3*A3

[0086] Find the length LNT of the summing vector NT:

[0087]

[0088] Then the XYZ components of the cutting vector N(xn, yn, zn) are:

[0089] xn = xnt / LNT

[0090] yn = ynt / LNT

[0091] zn = znt / LNT

[0092] (4) Combining the cutting depth S with S>0, the coordinates of the cutting point ps are calculated as follows:

[0093] ps = N*S + p0.

[0094] (5) Based on the cutting point ps and the cutting vector N, the equation of the cutting plane is: A*x + B*y + C*z + D = 0. Where A, B, C, and D are respectively:

[0095] A = nx

[0096] B = ny

[0097] C = nz

[0098] D = -A*x0 - B*y0 - C*z0.

[0099] 5. Based on the calculated cutting planes, cut all vertices of the original convex polyhedron that need to be cut. Specifically, in the 3D modeling software, the original convex polyhedron is cut according to the cutting plane set Set-Slice-Plane. The resulting set of cut polyhedra is denoted as Set-Poly.

[0100] If, in the set Set-Poly, the cut polyhedron is a convex polyhedron and its center point is within the polyhedron, then this cut polyhedron is recorded as the new convex polyhedron POLYHEDRON, and all other cut results are deleted. If the number of cuts in the loop T = 0, then proceed to step S5; if T > 0, then use the new convex polyhedron POLYHEDRON and execute steps two through five T times. Based on the new convex polyhedron POLYHEDRON after cutting, the reinforcing phase geometric model is obtained.

[0101] 6. Repeat steps 2 to 5 COUNT times to generate COUNT convex polyhedra POLYHEDRON, denoted as Set-Reins.

[0102] 7. Based on the reinforcing phase geometric model, establish a composite material model and perform mesh generation.

[0103] In a preferred embodiment of the present invention, the step is as follows:

[0104] (1) Construct a cube CUBE with side length L and volume L×L×L. Set the volume fraction VL.

[0105] (2) Sum the volumes of all convex polyhedra POLYHEDRON in the set-Reins, denoted as VOLUME_ALL. Calculate the target volume V = VL*L*L*L based on the side length L of the cube CUBE and the set volume fraction VL.

[0106] (3) Calculate the scaling factor. The volume of all convex polyhedra POLYHEDRON in the set Set-Reins is scaled using the proportional scaling factor SCALE.

[0107] (4) Create an empty set Set-Model.

[0108] (5) Select the center point of all convex polyhedra POLYHEDRON in the set Set-Reins, and denote it as CENTER. Move the convex polyhedra POLYHEDRON according to the vector -CENTER so that its center point is the origin O(0,0,0).

[0109] (6) Take three positive integers tx, ty, and tz within the range of 0 to 360 as the angle values ​​of rotation. Rotate the convex polyhedron POLYHEDRON by tx along the X-axis, rotate the convex polyhedron POLYHEDRON by ty along the Y-axis, and rotate the convex polyhedron POLYHEDRON by tz along the Z-axis.

[0110] (7) Take three positive integers px, py, and pz as coordinate values ​​within the range of 0 to L. The convex polyhedron POLYHEDRON moves according to the vector (px, py, pz).

[0111] (8) If the convex polyhedron POLYHEDRON does not intersect with the set Set-Model, and all vertices of the convex polyhedron POLYHEDRON are inside the cube CUBE, then put the convex polyhedron POLYHEDRON into the set Set-Model; otherwise, repeat steps S65 to S68.

[0112] (9) Assemble the Set-Model and CUBE into a solid to obtain the composite material model.

[0113] (10) Mesh the above composite material model. The mesh type is 4-node tetrahedral elements, and the total number of mesh elements is not less than 600,000.

[0114] 8. Input elastic properties and apply boundary conditions to the composite material model.

[0115] In a preferred embodiment of the present invention, the step is as follows:

[0116] Input the elastic properties of the reinforcing phase and the matrix into the Set-Model and CUBE sets respectively. The elastic properties include Young's modulus and Poisson's ratio.

[0117] The specific boundary conditions applied to the composite material model are as follows:

[0118] Apply principal strain ε along the X, Y, Z directions xx ε yy ε zz A shear strain ε is applied to the X=L surface of the model along the Y direction and the Y=L surface along the X direction. xyA shear strain ε is applied to the surface of model Y=L along the Z direction and the surface of model Z=L along the Y direction. yz A shear strain ε is applied to the Z=L surface of the model along the X direction and the X=L surface along the Z direction. zx Principal strain ε xx ε yy ε zz The value and shear strain ε xy ε yz ε zx The value is 0 or 2%.

[0119] 9. Based on the composite material model processed in the above steps, numerical simulations are performed on the mechanical behavior of the composite material to obtain the elastic modulus of the composite material at room temperature.

[0120] The practical application effects of the above-described method of the present invention will be illustrated below through specific embodiments.

[0121] Example 1

[0122] 1. Set the number of particles (COUNT) to 20 and the volume of the polyhedron (VOL) to 5000.

[0123] 2. Create 20 original convex polyhedra of arbitrary shapes in 3D modeling software, such as... Figure 1 As shown.

[0124] 3. Obtain the center points of all original convex polyhedra, such as Figure 2 As shown.

[0125] 4. Determine the area weight type value of the cut as 1, the cut depth as S as 5, and calculate the cut vector and cut point for all vertices of each original convex polyhedron. Based on the cut vector and cut point, obtain the set of cut planes, Set-Slice-Plane.

[0126] 5. In the 3D modeling software, the original convex polyhedron is cut according to the Set-Slice-Plane set of cutting surfaces. The convex polyhedron containing the center point CENTER is retained and defined as the new convex polyhedron POLYHEDRON; all other cutting results are deleted. The number of loop cuts is set to T = 0. The convex polyhedron has undergone one cut, yielding the enhanced phase geometry model REINFORCE.

[0127] VI. Repeat steps two through five 20 times to generate 20 convex polyhedra, denoted as Set-Reins, whose shapes are as follows: Figure 3 As shown.

[0128] Calculate the sum of volumes (VOLUME_ALL) of 20 convex polyhedra within the set Set-Reins. Set the side length of the cube CUBE to L = 100 and the volume fraction VL = 10%. Calculate the target volume V as:

[0129] V = VL * L * L * L

[0130] Calculate the scaling factor SCALE:

[0131]

[0132] The volume of the Set-Reins convex polyhedron is scaled according to the scaling factor SCALE.

[0133] Create an empty collection Set-Model.

[0134] Perform the following operation on the set-Reins convex polyhedron:

[0135] ① Obtain the center point CENTER of the convex polyhedron POLYHEDRON.

[0136] ② Move the convex polyhedron POLYHEDRON according to the vector -CENTER so that its center point is the origin O(0,0,0).

[0137] ③ Select three positive integers tx, ty, and tz within the range of 0 to 360 as the angle values ​​of the rotation angle.

[0138] ④ Rotate the convex polyhedron POLYHEDRON by tx along the X-axis, rotate the convex polyhedron POLYHEDRON by ty along the Y-axis, and rotate the convex polyhedron POLYHEDRON by tz along the Z-axis.

[0139] ⑤ Select three positive integers px, py, and pz as coordinate values ​​within the range of 0 to L.

[0140] ⑥ The convex polyhedron POLYHEDRON moves according to the vector (px, py, pz).

[0141] ⑦ If the convex polyhedron POLYHEDRON does not intersect with the set Set-Model, and all vertices of the convex polyhedron POLYHEDRON are inside the cube CUBE, then put the convex polyhedron POLYHEDRON into the set Set-Model; otherwise, repeat steps ③ to ⑥.

[0142] The convex polyhedral set-model represents the reinforcing phase particles of the composite material. The particle distribution before and after cutting is as follows: Figure 4 As shown.

[0143] The Set-Model and CUBE sets are assembled into a single solid to obtain the composite material model. This composite material model is then meshed using 4-node tetrahedral elements, with a minimum of 600,000 mesh elements. In this example, finite element models were created for both uncut and once-cut particles.

[0144] 8. Input material properties: The matrix is ​​a titanium alloy (Young's modulus 110 GPa, Poisson's ratio 0.33). The reinforcing phase is titanium boride (Young's modulus 560 GPa, Poisson's ratio 0.15). Apply principal strain ε simultaneously in the X, Y, and Z directions of the model. xx =2%, ε yy =ε zz =0%, shear strain ε xy ε yz ε zx The values ​​are all 0.

[0145] 9. Elastic deformation behavior of composite material models calculated using finite element software.

[0146] This embodiment establishes finite element models of titanium boride particle-reinforced composite materials, both uncut and cut once. Based on... Figure 4 It can be seen that after the particles are cut, the reinforcing phase particles change from sharp to rounded. Figure 5 It can be seen that the Young's modulus of the composite material decreases. This indicates that after cutting using the method provided by this invention, the stress concentration around the convex polyhedral particles is alleviated, and the total stress of the system decreases.

[0147] The embodiments described above are merely preferred embodiments of the present invention and are not intended to limit the invention. Those skilled in the art can make various changes and modifications without departing from the spirit and scope of the invention. Therefore, all technical solutions obtained through equivalent substitution or transformation fall within the protection scope of the present invention.

Claims

1. A method for predicting the elastic deformation behavior of composite materials based on the shape of a convex polyhedral reinforcing phase, characterized in that, The specific method is as follows: S1: Create the original convex polyhedron in 3D modeling software; S2: Select the center point of the original convex polyhedron created in step S1; S3: Select all vertices of the original convex polyhedron created in step S1 that need to be cut, determine the area weight of the cut, and calculate the cutting planes; create a set of planes for all cutting planes, denoted as Set-Slice-Plane; S4: Based on the calculated cutting plane, cut all the vertices of the original convex polyhedron that need to be cut to obtain the set of cut polyhedra Set-Poly; if the cut polyhedron in the set Set-Poly is a convex polyhedron and the center point is in the polyhedron, then the cut polyhedron is recorded as a new convex polyhedron POLYHEDRON, and the other cutting results are deleted. If the number of cyclic cuttings T = 0, then continue to step S5; if T > 0, then use a new convex polyhedron POLYHEDRON and execute steps S2~S4 T times. S5: Based on the new convex polyhedron POLYHEDRON after cutting in step S4, the geometric model of the enhanced phase is obtained; S6: Based on the reinforcing phase geometric model, establish the composite material model and perform mesh generation; S7: Input elastic properties and apply boundary conditions to the composite material model; S8: Based on the composite material model processed in step S7, numerical simulation of the mechanical behavior of the composite material is performed to obtain the elastic coefficient of the composite material at room temperature. In step S3, all vertices of the original convex polyhedron are selected and denoted as Set-Vertex; the specific method for calculating the cutting plane in step S3 is as follows: Select a vertex p0 (x0, y0, z0) within the set Set-Vertex, where vertex p0 is the intersection of the three edges and the three surfaces; Obtain the normal vectors and areas of the three surfaces containing vertex p0: the normal vector of the first surface is n1 (x1, y1, z1), and the area of ​​the first surface is A1; the normal vector of the second surface is n2 (x2, y2, z2), and the area of ​​the second surface is A2; the normal vector of the third surface is n3 (x3, y3, z3), and the area of ​​the third surface is A3; where n1, n2, and n3 are all directions away from the center point; The cutting vector N(xn, yn, zn) is calculated based on the area weight of the cut; combined with the cutting depth S, where S > 0, the coordinates of the cutting point ps are calculated as follows: ； Based on the cutting point ps and the cutting vector N, the equation of the cutting plane is obtained as follows: Where A, B, C, and D are respectively: ； ； ； ； The specific process of establishing the composite material model in step S6 is as follows: S61: Repeat steps S1~S5 several times to cut and obtain several convex polyhedra POLYHEDRON, denoted as Set-Reins; create a cube CUBE with side length L and volume L×L×L; set the volume fraction VL; S62: Sum the volumes of all convex polyhedra POLYHEDRON in the set-Reins, denoted as VOLUME_ALL; calculate the target volume based on the side length L of the cube CUBE and the set volume fraction VL. ; S63: Calculate the scaling factor ;Scale the volume of all convex polyhedra POLYHEDRON in the set Set-Reins according to the proportional scaling factor SCALE; S64: Create an empty Set-Model; S65: Select the center point of all convex polyhedra POLYHEDRON in the set Set-Reins, and denote it as CENTER; move the convex polyhedra POLYHEDRON according to the vector -CENTER so that its center point is the origin O (0, 0, 0). S66: Take three positive integers tx, ty, and tz in the range of 0 to 360 as the angle values ​​of rotation; rotate the convex polyhedron POLYHEDRON by tx along the X-axis, rotate the convex polyhedron POLYHEDRON by ty along the Y-axis, and rotate the convex polyhedron POLYHEDRON by tz along the Z-axis. S67: Select three positive integers px, py, and pz as coordinate values ​​within the range of 0 to L; The convex polyhedron POLYHEDRON moves according to the vector (px, py, pz); S68: If the convex polyhedron POLYHEDRON does not intersect with the empty set Set-Model, and all vertices of the convex polyhedron POLYHEDRON are inside the cube CUBE, then put the convex polyhedron POLYHEDRON into the empty set Set-Model; otherwise, repeat steps S65 to S68. S69: Assemble the empty set Set-Model and the cube CUBE into a single entity to obtain the composite material model.

2. The method for predicting the elastic deformation behavior of composite materials based on the shape of the convex polyhedral reinforcing phase according to claim 1, characterized in that, The 3D modeling software mentioned in step S1 includes AutoCAD, C4D, SolidWorks, NX, Rhinoceros, MAYA, 3DMAX, and Blender.

3. The method for predicting the elastic deformation behavior of composite materials based on the shape of the convex polyhedral reinforcing phase according to claim 1, characterized in that, The area weight of the cut is 0, 1, or 2.

4. The method for predicting the elastic deformation behavior of composite materials based on the shape of the convex polyhedral reinforcing phase according to claim 1, characterized in that, The calculation process of the cutting vector N(xn, yn, zn) is as follows: The vectors are summed based on the area weights of the cuts to obtain the summation vector NT(xnt, ynt, znt): 1) If the area weight of the cut is 0, the summation vector NT is: ； 2) If the area weight of the cut is 1, calculate the summation vector NT by combining the first surface area A1, the second surface area A2, and the third surface area A3: ； 3) If the area weight of the cut is 2, calculate the summation vector NT by combining the first surface area A1, the second surface area A2, and the third surface area A3: ； Find the length LNT of the summing vector NT: ； Then the XYZ components of the cutting vector N(xn, yn, zn) are: ； ； 。 5. The method for predicting the elastic deformation behavior of composite materials based on the shape of the convex polyhedral reinforcing phase according to claim 1, characterized in that, The composite material model is meshed with 4-node tetrahedral elements, and the total number of mesh elements is not less than 600,000.

6. The method for predicting the elastic deformation behavior of composite materials based on the shape of the convex polyhedral reinforcing phase according to claim 5, characterized in that, Input the elastic properties of the reinforcing phase and the matrix into the empty set Set-Model and the cube CUBE respectively; the elastic properties include Young's modulus and Poisson's ratio.

7. The method for predicting the elastic deformation behavior of composite materials based on the shape of the convex polyhedral reinforcing phase according to claim 1, characterized in that, The specific boundary conditions applied to the composite material model in step S7 are as follows: Apply principal strains along the X, Y, and Z directions , , Shear strain is applied to the X=L surface of the model along the Y direction and the Y=L surface along the X direction. Shear strain is applied to the Y=L surface of the model along the Z direction and the Z=L surface along the Y direction. Shear strain is applied to the Z=L surface of the model along the X direction and the X=L surface along the Z direction. Main strain , , and shear strain , , The value is 0 or 2%.

Images