Structural analysis method and device based on coupled micro-crack closure effect damage model

Abstract

The application provides a structure analysis method and device based on a coupling micro-crack closure effect damage model, and relates to the technical field of rocket engines.The method comprises the following steps: constructing a damage constitutive model based on a micro-crack closure effect; using material mechanical property test results to perform parameter identification on the damage constitutive model; and performing structure analysis on a thrust chamber based on the damage constitutive model after parameter identification.The method solves the technical problems of low accuracy and low efficiency of thrust chamber structure analysis, achieves the technical effects of improving structure analysis efficiency and precision, and improves the accuracy of thrust chamber life prediction.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of rocket engine, and in particular to a structure analysis method and device based on a damage model coupled with micro-crack closure effect. BACKGROUND

[0002] Reusability is a significant feature of the new generation of launch vehicles, and advanced liquid power is an important link to achieve reusability. In the working process of the reusable liquid rocket engine multiple start-stop, the thrust chamber as the core component of the reusable rocket engine, bears complex alternating thermal mechanical load, and the inner wall structure material shows the cyclic characteristics of kinematic hardening, isotropic hardening or softening, and is accompanied by time-dependent inelastic deformation and stress relaxation phenomena. Moreover, the material inevitably has micro-cracks, micro-pores and other micro-defects, and the micro-cracks and micro-pores expand under tension and close under compression, resulting in micro-crack closure effect. The micro-crack closure effect will affect the structural response of the thrust chamber in the subsequent working process, and further affect the accuracy of life prediction.

[0003] However, the existing technology often ignores the effect of micro-crack closure effect in the application of using the constitutive model to realize the structure analysis of the thrust chamber, resulting in low accuracy of the analysis results. Moreover, since the constitutive model contains numerous temperature-dependent material parameters, the parameter identification through model simulation is time-consuming and low in calculation efficiency, which greatly affects the efficiency of the structure analysis of the thrust chamber and the subsequent life prediction. SUMMARY

[0004] The purpose of the present application is to provide a structure analysis method and device based on a damage model coupled with micro-crack closure effect, to solve the technical problems of low accuracy and low efficiency of the structure analysis of the thrust chamber in the prior art.

[0005] In the first aspect, the embodiments of the present application provide a structure analysis method based on a damage model coupled with micro-crack closure effect, which comprises: constructing a damage constitutive model based on the micro-crack closure effect; identifying parameters of the damage constitutive model by using the material mechanical property test results; and performing structure analysis on the thrust chamber based on the damage constitutive model after parameter identification.

[0006] In some optional implementations, constructing a damage constitutive model based on the micro-crack closure effect comprises: determining a crack closure effect coefficient according to the cross-sectional area parameter of the material and the tension-compression state; determining a state potential function coupled with the micro-crack closure effect based on the crack closure effect coefficient and a basic state potential function; and modeling the material behavior based on the state variable, the state potential function coupled with the micro-crack closure effect, and the dissipation potential function to generate the damage constitutive model.

[0007] In some optional implementations, the damage constitutive model is a visco-plastic damage constitutive model considering micro-crack closure effect; the damage constitutive model includes initial independent variables; the initial independent variables include a plurality of initial material parameters; the material parameters include damage-related parameters and non-damage-related parameters; and the method further includes numerically calculating the damage constitutive model to determine target independent variables corresponding to the initial independent variables in a current time step.

[0008] In some optional implementations, the initial material parameters of the damage constitutive model are identified and optimized using the test values of the related parameters, including: obtaining stress and strain data of the material in different states through a type of material mechanical property test; determining test values of material mechanical property related parameters based on the stress and strain data; identifying and optimizing the initial material parameters of the damage constitutive model using the test values of the related parameters to determine optimized material parameters; and generating a visco-plastic damage constitutive model of the copper alloy based on the optimized material parameters.

[0009] In some optional implementations, the initial material parameters of the damage constitutive model are identified and optimized using the test values of the related parameters, including: performing global search on simulation results of the damage constitutive model using a random search algorithm to determine an initial parameter vector; and performing fine search on the initial parameter vector using a local optimization strategy to determine target material parameters.

[0010] In some optional implementations, the initial material parameters of the damage constitutive model are identified and optimized using the test values of the related parameters, including: simulating a single material unit using the damage constitutive model to determine first initial material parameters; the first initial material parameters are non-damage-related parameters; and simulating a full-size material sample model using the damage constitutive model to determine damage-related parameters of target material according to the non-damage-related parameters and the test values of the related parameters.

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

[0012] In a second aspect, an embodiment of the present application provides a structure analysis device based on a damage model of coupled micro-crack closure effect, which comprises: a model construction module, configured to construct a damage constitutive model based on micro-crack closure effect; a parameter identification module, configured to identify parameters of the damage constitutive model by using material mechanical property test results; and a structure analysis module, configured to perform structure analysis on the thrust chamber based on the damage constitutive model after parameter identification.

[0013] In a third aspect, an embodiment of the present application provides an electronic device, comprising a memory and a processor, wherein the memory stores a computer program capable of running on the processor, and the processor implements steps of the method of any one of the first aspect when executing the computer program.

[0014] In a fourth aspect, an embodiment of the present application provides a computer readable storage medium, which stores computer executable instructions, and the computer executable instructions cause a processor to execute the method of any one of the first aspect when the computer executable instructions are invoked and run by the processor.

[0015] The present application provides a structure analysis method and device based on a damage model of coupled micro-crack closure effect, which comprises: constructing a damage constitutive model based on micro-crack closure effect; identifying parameters of the damage constitutive model by using material mechanical property test results; and performing structure analysis on the thrust chamber based on the damage constitutive model after parameter identification, which solves the technical problems of low accuracy and low efficiency of thrust chamber structure analysis, and achieves the technical effects of improving structure analysis efficiency and precision and improving the accuracy of thrust chamber life prediction. BRIEF DESCRIPTION OF DRAWINGS

[0016] In order to more clearly illustrate the specific embodiments of the present application or the technical solutions in the prior art, the drawings needed in the specific embodiments or prior art description will be briefly introduced as follows. Obviously, the drawings in the following description are some embodiments of the present application, and those skilled in the art can also obtain other drawings according to these drawings without creative labor.

[0017] Figure 1 A flowchart of a structure analysis method based on a damage model of coupled micro-crack closure effect is provided for an embodiment of the present application.

[0018] Figure 2 A flowchart of another structure analysis method based on a damage model of coupled micro-crack closure effect is provided for an embodiment of the present application.

[0019] Figure 3A structural schematic diagram of a structure analysis device based on a coupling micro-crack closure effect damage model is provided for an embodiment of the present application.

[0020] Figure 4 A structural schematic diagram of an electronic device is provided for an embodiment of the present application. DETAILED DESCRIPTION

[0021] To make the objectives, technical solutions and advantages of the embodiments of the present application clearer, the following will clearly and completely describe the technical solutions in the embodiments of the present application with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some but not all of the embodiments of the present application. The components of the embodiments of the present application described and shown in the drawings can be arranged and designed in various different configurations.

[0022] Therefore, the following detailed description of the embodiments of the present application provided in the drawings is not intended to limit the scope of the claimed present application, but only represents selected embodiments of the present application. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the present application without creative labor are within the scope of protection of the present application.

[0023] It should be noted that: similar reference numbers and letters represent similar items in the following drawings, therefore, once an item is defined in one drawing, it does not need to be further defined and explained in subsequent drawings. The following will describe some embodiments of the present application in detail with reference to the drawings. The following embodiments and features in the embodiments can be combined with each other without conflict.

[0024] Under repeated loading conditions, such as multiple ignition and shutdown cycles of a rocket engine, micro-cracks can form and propagate in the inner wall of the thrust chamber. Micro-crack closure effect can slow down the growth rate of the cracks, thereby increasing the fatigue life of the thrust chamber. However, if the micro-cracks propagate to a certain size, it may cause irreversible damage, affecting the safety and reliability of the thrust chamber. That is, the micro-crack closure effect will affect the structural response of the thrust chamber in the subsequent working process, and thus affect the accuracy of life prediction. Therefore, the micro-crack closure effect needs to be considered, and a constitutive model capable of describing the mechanical behavior of the material in the inner wall of the thrust chamber is established to lay the foundation for the life prediction of the thrust chamber.

[0025] In the application of the prior art in realizing the thrust chamber structure analysis by using the constitutive model, the effect of the micro-crack closure effect is often ignored, resulting in low accuracy of the analysis results. And because the constitutive model contains numerous temperature-related material parameters, which need to be identified through material mechanical property tests and with the help of optimization algorithms, the simulation of the constitutive model coupled with material damage is time-consuming and low in calculation efficiency, thereby greatly affecting the efficiency of the thrust chamber structure analysis and subsequent life prediction.

[0026] Based on this, the embodiment of the present application provides a structure analysis method and device based on a coupled micro-crack closure effect damage model to solve the technical problems of low accuracy and low efficiency of the structure analysis of the existing technology.

[0027] In order to facilitate the understanding of the present embodiment, first of all, a structure analysis method based on a coupled micro-crack closure effect damage model is introduced in detail, as shown in Figure 1 The flowchart of the structure analysis method based on the coupled micro-crack closure effect damage model, which can be executed by an electronic device, mainly includes the following steps S110 to S130:

[0028] S110: constructing a damage constitutive model based on micro-crack closure effect;

[0029] In an embodiment, a damage constitutive model can be constructed for the specific application scenario of the rocket thrust chamber, which can be a visco-plastic damage constitutive model considering micro-crack closure effect. The construction of the damage constitutive model can be based on continuum damage mechanics, first define the crack closure effect coefficient, then derive the state potential function coupled with the micro-crack closure effect, and finally establish the visco-plastic damage constitutive model considering the micro-crack closure effect on the basis of the visco-plastic damage constitutive model. That is, constructing a damage constitutive model based on micro-crack closure effect can specifically include the following steps:

[0030] (S11) determining the crack closure effect coefficient according to the cross-sectional area parameters and the tension-compression state of the material;

[0031] Among them, the cross-sectional area parameters of the material include: the effective cross-sectional area of the micro-crack or micro-hole on the cross-section, and the total cross-sectional area of the cross-section. As a specific example, damage refers to the structural degradation process caused by atomic bond rupture and plastic expansion of micro-defects under external load. For a one-dimensional mean damage element, the damage can be expressed as:

[0032]

[0033] Among them, S D is the effective cross-sectional area of the micro-crack or micro-hole on the cross-section; S is the total cross-sectional area of the cross-section. According to the strain equivalence principle, the effective stress when stretching can be expressed as:

[0034]

[0035] Among them, is the cross-sectional force when stretching; is the effective cross-sectional area when stretching; The stress when not considering damage during stretching.

[0036] For most materials, the internal micro-defects will partially or completely close under compression load, and the closure of micro-defects can increase the effective load-bearing area of the material, and the material performance will be restored to a certain extent. The effective area during compression Satisfies:

[0037]

[0038] Define the crack closure effect coefficient h:

[0039]

[0040] The value range of this coefficient is [0, 1], when h=0, it is considered that the micro-cracks are completely closed; when h=1, it is considered that no micro-crack closure phenomenon occurs.

[0041] (S12) Based on the crack closure effect coefficient and the basic state potential function, determine the state potential function coupled with the micro-crack closure effect;

[0042] Wherein, the basic state potential function includes the following formulas (5)-(9). As a specific example, the state potential function can derive the state law equation of the material, and the viscoplastic damage constitutive model considering the micro-crack closure effect needs to couple the micro-crack closure effect with the state potential function. Based on the state dynamic coupling theory, without considering the micro-crack closure effect, according to the assumption of the local state method, the state potential function can be expressed as:

[0043]

[0044] Wherein, is the strain; r is the damage cumulative plastic strain; is the back strain; E is the elastic modulus; is the Poisson's ratio; is the material parameter related to temperature; tr(x) is the trace of tensor x, .

[0045] In one-dimensional case, the tensile or compression state of the material can be directly judged by the positive and negative of the strain or stress value. For three-dimensional case, the tensile or compression state of the material cannot be directly judged, according to the positive and negative of the principal strain component, the strain tensor is divided into positive and negative parts according to tensile and compression:

[0046]

[0047]

[0048] where the superscript "+" denotes the positive component and the superscript "-" denotes the negative component; is the eigenvalue of the strain tensor, i.e., the principal strain; is the unit eigenvector of the strain tensor, i.e., the principal direction of strain; is the Macaulay bracket, The positive and negative components of strain are related to the original strain satisfy the following operation rules:

[0049]

[0050]

[0051] Substituting equation (8) and equation (9) into equation (5) and considering the crack closure effect coefficient h, the state potential function coupled with micro-crack closure effect is obtained:

[0052]

[0053] (S13) Based on the state variable, the state potential function coupled with micro-crack closure effect, and the dissipation potential function, the material behavior is modeled, and a damage constitutive model is generated.

[0054] In the theory of irreversible thermodynamics, the material behavior can be modeled by state variables, state potential functions, and dissipation potential functions.

[0055] Under the small deformation assumption, the strain can be decomposed into elastic strain and inelastic strain :

[0056]

[0057] The state law equation can be derived from the state potential function, and the concomitant variable of the state variable is defined:

[0058] where is the stress; R is the isotropic strain hardening parameter; X is the back stress; is the damage concomitant variable; I is the second-order unit tensor.

[0059] Based on the state dynamic coupling theory, the dissipation potential function can be expressed as:

[0060]

[0061] where s is the deviatoric stress; is the initial yield strength of the material; S 1 、s 2 、md material parameters related to temperature; p d damage threshold represented by accumulated inelastic strain; p accumulated inelastic strain; H(x) is a unit step function.

[0062] The state variable evolution law related to dissipation mechanism can be derived from the dissipation potential function:

[0063]

[0064] where, plastic multiplier rate.

[0065] The first three terms of the dissipation potential function are Von Mises yield function For viscoplastic model, we have:

[0066]

[0067] where, viscous overstress.

[0068] According to the Norton power law formula, the relationship between viscous overstress and accumulated inelastic strain rate can be obtained:

[0069]

[0070] where, K and n material parameters related to temperature.

[0071] The plastic multiplier rate and the accumulated inelastic strain rate can be expressed as:

[0072]

[0073] Considering the partial recovery of metal hardening effect at high temperature, static recovery terms are added to the evolution equations of isotropic strain hardening parameter and back stress , and the form of superposition of three groups of back stress is used to simulate the follow-up hardening phenomenon of the material, that is:

[0074]

[0075] where, material parameters related to temperature; sign(x) is a sign function; Von Mises equivalent back stress.

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

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

[0078]

[0079] 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: .

[0080] 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.

[0081] 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:

[0082] The constitutive model constructed through the above steps contains 7 independent variables:

[0083]

[0084] The corresponding system of differential equations is:

[0085]

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

[0087]

[0088] 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:

[0089]

[0090] wherein, is the iteration step number; is the Jacobian matrix, denoted as:

[0091]

[0092] For three-dimensional analysis, the elastic strain and back stress contain 6 directional components, taking N = 6; for two-dimensional plane strain state and two-dimensional plane stress state, taking N = 4.

[0093] After the local residual equation system is calculated, the increment of each variable in the current time step is obtained , and each independent variable in the current time step can be calculated:

[0094]

[0095] By the above steps, the micro-crack closure effect is coupled into the viscoplastic damage constitutive model, and a viscoplastic damage constitutive model considering the micro-crack closure effect is successfully constructed, and the numericalization of the damage constitutive model is realized through the implicit Euler algorithm. The damage constitutive model coupled with the micro-crack closure effect can be used to describe the closure effect of micro-cracks and micro-holes in the structure under compression, and the accuracy of the constitutive model is improved.

[0096] In actual application, in order to apply the established constitutive model to the structure analysis and life prediction of the thrust chamber, material mechanical property tests can be carried out according to the material used for the inner wall of the thrust chamber, the working temperature range of the thrust chamber and the requirement of constitutive model parameter identification.

[0097] S120: parameter identification of the damage constitutive model by using the material mechanical property test results;

[0098] In one embodiment, in the above S120, the step of parameter identification of the damage constitutive model by using the material mechanical property test results includes:

[0099] (S21) obtaining stress-strain data of the material in different states through a kind of material mechanical property test;

[0100] (S22) determining the test values of the material mechanical property related parameters based on the stress-strain data;

[0101] The uniaxial tensile test can obtain the stress-strain data of the material under monotonic tension; the low-cycle fatigue test can obtain the stress-strain data of the material under cyclic loading; the uniaxial tensile test and the low-cycle fatigue test data are used to identify parameters simultaneously: S 1 、s 2 、m d 、p d。 The stress relaxation test can obtain the stress data of the material under high-temperature preservation; the stress relaxation test data are used to identify the viscous related parameters: Without considering the viscosity of the copper alloy below 700K, the crack closure effect coefficient is considered to be independent of temperature, and is selected as a constant value. The test scheme is shown in Table 2.

[0102] Table 2 Summary of material mechanical property tests

[0103]

[0104] (S23) The initial material parameters of the damage constitutive model are identified and optimized using the test values of the related parameters, and the optimized material parameters are determined.

[0105] (S24) Based on the optimized material parameters, a viscoplastic damage constitutive model of the copper alloy is generated.

[0106] In an embodiment, the step of identifying and optimizing the initial material parameters of the damage constitutive model using the test values of the related parameters in (S23) above can include: first using a random search algorithm to perform global search on the simulation results of the damage constitutive model to determine an initial parameter vector; and then using a local optimization strategy to perform fine search on the initial parameter vector to determine the target material parameters.

[0107] That is, the PRINO method can be used for parameter identification, which combines Price's global clustering-oriented method and local optimization strategy to identify the constitutive model parameters by minimizing the distance between the model response and the test data. The PRINO method first uses a random search algorithm for global search to find an initial parameter vector, and then uses a local optimization strategy for fine search. In this way, the accuracy and efficiency of parameter identification can be significantly improved by first performing global search and then performing local optimization.

[0108] In an embodiment, the step of identifying and optimizing the initial material parameters of the damage constitutive model using the test values of the related parameters in (S23) above can further include: first simulating a single material unit using the damage constitutive model to determine first initial material parameters; the first initial material parameters are non-damage related parameters; and then simulating the full-size material sample model using the damage constitutive model to determine the damage related parameters of the target material according to the non-damage related parameters and the test values of the related parameters.

[0109] Since the simulation using the constitutive model coupled with material damage is time-consuming and has low computational efficiency, to improve the parameter identification efficiency, the parameter identification process can be divided into two steps: first, ignore the damage and use a material unit for simulation to identify the material parameters related to elasticity, viscosity, hardening, isotropic hardening or softening; and then use the full-size material sample model for simulation to identify the material parameters related to damage and micro-crack closure effect according to the uniaxial tensile and low-cycle fatigue test data. By first simulating a material unit and then simulating the full-size material as a whole, the computational load can be greatly reduced, and the identification efficiency can be improved while ensuring the identification accuracy.

[0110] S130: performing structural analysis on the thrust chamber based on the damage constitutive model after parameter identification.

[0111] In an embodiment, the step of performing structural analysis on the thrust chamber based on the damage constitutive model after parameter identification in S130 above includes:

[0112] (S31) determining the temperature distribution result of the inner wall surface of the thrust chamber through three-dimensional flow-heat coupled simulation;

[0113] (S32) determining a target cross-section according to the temperature distribution result and performing heat transfer analysis on the target cross-section;

[0114] (S33) performing structural analysis on the thrust chamber using the damage constitutive model after parameter identification in combination with the results of the heat transfer analysis.

[0115] As a specific example, taking a rocket engine thrust chamber as the research object, three-dimensional flow-heat coupled simulation is carried out to obtain the temperature distribution result of the inner wall surface of the thrust chamber, and the highest temperature position is taken as the dangerous cross-section. The geometric parameters, temperature, pressure, and convective heat transfer coefficient of the dangerous cross-section are extracted to carry out thermal-structural analysis.

[0116] The thermal structural analysis adopts a sequential coupling method, and is first performed by heat transfer analysis, and then is performed by structural analysis on the basis of the heat transfer analysis, without considering the influence of structural deformation on temperature distribution. The established constitutive model is applied to structural analysis of multiple working cycles, to obtain stress-strain response and damage evolution law of the cooling channel of the dangerous section of the thrust chamber, and when the damage value D reaches 1, the structural failure is considered, and at this time, the number of working cycles experienced is the final structural life. It should be noted that in this example, for the rocket engine thrust chamber, three-dimensional flow heat transfer coupling simulation is adopted, compared with the traditional one-dimensional heat transfer simulation, the positioning is more accurate.

[0117] The application provides a structure analysis method based on a coupling micro-crack closure effect damage model, which comprises the following steps: constructing a damage constitutive model based on the micro-crack closure effect; performing parameter identification on the damage constitutive model by using material mechanical property test results; and performing structural analysis on a thrust chamber based on the damage constitutive model after parameter identification. The method solves the technical problems of low accuracy and low efficiency of thrust chamber structural analysis, improves the efficiency and accuracy of structural analysis, and improves the accuracy of thrust chamber life prediction.

[0118] As a specific example, in combination with Figure 2 As shown in the figure, the embodiment of the application first constructs a viscoplastic damage constitutive model considering the micro-crack closure effect (S210); then derives the numerical implementation process of the constitutive model based on the fully implicit Euler method (S220); then performs material mechanical property test (S230), and performs constitutive model parameter identification based on the PRINO method; finally, based on the constructed constitutive model, a calculation program is written and applied to thrust chamber structural analysis and life prediction (S240). The required model parameters to be identified and the required material mechanical property tests are shown in Table 1 and Table 2 above, respectively.

[0119] In another embodiment, the constitutive model constructed in the above method can also be programmed, so as to facilitate the analysis and application of the model. That is, according to the numerical implementation process of the constitutive model based on the fully implicit Euler method, a UMAT material subroutine is written by using Fortran language, and after being linked and compiled by the User Programmable Features module of ANSYS, the viscoplastic damage constitutive model considering the micro-crack closure effect established by the application can be called in ANSYS to perform structural analysis.

[0120] The structure analysis method based on the coupling micro-crack closure effect damage model provided in the embodiment improves the efficiency of parameter identification on the basis of ensuring the accuracy of parameter identification by selecting a suitable optimization algorithm and simplifying the parameter identification process.

[0121] In addition, the embodiment of the present application further provides a structure analysis device based on the damage model of coupled micro-crack closure effect, referring to Figure 3 As shown in the figure, the device comprises:

[0122] The model construction module 310 is configured to construct a damage constitutive model based on the micro-crack closure effect.

[0123] The parameter identification module 320 is configured to perform parameter identification on the damage constitutive model by using the material mechanical property test results.

[0124] The structure analysis module 330 is configured to perform structure analysis on the thrust chamber based on the damage constitutive model after parameter identification.

[0125] The device based on the damage model of coupled micro-crack closure effect provided by the embodiment of the present application can be specific hardware on the equipment or software or firmware installed on the equipment, etc. The device provided by the embodiment of the present application has the same implementation principle and technical effects as the foregoing method embodiments. For brevity and conciseness, the part of the device embodiment not mentioned in the foregoing method embodiments can be referred to the corresponding content in the foregoing method embodiments. Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process of the system, device and unit described above can be referred to the corresponding process in the foregoing method embodiments, which will not be described herein. The device based on the damage model of coupled micro-crack closure effect provided by the embodiment of the present application has the same technical features as the structure analysis method based on the damage model of coupled micro-crack closure effect provided by the foregoing embodiment, so it can also solve the same technical problems and achieve the same technical effects.

[0126] The embodiment of the present application further provides an electronic device, specifically, the electronic device comprises a processor and a storage device; the storage device stores a computer program, and the computer program performs the method described in any one of the above embodiments when executed by the processor.

[0127] Figure 4 A structural schematic diagram of an electronic device provided by the embodiment of the present application is shown in the figure, the electronic device 400 comprises 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 configured to execute the executable modules stored in the memory 41, such as a computer program.

[0128] The memory 41 can include a high-speed random access memory (RAM), and can also include a non-volatile memory, such as at least one disk memory. The communication connection between the system network element and at least one other network element is realized through at least one communication interface 43 (which can be wired or wireless), and the Internet, a wide area network, a local area network, a metropolitan area network, etc. can be used.

[0129] The bus 42 can be an ISA bus, a PCI bus, an EISA bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 4 Only one bidirectional arrow is used in the figure to represent the bus, but it does not mean that there is only one bus or only one type of bus.

[0130] The memory 41 is used to store programs, and the processor 40 executes the programs after receiving execution instructions. The method performed by the device defined by the flow of any embodiment of the present application can be applied to the processor 40 or implemented by the processor 40.

[0131] The processor 40 can be an integrated circuit chip with signal processing capability. In the implementation process, each step of the above method can be completed by the integrated logic circuit of hardware in the processor 40 or the instruction in the form of software. The processor 40 described above 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 gates or transistor logic devices, discrete hardware components. It can implement or execute the disclosed methods, steps and logic block diagrams in the embodiments of the present application. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor. The steps of the method disclosed in combination with the embodiments of the present application can be directly embodied as a hardware code processor for execution, or a combination of hardware and software modules in the code processor for execution. The software module can be located in a random access memory, a flash memory, a read-only memory, a programmable read-only memory or an electrically erasable programmable memory, a register, etc. The storage medium in the art. The storage medium is located in the memory 41, and the processor 40 reads the information in the memory 41 and combines the hardware to complete the steps of the above method.

[0132] Corresponding to the above method, the embodiment of the present application also provides a computer readable storage medium, the computer readable storage medium stores computer executable instructions, when the computer executable instructions are called and run by a processor, the computer executable instructions cause the processor to run the steps of the above method.

[0133] In the embodiments provided by the present application, it should be understood that the disclosed device and method can be implemented in other ways. The device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There can be another division manner for the actual implementation, for example, a plurality of units or components can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the displayed or discussed mutual couplings or direct couplings or communication connections between the units can be indirect couplings or communication connections through some interfaces, devices or units, and can be electrical, mechanical or in other forms.

[0134] The units described as separate components can or can not be physically separate, and the components displayed as units can or can not be physical units, that is, they can be located in one place, or can be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the embodiment scheme.

[0135] In addition, each functional unit in the embodiments provided by the present application can be integrated in one processing unit, or each unit can be a physical unit, or two or more units can be integrated in one unit.

[0136] If the functions are realized in the form of software function units and sold or used as independent products, they can be stored in a computer readable storage medium. Based on this understanding, the technical solutions of the present application essentially or the parts that make contributions to the prior art or parts of the technical solutions can be embodied in the form of a software product. The computer software product is stored in a storage medium, and includes a number of instructions for causing a computer device (which can be a personal computer, an electronic device, or a network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present application. The foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM, Read-Only Memory), a random access memory (RAM, Random Access Memory), a magnetic disk or an optical disk, and various media that can store program codes.

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

[0138] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, and are not limited thereto. Although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that the technical solutions recorded in the foregoing embodiments can still be modified, or some or all of the technical features can be replaced by equivalents. The modification or replacement does not make the essence of the corresponding technical solution deviate from the scope of the technical solutions of the embodiments of the present application.

Claims

1. A structural analysis method based on a coupled microcrack closure effect damage model, characterized by, The method comprises the following steps: constructing a damage constitutive model based on micro-crack closure effect; performing parameter identification on the damage constitutive model by using material mechanical property test results; performing structural analysis on a thrust chamber based on the damage constitutive model after parameter identification; constructing a damage constitutive model based on micro-crack closure effect, comprising: determining a crack closure effect coefficient according to a cross-sectional area parameter and a tension-compression state of a material; determining a state potential function coupled with micro-crack closure effect based on the crack closure effect coefficient and a basic state potential function; modeling material behavior based on a state variable, the state potential function coupled with micro-crack closure effect, and a dissipation potential function to generate a damage constitutive model.

2. The method of claim 1, wherein, The damage constitutive model is a visco-plastic damage constitutive model considering micro-crack closure effect; the damage constitutive model comprises initial independent variables; the initial independent variables comprise a plurality of initial material parameters; the material parameters comprise damage-related parameters and non-damage-related parameters; the method further comprises: performing numerical calculation on the damage constitutive model to determine target independent variables corresponding to the initial independent variables in a current time step.

3. The method of claim 1, wherein, performing parameter identification on the damage constitutive model by using material mechanical property test results, comprising: obtaining stress-strain data of a material in different states through a type of material mechanical property test; determining test values of material mechanical property related parameters based on the stress-strain data; identifying and optimizing initial material parameters of the damage constitutive model by using the test values of the related parameters to determine optimized material parameters; generating a visco-plastic damage constitutive model of a copper alloy based on the optimized material parameters.

4. The method of claim 3, wherein, identifying and optimizing the initial material parameters of the damage constitutive model by using the test values of the related parameters, comprising: performing global search on simulation results of the damage constitutive model by using a random search algorithm to determine an initial parameter vector; performing fine search on the initial parameter vector by using a local optimization strategy to determine target material parameters.

5. The method of claim 3, wherein, identifying and optimizing the initial material parameters of the damage constitutive model by using the test values of the related parameters, comprising: performing simulation on a single material unit by using the damage constitutive model to determine first initial material parameters; the first initial material parameters are non-damage-related parameters; performing simulation on a full-size material sample model by using the damage constitutive model, and determining damage-related parameters of a target material according to the non-damage-related parameters and the test values of the related parameters.

6. The method of claim 1, wherein, performing structural analysis on a thrust chamber based on the damage constitutive model after parameter identification, comprising: determining temperature distribution results of an inner wall surface of the thrust chamber through three-dimensional flow-heat coupled simulation; determining a target cross section according to the temperature distribution results, and performing heat transfer analysis on the target cross section; performing structural analysis on the thrust chamber by using the damage constitutive model after parameter identification in combination with results of the heat transfer analysis.

7. A structure analysis device based on a coupled microcracking closure effect damage model, characterized by, The method comprises the following steps: a model construction module configured to construct a damage constitutive model based on micro-crack closure effect; The damage constitutive model is constructed based on a micro-crack closure effect, including: determining a crack closure effect coefficient according to a cross-sectional area parameter and a tension-compression state of the material; determining a state potential function coupled with the micro-crack closure effect based on the crack closure effect coefficient and a basic state potential function; and modeling material behavior based on a state variable, the state potential function coupled with the micro-crack closure effect, and a dissipation potential function to generate the damage constitutive model; a parameter identification module configured to identify parameters of the damage constitutive model by using material mechanical performance test results; a structure analysis module configured to perform structure analysis on the thrust chamber based on the damage constitutive model after parameter identification.

8. An electronic device comprising a memory, a processor, the memory having stored therein a computer program executable on the processor, characterized in that, The processor executes the computer program to implement the steps of the method of any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that, The computer readable storage medium stores computer executable instructions, and when the computer executable instructions are called and executed by the processor, the computer executable instructions cause the processor to execute the method of any one of claims 1 to 6.

Images