CN117952023B - Design method for sectional flow regulation valve rod profile of solid attitude and orbit control engine - Google Patents
Design method for sectional flow regulation valve rod profile of solid attitude and orbit control engine Download PDFInfo
- Publication number
- CN117952023B CN117952023B CN202410349801.2A CN202410349801A CN117952023B CN 117952023 B CN117952023 B CN 117952023B CN 202410349801 A CN202410349801 A CN 202410349801A CN 117952023 B CN117952023 B CN 117952023B
- Authority
- CN
- China
- Prior art keywords
- throat area
- curve
- valve rod
- profile
- area curve
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 67
- 239000007787 solid Substances 0.000 title claims abstract description 24
- 238000009499 grossing Methods 0.000 claims abstract description 48
- 239000007921 spray Substances 0.000 claims abstract description 36
- 230000007704 transition Effects 0.000 claims abstract description 35
- 238000011144 upstream manufacturing Methods 0.000 claims abstract description 31
- 239000002245 particle Substances 0.000 claims description 24
- 238000005457 optimization Methods 0.000 claims description 22
- 230000008569 process Effects 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 10
- 239000006185 dispersion Substances 0.000 claims description 5
- 230000001105 regulatory effect Effects 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000001595 flow curve Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000003380 propellant Substances 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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)
Abstract
The invention discloses a design method for a valve rod profile of a gas valve by sectionally adjusting the flow of a solid attitude and orbit control engine, which comprises the following steps: taking the radius of the upstream transition arc of the throat part of the spray pipe, the convergence half angle and the smoothing times of the equivalent throat area curve as design variables, taking the root mean square error of the minimum target equivalent throat area curve and the design value of the equivalent throat area curve as an objective function, and optimizing to obtain the optimal radius of the upstream transition arc of the throat part of the spray pipe, the convergence half angle and the smoothing times of the equivalent throat area curve; and obtaining and outputting a final design value of the valve rod profile curve and the equivalent throat area curve based on the optimal upstream transition arc radius, the convergence half angle and the smoothing times of the equivalent throat area curve of the nozzle throat. The invention is applied to the field of solid attitude and orbit control engine design, can face any flow regulation requirement, obviously improves the matching degree of the design result and the requirement, and ensures the universality of the design result to any requirement.
Description
Technical Field
The invention relates to the technical field of solid attitude and orbit control engine design, in particular to a method for designing a valve rod profile of a gas valve by sectionally adjusting the flow of a solid attitude and orbit control engine.
Background
The gas type solid attitude and orbit control engine is used as an actuating mechanism of the kinetic energy interceptor, has simple structure, is safe and reliable, and is particularly suitable for systems which are powered by direct force, such as an air base, a sea base, a sky base and the like. For a pintle-type solid attitude and orbit control engine, thrust regulation is largely dependent on the geometry of the gas valve stem, the nozzle profile, and the stem travel (i.e., the axial relative position between the stem and the nozzle). With the diversity and complexity of the demands of the direct force control system, the demands of realizing multi-gear thrust adjustment by a set of gas valve system are becoming urgent, and the existing linear valve design method has difficulty in meeting the multi-gear thrust adjustment demands. Therefore, it is needed to develop the valve rod profile optimization design with the aim of realizing the sectional flow adjustment along with the valve rod travel by using a single valve opening under the condition of a given nozzle throat diameter, and the requirement of multi-gear thrust adjustment of a single set of mechanism is met to the greatest extent.
The valve profile design method commonly used at present comprises the following steps:
Based on the existing profile design (typically spherical or cubic), a design with a relatively close working profile is followed or trimmed. The method is simple, convenient and easy to realize, but needs the participation of experienced engineers, and the prior experience for the gas valve regulated by any section is lacked to be used as reference.
The valve profile optimization design is developed based on an optimization algorithm, the method generally uses the size parameters of a spray pipe and a throat bolt as design variables, uses the error between an obtained valve rod displacement-opening curve and a target curve as an objective function to establish an optimization model, and solves the problem based on an intelligent optimization algorithm. Although the method can realize the automation of the optimization process, additional information such as a valve profile curve equation is needed, the parameterization method of the profile and the selection of a design range are greatly depended, and unreasonable design range and parameterization method easily lead to the phenomenon that a solution meeting the set requirement does not exist. In addition, the method needs to reasonably plan the profile control parameters, and if the parameters are fewer, the profile control is limited, so that the linearity of the optimal design result is difficult to ensure; if the parameters are more, the optimization model is complicated, the calculation cost is unacceptable, and the method is difficult to be directly applied to engineering practice.
A linear valve design method based on recursive matching. According to the method, the valve rod molded surface data point recursion generation can be realized through the thought of piecewise linear approximation according to the maximum flow demand, and the linear flow regulation can be well realized. But there are a number of situations where no feasible solution exists due to the inability to backtrack during the recursion.
Disclosure of Invention
Aiming at the problems that in the prior art, the valve rod displacement of a gas valve of a solid attitude and track control engine and the valve flow are subjected to any sectional adjustment, the design result is poor or no solution exists, the invention provides a valve rod molded surface design method of the gas valve of the solid attitude and track control engine, and on the basis of a recursive matching design method, the invention provides a successive approximation method for dynamically smoothing flow requirements to realize infeasible solution treatment, and solves the design problem of the molded surface of an optimal spray pipe and a throat bolt under any flow adjustment requirement through integrated optimization of the spray pipe and the valve rod, thereby realizing the efficient design of the molded surface of the gas valve under the multistage thrust output of the solid attitude and track control engine, ensuring that the relation between the valve flow and the valve rod travel meets the design requirement to the greatest extent, obviously improving the coincidence degree of the design result and the requirement, and ensuring the universality of the design result to any requirement.
In order to achieve the purpose, the invention provides a design method for the valve rod profile of a gas valve by sectionally adjusting the flow of a solid attitude and orbit control engine, which comprises the following steps:
Taking the radius, convergence half angle and the smoothing times of the equivalent throat area curve of the upstream transition arc of the throat of the spray pipe as design variables, taking the root mean square error of the design values of the minimum equivalent throat area curve and the equivalent throat area curve as an objective function, and optimizing to obtain the optimal radius, convergence half angle and smoothing times of the equivalent throat area curve of the upstream transition arc of the throat of the spray pipe;
And obtaining and outputting a final design value of the valve rod profile curve and the equivalent throat area curve based on the optimal upstream transition arc radius, the convergence half angle and the smoothing times of the equivalent throat area curve of the nozzle throat.
In one embodiment, the calculating process of the equivalent throat area curve design value includes:
smoothing the target equivalent throat area curve based on the current equivalent throat area curve smoothing times;
Calculating a valve rod profile curve which currently meets the smooth target equivalent throat area curve by adopting a recursive matching method;
obtaining a current spray pipe profile curve based on the current spray pipe throat upstream transition arc radius and the convergence half angle;
and obtaining a current equivalent throat area curve design value based on the current valve rod profile curve and the spray pipe profile curve.
In one embodiment, the smoothing the target equivalent throat area curve based on the current equivalent throat area curve smoothing times is specifically:
Repeating the target equivalent throat area curve The secondary substitution is performed on the smoothing model to finish smoothing, wherein/>The current equivalent throat area curve smoothing times;
The smoothing model is as follows:
；
Wherein, Is the/>, on the target equivalent throat area curveDiscrete points/>Is the/>, on the target equivalent throat area curveDiscrete points/>Is the/>, on the target equivalent throat area curveDiscrete points/>Is the number of discrete points of the target equivalent throat area curve.
In one embodiment, the valve rod profile curve currently meeting the smoothed target equivalent throat area curve is calculated by adopting a recursive matching method specifically comprises the following steps:
step 201, according to the 1 st segment of the dispersion curve of the smoothed target equivalent throat area curve ToFor search interval, search/>Tangent point of the minimum hour throat plug and spout profile/>And will tangent point/>The nozzle profile coordinates at the location were used as the first point/>Wherein/>For the minimum throat area of the throat bolt profile and nozzle throat during the search of step 201,/>The abscissa value of the tangent point of the arc of the convergent section and the straight line of the convergent section of the spray pipe;
step 202, dot Stroke increase/>After that become/>As the unfolding point, and obtaining the point on the valve rod profile corresponding to the equivalent throat area according to the valve rod partial profile equation and the spray pipe profile equationPoint/>And dot/>Adding a valve rod molded surface coordinate data set, and initializing recurrence algebra/>Wherein/>Calculating step length;
Step 203, according to the smoothed target equivalent throat area curve Segment discrete curve/>To/>For search interval,/>To expand a point, search/>Optimal slope/>, of the local taylor expansion line of the valve stem profile at minimumAnd taking the valve rod minimum area point corresponding to the valve equivalent throat area calculated by the current valve rod partial profile equation and the spray pipe profile equation as a new point/>Wherein/>Is the minimum value of slope search,/>Maximum value of slope search,/>Searching for a minimum throat area of the throat bolt profile and the nozzle throat during step 203;
Step 204, reducing the coordinate travel of all points in the valve stem profile coordinate data set And the point obtained in step 203/>After the valve rod molded surface coordinate data set is added, judging/>Whether or not it is:
if yes, go to step 205;
Otherwise, let Then go to step 203 to continue iteration;
step 205, adding the circle center as the coordinate point data of the valve rod molded surface according to the existing coordinate point data of the valve rod molded surface in the valve rod molded surface coordinate data set Radius is/>Of (3), wherein,/>、For the point/>, at the last iteration of step 203Is the horizontal and vertical coordinate value of (2)/>For the point/>, at the last iteration of step 203Slope at;
And 206, adding the coordinate point set on the circular arc in the step 205 into the valve rod profile coordinate data set to obtain the valve rod profile curve which currently meets the smoothed target equivalent throat area curve.
In one embodiment, in step 203, the search interval。
In one embodiment, the optimal upstream transition arc radius, convergence half angle and equivalent throat area curve smoothing times of the nozzle throat are obtained by optimizing a continuous discrete particle swarm algorithm.
In one embodiment, in the process of optimizing by adopting a continuous discrete particle swarm optimization algorithm to obtain the optimal transition arc radius, convergence half angle and equivalent throat area curve smoothing times at the upstream of the throat of the spray pipe, the updating process of the particle speed and the position is as follows:
；
；
Wherein, For/>First/>, in the next iterationThe individual particles are at the/>Speed of dimension,/>、/>For/>First/>, in the next iterationThe individual particles are at the/>Speed and position of dimension,/>For/>First/>, in the next iterationThe individual particles are at the/>Dimension found historical optimal solution,/>For/>Historical optimal first/>, found by population in multiple iterationsDimensional coordinates,/>Representing inertial weights,/>、/>As learning factor,/>、/>For interval/>Random number between,/>For/>The first generation groupLocation of individual-For/>First/>, in the generation populationLocation of individual-For/>The first generation groupThe speed of the individual.
In one embodiment, in the process of optimizing by adopting a continuous discrete particle swarm optimization algorithm to obtain the optimal upstream transition arc radius, convergence half angle and equivalent throat area curve smoothing times of the nozzle throat, the convergence judging process is as follows:
If the optimal solution is not updated within 10 times, judging convergence, and outputting the upstream transition arc radius, the convergence half angle and the equivalent throat area curve smoothing times of the nozzle throat corresponding to the current optimal solution.
Compared with the prior art, the invention has the following beneficial technical effects:
1. Compared with a recursion matching design method, the invention realizes the integrated optimization of the spray pipe molded surface and the valve rod molded surface, and solves the problem that the recursion matching design is free from solution;
2. Compared with the conventional design scheme of the on-edge or fine-tuning, the multi-gear thrust and flow adjustment is realized, the design result is more stable and reliable, and the universality of the design result for any requirement is ensured;
3. the invention provides a method for smooth processing of sectional requirements, which realizes smooth transition between sectional flow regulation and reduces the design difficulty of a valve rod.
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 method for designing a valve rod profile of a gas valve by sectionally adjusting the flow of a solid attitude and orbit control engine in an embodiment of the invention;
FIG. 2 is a schematic view of linear gas valve profile calculation based on recursive optimization in an embodiment of the present invention, wherein: (a) A schematic diagram before recursion, and (b) a schematic diagram after recursion;
FIG. 3 is a flowchart of a continuous discrete particle swarm algorithm according to an embodiment of the invention;
FIG. 4 is a schematic diagram of the valve rod profile of a gas regulating valve of a solid attitude and orbit control engine in an embodiment of the invention;
FIG. 5 is a schematic diagram of a target flow trip curve for case 1 according to an embodiment of the present invention;
FIG. 6 is a graphical representation of the target equivalent throat area for case 1 in an embodiment of the present invention;
Fig. 7 is a schematic diagram of an objective function iteration curve of case 1 in an embodiment of the present invention;
FIG. 8 is a comparison of valve stem design profile results for case 1 in accordance with an embodiment of the present invention;
FIG. 9 is a schematic diagram of a target flow trip for case 2 according to an embodiment of the present invention;
FIG. 10 is a graphical representation of the target equivalent throat area for case 2 in an embodiment of the present invention;
Fig. 11 is a comparison of valve stem design profile results for case 2 in accordance with an embodiment of the present invention.
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 a solid attitude and track control engine flow sectionalized adjustment gas valve rod molded surface design method, which is based on a recursive matching design method, and is used for realizing smooth transition of arbitrary sectionalized flow adjustment by adopting data smoothing and realizing transition sharp point smooth range adjustment by selecting smooth times, realizing smooth flow matching design by adopting the recursive matching design method, and further ensuring that a travel flow curve obtained by design meets preset design requirements as far as possible by optimizing smooth times and a spray pipe convergence section molded surface.
Referring to fig. 1, the design method for the valve rod profile of the gas valve by adjusting the flow of the solid attitude and orbit control engine in a sectional manner specifically comprises the following steps:
Taking the upstream transition arc radius, convergence half angle and equivalent throat area curve smoothing times of the nozzle throat as design variables, and taking the root mean square error of the minimum target equivalent throat area curve and the design value of the equivalent throat area curve as an objective function;
Optimizing by adopting a continuous discrete particle swarm optimization to obtain the optimal upstream transition arc radius, convergence half angle and equivalent throat area curve smoothing times of the nozzle throat;
And obtaining and outputting a final design value of the valve rod profile curve and the equivalent throat area curve based on the optimal upstream transition arc radius, the convergence half angle and the smoothing times of the equivalent throat area curve of the nozzle throat.
In this embodiment, the target equivalent throat area curve may be based on a preset target flow stroke curveThe method is calculated and obtained in the specific implementation process:
Discretizing a target flow stroke curve within the total stroke as Segment, correspond to/>Individual point flow course data points，/>General/>100-200, Corresponding maximum travel is/>The stroke calculation step length is；
Calculating discrete points of a target equivalent throat area curve based on the discrete flow stroke curve, wherein the discrete points are as follows:
；
Wherein, Is the/>, on the target equivalent throat area curveDiscrete points, i.e. minimum passage area between nozzle and throat surfaces when the throat is moved to a specific position,/>On the target flow travel curve/>Flow of discrete points,/>Is the throat diameter of the spray pipe,/>For gasifier pressure,/>Is a characteristic speed of the propellant.
The method takes the upstream transition arc radius, convergence half angle and equivalent throat area curve smoothing times of the nozzle throat as design variables, and takes the root mean square error of the minimum target equivalent throat area curve and the design value of the equivalent throat area curve as an objective function, which comprises the following specific steps:
；
Wherein, Design value for equivalent throat area curveDiscrete points/>For transition arc radius upstream of nozzle throat/>For convergence half angle,/>For the number of times of smoothing the equivalent throat area curve,/>Root mean square error for the target equivalent throat area curve and the design value of the equivalent throat area curve, which is AND/>A function of the correlation.
In the present embodiment, for any givenThe calculation process of the equivalent throat area curve design value comprises the following steps:
Step 1, smoothing the target equivalent throat area curve based on the current number of times of smoothing the equivalent throat area curve, specifically, repeating each discrete point in the target equivalent throat area curve And substituting the target equivalent throat area curve with the sharp point into a smooth transition target equivalent throat area curve, wherein the smooth model is as follows:
；
Wherein, Is the/>, on the target equivalent throat area curveDiscrete points/>Is the/>, on the target equivalent throat area curveA plurality of discrete points;
Step 2, calculating to obtain a valve rod profile curve which currently meets the smooth target equivalent throat area curve by adopting a recursive matching method, wherein the specific implementation process is as follows:
step 201, according to the 1 st segment of the dispersion curve of the smoothed target equivalent throat area curve ToFor search interval, search/>Tangent point of the minimum hour throat plug and spout profile/>And will tangent point/>The nozzle profile coordinates at the location were used as the first point/>Wherein/>Is the abscissa value of the tangent point of the arc of the convergent section and the straight line of the convergent section of the spray pipe. Specifically, to/>For design variables,/>The first-order Taylor expansion of the spray pipe profile at the position is the local profile of the laryngeal plug, and the current profile stroke increase/> is obtained by adopting a dichotomy searchWhen the throat bolt profile and the minimum throat area/>, of the nozzle throatTangent point nearest to the target equivalent throat area/>Namely, the following problems are solved by adopting a dichotomy method:
；
；
step 202, dot Stroke increase/>After that become/>As the unfolding point, and obtaining the point on the valve rod profile corresponding to the equivalent throat area according to the valve rod partial profile equation and the spray pipe profile equationPoint/>As a second point of the laryngeal plug, shown in figure 2; and will be pointedAnd dot/>Adding valve rod profile coordinate data set and initializing recurrence algebra；
Step 203, according to the smoothed target equivalent throat area curveSegment dispersion curveCalculating the/>, based on a binary search methodLocal taylor expansion slope and/>, of segment dispersion curvePoints, i.e./>For search interval,/>To expand a point, searchOptimal slope/>, of the local taylor expansion line of the valve stem profile at minimumAnd taking the valve rod minimum area point corresponding to the valve equivalent throat area calculated by the current valve rod partial profile equation and the spray pipe profile equation as a new point/>Wherein/>Is the minimum value of slope search,/>Search Interval/>, maximum for slope search. As in step 201, to/>The minimum throat area/>, of the throat bolt molded surface and the nozzle throat, is obtained by adopting a dichotomy search for design variablesClosest to the tangent point of the target equivalent throat area, the following problem is solved by adopting a dichotomy:
；
；
Step 204, reducing the coordinate travel of all points in the valve stem profile coordinate data set And the point obtained in step 203/>After the valve rod molded surface coordinate data set is added, judging/>Whether or not it is:
if yes, go to step 205;
Otherwise, let Then go to step 203 to continue iteration;
step 205, adding the circle center as the coordinate point data of the valve rod molded surface according to the existing coordinate point data of the valve rod molded surface in the valve rod molded surface coordinate data set Radius is/>Of (3), wherein,/>For the point/>, at the last iteration of step 203Slope at/(I)、/>For the point/>, at the last iteration of step 203Is a horizontal and vertical coordinate value of (a);
step 206, adding the coordinate point set on the circular arc in the step 205 into the valve rod profile coordinate data set to obtain a valve rod profile curve which currently meets the smoothed target equivalent throat area curve;
step 3, obtaining a current spray pipe profile curve based on the current spray pipe throat upstream transition arc radius and the convergence half angle, wherein the calculation process is as follows:
；
Wherein, Is a point on the profile curve of the nozzle;
And 4, obtaining a current equivalent throat area curve design value based on the current valve rod profile curve and the spray pipe profile curve, namely calculating the deviation between a design result and a design requirement according to the current equivalent throat area curve design value and a target equivalent throat area curve, and obtaining the objective function of the optimization problem in the step 203.
In this embodiment, the continuous discrete particle swarm optimization is adopted to obtain the optimal transition arc radius, convergence half angle and equivalent throat area curve smoothing times at the upstream of the throat of the spray pipe, and referring to fig. 3, the specific implementation process is as follows:
First, initializing the upstream transition arc radius of the nozzle throat Convergence half angle/>Equivalent laryngeal area Curve smoothing times/>Designed interval and search population of/>、/>、Will/>Relaxing into continuous domains/>And generates 20 initial sample points.
Secondly, calculating the fitness of each sample point, namely the deviation between the design result of each sample point and the design requirement, namely the objective function; in the calculation, in each sample pointRounding to the last integer, i.e./>Wherein/>Is an upward rounding operation;
The velocity and position of each particle in the population is then updated using the following equation, i.e., assuming at the first In the next iteration, for the/>Individual particles/>Its/>Dimensional speed and position use/>And/>Representation of/>Historical optimal solution for particle finding/>First/>The dimensional coordinates are/>，/>Historical optimal first/>, found for populationDimensional coordinates, then at/>In the next iteration,/>The updating process of the particle speed and the position comprises the following steps:
；
；
Wherein, Representing inertial weights,/>As a learning factor, we generally get/>，/>For interval/>Random number between,/>For/>First/>, in the generation populationLocation of individual-For/>The first generation groupLocation of individual-For/>First/>, in the generation populationThe speed of the individual;
then, calculating the fitness of the new population, and updating the particle history optimum and the population history optimum;
Finally, convergence judgment is carried out: if the optimal solution is not updated in X times (for example, 10 times), judging convergence, and outputting the upstream transition circular arc radius of the nozzle throat corresponding to the current optimal solution Convergence half angle/>Equivalent laryngeal area Curve smoothing times/>Otherwise, the next iteration is carried out, and the population calculation is updated again.
At the upstream transition arc radius of the nozzle throat to obtain the optimumConvergence half angle/>Equivalent laryngeal area Curve smoothing times/>After that, can be based on/>The recursion matching method in the step 2 is carried out to obtain the final valve rod profile curve and is based on/>、/>And 3, obtaining a final spray pipe profile curve by performing the step, and finally obtaining a final equivalent throat area curve design value and a flow stroke curve design value.
The invention is further described below with reference to specific examples.
Taking two solid attitude and orbit control power gas valves as an example, a specific implementation scheme of the valve rod profile design is given, and the geometric configuration of the gas valve is shown in fig. 4.
Case 1: two sections are regulated in a segmented way, the maximum thrust flow is 1kg/s, and the maximum stroke is 20mm; the first section of travel is 10mm, and the flow rate is changed from 0 to 0.2kg/s; the second section has a stroke of 10mm, and the flow rate is changed from 0.2kg/s to 1kg/s; the corresponding specific target flow path curve is shown in fig. 5. The valve stem profile design embodiment is as follows:
1) Obtaining a target equivalent throat area curve according to the target flow stroke curve, namely an equivalent throat area and valve rod stroke demand curve shown in fig. 6;
2) Calculating the throat diameter to be 11.83mm according to the maximum equivalent throat area;
3) Optimizing the above-mentioned profile optimization problem by adopting a continuous-discrete particle swarm optimization algorithm by taking the radius of the transition circular arc at the upstream of the throat of the spray pipe, the convergence half angle and the smoothing times of the equivalent throat area curve as design variables, wherein the corresponding objective function iteration curve is shown in figure 7;
4) The obtained optimal solution is 、/>、/>；
5) Inputting the target equivalent throat area curve after smoothing for 13 times into an equivalent throat area recursion matching design algorithm, and calculating to obtain a corresponding valve rod profile bus, namely shown in fig. 8;
6) And calculating a change curve of the equivalent throat area along with the valve rod stroke based on the valve rod profile bus, and further calculating a flow stroke curve according to the relation between the flow and the equivalent throat area.
Simulation calculation results show that the method based on the combination of the optimal design and the recursive matching is more in accordance with the design requirement than the result obtained by direct recursive matching, and root mean square error calculation results of the design requirement and the design result are shown in a table 1, so that the method can obtain a reasonable design result which is more in accordance with the design requirement through continuous discrete mixing optimization on the premise of poor direct design effect.
Table 1 error comparison of different design methods
Case 2: three sections are regulated in a segmented way, the maximum thrust flow is 1kg/s, and the maximum stroke is 24mm; the first section of travel is 8mm, and the flow rate is changed from 0 to 0.1kg/s; the second section has a travel of 8mm and the flow rate is changed from 0.1kg/s to 0.8kg/s; the third section has a travel of 8mm and the flow rate is changed from 0.8kg/s to 1kg/s; the corresponding specific target flow trip curve is shown in fig. 9. The valve stem profile design embodiment is as follows:
1) Obtaining a target equivalent throat area curve according to the target flow stroke curve, namely an equivalent throat area and valve rod stroke demand curve shown in fig. 10;
2) Calculating the throat diameter to be 11.83mm according to the maximum equivalent throat area;
3) The optimization problem of the molded surface is optimized by adopting a continuous-discrete particle swarm optimization algorithm by taking the radius of an upstream transition arc of the throat part of the spray pipe, the convergence half angle and the smoothing frequency of an equivalent throat area curve as design variables, and the obtained optimal solution is 、/>、/>；
4) Inputting the target equivalent throat area curve after smoothing for 10 times into an equivalent throat area recursion matching design algorithm, and calculating to obtain a corresponding valve rod profile bus, namely shown in FIG. 11;
5) And calculating a change curve of the equivalent throat area along with the valve rod stroke based on the valve rod profile bus, and further calculating a flow stroke curve according to the relation between the flow and the equivalent throat area.
Simulation calculation results show that the method based on the combination of the optimal design and the recursive matching meets the design requirement better than the result obtained by direct recursive matching, and root mean square error calculation results of the design requirement and the design result are shown in a table 2, so that the method can obtain a reasonable design result which meets the design requirement better through continuous discrete mixing optimization on the premise of poor direct design effect.
Table 2 error comparison of different design methods
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 (7)
1. A design method for a valve rod profile of a gas valve by sectionally adjusting the flow of a solid attitude and orbit control engine is characterized by comprising the following steps:
Taking the radius, convergence half angle and the smoothing times of the equivalent throat area curve of the upstream transition arc of the throat of the spray pipe as design variables, taking the root mean square error of the design values of the minimum equivalent throat area curve and the equivalent throat area curve as an objective function, and optimizing to obtain the optimal radius, convergence half angle and smoothing times of the equivalent throat area curve of the upstream transition arc of the throat of the spray pipe;
based on the optimal upstream transition arc radius, convergence half angle and equivalent throat area curve smoothing times of the nozzle throat, obtaining and outputting a final valve rod profile curve and equivalent throat area curve design value;
The calculation process of the equivalent throat area curve design value comprises the following steps:
smoothing the target equivalent throat area curve based on the current equivalent throat area curve smoothing times;
Calculating a valve rod profile curve which currently meets the smooth target equivalent throat area curve by adopting a recursive matching method;
obtaining a current spray pipe profile curve based on the current spray pipe throat upstream transition arc radius and the convergence half angle;
and obtaining a current equivalent throat area curve design value based on the current valve rod profile curve and the spray pipe profile curve.
2. The method for designing the valve rod profile of the gas valve by sectionally adjusting the flow of the solid attitude and orbit control engine according to claim 1, wherein the smoothing of the target equivalent throat area curve based on the current smoothing times of the equivalent throat area curve is specifically as follows:
Repeating the target equivalent throat area curve The smoothing model is substituted for times to finish smoothing, wherein n is the current smoothing times of the equivalent throat area curve;
The smoothing model is as follows:
Wherein, Is the/>, on the target equivalent throat area curveDiscrete points/>Is the/>, on the target equivalent throat area curveDiscrete points/>Is the/>, on the target equivalent throat area curveDiscrete points/>+1 Is the number of discrete points of the target equivalent throat area curve.
3. The method for designing the valve rod profile of the solid attitude and orbit control engine flow sectionally-regulated gas valve according to claim 1, wherein the valve rod profile curve currently meeting the smoothed target equivalent throat area curve is calculated by adopting a recursive matching method, and is specifically as follows:
step 201, according to the 1 st segment of the dispersion curve of the smoothed target equivalent throat area curve To/>For search interval, search/>Tangent point of the minimum hour throat plug and spout profile/>And will tangent point/>The nozzle profile coordinates at the location were used as the first point/>Wherein/>For the minimum throat area of the throat bolt profile and nozzle throat during the search of step 201,/>The abscissa value of the tangent point of the arc of the convergent section and the straight line of the convergent section of the spray pipe;
step 202, dot Stroke increase/>After that become/>As the unfolding point, and according to the valve rod partial profile equation and the spray pipe profile equation, obtaining the point/> on the valve rod profile corresponding to the equivalent throat areaPoint of the dotAnd dot/>Adding a valve rod molded surface coordinate data set, and initializing recurrence algebra/>Wherein, the method comprises the steps of, wherein,Calculating step length;
Step 203, according to the smoothed target equivalent throat area curve Segment discrete curve/>ToFor search interval,/>To expand a point, search/>Optimal slope/>, of the local taylor expansion line of the valve stem profile at minimumAnd taking the valve rod minimum area point corresponding to the valve equivalent throat area calculated by the current valve rod partial profile equation and the spray pipe profile equation as a new point/>Wherein/>Is the minimum value of slope search,/>Maximum value of slope search,/>Searching for a minimum throat area of the throat bolt profile and the nozzle throat during step 203;
Step 204, reducing the coordinate travel of all points in the valve stem profile coordinate data set And the point obtained in step 203/>After the valve rod molded surface coordinate data set is added, judging/>Whether or not it is:
if yes, go to step 205;
Otherwise, let Then go to step 203 to continue iteration;
step 205, adding the circle center as the coordinate point data of the valve rod molded surface according to the existing coordinate point data of the valve rod molded surface in the valve rod molded surface coordinate data set Radius is/>Of (3), wherein,/>、/>For the point/>, at the last iteration of step 203Is the horizontal and vertical coordinate value of (2)/>For the point/>, at the last iteration of step 203Slope at;
And 206, adding the coordinate point set on the circular arc in the step 205 into the valve rod profile coordinate data set to obtain the valve rod profile curve which currently meets the smoothed target equivalent throat area curve.
4. The method for designing the valve stem profile of a gas valve for the staged flow adjustment of a solid state rail controlled engine as defined in claim 3, wherein in step 203, the search interval is set。
5. The method for designing the valve rod profile of the gas valve by sectionally adjusting the flow of the solid attitude and orbit control engine according to any one of claims 1 to 4, wherein the optimal upstream transition arc radius, convergence half angle and equivalent throat area curve smoothing times of the nozzle throat are obtained by optimizing a continuous discrete particle swarm algorithm.
6. The method for designing the valve rod profile of the gas valve by sectionally adjusting the flow of the solid attitude and orbit control engine according to claim 5, wherein in the process of optimizing by adopting a continuous discrete particle swarm algorithm to obtain the optimal upstream transition arc radius, convergence half angle and equivalent throat area curve smooth times of the nozzle throat, the updating process of the particle speed and the position is as follows:
Wherein, For/>First/>, in the next iterationThe individual particles are at the/>Speed of dimension,/>、/>For/>First/>, in the next iterationThe individual particles are at the/>Speed and position of dimension,/>For/>First/>, in the next iterationThe individual particles are at the/>The historical optimal solution found is maintained,For/>Historical optimal first/>, found by population in multiple iterationsDimensional coordinates,/>Representing inertial weights,/>、/>As learning factor,/>、/>For interval/>Random number between,/>For/>First/>, in the generation populationLocation of individual-For/>First/>, in the generation populationLocation of individual-For/>First/>, in the generation populationThe speed of the individual.
7. The method for designing the valve rod profile of the gas valve by sectionally adjusting the flow of the solid attitude and orbit control engine according to claim 6, wherein in the process of optimizing by adopting a continuous discrete particle swarm optimization to obtain the optimal upstream transition arc radius, convergence half angle and equivalent throat area curve smooth times of the nozzle throat, the convergence judgment process is as follows:
If the optimal solution is not updated within 10 times, judging convergence, and outputting the upstream transition arc radius, the convergence half angle and the equivalent throat area curve smoothing times of the nozzle throat corresponding to the current optimal solution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410349801.2A CN117952023B (en) | 2024-03-26 | 2024-03-26 | Design method for sectional flow regulation valve rod profile of solid attitude and orbit control engine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410349801.2A CN117952023B (en) | 2024-03-26 | 2024-03-26 | Design method for sectional flow regulation valve rod profile of solid attitude and orbit control engine |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117952023A CN117952023A (en) | 2024-04-30 |
| CN117952023B true CN117952023B (en) | 2024-06-07 |
Family
ID=90794698
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410349801.2A Active CN117952023B (en) | 2024-03-26 | 2024-03-26 | Design method for sectional flow regulation valve rod profile of solid attitude and orbit control engine |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117952023B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102362044A (en) * | 2009-02-17 | 2012-02-22 | 埃克森美孚研究工程公司 | Coated oil and gas well production devices |
| US10706361B1 (en) * | 2015-12-11 | 2020-07-07 | The Boeing Company | Hybrid feature selection for performance prediction of fluid control valves |
| CN113496062A (en) * | 2021-07-26 | 2021-10-12 | 中国人民解放军国防科技大学 | Method for designing motion time sequence of throat bolt of attitude and orbit control engine under stable pressure constraint |
| 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 |
| CN115859492A (en) * | 2022-09-23 | 2023-03-28 | 中国人民解放军国防科技大学 | Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine |
| CN117236025A (en) * | 2023-09-20 | 2023-12-15 | 河北工业大学 | Multi-valve adjustment thrust steady-state uncertainty analysis method for solid attitude and orbit control engine |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090038315A1 (en) * | 2007-08-06 | 2009-02-12 | Robert Paul Johnson | Ion pump and an electrochemical engine using same |
-
2024
- 2024-03-26 CN CN202410349801.2A patent/CN117952023B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102362044A (en) * | 2009-02-17 | 2012-02-22 | 埃克森美孚研究工程公司 | Coated oil and gas well production devices |
| US10706361B1 (en) * | 2015-12-11 | 2020-07-07 | The Boeing Company | Hybrid feature selection for performance prediction of fluid control valves |
| CN113496062A (en) * | 2021-07-26 | 2021-10-12 | 中国人民解放军国防科技大学 | Method for designing motion time sequence of throat bolt of attitude and orbit control engine under stable pressure constraint |
| 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 |
| CN115859492A (en) * | 2022-09-23 | 2023-03-28 | 中国人民解放军国防科技大学 | Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine |
| CN117236025A (en) * | 2023-09-20 | 2023-12-15 | 河北工业大学 | Multi-valve adjustment thrust steady-state uncertainty analysis method for solid attitude and orbit control engine |
Non-Patent Citations (2)
| Title |
|---|
| A surrogate-based optimization method for mixed-variable aircraft design;Zeping Wu等;《Engineering Optimization》;20210131;全文 * |
| 固体姿轨控发动机变工况过程推力不确定性分析;高经纬等;《国防科技大学学报》;20240125;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117952023A (en) | 2024-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102996511B (en) | CAD (Computer-Aided Design) aided design method of controllable diffused blade profile with curvature continuity | |
| CN102880052B (en) | Time scale function decomposition based hypersonic aircraft actuator saturation control method | |
| CN104776446A (en) | Combustion optimization control method for boiler | |
| CN113051775B (en) | Horizontal directional drilling track optimization method based on improved radial movement algorithm | |
| CN112115544A (en) | Rocket track optimization method based on improved genetic algorithm | |
| CN112903233B (en) | Optimization method for accurate positioning of two-stage series servo oil cylinder driving attack angle mechanism | |
| CN105977649A (en) | Shaped surface-oriented quick determination method for adjustment quantity of active panel of large parabolic antenna | |
| CN112115558B (en) | Steam turbine cascade molded line parameterized reconstruction method based on crowd searching algorithm | |
| CN112347579A (en) | Compressor blade profile design method and compressor blade profile | |
| CN106089801A (en) | A kind of compressor blade formative method | |
| CN104252132A (en) | Adaptive genetic algorithm-based interplanetary orbit control optimization method | |
| CN117952023B (en) | Design method for sectional flow regulation valve rod profile of solid attitude and orbit control engine | |
| CN108999845B (en) | Three-dimensional variable cross-section curved flow channel design method and device based on geometric fusion | |
| CN110589010A (en) | Hypersonic large-loading-space waverider design method | |
| CN113467241A (en) | Method for optimizing burn-up of convex curvature landing trajectory | |
| CN112231835B (en) | Thrust performance and deflection efficiency integrated vectoring nozzle outlet area optimization method | |
| CN115859492B (en) | Actual flow matching design method for molded surface of gas valve of solid attitude and orbit control engine | |
| CN112231836B (en) | Wing profile optimization method based on genetic algorithm and numerical simulation | |
| CN109606728A (en) | Method and system for designing precursor of hypersonic aircraft | |
| CN113741509A (en) | Hypersonic glide aircraft downward pressing section energy management method | |
| CN110245408B (en) | Design method for single-arc pressure surface blade profile of steam turbine | |
| CN116700013A (en) | Mixed optimization method for fastest climbing track of aircraft | |
| CN115358101A (en) | Jet pipe design method based on sound velocity solution and characteristic line reverse thrust | |
| CN109800519A (en) | A kind of aerofoil fan multiple spot Aerodynamic optimization design method | |
| CN117329147A (en) | Impeller integrated structure and design method thereof |
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 |