Multi-angle continuous fiber reinforced composite material parametric modeling method

Abstract

The invention provides a parametric modeling method for a multi-angle continuous fiber reinforced composite material, and relates to the technical field of composite material simulation, and the parametric modeling method comprises the following steps: carrying out internal matrix modeling under an original coordinate system; designing aggregates in three shapes, and modeling in the internal matrix model to generate three pores; performing fiber modeling to obtain a fiber body external connection cube, and performing external matrix modeling to obtain an external matrix model; performing Boolean operation based on the internal matrix model, the aggregate, the fiber body external connection cube and the external matrix model to obtain a final model; and performing parametric programming on the final model by utilizing Abaqus, and calculating to obtain the boundary dimension. According to the method, modeling can be carried out at any angle within the range, the cutting mechanism of the continuous fiber reinforced composite material with the pore defect at any fiber angle is researched, and more accurate workpiece modeling and a smoother simulation calculation process are achieved.

Description

Technical Field

[0001] The invention relates to the technical field of composite material simulation, and in particular to a parameterized modeling method of a multi-angle continuous fiber reinforced composite material. Background Art

[0002] Since the 1970s, continuous fiber reinforced ceramic matrix composites (FRCMCs) have attracted attention from various industries since their introduction. This excellent high-temperature structural material is widely used in special fields and occasions that are beyond the reach of existing traditional materials due to its excellent strength, hardness, thermal conductivity, high temperature resistance, wear resistance, chemical corrosion resistance and many other properties. FRCMCs, especially as high-tech strategic materials, have made great strides in the aerospace field. According to the production method and actual use, a variety of ceramic matrix composites have been born, the most common of which are fiber-reinforced ceramics, nano-ceramics, gradient functional ceramics, in-situ growth ceramics, heterogeneous particle dispersion strengthened ceramics, etc.

[0003] Fiber-reinforced ceramic matrix composites are mainly composed of three phases: matrix, interface and fiber. The matrix phase plays the main role of supporting and transferring loads, supporting the overall structure of the material, protecting and isolating the fibers, and transferring loads to the fibers. The fiber phase, as a reinforcing phase, mainly plays the role of strengthening and toughening. The interface phase, as a bridge, mainly combines the matrix and the fiber. If the interface phase is too strong, it will lead to brittle fracture or even overall planar fracture. If the interface phase is too weak, it will lead to interface debonding and fiber pullout. Therefore, appropriate interface strength is crucial in simulation calculations.

[0004] In view of the complex actual service site, FRCMCs must present specific service characteristics according to specific service conditions. Therefore, a variety of FRCMCs are derived according to the type of matrix phase. They are mainly silicon carbide ceramic matrix composites, ultra-high temperature ceramic matrix composites and oxide ceramic matrix composites. Taking ceramic matrix composites as an example, the most common ones are carbon fiber reinforced silicon carbide ceramic matrix composites (C f / SiC) and silicon carbide fiber reinforced silicon carbide ceramic matrix composites (SiC f / SiC). f / SiC and SiC f / SiC are both important materials in the aerospace field, but at the same time, both are highly brittle materials. The pore defects caused by their preparation process and the difficult processing characteristics of the materials themselves limit their processing methods, resulting in problems such as difficult to control the processing surface, severe tool wear, and low processing efficiency. Therefore, many scholars have conducted a lot of research on how to effectively process silicon carbide ceramic-based composites. In a systematic research program that combines traditional processing with special processing and verifies experimental results with simulation processing, in addition to analyzing the influence of processing parameters on cutting force and surface roughness, they also mainly study the cutting force, surface quality and cutting mechanism under different fiber angles.

[0005] However, most of the fiber angles in the current simulation of silicon carbide ceramic matrix composite cutting are concentrated on classic angles such as 0°, 45° and 90°. For example, the fiber only rotates around the x-axis, resulting in a single fiber angle direction. Therefore, the development of an efficient and accurate modeling method for arbitrary fiber orientation angles in continuous fiber reinforced composites has become an urgent need in the industry. Summary of the invention

[0006] In order to overcome the shortcomings of the prior art, the purpose of the present invention is to provide a multi-angle parametric modeling method for continuous fiber reinforced composites, which can be modeled at any angle within a range, and can randomly generate pore defects inside the matrix to study the cutting mechanism of continuous fiber reinforced composites containing pore defects at any fiber angle, thereby achieving more accurate workpiece modeling and smoother simulation calculation process, and effectively improving the accuracy and work efficiency of FRCMCs in finite element simulation.

[0007] To achieve the above object, the present invention provides the following scheme: a multi-angle parameterized modeling method of continuous fiber reinforced composite materials, comprising the following steps:

[0008] Constructing an original coordinate system, and performing internal matrix modeling in the original coordinate system to obtain an internal matrix model;

[0009] Designing three shapes of aggregates, namely, spherical, ellipsoidal and random polyhedral, and modeling in the internal matrix model based on the three shapes of aggregates to generate three shapes of pores, namely, spherical, ellipsoidal and random polyhedral;

[0010] In the original coordinate system, the fiber angle is set, the body diagonal of the internal matrix body model is set as the outer dimension of the cube inscribed outside the fiber body, the number of fibers is calculated and the fiber modeling is performed to obtain a fiber body model;

[0011] In the original coordinate system, the body diagonal of the cube inscribed outside the fiber body is set as the outer dimension of the external matrix, and the external matrix is ​​modeled to obtain an external matrix model;

[0012] Performing Boolean operations based on the internal matrix model, aggregate, the fiber body model and the external matrix model to obtain a final model;

[0013] Abaqus is used to perform parameter programming on the final model, and the external dimensions of the final model are calculated.

[0014] Optionally, three shapes of aggregates, namely, spherical, ellipsoidal and random polyhedral, are designed, and based on the three shapes of aggregates, modeling is performed in the internal matrix model to generate pores in three shapes, namely, spherical, ellipsoidal and random polyhedral, including:

[0015] Setting the outer dimensions of the internal matrix model to the region generation boundary of the aggregate center point;

[0016] Calculating the distance between the aggregate particle size and the region generation boundary, judging whether the aggregate interferes with the region generation boundary, and judging that the result is no when the distance between the aggregate particle size and the region generation boundary is greater than the aggregate particle size;

[0017] Calculate the distance between any two aggregate center points to determine whether there is interference between the aggregates. When the distance between any two aggregate center points is greater than the aggregate particle size, the judgment result is no.

[0018] Randomly setting the angle α and the angle β, and calculating the vertex coordinates of the aggregate based on the angle α and the angle β;

[0019] Among them, the particle size of spherical bone particles is the radius of the sphere, the particle size of ellipsoidal bone particles is the major axis size of the ellipsoid, and the particle size of random polyhedron-shaped bone particles is the radius of the circumscribed sphere of the random polyhedron.

[0020] Optionally, for random polyhedral bone particles, the inner product of the plane unit normal vector of the aggregate and all points in the plane is calculated, and it is determined whether the aggregate is a convex polyhedron. If the inner product is a non-negative value, the determination result is yes.

[0021] Optionally, performing Boolean operations based on the internal matrix model, aggregate, the fiber body model and the external matrix model to obtain a final model includes:

[0022] Circularly assembling the aggregate and the internal matrix model in the original coordinate system and performing Boolean operations to obtain an internal matrix containing pores, thereby completing a first assembly operation;

[0023] In the original coordinate system, the body center coordinates of the internal matrix model are set as the coordinate origin and a second coordinate system is constructed;

[0024] Assembling the internal matrix model and the external matrix model in the second coordinate system and performing Boolean operations to obtain a hollow external matrix, thereby completing a second assembly operation;

[0025] Assembling the porous internal matrix and the fiber body model to perform Boolean operation to obtain a porous hollow internal matrix, thereby completing the third assembly;

[0026] The fiber body model and the hollow external matrix are reassembled to perform Boolean operations to obtain a final model, thereby completing the fourth assembly.

[0027] Optionally, in the first assembly operation, the aggregate is a cutting body and the internal matrix model is a cut body.

[0028] Optionally, in the second assembly operation, the external matrix model is moved until the center of the external matrix model coincides with the coordinate origin of the second coordinate system to perform the second assembly; the internal matrix model is the cutting body, and the external matrix model is the cut body.

[0029] Optionally, in the third assembly operation, the fiber body model is moved until the center of the fiber body model coincides with the coordinate origin of the second coordinate system, and the fiber body model is rotated around the Y-axis and Z-axis directions by arbitrary angles to perform the third assembly; the fiber body model is the cutting body, and the porous internal matrix is ​​the cut body.

[0030] Optionally, in the fourth assembly operation, repeat the step of "moving the fiber body model until the center of the fiber body model coincides with the coordinate origin of the second coordinate system, and rotating the fiber body model successively around the Y-axis and Z-axis directions by any angles" to perform the fourth assembly; the hollow external matrix is ​​the cutting body, and the fiber body model is the cut body.

[0031] Optionally, it is characterized in that the final model is parametrically programmed using Abaqus, and the dimensions of the final model are calculated, including:

[0032] Define the length, width and height of the workpiece, and set the fiber diameter, fiber angle elevation, azimuth and porosity of the workpiece as input parameters;

[0033] The body diagonal of the internal matrix is ​​calculated using the length, width and height of the workpiece to obtain the outer dimensions of the cube inscribed outside the fiber body;

[0034] Defining the spacing between fibers and the spacing between fibers and the boundary of a cube inscribed outside the fiber body, and calculating the number of fibers according to the defined spacing, and then rounding up the number of fibers to obtain an integer-level number of fibers;

[0035] The actual dimensions of the cube inscribed outside the fiber body are calculated according to the integer-level fiber number and the input parameters, and the dimensions of the external matrix are calculated using the actual dimensions. The calculation formula for the dimensions is:

[0036]

[0037] Among them, Length, Width, and Height are the length, width, and height of the workpiece respectively. 1 is the body diagonal of the internal matrix, is the distance between fibers and the distance between fibers and the boundaries of the cube inscribed outside the fiber body, m is the number of fibers, m' is the rounded-up value of m, and L' 1 is the actual size of the cube circumscribing the fiber body, L 2 The dimensions of the external substrate.

[0038] The present invention discloses the following technical effects by providing a multi-angle parameterized modeling method for continuous fiber reinforced composite materials:

[0039] 1. The present invention can be based on the secondary development of Abaqus in Python language, and the parametric modeling of continuous fiber reinforced composite materials is taken as the research object. The fiber composite angles at two angles, azimuth (rotation around the Z axis) and elevation (rotation around the Y axis), are defined according to the three-dimensional coordinate system to meet the modeling of reinforced fibers at any angle within the range of 0° to 360°. At the same time, pore defects are randomly generated inside the matrix, which is convenient for studying the cutting mechanism of continuous fiber reinforced composite materials containing pore defects at any fiber angle.

[0040] 2. Due to the differences in the preparation processes of continuous fiber reinforced composite materials and the differences in the actual preparation environments, the present invention provides three types of pores, namely, spherical, ellipsoidal and random polyhedral shapes, to facilitate the comparison of simulation results.

[0041] 3. The present invention can combine Abaqus's rich built-in modules and the concise Python programming language to establish a more accurate workpiece model and a smoother simulation calculation process, effectively improving the accuracy of FRCMCs in finite element simulation, thereby reducing the number of experiments. Especially when multiple sets of simulation results are required for comparative analysis, the parametric modeling idea simplifies the traditional simulation modeling process, improves the efficiency of finite element simulation, shortens the simulation cycle, saves manpower and material resources, and enriches the researchability of the fiber cutting mechanism at any angle.

[0042] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0044] Figure 1 A schematic diagram of a method flow chart provided by an embodiment of the present invention;

[0045] Figure 2 Schematic diagram of three shapes of aggregates and internal matrix containing pores provided in embodiments of the present invention;

[0046] Figure 3 A schematic diagram of the fiber angle of a workpiece provided by an embodiment of the present invention;

[0047] Figure 4 A schematic diagram of workpiece assembly provided by an embodiment of the present invention;

[0048] Figure 5 A schematic diagram of a plug-in input interface provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0049] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0050] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0051] like Figure 1-4 As shown, the present invention provides a multi-angle parameterized modeling method for continuous fiber reinforced composite materials, comprising the following steps:

[0052] 1. Construct an original coordinate system, and perform internal matrix modeling in the original coordinate system, taking the external dimensions of the internal matrix as the final external dimensions of the workpiece matrix model to obtain the internal matrix model.

[0053] 2. Aggregate modeling is performed in the internal matrix area. Due to differences in manufacturing processes, the shape of the pores formed is also uncertain. For the authenticity of the results, three types of aggregates, spherical, ellipsoidal and random polyhedral, are designed. Based on the three types of aggregates, aggregate modeling is performed in the internal matrix model to generate pores (including spherical, ellipsoidal and random polyhedral shapes) that match the internal matrix. Specifically, it includes:

[0054] To ensure that the aggregate is accurately placed in the internal matrix area, the external dimensions of the internal matrix model are set to the area generation boundary of the center point of the aggregate; in order to avoid interference between the aggregate and the boundary, this embodiment is implemented through the following steps:

[0055] 1) Calculate the distance between the aggregate particle size and the region generation boundary to determine whether the aggregate interferes with the region generation boundary. When the distance between the aggregate particle size and the region generation boundary is greater than the aggregate particle size, the judgment result is no. 2) Calculate the distance between any two aggregate center points to determine whether the aggregates interfere with each other. When the distance between any two aggregate center points is greater than the aggregate particle size, the judgment result is no.

[0056] Generate a random coordinate P(X p ,Y p ,Z p ), then the judgment formula of boundary interference condition is as follows:

[0057] d h / 2+Length / 100≤X p ≤Length-d h / 2-Length / 100

[0058] d h / 2+Height / 100≤Y p ≤Height-d h / 2-Height / 100

[0059] d h / 2+Width / 100≤Z p ≤Widthd h / 2-Width / 100

[0060] Randomly set the angle α and the angle β, and calculate the spherical vertex coordinates V of the polyhedral aggregate based on the angle α and the angle β i (X i ,Y i ,Z i ):

[0061]

[0062] Among them, the particle size of spherical bone particles is the radius of the sphere, the particle size of ellipsoidal bone particles is the major axis size of the ellipsoid, and the particle size of random polyhedron-shaped bone particles is the radius of the circumscribed sphere of the random polyhedron.

[0063] For random polyhedral bone particles, in order to determine that the generated aggregates are all convex polyhedrons, the inner product of the plane unit normal vector of the aggregate and all points in the plane is calculated, and it is judged whether the aggregate is a convex polyhedron. If the inner product is a non-negative value, the judgment result is yes.

[0064] 3. In the original coordinate system, set the fiber angle, set the body diagonal of the internal matrix body model to the outer dimensions of the cube inscribed outside the fiber body, calculate the number of fibers and perform fiber modeling to obtain a fiber body model. Specifically include:

[0065] The fibers are modeled in the original coordinate system. To ensure that the fibers can be evenly distributed in the right-angle edge area of ​​the internal matrix after being assigned an angle, the fiber angle is set, and the body diagonal of the internal matrix model is set to the outer dimensions of the cube inscribed outside the fiber body so that the fibers are evenly distributed inside the internal matrix. The number of fibers is then calculated, and fiber modeling is performed to obtain a fiber body model.

[0066] 4. In the original coordinate system, ensure that the external matrix can completely cut off the excess part of the fiber, set the body diagonal of the cube inscribed outside the fiber body as the external dimension of the external matrix, and model the external matrix to obtain the external matrix model.

[0067] 5. Perform Boolean operations based on the internal matrix model, aggregate, the fiber body inscribed cube and the external matrix model to obtain the final model. Specifically including:

[0068] S1. The first assembly is to cyclically assemble the polyhedral aggregate and the internal matrix to perform Boolean operations. In the original coordinate system, the aggregate is used as the cutting body and the internal matrix is ​​used as the cut body to obtain the internal matrix containing pores.

[0069] S2. In the original coordinate system, the body center coordinates of the internal matrix model are set as the coordinate origin and a second coordinate system is constructed;

[0070] S3, the second assembly, assembling the internal matrix model and the external matrix model in the second coordinate system for Boolean operation. Move the external matrix model until the center of the external matrix model coincides with the origin of the second coordinate system; the internal matrix model is the cutting body, and the external matrix model is the cut body, to obtain a hollow external matrix;

[0071] S4, the third assembly, assembling the porous internal matrix and the fiber for Boolean operation. Move the fiber body model until the center of the fiber body model coincides with the coordinate origin of the second coordinate system, and rotate the fiber body model around the Y-axis and Z-axis directions by any angle; the fiber body model is the cutting body, and the porous internal matrix is ​​the cut body, to obtain a porous hollow internal matrix, and complete the third assembly;

[0072] S5, four times of assembly, reassembling the fiber body model and the hollow external matrix for Boolean operation, it is necessary to repeat the steps of "moving the fiber body model until the center of the fiber body model coincides with the coordinate origin of the second coordinate system, and rotating the fiber body model around the Y-axis and Z-axis directions in turn by any angle" The hollow external matrix is ​​the cutting body, and the fiber body model is the cut body to obtain the final model.

[0073] 6. Use Abaqus to perform parametric programming on the final model and calculate the dimensions of the final model. Specifically including:

[0074] Define the workpiece dimensions length (Length), width (Width), height (Height), fiber diameter (d), fiber angle elevation angle (Angle1), azimuth angle (Angle2), and porosity (Porosity) as input parameters;

[0075] The internal matrix body diagonal can be calculated

[0076]

[0077] That is, the calculated external dimensions of the cube inscribed outside the fiber body;

[0078] Defines the spacing between fibers in the x and y directions and the spacing between fibers and the cube boundaries surrounding the fiber bodies

[0079] And further calculate the number of fibers m according to the defined spacing:

[0080]

[0081] Since the obtained fiber number is not necessarily an integer, the fiber number m needs to be rounded up to obtain an integer fiber number m'. The program code for rounding m up is m'=int(math.ceil(m))+1;

[0082] Further, the internal matrix diagonal L' is calculated based on the integer fiber number and the input parameters. 1 , that is, the actual size L' of the fiber circumscribed cube 1 :

[0083]

[0084] The outer dimensions L of the external substrate are calculated using the actual outer dimensions. 2 :

[0085]

[0086] 7. Plug-in input interface after programming is completed

[0087] like Figure 5 As shown in the figure, the plug-in interface mainly displays the relevant information of the workpiece, so the main title of the plug-in interface is "Workpiece". Specifically, it can be divided into two first-level title boxes, "Workpieceparameter" and "Diagram", namely the workpiece parameters and model diagrams, and the contents of these two first-level title boxes correspond to each other.

[0088] The first secondary title box under the "Workpiece parameter" title box is "Porosity Option", i.e., the pore option; in this secondary title box, set the "Yes" and "No" options, i.e., select whether to generate pores inside the matrix of the workpiece, "Yes" means to generate pores, and "No" means not to generate pores. The second secondary title box under the "Workpiece parameter" title box is "PorosityType", i.e., select the type of pores to be generated; in this secondary title box, set the three pore options of "Polyhedron", "Sphere" and "Ellipe", which represent random polyhedron, circular and elliptical pores respectively, and the shape diagrams of the three pores of "Polyhedron", "Sphere" and "Ellipe" are respectively displayed in the "Diagram" title box. The third secondary title box under the "Workpieceparameter" title box is "Parameter", which is the specific input parameter of the workpiece; "Length", "Height", "Width", "Fiber_D", "Angle1", "Angle2" and "Porosity" are set in this secondary title box, which respectively represent the length, width and height of the workpiece's overall dimensions, fiber diameter, fiber angle elevation, azimuth and porosity. The workpiece's overall dimensions and fiber diameter are in millimeters (mm), the fiber angle is in degrees (o), and the porosity is expressed in percentage (%). In the model diagram of the "Diagram" title box, there are feature identifiers corresponding to each input parameter.

[0089] The three buttons at the bottom of the plug-in interface are "Continue", "Defaults" and "Cancel". The "Continue" button means to continue the next operation, specifically jumping to the next operation interface or submitting the content of the current operation interface; the "Defaults" button means to restore the default values ​​of input parameters. The default values ​​of input parameters are generally recommended values. When the user modifies the default values ​​of input parameters, they can be restored through the "Defaults" button; the "Cancel" button means to cancel the current operation and not save the current operation.

[0090] Therefore, the present invention provides a multi-angle parametric modeling method for continuous fiber reinforced composites, which can be modeled at any angle within a range, and can randomly generate pore defects inside the matrix to study the cutting mechanism of continuous fiber reinforced composites containing pore defects at any fiber angle, thereby achieving more accurate workpiece modeling and smoother simulation calculation process, and effectively improving the accuracy and work efficiency of FRCMCs in finite element simulation.

[0091] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0092] The principles and implementation methods of the present invention are described in this article using specific examples. The description of the above embodiments is only used to help understand the method and core idea of ​​the present invention. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A multi-angle parameterized modeling method for continuous fiber reinforced composite materials, characterized in that: The following steps are involved: Constructing an original coordinate system, and performing internal matrix modeling in the original coordinate system to obtain an internal matrix model; Designing three shapes of aggregates, namely, spherical, ellipsoidal and random polyhedral, and modeling in the internal matrix model based on the three shapes of aggregates to generate three shapes of pores, namely, spherical, ellipsoidal and random polyhedral; In the original coordinate system, the fiber angle is set, the body diagonal of the internal matrix body model is set as the outer dimension of the cube inscribed outside the fiber body, the number of fibers is calculated and the fiber modeling is performed to obtain a fiber body model; In the original coordinate system, the body diagonal of the cube inscribed outside the fiber body is set as the outer dimension of the external matrix, and the external matrix is ​​modeled to obtain an external matrix model; Performing Boolean operations based on the internal matrix model, aggregate, the fiber body model and the external matrix model to obtain a final model; Abaqus is used to perform parameter programming on the final model, and the external dimensions of the final model are calculated.

2. A multi-angle parameterized modeling method for continuous fiber reinforced composite materials according to claim 1, characterized in that: Aggregates in three shapes, spherical, ellipsoidal and random polyhedral, are designed, and based on the three shapes of aggregates, modeling is performed in the internal matrix model to generate pores in three shapes, spherical, ellipsoidal and random polyhedral, including: Setting the outer dimensions of the internal matrix model to the region generation boundary of the aggregate center point; Calculating the distance between the aggregate particle size and the region generation boundary, judging whether the aggregate interferes with the region generation boundary, and judging that the result is no when the distance between the aggregate particle size and the region generation boundary is greater than the aggregate particle size; Calculate the distance between any two aggregate center points to determine whether there is interference between the aggregates. When the distance between any two aggregate center points is greater than the aggregate particle size, the judgment result is no. Randomly setting the angle α and the angle β, and calculating the vertex coordinates of the aggregate based on the angle α and the angle β; Among them, the particle size of spherical bone particles is the radius of the sphere, the particle size of ellipsoidal bone particles is the major axis size of the ellipsoid, and the particle size of random polyhedron-shaped bone particles is the radius of the circumscribed sphere of the random polyhedron.

3. A multi-angle parameterized modeling method for continuous fiber reinforced composite materials according to claim 2, characterized in that: For random polyhedral bone particles, the inner product of the plane unit normal vector of the aggregate and all points in the plane is calculated, and it is determined whether the aggregate is a convex polyhedron. If the inner product is a non-negative value, the determination result is yes.

4. The multi-angle parameterized modeling method of continuous fiber reinforced composite materials according to claim 3, characterized in that: A Boolean operation is performed based on the internal matrix model, aggregate, the fiber body model and the external matrix model to obtain a final model, including: Circularly assembling the aggregate and the internal matrix model in the original coordinate system and performing Boolean operations to obtain an internal matrix containing pores, thereby completing a first assembly operation; In the original coordinate system, the body center coordinates of the internal matrix model are set as the coordinate origin and a second coordinate system is constructed; Assembling the internal matrix model and the external matrix model in the second coordinate system and performing Boolean operations to obtain a hollow external matrix, thereby completing a second assembly operation; Assembling the porous internal matrix and the fiber body model to perform Boolean operation to obtain a porous hollow internal matrix, thereby completing the third assembly; The fiber body model and the hollow external matrix are reassembled to perform Boolean operations to obtain a final model, thereby completing the fourth assembly.

5. The multi-angle parameterized modeling method of continuous fiber reinforced composite materials according to claim 4, characterized in that: In the first assembly operation, the aggregate is the cutting body and the internal matrix model is the cut body.

6. A multi-angle parameterized modeling method for continuous fiber reinforced composite materials according to claim 5, characterized in that: In the second assembly operation, the external matrix model is moved until the center of the external matrix model coincides with the coordinate origin of the second coordinate system to perform the second assembly; the internal matrix model is the cutting body, and the external matrix model is the cut body.

7. A multi-angle parameterized modeling method for continuous fiber reinforced composite materials according to claim 6, characterized in that: In the third assembly operation, the fiber body model is moved until the body center of the fiber body model coincides with the coordinate origin of the second coordinate system, and the fiber body model is rotated around the Y-axis and Z-axis directions by arbitrary angles to perform the third assembly; the fiber body model is the cutting body, and the porous internal matrix is ​​the cut body.

8. The multi-angle parameterized modeling method of continuous fiber reinforced composite materials according to claim 7, characterized in that: In the fourth assembly operation, the step of "moving the fiber body model until the center of the fiber body model coincides with the coordinate origin of the second coordinate system, and rotating the fiber body model successively around the Y-axis direction and the Z-axis direction by any angle" is repeated to perform the fourth assembly; the hollow external matrix is ​​the cutting body, and the fiber body model is the cut body.

9. A multi-angle parameterized modeling method for continuous fiber reinforced composite materials according to claim 8, characterized in that: Abaqus is used to perform parametric programming on the final model and calculate the dimensions of the final model, including: Define the length, width and height of the workpiece, and set the fiber diameter, fiber angle elevation, azimuth and porosity of the workpiece as input parameters; The body diagonal of the internal matrix is ​​calculated using the length, width and height of the workpiece to obtain the outer dimensions of the cube inscribed outside the fiber body; Defining the spacing between fibers and the spacing between fibers and the boundary of a cube inscribed outside the fiber body, and calculating the number of fibers according to the defined spacing, and then rounding up the number of fibers to obtain an integer-level number of fibers; The actual dimensions of the cube inscribed outside the fiber body are calculated according to the integer-level fiber number and the input parameters, and the dimensions of the external matrix are calculated using the actual dimensions. The calculation formula for the dimensions is: Among them, Length, Width, and Height are the length, width, and height of the workpiece respectively, and L1 is the body diagonal of the internal matrix. is the distance between fibers and between fibers and between fibers and the boundary of the cube bounded by the fiber body, m is the number of fibers, m′ is the rounded-up value of m, L′1 is the actual size of the cube bounded by the fiber body, and L2 is the size of the external matrix.

Images