A method for establishing RVE model of composite materials considering fiber misalignment

Abstract

The present invention discloses a method for establishing a composite material RVE model taking into account fiber dislocation, which relates to the technical field of composite materials and three-dimensional modeling. The method generates cross sections of multiple fibers in a two-dimensional RVE region and records the coordinates of the center points of the fibers. According to the center point of each fiber in the two-dimensional RVE region, a corresponding number of center control points are generated for each fiber along the Z-axis direction. The center control points of each fiber are fitted into a curve to obtain a distribution trajectory of the fiber along the Z-axis direction, which is recorded as the center point path of the fiber. According to the cross section of the fiber, a sweep operation is performed along the center point path of each fiber to generate a three-dimensional entity of all fibers, and the entity is assembled with a matrix to form a three-dimensional RVE model. The present invention makes up for the deficiency of existing composite material RVE modeling methods that can only generate straight fibers. The three-dimensional RVE model established by this method is more consistent with the actual composite material structural characteristics.

Description

Technical Field

[0001] The invention relates to the technical field of composite materials and three-dimensional modeling, in particular to a method for establishing a composite material RVE model taking fiber dislocation into consideration. Background Art

[0002] Fiber-reinforced composites have been widely used in hydrogen storage and transportation containers, aerospace, and other fields due to their excellent mechanical properties and lightweight characteristics. Fiber misalignment is a common defect in composite materials. It refers to the phenomenon that fibers randomly deviate from the designed direction. It is usually caused by manufacturing process limitations, material properties, or fluctuations in processing parameters. For example, uneven fiber winding tension, uneven mold surface, or thermal stress deformation caused by the curing process. Fiber misalignment can lead to local stress concentration, reduce the compressive strength and fatigue life of the material, and thus increase the risk of structural failure. Therefore, accurately simulating and analyzing the impact of fiber misalignment on the mechanical properties of composite materials is of great significance for the performance prediction and design application of composite materials.

[0003] Currently, the homogenized prediction and failure analysis of the macroscopic mechanical properties of fiber-reinforced composites, such as the elastic modulus and shear modulus, are mainly achieved based on the Representative Volume Element (RVE) model. This idealized assumption of straight fiber arrangement makes it impossible to consider the complex effects of fiber misalignment in the actual production process, making it difficult to achieve high-precision prediction and analysis of the mechanical properties of composite materials.

[0004] The Representative Volume Element (RVE) model is an important tool for studying the relationship between a material's macroscopic properties and microstructure in fields such as materials science and engineering. An RVE is a representative minimum volume element selected within a material's microstructure that reflects the overall macroscopic properties of the material. Its dimensions typically lie between the characteristic microstructural dimensions (such as fiber diameter and interfiber spacing) and the macrostructural dimensions. Summary of the Invention

[0005] In order to overcome the above-mentioned defects in the prior art, the present invention provides a method for establishing a composite material RVE model taking into account fiber dislocation, which makes up for the deficiency that the existing composite material RVE modeling method can only generate straight fibers. The three-dimensional RVE model established by this method is more consistent with the actual composite material structural characteristics.

[0006] To achieve the above object, the present invention adopts the following technical solutions, including:

[0007] A method for establishing a composite material RVE model considering fiber dislocation includes the following steps:

[0008] S1, generate cross sections of multiple fibers in the two-dimensional RVE region and record the coordinates of the center points of the fibers in the two-dimensional RVE region;

[0009] S2, based on the center point of each fiber in the two-dimensional RVE area obtained in step S1, generate corresponding central control points for each fiber along the Z-axis direction; fit the central control points of each fiber into a curve to obtain the distribution trajectory of the fiber along the Z-axis direction, which is recorded as the center point path of the fiber;

[0010] S3, based on the cross section of the fiber, a sweep operation is performed along the center point path of each fiber to generate a three-dimensional entity of all fibers, and assembled with a matrix to form a three-dimensional RVE model; the matrix is ​​obtained by extending the two-dimensional RVE area along the Z axis direction.

[0011] Preferably, step S1 is specifically as follows:

[0012] S11, define the area of ​​the two-dimensional RVE region S 0.Radius of fiber cross section R and area S , fiber volume percentage V f , the minimum value of fiber spacing L min and maximum value L max ;

[0013] S12, based on fiber volume percentage V f and fiber radius R The theoretical fiber number is calculated as S 0 V f / S ;

[0014] S13, randomly generate the center point coordinates of each fiber in the two-dimensional RVE area, and each generated center point coordinate is required to be within the two-dimensional RVE area and the distance from the center point of the adjacent fiber is 1 / 4. L To meet the fiber spacing requirements, L min ≤L≤L max If the generated center point coordinates do not meet the conditions, the center point coordinates are deleted and new center point coordinates are randomly generated. If the generated center point coordinates meet the conditions, the center point coordinates are retained and the center point coordinates of the next fiber are generated until the number of generated center point coordinates, that is, the number of fibers, reaches the theoretical number of fibers.

[0015] Preferably, step S2 is specifically as follows:

[0016] S21, defines that each fiber has N central control point, and the first i The central control points form a layer and are recorded as i layer, i =1,2,..., N ;

[0017] S22, for the first central control point of each fiber, the center point of each fiber in the two-dimensional RVE area is used as the corresponding first central control point. The coordinates of the first central control point are ( x 1, y 1, 0), where ( x 1, y 1) The coordinates of the center point of the fiber obtained in step S1 within the two-dimensional RVE area;

[0018] S23, for each fiber from the 2nd to the N Central control point, i +1 center control point coordinates ( x i+1 , y i+1 , z i+1 ) Based on the previous i The coordinates of the central control point ( x i , y i , z i ) is generated, and the specific generation algorithm is as follows:

[0019] Randomly generate an angle i , given a maximum moving distance r , calculate the single movement increment Δ of the center control point in the X-axis and Y-axis directions r x and Δ r y :

[0020] Δ r x = r ·cos i , Δ r y = r ·sin i ;

[0021] x i+1 = x i + α · Ninc ·Δ r x , y i+1 = y i + α · N inc ·Δ r y ;

[0022] in, N inc is the number of moves, c is the thickness of the three-dimensional RVE model, i.e., the length along the Z axis. α A random number between 0 and 1;

[0023] After each move, the remaining center control points of the current layer are traversed to determine whether the generated center control point interferes with the remaining center control points, that is, to determine whether the distance between the generated center control point coordinates and the remaining center control point coordinates is greater than the fiber diameter. If the distance is less than or equal to the fiber diameter, that is, interference occurs, the current move is canceled, and the center control point coordinates after the previous move are used as the second move. i +1 center control point coordinates; if the spacing is greater than the fiber diameter, that is, no interference occurs, the center control point coordinates after the current move are retained and the next move is continued until interference occurs;

[0024] S24, fitting a number of center control points of each fiber to generate a Bezier curve, and obtaining a distribution trajectory of the fiber along the Z-axis direction, which is recorded as the center point path of the fiber.

[0025] Preferably, when the fiber is located at the boundary or corner of the three-dimensional RVE model, that is, when the distance from the center control point of the fiber to the surface of the three-dimensional RVE model is less than the fiber radius, the fiber is mirrored to ensure that the generated three-dimensional RVE model meets the geometric periodicity conditions.

[0026] Preferably, in step S24, the Bessel function is:

[0027] ;

[0028] in, P i For the i The coordinates of the central control point, B ( t ) represents the Bezier curve, parameter t ∈[0,1], t The number of represents the number of curve points, is the Bernstein basis function;

[0029] ;

[0030] in, Indicates from N Select from the elements i The number of combinations of elements.

[0031] Preferably, the misalignment angle Γ of each fiber is defined as: P c ( x c , y c , z c )’s maximum angular offset; the specific calculation process is as follows:

[0032] ;

[0033] in, v for point P c The direction vector of x 1, y 1) is the center point coordinate of the fiber in the two-dimensional RVE region obtained in step S1, ( x 1, y 1, 0) is the coordinate of the first central control point;

[0034] ;

[0035] Among them, Г c for point P c The angular offset of the point P c Direction vector v Unit normal vector relative to the XY plane n The angle of (0,0,1);

[0036] .

[0037] Preferably, the method further includes step S4 of assigning material properties to the matrix and the fibers and dividing the meshes;

[0038] Among them, the specific ways to impart material properties to fibers are:

[0039] Any point on the path of the fiber center point P c Establish a local coordinate system at point P c is the origin of the local coordinate system, point P cThe tangent direction at is the z-axis direction of the local coordinate system, and any two mutually perpendicular directions in the fiber cross section are set as the x-axis direction and y-axis direction of the local coordinate system, thereby completing the assignment of transversely anisotropic material properties to the fiber.

[0040] The present invention also provides a computer program product, which includes a computer program / instruction, and when the computer program / instruction is executed by a processor, it implements the method for establishing a composite material RVE model considering fiber dislocation.

[0041] The advantages of the present invention are:

[0042] (1) The present invention provides a method for establishing a three-dimensional RVE model of composite materials taking into account fiber dislocation, which can efficiently and reasonably characterize the random fiber dislocation characteristics and ensure the authenticity of the model's geometric configuration. The three-dimensional RVE model established by this method can be more accurately used for simulation-based composite material performance prediction and failure analysis, providing a scientific basis for the design and application of composite material related research.

[0043] (2) Based on the random expansion algorithm, random movement of the central control point and Bezier curve fitting, the present invention realizes the efficient modeling of random fiber dislocation in the three-dimensional RVE model of the composite material, and quantifies the degree of fiber dislocation, which more realistically displays the fiber distribution of the composite material.

[0044] (3) The present invention efficiently completes the geometric generation, material property assignment, and meshing of a three-dimensional RVE model containing fiber misalignment through parameterization. The obtained models with different fiber volume fractions can be widely used in various RVE simulation problems and used to explore important issues such as the quantitative effect of fiber misalignment on the mechanical properties of composite materials. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] Figure 1 The present invention is a flow chart of a method for establishing an RVE model for composite materials taking fiber dislocation into consideration.

[0046] Figure 2 Schematic diagram of the modeling process of the three-dimensional RVE model of the present invention. DETAILED DESCRIPTION

[0047] 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.

[0048] Depend on Figure 1As shown, the present invention provides a method for establishing a composite material RVE model considering fiber misalignment, including the following process:

[0049] S1. Generate cross sections of multiple fibers within a two-dimensional RVE region using a random sequential expansion algorithm, and record the coordinates of the center points of the fibers within the two-dimensional RVE region; the two-dimensional RVE region is the region of the RVE (representative volume element) within the XY plane (the two-dimensional plane formed by the X-axis and the Y-axis).

[0050] Step S1 is specifically as follows:

[0051] S11, determine the basic geometric parameters of the model, including the size of the 2D RVE area a × b (width × height), radius of fiber cross section R and area S , fiber volume percentage V f , the minimum value of fiber spacing L min and maximum value L max In this embodiment, the fiber is a circular fiber, and the area of ​​the fiber cross section is S=πR 2 .

[0052] S12, based on fiber volume percentage V f The theoretical fiber number is calculated as abV f / πR 2 ; Among them, fiber volume percentage V f It refers to the ratio of the volume occupied by fibers to the volume of the entire composite material in a fiber-reinforced composite material, usually expressed as a percentage.

[0053] S13, randomly generate the center point coordinates of each fiber in the two-dimensional RVE area, and each generated center point coordinate is required to be within the two-dimensional RVE area and the distance from the center point of the adjacent fiber is 1 / 4. L To meet the fiber spacing requirements, L min ≤L≤L max If the generated center point coordinates do not meet the conditions, the center point coordinates are deleted and new center point coordinates are randomly generated. If the generated center point coordinates meet the conditions, the center point coordinates are retained and the center point coordinates of the next fiber are generated. The process is repeated until the number of generated center point coordinates, that is, the number of fibers, reaches the theoretical number of fibers. The center point coordinates of all randomly generated fibers in the two-dimensional RVE area are saved and each fiber is numbered.

[0054] S2, according to the center point of each fiber in the two-dimensional RVE area obtained in step S1, generate corresponding several center control points for each fiber along the Z-axis direction; based on the Bezier curve, fit the several center control points of each fiber to obtain the distribution trajectory of each fiber along the Z-axis direction, which is recorded as the center point path of the fiber.

[0055] Step S2 is specifically as follows:

[0056] S21, defines that each fiber has N central control point, and the first i The central control points form a layer and are recorded as i layer, i =1,2,..., N .

[0057] S22, for the first central control point of each fiber, the center point of each fiber in the two-dimensional RVE area is used as the corresponding first central control point. The coordinates of the first central control point are ( x 1, y 1, 0), where ( x 1, y 1) is the coordinate of the center point of the fiber in the two-dimensional RVE area obtained in step S1.

[0058] S23, for each fiber from the 2nd to the N Central control point, i +1 center control point coordinates ( x i+1 , y i+1 , z i+1 ) Based on the previous i The coordinates of the central control point ( x i , y i , z i ) is generated, and the specific generation algorithm is as follows:

[0059] Randomly generate an angle i ∈[0,2π), given a maximum moving distance r , calculate the single movement increment Δ of the center control point in the X-axis and Y-axis directions r x and Δ r y :

[0060] Δ r x= r ·cos i , Δ r y = r ·sin i ;

[0061] x i+1 = x i + α · N inc ·Δ r x , y i+1 = y i + α · N inc ·Δ r y ;

[0062] in, N inc The number of moves is used to gradually increase the distance the center control point moves in the XY plane; c is the thickness of the three-dimensional RVE model, i.e., the length along the Z axis. α A random number between 0 and 1; α It is a random number between 0 and 1, so that the moving path of the central control point is different each time, ensuring the randomness of the dislocation generation process.

[0063] After each move, the remaining center control points of the current layer are traversed to determine whether the generated center control point interferes with the remaining center control points, that is, to determine whether the distance between the generated center control point coordinates and the remaining center control point coordinates is greater than the fiber diameter. If the distance is less than or equal to the fiber diameter, that is, interference occurs, the current move is canceled, and the center control point coordinates after the previous move are used as the second move. i +1 center control point coordinates; if the spacing is greater than the fiber diameter, that is, no interference occurs, the center control point coordinates after the current movement are retained, and the next movement is continued based on this until interference occurs.

[0064] When a fiber is located at a boundary or corner of a 3D RVE model—that is, when the distance from the fiber's central control point to the four surfaces of the 3D RVE model (top, bottom, front, and back) is less than the fiber radius—the fiber needs to be geometrically periodicized to ensure that the generated RVE model satisfies the geometric periodicity condition as much as possible. Specifically, a mirror image fiber can be generated at the opposite or diagonal corner of the 3D RVE model. For example, if a fiber is located at the top boundary, a mirror image fiber is generated at the bottom boundary (in the –Y direction). Specifically, the fiber is copied and translated by –b units in the Y direction. If a fiber is located at the top left boundary, a mirror image fiber is generated at the bottom right boundary (in the +X and –Y directions). Specifically, the fiber is copied and translated by a and –b units in the X and Y directions, respectively. After generating the mirror image fiber, it is still necessary to check whether it interferes with other fibers.

[0065] S24, fitting a number of center control points of each fiber to generate a Bezier curve, and obtaining a distribution trajectory of the fiber along the Z-axis direction, which is recorded as the center point path of the fiber.

[0066] Among them, the Bezier curve function defines the parametric expression of the curve through a set of central control points. N The central control point is expressed as follows:

[0067] ;

[0068] in, P i For the i The coordinates of the central control point, B ( t ) represents a Bezier curve, is the Bernstein basis function;

[0069] ;

[0070] in, Indicates from N Select from the elements i The number of combinations of elements; parameter t ∈[0,1], t The number of represents the number of curve points, by changing t The number of can control the density of generated curve points, in order to ensure that each fiber N The smoothness of the curve generated by the central control point is t The number of is 100, that is, 100 values ​​evenly distributed from 0 to 1, which ensures the smoothness and accuracy of the curve.

[0071] S3. Based on the cross section of the fiber, a sweep operation is performed along the center point path of each fiber to generate a three-dimensional entity of all fibers, and assembled with a matrix to form a three-dimensional RVE model; the matrix is ​​obtained by extending the two-dimensional RVE area along the Z-axis direction.

[0072] In step S3, a circular cross section and interface of the fiber is created, and a sweep operation is performed along the center point path of the fiber, ensuring that the internal boundary and the normal direction of the cross section are maintained during the sweep process to obtain a three-dimensional entity of each fiber. a × b × c The rectangular parallelepiped is used as the matrix, and the fiber boundary is retained to merge the fiber and matrix into a complete three-dimensional RVE model.

[0073] The misalignment angle Г of each fiber is defined as: any point on the path of the fiber center point P c ( x c , y c , z c ) direction vector v Unit normal vector relative to the XY two-dimensional plane n The maximum value of the angle between (0,0,1). The specific calculation process is as follows:

[0074] Direction vector v The calculation formula is:

[0075] ;

[0076] in, v for point P c The direction vector of x 1, y 1) is the center point coordinate of the fiber in the two-dimensional RVE region obtained in step S1, ( x 1, y 1, 0) is the coordinate of the first central control point;

[0077] Any point P c The angular offset Г c Obtained by the angle calculation formula between two vectors:

[0078] ;

[0079] The misalignment angle Г of the fiber is obtained as:

[0080] .

[0081] The size of the fiber misalignment angle Г is mainly affected by factors such as the initial distribution of fibers on the XY plane, the random movement process of fibers, and the fiber spacing. It can be achieved by reducing the fiber volume fraction. V f , increase fiber spacing L Increase the fiber misalignment angle Г by other means, and improve the fiber Z The degree of deviation from the axis.

[0082] S4. Assign material properties to the matrix and fiber respectively and divide the mesh to complete the RVE modeling.

[0083] Determine the material properties of the matrix, including mechanical properties such as elastic modulus E , Poisson's ratio n and assign it to the base part.

[0084] Due to the transverse anisotropy of the fiber, it is necessary to determine the mechanical properties of the fiber, such as the elastic modulus ( E 11 、 E 22 、 E 33 ), shear modulus ( G 12 、 G 13 、 G 23 ), Poisson's ratio ( n 12 ,n 13 ,n 23 ) in all directions, so a local coordinate system is defined so that the fiber material direction is consistent with the geometric direction, and the material properties are assigned to the fiber part accordingly.

[0085] Regarding the establishment of the fiber local coordinate system, in order to correctly set the material direction of the fiber, first select the geometric area of ​​a single fiber. The main direction of the fiber is the fiber axis, which can be defined based on the center point path of the fiber. At any point on the center point path of the fiber P c Establish a local coordinate system at point P c is the origin of the local coordinate system, point P c The tangent direction at is the z-axis direction of the local coordinate system, and any two mutually perpendicular directions in the fiber cross section are set as the x-axis direction and y-axis direction of the local coordinate system, thereby completing the assignment of transversely anisotropic material properties to the fiber.

[0086] like Figure 2As shown, when the size of the 3D RVE model is 65 μm × 65 μm × 65 μm, the fiber volume percentage V f 50%, maximum moving distance r =0.01μm, radius of the fiber cross section R =4μm, fiber spacing 0.4μm≤ L ≤0.5μm, number of center control points N = 20. The meshing parameters are defined, with a mesh size of 2 μm and a 10-node quadratic tetrahedron element type (C3D10). The 3D RVE model is divided into 724,916 elements. Considering that the 3D RVE model has small local features, such as at boundaries, this meshing approach ensures both computational accuracy and efficiency.

[0087] The above are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for establishing a composite material RVE model considering fiber dislocation, characterized in that: The following steps are involved: S1, generate cross sections of multiple fibers in the two-dimensional RVE region and record the coordinates of the center points of the fibers in the two-dimensional RVE region; S2, based on the center point of each fiber in the two-dimensional RVE area obtained in step S1, generate corresponding central control points for each fiber along the Z-axis direction; fit the central control points of each fiber into a curve to obtain the distribution trajectory of the fiber along the Z-axis direction, which is recorded as the center point path of the fiber; S3, performing a sweep operation along the center point path of each fiber according to the cross section of the fiber to generate a three-dimensional solid of all fibers, and assembling them with a matrix to form a three-dimensional RVE model; the matrix is ​​obtained by extending the two-dimensional RVE area along the Z axis; Step S2 is specifically as follows: S21, defines that each fiber has N central control point, and the first i The central control points form a layer and are recorded as i layer, i =1,2,..., N ; S22, for the first central control point of each fiber, the center point of each fiber in the two-dimensional RVE area is used as the corresponding first central control point. The coordinates of the first central control point are ( x 1, y 1, 0), where ( x 1, y 1) The coordinates of the center point of the fiber obtained in step S1 within the two-dimensional RVE area; S23, for each fiber from the 2nd to the N Central control point, i +1 center control point coordinates ( x i+1 , y i+1 , z i+1 ) Based on the previous i The coordinates of the central control point ( x i , y i , z i ) is generated, and the specific generation algorithm is as follows: Randomly generate an angle θ , given a maximum moving distance ρ , calculate the single movement increment Δ of the center control point in the X-axis and Y-axis directions ρ x and Δ ρ y : D ρ x = ρ ·cos θ ,D ρ y = ρ ·sin θ ; x i+1 = x i + α · N inc ·D ρ x , y i+1 = y i + α · N inc ·D ρ y ; in, N inc is the number of moves, c is the thickness of the three-dimensional RVE model, i.e., the length along the Z axis. α A random number between 0 and 1; After each move, the remaining center control points of the current layer are traversed to determine whether the generated center control point interferes with the remaining center control points, that is, to determine whether the distance between the generated center control point coordinates and the remaining center control point coordinates is greater than the fiber diameter. If the distance is less than or equal to the fiber diameter, that is, interference occurs, the current move is canceled, and the center control point coordinates after the previous move are used as the second move. i +1 center control point coordinates; if the spacing is greater than the fiber diameter, that is, no interference occurs, the center control point coordinates after the current move are retained and the next move is continued until interference occurs; S24, fitting a number of center control points of each fiber to generate a Bezier curve, and obtaining a distribution trajectory of the fiber along the Z-axis direction, which is recorded as the center point path of the fiber.

2. The method for establishing a composite material RVE model considering fiber dislocation according to claim 1, characterized in that: Step S1 is specifically as follows: S11, define the area of ​​the two-dimensional RVE region S 0.Radius of fiber cross section R and area S , fiber volume percentage V f , the minimum value of fiber spacing L min and maximum value L max ; S12, based on fiber volume percentage V f and fiber radius R The theoretical fiber number is calculated as S 0 V f / S ; S13, randomly generate the center point coordinates of each fiber in the two-dimensional RVE area, and each generated center point coordinate is required to be within the two-dimensional RVE area and the distance from the center point of the adjacent fiber is 1 / 4. L To meet the fiber spacing requirements, L min ≤L≤ L max If the generated center point coordinates do not meet the conditions, the center point coordinates are deleted and new center point coordinates are randomly generated. If the generated center point coordinates meet the conditions, the center point coordinates are retained and the center point coordinates of the next fiber are generated until the number of generated center point coordinates, that is, the number of fibers, reaches the theoretical number of fibers.

3. The method for establishing a composite material RVE model considering fiber dislocation according to claim 1, characterized in that: When the fiber is located at the boundary or corner of the three-dimensional RVE model, that is, when the distance from the central control point of the fiber to the surface of the three-dimensional RVE model is less than the fiber radius, the fiber is mirrored to ensure that the generated three-dimensional RVE model meets the geometric periodicity condition.

4. The method for establishing a composite material RVE model considering fiber dislocation according to claim 1, characterized in that: In step S24, the Bessel function is: ； in, P i For the i The coordinates of the central control point, B ( t ) represents the Bezier curve, parameter t ∈[0,1], t The number of represents the number of curve points, is the Bernstein basis function; ； in, Indicates from N Select from the elements i The number of combinations of elements.

5. The method for establishing a composite material RVE model considering fiber dislocation according to claim 1, characterized in that: The misalignment angle Г of each fiber is defined as: any point on the path of the fiber center point P c ( x c , y c , z c )’s maximum angular offset; the specific calculation process is as follows: ； in, v for point P c The direction vector of x 1, y 1) is the center point coordinate of the fiber in the two-dimensional RVE region obtained in step S1, ( x 1, y 1, 0) is the coordinate of the first central control point; ； Among them, Г c for point P c The angular offset of the point P c Direction vector v Unit normal vector relative to the XY plane n The angle of (0,0,1); 。 6. The method for establishing a composite material RVE model considering fiber dislocation according to claim 1, characterized in that: The method further includes step S4 of assigning material properties to the matrix and the fibers and dividing the meshes; Among them, the specific ways to impart material properties to fibers are: Any point on the path of the fiber center point P c Establish a local coordinate system at point P c is the origin of the local coordinate system, point P c The tangent direction at is the z-axis direction of the local coordinate system, and any two mutually perpendicular directions in the fiber cross section are set as the x-axis direction and y-axis direction of the local coordinate system, thereby completing the assignment of transversely anisotropic material properties to the fiber.

7. A computer program product, characterized in that The invention comprises a computer program / instruction, which, when executed by a processor, realizes a method for establishing a composite material RVE model taking fiber dislocation into consideration as described in any one of claims 1 to 6.

Images