Structural analysis method and device based on coupling microcrack closing effect damage model

Abstract

The invention provides a structural analysis method and device based on a coupling microcrack closing effect damage model, and relates to the technical field of rocket engines, and the method comprises the steps: constructing a damage constitutive model based on a microcrack closing effect; carrying out parameter identification on the damage constitutive model by utilizing a material mechanical property test result; based on the damage constitutive model after parameter identification, structural analysis is conducted on the thrust chamber, the method solves the technical problems that structural analysis of the thrust chamber is low in accuracy and efficiency, and the technical effects of improving the structural analysis efficiency and precision and improving the accuracy of thrust chamber service life estimation are achieved.

Description

Technical Field

[0001] This invention relates to the field of rocket engine technology, and in particular to a structural analysis method and apparatus based on a damage model of coupled microcrack closure effect. Background Technology

[0002] Reusability is a significant feature of next-generation launch vehicles, and advanced liquid propulsion is a crucial element in achieving this. During the multiple start-up and shutdown cycles of a reusable liquid rocket engine, the thrust chamber, as the core component, endures complex alternating thermomechanical loads. Its internal wall material exhibits cyclical characteristics such as kinematic hardening, isotropic hardening, or softening, accompanied by time-dependent inelastic deformation and stress relaxation. Furthermore, the material inevitably possesses micro-defects such as microcracks and micropores. These microcracks and micropores expand under tension and close under compression, resulting in a microcrack closure effect. This microcrack closure effect affects the structural response of the thrust chamber during subsequent operation, thus impacting the accuracy of lifespan prediction.

[0003] In existing technologies, the use of constitutive models for thrust chamber structural analysis often neglects the effect of microcrack closure, resulting in low accuracy of analysis results. Furthermore, since constitutive models contain numerous temperature-related material parameters, parameter identification through model simulation is time-consuming and computationally inefficient, which also significantly affects the efficiency of thrust chamber structural analysis and subsequent life prediction. Summary of the Invention

[0004] The purpose of this invention is to provide a structural analysis method and apparatus based on a coupled microcrack closure effect damage model, so as to solve the technical problems of low accuracy and low efficiency in the analysis of thrust chamber structures in the prior art.

[0005] In a first aspect, embodiments of the present invention provide a structural analysis method based on a damage model of coupled microcrack closure effect. The method includes: constructing a damage constitutive model based on the microcrack closure effect; identifying parameters of the damage constitutive model using material mechanical property test results; and performing structural analysis on the thrust chamber based on the damage constitutive model after parameter identification.

[0006] In some optional implementations, a damage constitutive model is constructed based on the microcrack closure effect, including: determining the crack closure effect coefficient according to the cross-sectional area parameter and tensile / compressive state of the material; determining the state potential function coupled with the microcrack closure effect based on the above crack closure effect coefficient and the basic state potential function; and modeling the material behavior based on the state variables, the above-mentioned state potential function coupled with the microcrack closure effect, and the dissipation potential function to generate a damage constitutive model.

[0007] In some optional implementations, the damage constitutive model described above is a viscoplastic damage constitutive model that considers the microcrack closure effect; the damage constitutive model described above includes initial independent variables; the initial independent variables described above include multiple initial material parameters; the material parameters described above include: damage-related parameters and non-damage-related parameters; the method described above also includes: performing numerical calculations on the damage constitutive model described above to determine the target independent variables in the current time step corresponding to the initial independent variables described above.

[0008] In some optional implementations, the damage constitutive model is parameterized using the results of material mechanical property tests, including: obtaining stress-variable force data of the material under different states through a type of material mechanical property test; determining the experimental values ​​of relevant material mechanical property parameters based on the stress-variable force data; identifying and optimizing the initial material parameters of the damage constitutive model using the experimental values ​​of the relevant parameters, and determining the optimized material parameters; and generating a viscoplastic damage constitutive model of the copper alloy based on the optimized material parameters.

[0009] In some optional implementations, the initial material parameters of the aforementioned damage constitutive model are identified and optimized using the experimental values ​​of the relevant parameters. This includes: using a random search algorithm to perform a global search on the simulation results of the aforementioned damage constitutive model to determine the initial parameter vector; and using a local optimization strategy to perform a fine search on the aforementioned initial parameter vector to determine the target material parameters.

[0010] In some optional implementations, the initial material parameters of the damage constitutive model are identified and optimized using the experimental values ​​of the aforementioned relevant parameters, including: simulating a single material element using the aforementioned damage constitutive model to determine the first initial material parameter; the aforementioned first initial material parameter is a non-damage-related parameter; simulating a full-size material sample model using the aforementioned damage constitutive model, and determining the damage-related parameters of the target material based on the aforementioned non-damage-related parameters and the experimental values ​​of the aforementioned relevant parameters.

[0011] In some optional implementations, the thrust chamber is structurally analyzed based on the aforementioned damage constitutive model after parameter identification, including: determining the temperature distribution of the inner wall of the thrust chamber through three-dimensional flow-heat coupling simulation; determining the target cross section based on the aforementioned temperature distribution results, and performing heat transfer analysis on the aforementioned target cross section; and combining the results of the aforementioned heat transfer analysis, performing structural analysis on the aforementioned thrust chamber using the aforementioned damage constitutive model after parameter identification.

[0012] Secondly, embodiments of the present invention provide a structural analysis device based on a damage model of coupled microcrack closure effect. The device includes: a model construction module for constructing a damage constitutive model based on the microcrack closure effect; a parameter identification module for identifying parameters of the damage constitutive model using material mechanical property test results; and a structural analysis module for performing structural analysis on the thrust chamber based on the damage constitutive model after parameter identification.

[0013] Thirdly, embodiments of the present invention provide an electronic device, including a memory and a processor, wherein the memory stores a computer program that can run on the processor, and the processor executes the computer program to implement the steps of the method described in any of the first aspects above.

[0014] Fourthly, embodiments of the present invention provide a computer-readable storage medium storing computer-executable instructions, which, when invoked and executed by a processor, cause the processor to perform the method described in any of the first aspects above.

[0015] This invention provides a structural analysis method and apparatus based on a damage model with coupled microcrack closure effect. The method includes: constructing a damage constitutive model based on the microcrack closure effect; identifying parameters of the damage constitutive model using material mechanical property test results; and performing structural analysis on the thrust chamber based on the parameter-identified damage constitutive model. This method solves the technical problems of low accuracy and low efficiency in thrust chamber structural analysis, and achieves the technical effects of improving the efficiency and accuracy of structural analysis and improving the accuracy of thrust chamber life prediction. Attached Figure Description

[0016] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0017] Figure 1 A schematic flowchart of a structural analysis method based on a damage model of coupled microcrack closure effect provided in an embodiment of the present invention; Figure 2 A schematic flowchart of another structural analysis method based on a coupled microcrack closure effect damage model provided in an embodiment of the present invention; Figure 3 A schematic diagram of a structural analysis device based on a damage model of coupled microcrack closure effect provided in an embodiment of the present invention; Figure 4 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0018] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. The components of the embodiments of the present invention described and shown in the accompanying drawings can generally be arranged and designed in various different configurations.

[0019] Therefore, the following detailed description of the embodiments of the invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely to illustrate selected embodiments of the invention. All other embodiments obtained by those skilled in the art based on the embodiments of the invention without inventive effort are within the scope of protection of the invention.

[0020] It should be noted that similar reference numerals and letters in the following figures indicate similar items; therefore, once an item is defined in one figure, it does not need to be further defined and explained in subsequent figures. Some embodiments of the present invention will now be described in detail with reference to the accompanying drawings. Unless otherwise specified, the following embodiments and features thereof can be combined with each other.

[0021] Under repeated loading conditions, such as multiple ignition and shutdown cycles of a rocket engine, microcracks may form and propagate on the inner wall of the thrust chamber. The microcrack closure effect can slow crack growth, thereby increasing the fatigue life of the thrust chamber. However, if microcracks propagate to a certain size, they may cause irreversible damage, affecting the safety and reliability of the thrust chamber. In other words, the microcrack closure effect affects the structural response of the thrust chamber during subsequent operation, thus impacting the accuracy of life prediction. Therefore, it is necessary to consider the microcrack closure effect and establish a constitutive model that describes the mechanical behavior of the inner wall material of the thrust chamber, laying the foundation for thrust chamber life prediction.

[0022] In existing applications of constitutive models for thrust chamber structural analysis, the effect of microcrack closure is often overlooked, leading to low accuracy of the analysis results. Furthermore, since constitutive models contain numerous temperature-related material parameters, these parameters need to be identified through material mechanical property tests and optimization algorithms. However, simulations of constitutive models coupled with material damage are time-consuming and computationally inefficient, significantly impacting the efficiency of thrust chamber structural analysis and subsequent life prediction.

[0023] Based on this, embodiments of the present invention provide a structural analysis method and apparatus based on a coupled microcrack closure effect damage model to solve the technical problems of low accuracy and low efficiency in thrust chamber structural analysis in the prior art.

[0024] To facilitate understanding of this embodiment, a structural analysis method based on a coupled microcrack closure effect damage model disclosed in this invention will first be described in detail. (See [link to relevant documentation]). Figure 1 The diagram shows a structural analysis method based on a damage model of coupled microcrack closure effect. This method can be executed by electronic devices and mainly includes the following steps S110 to S130: S110: A damage constitutive model was constructed based on the microcrack closure effect; In one embodiment, a damage constitutive model can be constructed for the specific application scenario of a rocket thrust chamber. This damage constitutive model can be a viscoplastic damage constitutive model that considers the microcrack closure effect. The construction of this damage constitutive model can be based on continuum damage mechanics. First, the crack closure effect coefficient is defined, then the state potential function coupled with the microcrack closure effect is derived, and finally, a viscoplastic damage constitutive model considering the microcrack closure effect is established based on the viscoplastic damage constitutive model. That is, constructing a damage constitutive model based on the microcrack closure effect can specifically include the following steps: (S11) Determine the crack closure effect coefficient based on the cross-sectional area parameters and tensile / compressive state of the material; The cross-sectional area parameters of the material include: the effective cross-sectional area of ​​microcracks or micropores on the cross-section, and the total cross-sectional area. As a concrete example, damage refers to the structural deterioration process caused by the breaking of atomic bonds and the plastic expansion of micro-defects under external loads. For a one-dimensional mean damage element, damage... It can be represented as:

[0025] Among them, S D Let S be the effective cross-sectional area of ​​the microcracks or micropores on the cross-section; S is the total cross-sectional area of ​​the cross-section. According to the principle of strain equivalence, the effective stress during tension... It can be represented as:

[0026] in, The cross-section is subjected to force during tension; This is the effective cross-sectional area during stretching; The stress is calculated without considering damage during stretching.

[0027] For most materials, internal micro-defects will partially or completely close under compressive loads. This closure increases the material's effective load-bearing area, leading to a certain degree of recovery in material properties. (Effective area under compression) satisfy:

[0028] Define the crack closure effect coefficient h:

[0029] The coefficient ranges from [0,1]. When h=0, it is considered that all microcracks are closed; when h=1, it is considered that no microcrack closure has occurred.

[0030] (S12) Based on the crack closure effect coefficient and the basic state potential function, determine the state potential function of the coupled microcrack closure effect; The fundamental state potential function includes the following formulas (5)-(9). As a specific example, the state potential function can be used to derive the state law equation of the material. To construct a viscoplastic damage constitutive model considering the microcrack closure effect, it is necessary to couple the microcrack closure effect with the state potential function. Based on the state dynamic coupling theory, without considering the microcrack closure effect, according to the assumption of the local state method, the state potential function can be expressed as:

[0031] in, ρ is the strain; r is the cumulative plastic strain due to damage; For back strain; E It is the elastic modulus; Poisson's ratio; Here are the temperature-dependent material parameters; tr(x) is the trace of the tensor x. .

[0032] In one-dimensional cases, the sign of the strain or stress value directly indicates whether the material is in a tensile or compressive state. In three-dimensional cases, the tensile or compressive state cannot be directly determined. Therefore, based on the signs of the principal strain components, the strain tensor is divided into positive and negative parts for tension and compression.

[0033]

[0034] In this context, the superscript "+" indicates a positive component, and the superscript "-" indicates a negative component; These are the eigenvalues ​​of the strain tensor, i.e., the principal values ​​of strain; The unit eigenvector of the strain tensor, i.e., the principal strain direction; For Macaulay brackets, The positive and negative components of strain and the original strain The following operational rules must be met:

[0035]

[0036] Substituting equations (8) and (9) into equation (5), and considering the crack closure effect coefficient h, we can obtain the state potential function of the coupled microcrack closure effect:

[0037] (S13) The material behavior is modeled based on state variables, the state potential function coupled with microcrack closure effect and the dissipation potential function to generate a damage constitutive model.

[0038] In irreversible thermodynamics, material behavior can be modeled using state variables, state potential functions, and dissipation potential functions.

[0039] Under the assumption of small deformation, strain can be decomposed into elastic strain. With inelastic strain :

[0040] The state law equations can be derived from the state potential function, and the adjoint variables of the state variables can be defined: in, For stress; R It is an isotropic strain hardening parameter; X Back stress; is the variable associated with the damage; I is the second-order unit tensor.

[0041] Based on the state-dynamic coupling theory, the dissipative potential function can be expressed as:

[0042] in, s It is a deviatoric stress; The initial yield strength of the material; S 1 、s 2 、m d These are temperature-dependent material parameters; p d The damage threshold is represented by the cumulative inelastic strain; p H(x) represents the cumulative inelastic strain; H(x) is the unit step function.

[0043] From the dissipative potential function, the evolution law of state variables related to the dissipation mechanism can be derived:

[0044] in, It is the plastic multiplier law.

[0045] The first three terms of the dissipative potential function are the Von Mises yield function. For viscoplastic models, we have:

[0046] in, This is viscous overstress.

[0047] According to Norton's power law formula, the relationship between viscous overstress and cumulative inelastic strain rate can be obtained:

[0048] in, K and n These are temperature-dependent material parameters.

[0049] Plastic multiplier With cumulative inelastic strain rate The relationship can be represented as:

[0050] Considering the partial recovery of metal hardening effect at high temperatures, under isotropic strain hardening parameters... and back stress A static recovery term was added to the evolution equation for correction, and the kinematic hardening phenomenon of the material was simulated using a superposition of three sets of back stresses, namely:

[0051] in, These are temperature-dependent material parameters; sign(x) is the sign function. This is the Von Mises equivalent back stress.

[0052] Based on the above formulas, a viscoplastic damage constitutive model that considers the microcrack closure effect can be established. This model has 26 material parameters that need to be identified, as shown in Table 1.

[0053] Table 1. Required identification parameters in the constitutive model

[0054] In one embodiment, the constructed viscoplastic damage constitutive model considering the microcrack closure effect generally includes seven initial independent variables; the initial independent variables include multiple initial material parameters; the material parameters include: damage-related parameters and non-damage-related parameters. As a specific example, the initial material parameters are shown in Table 1, wherein the damage-related parameters include: S 1 、s 2、m d 、p d Non-damage-related parameters include: .

[0055] After constructing the damage constitutive model, it needs to be numerically calculated. In one embodiment, the above method may further include a numerical implementation process for the viscoplastic damage constitutive model that considers the microcrack closure effect, namely: performing numerical calculations on the damage constitutive model to determine the target independent variables within the current time step corresponding to the initial independent variables.

[0056] The numerical representation of the aforementioned damage constitutive model can be achieved using an implicit Euler algorithm. The implicit Euler method transforms the initial value problem of a differential equation into an initial value problem of a difference equation, thus providing a high degree of numerical stability. Taking the fully implicit Euler algorithm as an example, the specific implementation process is as follows: The constitutive model constructed through the above steps contains 7 independent variables:

[0057] The corresponding system of differential equations is:

[0058] Discretizing equation (27) using the implicit Euler method yields the local residual equations:

[0059] For a nonlinear system of equations, if the first... The approximate solution has been obtained in step 1. In order to solve the 2nd step... n An approximate solution with +1 step can be obtained using the Newton-Raphson iteration:

[0060] in, This represents the number of iteration steps. The Jacobian matrix is ​​denoted as:

[0061] For three-dimensional analysis, elastic strain and back stress contain components in six directions. N =6; For two-dimensional plane strain state and two-dimensional plane stress state, take N =4.

[0062] After the local residual equations are calculated, the increments of each variable within the current time step are obtained. It can calculate each independent variable within the current time step:

[0063] By coupling the microcrack closure effect into the viscoplastic damage constitutive model through the above steps, a viscoplastic damage constitutive model considering the microcrack closure effect was successfully constructed. The numerical model of this damage constitutive model was then realized using an implicit Euler algorithm. The damage constitutive model coupled with the microcrack closure effect can be used to describe the effect of microcrack and micropore closure under compression in a structure, improving the accuracy of the constitutive model.

[0064] In practical applications, in order to apply the established constitutive model to the structural analysis and life prediction of the thrust chamber, material mechanical property tests can be carried out based on the material used for the inner wall of the thrust chamber, the operating temperature range of the thrust chamber, and the requirements for identifying constitutive model parameters.

[0065] S120: Parameter identification of the damage constitutive model is performed using the test results of the material's mechanical properties; In one embodiment, the step of identifying parameters of the damage constitutive model using the material mechanical property test results in S120 includes: (S21) Obtain stress-variable force data of materials under different states through a type of material mechanical property test; (S22) Based on stress-variable force data, determine the experimental values ​​of relevant parameters of material mechanical properties; One type of material mechanical property test includes: uniaxial tensile test, low-cycle fatigue test, and stress relaxation test. The uniaxial tensile test obtains stress-strain data of the material under monotonic tension; the low-cycle fatigue test obtains stress-strain data of the material under cyclic loading; data from both the uniaxial tensile test and the low-cycle fatigue test are used to identify parameters. S 1 、s 2 、m d 、p d。 Stress relaxation tests can be used to obtain stress data of materials under high-temperature loading; stress relaxation test data are used to identify viscosity-related parameters. Ignoring the viscosity of copper alloys below 700K, and assuming that the crack closure effect coefficient is independent of temperature, a constant value was selected. The experimental scheme is shown in Table 2.

[0066] Table 2 Summary of Material Mechanical Property Tests

[0067] (S23) Using the experimental values ​​of relevant parameters, the initial material parameters of the damage constitutive model are identified and optimized, and the optimized material parameters are determined; (S24) Based on the optimized material parameters, a viscoplastic damage constitutive model of the copper alloy is generated.

[0068] In one embodiment, the step of identifying and optimizing the initial material parameters of the damage constitutive model using the experimental values ​​of relevant parameters in (S23) above may include: firstly, using a random search algorithm to perform a global search on the simulation results of the damage constitutive model to determine the initial parameter vector; and then using a local optimization strategy to perform a fine search on the initial parameter vector to determine the target material parameters.

[0069] In other words, the PRINO method can be used for parameter identification. This method combines Price's global clustering-guided approach with a local optimization strategy, identifying constitutive model parameters by minimizing the distance between the model response and the experimental data. The PRINO method first uses a random search algorithm to perform a global search to find the initial parameter vector, and then uses a local optimization strategy for a refined search. This approach of global search followed by local optimization can significantly improve the accuracy and efficiency of parameter identification.

[0070] In one embodiment, the step of identifying and optimizing the initial material parameters of the damage constitutive model using the experimental values ​​of relevant parameters in (S23) above may further include: firstly, simulating a single material unit using the damage constitutive model to determine the first initial material parameter; the first initial material parameter is a non-damage-related parameter; then, simulating a full-size material sample model using the damage constitutive model, and determining the damage-related parameters of the target material based on the experimental values ​​of the non-damage-related parameter and the relevant parameter.

[0071] Because constitutive models of coupled material damage are time-consuming and computationally inefficient, parameter identification can be divided into two steps to improve efficiency: First, neglecting damage, a single material element is used for simulation to identify material parameters related to elasticity, viscosity, kinematic hardening, isotropic hardening, or softening. Then, a full-size material specimen model is used for simulation, identifying material parameters related to damage and microcrack closure effects based on uniaxial tensile and low-cycle fatigue test data. By first simulating a single material element and then simulating the entire full-size material, the computational load can be significantly reduced, improving identification efficiency while maintaining accuracy.

[0072] S130: Based on the damage constitutive model after parameter identification, structural analysis of the thrust chamber is performed.

[0073] In one embodiment, the step of performing structural analysis of the thrust chamber based on the damage constitutive model after parameter identification in S130 above includes: (S31) The temperature distribution of the inner wall of the thrust chamber was determined by three-dimensional flow and heat transfer coupling simulation; (S32) Determine the target section based on the temperature distribution results, and perform heat transfer analysis on the target section; (S33) Based on the results of heat transfer analysis, the structural analysis of the thrust chamber is carried out using the damage constitutive model after parameter identification.

[0074] As a specific example, taking the thrust chamber of a rocket engine as the research object, a three-dimensional flow-heat transfer coupled simulation was conducted to obtain the temperature distribution on the inner wall of the thrust chamber, with the highest temperature point designated as the critical section. Geometric parameters, temperature, pressure, and convective heat transfer coefficient of the critical section were extracted for thermal structural analysis.

[0075] The thermal structure analysis employs a sequential coupling method, first performing heat transfer analysis, and then conducting structural analysis based on the heat transfer analysis, without considering the influence of structural deformation on temperature distribution. The established constitutive model is used for multi-cycle structural analysis to obtain the stress-strain response and damage evolution law of the cooling channel at the critical section of the thrust chamber. Structural failure is considered when the damage value D reaches 1, and the number of cycles completed at this point represents the final structural life. It should be noted that in this example, for the rocket engine thrust chamber, a three-dimensional flow-heat transfer coupled simulation is used, which provides more precise localization compared to traditional one-dimensional heat transfer simulation.

[0076] This invention provides a structural analysis method based on a damage model with coupled microcrack closure effect, comprising: constructing a damage constitutive model based on the microcrack closure effect; identifying parameters of the damage constitutive model using material mechanical property test results; and performing structural analysis of the thrust chamber based on the damage constitutive model with identified parameters. This method solves the technical problems of low accuracy and low efficiency in thrust chamber structural analysis, and achieves the technical effects of improving the efficiency and accuracy of structural analysis and improving the accuracy of thrust chamber life prediction.

[0077] As a concrete example, combined Figure 2 As shown, this embodiment of the invention first (S210) constructs a viscoplastic damage constitutive model considering the microcrack closure effect; then (S220) derives the numerical implementation process of the constitutive model based on the fully implicit Euler method; next, (S230) conducts material mechanical property tests and identifies constitutive model parameters based on the PRINO method; finally (S240) based on the constructed constitutive model, a calculation program is written and applied to the thrust chamber structure analysis and life prediction. The model parameters to be identified and the material mechanical property tests to be conducted are shown in Tables 1 and 2 above, respectively.

[0078] In another embodiment, the constitutive model constructed in the above method can also be programmed to facilitate model analysis and application. That is, based on the numerical implementation process of the constitutive model based on the fully implicit Euler method described above, a UMAT material subroutine is written in Fortran language. After linking and compiling through the User Programmable Features module of ANSYS, the viscoplastic damage constitutive model considering the microcrack closure effect established in this invention can be called in ANSYS to carry out structural analysis.

[0079] This embodiment provides a structural analysis method based on a coupled microcrack closure effect damage model. By selecting a suitable optimization algorithm and simplifying the parameter identification process, the efficiency of parameter identification is improved while ensuring the accuracy of parameter identification.

[0080] Furthermore, embodiments of the present invention also provide a structural analysis device based on a damage model of coupled microcrack closure effect, see [link to relevant documentation]. Figure 3 As shown, the device includes: Model building module 310 is used to build a damage constitutive model based on the microcrack closure effect; The parameter identification module 320 is used to identify parameters of the damage constitutive model using the test results of the material mechanical properties. The structural analysis module 330 is used to perform structural analysis on the thrust chamber based on the damage constitutive model after parameter identification.

[0081] The structural analysis device based on the coupled microcrack closure effect damage model provided in this application embodiment can be specific hardware on a device or software or firmware installed on the device. The implementation principle and technical effects of the device provided in this application embodiment are the same as those in the foregoing method embodiments. For the sake of brevity, any parts not mentioned in the device embodiment can be referred to the corresponding content in the foregoing method embodiments. Those skilled in the art will clearly understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can all be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here. The structural analysis device based on the coupled microcrack closure effect damage model provided in this application embodiment has the same technical features as the structural analysis method based on the coupled microcrack closure effect damage model provided in the foregoing embodiments, and therefore can solve the same technical problems and achieve the same technical effects.

[0082] This application also provides an electronic device, specifically, the electronic device includes a processor and a storage device; the storage device stores a computer program, and the computer program, when run by the processor, executes the method described in any of the above embodiments.

[0083] Figure 4This is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. The electronic device 400 includes: a processor 40, a memory 41, a bus 42, and a communication interface 43. The processor 40, the communication interface 43, and the memory 41 are connected through the bus 42. The processor 40 is used to execute executable modules, such as computer programs, stored in the memory 41.

[0084] The memory 41 may include high-speed random access memory (RAM) or non-volatile memory, such as at least one disk storage device. Communication between this system network element and at least one other network element is achieved through at least one communication interface 43 (which can be wired or wireless), such as the Internet, wide area network, local area network, metropolitan area network, etc.

[0085] Bus 42 can be an ISA bus, PCI bus, or EISA bus, etc. The bus can be divided into address bus, data bus, control bus, etc. For ease of representation, Figure 4 The symbol is represented by a single double-headed arrow, but this does not mean that there is only one bus or one type of bus.

[0086] The memory 41 is used to store programs. After receiving an execution instruction, the processor 40 executes the program. The method executed by the apparatus of the process definition disclosed in any of the foregoing embodiments of the present invention can be applied to the processor 40 or implemented by the processor 40.

[0087] Processor 40 may be an integrated circuit chip with signal processing capabilities. In implementation, each step of the above method can be completed by the integrated logic circuitry in the hardware of processor 40 or by instructions in software form. Processor 40 can be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), etc.; it can also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field-Programmable Gate Array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components. It can implement or execute the methods, steps, and logic block diagrams disclosed in the embodiments of this invention. The general-purpose processor can be a microprocessor or any conventional processor. The steps of the methods disclosed in the embodiments of this invention can be directly embodied in the execution of a hardware decoding processor, or executed by a combination of hardware and software modules in the decoding processor. The software modules can reside in random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, registers, or other mature storage media in the art. The storage medium is located in memory 41. The processor 40 reads the information in memory 41 and, in conjunction with its hardware, completes the steps of the above method.

[0088] Corresponding to the above method, this application embodiment also provides a computer-readable storage medium storing computer-executable instructions. When the computer-executable instructions are invoked and executed by a processor, the computer-executable instructions cause the processor to perform the steps of the above method.

[0089] In the embodiments provided in this application, it should be understood that the disclosed apparatus and methods can be implemented in other ways. The apparatus embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the displayed or discussed mutual couplings, direct couplings, or communication connections may be through some communication interfaces; indirect couplings or communication connections between devices or units may be electrical, mechanical, or other forms.

[0090] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0091] In addition, the functional units in the embodiments provided in this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0092] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention, essentially, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, electronic device, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0093] It should be noted that similar reference numerals and letters in the accompanying drawings indicate similar items. Therefore, once an item is defined in one accompanying drawing, it does not need to be further defined and explained in subsequent accompanying drawings. In addition, the terms "first," "second," "third," etc. are used only to distinguish descriptions and should not be construed as indicating or implying relative importance.

[0094] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A structural analysis method based on a damage model of coupled microcrack closure effect, characterized in that, include: A damage constitutive model was constructed based on the microcrack closure effect; The parameters of the damage constitutive model are identified using the results of material mechanical property tests. Based on the damage constitutive model after parameter identification, the thrust chamber is structurally analyzed.

2. The method according to claim 1, characterized in that, A damage constitutive model based on the microcrack closure effect is constructed, including: The crack closure effect coefficient is determined based on the cross-sectional area parameters and tensile / compressive state of the material. Based on the crack closure effect coefficient and the basic state potential function, the state potential function of the coupled microcrack closure effect is determined. Based on the state variables, the state potential function of the coupled microcrack closure effect, and the dissipation potential function, the material behavior is modeled to generate a damage constitutive model.

3. The method according to claim 2, characterized in that, The damage constitutive model is a viscoplastic damage constitutive model that considers the microcrack closure effect; the damage constitutive model includes initial independent variables; The initial independent variables include multiple initial material parameters; the material parameters include: damage-related parameters and non-damage-related parameters; the method further includes: The damage constitutive model is numerically calculated to determine the target independent variable within the current time step corresponding to the initial independent variable.

4. The method according to claim 1, characterized in that, Using the results of material mechanical property tests, parameter identification is performed on the damage constitutive model, including: Stress-variable force data of a material under different states are obtained through mechanical property tests of a class of materials. Based on the stress-variable force data, the experimental values ​​of relevant parameters of the material's mechanical properties are determined; Using the experimental values ​​of the relevant parameters, the initial material parameters of the damage constitutive model are identified and optimized, and the optimized material parameters are determined. Based on the optimized material parameters, a viscoplastic damage constitutive model for the copper alloy is generated.

5. The method according to claim 4, characterized in that, Using the experimental values ​​of the relevant parameters, the initial material parameters of the damage constitutive model are identified and optimized, including: The initial parameter vector is determined by globally searching the simulation results of the damage constitutive model using a random search algorithm. The initial parameter vector is finely searched using a local optimization strategy to determine the target material parameters.

6. The method according to claim 4, characterized in that, Using the experimental values ​​of the relevant parameters, the initial material parameters of the damage constitutive model are identified and optimized, including: The damage constitutive model is used to simulate a single material element to determine the first initial material parameters; the first initial material parameters are non-damage-related parameters. The damage constitutive model is used to simulate a full-size material specimen model. Based on the non-damage-related parameters and the experimental values ​​of the related parameters, the damage-related parameters of the target material are determined.

7. The method according to claim 1, characterized in that, Based on the damage constitutive model after parameter identification, structural analysis of the thrust chamber is performed, including: The temperature distribution on the inner wall of the thrust chamber was determined by three-dimensional flow-heat transfer coupling simulation. The target cross section is determined based on the temperature distribution results, and heat transfer analysis is performed on the target cross section. Based on the results of the heat transfer analysis, the thrust chamber is structurally analyzed using the damage constitutive model after parameter identification.

8. A structural analysis device based on a damage constitutive model of coupled microcrack closure effect, characterized in that, include: The model building module is used to construct a damage constitutive model based on the microcrack closure effect. The parameter identification module is used to identify parameters of the damage constitutive model using the test results of the material mechanical properties. The structural analysis module is used to perform structural analysis on the thrust chamber based on the damage constitutive model after parameter identification.

9. An electronic device comprising a memory and a processor, wherein the memory stores a computer program executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer-executable instructions that, when invoked and executed by a processor, cause the processor to perform the method according to any one of claims 1 to 7.

Images