Balanced rotating speed estimation method and device of projectile and storage medium

Abstract

The embodiment of the invention provides a projectile balance rotating speed estimation method and device and a storage medium, and the method comprises the steps: constructing a projectile geometric model according to preset projectile data and a preset empennage inclination angle; based on the projectile geometric model, performing steady calculation according to preset working condition data to obtain surface flow field data and space flow field data; according to the surface flow field data, the average value of pressure center or centroid coordinates on all empennages of the projectile and the average value of the radius from the pressure center or the centroid to a projectile shaft are obtained; according to the space flow field data, the average value of the pressure center or centroid coordinates and the average value of the radius from the pressure center or the centroid to the elastic shaft, the speed average value is obtained; and the balance rotating speed of the projectile is obtained according to the average value of the radius from the pressure center or the centroid to the projectile shaft, the speed average value and a preset empennage inclined angle. The calculation efficiency can be improved while the data acquisition precision is improved, so that the balance rotating speed calculation result precision is higher.

Description

Technical Field

[0001] This document relates to the field of aerodynamics, and in particular to a method, device and storage medium for estimating the equilibrium rotation speed of a projectile. Background Art

[0002] The balanced rotation speed of the projectile with inclined tail fins determines the stability and shooting accuracy of the projectile.

[0003] The existing technology mainly estimates the equilibrium rotation speed of the oblique tail projectile through an engineering estimation method based on the tail geometric parameters and the incoming flow parameters or an interpolation calculation method of the rolling moment coefficient at different rolling angular velocities obtained by CFD numerical calculation.

[0004] However, the engineering estimation method in the existing technology needs to use prior knowledge, such as nonlinear correction coefficients related to the shape, flow-related compressibility correction coefficients, friction resistance coefficients of the flat plate, etc. These prior data cannot cover all states, resulting in large errors in the estimated equilibrium speed. The interpolation calculation method uses unsteady flow field calculations, which require a long time to carry out unsteady flow field calculations under a constant roll angular velocity and have low calculation efficiency. Summary of the invention

[0005] In view of the above scheme, the present application aims to propose a method, device and storage medium for estimating the equilibrium rotation speed of a projectile to solve at least one of the above technical problems.

[0006] In a first aspect, one or more embodiments of this specification provide a method for estimating a balanced rotation speed of a projectile with an oblique tail fin, comprising:

[0007] Constructing a projectile geometry model according to preset projectile data and preset tail wing tilt angle;

[0008] Based on the projectile geometric model, a steady calculation is performed according to preset working condition data to obtain surface flow field data and spatial flow field data;

[0009] According to the surface flow field data, an average value of the coordinates of the pressure center or the centroid on each tail wing of the projectile and an average value of the radius of the pressure center or the centroid from the projectile axis are obtained;

[0010] Obtaining an average velocity value according to the spatial flow field data, an average value of the pressure center or the centroid coordinates, and an average value of the radius of the pressure center or the centroid from the projectile axis; and

[0011] The equilibrium rotation speed of the projectile is obtained according to the average value of the radius of the pressure center or the centroid from the projectile axis, the average speed and the preset tail fin tilt angle.

[0012] Further, meshing the projectile geometric model to obtain a computational fluid dynamics mesh;

[0013] Set the height and speed of the projectile according to the preset working condition data; and

[0014] According to the height and speed of the projectile, the surface flow field data and the spatial flow field data of the projectile are determined by using a computational fluid dynamics grid.

[0015] Further, according to the surface flow field data, the pressure center or centroid coordinates of each tail fin on the projectile are extracted; and

[0016] According to the pressure center or centroid coordinates of each tail wing, an average value of the pressure center or centroid coordinates and an average value of the radius of the pressure center or centroid from the missile axis are obtained.

[0017] Further, according to the spatial flow field data and the average value of the radius from the pressure center or the centroid to the projectile axis, an isosurface value is obtained;

[0018] Obtaining an isoline value according to the isosurface value and the average value of the pressure center or centroid coordinates; and

[0019] According to the isosurface value and the isoline value, an average value of the velocity on the isoline is obtained.

[0020] Furthermore, the calculation formula of the equilibrium rotation speed of the projectile is as follows:

[0021] ω x =u * *sin(δ) / r*

[0022] Among them, ω x Indicates the equilibrium speed;

[0023] x represents the coordinate axis;

[0024] u * Indicates the average speed;

[0025] δ represents the tail wing cant angle; and

[0026] r* represents the average value of the radius from the compression center or centroid to the projectile axis.

[0027] In a second aspect, an embodiment of the present application provides a device for estimating a balanced rotation speed of a projectile with an inclined tail fin, comprising:

[0028] A construction module, used for constructing a projectile geometric model according to preset projectile data and preset tail fin tilt angle;

[0029] A data acquisition module, used to perform steady calculation based on the projectile geometric model and preset working condition data to obtain surface flow field data and spatial flow field data;

[0030] A first calculation module is used to obtain the average value of the coordinates of the pressure center or the centroid on each tail wing of the projectile and the average value of the radius of the pressure center or the centroid from the projectile axis according to the surface flow field data;

[0031] A second calculation module is used to obtain a velocity average value according to the spatial flow field data, an average value of the pressure center or the centroid coordinates, and an average value of the radius of the pressure center or the centroid from the projectile axis; and

[0032] The third calculation module is used to obtain the equilibrium rotation speed of the projectile according to the average value of the radius of the pressure center or the centroid from the projectile axis, the average speed and the preset tail fin tilt angle.

[0033] Furthermore, the data acquisition module is configured as follows:

[0034] Meshing the projectile geometric model to obtain a computational fluid dynamics grid;

[0035] Set the height and speed of the projectile according to the preset working condition data; and

[0036] According to the height and speed of the projectile, the surface flow field data and the spatial flow field data of the projectile are determined by using a computational fluid dynamics grid.

[0037] Furthermore, the first computing module is configured to:

[0038] Extracting the pressure center or centroid coordinates of each tail fin on the projectile according to the surface flow field data; and

[0039] According to the pressure center or centroid coordinates of each tail wing, an average value of the pressure center or centroid coordinates and an average value of the radius of the pressure center or centroid from the missile axis are obtained.

[0040] Furthermore, the second computing module is configured to:

[0041] Obtaining an isosurface value according to the spatial flow field data and an average value of the radius of the pressure center or the centroid from the projectile axis;

[0042] Obtaining an isoline value according to the isosurface value and the average value of the pressure center or centroid coordinates; and

[0043] According to the isosurface value and the isoline value, an average value of the velocity on the isoline is obtained.

[0044] In a third aspect, an embodiment of the present application provides a storage medium for storing computer-executable instructions, characterized in that when the computer-executable instructions are executed, the steps of the method for estimating the equilibrium rotation speed of a projectile with an inclined tail fin described in any one of the first aspects are implemented.

[0045] Compared with the prior art, this application can at least achieve the following technical effects:

[0046] The present application can estimate the equilibrium rotation speed of the oblique tail fin projectile geometric model under any working condition, thereby improving the accuracy of data source acquisition; and only one steady-state calculation is required for the oblique tail fin projectile geometric model, thereby improving the calculation efficiency and achieving higher accuracy of the results. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate one or more embodiments of this specification or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments recorded in this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0048] Figure 1 A flow chart of a method for estimating the equilibrium rotation speed of a projectile with an oblique tail fin provided for one or more embodiments of this specification;

[0049] Figure 2 A geometric model diagram of a projectile with an oblique tail provided for one or more embodiments of this specification;

[0050] Figure 3 A schematic diagram of the structure of a balanced rotation speed estimation device for a projectile with an inclined tail fin provided in one or more embodiments of the present specification. DETAILED DESCRIPTION

[0051] In order to enable those skilled in the art to better understand the technical solutions in one or more embodiments of this specification, the following will be combined with the drawings in one or more embodiments of this specification to clearly and completely describe the technical solutions in one or more embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this specification, not all of the embodiments. Based on one or more embodiments of this specification, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of this document.

[0052] At present, there are two main methods for estimating the equilibrium speed of rotating missiles. One is the engineering estimation method based on the geometric parameters of the tail and the incoming flow parameters, and the other is the interpolation calculation method of the rolling moment coefficient at different roll angular velocities obtained based on CFD numerical calculations. The engineering estimation method is fast and only requires the shape parameters and the corresponding working condition parameters to carry out the calculation, but the error is large, about 10%, and it also requires prior knowledge, such as the nonlinear correction coefficient related to the shape, the compressibility correction coefficient related to the flow, the friction resistance coefficient of the plate, etc. These coefficients obtained from public data cannot cover all states. Therefore, the engineering estimation method is not available under some working conditions due to large errors; the interpolation method has a high calculation accuracy, and the error accuracy is less than 10%, but the interpolation method requires more unsteady flow field calculations at constant roll angular velocity. The interpolation method estimates the equilibrium speed in sequence. At least 2 unsteady calculations are required, and the calculation time required for an unsteady flow field is about 12 times that of a steady flow field, which is equivalent to more than 20 steady calculations.

[0053] In view of the above technical problems, the present application proposes a method for estimating the equilibrium rotation speed of a projectile with an inclined tail fin, such as Figure 1 As shown, the specific steps are as follows:

[0054] Step S1, constructing a projectile geometric model according to preset projectile data and preset tail wing tilt angle.

[0055] In an embodiment of the present application, the preset projectile data includes data such as the length, diameter, head shape, number and size of the tail fins of the projectile. Select a suitable three-dimensional coordinate system (such as a Cartesian coordinate system) to determine the origin (usually the geometric center or tail of the projectile); then use geometric modeling software to draw the main part of the projectile based on the basic size data of the projectile. According to the preset number and size of the tail fins and the tail fin tilt angle and other data, add a tail fin to the tail of the projectile. Among them, the tilt angle of the tail fin can ensure that each tail fin is tilted at a preset angle to ensure the stability of the projectile during flight. Finally, use the measurement tool of the geometric modeling software to verify whether the size of the projectile and the tilt angle of the tail fin meet the preset data, and adjust as needed until the design requirements are met. Finally, a constructed projectile geometric model is obtained.

[0056] Step S2, based on the projectile geometric model, perform steady calculation according to preset working condition data to obtain surface flow field data and spatial flow field data.

[0057] In an embodiment of the present application, the geometric model of the projectile is gridded to obtain a computational fluid dynamics grid; the height and speed of the projectile are set according to preset working condition data; and based on the height and speed of the projectile, the surface flow field data and spatial flow field data of the projectile are determined using the computational fluid dynamics grid.

[0058] Specifically, mesh the geometric model of the projectile using mesh generation software and divide it into a computational fluid dynamics grid file. Then select a suitable CFD solver and import the generated mesh. Collect and prepare all necessary working condition data, set working condition data for the projectile, and the working condition data includes data such as height and speed. For example, choose to estimate the equilibrium speed at an altitude of 5000 meters and a speed of 1000m / s. Perform CFD steady-state calculations under given working conditions, start the solution process, monitor convergence, and adjust the iteration parameters if necessary until a stable solution is obtained. After the calculation converges, extract surface flow field data (such as pressure distribution, temperature, etc.) and spatial flow field data (such as velocity field, temperature field, pressure field, velocity components in each direction, etc.).

[0059] Step S3, obtaining the average value of the coordinates of the pressure center or the centroid on each tail wing of the projectile and the average value of the radius of the pressure center or the centroid from the projectile axis according to the surface flow field data.

[0060] In an embodiment of the present application, the pressure center or centroid coordinates of each tail wing on the projectile are extracted based on the surface flow field data; and based on the pressure center or centroid coordinates of each tail wing, the average value of the pressure center or centroid coordinates and the average value of the radius of the pressure center or centroid from the projectile axis are obtained.

[0061] Specifically, in the post-processing software Tecplot, first, the tail area is identified, and the pressure center or centroid coordinates (X, Y, Z) on each tail are extracted. The coordinate X axis coincides with the missile axis, such as Figure 2 As shown, it can be forward or backward; the coordinate Y axis and Z axis are perpendicular to the X axis respectively. Among them, the pressure center is the equivalent point of action of the aerodynamic force on the object, which is usually calculated by integrating the pressure distribution; the centroid is the "geometric center" of the geometric body, which is obtained by calculating the center of mass of the geometric shape. Secondly, the X-axis of the pressure center or centroid coordinates of each tail wing are summed, and then divided by the number of tail wing, to obtain the average value of the pressure center or centroid coordinates. Then, according to the coordinates or direction of the missile axis, calculate the distance from the pressure center or centroid of each tail wing to the missile axis; sum these distances, and then according to the tail wing data, obtain the average value of the radius of the pressure center or centroid from the missile axis.

[0062] For example: If you have 4 tail fins, and the coordinates of the center of pressure of each tail fin are (x1, y1, z1), (x2, y2, z2), (x3, y3, z3) and (x4, y4, z4), then the average coordinate is ((x1+x2+x3+x4) / 4. If the distances from the center of pressure of the 4 tail fins to the missile axis are d1, d2, d3 and d4 respectively, then the average radius is (d1+d2+d3+d4) / 4.

[0063] Step S4, obtaining an average velocity value according to the spatial flow field data, the average value of the pressure center or the centroid coordinates and the average value of the radius of the pressure center or the centroid from the projectile axis.

[0064] In an embodiment of the present application, the isosurface value is obtained based on the spatial flow field data and the average value of the radius of the pressure center or centroid from the elastic axis; the isoline value is obtained based on the isosurface value and the average value of the pressure center or centroid coordinates; and the average value of the velocity on the isovalue line is obtained based on the isosurface value and the isoline value.

[0065] Specifically, using spatial flow field data, using the post-processing software Tecplot, set the distance from the elastic axis to be equal to the average value of the radius of the pressure center or the centroid, thereby generating an isosurface. Among them, an isosurface refers to a surface with equal values ​​on a certain specific variable (such as velocity, pressure, etc.), that is, an isosurface formed by variables such as velocity and pressure in the flow field. Then set the x coordinate to be equal to the average value of the x coordinate of the pressure center or the centroid, thereby generating an isoline on the isosurface. Among them, the isoline value refers to a line with equal variable values ​​in a certain specific direction (such as a direction perpendicular to the isosurface) on the isosurface. Finally, extract the velocity data of each point on the isosurface, and obtain the velocity average value based on the velocity data of each point.

[0066] Step S5, obtaining the equilibrium rotation speed of the projectile according to the average value of the radius of the pressure center or the centroid from the projectile axis, the average speed and the preset tail fin tilt angle.

[0067] In the embodiment of the present application, the average speed u is used * Multiply by the sine value of the tail wing tilt angle δ, divide the product by the average value r* of the radius from the pressure center or the centroid to the missile axis, and calculate the equilibrium speed;

[0068] ω x =u * *sin(δ) / r*

[0069] Among them, ω x Indicates the equilibrium speed;

[0070] x represents the coordinate;

[0071] u * Indicates the average speed;

[0072] δ represents the tail wing tilt angle;

[0073] r* represents the average value of the radius from the compression center or centroid to the projectile axis;

[0074] The present application can estimate the equilibrium speed without the need for prior knowledge, and is applicable to any inclined tail fin projectile model and any working condition, with a wider range of applicable working conditions than the engineering estimation method. The present application uses the processing results of the static flow field in the inclined tail fin projectile model, which has a more reliable flow field data basis and higher result accuracy than the engineering estimation method. In the present application, only one steady-state CFD calculation is required for the inclined tail fin projectile model, which requires much less calculation time than the interpolation calculation method, and can also estimate the equilibrium speed more quickly.

[0075] The embodiment of the present application provides a device for estimating the equilibrium rotation speed of a projectile with an inclined tail fin, such as Figure 3 As shown, including:

[0076] A construction module 101 is used to construct a projectile geometric model according to preset projectile data and preset tail fin tilt angle;

[0077] A data acquisition module 102 is used to perform steady-state calculations based on the projectile geometric model and preset working condition data to obtain surface flow field data and spatial flow field data;

[0078] The first calculation module 103 is used to obtain the average value of the pressure center or the centroid coordinates on each tail wing of the projectile and the average value of the radius of the pressure center or the centroid from the projectile axis according to the surface flow field data;

[0079] A second calculation module 104 is used to obtain a velocity average value according to the spatial flow field data, an average value of the pressure center or the centroid coordinates, and an average value of the radius of the pressure center or the centroid from the projectile axis; and

[0080] The third calculation module 105 is used to obtain the equilibrium rotation speed of the projectile according to the average value of the radius of the pressure center or the centroid from the projectile axis, the average speed and the preset tail fin tilt angle.

[0081] Furthermore, the data acquisition module is configured as follows:

[0082] Meshing the projectile geometric model to obtain a computational fluid dynamics grid;

[0083] Set the height and speed of the projectile according to the preset working condition data; and

[0084] According to the height and speed of the projectile, the surface flow field data and the spatial flow field data of the projectile are determined by using a computational fluid dynamics grid.

[0085] Furthermore, the first computing module is configured to:

[0086] Extracting the pressure center or centroid coordinates of each tail fin on the projectile according to the surface flow field data; and

[0087] According to the pressure center or centroid coordinates of each tail wing, an average value of the pressure center or centroid coordinates and an average value of the radius of the pressure center or centroid from the missile axis are obtained.

[0088] Furthermore, the second computing module is configured to:

[0089] Obtaining an isosurface value according to the spatial flow field data and an average value of the radius of the pressure center or the centroid from the projectile axis;

[0090] Obtaining an isoline value according to the isosurface value and the average value of the pressure center or centroid coordinates; and

[0091] According to the isosurface value and the isoline value, an average value of the velocity on the isoline is obtained.

[0092] An embodiment of the present application provides a storage medium for storing computer executable instructions, which, when executed, implements the steps of the method for estimating the equilibrium rotation speed of a projectile with an oblique tail fin described in any one of the above embodiments.

[0093] It should be noted that the embodiments of the storage medium in this specification and the embodiments of the blockchain-based service provision method in this specification are based on the same inventive concept. Therefore, the specific implementation of this embodiment can refer to the aforementioned corresponding implementation of the blockchain-based service provision method, and the repeated parts will not be repeated.

[0094] The above is a description of a specific embodiment of the specification. Other embodiments are within the scope of the appended claims. In some cases, the actions or steps recorded in the claims can be performed in an order different from that in the embodiments and still achieve the desired results. In addition, the processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.

[0095] In the 1930s, improvements to a technology could be clearly distinguished as hardware improvements (for example, improvements to the circuit structure of diodes, transistors, switches, etc.) or software improvements (improvements to the method flow). However, with the development of technology, many improvements to the method flow today can be regarded as direct improvements to the hardware circuit structure. Designers almost always obtain the corresponding hardware circuit structure by programming the improved method flow into the hardware circuit. Therefore, it cannot be said that an improvement in a method flow cannot be implemented using a hardware entity module. For example, a programmable logic device (PLD) (such as a field programmable gate array (FPGA)) is such an integrated circuit whose logical function is determined by the user's programming of the device. Designers can "integrate" a digital system on a PLD by programming it themselves, without having to ask a chip manufacturer to design and produce a dedicated integrated circuit chip. Moreover, nowadays, instead of manually making integrated circuit chips, this kind of programming is mostly implemented by "logic compiler" software, which is similar to the software compiler used when developing and writing programs, and the original code before compilation must also be written in a specific programming language, which is called hardware description language (HDL). There is not only one HDL, but many kinds, such as ABEL (Advanced Boolean Expression Language), AHDL (Altera Hardware Description Language), Confluence, CUPL (Cornell University Programming Language), HDCal, JHDL (Java Hardware Description Language), Lava, Lola, MyHDL, PALASM, RHDL (Ruby Hardware Description Language), etc. The most commonly used ones are VHDL (Very-High-Speed ​​Integrated Circuit Hardware Description Language) and Verilog. Those skilled in the art should also know that it is only necessary to program the method flow slightly in the above-mentioned hardware description languages ​​and program it into the integrated circuit, and then it is easy to obtain the hardware circuit that implements the logic method flow.

[0096] The controller can be implemented in any appropriate manner, for example, the controller can take the form of a microprocessor or processor and a computer-readable medium storing a computer-readable program code (such as software or firmware) that can be executed by the (micro)processor, a logic gate, a switch, an application-specific integrated circuit (ASIC), a programmable logic controller, and an embedded microcontroller. Examples of controllers include, but are not limited to, the following microcontrollers: ARC 625D, Atmel AT91SAM, Microchip PIC18F26K20, and Silicone Labs C8051F320. The memory controller can also be implemented as part of the control logic of the memory. Those skilled in the art also know that in addition to implementing the controller in a purely computer-readable program code manner, the controller can be implemented in the form of a logic gate, a switch, an application-specific integrated circuit, a programmable logic controller, and an embedded microcontroller by logically programming the method steps. Therefore, this controller can be considered as a hardware component, and the devices included therein for implementing various functions can also be regarded as structures within the hardware component. Or even, the devices for implementing various functions can be regarded as both software modules for implementing the method and structures within the hardware component.

[0097] The systems, devices, modules or units described in the above embodiments may be implemented by computer chips or entities, or by products with certain functions. A typical implementation device is a computer. Specifically, the computer may be, for example, a personal computer, a laptop computer, a cellular phone, a camera phone, a smart phone, a personal digital assistant, a media player, a navigation device, an email device, a game console, a tablet computer, a wearable device, or a combination of any of these devices.

[0098] For the convenience of description, the above devices are described in terms of functions and are divided into various units. Of course, when implementing the embodiments of this specification, the functions of each unit can be implemented in the same or multiple software and / or hardware.

[0099] It should be understood by those skilled in the art that one or more embodiments of this specification may be provided as a method, system or computer program product. Therefore, one or more embodiments of this specification may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, this specification may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0100] This specification is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of this specification. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0101] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0102] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0103] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0104] The memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0105] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer readable media does not include temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0106] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0107] One or more embodiments of the present specification may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform specific tasks or implement specific abstract data types. One or more embodiments of the present specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules may be located in local and remote computer storage media, including storage devices.

[0108] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0109] The above description is only an embodiment of this document and is not intended to limit this document. For those skilled in the art, this document may have various changes and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of this document should be included in the scope of the claims of this document.

Claims

1. A method for estimating the equilibrium rotation speed of a projectile with an inclined tail fin, characterized in that include: Constructing a projectile geometry model according to preset projectile data and preset tail wing tilt angle; Based on the projectile geometric model, a steady calculation is performed according to preset working condition data to obtain surface flow field data and spatial flow field data; According to the surface flow field data, an average value of the coordinates of the pressure center or the centroid on each tail wing of the projectile and an average value of the radius of the pressure center or the centroid from the projectile axis are obtained; Obtaining an average velocity value according to the spatial flow field data, an average value of the pressure center or the centroid coordinates, and an average value of the radius of the pressure center or the centroid from the projectile axis; and The equilibrium rotation speed of the projectile is obtained according to the average value of the radius of the pressure center or the centroid from the projectile axis, the average speed and the preset tail fin tilt angle.

2. The method according to claim 1, characterized in that The surface flow field data and the spatial flow field data are obtained by performing steady calculation based on the projectile geometric model according to the preset working condition data, including: Meshing the projectile geometric model to obtain a computational fluid dynamics grid; Set the height and speed of the projectile according to the preset working condition data; and According to the height and speed of the projectile, the surface flow field data and the spatial flow field data of the projectile are determined by using a computational fluid dynamics grid.

3. The method according to claim 1, characterized in that The average value of the pressure center or centroid coordinates on each tail wing of the projectile and the average value of the radius of the pressure center or centroid from the projectile axis obtained according to the surface flow field data include: Extracting the pressure center or centroid coordinates of each tail fin on the projectile according to the surface flow field data; and According to the pressure center or centroid coordinates of each tail wing, an average value of the pressure center or centroid coordinates and an average value of the radius of the pressure center or centroid from the missile axis are obtained.

4. The method according to claim 1, characterized in that The velocity average value obtained according to the spatial flow field data, the average value of the pressure center or the centroid coordinates and the average value of the radius of the pressure center or the centroid from the missile axis includes: Obtaining an isosurface value according to the spatial flow field data and an average value of the radius from the pressure center or the centroid to the projectile axis; Obtaining an isoline value according to the isosurface value and the average value of the pressure center or centroid coordinates; and According to the isosurface value and the isoline value, an average value of the velocity on the isoline is obtained.

5. The method according to claim 1, characterized in that The calculation formula of the equilibrium rotation speed of the projectile is as follows: oh x =u * *sin(δ) / r* Among them, ω x Indicates the equilibrium speed; x represents the coordinate axis; u * Indicates the average speed; δ represents the tail wing cant angle; and r* represents the average value of the radius from the compression center or centroid to the projectile axis.

6. A device for estimating the equilibrium rotation speed of a projectile with an inclined tail fin, characterized in that include: A construction module, used for constructing a projectile geometric model according to preset projectile data and preset tail fin tilt angle; A data acquisition module, used to perform steady calculation based on the projectile geometric model and preset working condition data to obtain surface flow field data and spatial flow field data; A first calculation module is used to obtain the average value of the coordinates of the pressure center or the centroid on each tail wing of the projectile and the average value of the radius of the pressure center or the centroid from the projectile axis according to the surface flow field data; A second calculation module is used to obtain an average velocity value according to the spatial flow field data, an average value of the pressure center or the centroid coordinates, and an average value of the radius of the pressure center or the centroid from the missile axis; as well as The third calculation module is used to obtain the equilibrium rotation speed of the projectile according to the average value of the radius of the pressure center or the centroid from the projectile axis, the average speed and the preset tail fin tilt angle.

7. The device according to claim 6, characterized in that The data acquisition module is configured as follows: Meshing the projectile geometric model to obtain a computational fluid dynamics grid; Set the height and speed of the projectile according to the preset working condition data; and According to the height and speed of the projectile, the surface flow field data and the spatial flow field data of the projectile are determined by using a computational fluid dynamics grid.

8. The device according to claim 6, characterized in that The first computing module is configured to: Extracting the pressure center or centroid coordinates of each tail fin on the projectile according to the surface flow field data; and According to the pressure center or centroid coordinates of each tail wing, an average value of the pressure center or centroid coordinates and an average value of the radius of the pressure center or centroid from the missile axis are obtained.

9. The device according to claim 6, characterized in that The second computing module is configured to: Obtaining an isosurface value according to the spatial flow field data and an average value of the radius from the pressure center or the centroid to the projectile axis; Obtaining an isoline value according to the isosurface value and the average value of the pressure center or centroid coordinates; and According to the isosurface value and the isoline value, an average value of the velocity on the isoline is obtained.

10. A storage medium for storing computer executable instructions, characterized in that: When the computer executable instructions are executed, the steps of the method for estimating the equilibrium rotation speed of a projectile with an oblique tail fin according to any one of claims 1 to 5 are implemented.

Images