CN115859492B - Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine - Google Patents

Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine Download PDF

Info

Publication number
CN115859492B
CN115859492B CN202211161938.2A CN202211161938A CN115859492B CN 115859492 B CN115859492 B CN 115859492B CN 202211161938 A CN202211161938 A CN 202211161938A CN 115859492 B CN115859492 B CN 115859492B
Authority
CN
China
Prior art keywords
curve
valve
profile
valve rod
flow
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.)
Active
Application number
CN202211161938.2A
Other languages
Chinese (zh)
Other versions
CN115859492A (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 CN202211161938.2A priority Critical patent/CN115859492B/en
Publication of CN115859492A publication Critical patent/CN115859492A/en
Application granted granted Critical
Publication of CN115859492B publication Critical patent/CN115859492B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Lift Valve (AREA)

Abstract

The invention discloses an actual flow matching design method of a solid attitude and orbit control engine gas valve profile, which is based on a theoretical target flow curve, combines a binary search method, taylor local expansion and the like to recursively obtain a valve rod profile curve coordinate point set, corrects the theoretical target flow curve based on the actual maximum flow of a valve to obtain an actual target flow curve, adopts the same valve rod profile design method to obtain a corresponding valve rod profile coordinate point set, obtains a valve rod displacement-valve flow curve through computational fluid dynamics simulation, corrects the valve rod profile curve coordinate point set according to the error between the valve rod displacement-valve flow curve and the target flow curve, and circularly carries out computational fluid dynamics simulation-valve rod profile curve coordinate point set correction until the error meets the design requirement. The invention is applied to the field of solid attitude and orbit control engine design, can obviously improve the effectiveness of a design result, and ensures that the linearity of the valve rod displacement and the actual flow of the valve meets the design requirement.

Description

Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine
Technical Field
The invention relates to the technical field of solid attitude and orbit control engine design, in particular to a method for matching and designing the actual flow of a gas valve molded surface of a solid attitude and orbit control engine.
Background
The gas type solid attitude control engine is used as an actuating mechanism of the kinetic energy interceptor, has the characteristics of simple structure, safety, reliability and the like, and is particularly suitable for intercepting weapon systems with high safety requirements on space bases, sea bases and the like. Solid attitude control engine thrust adjustment is largely dependent on the valve stem geometry in the valve, the nozzle profile, and the relative position between the valve stem and the nozzle. If the flow of a single valve is nonlinear along with the stroke of the valve rods, the constant total flow of the valve is difficult to ensure in the motion process of each valve rod, so that the pressure of a combustion chamber can be changed along with the valve, and the difficulty of thrust adjustment of a solid attitude and orbit control engine and the complexity of a control system are greatly improved. Therefore, there is a need to develop valve profile optimization designs that target the high degree of linear variation of individual valve actual flow with valve stem displacement for a given nozzle profile, thereby simplifying control system complexity as much as possible.
The valve profile design method commonly used at present comprises the following steps:
valve profile design based on existing schemes: based on the valve profile design result which is relatively close to the design requirement, the valve profile is manually modified by a trial and error method to obtain a design scheme which can meet the requirement. The method is simple and efficient, is easy to realize, effectively avoids complex operation, but needs the participation of engineers with higher engineering project practice experience, and is difficult to ensure that the design result is an optimal scheme, and the linearity is difficult to meet the requirement.
Valve profile optimization design based on intelligent optimization algorithm: the size parameters of the spray pipe molded surface and the valve rod molded surface are used as design variables, the error between the valve rod displacement-valve equivalent throat area curve and the target curve is used as an objective function, and the valve molded surface optimal design is realized based on an intelligent optimization algorithm. The method can realize automatic optimal design of the valve profile, and the design result is accurate. However, a spray pipe and valve rod profile parameterized model needs to be constructed, profile control parameters need to be reasonably planned, the control on the profile shape is limited due to too few parameters, and a design result meeting linearity requirements is difficult to obtain; too many parameters can lead to high complexity of the optimization model, unacceptable time cost of the optimization design and great difficulty in application to engineering practice.
Disclosure of Invention
Aiming at the problems that in the prior art, the gas valve of the solid attitude and track control engine has obvious difference between the theoretical flow and the actual flow of the valve and the maximum theoretical flow of the valve cannot be achieved in the linearization process of the valve flow, and the like, the invention provides the actual flow matching design method of the gas valve profile of the solid attitude and track control engine, which can realize the valve profile design with higher efficiency, correct the difference between the theoretical flow and the actual flow of the valve, ensure the linearity of the valve rod displacement-the valve flow to meet the design requirement, and avoid the problems that the valve rod displacement-the theoretical flow curve has high linearity but the valve rod displacement-the valve flow curve cannot achieve the same linearity.
In order to achieve the purpose, the invention provides a method for matching and designing the actual flow of the molded surface of a gas valve of a solid attitude and orbit control engine, which comprises the following steps:
step 1, acquiring a coordinate point set of a spray pipe profile curve, and initializing a flow curve into a theoretical target flow curve;
step 2, determining a corresponding equivalent throat area curve based on the flow curve, and obtaining a valve rod profile curve coordinate point set meeting the equivalent throat area curve based on a recursion matching design method;
step 3, performing computational fluid dynamics simulation based on a spray pipe profile curve coordinate point set and a valve rod profile curve coordinate point set to obtain actual maximum flow of the valve, correcting a theoretical target flow curve according to the actual maximum flow of the valve to obtain an actual target flow curve, and updating the flow curve into the actual target flow curve;
step 4, obtaining an updated valve rod profile curve coordinate point set according to the updated flow curve, and then carrying out computational fluid dynamics simulation according to the spray pipe profile curve coordinate point set and the updated valve rod profile curve coordinate point set to obtain a valve rod displacement-valve flow curve;
step 5, judging whether the current valve rod displacement-valve flow curve meets the convergence condition:
if yes, outputting a current valve rod profile curve coordinate point set;
otherwise, updating the flow curve according to the error between the current valve rod displacement-valve flow curve and the actual target flow curve, and returning to the value step 4 for iteration.
Compared with the prior art, the invention has the advantages that:
1. compared with the valve profile design based on the existing scheme, the design result of the invention has higher accuracy;
2. compared with the valve profile optimization design based on the intelligent optimization algorithm, the method has the advantages that the design flow is simpler, the design efficiency is higher, and the linearity of the design result can be ensured under the condition of low calculation cost.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to the structures shown in these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a design method in an embodiment of the invention;
FIG. 2 is a schematic view of a coordinate point corresponding to a minimum throat area formed by a tangent line of a nozzle profile and a valve stem profile and a translated tangent line and the nozzle profile when a valve is completely closed in an embodiment of the present invention;
FIG. 3 is a schematic diagram of Taylor local expansion slope search in accordance with an embodiment of the present invention;
FIG. 4 is a schematic diagram of a gas valve profile of a solid attitude and orbit control engine in an example of an embodiment of the invention;
FIG. 5 is a schematic diagram showing the comparison of flow matching design results when theoretical target flow is not corrected in the example of the embodiment of the present invention;
FIG. 6 is a schematic diagram showing the comparison of flow matching design results after correcting the theoretical target flow in the example of the embodiment of the invention;
FIG. 7 is a comparison of valve stem profile design results for an example of an embodiment of the present invention;
FIG. 8 is a graph showing the comparison of flow matching design results for a valve stem with a relative stroke of 8mm in an example of an embodiment of the present invention;
FIG. 9 is a graphical representation of a comparison of valve stem profile design results for an example of an embodiment of the present invention with a valve stem relative travel of 8 mm;
the achievement of the objects, functional features and advantages of the present invention will be further described with reference to the accompanying drawings, in conjunction with the embodiments.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
It should be noted that all directional indicators (such as up, down, left, right, front, and rear … …) in the embodiments of the present invention are merely used to explain the relative positional relationship, movement, etc. between the components in a particular posture (as shown in the drawings), and if the particular posture is changed, the directional indicator is changed accordingly.
Furthermore, descriptions such as those referred to as "first," "second," and the like, are provided for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implying an order of magnitude of the indicated technical features in the present disclosure. 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 invention, the meaning of "plurality" means at least two, for example, two, three, etc., unless specifically defined otherwise.
In the present invention, unless specifically stated and limited otherwise, the terms "connected," "affixed," and the like are to be construed broadly, and for example, "affixed" may be a fixed connection, a removable connection, or an integral body; 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 above terms in the present invention can be understood by those of ordinary skill in the art according to the specific circumstances.
In addition, the technical solutions of the embodiments of the present invention 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 as not existing, and not falling within the scope of protection claimed by the present invention.
The embodiment discloses an actual flow matching design method of a solid attitude and orbit control engine gas valve profile, which is based on a theoretical target flow curve, combines a binary search method, a Taylor local expansion and other recursions to obtain a valve rod profile curve coordinate point set, corrects the theoretical target flow curve based on the actual maximum flow of a valve to obtain an actual target flow curve, adopts the same valve rod profile design method to obtain a corresponding valve rod profile coordinate point set, obtains a valve rod displacement-valve flow curve through computational fluid dynamics simulation, corrects the valve rod profile curve coordinate point set according to root mean square error between the valve flow curve and the target flow curve, and circularly carries out the computational fluid dynamics simulation-valve rod profile curve coordinate point set correction process until the error meets the design requirement. Compared with other implementation schemes, the method can obviously improve the effectiveness of the design result, and ensures that the linearity of the valve rod displacement and the valve actual flow meets the design requirement.
Referring to fig. 1, the design method for matching the actual flow of the gas valve profile of the solid attitude and orbit control engine in the embodiment specifically comprises the following steps 1-5.
Step 1, obtaining the throat diameter of a spray pipe, the radius of transitional circular arcs at the upstream of the throat, the convergence half angle and the working pressure
Figure 604704DEST_PATH_IMAGE001
Characteristic speed of propellant->
Figure 397080DEST_PATH_IMAGE002
Theoretical target flow curve->
Figure 597117DEST_PATH_IMAGE003
Convergence error satisfying design requirements +.>
Figure 324902DEST_PATH_IMAGE004
. According to known throat diameter and throat position of spray pipeThe curve coordinate point set of the spray pipe profile is determined by the radius of the transitional arc and the convergence half angle, and different curve coordinate point sets of the spray pipe profile can be obtained by changing the throat diameter of the spray pipe and/or the radius of the transitional arc upstream of the throat part and/or the convergence half angle, so that the design requirements of different solid attitude rail control engines are met. Initializing a flow curve->
Figure 938286DEST_PATH_IMAGE005
For a theoretical target flow curve, i.e.
Figure 144139DEST_PATH_IMAGE006
。
And 2, determining a corresponding equivalent throat area curve based on the flow curve, and obtaining a valve rod profile curve coordinate point set meeting the equivalent throat area curve based on a recursion matching design method.
In the specific implementation process, the flow formula of the one-dimensional isentropic flow of the ideal gas can be adopted
Figure 401945DEST_PATH_IMAGE007
Calculating to obtain a flow curve +.>
Figure 159685DEST_PATH_IMAGE008
Corresponding equivalent throat area curve->
Figure 732137DEST_PATH_IMAGE009
I.e. +.>
Figure 741681DEST_PATH_IMAGE010
。
In this embodiment, the specific process of obtaining the valve rod profile curve coordinate point set meeting the equivalent throat area curve based on the recursive matching design method is as follows:
step 2.1, calculating to obtain the tangent point of the spray pipe molded surface and the valve rod molded surface when the valve is completely closed based on a dichotomy
Figure 978627DEST_PATH_IMAGE011
The specific implementation mode is as follows:
step 2.1.1, obtaining the valve rod displacement in the equivalent throat area curve as
Figure 782635DEST_PATH_IMAGE012
Equivalent throat area at that time
Figure 105032DEST_PATH_IMAGE013
, wherein ,/>
Figure 918267DEST_PATH_IMAGE014
For the relative travel of the valve stem>
Figure 885086DEST_PATH_IMAGE015
Calculating step length;
step 2.1.2, searching the tangent point of the spray pipe and the valve rod when the valve is completely closed on the profile curve of the spray pipe by adopting a dichotomy method
Figure 719050DEST_PATH_IMAGE016
So that it satisfies the following equation:
Figure 528743DEST_PATH_IMAGE017
wherein ,
Figure 145669DEST_PATH_IMAGE018
to be at the spot->
Figure 826050DEST_PATH_IMAGE019
Straight line edge tangent to the nozzle>
Figure 237439DEST_PATH_IMAGE020
Negative axis translation +.>
Figure 409795DEST_PATH_IMAGE021
The minimum throat area formed by the rear and the spray pipe profile is that:
Figure 423887DEST_PATH_IMAGE022
wherein ,
Figure 365298DEST_PATH_IMAGE023
is the curve coordinate point set of the spray pipe profile>
Figure 806644DEST_PATH_IMAGE024
Respectively a valve rod profile curve coordinate point set, +.>
Figure 466295DEST_PATH_IMAGE025
Concentrate the curve coordinate point of the spray pipe profilenThe abscissa value of the dot, +.>
Figure 304586DEST_PATH_IMAGE026
Concentrate the curve coordinate point of the valve rod profilepThe abscissa value of the point;
step 2.2, the molded surface of the spray pipe is dotted
Figure 366083DEST_PATH_IMAGE027
The tangential line edge->
Figure 978330DEST_PATH_IMAGE028
Negative axis translation +.>
Figure 859699DEST_PATH_IMAGE029
The rear straight line is used as a valve rod partial profile curve, and a coordinate point corresponding to the minimum throat area formed by the spray pipe profile on the valve rod partial profile curve is obtained by solving>
Figure 215594DEST_PATH_IMAGE030
Point->
Figure 866018DEST_PATH_IMAGE031
And (4) point->
Figure 914745DEST_PATH_IMAGE032
Sequentially adding a valve rod profile curve coordinate point set and initializing recurrence algebra +.>
Figure 283410DEST_PATH_IMAGE033
Namely, as shown in FIG. 2, the transition arc in FIG. 2 is the transition arc of the curve of the spray pipe profile, and the arc tangent is the upper point of the spray pipe profile +.>
Figure 442996DEST_PATH_IMAGE034
The tangent line is located, and the local molded surface of the valve rod is a point on the molded surface of the spray pipe
Figure 947926DEST_PATH_IMAGE035
The tangential line edge->
Figure 308500DEST_PATH_IMAGE036
Negative axis translation +.>
Figure 289095DEST_PATH_IMAGE037
A rear straight line;
step 2.3, obtaining the valve rod displacement in the equivalent throat area curve as
Figure 127738DEST_PATH_IMAGE038
Equivalent throat area ∈>
Figure 877388DEST_PATH_IMAGE039
Searching for satisfying equivalent throat area based on dichotomy>
Figure 408863DEST_PATH_IMAGE040
Slope of valve stem profile local taylor expansion line
Figure 486541DEST_PATH_IMAGE041
And the equivalent throat area of the valve rod profile and the spray pipe profile on the local Taylor expansion straight line of the current valve rod profile is formed +.>
Figure 253509DEST_PATH_IMAGE042
Is a new point +.>
Figure 733031DEST_PATH_IMAGE043
And new point +.>
Figure 297392DEST_PATH_IMAGE044
Adding a valve rod profile curve coordinate point set, referring to fig. 3, the specific implementation mode is as follows:
step 2.3.1, initializing a slope search interval
Figure 127945DEST_PATH_IMAGE045
And obtain the valve rod displacement in the equivalent throat area curve as +.>
Figure 167445DEST_PATH_IMAGE046
Equivalent throat area ∈>
Figure 32633DEST_PATH_IMAGE047
, wherein ,/>
Figure 640332DEST_PATH_IMAGE048
For maximum value of slope search, +.>
Figure 82815DEST_PATH_IMAGE049
Minimum value for slope search; />
Step 2.3.2, let point
Figure 801372DEST_PATH_IMAGE050
For the local expansion point, search slope +.>
Figure 521066DEST_PATH_IMAGE051
And let the slope be +.>
Figure 424300DEST_PATH_IMAGE052
And go through the local expansion point->
Figure 963866DEST_PATH_IMAGE053
As a current local expansion line, wherein the slope +.>
Figure 876327DEST_PATH_IMAGE054
The method meets the following conditions: current officeThe minimum throat area formed by the straight line of expansion of the part and the spray pipe profile is equal to the target equivalent throat area +.>
Figure 184949DEST_PATH_IMAGE055
；
Step 2.3.3, the minimum throat area formed by the local Taylor expansion straight line of the current valve rod profile and the spray pipe profile is equal to the equivalent throat area
Figure 400030DEST_PATH_IMAGE056
Is a new point +.>
Figure 285946DEST_PATH_IMAGE057
And->
Figure 611885DEST_PATH_IMAGE058
Adding a valve rod profile curve coordinate point set;
step 2.4, convergence judgment: order the
Figure 165226DEST_PATH_IMAGE059
Judging->
Figure 285629DEST_PATH_IMAGE060
Whether or not it is:
if yes, enter step 2.5;
otherwise, returning to the step 2.3 to continue iteration;
step 2.5, adding the circle center as the coordinate point data of the valve rod profile according to the coordinate point set of the valve rod profile curve
Figure 924421DEST_PATH_IMAGE061
Radius is->
Figure 54051DEST_PATH_IMAGE062
Is a circular arc of (2), wherein>
Figure 71685DEST_PATH_IMAGE063
For the valve rod profile local Taylor at the last iteration of step 2.3Slope of the expansion line>
Figure 750273DEST_PATH_IMAGE064
、/>
Figure 751727DEST_PATH_IMAGE065
For the point at the last iteration of step 2.3
Figure 278523DEST_PATH_IMAGE066
Is a horizontal and vertical coordinate value of (a);
and 2.6, adding the coordinate point set on the circular arc in the step 2.5 into the valve rod profile curve coordinate point set, and outputting the valve rod profile curve coordinate point set.
In this embodiment, the slope search interval in step 2.3.1 is
Figure 681822DEST_PATH_IMAGE067
。
In the specific implementation process, the specific searching process of performing the binary search in the steps 2.1.2 and 2.3.2 is a conventional technical means in the field, so that a detailed description thereof is omitted in this embodiment.
Step 3, performing computational fluid dynamics simulation based on a spray pipe profile curve coordinate point set and a valve rod profile curve coordinate point set to obtain actual maximum flow of the valve, correcting a theoretical target flow curve according to the actual maximum flow of the valve to obtain an actual target flow curve, and updating the flow curve into the actual target flow curve, wherein the specific implementation process is as follows:
step 3.1, importing a spray pipe profile curve coordinate point set and a valve rod profile curve coordinate point set into grid division software to finish unstructured grid division and grid file output;
step 3.2, importing the grid file output by the grid division software into computational fluid dynamics simulation software, and obtaining the actual maximum flow of the valve by solving a two-dimensional axisymmetric conservation type Navier-Stokes equation of compressible flow
Figure 675186DEST_PATH_IMAGE068
；
And 3.3, correcting a theoretical target flow curve according to the actual maximum flow of the valve to obtain an actual target flow curve, wherein the actual target flow curve is specifically as follows:
Figure 757412DEST_PATH_IMAGE069
Figure 494424DEST_PATH_IMAGE070
in the formula ,
Figure 486650DEST_PATH_IMAGE071
for the actual target flow curve, +.>
Figure 244391DEST_PATH_IMAGE072
For the theoretical target flow curve, +.>
Figure 954858DEST_PATH_IMAGE073
For the actual maximum flow of the valve>
Figure 354615DEST_PATH_IMAGE074
The theoretical maximum flow of the valve;
and 3.4, updating the flow curve into an actual target flow curve.
Step 4, obtaining an updated valve rod profile curve coordinate point set according to the updated flow curve, and then carrying out computational fluid dynamics simulation according to the spray pipe profile curve coordinate point set and the updated valve rod profile curve coordinate point set to obtain a valve rod displacement-valve flow curve, wherein the specific implementation process is as follows:
step 4.1, obtaining an updated valve rod profile curve coordinate point set according to the updated flow curve by adopting the same method as the step 2;
step 4.2, importing the spray pipe profile curve coordinate point set and the updated valve rod profile curve coordinate point set into grid division software to complete unstructured grid division and new grid file output;
step 4.3, the new grid file output by the grid division software is imported into computational fluid dynamics simulation software, and a valve rod displacement-valve flow curve is obtained by solving a two-dimensional axisymmetric conservation type Navie-Stokes equation of compressible flow
Figure 466928DEST_PATH_IMAGE075
。
Step 5, judging whether the current valve rod displacement-valve flow curve meets the convergence condition, namely judging
Figure 536515DEST_PATH_IMAGE076
Whether or not it is, wherein->
Figure 858912DEST_PATH_IMAGE077
For the number of valve stem steps, < >>
Figure 672147DEST_PATH_IMAGE078
Is the first in the valve flow curveiValve flow corresponding to the displacement of the valve rod +.>
Figure 904545DEST_PATH_IMAGE079
Is the first in the actual target flow curveiValve flow corresponding to the displacement of the valve rod +.>
Figure 4088DEST_PATH_IMAGE080
The convergence error is:
if true, the current valve stem displacement-valve flow curve
Figure 689148DEST_PATH_IMAGE081
The corresponding valve rod profile curve coordinate point set is output as a design result;
if not, after updating the flow curve based on the error between the current valve rod displacement-valve flow curve and the actual target flow curve, returning to the value step 4 for iteration, wherein the updating process of the flow curve is as follows
Figure 430707DEST_PATH_IMAGE082
。
The invention is further described below with reference to specific examples.
Fig. 4 is a geometric configuration of a gas valve of a solid attitude and orbit control engine, and part of design parameters and working parameters are shown in table 1, and then a specific implementation scheme of valve rod profile design of the gas valve of the solid attitude and orbit control engine is as follows:
(1) A valve rod profile coordinate point set is obtained by adopting a binary search method and Taylor local expansion recursion;
1.1 Initializing a profile tangent point search interval to be a transition arc boundary, calculating a tangent point between a spray pipe profile and a valve rod profile based on a binary search method, and updating a valve rod profile coordinate point set;
1.2 Updating the existing valve rod profile coordinate point set based on the calculation step length;
1.3 Calculating the Taylor local expansion slope based on a binary search method, solving the equivalent throat area, and recursively obtaining the next local expansion point calculation point;
1.4 Convergence determination: if the target curve is calculated, the step 1.5) is switched to, otherwise, the step 1.2) is switched to continue the recursive calculation;
1.5 Increasing the arc of the head of the valve rod molded surface;
1.6 Outputting a valve rod profile curve coordinate point set corresponding to the target flow curve;
(2) Calling grid dividing software to output an unstructured grid file;
(3) Calling computational fluid dynamics simulation software to obtain the actual maximum flow of the valve;
(4) Correcting the theoretical target flow curve according to the actual maximum flow of the valve to obtain an actual target flow curve;
(5) Obtaining a valve rod profile curve and a spray pipe profile curve coordinate point set meeting the flow curve requirement according to the method of the step (1);
(6) Calling grid dividing software to output an unstructured grid file;
(7) Invoking computational fluid dynamics simulation software to obtain a valve rod displacement-valve flow curve;
(8) And (3) convergence judgment: outputting a valve rod profile curve and a spray pipe profile curve coordinate point set meeting linearity requirements if the root mean square error of the valve flow curve and the target flow curve meets design requirements, otherwise updating flow curve requirements according to the error of the flow curve and the actual target flow curve, and turning to the step (5) to continue iteration.
Table 1 design parameters and operating parameters of gas valve section of solid attitude and orbit control engine
Figure 517612DEST_PATH_IMAGE083
Based on the parameters, aiming at the state that the relative stroke of the valve rod is 4mm, the valve rod profile design under the condition that the theoretical target flow is not corrected and the theoretical target flow is corrected is developed. As shown in fig. 5, due to flow loss, the matched design flow obtained based on the recursive optimization method cannot completely coincide with the theoretical target flow, and the maximum theoretical flow of the valve cannot be achieved due to the flow loss, and the linearity of the gas valve near the maximum flow cannot be ensured if the theoretical target flow is not corrected;
the flow matching design result after correcting the theoretical target flow is shown in fig. 6, and the flow matching design result is basically identical with the target flow at the moment, meanwhile, because the difference exists between the theoretical flow and the actual flow of the valve, the high linear change of the valve flow along with the displacement of the valve rod cannot be truly realized based on the profile design result of the recursive optimization method;
the valve rod profile design result is shown in fig. 7, and the gas valve profile design result based on numerical simulation is obviously different from the gas valve profile design result based on recursive optimization;
aiming at the state that the relative stroke of the valve rod is 8mm, valve rod profile design under the correction theoretical target flow is carried out, as shown in fig. 8 and 9, simulation results show that the linearity of the flow matching design result after the correction theoretical target flow is better, and the gas valve profile design result based on numerical simulation is obviously different from the gas valve profile design result based on recursive optimization;
in summary, according to the valve rod molded surface obtained by the actual flow matching design method of the solid attitude and orbit control engine gas valve molded surface based on numerical simulation, the valve rod displacement-valve flow high-linearity change can be ensured, and the design result meets the design requirement.
The foregoing description is only of the preferred embodiments of the present invention and is not intended to limit the scope of the invention, and all equivalent structural changes made by the description of the present invention and the accompanying drawings or direct/indirect application in other related technical fields are included in the scope of the invention.

Claims (8)

1. The actual flow matching design method of the solid attitude and orbit control engine gas valve profile is characterized by comprising the following steps of:
step 1, acquiring a coordinate point set of a spray pipe profile curve, and initializing a flow curve into a theoretical target flow curve;
step 2, determining a corresponding equivalent throat area curve based on the flow curve, and obtaining a valve rod profile curve coordinate point set meeting the equivalent throat area curve based on a recursion matching design method;
step 3, performing computational fluid dynamics simulation based on a spray pipe profile curve coordinate point set and a valve rod profile curve coordinate point set to obtain actual maximum flow of the valve, correcting a theoretical target flow curve according to the actual maximum flow of the valve to obtain an actual target flow curve, and updating the flow curve into the actual target flow curve;
step 4, obtaining an updated valve rod profile curve coordinate point set according to the updated flow curve, and then carrying out computational fluid dynamics simulation according to the spray pipe profile curve coordinate point set and the updated valve rod profile curve coordinate point set to obtain a valve rod displacement-valve flow curve;
step 5, judging whether the current valve rod displacement-valve flow curve meets the convergence condition:
if yes, outputting a current valve rod profile curve coordinate point set;
otherwise, updating the flow curve according to the error between the current valve rod displacement-valve flow curve and the actual target flow curve, and returning to the step 4 to continue iteration;
in step 2, the valve rod profile curve coordinate point set meeting the equivalent throat area curve is obtained based on a recursion matching design method, specifically:
step 2.1, calculating to obtain the tangent point of the spray pipe molded surface and the valve rod molded surface when the valve is completely closed based on a dichotomy
Figure QLYQS_1
；
Step 2.2, solving to obtain the upper point of the spray pipe molded surface
Figure QLYQS_3
The tangential line edge->
Figure QLYQS_6
Negative axis translation +.>
Figure QLYQS_8
Coordinate point corresponding to minimum throat area formed by spray pipe profile>
Figure QLYQS_4
Point->
Figure QLYQS_5
And (4) point->
Figure QLYQS_7
Sequentially adding a valve rod profile curve coordinate point set and initializing recurrence algebra +.>
Figure QLYQS_9
, wherein ,/>
Figure QLYQS_2
Calculating step length;
step 2.3, obtaining the valve rod displacement in the equivalent throat area curve as
Figure QLYQS_10
Equivalent throat area ∈>
Figure QLYQS_14
Searching for satisfying equivalent throat area based on dichotomy>
Figure QLYQS_16
Slope of valve stem profile local taylor expansion line +.>
Figure QLYQS_11
And the minimum throat area formed by the current valve rod profile local Taylor expansion straight line and the spray pipe profile is equal to the equivalent throat area +.>
Figure QLYQS_13
Is a new point +.>
Figure QLYQS_15
And->
Figure QLYQS_17
Adding a valve rod profile curve coordinate point set, wherein,
Figure QLYQS_12
is the relative stroke of the valve rod;
step 2.4, order
Figure QLYQS_18
Judging->
Figure QLYQS_19
Whether or not it is:
if yes, enter step 2.5;
otherwise, returning to the step 2.3 to continue iteration;
step 2.5, adding the circle center as the coordinate point data of the valve rod profile according to the coordinate point set of the valve rod profile curve
Figure QLYQS_20
Radius is->
Figure QLYQS_21
Is a circular arc of (2), wherein>
Figure QLYQS_22
For the slope of the partial taylor expansion line of the valve stem profile at the last iteration of step 2.3 +.>
Figure QLYQS_23
、/>
Figure QLYQS_24
For the point at the last iteration of step 2.3 +.>
Figure QLYQS_25
Is a horizontal and vertical coordinate value of (a);
and 2.6, adding the coordinate point set on the circular arc in the step 2.5 into the valve rod profile curve coordinate point set, and outputting the valve rod profile curve coordinate point set.
2. The method for matching design of actual flow of gas valve profile of solid attitude and orbit control engine according to claim 1, wherein in step 2.1, the tangent point of the spray pipe profile and the valve rod profile when the valve is completely closed is calculated based on the dichotomy
Figure QLYQS_26
The method specifically comprises the following steps:
step 2.1.1, obtaining the valve rod displacement in the equivalent throat area curve as
Figure QLYQS_27
Equivalent throat area ∈>
Figure QLYQS_28
；
Step 2.1.2, searching the valve on the spray pipe profile curve by adopting a dichotomyTangential point of nozzle and valve rod when fully closed
Figure QLYQS_29
So that it satisfies the following equation:
Figure QLYQS_30
，
wherein ,
Figure QLYQS_31
to be at the spot->
Figure QLYQS_32
Straight line edge tangent to the nozzle>
Figure QLYQS_33
Negative axis translation +.>
Figure QLYQS_34
The minimum throat area formed by the rear and the spray pipe profile is that:
Figure QLYQS_35
，
wherein ,
Figure QLYQS_36
is the curve coordinate point set of the spray pipe profile>
Figure QLYQS_37
Respectively a valve rod profile curve coordinate point set, +.>
Figure QLYQS_38
Concentrate the curve coordinate point of the spray pipe profilenThe abscissa value of the dot, +.>
Figure QLYQS_39
Concentrate the curve coordinate point of the valve rod profilepThe abscissa value of the dot.
3. The method for matching and designing the actual flow of the gas valve profile of the solid attitude and orbit control engine according to claim 1, wherein the step 2.3 specifically comprises:
step 2.3.1, initializing a slope search interval
Figure QLYQS_40
And obtain the valve rod displacement in the equivalent throat area curve as
Figure QLYQS_41
Equivalent throat area ∈>
Figure QLYQS_42
, wherein ,/>
Figure QLYQS_43
For maximum value of slope search, +.>
Figure QLYQS_44
Minimum value for slope search;
step 2.3.2, let point
Figure QLYQS_45
For the local expansion point, search slope +.>
Figure QLYQS_46
And let the slope be +.>
Figure QLYQS_47
And go through the local expansion point->
Figure QLYQS_48
As a current local expansion line, wherein the slope +.>
Figure QLYQS_49
The method meets the following conditions: current localityThe minimum throat area formed by the expansion line and the spray pipe profile is equal to the target equivalent throat area under the current valve rod displacement +.>
Figure QLYQS_50
；
Step 2.3.3, the minimum throat area formed by the local Taylor expansion straight line of the current valve rod profile and the spray pipe profile is equal to the equivalent throat area
Figure QLYQS_51
Is a new point +.>
Figure QLYQS_52
And->
Figure QLYQS_53
And adding the valve rod profile curve coordinate point set.
4. The method for matching and designing the actual flow of the gas valve profile of the solid attitude and orbit control engine according to any one of claims 1 to 3, wherein in the step 3, the actual maximum flow of the valve is obtained by performing computational fluid dynamics simulation based on a spray pipe profile curve coordinate point set and a valve rod profile curve coordinate point set, specifically:
importing the spray pipe profile curve coordinate point set and the valve rod profile curve coordinate point set into grid division software to complete unstructured grid division and grid file output;
the grid file is imported into computational fluid dynamics simulation software, and the actual maximum flow of the valve is obtained by solving a two-dimensional axisymmetric conservation type Navie-Stokes equation of compressible flow
Figure QLYQS_54
。/>
5. The method for matching and designing the actual flow of the gas valve profile of the solid attitude and orbit control engine according to claim 4, wherein in the step 3, the theoretical target flow curve is corrected according to the actual maximum flow of the valve and the actual target flow curve is obtained, specifically:
Figure QLYQS_55
，
Figure QLYQS_56
，
in the formula ,
Figure QLYQS_57
for the actual target flow curve, +.>
Figure QLYQS_58
For the theoretical target flow curve, +.>
Figure QLYQS_59
For the actual maximum flow of the valve>
Figure QLYQS_60
Is the theoretical maximum flow of the valve.
6. The method for matching and designing the actual flow of the gas valve profile of the solid attitude and orbit control engine according to any one of claims 1 to 3, wherein in the step 4, the valve stem displacement-valve flow curve is obtained by performing computational fluid dynamics simulation according to a spray pipe profile curve coordinate point set and an updated valve stem profile curve coordinate point set, specifically:
importing the spray pipe profile curve coordinate point set and the updated valve rod profile curve coordinate point set into grid division software to complete unstructured grid division and new grid file output;
the new grid file is imported into computational fluid dynamics simulation software, and a valve rod displacement-valve flow curve is obtained by solving a two-dimensional axisymmetric conservation type Navier-Stokes equation of compressible flow
Figure QLYQS_61
。
7. A method for matching design of actual flow of a gas valve profile of a solid state rail controlled engine according to any one of claims 1 to 3, wherein in step 5, the convergence condition is:
Figure QLYQS_62
，
in the formula ,
Figure QLYQS_63
for the number of valve stem steps, < >>
Figure QLYQS_64
Is the first in the valve flow curveiThe valve flow corresponding to the displacement of the valve rod,
Figure QLYQS_65
is the first in the actual target flow curveiValve flow corresponding to the displacement of the valve rod +.>
Figure QLYQS_66
Is the convergence error.
8. A method for matching design of actual flow of a gas valve profile of a solid state rail controlled engine according to any one of claims 1 to 3, wherein in step 5, the flow curve is updated based on an error between a current valve stem displacement-valve flow curve and an actual target flow curve, specifically:
Figure QLYQS_67
，
in the formula ,
Figure QLYQS_68
for the flow curve +.>
Figure QLYQS_69
For the actual target flow curve, +.>
Figure QLYQS_70
Is a valve stem displacement-valve flow curve. />
CN202211161938.2A 2022-09-23 2022-09-23 Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine Active CN115859492B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211161938.2A CN115859492B (en) 2022-09-23 2022-09-23 Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211161938.2A CN115859492B (en) 2022-09-23 2022-09-23 Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine

Publications (2)

Publication Number Publication Date
CN115859492A CN115859492A (en) 2023-03-28
CN115859492B true CN115859492B (en) 2023-05-05

Family

ID=85661096

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211161938.2A Active CN115859492B (en) 2022-09-23 2022-09-23 Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine

Country Status (1)

Country Link
CN (1) CN115859492B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117952023B (en) * 2024-03-26 2024-06-07 中国人民解放军国防科技大学 Design method for sectional flow regulation valve rod profile of solid attitude and orbit control engine

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110645118B (en) * 2019-10-31 2021-01-05 湖北航天技术研究院总体设计所 Equal-flow solid attitude control thrust device and design method
CN114861352A (en) * 2022-05-05 2022-08-05 中国人民解放军国防科技大学 Method and device for designing valve rod profile of solid attitude and orbit control engine gas valve

Also Published As

Publication number Publication date
CN115859492A (en) 2023-03-28

Similar Documents

Publication Publication Date Title
CN115859492B (en) Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine
CN111191359B (en) Method for constructing space killer zone model of air-defense missile and calculating shooting data
CN103486905B (en) Determining method for terminal guidance shift-exchange conditions of reenter vehicle
CN114233513B (en) Multi-dimensional deflection pea-shaped S-bend spray pipe structure
CN114440711B (en) Four-stage solid carrier rocket trajectory optimization method based on particle swarm optimization
CN106708104B (en) Digital valve flow segmented control algorithm
CN114861352A (en) Method and device for designing valve rod profile of solid attitude and orbit control engine gas valve
CN112199853A (en) Winged missile with steering engine bulge and bulge optimization design method thereof
CN110589010A (en) Hypersonic large-loading-space waverider design method
CN116679743A (en) Nonlinear optimal flight time control guidance method for intercepting maneuvering target
CN115828412A (en) Single-chamber double-thrust solid rocket engine proportioning parameter design method considering outer ballistic constraint
CN113609615B (en) Turbine blade multi-exhaust gas film cold efficiency correction calculation method
CN113189873B (en) Rapid high-precision calculation method for rotating biprism based on symmetric error fitting
CN113359819B (en) Optimal guidance law with collision angle constraint and acceleration limitation
CN112231835B (en) Thrust performance and deflection efficiency integrated vectoring nozzle outlet area optimization method
CN117556537A (en) Turbine blade profile modeling method based on suction surface molded lines and channel width distribution
CN110728010A (en) Method for calculating firing data of antiaircraft gun fire control system
CN110220414B (en) Coincidence method in terminal guided projectile firing plan
CN116596977A (en) Registration method for low-overlapping-rate point cloud
CN109165444B (en) Design method for simultaneously and cooperatively intercepting number of aircrafts and spatial position distribution by multiple aircrafts
Zhang et al. Closed-form time-to-go estimation for proportional navigation guidance considering drag
CN114935277B (en) Online planning method for ideal trajectory of gliding Cheng Zhidao shell
CN111737815B (en) Flight system drop point error correction method based on dimensionless form
CN117952023B (en) Design method for sectional flow regulation valve rod profile of solid attitude and orbit control engine
CN113076610A (en) Design method of binary adjustable air inlet channel

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