A construction method for a particle-fiber hybrid reinforcement model

Abstract

The present invention discloses a method for constructing a particulate fiber hybrid reinforcement model, which solves the problem that there is currently a lack of parametric and rapid establishment of numerical models for the microscopic configuration of particulate / fiber hybrid reinforced composites. The key points of the technical solution are to determine the representative unit type and corresponding size parameters, determine the particulate parameters and fiber parameters, establish a matrix model, calculate the initial number of particles according to the particulate parameters, establish the corresponding Thiessen polygon structure, and perform proportional scaling to form an initial particulate model. A fiber model is established by random placement according to the representative unit size parameters and fiber parameters. The fiber model minus the initial particulate model gives the target particulate model. The fiber model, the target particulate model and the matrix model are merged to form a particulate / fiber hybrid reinforcement model that meets the set parameters. The method for constructing a particulate fiber hybrid reinforcement model of the present invention can fill the gap in the prior art and can rapidly and effectively construct a particulate / fiber hybrid reinforcement model.

Description

Technical Field

[0001] The present invention relates to a finite element modeling method for composite materials, and particularly to a method for constructing a particle-fiber hybrid reinforcement model. Background Art

[0002] Composite materials made by adding a second high-strength reinforcing phase to a matrix have advantages such as high strength, high elastic modulus, wear resistance, good electrical and thermal conductivity compared to single materials, and are widely used in industries such as aerospace, electronics, automotive, and construction.

[0003] In recent years, the research on composite materials has shown a development trend from single reinforcement to multi-reinforcement. Especially for particle / fiber hybrid reinforced composite materials, their toughness has been improved compared to particle or fiber reinforced composite materials. Numerical analysis is one of the important means for the design of advanced composite materials. Compared with experimental analysis, it consumes less materials and time. However, for particle / fiber hybrid reinforced composite materials, there is currently no relatively convenient numerical modeling method proposed. Summary of the Invention

[0004] The purpose of the present invention is to provide a method for constructing a particle-fiber hybrid reinforcement model, which can fill the gaps in the existing technology and can quickly and effectively construct a particle / fiber hybrid reinforcement model.

[0005] The above technical purpose of the present invention is achieved through the following technical solutions:

[0006] A method for constructing a particle-fiber hybrid reinforcement model includes the following steps:

[0007] S1. Determine the representative unit type and its corresponding size parameters, and determine the particle parameters and fiber parameters;

[0008] S2. In the finite element software Abaqus, establish a matrix model;

[0009] S3. Calculate the initial number of particles according to the particle parameters;

[0010] S4. According to the representative unit size parameters and the number of particles, establish a corresponding Voronoi polygon structure in the finite element software Abaqus, and remove the outermost cells of the polygon structure;

[0011] S5. Scale the Voronoi polygon structure proportionally to form an initial particle model;

[0012] S6. According to the representative unit size parameters and the fiber parameters, use the random placement method to establish a fiber model in the finite element software Abaqus;

[0013] S7. In the finite element software Abaqus, subtract the initial particle model from the fiber model, and remove the broken particles and redundant particles to obtain a target particle model that meets the set parameters.

[0014] S8. Combine the fiber model, the target particle model, and the matrix model to form a particle / fiber hybrid reinforced model that meets the set parameters.

[0015] In summary, the present invention has the following beneficial effects:

[0016] Aiming at the problem that there is a lack of a parametric and rapid method for establishing a numerical model for the microscopic configuration of particle / fiber hybrid reinforced composites, by utilizing the feature that each cell of the Voronoi polygon has a convex polygon, a single cell after the discretization of the Voronoi polygon is used to represent discrete particles with complex morphologies; a random placement method is adopted to establish a fiber model, and by writing a plug-in for the Abaqus finite element software, a method for rapidly and effectively constructing a particle / fiber hybrid reinforced model is established. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] Figure 1 is a flowchart showing the process of this method;

[0018] Figure 2 is a schematic diagram showing the process of obtaining the initial particle model of a two-dimensional model;

[0019] Figure 3 is a schematic diagram showing the process of obtaining the particle / fiber hybrid reinforced model of a two-dimensional model;

[0020] Figure 4 is a schematic diagram showing the process of obtaining the initial particle model of a three-dimensional model;

[0021] Figure 5 is a schematic diagram showing the process of obtaining the particle / fiber hybrid reinforced model of a three-dimensional model. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0022] The present invention will be further described in detail below with reference to the accompanying drawings.

[0023] According to one or more embodiments, a method for constructing a particle / fiber hybrid reinforced model is disclosed. As Figure 1 shown, it includes the following steps:

[0024] S1. Determine the representative unit type and its corresponding size parameters, and determine the particle parameters and fiber parameters.

[0025] The representative unit type is two-dimensional or three-dimensional, and the size parameters correspond to the representative unit type. The particle parameters include the equivalent radius r p of the particle, and the particle volume fraction V f1 ; the fiber parameters include the equivalent radius rf , the fiber volume fraction V f2 , the fiber length range (l min , l max ); the representative unit size parameters include the length L, width W, and height H. Correspondingly, for the two-dimensional representative unit, it is [(L + 12r p ), (W + 12r p )], or for the three-dimensional representative unit, it is [(L + 8r p ), (W + 8r p ), (H + 8r p )].

[0026] S2. In the finite element software Abaqus, establish a matrix model.

[0027] The matrix model is a two-dimensional model centered at (0, 0) with length and width (L, W) corresponding to the two-dimensional representative unit, or a three-dimensional model centered at (0, 0) with length, width, and height (L, W, H) corresponding to the three-dimensional representative unit.

[0028] S3. Calculate the initial number of particles according to the particle parameters.

[0029] According to the equivalent radius r p of the particle and the particle volume fraction V f1 , calculate the initial number of particles:

[0030] For the two-dimensional model, the number of particles n p-2 is rounded down;

[0031] For the three-dimensional model, the number of particles n p-3 is rounded down.

[0032] where f s is the multiplication coefficient, and the value range is preferably set to (1.8, 2).

[0033] S4. According to the representative unit size parameters [(L + 12r p ), (W + 12r p )] (two-dimensional) or [(L + 8r p ), (W + 8r p ), (H + 8r p )] (three-dimensional), the number of particles n p-2 (two-dimensional) or n p-3 (three-dimensional), use Python language for secondary development, establish the corresponding two-dimensional or three-dimensional Voronoi polygon structure in the finite element software Abaqus, and remove the outermost cells of the polygon structure.

[0034] S5. Scale the Thiessen polygon structure proportionally to form an initial particle model.

[0035] Scale down the Thiessen polygon structure to form a particle group.

[0036] For a two-dimensional model, the scaling factor

[0037] For a three-dimensional model, the scaling factor

[0038] Moreover, remove the particles outside the matrix, and remove the particles whose 75% area at the matrix boundary is outside the matrix; remove the part of the particle group outside the matrix range to obtain the initial particle model. Generate an initial particle model with the same distribution and the same volume, where the particles are circular / spherical, according to the centroid position and volume of the particles.

[0039] S6. According to the representative unit size parameters: (L, W) (two-dimensional) or (L, W, H) (three-dimensional); fiber parameters: fiber equivalent radius r f , fiber volume fraction V f2 , fiber length range (l min , l max ); Use Python language for secondary development, and adopt the random placement method to establish a fiber model in the finite element software Abaqus.

[0040] S7. In the finite element software Abaqus, subtract the initial particle model from the fiber model, and remove the broken particles and redundant particles to obtain a target particle model that meets the set parameters.

[0041] S8. Merge the fiber model, the target particle model and the matrix model to form a particle / fiber hybrid reinforced model that meets the set parameters.

[0042] For clarity, the following examples with representative unit types of two-dimensional and three-dimensional are given:

[0043] First, for a two-dimensional model, such as Figure 2 and Figure 3 shown as follows:

[0044] Step 1. Determine that the representative unit type is a two-dimensional model with length L = 100 and width W = 100; determine the particle equivalent radius r p = 1.5, particle volume fraction V f1 = 15%; fiber equivalent radius r f = 0.2, fiber volume fraction V f2 = 0.03, fiber length range (l min = 10, l max = 40).

[0045] Step 2: In the finite element software Abaqus, establish a two-dimensional matrix model part-base with the center at (0, 0) and the length and width of (100, 100), as shown in Figure 2 (a) in

[0046] Step 3: According to the equivalent radius of the particle 1.5 and the particle volume fraction 0.15, calculate the initial number of particles. For the two-dimensional model, the number of particles is

[0047]

[0048] where the multiplication coefficient is taken as 1.8.

[0049] Step 4: According to the representative unit size parameters (100, 100) and the number of particles 548, use Python language for secondary development to establish a Voronoi polygon structure with the center at (0, 0), the length and width of (118, 118), and the number of cells 548 in the finite element software Abaqus, as shown in Figure 2 (b) in, and generate 548 cells in the area of (118×118).

[0050] Step 5: Scale down the Voronoi polygon structure proportionally to form a particle group. For the two-dimensional model, the scaling coefficient of the side length is

[0051]

[0052] And:

[0053] 5.1 Remove the particles outside the matrix, as shown in Figure 2 (c) in;

[0054] 5.2 Remove the particles with most of their area outside the matrix at the boundary;

[0055] 5.2 Fine-tune the positions of other particles at the matrix boundary to prevent fine particle fragments from appearing at the boundary;

[0056] 5.3 Remove the part of the particle group outside the matrix range to make it all inside the matrix, and cut off the part of the particle group outside the matrix to obtain the initial particle group part-particles, as shown in Figure 2 (d) in.

[0057] 5.4 Generate an initial particle model with the same distribution and the same area of particles as circular according to the centroid position and area of the particles, as shown in Figure 3 (b) in.

[0058] Step 6: According to the representative unit size parameters, (100, 100), fiber equivalent radius 0.2, fiber volume fraction 0.03, fiber length range (10, 40); use Python language for secondary development, and adopt the random placement method to establish a fiber group model part-fibers centered at (0, 0) in the finite element software Abaqus, as Figure 3 shown in (a) of

[0059] Step 7: In the finite element software, subtract the fiber model from the initial particle model, that is, remove the particles that intersect with the fibers, as Figure 3 shown in (c) of Figure 3 ; and remove the broken particles and redundant particles to obtain the target particle model part-particles-re that meets the set parameters, as shown in (d) of

[0060] Figure 3 Set the target value to 15% to obtain the target particle model. Figure 3 Step 8: Combine the fiber group model part-fibers, the target particle model part-particles-re and the matrix model part-base to form a particle / fiber hybrid reinforced model that meets the set parameters, as Figure 3 shown in (e) of

[0061] For the two-dimensional and three-dimensional models, as Figure 4 and Figure 5 shown:

[0062] Step 1: Determine that the representative unit type is a three-dimensional model, with length L = 50, width W = 50, and height H = 50; determine the particle equivalent radius r p = 4, particle volume fraction V f1 = 15%; fiber equivalent radius r f = 1, fiber volume fraction V f2 = 0.01, fiber length range (l min = 10, l max = 35).

[0063] Step 2: In the finite element software Abaqus, establish a two-dimensional matrix model part-base centered at (0, 0, 0) with length, width, and height of (50, 50, 50), as Figure 4 shown in (a) of

[0064] Step 3: According to the particle equivalent radius 4 and the particle volume fraction 0.15; calculate the number of initial particles. The number of particles in the three-dimensional model is

[0065]

[0066] Among them, the multiplication coefficient takes the value of 1.8.

[0067] Step 4: According to the representative unit size parameters (100, 100, 100); the number of particles is 561. Use Python language for secondary development to establish a Voronoi polygon structure centered at (0, 0, 0) with a length, width, and height of (82, 82, 82) and 561 cells in the finite element software Abaqus.

[0068] Step 5: Scale down the Voronoi polygon structure proportionally to form a particle group. As shown in (b) of [], generate 561 cells in the area of (82×82×82); remove the outermost cells, and scale down the original cells by a factor of 0.653 to form a particle group. For a two-dimensional model, the scaling factor for the side length is Figure 4 shown in [], and:

[0069]

[0070] And:

[0071] 5.1 Remove the particles outside the matrix.

[0072] 5.2 Remove the particles with most of their area outside the matrix at the boundary;

[0073] 5.2 Fine-tune the positions of other particles at the matrix boundary to prevent small particle fragments from appearing at the boundary;

[0074] 5.3 Remove the part of the particle group outside the matrix range to obtain the initial particle group part-particles, as shown in (c) of []. Figure 4 shown in [].

[0075] 5.4 Generate an initial particle model with the same distribution and the same area of particles as circular according to the centroid position and area of the particles, as shown in (d) of []. Figure 4 shown in [].

[0076] Step 6: According to the representative unit size parameters, (50, 50, 50), the equivalent radius of the fiber is 1, the fiber volume fraction is 0.01, and the fiber length range is (10, 35); use Python language for secondary development and adopt the random placement method to establish a fiber group model part-fibers centered at (0, 0, 0) in the finite element software Abaqus, as shown in (a) of []. Figure 5 shown in [].

[0077] Step 7. In the finite element software, subtract the initial particle model from the fiber model, and remove the broken particles and redundant particles to obtain the target particle model part-particles-re that meets the set parameters. That is, remove the particles that intersect with the fibers and randomly remove particles until the particle volume fraction reaches the set target value of 15%, to obtain the target particle model, as shown in Figure 5 (b) in

[0078] Step 8. Combine the fiber group model part-fibers, the target particle model part-particles-re, and the matrix model part-base to form a particle / fiber hybrid reinforced model that meets the set parameters, as shown in Figure 5 (c) in. Perform the processing steps from Figure 5 (b) to (c) in on the initial particle model with circular particles to obtain a particle / fiber hybrid reinforced model with polyhedral particles, as shown in Figure 5 (d) in

[0079] This specific embodiment is only an interpretation of the present invention and is not a limitation thereof. Those skilled in the art can make modifications without creative contributions to this embodiment according to needs after reading this specification, but as long as it is within the scope of the claims of the present invention, it is protected by the patent law.

Claims

1. A method for constructing a particle-fiber hybrid reinforcement model, characterized in that, It includes the following steps: S1. Determine the representative unit type and its corresponding dimensional parameters, and determine the particle parameters and fiber parameters. Among them, the representative unit type is two-dimensional or three-dimensional; the particle parameters include the particle equivalent radius r p , the particle volume fraction V f1 ; the fiber parameters include the fiber equivalent radius r f , the fiber volume fraction V f2 , the fiber length range (l min , l max ); the representative unit dimensional parameters include the length L, width W, and height H, which are [(L + 12r p ), (W + 12r p )] corresponding to the two-dimensional representative unit or [(L + 8r p ), (W + 8r p ), (H + 8r p )] corresponding to the three-dimensional representative unit; S2. Establish a matrix model in the finite element software Abaqus; S3. Calculate the initial number of particles according to the particle parameters. Specifically, the calculation of the initial number of particles is as follows: for a two-dimensional model, The number of particles n p-2 is rounded; for a three-dimensional model, The number of particles n p-3 is rounded; f s is a multiplication factor, and its value range is (1.8, 2); S4. According to the representative unit size parameters and the number of particles, establish the corresponding Thiessen polygon structure in the finite element software Abaqus, and remove the outermost cells of the polygon structure; S5. Scale the Thiessen polygon structure proportionally to form an initial particle model. Specifically, scaling the Thiessen polygon structure proportionally to form an initial particle model means: scaling down the Thiessen polygon structure to form a particle swarm. For a two-dimensional model, the scaling factor For a three-dimensional model, the scaling factor Remove the particles outside the matrix, and remove the particles with 75% of their area at the matrix boundary outside the matrix; remove the part of the particle swarm outside the matrix range to obtain the initial particle model. Generate an initial particle model with particles of the same distribution and the same volume and in the shape of circles / spheres according to the centroid position and volume of the particles. S6. According to the representative unit size parameters and fiber parameters, establish a fiber model in the finite element software Abaqus by using the random placement method; S7. In the finite element software Abaqus, subtract the initial particle model from the fiber model, and remove the broken particles and redundant particles to obtain a target particle model that meets the set parameters; S8. Merge the fiber model, the target particle model and the matrix model to form a particle / fiber hybrid reinforced model that meets the set parameters.

2. The construction method of the particle-fiber hybrid reinforcement model according to claim 1, characterized in that: The matrix model is a two-dimensional model with the center at (0, 0), length and width of (L, W) established corresponding to the two-dimensional representative unit, or a three-dimensional model with the center at (0, 0), length, width and height of (L, W, H) established corresponding to the three-dimensional representative unit.

Images