CN117556551A - Lightweight design method, device, equipment and medium for engine combustion chamber shell - Google Patents

Lightweight design method, device, equipment and medium for engine combustion chamber shell Download PDF

Info

Publication number
CN117556551A
CN117556551A CN202410042685.XA CN202410042685A CN117556551A CN 117556551 A CN117556551 A CN 117556551A CN 202410042685 A CN202410042685 A CN 202410042685A CN 117556551 A CN117556551 A CN 117556551A
Authority
CN
China
Prior art keywords
combustion chamber
file
elite
external elite
chamber shell
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202410042685.XA
Other languages
Chinese (zh)
Other versions
CN117556551B (en
Inventor
武泽平
张生全
马帅超
张为华
高经纬
王东辉
闫野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202410042685.XA priority Critical patent/CN117556551B/en
Publication of CN117556551A publication Critical patent/CN117556551A/en
Application granted granted Critical
Publication of CN117556551B publication Critical patent/CN117556551B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Combined Controls Of Internal Combustion Engines (AREA)

Abstract

The application belongs to the technical field of engine design, and relates to a lightweight design method, device, equipment and medium for an engine combustion chamber shell. The method comprises the following steps: sampling by adopting an optimized Latin hypercube to obtain a plurality of initial sample points; setting fitting points of the inner shape surface of the sealing head of the combustion chamber shell according to initial sample points, taking coordinates of the fitting points as design variables, taking the number of the fitting points as dimensions of the design variables, taking the mass of the combustion chamber shell as an objective function, taking the stress or deformation of the sealing head of the combustion chamber shell as a constraint function, and establishing a constraint optimization model of the combustion chamber shell; performing high-precision simulation according to the initial sample points, establishing external elite files according to the initial sample points and corresponding high-precision simulation values, and setting the number of the external elite files; when the external elite file converges, outputting the current sample point as the optimal solution of the constraint optimization model of the combustion chamber shell. By adopting the method and the device, the optimal light fitting point can be obtained.

Description

Lightweight design method, device, equipment and medium for engine combustion chamber shell
Technical Field
The present application relates to the technical field of engine design, and in particular, to a lightweight design method, device, equipment and medium for an engine combustion chamber housing.
Background
The combustion chamber of a solid rocket engine is the main component of a missile or rocket and is the means for storing and providing combustion of a propellant charge. Among them, the combustion chamber housing is an important part of the rocket engine, and is mainly used for fixing the combustion state inside the combustion chamber and resisting the impact and erosion of high-temperature and high-pressure gas. The typical combustor casing consists of a cylinder, a front end enclosure, a rear end enclosure, a front skirt, a rear skirt, a front joint and a rear joint. In order to meet the high-index design requirements of rocket engines on the combustion chamber shell, the lightweight design of the combustion chamber shell is an important task of technical implementation.
In the prior art, the common lightweight design method of the shell closure head of the combustion chamber comprises the following steps:
(1) The design method for establishing the geometric model comprises the following steps: the existing geometric model for optimizing the quality of the shell seal head of the combustion chamber has the advantages that the inner seal head curve and the outer seal head curve are conventional curve equations, such as ellipses, circles and the like.
(2) The design method based on experience and CAE (computer aided engineering) comprises the following steps: the engineer relies on the past engineering design experience to design the primary combustion chamber shell structure scheme, and then combines CAE simulation analysis and design indexes to adjust the design of the structural surface of the shell closure head.
(3) The design method based on the intelligent optimization algorithm comprises the following steps: and an intelligent optimization algorithm is utilized to construct iterative search optimization under constraint restriction, and an engineer is not required to have engineering structure design experience.
(4) The agent model-based optimization algorithm comprises the following implementation steps: establishing an optimal design model: the method comprises the steps of defining design variables, optimization targets and constraint variables of optimization problems, establishing a parameterized model of structural design of the combustion chamber shell, taking curves of inner and outer surfaces and transition curves as the design variables, and once the curves of the inner and outer surfaces and the transition section are determined, only determining the structure of the combustion chamber shell, selecting the parameters as the constraint variables according to stress indexes or deformation indexes provided by the whole body, wherein the optimization targets are the mass of the combustion chamber shell; constructing a proxy model of the target: constructing a target model and a proxy model of a constraint model by using sample points of experimental design, wherein the common proxy model methods include a radial basis method, a kriging method, a polynomial method and the like; sampling and iterative optimization: searching the agent model through the evolution algorithm, finding out sample points meeting the sampling criteria, carrying out high-precision simulation, stopping if the algorithm converges, and otherwise, updating the agent model to enter the next sampling.
However, the current lightweight design methods of these combustor casing heads have the following drawbacks:
(1) The design method for establishing the geometric model comprises the following steps: this geometric modeling approach is limited by the curve equation, resulting in a lightweight design result with a significant amount of room for improvement.
(2) The design method based on experience and CAE (computer aided engineering) comprises the following steps: first, this approach requires a certain amount of engineering experience in designing the combustor casing, and the degree of optimization is not well understood by a typical novice engineer. Secondly, the lightweight design by the method needs to be repeatedly combined with CAE simulation, and the simulation result of each time cannot be guaranteed to meet the design requirement, and continuous trial and error is needed, so that the whole design period consumes longer time, and the cost is increased. Finally, the result obtained by the lightweight design realized by the method can achieve a certain quality optimization effect, but is difficult to achieve the optimal effect.
(3) The design method based on the intelligent optimization algorithm comprises the following steps: the problem of lightweight design of the combustion chamber shell can be effectively solved, but thousands of iterations are usually needed, and the cost is high.
(4) Agent model-based optimization algorithm: the optimization efficiency can be effectively improved, but for the optimization problem with constraint, a simple punishment function method is mostly adopted, so that convergence accuracy is difficult to ensure.
Disclosure of Invention
Based on the above, it is necessary to provide a lightweight design method, device, equipment and medium for an engine combustion chamber shell, which can obtain an optimal lightweight fitting point, greatly improve the optimization efficiency and the optimization precision, and meanwhile, the method is not limited by a curve equation, does not need related design experience, shortens the period and reduces the cost.
A lightweight design method for an engine combustion chamber shell comprises the following steps:
sampling by adopting an optimized Latin hypercube to obtain a plurality of initial sample points;
setting fitting points of the inner shape surface of the sealing head of the combustion chamber shell according to initial sample points, taking coordinates of the fitting points as design variables, taking the number of the fitting points as dimensions of the design variables, taking the mass of the combustion chamber shell as an objective function, taking the stress or deformation of the sealing head of the combustion chamber shell as a constraint function, and establishing a constraint optimization model of the combustion chamber shell;
performing high-precision simulation according to the initial sample points, establishing external elite files according to the initial sample points and corresponding high-precision simulation values, and setting the number of the external elite files;
when the external elite file converges, outputting the current sample point as the optimal solution of the constraint optimization model of the combustion chamber shell.
In one embodiment, when the external elite file does not converge:
constructing an initial proxy model of the combustion shell quality optimization problem according to the initial sample points and the external elite file;
according to the initial proxy model, a three-stage constraint sampling method is adopted to obtain a plurality of next sample points;
acquiring a filling sample, optimizing an external elite file, and updating the external elite file;
when the updated external elite file is judged to be converged, outputting a current sample point;
and when the updated external elite file is judged not to be converged, updating the initial proxy model, regenerating the next sample point, updating the external elite file again, and judging until the updated external elite file is converged.
In one embodiment, building a constrained optimization model of a combustor casing includes:
in the method, in the process of the invention,for the purpose of +.>For design variables +.>For the dimension of the design variable +.>In order to constrain the function of the signal,is the number of constraint functions.
In one embodiment, setting the number of external elite files includes:
in the method, in the process of the invention,for the number of external elite files,min order to design the dimensions of the variables,Nis the number of initial sample points.
In one embodiment, constructing an initial proxy model of the combustion casing quality optimization problem from the initial sample points and the external elite archive, comprises:
in the method, in the process of the invention,for the initial proxy model, ++>For design variables +.>Is the coefficient of the basis function>As a function of the basis function,is the distance from the sample point to the center point.
In one embodiment, according to an initial proxy model, a three-stage constraint sampling method is adopted to obtain a plurality of next sample points, including:
in the method, in the process of the invention,obtaining the next sample point for the processing of the initial proxy model and the existing sample points, +.>Is->An initial proxy model after sub-sampling, +.>Is a certain sample point.
In one embodiment, obtaining a fill sample, optimizing an external elite archive, and updating the external elite archive includes:
when all external elite files are not feasible, the feasibility criteria for optimizing the external elite files with the filling sample meet:
in the method, in the process of the invention,normalized constraint violation for the current sample point, +.>For minimum normalized constraint violation in external elite archive,/->Is an external elite file;
when part of the external elite file is not feasible, the feasibility criterion for optimizing the external elite file by using the filling sample meets the following conditions:
in the method, in the process of the invention,for a certain sample pointXIs>For the optimal target value in the external elite file, +.>Normalized constraint violation for the current sample point, +.>Maximum constraint violation for a single constraint in the population;
when all external elite files are viable, the feasibility criteria for optimizing the external elite files with the fill samples satisfy:
according to the feasibility criteria, if the filling sample is better than the worst sample in the external elite file, the filling sample is used for replacing the worst sample in the external elite file, otherwise the external elite file is kept unchanged.
A lightweight design device for an engine combustion chamber housing, comprising:
the sampling module is used for sampling by adopting an optimized Latin hypercube to obtain a plurality of initial sample points;
the modeling module is used for setting fitting points of an inner shape surface of the sealing head of the combustion chamber shell according to the initial sample points, taking coordinates of the fitting points as design variables, taking the number of the fitting points as dimensions of the design variables, taking the mass of the combustion chamber shell as an objective function, and taking the stress or deformation of the sealing head of the combustion chamber shell as a constraint function to establish a constraint optimization model of the combustion chamber shell;
the file module is used for carrying out high-precision simulation according to the initial sample points, establishing external elite files according to the initial sample points and corresponding high-precision simulation values, and setting the number of the external elite files;
and the output module is used for outputting the current sample point as the optimal solution of the constraint optimization model of the combustion chamber shell when the external elite file converges.
A computer device comprising a memory storing a computer program and a processor which when executing the computer program performs the steps of:
sampling by adopting an optimized Latin hypercube to obtain a plurality of initial sample points;
setting fitting points of the inner shape surface of the sealing head of the combustion chamber shell according to initial sample points, taking coordinates of the fitting points as design variables, taking the number of the fitting points as dimensions of the design variables, taking the mass of the combustion chamber shell as an objective function, taking the stress or deformation of the sealing head of the combustion chamber shell as a constraint function, and establishing a constraint optimization model of the combustion chamber shell;
performing high-precision simulation according to the initial sample points, establishing external elite files according to the initial sample points and corresponding high-precision simulation values, and setting the number of the external elite files;
when the external elite file converges, outputting the current sample point as the optimal solution of the constraint optimization model of the combustion chamber shell.
A computer readable storage medium having stored thereon a computer program which when executed by a processor performs the steps of:
sampling by adopting an optimized Latin hypercube to obtain a plurality of initial sample points;
setting fitting points of the inner shape surface of the sealing head of the combustion chamber shell according to initial sample points, taking coordinates of the fitting points as design variables, taking the number of the fitting points as dimensions of the design variables, taking the mass of the combustion chamber shell as an objective function, taking the stress or deformation of the sealing head of the combustion chamber shell as a constraint function, and establishing a constraint optimization model of the combustion chamber shell;
performing high-precision simulation according to the initial sample points, establishing external elite files according to the initial sample points and corresponding high-precision simulation values, and setting the number of the external elite files;
when the external elite file converges, outputting the current sample point as the optimal solution of the constraint optimization model of the combustion chamber shell.
According to the lightweight design method, device, equipment and medium for the engine combustion chamber shell, firstly, parameterization modeling is conducted on the combustion chamber shell, fitting points are used for driving and fitting to generate an internal surface curve, and coordinates of the fitting points are selected as design variables. Then, experimental design is carried out on the design variables of the shell head of the combustion chamber of the rocket engine by adopting an optimized surging Ding Lifang design method, a group of initial sample points are obtained, and an initial proxy model is built based on the initial sample points. In the subsequent iteration process, constraint sampling is carried out by using an elite file enhanced constraint sampling method and the agent model is continuously updated, so that the accuracy of the agent model and the constraint model is continuously improved, and the accurate and efficient quality optimization design in the true sense is realized.
Drawings
FIG. 1 is an application scenario diagram of a lightweight design approach for an engine combustor casing in one embodiment;
FIG. 2 is a flow diagram of a method of lightweight design of an engine combustor casing in one embodiment;
FIG. 3 is a three-dimensional model of a combustor casing in one embodiment;
FIG. 4 is a two-dimensional model of a combustor casing in one embodiment;
FIG. 5 is a schematic frame diagram of a method of lightweight design of an engine combustor casing in one embodiment;
FIG. 6 is a diagram of iterative data of an fit point generation inlier scheme in one embodiment;
FIG. 7 is a stress cloud at one perspective in one embodiment;
FIG. 8 is a stress cloud at another perspective in one embodiment;
FIG. 9 is a graph of intermediate wall thickness plan iteration data for one embodiment;
FIG. 10 is a stress cloud of a cross-section of an embodiment;
FIG. 11 is a stress cloud of another cross-section in one embodiment;
FIG. 12 is a block diagram of a lightweight design device for an engine combustion chamber housing in one embodiment;
fig. 13 is an internal structural view of a computer device in one embodiment.
Reference numerals:
an end socket A and a cylinder B;
1 external surface, 2 internal surface, 3 end socket section, and 4 transition section.
Detailed Description
In order to make the objects, technical solutions and advantages of the present application more apparent, the present application will be further described in detail with reference to the accompanying drawings and examples. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
It should be noted that all directional indicators (such as up, down, left, right, front, and rear … …) in the embodiments of the present application are merely used to explain the relative positional relationship, movement, etc. between the components in a specific posture (as shown in the drawings), and if the specific posture is changed, the directional indicator is correspondingly changed.
In addition, descriptions such as those related to "first," "second," and the like, are provided for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated in this application. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include at least one such feature. In the description of the present application, "plurality of sets" means at least two sets, e.g., two sets, three sets, etc., unless specifically defined otherwise.
In the present application, unless explicitly specified and limited otherwise, the terms "coupled," "secured," and the like are to be construed broadly, and for example, "secured" may be either permanently attached or removably attached, or integrally formed; the device can be mechanically connected, electrically connected, physically connected or wirelessly connected; either directly or indirectly, through intermediaries, or both, may be in communication with each other or in interaction with each other, unless expressly defined otherwise. The specific meaning of the terms in this application will be understood by those of ordinary skill in the art as the case may be.
In addition, the technical solutions of the embodiments of the present application may be combined with each other, but it is necessary to be based on the fact that those skilled in the art can implement the technical solutions, and when the technical solutions are contradictory or cannot be implemented, the combination of the technical solutions should be considered to be absent, and is not within the scope of protection claimed in the present application.
The lightweight design method of the engine combustion chamber shell can be applied to an application scene graph shown in fig. 1. The terminal 102 communicates with the server 104 through a network, where the terminal 102 may include, but is not limited to, various personal computers, notebook computers, smartphones, tablet computers, and portable wearable devices, and the server 104 may be various portal sites, servers corresponding to a background of a working system, and the like.
The application provides a lightweight design method for an engine combustion chamber shell, as shown in fig. 2, in an embodiment, taking a terminal in fig. 1 as an example, the method is applied to the description, and the method includes:
and 202, sampling by adopting an optimized Latin hypercube to obtain a plurality of initial sample points.
In this step, how to sample and generate the initial sample points in the design space belongs to the prior art, and is not described herein.
And 204, setting fitting points of an inner shape surface of the sealing head of the combustion chamber shell according to the initial sample points, taking coordinates of the fitting points as design variables, taking the number of the fitting points as dimensions of the design variables, taking the mass of the combustion chamber shell as an objective function, and taking the stress or deformation of the sealing head of the combustion chamber shell as a constraint function, so as to establish a constraint optimization model of the combustion chamber shell.
In this step, as shown in fig. 3, a three-dimensional model diagram of a combustion chamber housing including: cylinder B and seal head A at the front end.
Due to the density of the combustion chamber housing material once selected (in the case of a single material)ρThe mass of the combustion chamber housing is determined as a function of the volume V, the size of which is in turn determined in particular by the geometric parameters, so that, depending on the geometric parameters of the modeling of the combustion chamber housing, the design variables of the optimization model can be determined, the dimensions of which are determined by the geometric parameters specifically selected.
In one embodiment, under the condition that the curve equation of the profile surface of the shell head of the combustion chamber is determined, an elite file enhancement optimization algorithm is utilized to drive coordinate data of an inner profile fitting point, and an optimal inner profile curve is generated through fitting.
As shown in the two-dimensional model diagram of the combustion chamber shell in fig. 4, the structure of the seal head is determined by the curve of the inner shape surface 2 and the curve of the outer shape surface 1,is fit point (i.e.)>,/>Comprising the following steps: fitting point of end socket section 3 and fitting point of transition section 4, < ->Representing the number of fitting points of the end socket segment>Representing the number of fitting points of the transition section>For the length of the outline vertex from the origin of coordinates,is the radius of revolution of the cylinder +.>Is the wall thickness of the cylinder. Selecting fitting point of end socket section->And transition end pointIs->Coordinates and transition piece->Is->The coordinates are used as design variables of the optimization model, and the dimension of the design variables is the total number of fitting points +.>. The fitting points of the inner shape surfaces of the seal heads are shown in table 1.
TABLE 1 fitting points of inner shape surfaces of seal heads
Specifically, establishing a constrained optimization model of the combustor casing includes:
in the method, in the process of the invention,for the objective function, i.e. the mass of the combustion chamber housing (combustion chamber housing structural mass), +.>For designing variables, i.e. fitting pointsPA kind of electronic devicexOr alternativelyyCoordinates of->Fitting points for the dimensions of the design variablesPTotal number of->For constraint functions, i.e. stress or deformation>Is the number of constraint functions.
Step 206, performing high-precision simulation according to the initial sample points, establishing external elite files according to the initial sample points and corresponding high-precision simulation values, and setting the number of the external elite files.
In this step, the concept of an external elite archive is introduced for storing several optimal individuals up to nowWherein->Is->Sample spots->For corresponding high-precision simulation (i.e., high-fidelity simulation) values. Through continuous updating in the algorithm iteration process, the performance of the individuals in the external elite file is always kept superior to that of other individuals, and the number of the individuals in the external elite file is always kept unchanged in the iteration and is the same as the number of the initial experimental design.
Establishing an external elite file and setting the number of the external elite files, wherein the method comprises the following steps:
in the method, in the process of the invention,for the number of external elite files, +.>For the dimension of the design variable +.>Is the number of initial sample points.
Step 208, when the external elite file converges, outputting the optimal solution of the constraint optimization model with the current sample point as the combustion chamber shell.
In this step, when the external elite file does not converge:
constructing an initial proxy model of the combustion shell quality optimization problem according to the initial sample points and the external elite file;
according to the initial proxy model, a three-stage constraint sampling method is adopted to obtain a plurality of next sample points;
acquiring a filling sample, optimizing an external elite file, and updating the external elite file;
when the updated external elite file is judged to be converged, outputting a current sample point;
and when the updated external elite file is judged not to be converged, updating the initial proxy model, regenerating the next sample point, updating the external elite file again, and judging until the updated external elite file is converged.
In the present embodiment, specifically:
according to the initial sample points and the external elite file, an initial proxy model of the combustion shell quality optimization problem is constructed by adopting a radial basis interpolation method, and the method comprises the following steps:
in the method, in the process of the invention,for the initial proxy model, ++>For design variables +.>Is the coefficient of the basis function>As basis functions, usually taking Gaussian basis functions, < ->Is the distance from the sample point to the center point.
According to the initial proxy model, a three-stage constraint sampling method is adopted to obtain a plurality of next sample points, including:
in the method, in the process of the invention,obtaining the next sample point for the processing of the initial proxy model and the existing sample points, +.>Is->An initial proxy model after sub-sampling, +.>Is a certain sample point.
Obtaining a filling sample according to an optimization algorithm (the algorithm is in the prior art), optimizing the external elite file, and updating the external elite file, including:
when all external elite files are not viable (i.e., there is at least one strong constraint present, the purpose of the next fill sample is to explore the boundaries of the constraint and find a viable region), the feasibility criteria for optimizing the external elite file with the fill sample satisfies:
in the method, in the process of the invention,normalized constraint violation for the current sample point, +.>For minimum normalized constraint violation in external elite archive,/->Is an external elite file;
when part of the external elite file is not feasible (i.e. the purpose of filling is to improve the target performance and to continuously approach the constraint boundary), the feasibility criterion for optimizing the external elite file with the filling sample is satisfied:
in the method, in the process of the invention,for a certain sample point->Is>For the optimal target value in the external elite file, +.>Normalized constraint violation for the current sample point, +.>Maximum constraint violation for a single constraint in the population;
when all external elite files are feasible (i.e., find globally optimal solution), the feasibility criteria for optimizing the external elite files with the fill samples satisfy:
according to the feasibility criterion, if the filling sample is better than the worst sample in the external elite file, the filling sample is used for replacing the worst sample in the external elite file, so that the external elite file is updated, otherwise, the external elite file is kept unchanged.
It is necessary to explain that: the convergence of the external elite file means that the preset iteration times are reached.
As shown in a frame schematic diagram of a lightweight design method of an engine combustion chamber shell in fig. 5, firstly, parameterizing and modeling the combustion chamber shell to obtain a parameterized geometric model, generating an inner shape surface curve by using fitting point driving fitting, and selecting coordinates of fitting points of the inner shape surface as design variables. Then, experimental design is carried out on design variables of a shell head of a combustion chamber of the rocket engine by adopting an optimized Latin hypercube design method, a group of initial sample points are obtained, high-fidelity simulation aiming at target constraint is carried out based on the initial sample points, and the number of elite files is set. And when judging that the elite file is not converged, selecting the elite file and constructing an initial proxy model. In the subsequent iteration process, the next sample point is selected based on a three-stage sampling method, the elite file is updated, constraint sampling is carried out by using a elite file enhanced constraint sampling method, and training samples are continuously updated to update the proxy model, so that the precision of the proxy model and the constraint model is continuously improved, and the truly accurate and efficient quality optimization design is realized.
According to the lightweight design method of the engine combustion chamber shell, aiming at the problems of lightweight, precise quality and efficient optimization of a solid rocket engine combustion chamber shell structure, a geometric model building method for generating a combustion chamber shell head inner shape surface curve by multiple fitting points is provided, an optimization model taking the coordinates of the fitting points of the combustion chamber shell head inner shape surface as design variables is built, the design target is that the mass of the combustion chamber shell is minimum under the condition of meeting constraint conditions through optimization design, a modeling method for generating the inner shape surface curve by the fitting points and an elite file enhancement constraint sampling method are adopted for searching, and an end socket shell structure quality optimization method based on elite file enhancement is designed. The combustion chamber shell is analyzed by the CAE method, and the final result shows that the combustion chamber shell lightweight design method can effectively reduce the quality of the combustion chamber shell while ensuring the strength of the combustion chamber shell, thereby achieving the original lightweight design purpose.
The method is based on a proxy model technology, and provides a combustor shell lightweight optimization modeling method for generating an inner shape surface based on data point fitting, wherein a certain number of fitting points are combined and arranged, different head inner shape surface spline curves are generated through fitting, different head inner shape surface curves are generated through changing coordinate data fitting of the fitting points, structural limitation of the inner shape surface and the outer shape surface head formed by a conventional single curve equation can be broken through, the method has the advantages of high freedom degree and strong universality, and lightweight design can be performed to the greatest extent; by providing the elite file enhanced constraint sampling method based on the inaccurate search idea, the feasibility and the optimality of balanced optimization are improved on the premise that the optimization result meets constraint condition setting, the algorithm efficiency is improved, the true feasible global optimal solution is positioned, the efficient design of quality optimization is realized, the design performance is improved while the simulation times are reduced, the overall quality of the combustion chamber shell tends to an optimal state, and an accurate and efficient design method is provided for the lightweight design of the rocket engine combustion chamber shell. The method can break through the limitation of the conventional inner and outer seal head curve parameter equation, and optimize the structural quality of the seal head of the combustion chamber to the greatest extent; when the dimension of the design variable is fixed, structural quality optimization efficiency of the shell head of the combustion chamber can be greatly improved based on elite file enhancement; the design variables and the quantity of the variables can be freely selected, so that the combustor shell has the characteristics of large freedom degree and good universality in the light design direction.
In one specific embodiment, the head is fitted from +12 points of the oval profile to create an internally contoured combustor casing, with the design variables selected as shown in Table 2. The curve of the external surface is taken as an elliptic equation, the internal surface is generated by 12 points in a fitting way (8 points of the seal head and 4 points of the transition section),,/>,/>for example, the fitting point of the closure head section is +.>End of transition piece->Is->Coordinates as design variables, individual fitting points +.>The coordinates are set to a certain different constant, specifically chosen +.>The values are related to the profile parameter equation, i.e. +.>Point of transition->To->Coordinates as design variables, each of which +.>The coordinates are set to a different constant, specifically chosen +.>Value and->、/>、/>Related to equal parameters, i.e.P/>The design variable dimension is 12.
Table 2 design variable selection
Constraint that maximum stress is less than or equal to 1080Mpa, input is internal pressure of 50Mpa, initial sampleN40 were selected. The total number of iterations was 460, as shown in fig. 6, which had begun to converge at 400 f, the final optimized mass was 3.268kg, and the coordinates of the fitting points for each internal surface were shown in table 3. Based on the upper partThe finite element analysis stress cloud diagram of the three-dimensional model of the combustion chamber shell, which is generated by the method, under the input of 50Mpa internal pressure is shown in fig. 7 and 8, and according to the information of fig. 7 and 8, the maximum stress of the inner surface and the outer surface is 1080Mpa, and red areas of the inner surface and the outer surface are overlapped, and almost the whole seal head structure is contained. It can be concluded that the combustion chamber housing head reaches a nearly constant stiffness, i.e. mass optimized, condition.
TABLE 3 coordinates of fitting points for each inner surface
The profile curve is used as an elliptic equation,,/>,/>the design of the combustor shell requires the end socket to be designed into a structure with equal wall thickness, and an iteration data diagram of the scheme with equal wall thickness calculated by using a particle swarm algorithm is shown in fig. 9. The optimal mass was 3.779kg, and convergence began at 70 iterations. The stress cloud diagrams of finite element analysis under the input of 50Mpa are shown in fig. 10 and 11, and according to the information of fig. 10 and 11, the maximum stress of the external surface does not reach 1080Mpa, only the internal surface reaches 1080Mpa in a certain rotation area, and the internal surface shows excessive redundancy of the rigidity of the top area of the seal head structure, so that the quality has a large optimization space.
The results of the two schemes are shown in Table 4, and the method has obvious advantages.
Table 4 comparison of the results of the two schemes
Compared with the prior art, the method has the advantages of less manual participation process, high design speed, capability of automatically processing constraint conditions in the optimization process, high design universality, capability of avoiding design optimization from separating from a feasible region, excellent performance of a design result, high design automation degree and high efficiency, and capability of effectively meeting the lightweight design of the combustion chamber shell.
It should be understood that, although the steps in the flowchart of fig. 2 are shown in sequence as indicated by the arrows, the steps are not necessarily performed in sequence as indicated by the arrows. The steps are not strictly limited to the order of execution unless explicitly recited herein, and the steps may be executed in other orders. Moreover, at least some of the steps in fig. 2 may include multiple sub-steps or stages that are not necessarily performed at the same time, but may be performed at different times, nor do the order in which the sub-steps or stages are performed necessarily performed in sequence, but may be performed alternately or alternately with at least a portion of the sub-steps or stages of other steps or other steps.
The present application also provides a lightweight design device for an engine combustion chamber housing, as shown in fig. 12, comprising, in one embodiment: a sampling module 1202, a modeling module 1204, an archive module 1206, and an output module 1208, wherein:
the sampling module 1202 is configured to sample with an optimized latin hypercube to obtain a plurality of initial sample points;
the modeling module 1204 is configured to set fitting points of an inner shape surface of the seal head of the combustor casing according to the initial sample points, set coordinates of the fitting points as design variables, set the number of the fitting points as dimensions of the design variables, set mass of the combustor casing as an objective function, and set stress or deformation of the seal head of the combustor casing as a constraint function, so as to establish a constraint optimization model of the combustor casing;
the archive module 1206 is configured to perform high-precision simulation according to the initial sample point, establish an external elite archive according to the initial sample point and the corresponding high-precision simulation value, and set the number of external elite archives;
the output module 1208 is configured to output the current sample point as an optimal solution of the constrained optimization model of the combustor casing when the external elite file converges.
For specific limitations on the lightweight design device of the engine combustion chamber housing, reference may be made to the above limitations on the lightweight design method of the engine combustion chamber housing, and no further description is given here. Each of the modules in the above-described apparatus may be implemented in whole or in part by software, hardware, and combinations thereof. The above modules may be embedded in hardware or may be independent of a processor in the computer device, or may be stored in software in a memory in the computer device, so that the processor may call and execute operations corresponding to the above modules.
In one embodiment, a computer device is provided, which may be a terminal, and the internal structure thereof may be as shown in fig. 13. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected by a system bus. Wherein the processor of the computer device is configured to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The network interface of the computer device is used for communicating with an external terminal through a network connection. The computer program, when executed by a processor, implements a method for lightweight design of an engine combustion chamber housing. The display screen of the computer equipment can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer equipment can be a touch layer covered on the display screen, can also be keys, a track ball or a touch pad arranged on the shell of the computer equipment, and can also be an external keyboard, a touch pad or a mouse and the like.
It will be appreciated by those skilled in the art that the structure shown in fig. 13 is merely a block diagram of a portion of the structure associated with the present application and is not limiting of the computer device to which the present application applies, and that a particular computer device may include more or fewer components than shown, or may combine some of the components, or have a different arrangement of components.
In an embodiment a computer device is provided comprising a memory storing a computer program and a processor implementing the steps of the method of the above embodiments when the computer program is executed.
In one embodiment, a computer readable storage medium is provided, on which a computer program is stored which, when executed by a processor, implements the steps of the method of the above embodiments.
Those skilled in the art will appreciate that implementing all or part of the above described methods may be accomplished by way of a computer program stored on a non-transitory computer readable storage medium, which when executed, may comprise the steps of the embodiments of the methods described above. Any reference to memory, storage, database, or other medium used in the various embodiments provided herein may include non-volatile and/or volatile memory. The nonvolatile memory can include Read Only Memory (ROM), programmable ROM (PROM), electrically Programmable ROM (EPROM), electrically Erasable Programmable ROM (EEPROM), or flash memory. Volatile memory can include Random Access Memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in a variety of forms such as Static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous Link DRAM (SLDRAM), memory bus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), among others.
The technical features of the above embodiments may be arbitrarily combined, and all possible combinations of the technical features in the above embodiments are not described for brevity of description, however, as long as there is no contradiction between the combinations of the technical features, they should be considered as the scope of the description.
The above examples only represent a few embodiments of the present application, which are described in more detail and are not to be construed as limiting the scope of the present application. It should be noted that it would be apparent to those skilled in the art that various modifications and improvements could be made without departing from the spirit of the present application, which would be within the scope of the present application. Accordingly, the scope of protection of the present application is to be determined by the claims appended hereto.

Claims (10)

1. The method for designing the engine combustion chamber shell in a lightweight way is characterized by comprising the following steps:
sampling by adopting an optimized Latin hypercube to obtain a plurality of initial sample points;
setting fitting points of the inner shape surface of the sealing head of the combustion chamber shell according to initial sample points, taking coordinates of the fitting points as design variables, taking the number of the fitting points as dimensions of the design variables, taking the mass of the combustion chamber shell as an objective function, taking the stress or deformation of the sealing head of the combustion chamber shell as a constraint function, and establishing a constraint optimization model of the combustion chamber shell;
performing high-precision simulation according to the initial sample points, establishing external elite files according to the initial sample points and corresponding high-precision simulation values, and setting the number of the external elite files;
when the external elite file converges, outputting the current sample point as the optimal solution of the constraint optimization model of the combustion chamber shell.
2. The method of claim 1, wherein when the external elite file does not converge:
constructing an initial proxy model of the combustion shell quality optimization problem according to the initial sample points and the external elite file;
according to the initial proxy model, a three-stage constraint sampling method is adopted to obtain a plurality of next sample points;
acquiring a filling sample, optimizing an external elite file, and updating the external elite file;
when the updated external elite file is judged to be converged, outputting a current sample point;
and when the updated external elite file is judged not to be converged, updating the initial proxy model, regenerating the next sample point, updating the external elite file again, and judging until the updated external elite file is converged.
3. The method of lightweight design of an engine combustor casing of claim 2, wherein building a constrained optimization model of the combustor casing comprises:
in the method, in the process of the invention,for the purpose of +.>For design variables +.>For the dimension of the design variable +.>For constraint function, ++>Is the number of constraint functions.
4. The method of designing a lightweight engine combustor casing according to claim 3, wherein setting the number of external elite files comprises:
in the method, in the process of the invention,for the number of external elite files,min order to design the dimensions of the variables,Nis the number of initial sample points.
5. The method of designing a combustion chamber housing for an engine according to any one of claims 2 to 4, wherein constructing an initial proxy model of a combustion chamber housing quality optimization problem from the initial sample points and the external elite file comprises:
in the method, in the process of the invention,for the initial proxy model, ++>For design variables +.>Is the coefficient of the basis function>As a basis function +.>Is the distance from the sample point to the center point.
6. The method for lightweight design of engine combustor casing according to claim 5, wherein the three-stage constrained sampling method is adopted according to an initial proxy model to obtain a plurality of next sample points, comprising:
in the method, in the process of the invention,obtaining the next sample point for the processing of the initial proxy model and the existing sample points, +.>Is->An initial proxy model after sub-sampling, +.>Is a certain sample point.
7. The method of claim 6, wherein obtaining a fill sample, optimizing an external elite file, and updating the external elite file comprises:
when all external elite files are not feasible, the feasibility criteria for optimizing the external elite files with the filling sample meet:
in the method, in the process of the invention,normalized constraint violation for the current sample point, +.>For minimum normalized constraint violation in external elite archive,/->Is an external elite file;
when part of the external elite file is not feasible, the feasibility criterion for optimizing the external elite file by using the filling sample meets the following conditions:
in the method, in the process of the invention,for a certain sample pointXIs>For the optimal target value in the external elite file, +.>Normalized constraint violation for the current sample point, +.>Maximum constraint violation for a single constraint in the population;
when all external elite files are viable, the feasibility criteria for optimizing the external elite files with the fill samples satisfy:
according to the feasibility criteria, if the filling sample is better than the worst sample in the external elite file, the filling sample is used for replacing the worst sample in the external elite file, otherwise the external elite file is kept unchanged.
8. The lightweight design device of engine combustion chamber casing, characterized by comprising:
the sampling module is used for sampling by adopting an optimized Latin hypercube to obtain a plurality of initial sample points;
the modeling module is used for setting fitting points of an inner shape surface of the sealing head of the combustion chamber shell according to the initial sample points, taking coordinates of the fitting points as design variables, taking the number of the fitting points as dimensions of the design variables, taking the mass of the combustion chamber shell as an objective function, and taking the stress or deformation of the sealing head of the combustion chamber shell as a constraint function to establish a constraint optimization model of the combustion chamber shell;
the file module is used for carrying out high-precision simulation according to the initial sample points, establishing external elite files according to the initial sample points and corresponding high-precision simulation values, and setting the number of the external elite files;
and the output module is used for outputting the current sample point as the optimal solution of the constraint optimization model of the combustion chamber shell when the external elite file converges.
9. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor implements the steps of the method of any of claims 1 to 7 when the computer program is executed.
10. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method of any of claims 1 to 7.
CN202410042685.XA 2024-01-11 2024-01-11 Lightweight design method, device, equipment and medium for engine combustion chamber shell Active CN117556551B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202410042685.XA CN117556551B (en) 2024-01-11 2024-01-11 Lightweight design method, device, equipment and medium for engine combustion chamber shell

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202410042685.XA CN117556551B (en) 2024-01-11 2024-01-11 Lightweight design method, device, equipment and medium for engine combustion chamber shell

Publications (2)

Publication Number Publication Date
CN117556551A true CN117556551A (en) 2024-02-13
CN117556551B CN117556551B (en) 2024-04-09

Family

ID=89811473

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202410042685.XA Active CN117556551B (en) 2024-01-11 2024-01-11 Lightweight design method, device, equipment and medium for engine combustion chamber shell

Country Status (1)

Country Link
CN (1) CN117556551B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117932824A (en) * 2024-03-22 2024-04-26 北京理工大学 Rotary inclusion type engine combustion chamber configuration parameterization design method

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106096127A (en) * 2016-06-07 2016-11-09 浙江大学 Robust error estimator method containing interval parameter uncertainty structure
CN109759897A (en) * 2019-03-14 2019-05-17 西安航天动力机械有限公司 A kind of measurement of the horizontal assembly of large size shell and aligning method
CN112434448A (en) * 2021-01-27 2021-03-02 中国人民解放军国防科技大学 Proxy model constraint optimization method and device based on multipoint adding
CN112528441A (en) * 2021-02-18 2021-03-19 中国人民解放军国防科技大学 Throat-plug type variable thrust engine overall parameter design method, device and equipment
CN113343545A (en) * 2021-08-02 2021-09-03 中国人民解放军国防科技大学 Structure adaptive optimization design method, device, equipment and medium
CN113761645A (en) * 2021-05-26 2021-12-07 青岛海洋地质研究所 Method, device and equipment for designing underwater vehicle shell
CN113806991A (en) * 2021-11-17 2021-12-17 天津仁爱学院 Engine combustion noise optimization prediction method and device and storage medium
CN115730404A (en) * 2022-11-25 2023-03-03 陕西法士特齿轮有限责任公司 Method, system, equipment and medium for designing molded surface of outer ring of one-way clutch
WO2023045244A1 (en) * 2021-09-26 2023-03-30 中国华能集团清洁能源技术研究院有限公司 Offshore wind turbine support structure optimization design method and system based on proxy model
CN116541963A (en) * 2023-06-02 2023-08-04 合肥工业大学 Solid rocket engine shell modeling method and system
US20230274047A1 (en) * 2022-05-06 2023-08-31 Nanjing University Of Aeronautics And Astronautics Method for analyzing fuselage profile based on measurement data of whole aircraft
CN116738583A (en) * 2023-08-16 2023-09-12 中国人民解放军国防科技大学 Solid rocket engine charging configuration constraint design method
CN117313250A (en) * 2023-11-30 2023-12-29 陕西空天信息技术有限公司 Design method, device, storage medium and equipment for solid rocket engine medicine column

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106096127A (en) * 2016-06-07 2016-11-09 浙江大学 Robust error estimator method containing interval parameter uncertainty structure
CN109759897A (en) * 2019-03-14 2019-05-17 西安航天动力机械有限公司 A kind of measurement of the horizontal assembly of large size shell and aligning method
CN112434448A (en) * 2021-01-27 2021-03-02 中国人民解放军国防科技大学 Proxy model constraint optimization method and device based on multipoint adding
CN112528441A (en) * 2021-02-18 2021-03-19 中国人民解放军国防科技大学 Throat-plug type variable thrust engine overall parameter design method, device and equipment
CN113761645A (en) * 2021-05-26 2021-12-07 青岛海洋地质研究所 Method, device and equipment for designing underwater vehicle shell
CN113343545A (en) * 2021-08-02 2021-09-03 中国人民解放军国防科技大学 Structure adaptive optimization design method, device, equipment and medium
WO2023045244A1 (en) * 2021-09-26 2023-03-30 中国华能集团清洁能源技术研究院有限公司 Offshore wind turbine support structure optimization design method and system based on proxy model
CN113806991A (en) * 2021-11-17 2021-12-17 天津仁爱学院 Engine combustion noise optimization prediction method and device and storage medium
US20230274047A1 (en) * 2022-05-06 2023-08-31 Nanjing University Of Aeronautics And Astronautics Method for analyzing fuselage profile based on measurement data of whole aircraft
CN115730404A (en) * 2022-11-25 2023-03-03 陕西法士特齿轮有限责任公司 Method, system, equipment and medium for designing molded surface of outer ring of one-way clutch
CN116541963A (en) * 2023-06-02 2023-08-04 合肥工业大学 Solid rocket engine shell modeling method and system
CN116738583A (en) * 2023-08-16 2023-09-12 中国人民解放军国防科技大学 Solid rocket engine charging configuration constraint design method
CN117313250A (en) * 2023-11-30 2023-12-29 陕西空天信息技术有限公司 Design method, device, storage medium and equipment for solid rocket engine medicine column

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
ZHIXIANG WANG等: "Lightweight design of cylindrical stiffened shells in launch vehicles by a dual-elite population sequential approximation optimization approach", ENGINEERING OPTIMIZATION, 9 June 2020 (2020-06-09), pages 1 - 22 *
彭莎莎 等: "薄壁大开口分段固体发动机卧式对接装配研究", 航天制造技术, no. 06, 25 December 2014 (2014-12-25), pages 11 - 15 *
杨家伟 等: "基于非精确搜索的固体火箭发动机序列近似约束优化设计方法", 《机械设计》, 6 December 2023 (2023-12-06), pages 1 - 6 *
武泽平 等: "应用径向基代理模型实现序列自适应再采样优化策略", 国防科技大学学报, no. 06, 28 December 2014 (2014-12-28), pages 18 - 24 *
胡宽 等: "基于网络理论的固体火箭发动机纤维缠绕壳体优化设计", 兵工学报, no. 09, 15 September 2008 (2008-09-15), pages 1099 - 1103 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117932824A (en) * 2024-03-22 2024-04-26 北京理工大学 Rotary inclusion type engine combustion chamber configuration parameterization design method
CN117932824B (en) * 2024-03-22 2024-06-07 北京理工大学 Rotary inclusion type engine combustion chamber configuration parameterization design method

Also Published As

Publication number Publication date
CN117556551B (en) 2024-04-09

Similar Documents

Publication Publication Date Title
Choi et al. Gradient-based constrained optimization using a database of linear reduced-order models
CN117556551B (en) Lightweight design method, device, equipment and medium for engine combustion chamber shell
US11734471B2 (en) Topology optimization for subtractive manufacturing techniques
Zhou et al. A two-stage adaptive multi-fidelity surrogate model-assisted multi-objective genetic algorithm for computationally expensive problems
US8248407B2 (en) System and method of computer-aided design of a 3D object modeled by geometries
CN101877021B (en) Method for automatically converting data among kinetic analysis, three-dimensional modeling and finite-element analysis software
JP5045853B2 (en) Calculation data generation apparatus, calculation data generation method, and calculation data generation program
Zhang et al. An improved two-stage framework of evidence-based design optimization
Talebi et al. Stress analysis of 3D complex geometries using the scaled boundary polyhedral finite elements
Tromme et al. Structural optimization of multibody system components described using level set techniques
Daxini et al. Parametric shape optimization techniques based on Meshless methods: A review
Chen et al. Tetrahedral mesh improvement by shell transformation
Koziel et al. Rapid multi-objective aerodynamic design using co-kriging and space mapping
Liu et al. A sequential sampling generation method for multi-fidelity model based on Voronoi region and sample density
Zheng et al. An improved local remeshing algorithm for moving boundary problems
CN113821877B (en) Hypersonic aerodynamic thermal numerical simulation-oriented regional decomposition calculation method and hypersonic aerodynamic thermal numerical simulation-oriented regional decomposition calculation device
Schneider et al. A new approach for force‐displacement co‐simulation using kinematic coupling constraints
WO2020023811A1 (en) 3d object design synthesis and optimization using existing designs
Ulu et al. Structural design using laplacian shells
CN112926153B (en) Continuous body structure fatigue topology optimization design method combining K-S function
Ren et al. Multi-fidelity aerodynamic shape optimization using manifold mapping
Meng et al. A nurbs-enhanced finite volume method for steady euler equations with goal-oriented h-adaptivity
US20080052048A1 (en) Analytical Model Creation System
JP5164808B2 (en) Analysis mesh generator
Chen et al. Surface triangular mesh and volume tetrahedral mesh generations for biomolecular modeling

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant