OPTIMIZED DESIGNING METHOD FOR LAMINATE PREFORM OF CERAMIC MATRIX COMPOSITE

Abstract

The present disclosure provides an optimized designing method for a laminate preform of a ceramic matrix composite e. With overall consideration of a strength requirement of a component, a geometric shape and properties of the laminate preform, the method includes: optimizing, based on a corresponding mechanical formula, a fiber volume fraction of each laminate constituting the preform and a fiber direction in the laminate, and selecting a preferred microscopic structure for each laminate, thereby taking full advantage of material performance. The method is applicable for optimized design of various components of the ceramic matrix composite.

Claims

1. An optimized designing method for a laminate preform of a ceramic matrix composite, comprising: step 1: creating, based on a strength analysis result of a macroscopic model of a ceramic matrix composite component, a preliminary design solution to a laminate preform structure; step 2: carrying out finite element strength analysis on the laminate preform structure in step 1, to obtain corresponding first principal stress, normal stress and shear stress of each finite element unit; step 3: calculating, according to a stress result obtained in step 2, an included angle between the first principal stress of each finite element unit and a reference rectangular coordinate system; step 4: utilizing the included angle of each finite element unit obtained in step 3 to change an element coordinate system of each finite element unit, such that a material principal direction of each finite element unit is consistent with a direction of the first principal stress; step 5, utilizing a rule of mixture for a composite to calculate a material parameter of each finite element unit; step 6: carrying out, according to the material parameter of each finite element unit obtained in step 5, strength analysis on the laminate preform structure in step 1 again, and returning to step 1 to redesign the laminate preform structure when a strength requirement is not met; and step 7: selecting, according to a transverse strength calculation formula of the rule of mixture for a composite in conjunction with a strength analysis result obtained in step 6, an appropriate microscopic structure type for each laminate.

2. The optimized designing method for a laminate preform of a ceramic matrix composite according to claim 1, wherein in the step 1, the method further comprises: carrying out finite element strength analysis on the macroscopic model of the ceramic matrix composite component to extract a position of a unit in a stress distribution boundary region, and importing the position of the unit into the macroscopic model for reflection of the position, whereby the laminate preform structure of the ceramic matrix composite component is preliminarily designed.

3. The optimized designing method for a laminate preform of a ceramic matrix composite according to claim 1, wherein in the step 2, the corresponding first principal stress σ.sub.1, normal stresses σ.sub.x and σ.sub.y, and shear stresses τ.sub.yx, τ.sub.zx and τ.sub.zy of each finite element unit are obtained, σ.sub.x and σ.sub.y denoting normal stresses in an X-axis direction and a Y-axis direction respectively, τ.sub.yx denoting a shear stress in a direction perpendicular to a Y axis and parallel to an X axis, τ.sub.zx denoting a shear stress in a direction perpendicular to a Z axis and parallel to the X axis, τ.sub.zy denoting a shear stress in a direction perpendicular to the Z axis and parallel to the Y axis, and X, Y and Z being three coordinate axes in the reference rectangular coordinate system.

4. The optimized designing method for a laminate preform of a ceramic matrix composite according to claim 3, wherein in the step 3, the included angle is one between a first principal stress direction of each finite element unit and each of the three coordinate axes of the reference rectangular coordinate system.

5. The optimized designing method for a laminate preform of a ceramic matrix composite according to claim 3, wherein in the step 3, the following formula of a state of stress at a given point is utilized to obtain the included angle between the first principal stress direction of each finite element unit and each of the three coordinate axes of the reference rectangular coordinate system: $l_{1} (σ_{x} - σ_{1}) + m_{1} τ_{y x} + n_{1} τ_{z x} = 0 l_{1} τ_{x y} + m_{1} (σ_{y} - σ_{1}) + n_{1} τ_{{zy}^{=}} 0 l_{1} = \frac{1}{\sqrt{1 + {(\frac{m_{1}}{l_{1}})}^{2} + {(\frac{n_{1}}{l_{1}})}^{2}}} m_{1} = \frac{1}{\sqrt{1 + {(\frac{l_{1}}{m_{1}})}^{2} + {(\frac{n_{1}}{m_{1}})}^{2}}} n_{1} = \frac{1}{\sqrt{1 + {(\frac{l_{1}}{n_{1}})}^{2} + {(\frac{m_{1}}{n_{1}})}^{2}}}$ wherein l.sub.1, m.sub.1 and n.sub.1 are the included angles between the first principal stress and the X axis, the first principal stress and the Y axis, and the first principal stress and the Z axis of the reference rectangular coordinate system, respectively.

6. The optimized designing method for a laminate preform of a ceramic matrix composite according to claim 4, wherein in the step 3, the following formula of a state of stress at a given point is utilized to obtain the included angle between the first principal stress direction of each finite element unit and each of the three coordinate axes of the reference rectangular coordinate system: $l_{1} (σ_{x} - σ_{1}) + m_{1} τ_{y x} + n_{1} τ_{z x} = 0 l_{1} τ_{x y} + m_{1} (σ_{y} - σ_{1}) + n_{1} τ_{{zy}^{=}} 0 l_{1} = \frac{1}{\sqrt{1 + {(\frac{m_{1}}{l_{1}})}^{2} + {(\frac{n_{1}}{l_{1}})}^{2}}} m_{1} = \frac{1}{\sqrt{1 + {(\frac{l_{1}}{m_{1}})}^{2} + {(\frac{n_{1}}{m_{1}})}^{2}}} n_{1} = \frac{1}{\sqrt{1 + {(\frac{l_{1}}{n_{1}})}^{2} + {(\frac{m_{1}}{n_{1}})}^{2}}}$ wherein l.sub.1, m.sub.1 and n.sub.1 are the included angles between the first principal stress and the X axis, the first principal stress and the Y axis, and the first principal stress and the Z axis of the reference rectangular coordinate system, respectively.

7. The optimized designing method for a laminate preform of a ceramic matrix composite according to claim 1, wherein in the step 5, the material parameters comprise a fiber volume fraction, an elastic modulus, a Poisson's ratio, a shear modulus, and transverse strength.

8. The optimized designing method for a laminate preform of a ceramic matrix composite according to claim 1, wherein in the step 5, on the premise that the first principal stress, fiber strength and matrix strength are known, the following rule of mixture for a composite is utilized to obtain the material parameter of each finite element unit: a fiber volume fraction: $V_{f} = \frac{σ_{1} - σ_{m}}{σ_{f} - σ_{m}}$ wherein V.sub.f is the fiber volume fraction, σ.sub.1 is the first principal stress, σ.sub.f is a maximum stress borne by a reinforcement fiber, and σ.sub.m is a maximum stress borne by a matrix; an elastic modulus: $E_{1} = E_{f} V_{f} + E_{m} (1 - V_{f}) E_{2} = E_{3} = \frac{E_{m} E_{f}}{E_{m} V_{f} + E_{f} (1 - V_{f})}$ wherein E.sub.1, E.sub.2 and E.sub.3 are elastic moduli in three directions 1, 2 and 3 of the composite, a direction 1 being a fiber reinforcement direction, and directions 2 and 3 being two directions perpendicular to the fiber reinforcement direction, E.sub.f is an elastic modulus of the reinforcement fiber, and E.sub.m is an elastic modulus of the matrix; a Poisson's ratio: $v_{1 2} = v_{1 3} = v_{f} V_{f} + v_{m} (1 - V_{f}) v_{2 3} = \frac{E_{2}}{2 G_{2 3}} - 1$ wherein v.sub.12, v.sub.13 and v.sub.23 are Poisson's ratios in the three directions of the composite, v.sub.12 denoting strain, caused by unit strain in the direction 1, in the direction 2, v.sub.13 denoting strain, caused by unit strain in the direction 1, in the direction 3, v.sub.23 denoting strain, caused by unit strain in the direction 2, in the direction 3, of is a Poisson's ratio of the reinforcement fiber, and v.sub.m is a Poisson's ratio of the matrix; a shear modulus: $G_{1 2} = G_{1 3} = \frac{G_{m} G_{f}}{G_{m} V_{f} + G_{f} (1 - V_{f})} G_{2 3} = \frac{G_{f 2 3} G_{m}}{G_{f 2 3} (1 - V_{f}) + G_{m} V_{f}}$ wherein G.sub.12, G.sub.13 and G.sub.23 are shear moduli in three directions of the composite, and denote ratios of a shear stress to shear strain in the directions 1 to 2, the directions 1 to 3 and the directions 2 to 3, G.sub.f is a shear modulus of the reinforcement fiber, G.sub.m is a shear modulus of the matrix, and G.sub.f23 is a shear modulus in the directions 2 to 3 of the reinforcement fiber; and transverse strength:
σ.sub.2=σ.sub.m[1−2(V.sub.f/π).sup.1/2] wherein σ.sub.2 is the transverse strength and σ.sub.m is matrix strength.

9. The optimized designing method for a laminate preform of a ceramic matrix composite according to claim 7, wherein in the step 5, on the premise that the first principal stress, fiber strength and matrix strength are known, the following rule of mixture for a composite is utilized to obtain the material parameter of each finite element unit: a fiber volume fraction: $V_{f} = \frac{σ_{1} - σ_{m}}{σ_{f} - σ_{m}}$ wherein V.sub.f is the fiber volume fraction, σ.sub.1 is the first, principal stress, σ.sub.f is a maximum stress borne by a reinforcement fiber, and σ.sub.m is a maximum stress borne by a matrix; an elastic modulus: $E_{1} = E_{f} V_{f} + E_{m} (1 - V_{f}) E_{2} = E_{3} = \frac{E_{m} E_{f}}{E_{m} V_{f} + E_{f} (1 - V_{f})}$ wherein E.sub.1, E.sub.2 and E.sub.3 are elastic moduli in three directions 1, 2 and 3 of the composite, a direction 1 being a fiber reinforcement direction, and directions 2 and 3 being two directions perpendicular to the fiber reinforcement direction, E.sub.f is an elastic modulus of the reinforcement fiber, and E.sub.m is an elastic modulus of the matrix; a Poisson's ratio: $v_{1 2} = v_{1 3} = v_{f} V_{f} + v_{m} (1 - V_{f}) v_{2 3} = \frac{E_{2}}{2 G_{2 3}} - 1$ wherein v.sub.12, v.sub.13 and v.sub.23 are Poisson's ratios in the three directions of the composite, v.sub.12 denoting strain, caused by unit strain in the direction 1, in the direction 2, v.sub.13 denoting strain, caused by unit strain in the direction 1, in the direction 3, v.sub.23 denoting strain, caused by unit, strain in the direction 2, in the direction 3, v.sub.f is a Poisson's ratio of the reinforcement fiber, and v.sub.m is a Poisson's ratio of the matrix: a shear modulus: $G_{1 2} = G_{1 3} = \frac{G_{m} G_{f}}{G_{m} V_{f} + G_{f} (1 - V_{f})} G_{2 3} = \frac{G_{f 2 3} G_{m}}{G_{f 2 3} (1 - V_{f}) + G_{m} V_{f}}$ wherein G.sub.12, G.sub.13 and G.sub.23 are shear moduli in three directions of the composite, and denote ratios of a shear stress to shear strain in the directions 1 to 2, the directions 1 to 3 and the directions 2 to 3, G.sub.f is a shear modulus of the reinforcement fiber, G.sub.m is a shear modulus of the matrix, and G.sub.f23 is a shear modulus in the directions 2 to 3 of the reinforcement fiber; and transverse strength:
σ.sub.2=σ.sub.m[1−2(V.sub.f/π).sup.1/2] wherein σ.sub.2 is the transverse strength and am is matrix strength.

10. The optimized designing method for a laminate preform of a ceramic matrix composite according to claim 1, wherein in the step 7, a microscopic structure type comprises a unidirectional reinforcement structure, a plain weave structure and a three-dimensional weave structure.

11. A turbine rotor blade of an engine, wherein a laminate preform structure of a blade body of the turbine rotor blade of the engine is formed by stacking a SiC fiber cloth, the blade body is divided into three parts of a middle, a concave and a convex, according to a thickness of a single layer of SiC fiber cloth, a solution to a laminate of each part is refined, a blade body structure is obtained, and the laminate preform structure of the blade body is formed by stacking 14 SiC fiber cloths of different shapes, laminates 1-4 using a plain weave structure, laminates 5-8 using a unidirectional reinforcement structure, fiber volume fractions in laminates 1-6 being 10.4%, 10.0%, 8.5%, 8.1%, 5.0% and 6.8%, respectively, and fiber volume fractions in laminates 7-14 being all 5.0%.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0031] FIG. 1 is a flowchart of a design solution for optimizing a laminate preform;

[0032] FIG. 2 is a diagram showing a position distribution of a stress boundary unit at a blade body portion of a turbine rotor blade; and

[0033] FIG. 3 is a schematic diagram of a design solution to a laminate preform of a blade body of the turbine rotor blade.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0034] The present disclosure is described in further detail below with reference to the accompanying drawings.

[0035] As shown in FIG. 1, an optimized designing method for a laminate preform particularly includes the steps as follows.

[0036] In step 1: a preliminary design solution is created based on a strength analysis result of a macroscopic model of a ceramic matrix composite component. Finite element strength analysis is carried out on the macroscopic model of the ceramic matrix composite component to extract a position of a unit in a stress distribution boundary region, and the position is reflected in a component model to complete preliminary design of a laminate preform structure of the component.

[0037] In step 2: a stress result after a finite element analysis is extracted. Finite element strength analysis is carried out on a designed laminate preform model of the component in step 1, to obtain the corresponding first principal stress σ.sub.1, normal stresses σ.sub.x and σ.sub.y, and shear stresses τ.sub.yx, τ.sub.zx and τ.sub.zy of each finite element unit.

[0038] In step 3: an included angle between the first principal stress of each finite element unit and a reference coordinate system is derived. A formula of a state of stress at a given point (where, l.sub.1, m.sub.1 and n.sub.1 are included angles between the first principal stress and an X axis, the first principal stress and a Y axis, and the first principal stress and a Z axis of the reference coordinate system, respectively) is utilized to obtain an included angle between a first principal stress direction of each finite element unit and each of three coordinate axes of the reference rectangular coordinate system:

[00006] $l_{1} (σ_{x} - σ_{1}) + m_{1} τ_{y x} + n_{1} τ_{z x} = 0 l_{1} τ_{x y} + m_{1} (σ_{y} - σ_{1}) + n_{1} τ_{{zy}^{=}} 0 l_{1} = \frac{1}{\sqrt{1 + {(\frac{m_{1}}{l_{1}})}^{2} + {(\frac{n_{1}}{l_{1}})}^{2}}} m_{1} = \frac{1}{\sqrt{1 + {(\frac{l_{1}}{m_{1}})}^{2} + {(\frac{n_{1}}{m_{1}})}^{2}}} n_{1} = \frac{1}{\sqrt{1 + {(\frac{l_{1}}{n_{1}})}^{2} + {(\frac{m_{1}}{n_{1}})}^{2}}}$

[0039] In step 4: a material principal direction changes with an element coordinate system in finite element software. The included angle of each finite element unit obtained in step 3 is utilized to change an element coordinate system of each finite element unit, so that a material principal direction of each unit is consistent with the direction of the first principal stress.

[0040] In step 5: given that bearing strength of the composite considered in the present method refers to a maximum stress borne by both a fiber and a matrix, a rule of mixture for a composite may be utilized to obtain a volume fraction (fiber volume fraction) of each finite element unit on the premise that the first principal stress, fiber strength σ.sub.f (that is, a maximum stress borne by a reinforcement fiber) and matrix strength σ.sub.m (that is, a maximum stress borne by the matrix). An influence, on other material parameters, from the volume fraction may also be derived from the rule of mixture for a composite: [0041] a fiber volume fraction is calculated by:

[00007] $V_{f} = \frac{σ_{1} - σ_{m}}{σ_{f} - σ_{m}}$

where, V.sub.f is the fiber volume fraction, σ.sub.1 is the first principal stress, σ.sub.f is the maximum stress borne by the reinforcement fiber, and σ.sub.m is the maximum stress borne by the matrix; [0042] an elastic modulus is calculated by:

[00008] $E_{1} = E_{f} V_{f} + E_{m} (1 - V_{f}) E_{2} = E_{3} = \frac{E_{m} E_{f}}{E_{m} V_{f} + E_{f} (1 - V_{f})}$

where, E.sub.1, E.sub.2 and E.sub.3 are elastic moduli in three directions 1, 2 and 3 of the composite, a direction 1 being a fiber reinforcement direction, and directions 2 and 3 being two directions perpendicular to the fiber reinforcement direction, E.sub.f is an elastic modulus of the reinforcement fiber, and E.sub.m is an elastic modulus of the matrix; [0043] a Poisson's ratio is calculated by:

[00009] $v_{12} = v_{13} = v_{f} V_{f} + v_{m} (1 - V_{f}) v_{23} = \frac{E_{2}}{2 G_{2 3}} - 1$

where, v.sub.12, v.sub.13 and v.sub.23 are Poisson's ratios in the three directions of the composite, of is a Poisson's ratio of the reinforcement fiber, and v.sub.m is a Poisson's ratio of the matrix; [0044] a shear modulus is calculated by:

[00010] $G_{1 2} = G_{1 3} = \frac{G_{m} G_{f}}{G_{m} V_{f} + G_{f} (1 - V_{f})} G_{2 3} = \frac{G_{f 2 3} G_{m}}{G_{f 2 3} (1 - V_{f}) + G_{m} V_{f}}$

where, G.sub.12, G.sub.13 and G.sub.23 are shear moduli in three directions of the composite, G.sub.f is a shear modulus of the reinforcement fiber, G.sub.m is a shear modulus of the matrix, and G.sub.f23 is a shear modulus in a direction 2 to 3 of the reinforcement fiber; and [0045] transverse strength is calculated by:

σ.sub.2=σ.sub.m[1−2(V.sub.f/π).sup.1/2]

where, σ.sub.2 is the transverse strength, and σ.sub.m is matrix strength, the transverse strength being perpendicular to a reinforcement direction.

[0046] In step 6: according to a new material parameter of each finite element unit obtained in step 5, strength analysis is carried out on the laminate preform structure of the ceramic matrix composite component in step 1 again, returning to step 1 to redesign the laminate preform structure when a strength requirement is not met.

[0047] In step 7: a final optimization design solution is determined. A microscopic structure of the composite includes a unidirectional reinforcement structure, a plain weave structure and a three-dimensional weave structure. The unidirectional reinforcement structure has only one reinforcement direction, and has weaker strength in the other two directions; the plain weave structure has two reinforcement directions, and has weaker strength in a thickness direction; and the three-dimensional weave structure has excellent mechanical performance in three directions. According to the transverse strength calculation formula of the rule of mixture for a composite, a calculation result obtained in step 6 is analyzed, and an appropriate microscopic structure type (unidirectional reinforcement, plain weave and three-dimensional weave) is selected for each laminate to meet the strength requirement. To this end, the method provides a macroscopic shape of each laminate of a component preform, selects a preferred microscopic structure type for each laminate, and optimizes a fiber volume fraction of the laminate and a fiber direction (the first principal stress direction) in the laminate.

[0048] Next, the present method is illustrated by taking a complex component, a turbine rotor blade of a certain aero-engine, as an example. Emphasis is put on an optimization result of a laminate preform of its blade body portion. The method particularly includes:

[0049] (1) A laminate preform structure is formed by stacking a SiC fiber cloth, firstly, finite element analysis software is utilized to analyze strength of a macroscopic model of a ceramic matrix turbine rotor blade to extract a coordinate of stress boundary points and further import these points into the macroscopic model of the blade to make them reflected in a component model, as shown in FIG. 2. In consideration of a geometric shape of a blade body and the coordinates of a series of imported stress boundary points, it is found that the blade body may be divided into three parts of a middle, a concave and a convex. Then, a solution to a laminate of each part is refined according to a thickness of a single layer of SiC fiber cloth to finally obtain a blade body structure. As shown in FIG. 3, a laminate preform structure of the blade body portion is formed by stacking 14 SiC fiber cloths of different shapes.

[0050] (2) The finite element analysis software is utilized to analyze strength of the laminate preform structure of the blade body designed in (1), to obtain the corresponding first principal stress σ.sub.1, normal stresses σ.sub.x and σ.sub.y, and shear stresses τ.sub.yx, τ.sub.zx and τ.sub.zy of each finite element unit.

[0051] (3) A formula of a state of stress at a given point in step 3 is utilized to obtain an included angle between a first principal stress direction of each finite element unit and each of three coordinate axes of a reference rectangular coordinate system. In the finite element analysis software, elastic moduli in three directions are defined when a material parameter of a silicon carbide fiber is defined, where a direction with a maximum modulus is a material principal direction. A material principal direction in each finite element unit changes as a direction of a coordinate axis of an element coordinate system changes. A material principal direction of the unit may be changed by changing the direction of the coordinate axis of the element coordinate system, and the element coordinate system may be changed by creating a local coordinate system. In the finite element analysis software, the local coordinate system is created by sequentially rotating X, Y and Z axes of the reference coordinate system, so it is required to know three included angles between the principal direction of each unit and a global coordinate system. A fiber direction is optimized based on the method.

[0052] (4) The formula of the rule of mixture for a composite in step 5 is utilized to determine an optimal volume fraction of the laminate preform structure, subjected to fiber optimization in (3), of the blade body and modify a corresponding material parameter.

[0053] (5) Strength analysis is carried out for a solution after volume fraction optimization in (4) again, proceeding to a next step when a strength requirement is met, otherwise returning to (1).

[0054] (6) A preferred microscopic structure of each laminate is selected according to step 7 to obtain a final solution to optimized design of the laminate preform structure of the blade body. As shown in FIG. 3, after a preliminary design solution to the laminate preform of a blade body component is optimized, laminates 1-4 use the plain weave structure, laminates 5-8 use the unidirectional reinforcement structure, fiber volume fractions in laminates 1-6 are 10.4%, 10.0%, 8.5%, 8.1%, 5.0% and 6.8%, respectively, and fiber volume fractions in laminates 7-14 are all 5.0%. [92] What is described above is merely the preferred implementation of the present disclosure, the scope of protection of the present disclosure is not limited to the above embodiments, and all technical solutions following the idea of the present disclosure fall within the scope of protection of the present disclosure. It should be noted that several modifications and adaptations made by those of ordinary skill in the art without departing from the principle of the present disclosure should fall within the scope of protection of the present disclosure.

OPTIMIZED DESIGNING METHOD FOR LAMINATE PREFORM OF CERAMIC MATRIX COMPOSITE

Inventors

Cpc classification

Classification Explorer

G06F2119/14

PHYSICS

Classification Explorer

G06F2111/10

PHYSICS

Classification Explorer

G06F30/23

PHYSICS

Classification Explorer

G06F2113/26

PHYSICS

International classification

Classification Explorer

G06F30/23

PHYSICS

Abstract

Claims

Description