A multi-angle parametric modeling method of continuous fiber reinforced composites

Abstract

The application provides a multi-angle parametric modeling method of continuous fiber reinforced composite material, and relates to the technical field of composite material simulation, and comprises the following steps: internal matrix modeling is carried out in an original coordinate system; three kinds of aggregate shapes are designed, and three kinds of pores are generated by modeling in the internal matrix model; fiber modeling is carried out to obtain a fiber body circumscribed cube, external matrix modeling is carried out to obtain an external matrix model; Boolean operation is carried out based on the internal matrix model, the aggregate, the fiber body circumscribed cube and the external matrix model to obtain a final model; Abaqus is used for parametric programming of the final model, and the size of the final model is calculated. The application can model at any angle within a range, study the cutting mechanism of the continuous fiber reinforced composite material containing a pore defect at any fiber angle, and realize more accurate workpiece modeling and more smooth simulation calculation process.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of composite material simulation, in particular to a multi-angle parameterized modeling method of continuous fiber reinforced composite material. BACKGROUND

[0002] Since the 1970s, continuous fiber reinforced ceramic matrix composites (FRCMCs) have been concerned by all industries since their advent. This excellent high-temperature structural material is widely used in special fields and occasions that cannot be achieved by existing traditional materials due to its excellent strength, hardness, thermal conductivity, high-temperature resistance, wear resistance, and chemical corrosion resistance. FRCMCs, especially as high-tech strategic materials, have made great strides in the field of aerospace. According to the production method and actual use, various ceramic matrix composites are born, the most common of which are fiber reinforced ceramics, nano ceramics, gradient functional ceramics, in-situ grown ceramics, and heterogeneous particle dispersion strengthened ceramics.

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

[0004] For complex actual service sites, FRCMCs must exhibit specific service characteristics according to specific service conditions, so various FRCMCs are derived according to the type of the matrix phase. Mainly including 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). C f / SiC and SiC fBoth SiC and Si3N4 are important materials in the field of aerospace, but both of them are high brittle materials. The preparation process of the two materials produces pore defects and the characteristics of the materials themselves are difficult to process, which limits the processing methods, resulting in problems such as uncontrolled processing surface, serious tool wear and low processing efficiency. Therefore, many scholars have made a lot of research on how to effectively process silicon carbide ceramic matrix composites. In the system research scheme combining traditional processing and special processing, experimental results and simulation processing, in addition to analyzing the influence of processing parameters on cutting force and surface roughness, the cutting force, surface quality and cutting mechanism under different fiber angles are also studied.

[0005] However, the fiber angles in most current simulations of cutting of silicon carbide ceramic matrix composites are concentrated on classic angles such as 0°, 45° and 90°, such as the fiber only rotating around the x-axis, resulting in a single fiber angle direction. Therefore, it is an urgent need in the industry to develop a high-efficiency and accurate modeling method for continuous fiber reinforced composites with arbitrary fiber orientation angles. SUMMARY

[0006] In order to overcome the shortcomings of the prior art, the purpose of the present application is to provide a multi-angle parameterized modeling method for continuous fiber reinforced composites, which can model at any angle within a range, and can randomly generate pore defects in the matrix to study the cutting mechanism of continuous fiber reinforced composites with pore defects under any fiber angle, realize more accurate workpiece modeling and more smooth simulation calculation process, and effectively improve the precision and work efficiency of FRCMCs in finite element simulation.

[0007] To achieve the above purpose, the present application provides the following scheme: a multi-angle parameterized modeling method for continuous fiber reinforced composites, comprising the following steps:

[0008] An original coordinate system is constructed, and internal matrix modeling is performed in the original coordinate system to obtain an internal matrix model;

[0009] Three kinds of aggregate shapes of spherical, ellipsoidal and random polyhedral are designed, and based on the three kinds of aggregate shapes, modeling is performed in the internal matrix model to generate pores of spherical, ellipsoidal and random polyhedral shapes;

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

[0011] In the original coordinate system, the body diagonal of the fiber body circumscribed cube is set as the size of the outer shape of the external matrix, and external matrix modeling is performed to obtain an external matrix model;

[0012] performing a Boolean operation based on the internal matrix model, the aggregate, the fibrous body model and the external matrix model to obtain a final model;

[0013] performing parameterized programming on the final model by using Abaqus and calculating the size of the final model.

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

[0015] setting the size of the internal matrix model as a region generation boundary of a center point of the aggregate;

[0016] calculating the distance between the size of the aggregate and the region generation boundary to determine whether the aggregate interferes with the region generation boundary, and the result is no when the distance is greater than the size of the aggregate;

[0017] calculating the distance between any two center points of the aggregate to determine whether the aggregate interferes with each other, and the result is no when the distance is greater than the size of the aggregate;

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

[0019] wherein the size of the spherical aggregate is the radius of the spherical shape, the size of the ellipsoidal aggregate is the major axis size of the ellipsoidal shape, and the size of the random polyhedral aggregate is the radius of the circumscribed spherical shape of the random polyhedral shape.

[0020] Optionally, for the random polyhedral aggregate, 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, and the result is yes if the inner product is a non-negative value.

[0021] Optionally, performing a Boolean operation based on the internal matrix model, the aggregate, the fibrous body model and the external matrix model to obtain a final model, including:

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

[0023] under the original coordinate system, setting the body center coordinates of the internal matrix model as the coordinate origin and constructing a second coordinate system;

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

[0025] Assembling and performing Boolean operation on the internal matrix model containing pores and the fiber body model to obtain a hollow internal matrix containing pores, and completing a third assembling operation;

[0026] Re-assembling and performing Boolean operation on the fiber body model and the hollow external matrix to obtain a final model, and completing a fourth assembling operation.

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

[0028] Optionally, in the second assembling operation, the external matrix model is moved to the body center of the external matrix model, and the body center of the external matrix model is coincided with the coordinate origin of the second coordinate system, so as to perform the second assembling operation; the internal matrix model is a cutting body, and the external matrix model is a cut body.

[0029] Optionally, in the third assembling operation, the fiber body model is moved to the body center of the fiber body model, and the body center of the fiber body model is coincided with the coordinate origin of the second coordinate system, and the fiber body model is rotated by an arbitrary angle around the Y-axis direction and the Z-axis direction in sequence, so as to perform the third assembling operation; the fiber body model is a cutting body, and the internal matrix containing pores is a cut body.

[0030] Optionally, in the fourth assembling operation, the step of "moving the fiber body model to the body center of the fiber body model, and the body center of the fiber body model is coincided with the coordinate origin of the second coordinate system, and the fiber body model is rotated by an arbitrary angle around the Y-axis direction and the Z-axis direction in sequence" is repeated, so as to perform the fourth assembling operation; the hollow external matrix is a cutting body, and the fiber body model is a cut body.

[0031] Optionally, the final model is parameterized programmed by Abaqus, and the size of the final model is calculated, including:

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

[0033] Calculating the body diagonal of the internal matrix by using the length, width and height of the workpiece, to obtain the size of the circumscribed cube of the fiber body;

[0034] Defining the fiber and the distance between the fibers, and the distance between the fibers and the boundary of the circumscribed cube of the fiber body, and calculating the number of fibers according to the defined distance, and then rounding up the number of fibers to obtain the integer level number of fibers;

[0035] According to the integer level fiber quantity and the input parameter, the actual size of the fiber body circumscribed cube is calculated, and the size of the external matrix is calculated by using the actual size, and the calculation formula of the size is:

[0036]

[0037] Wherein, Length, Width, Height are the length, width and height of the workpiece, L1 is the body diagonal of the internal matrix, The distance between the fibers and the distance between the fibers and the boundary of the fiber body circumscribed cube, m is the fiber quantity, m' is the upward integer value of m, L'1 is the actual size of the fiber body circumscribed cube, L2 is the size of the external matrix.

[0038] The present application provides a kind of multi-angle continuous fiber reinforced composite parameterized modeling method, and the following technical effects are disclosed:

[0039] 1, the present application can be based on Python language Abaqus secondary development, to continuous fiber reinforced composite parameterized modeling as research object, according to three-dimensional coordinate system defines the fiber composite angle under two angles of azimuth angle (rotation around Z axis) and elevation angle (rotation around Y axis), meet the modeling of any angle of 0°~360° range of reinforcing fiber, while generating pore defects randomly in matrix, it is convenient to study the cutting mechanism of continuous fiber reinforced composite material containing pore defects under any fiber angle.

[0040] 2, due to the difference of continuous fiber reinforced composite material preparation process and the difference of actual preparation environment, the present application provides three kinds of pores of spherical, ellipsoidal and random polyhedral shape, which is convenient for comparison of simulation results.

[0041] 3, the present application can combine the rich built-in modules of Abaqus and the simple Python programming language, establish more accurate workpiece model, more smooth simulation calculation process, effectively improve the precision of FRCMCs in finite element simulation, thereby reducing the number of tests, especially when multiple simulation results are compared and analyzed, the parameterized 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 fiber cutting mechanism under any angle.

[0042] The technical solutions of the present application will be further described in detail below with the help of drawings and examples. DETAILED DESCRIPTION

[0043] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments. 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 any creative work.

[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 a porous internal matrix provided in an embodiment 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 the plug-in input interface provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0049] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. 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 making creative efforts 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 Figures 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. The external dimensions of the internal matrix are used as the final external dimensions of the workpiece matrix model to obtain the internal matrix model.

[0053] 2. In the internal matrix area, the shape of the formed pores is uncertain due to the difference of manufacturing process, in order to ensure the authenticity of the results, three kinds of shapes of aggregate, i.e. spherical, ellipsoidal and random polyhedral, are designed, and based on the three kinds of aggregate, the aggregate modeling is carried out in the internal matrix model to generate pores (including three kinds of shapes of spherical, ellipsoidal and random polyhedral) matching the internal matrix. Specifically, it comprises:

[0054] In order to ensure the accurate placement of aggregate in the internal matrix area, the size of the internal matrix model is set as the region generating boundary of the center point of the aggregate; in order to avoid the interference of the aggregate with the boundary, the embodiment is realized by the following steps:

[0055] 1) Calculate the distance between the aggregate particle size and the region generating boundary, and judge whether the aggregate interferes with the region generating boundary, when the distance between the aggregate particle size and the region generating boundary is greater than the aggregate particle size, the result is no; 2) Calculate the distance between any two aggregate center points, and judge whether the aggregate interferes with each other, when the distance between any two aggregate center points is greater than the aggregate particle size, the result is no.

[0056] A random coordinate P (X p ,Y p ,Z p ) is generated in the placement area, and the judgment formula of the 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 ≤Width-d h / 2-Width / 100

[0060] Wherein, d h represents the particle size of the aggregate, wherein the particle size of the spherical aggregate is the radius of the spherical shape, the particle size of the ellipsoidal aggregate is the long axis size of the ellipsoidal shape, and the particle size of the random polyhedral aggregate is the radius of the circumscribed spherical shape of the random polyhedral shape, Length is the length of the workpiece size, Width is the width of the workpiece size, and Height is the height of the workpiece size.

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

[0062]

[0063] Wherein, the particle size of the spherical bone particle is the radius of the spherical shape, the particle size of the ellipsoidal bone particle is the long axis size of the ellipsoidal shape, and the particle size of the random polyhedral bone particle is the radius of the circumscribed spherical shape of the random polyhedral shape.

[0064] For the random polyhedral bone particle, in order to determine that the generated aggregate is a convex polyhedron, 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.

[0065] 3. In the original coordinate system, set the fiber angle, set the body diagonal of the internal matrix model as the size of the circumscribed cube of the fiber body, calculate the number of fibers and perform fiber modeling to obtain a fiber body model. Specifically, it includes:

[0066] In the original coordinate system, the fiber is modeled, in order to ensure that the fiber can be uniformly distributed in the right-angle edge area of the internal matrix after being given an angle, by: setting the fiber angle, and setting the body diagonal of the internal matrix model as the size of the circumscribed cube of the fiber body, so that the fiber is uniformly distributed in the internal matrix, then calculating the number of fibers and performing fiber modeling to obtain a fiber body model.

[0067] 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 circumscribed cube of the fiber body as the size of the external matrix, and perform external matrix modeling to obtain an external matrix model.

[0068] 5. Perform Boolean operation based on the internal matrix model, aggregate, circumscribed cube of the fiber body, and external matrix model to obtain a final model. Specifically, it includes:

[0069] S1, first assembly, cyclically assemble the polyhedral aggregate and the internal matrix to perform Boolean operation. In the original coordinate system, take the aggregate as the cutting body and the internal matrix as the cut body to obtain an internal matrix containing pores.

[0070] S2, in the original coordinate system, set the body center coordinates of the internal matrix model as the coordinate origin and construct a second coordinate system;

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

[0072] S4, third assembly, assemble the internal matrix containing pores and the fiber to make Boolean operation. Move the fiber body model to the body center of the fiber body model, which coincides with the coordinate origin of the second coordinate system, and rotate the fiber body model by an arbitrary angle around the Y-axis direction and the Z-axis direction in turn; the fiber body model is a cutting body, and the internal matrix containing pores is a cut body, to obtain a hollow internal matrix containing pores, and complete the third assembly;

[0073] S5, fourth assembly, reassemble the fiber body model and the hollow external matrix to make Boolean operation, which needs to repeat the step of "moving the fiber body model to the body center of the fiber body model, which coincides with the coordinate origin of the second coordinate system, and rotating the fiber body model by an arbitrary angle around the Y-axis direction and the Z-axis direction in turn"; the hollow external matrix is a cutting body, and the fiber body model is a cut body, to obtain a final model.

[0074] 6, parameterize the final model by using Abaqus, and calculate the size of the final model. Specifically, it includes:

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

[0076] Calculate the body diagonal of the internal matrix

[0077]

[0078] That is, the calculated size of the circumscribed cube of the fiber body;

[0079] Define the distance between the fibers and the distance between the fibers and the boundary of the circumscribed cube of the fiber body in the x direction and the y direction

[0080] And further calculate the number of fibers m according to the defined distance:

[0081]

[0082] Wherein, d represents the fiber diameter.

[0083] 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;

[0084] The internal matrix diagonal L'1, that is, the actual outer dimension L'1 of the fiber circumscribed cube, is further calculated based on the integer-level fiber number and the input parameters:

[0085]

[0086] Here, d represents the fiber diameter.

[0087] The outer dimension L2 of the outer base is calculated using the actual outer dimension:

[0088]

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

[0090] like Figure 5 As shown in the figure, the plug-in interface mainly displays relevant information about the workpiece, so the main title of the plug-in interface is "Workpiece". Specifically, it can be divided into two first-level title boxes: "Workpiece parameter" and "Diagram", which correspond to the workpiece parameters and model diagram.

[0091] The first sub-title box under the title box of "Workpiece parameter" is "Porosity Option", that is, whether to generate pores in the interior of the base body of the workpiece; "Yes" means to generate pores, and "No" means not to generate pores. The second sub-title box under the title box of "Workpiece parameter" is "Porosity Type", that is, the type of generated pores; three pore options of "Polyhedron", "Sphere" and "Ellipe" are set in the sub-title box, which respectively represent random polyhedral, circular and elliptical pores, and the shape diagrams of the three pores of "Polyhedron", "Sphere" and "Ellipe" are respectively displayed in the "Diagram" title box. The third sub-title box under the title box of "Workpiece parameter" is "Parameter", that is, the specific input parameters of the workpiece; "Length", "Height", "Width", "Fiber_D", "Angle1", "Angle2" and "Porosity" are set in the sub-title box, which respectively represent the length, width and height of the workpiece, the fiber diameter, the fiber angle elevation angle, the azimuth angle and the porosity, the length and width of the workpiece and the fiber diameter are in millimeters (mm), the fiber angle is in degrees (o), and the porosity is in percentage (%). In the model diagram of the "Diagram" title box, there are feature marks corresponding to each input parameter.

[0092] 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, and the specific operations include jumping to the next operation interface or submitting the content of the current operation interface; the "Defaults" button means to restore the default value of the input parameter, and the default value of the input parameter is generally the recommended value; when the user modifies the default value of the input parameter, the "Defaults" button can be used to restore it; the "Cancel" button means to cancel the current operation and not to save the current operation.

[0093] Therefore, the present application can model at any angle within a range, and can randomly generate pore defects in the interior of the base body to study the cutting mechanism of continuous fiber reinforced composite materials containing pore defects at any fiber angle, so as to realize more accurate workpiece modeling and more smooth simulation calculation process, and effectively improve the precision and work efficiency of FRCMCs in finite element simulation.

[0094] The various embodiments described in this specification are presented for the purpose of illustrating the principles of the present application and its best mode of operation. Each of the embodiments described in this specification has been provided for the purpose of illustration and is not intended to limit the application.

[0095] The principles and implementations of the present application have been described in the specification with specific examples. The above description of the embodiments is only for the purpose of helping to understand the method of the present application and its core idea. Meanwhile, for those skilled in the art, the specific implementation and application range of the present application can be changed according to the idea of the present application. In summary, the content of the specification should not be understood as a limitation of the present application.

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 types of aggregates: spherical, ellipsoidal, and random polyhedron-shaped; and modeling the three types of aggregates in the internal matrix model to generate pores in the three types of spherical, ellipsoidal, and random polyhedron-shaped pores; In the original coordinate system, the fiber angle is set, the body diagonal of the internal matrix model is set as the outer dimensions of the cube inscribed outside the fiber body, the number of fibers is calculated and fiber modeling is performed to obtain a fiber body model; In the original coordinate system, the diagonal of the cube inscribed outside the fiber body is set as the outer dimensions 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; Performing parametric programming on the final model using Abaqus, and calculating the dimensions of the final model; A Boolean operation is performed based on the internal matrix model, the 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, wherein the aggregate is a cutting body and the internal matrix model is a cut body; 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 a Boolean operation to obtain a hollow external matrix, thereby completing a second assembly operation. 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, thereby performing the second assembly; the internal matrix model is a cutting body, and the external matrix model is a cut body; Assembling the porous internal matrix and the fiber body model and performing a Boolean operation to obtain a porous hollow internal matrix, thereby completing the third assembly. 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 the Z-axis by any angle to perform the third assembly; the fiber body model is the cutting body, and the porous internal matrix is ​​the cut body; The fiber body model and the hollow outer matrix are reassembled to perform Boolean operations to obtain a final model, thereby completing the fourth assembly.

2. A multi-angle parameterized modeling method for continuous fiber reinforced composite materials according to claim 1, characterized in that: For a random polyhedral aggregate, the inner product of the 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.

3. The multi-angle parameterized modeling method of continuous fiber reinforced composite materials according to claim 2, 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 and Z-axis directions 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.

4. The multi-angle parameterized modeling method of continuous fiber reinforced composite materials according to claim 3, 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; Calculating the body diagonal of the internal matrix using the length, width and height of the workpiece to obtain the outer dimensions of the cube inscribed on the fiber body; defining the spacing between fibers and the spacing between fibers and the boundary of a cube inscribed outside the fiber body, calculating the number of fibers based on 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 based on the integer-level fiber number and the input parameters. The dimensions of the external matrix are calculated using the actual dimensions. The calculation formula for the dimensions is: Where d represents the fiber diameter, Length, Width, and Height are the length, width, and height of the workpiece, respectively, and L1 is the 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, L ' 1 is the actual size of the cube inscribed outside the fiber body, and L2 is the size of the external matrix.

Images