Three-dimensional parametric modeling method for multi-layer heterogeneous composite cladding tube with two-dimensional fiber woven structure

Abstract

The invention relates to a three-dimensional parametric modeling method for a multi-layer heterogeneous composite cladding tube with a two-dimensional fiber woven structure, which comprises the following steps of: 1, setting a parameter equation of a geometric profile curve of a fiber bundle, and materializing a space coated by the parameter equation; 2, establishing a hexahedron matrix model completely covering a two-dimensional (2D) woven fiber bundle and a hexahedron interface and lining metal model, and carrying out union operation on each component model by utilizing a Boolean operation function; and 3, according to the inner diameter of the actual cladding tube, bending the hexahedron model containing the two-dimensional braided structure into a corresponding radian, and performing circumferential array operation according to the radian to complete modeling of the multi-layer heterogeneous composite cladding tube. According to the parametric modeling method, the limitation that the component model of the multi-layer heterogeneous composite cladding tube is not easy to build completely and the number of layers of braided preforms is uncontrollable can be solved, meanwhile, flexible change can be achieved according to the geometric dimension of the real component model, the modeling time is saved, and accurate simulation of the multi-layer heterogeneous composite cladding tube is achieved.

Description

Technical Field

[0001] The invention belongs to the field of nuclear fuel composite cladding tubes and relates to a three-dimensional parameterized modeling method of a multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braiding structure. Background Art

[0002] The advantage of using simulation technology in computational materials science to study materials is that it can not only simulate various experimental processes to understand the internal microscopic properties of materials and their macroscopic mechanical behavior, but also predict their performance before these new materials are actually prepared, providing strong theoretical guidance for designing new structural materials with excellent performance.

[0003] SiC f / SiC ceramic matrix composite is one of the candidate materials for nuclear cladding tube. The continuous SiC fiber toughening SiC matrix has two ways: winding and weaving. However, under the same fiber volume fraction and fiber angle, the mechanical properties of the SiC ceramic matrix composite with braided structure are better than those of the SiC ceramic matrix composite with wound structure. Therefore, for SiC f / SiC ceramic-based composite cladding tubes in the field of nuclear material applications are mainly focused on the research of SiC fiber braided structures.

[0004] At the same time, in order to improve the current SiC f / SiC composites crack at low stress levels and do not meet the airtightness standards. Therefore, a composite made of refractory metals and SiC was introduced. f The new cladding structure composed of SiC is designed to improve airtightness and increase accident tolerance.

[0005] At present, in the field of simulation design of ceramic matrix composites, finite element methods are mostly used for multi-scale simulation design to predict and analyze the mechanical properties and damage behavior of ceramic matrix composites. Multi-scale simulation design is to use microscopic technology, digital imaging technology, and other methods to use finite element simulation software to reproduce the structure of the material at different scales such as micro, meso and macro, and then simulate the design of the material under different working conditions and physical fields. Summary of the invention

[0006] Technical issues to be solved

[0007] In order to avoid the shortcomings of the prior art, the present invention proposes a three-dimensional parametric modeling method for a multi-layer heterogeneous composite cladding tube with a two-dimensional fiber weaving structure, with the purpose of solving the problem that macroscopic composite cladding tubes with complex 2D fiber weaving structures are difficult to model, and provides a reliable calculation model for accurate finite element simulation analysis of multi-layer heterogeneous composite cladding tubes.

[0008] Technical Solution

[0009] A three-dimensional parametric modeling method for a multilayer heterogeneous composite cladding tube with a two-dimensional (2D) fiber braided structure is characterized by the following steps:

[0010] Step 1: Fiber bundle outline solidification and component models:

[0011] Fiber bundle contour solidification: The geometric contour curves of the 2D braided fiber bundle in different directions are described by the parametric equations of the geometric contour curves. The filling function in the COMSOL finite element software is used to solidify the space covered by the parametric equations to complete the modeling of the 2D braided fiber bundle and obtain the 2D braided fiber bundle model.

[0012] The parametric equation is described using the geometric dimension information of the thickness t, length l, width w of a single fiber bundle and the distance g between two adjacent fiber bundles, and the unit is mm;

[0013] Component modeling: The components include matrix, interface and lining metal. After hexahedron modeling, matrix model, interface model and lining metal model are obtained. The lengths of the three component models are l m and width w m Satisfy m =l=2w+2g,w m =l=2w+2g;

[0014] The number of matrix models is twice that of the 2D braided fiber bundle model, and one side of the hexahedron matches the 2D braided fiber bundle model and completely covers the 2D braided fiber bundle;

[0015] Step 2: Use the Boolean operation function in COMSOL finite element software to perform a union operation on the fiber bundle, matrix, interface and lining metal to complete the hexahedral modeling of the single-layer fiber preform;

[0016] During the union operation, two matrix models completely cover a 2D braided fiber bundle model to form model A. One end of the multiple models A after being stacked and united is the interface model and the lining metal model in turn, and the hexahedral modeling with a single-layer fiber preform is obtained, which is a heterogeneous cladding tube model of a multi-layer fiber preform.

[0017] Step 3, modeling of multi-layer heterogeneous cladding tube: according to the inner diameter of the actual cladding tube, the large deformation function in the COMSOL finite element software is used to bend the multi-layer heterogeneous cladding tube containing multiple layers of two-dimensional braided structure into a designed angle, and a circular array operation is performed according to the arc to complete the modeling of the multi-layer heterogeneous composite cladding tube.

[0018] The parameter equations of the geometric profile curves of the fiber bundle in different directions described in step 1 include four curves of the longitudinal section and the transverse section of the fiber bundle, that is, four parameter equations are set:

[0019]

[0020]

[0021]

[0022]

[0023] Among them, y1(x) and y2(x) are the two curve parameter equations of the longitudinal section of the fiber bundle; y3(x) and y4(x) are the two curve parameter equations of the transverse section of the fiber bundle; x is the independent variable of the parameter equation, the unit is mm; t is the thickness of the fiber bundle, the unit is mm; l is the length of the fiber bundle, the unit is mm; g is the distance between two adjacent fiber bundles, the unit is mm, that is, g=l / 2-w; w is the width of the fiber bundle, the unit is mm.

[0024] When the 2D woven fiber bundle is spatially solidified in step 1, the cross section of the fiber bundle is first solidified using the fill function in the COMSOL software, and then the cross section is used as the source surface and the longitudinal cross-sectional curve of the fiber bundle is used as the path for sweeping and filling to realize the spatial solidification of the 2D woven fiber bundle.

[0025] The Boolean union function described in step 2 is to remove the space occupied by the 2D woven fibers wrapped in the hexahedral matrix to form an irregularly shaped matrix, that is, the matrix and the 2D woven fibers become two independent geometric models.

[0026] The multilayer fiber preform described in step 3 has three layers, and the modeling of the three-layer fiber preform is achieved by arranging two layers of fiber preforms in the thickness direction of the single-layer fiber preform, for a total of three layers of fiber preforms.

[0027] The thicknesses of the models described in step 1 are: 0.205 mm for a single fiber bundle, 0.475 mm for a matrix, 0.100 mm for an interface, and 0.399 mm for a metal lining. The width w of the single fiber bundle model is 0.925 mm, and the length l is 2.159 mm. The width and length of the other component models are the same, both 2.159 mm.

[0028] The angle corresponding to the multi-layer heterogeneous composite cladding tube described in step 3 after large deformation is:

[0029]

[0030] The inner diameter of the cladding tube is 11 mm, the thickness of the inner lining metal is 2.159 mm, and the angle corresponding to the multi-layer heterogeneous composite cladding tube after large deformation is 22.517°.

[0031] The circular array described in step 3 is based on the straight line where the origin is located as the rotation axis, and the array angle is the angle corresponding to the multi-layer heterogeneous composite cladding tube after large deformation.

[0032] The circular array described in step 3 is based on the straight line where the origin is located as the rotation axis, the array angle is 22.517°, the number of arrays is 16, and the establishment of the three-dimensional model of the multi-layer heterogeneous composite cladding tube is completed.

[0033] Beneficial Effects

[0034] The present invention proposes a three-dimensional parametric modeling method for a multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braiding structure, comprising the following steps: 1: setting a parametric equation of the geometric contour curve of the fiber bundle, and materializing the space covered by the parametric equation; 2: establishing a hexahedral matrix model that completely covers the two-dimensional (2D) braided fiber bundle, as well as a hexahedral interface and lining metal model, and using the Boolean operation function to perform a union operation on each component model; 3: according to the inner diameter of the actual cladding tube, the hexahedral model containing the two-dimensional braiding structure is bent into a corresponding arc, and a circular array operation is performed according to the arc to complete the modeling of the multi-layer heterogeneous composite cladding tube. The parametric modeling method can solve the limitations that the component model of the multi-layer heterogeneous composite cladding tube is not easy to establish completely and the number of braided preform layers is uncontrollable. At the same time, it can be flexibly changed according to the geometric size of the real component model, saving modeling time, and realizing accurate simulation and simulation of the multi-layer heterogeneous composite cladding tube.

[0035] The advantages and beneficial effects of the present invention are as follows: the three-dimensional parametric modeling method is used to solve the limitations that the component model of the multi-layer heterogeneous composite cladding tube is difficult to establish completely and the number of braided preform layers is uncontrollable. At the same time, it can be flexibly changed according to the geometric dimensions of the real component model, saving modeling time, realizing accurate simulation of the multi-layer heterogeneous composite cladding tube, and providing a reliable calculation model for realizing finite element simulation analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] Figure 1 : A flow chart of a multi-layer heterogeneous composite cladding tube model established for an embodiment of the present invention;

[0037] Figure 2 : A schematic diagram of a matrix and fiber model after a Boolean operation union established in an embodiment of the present invention;

[0038] Figure 3 : Schematic diagram of the changes of the multi-layer heterogeneous hexahedron model established in the embodiment of the present invention before and after large deformation;

[0039] Figure 4 : Schematic diagram of a multi-layer heterogeneous composite cladding tube model established in an embodiment of the present invention. DETAILED DESCRIPTION

[0040] The present invention will now be further described with reference to the embodiments and the accompanying drawings:

[0041] Embodiment 1:

[0042] S1: Fiber bundle contour solidification: set the parametric equations of the fiber bundle geometric contour curves in different directions respectively. The parametric equations can describe the geometric dimension information of the fiber bundle, such as thickness (t), length (l), width (w), etc., and then use the filling function in COMSOL finite element software to solidify the space covered by the parametric equation to complete the modeling of the 2D braided fiber bundle;

[0043] The parameter equations of the longitudinal section and transverse section of the fiber bundle, a total of 4 curves, that is, 4 parameter equations are set:

[0044]

[0045]

[0046]

[0047]

[0048] Among them, y1(x) and y2(x) are two curve parameter equations of the longitudinal section of the fiber bundle, y3(x) and y4(x) are two curve parameter equations of the transverse section of the fiber bundle, x is the independent variable of the parameter equation (mm), t is the thickness of the fiber bundle (mm), l is the length of the fiber bundle (mm), g is the distance between two adjacent fiber bundles (mm), that is, g=l / 2-w, w is the width of the fiber bundle (mm).

[0049] The fiber bundle cross section is then solidified using the filling function in COMSOL finite element software. The cross section is then used as the source surface and the longitudinal cross-sectional curve of the fiber bundle is used as the path for sweeping and filling to realize the spatial solidification of the 2D braided fiber bundle.

[0050] S2: Modeling of other component models (matrix, interface, lining metal): According to the length and thickness of the single-layer 2D woven fiber bundle established in step S1, a hexahedral matrix model that completely covers the fiber bundle and a hexahedral interface and lining metal model are established.

[0051] The thicknesses of the fiber bundle, matrix, interface and metal lining are respectively 0.205 mm for the fiber bundle, 0.475 mm for the matrix, 0.100 mm for the interface and 0.399 mm for the metal lining. The width of a single fiber bundle model is 0.925 mm and the length is 2.159 mm. The width and length of the other component models are the same, both 2.159 mm.

[0052] Then, the Boolean operation function in COMSOL finite element software is used to remove the space occupied by the 2D woven fibers wrapped in the hexahedral matrix to form an irregularly shaped matrix, that is, the matrix and the 2D woven fibers become two independent geometric models, such as Figure 2 As shown, the hexahedral modeling with a single-layer fiber preform is completed;

[0053] S3: Modeling of multi-layer heterogeneous cladding tube: Based on the heterogeneous cladding tube model of the single-layer fiber preform established in step S2, two layers of fiber preforms are arrayed in the thickness direction, for a total of three layers of fiber preforms, to achieve modeling of a multi-layer fiber preform.

[0054] According to the actual inner diameter of the cladding tube (11 mm), the angle corresponding to the multi-layer heterogeneous composite cladding tube with lining metal (l = 2.159 mm) as the inner wall after large deformation is:

[0055]

[0056] The established model is rotated with the origin line as the axis of rotation and the array angle is 22.517°, such as Figure 3 As shown, the number of arrays is 16, and the three-dimensional model of the multi-layer heterogeneous composite cladding tube is established. Figure 4 shown.

[0057] Referring to Example 1 of the present invention, the geometric parameters of different component models can be set according to actual observation results, and the multi-layer heterogeneous composite cladding tube model established by the present invention can be realized. Therefore, the present invention will no longer list and explain them one by one.

Claims

1. A three-dimensional parametric modeling method for a multilayer heterogeneous composite cladding tube with a two-dimensional (2D) fiber braided structure, characterized in that Here are the steps: Step 1: Fiber bundle outline solidification and component models: Fiber bundle contour solidification: The geometric contour curves of the 2D braided fiber bundle in different directions are described by the parametric equations of the geometric contour curves. The filling function in the COMSOL finite element software is used to solidify the space covered by the parametric equations to complete the modeling of the 2D braided fiber bundle and obtain the 2D braided fiber bundle model. The parametric equation is described using the geometric dimension information of the thickness t, length l, width w of a single fiber bundle and the distance g between two adjacent fiber bundles, and the unit is mm; Component modeling: The components include matrix, interface and lining metal. After hexahedron modeling, matrix model, interface model and lining metal model are obtained. The lengths of the three component models are l m and width w m Satisfy m =l=2w+2g,w m =l=2w+2g; The number of matrix models is twice that of the 2D braided fiber bundle model, and one side of the hexahedron matches the 2D braided fiber bundle model and completely covers the 2D braided fiber bundle; Step 2: Use the Boolean operation function in COMSOL finite element software to perform a union operation on the fiber bundle, matrix, interface and lining metal to complete the hexahedral modeling of the single-layer fiber preform; During the union operation, two matrix models completely cover a 2D braided fiber bundle model to form model A. One end of the multiple models A after being stacked and united is the interface model and the lining metal model in turn, and the hexahedral modeling with a single-layer fiber preform is obtained, which is a heterogeneous cladding tube model of a multi-layer fiber preform. Step 3, modeling of multi-layer heterogeneous cladding tube: according to the inner diameter of the actual cladding tube, the large deformation function in the COMSOL finite element software is used to bend the multi-layer heterogeneous cladding tube containing multiple layers of two-dimensional braided structure into a designed angle, and a circular array operation is performed according to the arc to complete the modeling of the multi-layer heterogeneous composite cladding tube.

2. The three-dimensional parametric modeling method of the multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braided structure according to claim 1 is characterized in that: The parameter equations of the geometric profile curves of the fiber bundle in different directions described in step 1 include four curves of the longitudinal section and the transverse section of the fiber bundle, that is, four parameter equations are set: Among them, y1(x) and y2(x) are the two curve parameter equations of the longitudinal section of the fiber bundle; y3(x) and y4(x) are the two curve parameter equations of the transverse section of the fiber bundle; x is the independent variable of the parameter equation, the unit is mm; t is the thickness of the fiber bundle, the unit is mm; l is the length of the fiber bundle, the unit is mm; g is the distance between two adjacent fiber bundles, the unit is mm, that is, g=l / 2-w; w is the width of the fiber bundle, the unit is mm.

3. The three-dimensional parametric modeling method of the multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braided structure according to claim 1 is characterized in that: When the 2D woven fiber bundle is spatially solidified in step 1, the cross section of the fiber bundle is first solidified using the fill function in the COMSOL software, and then the cross section is used as the source surface and the longitudinal cross-sectional curve of the fiber bundle is used as the path for sweeping and filling to realize the spatial solidification of the 2D woven fiber bundle.

4. The three-dimensional parametric modeling method of the multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braided structure according to claim 1, characterized in that: The Boolean union function described in step 2 is to remove the space occupied by the 2D woven fibers wrapped in the hexahedral matrix to form an irregularly shaped matrix, that is, the matrix and the 2D woven fibers become two independent geometric models.

5. The three-dimensional parametric modeling method of the multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braided structure according to claim 1, characterized in that: The multilayer fiber preform described in step 3 has three layers, and the modeling of the three-layer fiber preform is achieved by arranging two layers of fiber preforms in the thickness direction of the single-layer fiber preform, for a total of three layers of fiber preforms.

6. The three-dimensional parametric modeling method of the multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braided structure according to claim 1, characterized in that: The thicknesses of the models described in step 1 are: 0.205 mm for a single fiber bundle, 0.475 mm for a matrix, 0.100 mm for an interface, and 0.399 mm for a metal lining. The width w of the single fiber bundle model is 0.925 mm, and the length l is 2.159 mm. The width and length of the other component models are the same, both 2.159 mm.

7. The three-dimensional parametric modeling method of the multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braided structure according to claim 1, characterized in that: The angle corresponding to the multi-layer heterogeneous composite cladding tube described in step 3 after large deformation is:

8. The three-dimensional parametric modeling method of the multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braided structure according to claim 7, characterized in that: The inner diameter of the cladding tube is 11 mm, the thickness of the inner lining metal is 2.159 mm, and the angle corresponding to the multi-layer heterogeneous composite cladding tube after large deformation is 22.517°.

9. The three-dimensional parametric modeling method of the multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braided structure according to claim 1 or 8, characterized in that: The circular array described in step 3 is based on the straight line where the origin is located as the rotation axis, and the array angle is the angle corresponding to the multi-layer heterogeneous composite cladding tube after large deformation.

10. The three-dimensional parametric modeling method of the multi-layer heterogeneous composite cladding tube with a two-dimensional fiber braided structure according to claim 7 or 8, characterized in that: The circular array described in step 3 is based on the straight line where the origin is located as the rotation axis, the array angle is 22.517°, the number of arrays is 16, and the establishment of the three-dimensional model of the multi-layer heterogeneous composite cladding tube is completed.

Images