Cell modeling method, system, equipment and medium for two-dimensional braided composite materials
By obtaining the microscopic cell parameters of the woven composite material and generating the microscopic cell model, the problem of jagged discontinuity in the fiber bundle/matrix interface in the existing technology is solved, and the accuracy of the simulation analysis is significantly improved.
Patent Information
- Application Number
- CN202411287706.0
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2024-09-14
- Publication Date
- 2025-01-21
- Estimated Expiration
- 2044-09-14
AI Technical Summary
In the prior art, the mesoscopic cell element modeling of braided composite materials uses rectangular units, resulting in serrated discontinuous features appearing at the fiber bundle/matrix interface, affecting the accuracy of the numerical simulation analysis results.
By obtaining the cellular cells of the woven composite material and determining its size parameters, warp parameters and weft parameters, these parameters are input into the finite element software Abaqus to generate a cellular cell model. The model includes creating a three-dimensional solid part, determining fiber bundle outlines and sweeping curves, and sweeping along these curves to divide into warp fiber bundles, weft fiber bundles and matrix, and finally generate a microscopic cell model of the woven composite material.
This method does not require model import through pre-processing software, the geometric model has high accuracy, and the fiber bundle/matrix interface does not have jagged discontinuous features, which significantly improves the accuracy and precision of the calculation results.
Smart Images

Figure CN118798010B_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of three-dimensional modeling, and in particular to a cell modeling method, system, equipment and medium for two-dimensional braided composite materials. Background Art
[0002] Braided composite materials have excellent mechanical properties such as high specific strength, high specific stiffness, and low density, and have broad application prospects in engineering fields such as aerospace, vehicles, and construction. Composite materials are highly designable at multiple levels of materials, processes, and structures. Changes in their microscopic braiding parameters have a significant impact on the mechanical behavior of macroscopic materials / components. Therefore, mesoscopic cell modeling and analysis is a key means to fully explore the mechanical properties of braided composite materials and design composite components.
[0003] In the prior art, the micro-cell modeling of woven composite materials mostly adopts pre-processing modeling software, and then imports finite element software for numerical simulation analysis. The specific process includes: first, use the import function of the software to open the dialog box and import the geometric model of the woven composite material. For woven composite materials, it is necessary to understand its weaving process in detail, especially the movement mode of the yarn carrier and the spatial curve of the yarn. Based on the variable cross-section theory of the yarn and the analysis of the weaving process, the parametric modeling function in the pre-processing software is used to establish the spatial entity model of the yarn. The yarn model is combined with the matrix model to create a three-dimensional woven composite material entity model containing the micro-structure characteristics of the yarn. The composite material model is divided into a series of small units, and corresponding material properties are assigned to each unit. Boundary conditions are assigned to the model to obtain a micro-cell model. After the modeling is completed, the model is checked and optimized. Finally, the built model is exported to a file format that can be recognized by the finite element analysis software for subsequent finite element analysis.
[0004] The defects of the above-mentioned prior art are: the pre-processing modeling software adopts rectangular parallelepiped units, which show obvious jagged discontinuity characteristics at the fiber bundle / matrix interface, affecting the accuracy of the numerical simulation analysis results. Summary of the invention
[0005] Based on this, it is necessary to provide a cellular modeling method, system, equipment and medium for two-dimensional woven composite materials to address the above-mentioned technical problems.
[0006] An embodiment of the present invention provides a cell modeling method for a two-dimensional braided composite material, comprising:
[0007] Obtaining the mesoscopic cells of the woven composite material, and determining the size parameters, warp yarn parameters and weft yarn parameters of the mesoscopic cells; wherein the size parameters include: cell length Length, cell width Width and cell thickness Thickness;
[0008] The size parameters, warp yarn parameters and weft yarn parameters of the meso-cell are input into the finite element software Abaqus to generate a meso-cell model. Generating the meso-cell model includes the following steps:
[0009] Create a three-dimensional solid part according to the size parameters of the mesoscopic cell; the length of the three-dimensional solid part is 2×Length, the width is 2×Width, and the thickness is 2×Thickness;
[0010] Taking the YZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, the contour of the warp fiber bundle is determined according to the size parameters and the warp parameters; taking the XZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, the contour of the weft fiber bundle is determined according to the size parameters and the weft parameters; taking the x=2×Length plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, the weft sweep curve is determined according to the size parameters and the weft parameters; taking the y=2×Width plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, the warp sweep curve is determined according to the size parameters and the warp parameters;
[0011] Sweep the warp fiber bundle contour and the weft fiber bundle contour along the warp sweep curve and the weft sweep curve respectively, and divide the three-dimensional solid part into the warp fiber bundle, the weft fiber bundle and the matrix;
[0012] The warp fiber bundles, weft fiber bundles and the parts of the matrix with 0≤x≤Length, 0≤y≤Width, and 0≤z≤Thickness are retained to obtain a mesoscopic cell model of the woven composite material; wherein x is the length of the mesoscopic cell model, y is the width of the mesoscopic cell model, and z is the thickness of the mesoscopic cell model.
[0013] In addition, the warp parameters of the meso-cell include: warp major axis WarpL, warp minor axis WarpS, and warp curvature radius WarpR;
[0014] The weft yarn parameters of the micro-cell include: weft yarn major axis WeftL, weft yarn minor axis WeftS, and weft yarn curvature radius WeftR.
[0015] In addition, the YZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus is used as a sketch, and the warp fiber bundle profile is determined according to the size parameters and the warp yarn parameters, which specifically includes:
[0016] The side point (0, Width / 2, Thickness / 2) is used as the origin of the sketch coordinates, (0, -Thickness / 4) is used as the center coordinates of the ellipse contour, (WarpL / 2, -Thickness / 4) is used as the endpoint of the major axis of the ellipse contour, and (0, WarpS / 2-Thickness / 4) is used as the endpoint of the minor axis of the ellipse contour to obtain the contour of the first warp fiber bundle;
[0017] The side point (0, Width / 2, Thickness / 2) is used as the origin of the sketch coordinates, (-Width / 2, Thickness / 4) is used as the center coordinates of the ellipse contour, (WarpL / 2-Width / 2, Thickness / 4) is used as the endpoint of the major axis of the ellipse contour, and (-Width / 2, WarpS / 2+Thickness / 4) is used as the endpoint of the minor axis of the ellipse contour, to obtain the contour of the second warp fiber bundle;
[0018] Take the side point (0, Width / 2, Thickness / 2) as the sketch coordinate origin, (Width / 2, Thickness / 4) as the center coordinates of the ellipse contour, (WarpL / 2+Width / 2, Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (Width / 2, WarpS / 2+Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the third warp fiber bundle.
[0019] In addition, the XZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus is used as a sketch, and the weft fiber bundle profile is determined according to the size parameters and the weft yarn parameters, which specifically includes:
[0020] The side point (Length / 2, 0, Thickness / 2) is used as the origin of the sketch coordinates, (0, Thickness / 4) is used as the center coordinates of the ellipse contour, (WeftL / 2, Thickness / 4) is used as the endpoint of the major axis of the ellipse contour, and (0, WeftS / 2+Thickness / 4) is used as the endpoint of the minor axis of the ellipse contour to obtain the contour of the first weft fiber bundle;
[0021] Take the side point (Length / 2, 0, Thickness / 2) as the origin of the sketch coordinates, (-Length / 2, -Thickness / 4) as the center coordinates of the ellipse contour, (WeftL / 2-Length / 2, -Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (-Length / 2, WeftS / 2-Thickness / 4) as the endpoint of the minor axis of the ellipse contour, and obtain the contour of the second weft fiber bundle;
[0022] Take the side point (Length / 2, 0, Thickness / 2) as the sketch coordinate origin, (Length / 2, -Thickness / 4) as the center coordinate of the ellipse contour, (WeftL / 2+Length / 2, -Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (Length / 2, WeftS / 2-Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the third weft fiber bundle.
[0023] In addition, the part of the x=2×Length plane of the three-dimensional coordinate system of the finite element software Abaqus is used as a sketch, and the weft yarn sweep curve is determined according to the size parameters and the weft yarn parameters, which specifically includes:
[0024] Take the side point (2×Length, Width / 2, Thickness / 2) as the origin of the sketch coordinates, and draw two circles with (-Width / 2, WeftR-Thickness / 4) and (0, Thickness / 4- WeftR) as the centers and WeftR as the radius.
[0025] Draw a common tangent line between the two circles, cut the redundant curve, and mirror-copy the remaining curve with the Y axis of the three-dimensional coordinate system of the finite element software Abaqus as the symmetry axis to obtain the first weft sweep curve;
[0026] Then, the X-axis of the three-dimensional coordinate system of the finite element software Abaqus is used as the symmetry axis, and the first weft yarn sweep curve is mirror-copied to obtain the first weft yarn sweep curve and the second weft yarn sweep curve.
[0027] In addition, the part of the y=2×cell width Width plane of the three-dimensional coordinate system of the finite element software Abaqus is used as a sketch, and the warp sweep curve is determined according to the size parameters and the warp parameters, which specifically includes:
[0028] Take the side point (Length / 2, 2×Width, Thickness / 2) as the origin of the sketch coordinate system, and draw two circles with (-Length / 2, WarpR-Thickness / 4) and (0, Thickness / 4-WarpR) as the centers and WarpR as the radius.
[0029] Draw a common tangent line between the two circles, cut the redundant curve, and mirror-copy the remaining curve with the Y axis of the three-dimensional coordinate system of the finite element software Abaqus as the symmetry axis to obtain the first warp sweep curve;
[0030] The X-axis of the three-dimensional coordinate system of the finite element software Abaqus is taken as the symmetry axis, and the first warp sweep curve is mirror-copied to obtain the first warp sweep curve and the second warp sweep curve.
[0031] In addition, the cell modeling system for two-dimensional braided composite materials includes:
[0032] A cell acquisition module is used to acquire the microscopic cells of the woven composite material and determine the size parameters, warp parameters and weft parameters of the microscopic cells; wherein the size parameters include: cell length Length, cell width Width, cell thickness Thickness;
[0033] The model building module is used to input the size parameters, warp yarn parameters and weft yarn parameters of the meso-cell into the finite element software Abaqus to generate a meso-cell model; the model building module includes:
[0034] The entity creation unit is used to create a three-dimensional entity part according to the size parameters of the mesoscopic cell; the length of the three-dimensional entity part is 2×Length, the width is 2×Width, and the thickness is 2×Thickness;
[0035] The curve construction unit is used to use the portion of the YZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch to determine the contour of the warp fiber bundle according to the size parameters and the warp parameters; use the portion of the XZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch to determine the contour of the weft fiber bundle according to the size parameters and the weft parameters; use the portion of the x=2×Length plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch to determine the weft sweep curve according to the size parameters and the weft parameters; use the portion of the y=2×Width plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch to determine the warp sweep curve according to the size parameters and the warp parameters;
[0036] A sweeping unit is used to sweep the warp fiber bundle contour and the weft fiber bundle contour along the warp sweeping curve and the weft sweeping curve respectively, so as to divide the three-dimensional solid part into the warp fiber bundle, the weft fiber bundle and the matrix;
[0037] The model building unit is used to retain the parts of the warp fiber bundles, the weft fiber bundles and the matrix with 0≤x≤Length, 0≤y≤Width, and 0≤z≤Thickness, so as to obtain a mesoscopic cell model of the woven composite material; wherein x is the length of the mesoscopic cell model, y is the width of the mesoscopic cell model, and z is the thickness of the mesoscopic cell model.
[0038] In addition, a computer device includes a memory and a processor, wherein a computer program is stored in the memory, and when the processor executes the computer program, the steps of the above-mentioned two-dimensional woven composite material cell modeling method are implemented.
[0039] In addition, a storage medium stores a computer program, which implements the steps of the above-mentioned two-dimensional woven composite material cell modeling method when executed by a processor.
[0040] Compared with the prior art, the cell modeling method, system, device and medium of the two-dimensional braided composite material provided by the embodiments of the present invention have the following beneficial effects:
[0041] The meso-cells of the woven composite material are obtained, and the size parameters, warp yarn parameters and weft yarn parameters of the meso-cells are determined; wherein the size parameters include: cell length Length, cell width Width and cell thickness Thickness; the size parameters, warp yarn parameters and weft yarn parameters of the meso-cells are input into the finite element software Abaqus to generate a meso-cell model.
[0042] The present invention does not need to import the model through pre-processing software, and the geometric model has high accuracy. In the process of sweeping the warp fiber bundle contour and the weft fiber bundle contour along the warp yarn sweep curve and the weft yarn sweep curve respectively, no jagged discontinuity features will appear at the fiber bundle / matrix interface, which can significantly improve the accuracy and precision of the calculation results. BRIEF DESCRIPTION OF THE DRAWINGS
[0043] Figure 1 A flowchart of a method for cellular modeling of a two-dimensional braided composite material provided in one embodiment. DETAILED DESCRIPTION
[0044] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention.
[0045] In one embodiment, a cell modeling method for a two-dimensional braided composite material is provided, such as Figure 1 As shown, the method includes:
[0046] Step 1 Obtain the mesocell of the woven composite material and determine the size parameters, warp parameters, and weft parameters of the mesocell. The size parameters of the mesocell include: cell length Length, cell width Width, and cell thickness Thickness. The warp parameters of the mesocell include: warp major axis WarpL, warp minor axis WarpS, and warp curvature radius WarpR. The weft parameters of the mesocell include: weft major axis WeftL, weft minor axis WeftS, and weft curvature radius WeftR.
[0047] Step 2 Input the size parameters, warp yarn parameters and weft yarn parameters of the meso-cell into the finite element software Abaqus to generate a meso-cell model. Generating a meso-cell model includes the following steps:
[0048] Step 2.1 Create a three-dimensional solid part according to the size parameters of the mesoscopic cell; the length of the three-dimensional solid part is 2×Length, the width is 2×Width, and the thickness is 2×Thickness.
[0049] Step 2.2 According to the size parameters, warp yarn parameters and weft yarn parameters, the warp yarn fiber bundle profile, weft yarn fiber bundle profile, weft yarn sweep curve and warp yarn sweep curve are respectively made. The specific process is as follows:
[0050] (1) The YZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus is used as a sketch, and the warp fiber bundle profile is determined according to the size parameters and warp yarn parameters, which specifically includes:
[0051] Take the side point (0, Width / 2, Thickness / 2) as the sketch coordinate origin, (0, -Thickness / 4) as the center coordinate of the ellipse contour, (WarpL / 2, -Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (0, WarpS / 2-Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the first warp fiber bundle.
[0052] Take the side point (0, Width / 2, Thickness / 2) as the sketch coordinate origin, (-Width / 2, Thickness / 4) as the center coordinates of the ellipse contour, (WarpL / 2-Width / 2, Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (-Width / 2, WarpS / 2+Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the second warp fiber bundle.
[0053] Take the side point (0, Width / 2, Thickness / 2) as the sketch coordinate origin, (Width / 2, Thickness / 4) as the center coordinates of the ellipse contour, (WarpL / 2+Width / 2, Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (Width / 2, WarpS / 2+Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the third warp fiber bundle.
[0054] (2) Take the XZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, and determine the weft fiber bundle contour based on the size parameters and weft yarn parameters:
[0055] Take the side point (Length / 2, 0, Thickness / 2) as the sketch coordinate origin, (0, Thickness / 4) as the center coordinate of the ellipse contour, (WeftL / 2, Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (0, WeftS / 2+Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the first weft fiber bundle.
[0056] Take the side point (Length / 2, 0, Thickness / 2) as the sketch coordinate origin, (-Length / 2, -Thickness / 4) as the center coordinates of the ellipse contour, (WeftL / 2-Length / 2, -Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (-Length / 2, WeftS / 2-Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the second weft fiber bundle.
[0057] Take the side point (Length / 2, 0, Thickness / 2) as the sketch coordinate origin, (Length / 2, -Thickness / 4) as the center coordinate of the ellipse contour, (WeftL / 2+Length / 2, -Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (Length / 2, WeftS / 2-Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the third weft fiber bundle.
[0058] (3) Take the part of the three-dimensional solid part in the x=2×cell length Length plane of the three-dimensional coordinate system of the finite element software Abaqus as the sketch, and determine the weft sweep curve according to the size parameters and weft parameters:
[0059] Take the side point (2×Length, Width / 2, Thickness / 2) as the origin of the sketch coordinates, and draw two circles with (-Width / 2, WeftR-Thickness / 4) and (0, Thickness / 4- WeftR) as the centers and WeftR as the radius.
[0060] Draw a common tangent line between the two circles, cut the redundant curve, and mirror-copy the remaining curve with the Y axis of the three-dimensional coordinate system of the finite element software Abaqus as the symmetry axis to obtain the first weft sweep curve;
[0061] Then, the X-axis of the three-dimensional coordinate system of the finite element software Abaqus is used as the symmetry axis, and the first weft yarn sweep curve is mirror-copied to obtain the first weft yarn sweep curve and the second weft yarn sweep curve.
[0062] (4) Take the part of the y=2×Width plane of the three-dimensional coordinate system of the finite element software Abaqus as the sketch, and determine the warp sweep curve according to the size parameters and warp parameters:
[0063] Take the side point (Length / 2, 2×Width, Thickness / 2) as the sketch coordinate origin, draw two circles with (-Length / 2, WarpR-Thickness / 4) and (0, Thickness / 4-WarpR) as the center and WarpR as the radius.
[0064] Draw the common tangent of the two circles, cut the redundant curve, mirror-copy the remaining curve with the Y axis of the three-dimensional coordinate system of the finite element software Abaqus as the symmetry axis, and obtain the first warp sweep curve.
[0065] The X-axis of the three-dimensional coordinate system of the finite element software Abaqus is used as the symmetry axis, and the first warp sweep curve is mirror-copied to obtain the first warp sweep curve and the second warp sweep curve.
[0066] Step 2.3 Sweep the warp fiber bundle contour and the weft fiber bundle contour along the warp sweep curve and the weft sweep curve respectively, and divide the three-dimensional solid part into the warp fiber bundle, the weft fiber bundle and the matrix.
[0067] Step 2.4: Retain the warp fiber bundles, weft fiber bundles and the part of the matrix with 0≤x≤Length, 0≤y≤Width, 0≤z≤Thickness to obtain the mesoscopic cell model of the woven composite material, where x is the length of the mesoscopic cell model, y is the width of the mesoscopic cell model, and z is the thickness of the mesoscopic cell model.
[0068] Based on the same inventive concept, the present invention also includes:
[0069] Extract the Python command stream generated by the above modeling process, and represent the specific size parameters in the command stream with key size parameters. Encapsulate the modified command stream, generate an Abaqus plug-in, and use the key size parameters of the mesoscopic cell as input parameters in the plug-in interface.
[0070] After completing the above operations for the first time, the developed Abaqus plug-in can be used to input key size parameters to automatically generate a mesoscopic cellular model of a two-dimensional woven composite material with one click.
[0071] In the Python command flow of Abaqus, some operations use masks to represent the selected objects by default. To facilitate indexing and selecting objects during secondary development, you can enter the following command statement in the command bar: session.journalOptions.setValues(replayGeometry=INDEX,recoverGeometry=INDEX). At this time, you can use the number of the collection to be selected in the geometry to index.
[0072] In one embodiment, a cellular modeling system for a two-dimensional braided composite material is provided, the system comprising:
[0073] The cell acquisition module is used to obtain the microscopic cells of the woven composite material and determine the size parameters, warp parameters and weft parameters of the microscopic cells; wherein the size parameters include: cell length Length, cell width Width, cell thickness Thickness.
[0074] The model building module is used to input the size parameters, warp yarn parameters and weft yarn parameters of the meso-cell into the finite element software Abaqus to generate a meso-cell model; the model building module includes:
[0075] The entity creation unit is used to create a three-dimensional entity part according to the size parameters of the mesoscopic cell; the length of the three-dimensional entity part is 2×Length, the width is 2×Width, and the thickness is 2×Thickness.
[0076] The curve construction unit is used for determining the contour of the warp fiber bundle according to the size parameters and the warp parameters by taking the portion of the YZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch; determining the contour of the weft fiber bundle according to the size parameters and the weft parameters by taking the portion of the XZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch; determining the weft sweep curve according to the size parameters and the weft parameters by taking the portion of the x=2×Length plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch; determining the warp sweep curve according to the size parameters and the warp parameters by taking the portion of the y=2×Width plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch.
[0077] The sweeping unit is used to sweep the warp fiber bundle contour and the weft fiber bundle contour along the warp sweeping curve and the weft sweeping curve respectively, so as to divide the three-dimensional solid part into the warp fiber bundle, the weft fiber bundle and the matrix.
[0078] The model building unit is used to retain the warp fiber bundle, the weft fiber bundle and the part of 0≤x≤Length, 0≤y≤Width, 0≤z≤Thickness in the matrix to obtain the mesoscopic cell model of the woven composite material, where x is the length of the mesoscopic cell model, y is the width of the mesoscopic cell model, and z is the thickness of the mesoscopic cell model.
[0079] In one embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and when the processor executes the computer program, steps of a method for cell modeling of a two-dimensional woven composite material are implemented.
[0080] In one embodiment, a storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the steps of a method for cell modeling of a two-dimensional woven composite material are implemented.
[0081] The above-mentioned embodiments only express several implementation methods of the present invention, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for ordinary technicians in this field, several variations and improvements can be made without departing from the concept of the present invention, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the patent of the present invention shall be subject to the attached claims.
Claims
1. A cell modeling method for a two-dimensional braided composite material, characterized in that: include: Obtaining a mesoscopic cell of a woven composite material, and determining the size parameters, warp yarn parameters, and weft yarn parameters of the mesoscopic cell; wherein the size parameters include: cell length Length, cell width Width, and cell thickness Thickness; The size parameters, warp yarn parameters and weft yarn parameters of the meso-cell are input into the finite element software Abaqus to generate a meso-cell model; the generation of the meso-cell model includes the following steps: Create a three-dimensional solid part according to the size parameters of the mesoscopic cell; the length of the three-dimensional solid part is 2×Length, the width is 2×Width, and the thickness is 2×Thickness; Taking the YZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, the contour of the warp fiber bundle is determined according to the size parameters and the warp parameters; taking the XZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, the contour of the weft fiber bundle is determined according to the size parameters and the weft parameters; taking the x=2×Length plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, the weft sweep curve is determined according to the size parameters and the weft parameters; taking the y=2×Width plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, the warp sweep curve is determined according to the size parameters and the warp parameters; Sweep the warp fiber bundle contour and the weft fiber bundle contour along the warp sweep curve and the weft sweep curve respectively, and divide the three-dimensional solid part into the warp fiber bundle, the weft fiber bundle and the matrix; The warp fiber bundles, weft fiber bundles and the parts of the matrix with 0≤x≤Length, 0≤y≤Width, and 0≤z≤Thickness are retained to obtain a mesoscopic cell model of the woven composite material; wherein x is the length of the mesoscopic cell model, y is the width of the mesoscopic cell model, and z is the thickness of the mesoscopic cell model.
2. The cell modeling method of the two-dimensional braided composite material according to claim 1, characterized in that: The warp yarn parameters of the mesoscopic cell include: the warp yarn major axis WarpL, the warp yarn minor axis WarpS, and the warp yarn curvature radius WarpR; The weft yarn parameters of the micro-cell include: the weft yarn major axis WeftL, the weft yarn minor axis WeftS, and the weft yarn curvature radius WeftR.
3. The cell modeling method of the two-dimensional braided composite material according to claim 2, characterized in that: The method uses the YZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, and determines the warp fiber bundle contour according to the size parameters and the warp yarn parameters, which specifically includes: The side point (0, Width / 2, Thickness / 2) is used as the origin of the sketch coordinates, (0, -Thickness / 4) is used as the center coordinates of the ellipse contour, (WarpL / 2, -Thickness / 4) is used as the endpoint of the major axis of the ellipse contour, and (0, WarpS / 2-Thickness / 4) is used as the endpoint of the minor axis of the ellipse contour to obtain the contour of the first warp fiber bundle; The side point (0, Width / 2, Thickness / 2) is used as the origin of the sketch coordinates, (-Width / 2, Thickness / 4) is used as the center coordinates of the ellipse contour, (WarpL / 2-Width / 2, Thickness / 4) is used as the endpoint of the major axis of the ellipse contour, and (-Width / 2, WarpS / 2+Thickness / 4) is used as the endpoint of the minor axis of the ellipse contour, to obtain the contour of the second warp fiber bundle; Take the side point (0, Width / 2, Thickness / 2) as the sketch coordinate origin, (Width / 2, Thickness / 4) as the center coordinates of the ellipse contour, (WarpL / 2+Width / 2, Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (Width / 2, WarpS / 2+Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the third warp fiber bundle.
4. The cell modeling method of the two-dimensional braided composite material according to claim 2, characterized in that: The method uses the part of the XZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch, and determines the contour of the weft fiber bundle according to the size parameters and the weft yarn parameters, which specifically includes: The side point (Length / 2, 0, Thickness / 2) is used as the origin of the sketch coordinates, (0, Thickness / 4) is used as the center coordinates of the ellipse contour, (WeftL / 2, Thickness / 4) is used as the endpoint of the major axis of the ellipse contour, and (0, WeftS / 2+Thickness / 4) is used as the endpoint of the minor axis of the ellipse contour to obtain the contour of the first weft fiber bundle; Take the side point (Length / 2, 0, Thickness / 2) as the origin of the sketch coordinates, (-Length / 2, -Thickness / 4) as the center coordinates of the ellipse contour, (WeftL / 2-Length / 2, -Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (-Length / 2, WeftS / 2-Thickness / 4) as the endpoint of the minor axis of the ellipse contour, and obtain the contour of the second weft fiber bundle; Take the side point (Length / 2, 0, Thickness / 2) as the sketch coordinate origin, (Length / 2, -Thickness / 4) as the center coordinate of the ellipse contour, (WeftL / 2+Length / 2, -Thickness / 4) as the endpoint of the major axis of the ellipse contour, and (Length / 2, WeftS / 2-Thickness / 4) as the endpoint of the minor axis of the ellipse contour to obtain the contour of the third weft fiber bundle.
5. The cell modeling method of a two-dimensional braided composite material according to claim 2, characterized in that: The method uses the part of the x=2×Length plane of the three-dimensional coordinate system of the finite element software Abaqus as a sketch, and determines the weft sweep curve according to the size parameters and the weft parameters, which specifically includes: Take the side point (2×Length, Width / 2, Thickness / 2) as the origin of the sketch coordinates, and draw two circles with (-Width / 2, WeftR-Thickness / 4) and (0, Thickness / 4- WeftR) as the centers and WeftR as the radius. Draw a common tangent line between the two circles, cut the redundant curve, and mirror-copy the remaining curve with the Y axis of the three-dimensional coordinate system of the finite element software Abaqus as the symmetry axis to obtain the first weft sweep curve; Then, the X-axis of the three-dimensional coordinate system of the finite element software Abaqus is used as the symmetry axis, and the first weft yarn sweep curve is mirror-copied to obtain the first weft yarn sweep curve and the second weft yarn sweep curve.
6. The cell modeling method of the two-dimensional braided composite material according to claim 2, characterized in that: The method uses the part of the three-dimensional solid part in the y=2×Width plane of the three-dimensional coordinate system of the finite element software Abaqus as a sketch, and determines the warp sweep curve according to the size parameters and the warp parameters, which specifically includes: Take the side point (Length / 2, 2×Width, Thickness / 2) as the origin of the sketch coordinate system, and draw two circles with (-Length / 2, WarpR-Thickness / 4) and (0, Thickness / 4-WarpR) as the centers and WarpR as the radius. Draw a common tangent line between the two circles, cut the redundant curve, and mirror-copy the remaining curve with the Y axis of the three-dimensional coordinate system of the finite element software Abaqus as the symmetry axis to obtain the first warp sweep curve; The X-axis of the three-dimensional coordinate system of the finite element software Abaqus is used as the symmetry axis, and the first warp sweep curve is mirror-copied to obtain the first warp sweep curve and the second warp sweep curve.
7. A cell modeling system for two-dimensional braided composite materials, characterized in that: include: A cell acquisition module is used to acquire the microscopic cells of the woven composite material and determine the size parameters, warp parameters and weft parameters of the microscopic cells; wherein the size parameters include: cell length Length, cell width Width, cell thickness Thickness; A model building module is used to input the size parameters, warp yarn parameters and weft yarn parameters of the mesoscopic cell into the finite element software Abaqus to generate a mesoscopic cell model; the model building module includes: The entity creation unit is used to create a three-dimensional entity part according to the size parameters of the mesoscopic cell; the length of the three-dimensional entity part is 2×Length, the width is 2×Width, and the thickness is 2×Thickness; The curve construction unit is used to use the portion of the YZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch to determine the contour of the warp fiber bundle according to the size parameters and the warp parameters; use the portion of the XZ plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch to determine the contour of the weft fiber bundle according to the size parameters and the weft parameters; use the portion of the x=2×Length plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch to determine the weft sweep curve according to the size parameters and the weft parameters; use the portion of the y=2×Width plane of the three-dimensional solid part in the three-dimensional coordinate system of the finite element software Abaqus as a sketch to determine the warp sweep curve according to the size parameters and the warp parameters; A sweeping unit is used to sweep the warp fiber bundle contour and the weft fiber bundle contour along the warp sweeping curve and the weft sweeping curve respectively, so as to divide the three-dimensional solid part into the warp fiber bundle, the weft fiber bundle and the matrix; The model building unit is used to retain the parts of the warp fiber bundles, the weft fiber bundles and the matrix with 0≤x≤Length, 0≤y≤Width, and 0≤z≤Thickness, so as to obtain a mesoscopic cell model of the woven composite material; wherein x is the length of the mesoscopic cell model, y is the width of the mesoscopic cell model, and z is the thickness of the mesoscopic cell model.
8. A computer device comprising a memory and a processor, wherein a computer program is stored in the memory, wherein: When the processor executes the computer program, the steps of the two-dimensional woven composite material cell modeling method described in any one of claims 1-6 are implemented.
9. A storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the two-dimensional woven composite material cell modeling method described in any one of claims 1 to 6 are implemented.
Citation Information
Patent Citations
Design method of woven fabric composite material properties based on controllable microstructure
CN106649984A
Modeling method of 2.5-dimensional woven composite material microstructure fiber bundle
CN114239276A