Simulation analysis method for force and heat intensity of porous thermal protection structure in high-speed flow field environment

Abstract

The invention provides a force and heat intensity simulation analysis method for a porous thermal protection structure in a high-speed flow field environment. The method comprises the following steps: 1, carrying out aerodynamic thermal simulation calculation by adopting a CFD (Computational Fluid Dynamics) method to obtain cold wall heat flux density and recovery enthalpy of the porous thermal protection structure at a typical moment; 2, performing transient heat transfer CAE simulation analysis on the porous thermal protection structure according to a result obtained in the step 1 to obtain the temperature of the porous thermal protection structure; 3, the distributed aerodynamic force of the porous thermal protection structure is obtained; 4, respectively carrying out stretching and shearing tests on the porous thermal protection material at different temperatures to obtain the minimum values of rigidity in the stretching and shearing directions and the minimum value of nominal stress of initial damage at different temperatures; and 5, completing the force and heat simulation analysis of the porous thermal protection structure based on parametric modeling through Abaqus secondary development. The method solves the problem that in the prior art, local block-shaped stripping of the porous thermal protection structure with random defects in the high-speed flow field environment is difficult to predict.

Description

Technical Field

[0001] The present invention belongs to the technical field of structural thermal strength analysis, and in particular relates to a method for simulating and analyzing the thermal strength of a porous thermal protection structure in a high-speed flow field environment. Background Art

[0002] Structural thermal strength analysis is an important method for the design and verification of aircraft structures and thermal protection systems. As the flight speed of modern aircraft continues to increase, the thermal loads faced by thermal protection structures extending out of the aircraft surface are becoming increasingly severe due to the strong compression of the airflow behind the shock wave and the friction between the airflow and the structure in the boundary layer. For porous thermal protection structures extending out of the aircraft surface, due to the existence of random pore defects, existing conventional thermal simulation analysis technologies often find it difficult to predict the local block peeling of thermal protection structures in high-speed flow fields, and thus cannot meet the needs of thermal strength analysis of porous thermal protection structures under high-speed flow fields. Summary of the invention

[0003] The present invention aims to solve at least one of the technical problems existing in the prior art.

[0004] To this end, the present invention provides a method for simulating and analyzing the mechanical and thermal strength of porous thermal protection structures under a high-speed flow field environment.

[0005] The technical solution of the present invention is as follows:

[0006] According to one aspect, a method for simulating and analyzing the mechanical and thermal strength of a porous thermal protection structure in a high-speed flow field environment is provided, the method comprising:

[0007] Step 1: Use CFD method to perform aerodynamic thermal simulation calculation to obtain the cold wall heat flux density and recovery enthalpy of the porous thermal protection structure at typical moments;

[0008] Step 2: Perform transient heat transfer CAE simulation analysis of the porous thermal protection structure according to the result obtained in step 1 to obtain the temperature of the porous thermal protection structure;

[0009] Step 3: Use CFD method to perform aerodynamic simulation calculation to obtain the distributed aerodynamic force of the porous thermal protection structure;

[0010] Step 4: Conduct tensile and shear tests on porous thermal protection materials at different temperatures to obtain the minimum values ​​of tensile and shear stiffness and the minimum value of nominal stress at the beginning of damage at different temperatures;

[0011] Step 5: Complete the mechanical and thermal simulation analysis of porous thermal protection structures based on parametric modeling through Abaqus secondary development, including:

[0012] 5.1) Prefabricate the peeling path and establish a parameterized model;

[0013] 5.2) Set up constitutive model and units;

[0014] 5.3) Based on steps 5.1) and 5.2), the distributed aerodynamic force calculated in step 3 and the temperature field of the porous thermal protection structure calculated in step 2 are used as the boundary conditions of the mechanical and thermal loads to carry out mechanical and thermal coupling analysis;

[0015] 5.4) Based on step 5.3), parameter optimization is performed to obtain the most favorable and most unfavorable positions for peeling of the porous thermal protection structure.

[0016] Furthermore, the step 2 specifically includes:

[0017] 2.1) Construct a finite element model of a porous thermal protection structure; wherein a three-dimensional rectangular coordinate system is established, the origin O is the vertex of the edge of the windward surface of the thermal protection structure, the X-axis is in the same direction as the airflow direction, the Y-axis is vertically upward, and the Z-axis is perpendicular to the OXY plane, and the direction is determined by the right-hand rule;

[0018] 2.2) Taking the cold wall heat flux density and recovery enthalpy calculated at the typical time in step 1 as input, the cold wall heat flux density and recovery enthalpy of each unit integral point on the surface of the porous thermal protection structure at any calculation time are obtained through an interpolation algorithm;

[0019] 2.3) In combination with the flow field parameters, the cold wall heat flux density and recovery enthalpy of each unit integral point on the surface of the porous thermal protection structure obtained in step 2.2) at the calculation time are used as input to convert the cold and hot wall heat flux densities, and obtain the hot wall heat flux density loaded on the surface of the porous thermal protection structure at any calculation time;

[0020] 2.4) Using the heat flux density of the hot wall at any calculation time obtained in step 2.3) as input, a transient heat transfer CAE simulation analysis of the porous thermal protection structure is performed to obtain the temperature of the porous thermal protection structure.

[0021] Further, step 5.1) specifically includes:

[0022] Two parallel reference surfaces with a fixed distance are used as two peeling surfaces, which are counted as A and B surfaces respectively. The peeling surface A takes the XOZ surface in the coordinate system of step 2.1) as the starting plane, and takes any point M on the OY axis with a distance d from the origin O as the reference, ensuring that the peeling surface A always passes through point M, and then takes M as the rotation center, rotates around the Z axis with an angle of θ1, then rotates around the Y axis with an angle of θ2, and then rotates around the X axis with an angle of θ3, so as to control the parameters d, θ1, θ2, and θ3 of the A surface to achieve arbitrary angle cutting of the porous thermal protection structure model, and define the area between the A and B surfaces after cutting as the peeling area.

[0023] Further, step 5.2) specifically includes:

[0024] The bilinear cohesion model is used to design the peeling area, and the minimum stiffness values ​​in the tensile and shear directions at different temperatures are input for the elastic segment. The damage initiation is judged based on the quadratic nominal stress criterion, and the minimum nominal stress values ​​at which damage starts in the tensile and shear directions at different temperatures are input. The cohesive unit is selected as the unit type.

[0025] The remaining area except the peeling area is designed as isotropic material, and the three-dimensional stress element is selected as the element type.

[0026] Further, step 5.4) specifically includes:

[0027] Taking the maximum or minimum of the damage variable QUADSCRT of each working condition as the goal, and the corresponding d, θ1, θ2, θ3 parameter combination as the variable, the global or local algorithm is used to optimize and obtain the most favorable and unfavorable positions for peeling of the porous thermal protection structure.

[0028] Furthermore, the algorithm in step 5.4) includes a genetic algorithm and a sequential quadratic programming.

[0029] Furthermore, the step 4 specifically includes:

[0030] 4.1) Conducting tensile and shear tests on the porous thermal protection structure at different typical temperatures, where the typical temperatures can cover the temperature range of the porous thermal protection structure calculated in step 2, and obtaining tensile and shear force-displacement curves at different typical temperatures;

[0031] 4.2) The minimum values ​​of stiffness in the tensile and shear directions and the minimum values ​​of nominal stress at the onset of damage for all samples at different temperatures were counted respectively.

[0032] Furthermore, in step 4.1), there are no less than 5 valid samples at each temperature.

[0033] According to another aspect, a computer-readable storage medium is provided, storing a computer program, wherein the computer program implements the simulation analysis method as described above when executed by a processor.

[0034] According to yet another aspect, a computer program product is provided, comprising computer instructions, wherein when the computer instructions are executed by a processor, the simulation analysis method as described above is implemented.

[0035] Compared with the prior art, the present invention provides a method for simulating and analyzing the mechanical and thermal strength of porous thermal protection structures under a high-speed flow field environment, which can effectively predict the situation of local block peeling, and solves the problem that conventional mechanical and thermal simulation analysis technology is difficult to predict the local block peeling of porous thermal protection structures containing random defects under a high-speed flow field environment. The prefabricated peeling path and parameterized search method proposed in the present invention can quickly search out the peeling surface of the porous thermal protection structure under mechanical and thermal loads, and can provide a solution to related similar problems. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] The included drawings are used to provide a further understanding of the embodiments of the present invention, which constitute a part of the specification, are used to illustrate the embodiments of the present invention, and together with the text description, explain the principles of the present invention. Obviously, the drawings in the following description are only some embodiments of the present invention, and for ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0037] Figure 1 A schematic flow chart of a method for simulating and analyzing the mechanical and thermal strength of a porous thermal protection structure in a high-speed flow field environment is shown according to a specific embodiment of the present invention;

[0038] Figure 2 A schematic diagram of a simplified model of parameterized modeling provided according to a specific embodiment of the present invention is shown. DETAILED DESCRIPTION

[0039] It should be noted that, in the absence of conflict, the embodiments in this application and the features in the embodiments can be combined with each other. The technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. The following description of at least one exemplary embodiment is actually only illustrative and is by no means intended to limit the present invention and its application or use. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0040] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit the exemplary embodiments according to the present application. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprise" and / or "include" are used in this specification, it indicates the presence of features, steps, operations, devices, components and / or combinations thereof.

[0041] Unless otherwise specifically stated, the relative arrangement of the parts and steps described in these embodiments, numerical expressions and numerical values ​​do not limit the scope of the present invention. At the same time, it should be understood that, for ease of description, the sizes of the various parts shown in the accompanying drawings are not drawn according to the actual proportional relationship. The technology, method and equipment known to ordinary technicians in the relevant field may not be discussed in detail, but in appropriate cases, the technology, method and equipment should be regarded as a part of the authorization specification. In all examples shown and discussed here, any specific value should be interpreted as being merely exemplary, rather than as a limitation. Therefore, other examples of exemplary embodiments may have different values. It should be noted that similar reference numerals and letters represent similar items in the following drawings, so once a certain item is defined in an accompanying drawing, it does not need to be further discussed in subsequent drawings.

[0042] like Figure 1 As shown, in one embodiment of the present invention, a method for simulating and analyzing the mechanical and thermal strength of a porous thermal protection structure in a high-speed flow field environment is provided, the method comprising:

[0043] Step 1: Use CFD method to perform aerodynamic thermal simulation calculation to obtain the cold wall heat flux density and recovery enthalpy of the porous thermal protection structure at typical moments;

[0044] Step 2: Perform transient heat transfer CAE simulation analysis of the porous thermal protection structure according to the result obtained in step 1 to obtain the temperature of the porous thermal protection structure;

[0045] Step 3: Use CFD method to perform aerodynamic simulation calculation to obtain the distributed aerodynamic force of the porous thermal protection structure;

[0046] Step 4: Conduct tensile and shear tests on porous thermal protection materials at different temperatures to obtain the minimum values ​​of tensile and shear stiffness and the minimum value of nominal stress at the beginning of damage at different temperatures;

[0047] Step 5: Complete the mechanical and thermal simulation analysis of porous thermal protection structures based on parametric modeling through Abaqus secondary development, including:

[0048] 5.1) Prefabricate the peeling path and establish a parameterized model;

[0049] 5.2) Set up constitutive model and units;

[0050] 5.3) Based on steps 5.1) and 5.2), the distributed aerodynamic force calculated in step 3 and the temperature field of the porous thermal protection structure calculated in step 2 are used as the boundary conditions of the mechanical and thermal loads to carry out mechanical and thermal coupling analysis;

[0051] 5.4) Based on step 5.3), parameter optimization is performed to obtain the most favorable and most unfavorable positions for peeling of the porous thermal protection structure.

[0052] In the embodiment of the present invention, the step 2 specifically includes:

[0053] 2.1) Construct a finite element model of a porous thermal protection structure; wherein a three-dimensional rectangular coordinate system is established, the origin O is the vertex of the edge of the windward surface of the thermal protection structure, the X-axis is in the same direction as the airflow direction, the Y-axis is vertically upward, and the Z-axis is perpendicular to the OXY plane, and the direction is determined by the right-hand rule;

[0054] 2.2) Taking the cold wall heat flux density and recovery enthalpy calculated at the typical time in step 1 as input, the cold wall heat flux density and recovery enthalpy of each unit integral point on the surface of the porous thermal protection structure at any calculation time are obtained through an interpolation algorithm;

[0055] 2.3) In combination with the flow field parameters, the cold wall heat flux density and recovery enthalpy of each unit integral point on the surface of the porous thermal protection structure obtained in step 2.2) at the calculation time are used as input to convert the cold and hot wall heat flux densities, and obtain the hot wall heat flux density loaded on the surface of the porous thermal protection structure at any calculation time;

[0056] 2.4) Using the heat flux density of the hot wall at any calculation time obtained in step 2.3) as input, a transient heat transfer CAE simulation analysis of the porous thermal protection structure is performed to obtain the temperature of the porous thermal protection structure.

[0057] In the embodiment of the present invention, step 5.1) specifically includes:

[0058] Two parallel reference surfaces with a fixed distance are used as two peeling surfaces, which are counted as A and B surfaces respectively. The peeling surface A takes the XOZ surface in the coordinate system of step 2.1) as the starting plane, and takes any point M on the OY axis with a distance d from the origin O as the reference, ensuring that the peeling surface A always passes through point M, and then takes M as the rotation center, rotates around the Z axis with an angle of θ1, then rotates around the Y axis with an angle of θ2, and then rotates around the X axis with an angle of θ3, so as to control the parameters d, θ1, θ2, and θ3 of the A surface to achieve arbitrary angle cutting of the porous thermal protection structure model, and define the area between the A and B surfaces after cutting as the peeling area.

[0059] In the embodiment of the present invention, step 5.2) specifically includes:

[0060] The bilinear cohesion model is used to design the peeling area, and the minimum stiffness values ​​in the tensile and shear directions at different temperatures are input for the elastic segment. The damage initiation is judged based on the quadratic nominal stress criterion, and the minimum nominal stress values ​​at which damage starts in the tensile and shear directions at different temperatures are input. The cohesive unit is selected as the unit type.

[0061] The remaining area except the peeling area is designed as isotropic material, and the three-dimensional stress element is selected as the element type.

[0062] In the embodiment of the present invention, step 5.4) specifically includes:

[0063] Taking the maximum or minimum of the damage variable QUADSCRT of each working condition as the goal, and the corresponding d, θ1, θ2, θ3 parameter combination as the variable, the global or local algorithm is used to optimize and obtain the most favorable and unfavorable positions for peeling of the porous thermal protection structure.

[0064] It can be seen that the existing conventional mechanical and thermal simulation analysis technology is difficult to predict the problem of local block peeling of porous thermal protection structures with random defects in high-speed flow field environments. The embodiment of the present invention proposes a prefabricated peeling path, using the idea of ​​the measured minimum material parameters in the peeling area; and proposes a parameterized search method to quickly search for the peeling surface of the porous thermal protection structure under mechanical and thermal loads, which can provide a solution to related similar problems. Compared with the existing technology, the method of the present invention can effectively predict the situation of local block peeling.

[0065] According to an optional embodiment of the present invention, the algorithm in step 5.4) includes a genetic algorithm and a sequential quadratic programming.

[0066] In the above embodiment, the step 4 specifically includes:

[0067] 4.1) Conducting tensile and shear tests on the porous thermal protection structure at different typical temperatures, where the typical temperatures can cover the temperature range of the porous thermal protection structure calculated in step 2, and obtaining tensile and shear force-displacement curves at different typical temperatures;

[0068] 4.2) The minimum values ​​of stiffness in the tensile and shear directions and the minimum values ​​of nominal stress at the onset of damage for all samples at different temperatures were counted respectively.

[0069] Preferably, in step 4.1), there are no less than 5 valid samples at each temperature.

[0070] According to another embodiment, a computer-readable storage medium is provided, storing a computer program, wherein the computer program implements the simulation analysis method as described above when executed by a processor.

[0071] According to another embodiment, a computer program product is provided, including computer instructions, which implement the simulation analysis method described above when executed by a processor.

[0072] In order to further understand the simulation analysis method of the embodiment of the present invention, a specific embodiment is described in detail below:

[0073] like Figure 1 to Figure 2 As shown, according to a specific embodiment of the present invention, a method for simulating and analyzing the mechanical and thermal strength of a porous thermal protection structure in a high-speed flow field environment is provided, and the method includes the following steps.

[0074] 1) The CFD method is used to perform aerodynamic thermal simulation calculations to obtain the cold wall heat flux density and recovery enthalpy of the porous thermal protection structure at typical moments.

[0075] 2) According to the result of the aerodynamic thermal simulation calculation in step 1), a transient heat transfer CAE simulation analysis of the porous thermal protection structure is performed to obtain the temperature of the porous thermal protection structure:

[0076] 2.1) Construct a finite element model of a porous thermal protection structure, such as Figure 2 As shown. A three-dimensional rectangular coordinate system is established, with the origin O being the vertex of the edge of the windward surface of the thermal protection structure, the X axis being in the same direction as the airflow direction, the Y axis being vertically upward, and the Z axis being perpendicular to the OXY plane, and the direction being determined by the right-hand rule;

[0077] 2.2) Using the typical cold wall heat flux density and recovery enthalpy obtained in step 1) as input, the cold wall heat flux density and recovery enthalpy of each unit integral point on the surface of the porous thermal protection structure at any calculation time are obtained through an interpolation algorithm;

[0078] 2.3) In combination with the flow field parameters, the cold wall heat flux density and recovery enthalpy of each unit integral point on the surface of the porous thermal protection structure obtained in step 2.2) at the calculation time are used as input, and the hot wall heat flux density loaded on the surface of the porous thermal protection structure at any calculation time is obtained through the cold and hot wall heat flux density conversion formula;

[0079] 2.4) Using the heat flux density of the hot wall at any calculation time obtained in step 2.3) as input, a transient heat transfer CAE simulation analysis of the porous thermal protection structure is performed to obtain the temperature of the porous thermal protection structure.

[0080] 3) Use CFD method to perform aerodynamic simulation calculation to obtain the distributed aerodynamic force of porous thermal protection structure;

[0081] 4) Carry out tensile and shear tests on porous thermal protection materials at different temperatures to obtain the minimum values ​​of tensile and shear stiffness and the minimum value of nominal stress at the beginning of damage at different temperatures:

[0082] 4.1) Conduct tensile and shear tests on the porous thermal protection structure at different typical temperatures, where the typical temperatures can cover the temperature range of the porous thermal protection structure calculated in step 2), and at least 5 valid samples are used at each temperature to obtain tensile and shear force-displacement curves at different typical temperatures;

[0083] 4.2) Count the lowest values ​​of tensile and shear stiffness and the lowest nominal stress at the beginning of damage for all samples at different temperatures;

[0084] 5) Complete the mechanical and thermal simulation analysis of porous thermal protection structures based on parametric modeling through Abaqus secondary development:

[0085] 5.1) Prefabricate the peeling path and establish a parameterized model: Figure 2 , using two parallel reference surfaces with a fixed distance as two peeling surfaces, counted as A and B surfaces respectively, where peeling surface A takes the XOZ surface as the starting plane, and takes any point M on the OY axis with a distance d from the origin O as the reference, ensuring that peeling surface A always passes through point M, and then takes M as the rotation center, rotates around the Z axis at an angle of θ1, then around the Y axis at an angle of θ2, and then around the X axis at an angle of θ3. By controlling the parameters d, θ1, θ2, and θ3 of surface A, the porous thermal protection structure model can be cut at any angle, and the area between surfaces A and B after cutting is defined as the peeling area;

[0086] 5.2) Set constitutive model and unit: The bilinear cohesion model is used in the peeling area, and the minimum stiffness values ​​in the tensile and shear directions at different temperatures are input for the elastic segment. The damage initiation is judged based on the quadratic nominal stress criterion, and the minimum nominal stress value at which the damage starts in the tensile and shear directions at different temperatures is input. The damaged segment is ignored, and the unit type is cohesive unit. The remaining area except the peeling area is regarded as isotropic material, and conventional engineering material parameters are used. The unit type is selected as three-dimensional stress unit.

[0087] 5.3) Use the distributed aerodynamic force calculated in step 3) and the temperature field of the porous thermal protection structure calculated in step 2) as the boundary conditions of the mechanical and thermal loads to carry out mechanical and thermal coupling analysis;

[0088] 5.4) Parameter optimization: Taking the maximum / minimum of the damage variable QUADSCRT of each working condition as the goal, and the corresponding d, θ1, θ2, θ3 parameter combination as the variable, the variation range of d, θ1, θ2, θ3 parameters is given respectively, and the global / local algorithms such as genetic algorithm and sequential quadratic programming are used to optimize, so as to obtain the most favorable and most unfavorable positions for the peeling of the porous thermal protection structure, thereby effectively guiding the structural design.

[0089] Features described and / or illustrated above for one embodiment may be used in the same or similar manner in one or more other embodiments, and / or combined with or used in place of features in other embodiments.

[0090] It should be emphasized that the term "include / comprises" when used herein refers to the presence of features, integers, steps or components, but does not exclude the presence or addition of one or more other features, integers, steps, components or combinations thereof.

[0091] The above method of the present invention can be implemented by hardware, or by hardware combined with software. The present invention relates to such a computer-readable program, which, when executed by a logic component, enables the logic component to implement the above-mentioned device or component, or enables the logic component to implement the above-mentioned various methods or steps. The present invention also relates to a storage medium for storing the above program, such as a hard disk, a magnetic disk, an optical disk, a DVD, a flash memory, etc.

[0092] The many features and advantages of these embodiments are apparent from this detailed description, and thus the appended claims are intended to cover all such features and advantages of these embodiments that fall within their true spirit and scope. Furthermore, since numerous modifications and changes will readily occur to those skilled in the art, it is not intended that the embodiments of the invention be limited to the exact construction and operation illustrated and described, but rather all suitable modifications and equivalents falling within the scope thereof are intended to be covered.

[0093] Parts of the present invention that are not described in detail are well known to those skilled in the art.

Claims

1. A method for simulating and analyzing the mechanical and thermal strength of porous thermal protection structures under high-speed flow field environment, characterized in that: The simulation analysis method comprises: Step 1: Use CFD method to perform aerodynamic thermal simulation calculation to obtain the cold wall heat flux density and recovery enthalpy of the porous thermal protection structure at typical moments; Step 2: Perform transient heat transfer CAE simulation analysis of the porous thermal protection structure according to the result obtained in step 1 to obtain the temperature of the porous thermal protection structure; Step 3: Use CFD method to perform aerodynamic simulation calculation to obtain the distributed aerodynamic force of the porous thermal protection structure; Step 4: Conduct tensile and shear tests on porous thermal protection materials at different temperatures to obtain the minimum values ​​of tensile and shear stiffness and the minimum value of nominal stress at the beginning of damage at different temperatures; Step 5: Complete the mechanical and thermal simulation analysis of porous thermal protection structures based on parametric modeling through Abaqus secondary development, including: 5.1) Prefabricate the peeling path and establish a parameterized model; 5.2) Set up constitutive model and units; 5.3) Based on steps 5.1) and 5.2), the distributed aerodynamic force calculated in step 3 and the temperature field of the porous thermal protection structure calculated in step 2 are used as the boundary conditions of the mechanical and thermal loads to carry out mechanical and thermal coupling analysis; 5.4) Based on step 5.3), parameter optimization is performed to obtain the most favorable and most unfavorable positions for peeling of the porous thermal protection structure.

2. The method for simulating and analyzing the mechanical and thermal strength of porous thermal protection structures under a high-speed flow field environment according to claim 1 is characterized in that: The step 2 specifically includes: 2.1) Construct a finite element model of a porous thermal protection structure; wherein a three-dimensional rectangular coordinate system is established, the origin O is the vertex of the edge of the windward surface of the thermal protection structure, the X-axis is in the same direction as the airflow direction, the Y-axis is vertically upward, and the Z-axis is perpendicular to the OXY plane, and the direction is determined by the right-hand rule; 2.2) Taking the cold wall heat flux density and recovery enthalpy calculated at the typical time in step 1 as input, the cold wall heat flux density and recovery enthalpy of each unit integral point on the surface of the porous thermal protection structure at any calculation time are obtained through an interpolation algorithm; 2.3) In combination with the flow field parameters, the cold wall heat flux density and recovery enthalpy of each unit integral point on the surface of the porous thermal protection structure obtained in step 2.2) at the calculation time are used as input to convert the cold and hot wall heat flux densities, and obtain the hot wall heat flux density loaded on the surface of the porous thermal protection structure at any calculation time; 2.4) Using the hot wall heat flux density at any calculation time obtained in step 2.3) as input, a transient heat transfer CAE simulation analysis of the porous thermal protection structure is performed to obtain the temperature of the porous thermal protection structure.

3. The method for simulating and analyzing the mechanical and thermal strength of porous thermal protection structures under a high-speed flow field environment according to claim 2 is characterized in that: Step 5.1) specifically includes: Two parallel reference surfaces with a fixed distance are used as two peeling surfaces, which are counted as A and B surfaces respectively. The peeling surface A takes the XOZ surface in the coordinate system of step 2.1) as the starting plane, and takes any point M on the OY axis with a distance d from the origin O as the reference, ensuring that the peeling surface A always passes through point M, and then takes M as the rotation center, rotates around the Z axis with an angle of θ1, then rotates around the Y axis with an angle of θ2, and then rotates around the X axis with an angle of θ3, so as to control the parameters d, θ1, θ2, and θ3 of the A surface to achieve arbitrary angle cutting of the porous thermal protection structure model, and define the area between the A and B surfaces after cutting as the peeling area.

4. The method for simulating and analyzing the mechanical and thermal strength of porous thermal protection structures under a high-speed flow field environment according to claim 3 is characterized in that: Step 5.2) specifically includes: The bilinear cohesion model is used to design the peeling area, and the minimum stiffness values ​​in the tensile and shear directions at different temperatures are input for the elastic segment. The damage initiation is judged based on the quadratic nominal stress criterion, and the minimum nominal stress values ​​at which damage starts in the tensile and shear directions at different temperatures are input. The cohesive unit is selected as the unit type. The remaining area except the peeling area is designed as isotropic material, and the three-dimensional stress element is selected as the element type.

5. The method for simulating and analyzing the mechanical and thermal strength of porous thermal protection structures under a high-speed flow field environment according to claim 4 is characterized in that: Step 5.4) specifically includes: Taking the maximum or minimum of the damage variable QUADSCRT of each working condition as the goal, and the corresponding d, θ1, θ2, θ3 parameter combination as the variable, the global or local algorithm is used to optimize and obtain the most favorable and unfavorable positions for peeling of the porous thermal protection structure.

6. A method for simulating and analyzing the mechanical and thermal strength of porous thermal protection structures under a high-speed flow field environment according to any one of claims 1 to 5, characterized in that: The step 4 specifically includes: 4.1) Conducting tensile and shear tests on the porous thermal protection structure at different typical temperatures, where the typical temperatures can cover the temperature range of the porous thermal protection structure calculated in step 2, and obtaining tensile and shear force-displacement curves at different typical temperatures; 4.2) The minimum values ​​of stiffness in the tensile and shear directions and the minimum values ​​of nominal stress at the onset of damage for all samples at different temperatures were counted respectively.

7. The method for simulating and analyzing the mechanical and thermal strength of porous thermal protection structures under high-speed flow field environment according to claim 6 is characterized in that: In step 4.1), there are no less than 5 valid samples at each temperature.

8. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the simulation analysis method described in claims 1-8 is implemented.

9. A computer program product, characterized in that It includes computer instructions, and when the computer instructions are executed by a processor, the simulation analysis method described in claims 1-8 is implemented.

Images