A matrix crack density prediction method, device, equipment, medium and product

Abstract

The present application discloses a matrix crack density prediction method, device, equipment, medium and product, which relate to the field of mechanical property calculation of composite materials. The method comprises the following steps: establishing a time- and temperature-dependent interface phase oxidation model of a fiber-reinforced ceramic matrix composite material based on a thermal-mechanical-oxygen coupling mechanism; establishing a composite material unit cell model based on a damage state, and determining a time- and temperature-dependent internal fiber and matrix microscopic stress field of the composite material; using a fracture mechanics interface debonding criterion to determine a temperature-, stress- and time-dependent interface debonding length; determining a temperature-, stress- and time-dependent matrix strain energy, and using an initial matrix cracking stress and crack spacing to determine a critical matrix strain energy; using a critical matrix strain energy criterion to predict a temperature-, stress- and time-dependent matrix crack density based on the critical matrix strain energy. The present application can accurately predict the matrix crack density of a fiber-reinforced ceramic matrix composite material under thermal-mechanical-oxygen coupling.

Description

Technical Field

[0001] The present application relates to the field of mechanical property calculation of composite materials, and in particular to a matrix crack density prediction method, device, equipment, medium and product. Background Art

[0002] Fiber-reinforced ceramic matrix composites can withstand complex stress loads and serve for a long time in extreme environments such as high temperature and chemical corrosion. With the rapid development of high-performance engines, hypersonic aircraft, nuclear power and other fields, fiber-reinforced ceramic matrix composites have become the key material basis for various high-end equipment and have received extensive attention from countries around the world.

[0003] The service environment of fiber-reinforced ceramic matrix composites is extremely harsh, including complex stress states, high temperatures, chemical corrosion, etc. With the increase of ambient temperature, the increase of use load, and the effect of the oxidizing environment, the modulus and strength gradually decline over time, and the ability to resist deformation becomes worse. Under the action of heat-mechanics-oxygen coupling, the matrix crack density gradually increases, allowing oxidizing gas components to enter the interior of the composite material, resulting in increased damage. In order to ensure the reliability and safety of fiber-reinforced ceramic matrix composites during use, it is necessary to establish a matrix crack density prediction model for fiber-reinforced ceramic matrix composites under the action of heat-mechanics-oxygen coupling, and predict the evolution relationship between crack surface density and temperature, stress and time.

[0004] However, there are few studies on the matrix crack density of fiber-reinforced ceramic matrix composites under the action of thermal-mechanical-oxygen coupling, and it is impossible to accurately predict the matrix crack density of fiber-reinforced ceramic matrix composites under the action of thermal-mechanical-oxygen coupling. How to consider the coupling effects of multiple factors such as temperature, load and oxidation, and predict the matrix crack density of fiber-reinforced ceramic matrix composites, is a key technical problem that needs to be solved to ensure the actual engineering application of fiber-reinforced ceramic matrix composites. Summary of the invention

[0005] The purpose of this application is to provide a matrix crack density prediction method, device, equipment, medium and product to solve the problem of being unable to accurately predict the matrix crack density of fiber-reinforced ceramic matrix composites under the action of heat-mechanics-oxygen coupling.

[0006] To achieve the above objectives, this application provides the following solutions.

[0007] In a first aspect, the present application provides a method for predicting matrix crack density, comprising the following steps.

[0008] Based on the thermal-mechanical-oxygen coupling mechanism, a time- and temperature-dependent interfacial phase oxidation model of fiber-reinforced ceramic matrix composites was established.

[0009] Based on the damage state of fiber-reinforced ceramic matrix composites under the coupling of heat, force and oxygen, a composite unit cell model is established according to the interface phase oxidation model, and the time and temperature dependent microscopic stress field of the internal fibers and the matrix of the composite material is determined according to the composite unit cell model.

[0010] Based on the microscopic stress field between the fibers and the matrix inside the composite material, the interface debonding criterion of fracture mechanics is adopted to determine the temperature, stress and time dependent interface debonding length.

[0011] Based on the microscopic stress field between the fibers and the matrix inside the composite material and the interface debonding length, the temperature, stress and time-dependent matrix strain energy is determined, and based on the temperature, stress and time-dependent matrix strain energy, the critical matrix strain energy is determined using the initial matrix cracking stress and crack spacing.

[0012] A critical matrix strain energy criterion is adopted to predict the temperature, stress and time dependent matrix crack density based on the critical matrix strain energy.

[0013] In a second aspect, the present application provides a matrix crack density prediction device, comprising the following modules.

[0014] The interface phase oxidation model building module is used to establish the time and temperature dependent interface phase oxidation model of fiber reinforced ceramic matrix composites based on the thermal-mechanical-oxygen coupling mechanism.

[0015] The module for determining the microscopic stress field of the fibers and the matrix inside the composite material is used to establish a composite material unit cell model based on the damage state of the fiber-reinforced ceramic matrix composite material under the action of heat-mechanics-oxygen coupling and the interface phase oxidation model, and to determine the time and temperature dependent microscopic stress field of the fibers and the matrix inside the composite material based on the composite material unit cell model.

[0016] The interface debonding length determination module is used to determine the temperature, stress and time-dependent interface debonding length based on the microscopic stress field between the fiber and the matrix inside the composite material and the fracture mechanics interface debonding criterion.

[0017] The critical matrix strain energy determination module is used to determine the temperature, stress and time-dependent matrix strain energy based on the microscopic stress field between the fibers and the matrix inside the composite material and the interface debonding length, and to determine the critical matrix strain energy based on the temperature, stress and time-dependent matrix strain energy using the initial matrix cracking stress and crack spacing.

[0018] The matrix crack density prediction module is used to predict the temperature, stress and time dependent matrix crack density based on the critical matrix strain energy using a critical matrix strain energy criterion.

[0019] In a third aspect, the present application provides a computer device, comprising: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement any of the above-described matrix crack density prediction methods.

[0020] In a fourth aspect, the present application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements any of the above-described matrix crack density prediction methods.

[0021] In a fifth aspect, the present application provides a computer program product, including a computer program, which, when executed by a processor, implements any of the above-mentioned matrix crack density prediction methods.

[0022] According to the specific embodiments provided by the present application, the present application discloses the following technical effects: based on the heat-mechanics-oxygen coupling mechanism, the present application establishes a time- and temperature-dependent interface phase oxidation model for fiber-reinforced ceramic matrix composites, based on the damage state of the fiber-reinforced ceramic matrix composites under the heat-mechanics-oxygen coupling, and according to the interface phase oxidation model, a composite unit cell model is established to determine the time- and temperature-dependent internal fiber and matrix microscopic stress field of the composite material, thereby using the fracture mechanics interface debonding criterion based on the time- and temperature-dependent internal fiber and matrix microscopic stress field of the composite material under the heat-mechanics-oxygen coupling to determine the temperature-, stress- and time-dependent interface debonding length, obtain the temperature-, stress- and time-dependent matrix strain energy, use the initial matrix cracking stress and crack spacing to determine the critical matrix strain energy, and finally, use the critical matrix strain energy criterion based on the critical matrix strain energy. The present application is based on the mechanism of thermal-mechanical-oxygen coupling and the time- and temperature-dependent microscopic stress field of the fibers and matrix inside the composite material under the thermal-mechanical-oxygen coupling, so as to predict the matrix crack density of the fiber-reinforced ceramic matrix composite material under the thermal-mechanical-oxygen coupling. In addition, the present application also involves a composite unit cell model established based on the damage state of the fiber-reinforced ceramic matrix composite material under the thermal-mechanical-oxygen coupling, a fracture force interface debonding criterion, and a critical matrix strain energy criterion, etc., so as to more accurately predict the matrix crack density of the fiber-reinforced ceramic matrix composite material under the thermal-mechanical-oxygen coupling. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0024] Figure 1It is a schematic flow chart of a method for predicting matrix crack density in one embodiment of the present application.

[0025] Figure 2 A schematic diagram of a crack density curve of a SiC / SiC composite material matrix under the action of heat-mechanics-oxygen coupling provided in one embodiment of the present application. DETAILED DESCRIPTION

[0026] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0027] In order to make the above-mentioned objects, features and advantages of the present application more obvious and easy to understand, the present application is further described in detail below with reference to the accompanying drawings and specific implementation methods.

[0028] The present application embodiment provides a matrix crack density prediction method, which is executed by a computer device, and can be executed by a computer device such as a terminal or a server alone, or by a terminal and a server together. In the present application embodiment, Figure 1 As shown, the method includes the following steps.

[0029] S1: Based on the thermal-mechanical-oxygen coupling mechanism, a time- and temperature-dependent interfacial phase oxidation model for fiber-reinforced ceramic matrix composites is established.

[0030] S2: Based on the damage state of fiber-reinforced ceramic matrix composites under the coupling of heat, force and oxygen, a composite unit cell model is established according to the interface phase oxidation model, and the time and temperature dependent microscopic stress field of the internal fibers and the matrix of the composite material is determined according to the composite unit cell model.

[0031] S3: Based on the microscopic stress field between the fibers and the matrix inside the composite material, the interface debonding criterion of fracture mechanics is adopted to determine the temperature, stress and time dependent interface debonding length.

[0032] S4: Based on the microscopic stress field between the fibers and the matrix inside the composite material and the interface debonding length, the temperature, stress and time-dependent matrix strain energy is determined, and based on the temperature, stress and time-dependent matrix strain energy, the critical matrix strain energy is determined using the initial matrix cracking stress and crack spacing.

[0033] S5: Using a critical matrix strain energy criterion, the temperature, stress and time dependent matrix crack density is predicted based on the critical matrix strain energy.

[0034] In an exemplary embodiment, S1 specifically includes: under the action of heat-mechanics-oxygen coupling, the carbon interface phase on the fiber surface of the fiber-reinforced ceramic matrix composite material reacts with oxygen, so that the carbon interface phase on the fiber surface is oxidized, and a time and temperature dependent interface phase oxidation model of the fiber-reinforced ceramic matrix composite material is established.

[0035] The interface phase oxidation model is: ;in, l o is the oxidation length, the oxidation length l o is a function of temperature and time; α and β is the temperature-dependent oxidation parameter, , , T is the Kelvin temperature; γ is a constant, t is the oxidation time.

[0036] In an exemplary embodiment, the microscopic stress field between the fiber and the matrix inside the composite material is:

[0037] .

[0038] .

[0039] in, is the fiber stress, is stress, x is the length along the fiber axis, x=0 at the matrix crack; Ф f The fibers on the matrix crack plane bear the stress. r f is the fiber radius; is the shear stress in the interface oxidation zone; is the temperature and time dependent interfacial oxidation length; is the shear stress in the interface debonding zone; is the temperature and time dependent interfacial debonding length; is the fiber axial stress in the interface bonding zone; ρ is the shear lag model parameter; is the stress in the bonding area of ​​the matrix; is the volume content ratio of fiber to matrix; λ is the model parameter; is the axial stress of the matrix in the interface bonding area; is the temperature and time dependent matrix crack spacing.

[0040] In an exemplary embodiment, the fracture mechanics interface debonding criterion is: ; Among them, Г i is the interfacial debonding energy, df ( x =0, t , σ , T ) is the temperature, stress and time dependent fiber axial displacement in the matrix crack plane, v ( x,t , σ , T ) is the temperature, stress and time dependent fiber relative matrix displacement; is the interface debonding length.

[0041] .

[0042] .

[0043] in, E f , E m , E c are the Young's modulus of fiber, matrix and composite material, respectively. is the axial displacement of the fiber relative to the matrix; V m is the matrix volume fraction.

[0044] The interface debonding length satisfies: .

[0045] In an exemplary embodiment, the temperature, stress and time dependent matrix strain energy as follows.

[0046] .

[0047] in, is the substrate Young’s modulus; is the temperature, stress and time dependent interfacial debonding length; is the temperature, stress and time dependent matrix crack spacing.

[0048] The critical matrix strain energy for: ;in, is the critical matrix strain energy parameter; is the initial matrix cracking stress.

[0049] In an exemplary embodiment, the temperature, stress and time dependent matrix crack density ρ m for: .

[0050] Taking the SiC / SiC composite material under heat-mechanical-oxygen coupling as an example, the matrix crack density prediction method provided in the present application is used to predict the matrix crack density of the SiC / SiC composite material under heat-mechanical-oxygen coupling. As the temperature and stress change, the matrix crack density is as follows: Figure 2 shown.

[0051] Based on the same inventive concept, the embodiment of the present application also provides a substrate crack density prediction device for implementing the above-mentioned substrate crack density prediction method. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme recorded in the above-mentioned method, so the specific limitations of one or more substrate crack density prediction device embodiments provided below can refer to the limitations of the substrate crack density prediction method above, and will not be repeated here.

[0052] In an exemplary embodiment, a matrix crack density prediction device is provided, which includes the following modules.

[0053] The interface phase oxidation model building module is used to establish the time and temperature dependent interface phase oxidation model of fiber reinforced ceramic matrix composites based on the thermal-mechanical-oxygen coupling mechanism.

[0054] The module for determining the microscopic stress field of the fibers and the matrix inside the composite material is used to establish a composite material unit cell model based on the damage state of the fiber-reinforced ceramic matrix composite material under the action of heat-mechanics-oxygen coupling and the interface phase oxidation model, and to determine the time and temperature dependent microscopic stress field of the fibers and the matrix inside the composite material based on the composite material unit cell model.

[0055] The interface debonding length determination module is used to determine the temperature, stress and time-dependent interface debonding length based on the microscopic stress field between the fiber and the matrix inside the composite material and the fracture mechanics interface debonding criterion.

[0056] The critical matrix strain energy determination module is used to determine the temperature, stress and time-dependent matrix strain energy based on the microscopic stress field between the fibers and the matrix inside the composite material and the interface debonding length, and to determine the critical matrix strain energy based on the temperature, stress and time-dependent matrix strain energy using the initial matrix cracking stress and crack spacing.

[0057] The matrix crack density prediction module is used to predict the temperature, stress and time dependent matrix crack density based on the critical matrix strain energy using a critical matrix strain energy criterion.

[0058] In an exemplary embodiment, a computer device is provided, which may be a server or a terminal. The computer device includes a processor, a memory, an input / output interface (I / O for short) and a communication interface. The processor, the memory and the input / output interface are connected via a system bus, and the communication interface is connected to the system bus via the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store matrix crack density prediction data. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal via a network connection. When the computer program is executed by the processor, a matrix crack density prediction method is implemented.

[0059] In an exemplary embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and the above method is implemented when the processor executes the computer program.

[0060] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, which implements the above method when executed by a processor.

[0061] In an exemplary embodiment, a computer program product is provided, including a computer program, which implements the above method when executed by a processor.

[0062] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetoresistive random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM may be in various forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0063] In this application, all actions to obtain signals, information or data are carried out in compliance with the relevant data protection laws and policies of the country where they are located and with the authorization given by the owner of the corresponding device.

[0064] The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. The non-relational database may include a distributed database based on blockchain, etc., but is not limited thereto. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., but is not limited thereto.

[0065] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0066] This article uses specific examples to illustrate the principles and implementation methods of this application. The description of the above embodiments is only used to help understand the method and core ideas of this application. At the same time, for those skilled in the art, according to the ideas of this application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting this application.

Claims

1. A matrix crack density prediction method, characterized in that: The matrix crack density prediction method comprises: Based on the thermal-mechanical-oxygen coupling mechanism, a time- and temperature-dependent interfacial phase oxidation model of fiber-reinforced ceramic matrix composites is established, including: Under the action of heat-mechanics-oxygen coupling, the carbon interface phase on the fiber surface of the fiber-reinforced ceramic matrix composite reacts with oxygen, causing the carbon interface phase on the fiber surface to oxidize, and a time- and temperature-dependent interface phase oxidation model of the fiber-reinforced ceramic matrix composite is established; The interface phase oxidation model is: ;in, l o is the oxidation length; α and β is the temperature-dependent oxidation parameter, , , T is the Kelvin temperature; γ is a constant, t is the oxidation time; Based on the damage state of fiber-reinforced ceramic matrix composites under the action of heat-mechanics-oxygen coupling, a composite unit cell model is established according to the interface phase oxidation model, and the time- and temperature-dependent internal fiber and matrix microscopic stress field of the composite material is determined according to the composite unit cell model; Based on the microscopic stress field between the fiber and the matrix inside the composite material, the interface debonding criterion of fracture mechanics is adopted to determine the temperature, stress and time dependent interface debonding length; The fracture mechanics interface debonding criterion is: ； Among them, i is the interfacial debonding energy, d f ( x =0, t , σ , T ) is the temperature, stress and time dependent fiber axial displacement in the matrix crack plane, v ( x,t , σ , T ) is the temperature, stress and time dependent fiber relative matrix displacement; is the interface debonding length; r f is the fiber radius; Ф f The fibers in the matrix crack plane bear the stress; is the shear stress in the interface debonding zone; is the stress, x is the length along the fiber axis; Determining temperature, stress and time-dependent matrix strain energy based on the microscopic stress field between the fibers and the matrix inside the composite material and the interface debonding length, and determining critical matrix strain energy based on the temperature, stress and time-dependent matrix strain energy using initial matrix cracking stress and crack spacing; A critical matrix strain energy criterion is adopted to predict the temperature, stress and time dependent matrix crack density based on the critical matrix strain energy.

2. The matrix crack density prediction method according to claim 1, characterized in that: The microscopic stress field of the fiber and matrix inside the composite material is: ； ； in, is the fiber stress, is stress, x is the length along the fiber axis, x=0 at the matrix crack; Ф f The fibers on the matrix crack plane bear the stress. r f is the fiber radius; is the shear stress in the interface oxidation zone; is the temperature and time dependent interfacial oxidation length; is the shear stress in the interface debonding zone; is the temperature and time dependent interfacial debonding length; is the fiber axial stress in the interface bonding zone; ρ is the shear lag model parameter; is the stress in the bonding area of ​​the matrix; is the volume content ratio of fiber to matrix; λ is the model parameter; is the axial stress of the matrix in the interface bonding area; is the temperature and time dependent matrix crack spacing.

3. The matrix crack density prediction method according to claim 2, characterized in that: The temperature, stress and time dependent matrix strain energy for: ； in, is the substrate Young’s modulus; is the temperature, stress and time dependent interfacial debonding length; is the temperature, stress and time dependent matrix crack spacing; The critical matrix strain energy for: ； in, is the critical matrix strain energy parameter; is the initial matrix cracking stress.

4. The matrix crack density prediction method according to claim 3, characterized in that: The temperature, stress and time dependence of matrix crack density ρ m for: .

5. A matrix crack density prediction device using the matrix crack density prediction method according to any one of claims 1 to 4, characterized in that: The matrix crack density prediction device comprises: The interface phase oxidation model building module is used to establish the time-temperature dependent interface phase oxidation model of fiber reinforced ceramic matrix composites based on the thermal-mechanical-oxygen coupling mechanism; A module for determining the microscopic stress field of the fiber and matrix inside the composite material, which is used to establish a composite material unit cell model based on the damage state of the fiber-reinforced ceramic matrix composite material under the action of heat-mechanics-oxygen coupling and the interface phase oxidation model, and to determine the time- and temperature-dependent microscopic stress field of the fiber and matrix inside the composite material based on the composite material unit cell model; An interface debonding length determination module is used to determine the temperature, stress and time-dependent interface debonding length based on the microscopic stress field between the fiber and the matrix inside the composite material and the fracture mechanics interface debonding criterion; A critical matrix strain energy determination module is used to determine the temperature, stress and time-dependent matrix strain energy based on the microscopic stress field between the fiber and the matrix inside the composite material and the interface debonding length, and to determine the critical matrix strain energy based on the temperature, stress and time-dependent matrix strain energy by using the initial matrix cracking stress and crack spacing; The matrix crack density prediction module is used to predict the temperature, stress and time dependent matrix crack density based on the critical matrix strain energy using a critical matrix strain energy criterion.

6. A computer device comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the matrix crack density prediction method according to any one of claims 1 to 4.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the matrix crack density prediction method described in any one of claims 1 to 4 is implemented.

8. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the matrix crack density prediction method described in any one of claims 1 to 4 is implemented.

Images