CN109781095B - Method for predicting motion trail of underwater vehicle based on pressure distribution - Google Patents
Method for predicting motion trail of underwater vehicle based on pressure distribution Download PDFInfo
- Publication number
- CN109781095B CN109781095B CN201910071175.4A CN201910071175A CN109781095B CN 109781095 B CN109781095 B CN 109781095B CN 201910071175 A CN201910071175 A CN 201910071175A CN 109781095 B CN109781095 B CN 109781095B
- Authority
- CN
- China
- Prior art keywords
- navigation body
- formula
- deflection angle
- pressure
- absolute
- 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
Images
Landscapes
- Navigation (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
The invention relates to a method for predicting a motion trail of an underwater vehicle based on pressure distribution, belonging to the technical field of ship and underwater vehicle engineering and underwater vehicle motion trail estimation. The method is based on an underwater vehicle test, pressure measuring points are arranged on the surface of the vehicle, pressure data of the pressure measuring points on the surface of the vehicle are obtained through the test, and the motion rule and the track of the underwater vehicle are finally obtained through further research and calculation based on the data. The invention can provide scientific basis for the technologies of attitude control, motion mode control and the like of the underwater vehicle and provide technical support for the pressure measuring point arrangement scheme of the developed underwater vehicle test.
Description
Technical Field
The invention relates to a method for predicting a motion trail of an underwater vehicle based on pressure distribution, belonging to the technical field of ship and underwater vehicle engineering and underwater vehicle motion trail estimation.
Background
With the development of modern ocean engineering and shipbuilding industry, the navigation safety of underwater vehicles has been more and more generally concerned and emphasized, wherein the prediction of the motion trail of the underwater vehicle is one of the important points of concern. The motion track of the navigation body is restricted by the complex underwater environment, so that the water outlet position, the motion area and the like of the navigation body are changed, and a decision maker is difficult to control. Therefore, the method for predicting the motion trail of the underwater vehicle is a brand new technical problem and has more practical engineering value and scientific significance.
Since the 80's in the 20 world, with the development of computer equipment and the advancement of computing technology, the development of computational fluid dynamics has further pushed the research on underwater vehicle problems. From traditional experimental research, students turn to a research mode combining numerical values and experiments, and mainly research the aspects of the drag reduction characteristics of the navigation body, the mechanical properties of the navigation body and the like. For the research on the prediction of the motion track of the underwater vehicle, a scholars proposes to estimate the track according to the factors of sea waves and ocean currents. However, the method for estimating the motion trail of the navigation body according to the external change rule is not strict, on one hand, the method heavily researches the change rule of the water body and ignores the influence on the navigation body; on the other hand, the method only considers the influence factors of sea waves and ocean currents, and ignores other restriction conditions generated by the complex ocean environment.
Disclosure of Invention
The invention aims to provide a method for predicting the motion trail of an underwater vehicle based on pressure distribution, which can realize the prediction of the motion trail of the underwater vehicle and further solve the problem of the prediction of the motion trail of the underwater vehicle.
The purpose of the invention is realized by the following technical scheme.
A method for predicting a motion track of an underwater vehicle based on pressure distribution specifically comprises the following steps:
step one, arranging pressure measuring points on the surface of an aeronautical body;
the pressure measurement points are shown in table 1;
step two, acquiring pressure data P of the pressure measuring points of the upstream surface and the downstream surface of the navigation body in the direction of a relative coordinate system x' in the step oneaiAnd Pbi;
Step three, according to the pressure data of the upstream surface and the back surface of the navigation body in the step two, calculating the pressure difference delta P of the pressure measuring point data of the upstream surface and the back surface of the navigation body by using a formula (1)x’i(ii) a Carrying out differential pressure surface integral on the measuring points by using a formula (2), and then calculating to obtain the resultant force F of the measuring points in the direction of a relative coordinate system xx’iFinally, the resultant force of each measuring point in the direction of the relative coordinate system x' is summed by using a formula (3) to obtain the resultant force F of the navigation body in the direction of the relative coordinate system xti;
ΔPx,i=Pai-Pbi(1)
Fx,i=0.6ΔPx,iΔSx,(2)
Wherein the constant 0.6 is the pressure integral correction coefficient, Delta Sx’Is the cross-sectional area of the navigation body.
Fourthly, based on the absolute coordinate system, according to the surface resultant force F of the sailing body obtained in the third steptiCombining the rotary inertia J of the navigation body in the formula (4) and the torque M of the navigation body in the formula (5), the angular acceleration alpha of the navigation body is calculated by using the formula (6)i;
Wherein m is the navigation body mass, l is the navigation body length, diThe distance between the pressure measuring point and the centroid is shown.
Step five, obtaining the angular acceleration alpha through the step fouriObtaining the relative deflection angle delta theta of the navigation body by using the formula (7)iThe absolute deflection angle theta can be calculated by combining the relative deflection angle and using the formula (8)iThe absolute deflection angle is the sum of the relative deflection angle at the moment and the absolute deflection angle at the last moment;
Δθi=αiΔt2(7)
θi=Δθi+θi-1(8)
where Δ t is the navigation body calculation time step, θi-1Is the absolute deflection angle of the vehicle at a moment.
Step six, forecasting the speed of the navigation body at the next moment by combining the absolute deflection angle of the navigation body obtained in the step five through the initial speeds of the absolute coordinate system of the navigation body along the y direction and the x direction and by using a formula (9) and a formula (10);
wherein v isx(ti)And vy(ti)Respectively the last moment speed of the navigation body.
Step seven: integrating the speed of the navigation body predicted in the sixth step in time by using a formula (11) and a formula (12) to obtain the relative displacement of the navigation body in the time step and the absolute displacement of the navigation body in the last momentAdding up to obtain the absolute displacement of the predicted time of the navigation body
Repeating the fifth step to the seventh step, and carrying out time iteration on the deflection angle and the speed of the navigation body to obtain the displacement of the next moment;
step eight: and fitting the motion trail of the navigation body according to the displacement obtained by the navigation body in the step seven, namely realizing the prediction of the motion trail of the navigation body.
TABLE 1 arrangement of pressure measuring point combinations of a navigation body
Advantageous effects
1. According to the method for predicting the motion trail of the underwater vehicle based on the pressure distribution, the mechanical properties of the vehicle are fully considered, and the prediction accuracy is improved;
2. the method for predicting the motion trail of the underwater vehicle based on the pressure distribution can provide an arrangement scheme for the arrangement of pressure measuring points in a subsequent underwater vehicle test, and saves the test cost and time;
drawings
FIG. 1 is a flow chart of a calculation method for predicting a motion trajectory of an underwater vehicle based on pressure distribution according to the present invention;
FIG. 2 is a schematic diagram of a navigation body coordinate system, upstream surface and downstream surface pressure measurement point distribution and related dimensions;
FIG. 3 is a graph of the variation of the upstream surface pressure of the vehicle in case 2;
FIG. 4 is a graph of the variation of the pressure of the back surface of the vehicle in case 2;
FIG. 5 is a graph showing the variation of the pressure difference between the upstream surface and the downstream surface of the vehicle in the case of the scheme 2;
FIG. 6 is a graph showing the resultant force variation of the surface measurement points of the sailing body in the embodiment 2;
FIG. 7 is a graph of variation of angular acceleration of the sailing body according to scheme 2;
FIG. 8 is a graph of absolute yaw angle change for a solution 2 vehicle;
FIG. 9 is a graph of variation in speed of vehicle movement for scenario 2;
FIG. 10 is a diagram comparing the prediction and experiment of the motion trail of the navigation body in the scheme 2;
FIG. 11 is a diagram of the prediction and experimental comparison of the trajectory of the vehicle in case 10;
fig. 12 is a diagram for comparing the prediction and experiment of the motion trail of the navigation body with 13 schemes.
Detailed Description
The following detailed description of embodiments of the invention refers to the accompanying drawings.
Example 1
The navigation body is positioned in a 5-level wave water body environment with the water depth of 20m, the initial speeds in the x direction and the y direction are respectively 2m/s and 20m/s, and the adopted pressure measuring point arrangement scheme is scheme 2.
A specific flow of a calculation method for predicting a motion trajectory of an underwater vehicle based on pressure distribution is mainly realized by the following steps as shown in FIG. 1:
the method comprises the following steps: arranging a surface pressure measuring point of the navigation body;
based on an underwater vehicle test, pressure measuring points are arranged on the surface of a vehicle, the number of single-side pressure measuring points of the vehicle is 52, the numbers are gradually increased from the head to the tail, the numbers of the pressure measuring points on the upstream surface and the downstream surface are the same, and the total number of the measuring points is 102, as shown in FIG. 2; the distribution scheme of the combination of the pressure measuring points of the navigation body is shown in the table 1;
step two: acquiring surface pressure measuring point data of a navigation body;
establishing a relative coordinate system and an absolute coordinate system of the navigation body, as shown in fig. 2; acquiring pressure data of a pressure measuring point on the surface of a navigation body relative to the upstream surface and the back flow surface in the x' direction of a coordinate system, wherein the pressure data are PaiAnd PbiAs shown in fig. 3 and 4;
step three: calculating the surface measuring point resultant force of the sailing body;
according to the pressure data of the upstream surface and the back surface of the navigation body in the step two, the pressure difference delta P of the pressure measuring point data of the upstream surface and the back surface of the navigation body is obtained by using the formula (1)x’iAs shown in fig. 5; carrying out differential pressure surface integral on the measuring points by using a formula (2), and then calculating to obtain the resultant force F of the measuring points in the direction of a relative coordinate system xx’i(ii) a Finally, the resultant force of each measuring point in the direction of the relative coordinate system x' is summed by using a formula (3) to obtain the resultant force F of the navigation body in the direction of the relative coordinate system xtiAs shown in fig. 6;
ΔPx,i=Pai-Pbi(1)
Fx,i=0.6ΔPx,iΔSx,(2)
wherein the constant 0.6 is the integral correction of pressurePositive coefficient, Δ Sx’Is the cross-sectional area of the navigation body.
Step four: calculating the angular acceleration of the navigation body;
based on the absolute coordinate system, according to the surface resultant force F of the sailing body in the third steptiCombining the rotary inertia J of the navigation body in the formula (4) and the torque M of the navigation body in the formula (5), the angular acceleration alpha of the navigation body is calculated by using the formula (6)iAs shown in fig. 7;
wherein m is the navigation body mass, l is the navigation body length, diThe distance between the pressure measuring point and the centroid is shown.
Step five: calculating the absolute deflection angle of the navigation body;
angular acceleration alpha obtained by the fourth stepiObtaining the relative deflection angle delta theta of the navigation body by using the formula (7)iThe absolute deflection angle theta can be calculated by combining the relative deflection angle and using the formula (8)iAs shown in fig. 8, the absolute deflection angle is the sum of the relative deflection angle at the present moment and the absolute deflection angle at the previous moment;
Δθi=αiΔt2(7)
θi=Δθi+θi-1(8)
step six: calculating the movement speed of the navigation body;
knowing that the initial speeds of the absolute coordinate system of the navigation body along the y direction and the x direction are 20m/s and 2m/s respectively, predicting the speed of the navigation body at the next moment by using a formula (9) and a formula (10) in combination with the absolute deflection angle of the navigation body obtained in the step five, and showing in fig. 9;
wherein v isx(ti)And vy(ti)Respectively the last moment speed of the navigation body.
Step seven: calculating the displacement of the navigation body;
integrating the speed of the navigation body obtained in the step six in time by using a formula (11) and a formula (12), so as to obtain the relative displacement of the navigation body in the time step and the absolute displacement of the navigation body in the last momentAdding up to obtain the absolute displacement of the navigation body at the moment
Repeating the fifth step to the seventh step, and carrying out time iteration on the deflection angle and the speed of the navigation body to obtain the displacement of the next moment;
step eight: calculating the motion trail of the navigation body;
and according to the displacement of the navigation body in the sixth step, fitting to obtain a motion track of the navigation body, namely realizing the prediction of the motion track of the navigation body, as shown in fig. 10.
Example 2
The navigation body is positioned in a 5-level wave water body environment with the water depth of 20m, the initial speeds in the x direction and the y direction are respectively 2m/s and 20m/s, and the scheme of the arrangement scheme of the pressure measuring points is a scheme 10.
The specific flow of the calculation method for predicting the motion trail of the underwater vehicle based on the pressure distribution is consistent with that of case 1, and the predicted motion trail is shown in fig. 11.
Based on the pressure measuring point combination arrangement scheme of the table 1, the motion trail prediction of 13 groups of schemes is given, as shown in fig. 12. The deviation error of the combined distribution of the 13 pressure measuring points in the table 1 is given in the table 2 by taking the distance measurement 1m at the position of 20m of the water surface as the judgment standard of the deviation error.
Deviation error of distribution scheme of 213 pressure measuring points in table
The above description is only for the purpose of illustrating the present invention, and modifications and equivalents thereof may be made by those skilled in the art. All changes, equivalents, modifications and the like which come within the spirit and principle of the invention are desired to be protected.
Claims (1)
1. A method for predicting the motion trail of an underwater vehicle based on pressure distribution is characterized by comprising the following steps: the method comprises the following specific steps:
step one, arranging pressure measuring points on the surface of an aeronautical body;
step two, acquiring pressure data P of the pressure measuring points of the upstream surface and the downstream surface of the navigation body in the direction of a relative coordinate system x' in the step oneaiAnd Pbi;
Step three, according to the pressure data of the upstream surface and the back surface of the navigation body in the step two, calculating the pressure difference delta P of the pressure measuring point data of the upstream surface and the back surface of the navigation body by using a formula (1)x’i(ii) a Carrying out differential pressure surface integral on the measuring points by using a formula (2), and then calculating to obtain the resultant force F of the measuring points in the direction of a relative coordinate system xx’iFinally, the resultant force of each measuring point in the direction of the relative coordinate system x' is summed by using a formula (3) to obtain the resultant force F of the navigation body in the direction of the relative coordinate system xti;
ΔPx,i=Pai-Pbi(1)
Fx’i=0.6ΔPx’iΔSx’(2)
Wherein the constant 0.6 is the pressure integral correction coefficient, Delta Sx’Is the cross-sectional area of the navigation body;
fourthly, based on the absolute coordinate system, according to the surface resultant force F of the sailing body obtained in the third steptiCombining the rotary inertia J of the navigation body in the formula (4) and the torque M of the navigation body in the formula (5), the angular acceleration alpha of the navigation body is calculated by using the formula (6)i;
Wherein m is the navigation body mass, l is the navigation body length, diThe distance between the pressure measuring point and the centroid is taken as the distance;
step five, obtaining the angular acceleration alpha through the step fouriObtaining the relative deflection angle delta theta of the navigation body by using the formula (7)iThe absolute deflection angle theta can be calculated by combining the relative deflection angle and using the formula (8)iThe absolute deflection angle is the sum of the relative deflection angle at the moment and the absolute deflection angle at the last moment;
Δθi=αiΔt2(7)
θi=Δθi+θi-1(8)
where Δ t is the navigation body calculation time step, θi-1The absolute deflection angle of the navigation body at the last moment;
step six, forecasting the speed of the navigation body at the next moment by combining the absolute deflection angle of the navigation body obtained in the step five through the initial speeds of the absolute coordinate system of the navigation body along the y direction and the x direction and by using a formula (9) and a formula (10);
wherein v isx(ti)And vy(ti)Respectively the speed of the navigation body at the last moment;
step seven: integrating the speed of the navigation body predicted in the sixth step in time by using a formula (11) and a formula (12) to obtain the relative displacement of the navigation body in the time step and the absolute displacement of the navigation body in the last momentAdding up to obtain the absolute displacement of the predicted time of the navigation body
Repeating the fifth step to the seventh step, and carrying out time iteration on the deflection angle and the speed of the navigation body to obtain the displacement of the next moment;
step eight: and fitting the motion trail of the navigation body according to the displacement obtained by the navigation body in the step seven, namely realizing the prediction of the motion trail of the navigation body.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910071175.4A CN109781095B (en) | 2019-01-25 | 2019-01-25 | Method for predicting motion trail of underwater vehicle based on pressure distribution |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910071175.4A CN109781095B (en) | 2019-01-25 | 2019-01-25 | Method for predicting motion trail of underwater vehicle based on pressure distribution |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109781095A CN109781095A (en) | 2019-05-21 |
CN109781095B true CN109781095B (en) | 2020-10-27 |
Family
ID=66502525
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910071175.4A Active CN109781095B (en) | 2019-01-25 | 2019-01-25 | Method for predicting motion trail of underwater vehicle based on pressure distribution |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109781095B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112363394B (en) * | 2020-11-02 | 2022-04-15 | 中国科学院力学研究所 | Ventilation supercavitation form control method based on machine learning |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2442718C1 (en) * | 2010-09-13 | 2012-02-20 | Федеральное государственное образовательное учреждение высшего профессионального образования "Мурманский государственный технический университет" (ФГОУВПО "МГТУ") | Determination method for dynamic parameters of marine movement mathematical model |
CN105333988A (en) * | 2015-11-25 | 2016-02-17 | 中国海洋大学 | Artificial lateral line pressure detection method |
CN106564577A (en) * | 2016-11-02 | 2017-04-19 | 中国海洋大学 | Multifunctional AUV based on bionic lateral line |
CN107084723A (en) * | 2017-05-12 | 2017-08-22 | 中国人民解放军91550部队 | A kind of underwater sailing body movement locus method of estimation under marine environment |
-
2019
- 2019-01-25 CN CN201910071175.4A patent/CN109781095B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2442718C1 (en) * | 2010-09-13 | 2012-02-20 | Федеральное государственное образовательное учреждение высшего профессионального образования "Мурманский государственный технический университет" (ФГОУВПО "МГТУ") | Determination method for dynamic parameters of marine movement mathematical model |
CN105333988A (en) * | 2015-11-25 | 2016-02-17 | 中国海洋大学 | Artificial lateral line pressure detection method |
CN106564577A (en) * | 2016-11-02 | 2017-04-19 | 中国海洋大学 | Multifunctional AUV based on bionic lateral line |
CN107084723A (en) * | 2017-05-12 | 2017-08-22 | 中国人民解放军91550部队 | A kind of underwater sailing body movement locus method of estimation under marine environment |
Non-Patent Citations (1)
Title |
---|
基于扰动压力与速度测量水中目标航向与航速的新方法;顾建农等;《舰船科学技术》;20120531;第34卷(第5期);47-50 * |
Also Published As
Publication number | Publication date |
---|---|
CN109781095A (en) | 2019-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108846207B (en) | Ship motion forecasting device and forecasting method based on nonlinear model identification | |
CN110032197B (en) | Unmanned ship path tracking control method based on limited time guidance and control | |
CN105807789B (en) | UUV control methods based on the compensation of T-S Fuzzy Observers | |
CN109901402B (en) | Autonomous underwater robot path tracking method based on course smoothing technology | |
CN109946976B (en) | Wide-navigational-speed AUV motion control method | |
CN109781095B (en) | Method for predicting motion trail of underwater vehicle based on pressure distribution | |
Ueno et al. | Rudder effectiveness correction for scale model ship testing | |
CN112051742A (en) | MPC-based full-drive ship track tracking method | |
CN111798702A (en) | Unmanned ship path tracking control method, system, storage medium and terminal | |
CN112327883B (en) | Under-actuated surface vessel path tracking control method based on extended state observer | |
CN117389312A (en) | Model-based three-dimensional tracking control method for counter roll of underwater vehicle | |
CN117369252A (en) | Heterogeneous collaborative track tracking method for unmanned ship on water surface and unmanned underwater ship | |
Kobets et al. | Control of the Pivot Point Position of a Conventional Single-Screw Vessel. | |
CN115480580A (en) | NMPC-based underwater robot path tracking and obstacle avoidance control method | |
CN114706320A (en) | DP ship thrust optimal distribution simulation method under multi-constraint condition | |
CN110837254B (en) | Longitudinal motion control method for hydrofoil craft | |
Ahmed et al. | Nonlinear dynamics of novel flight-style autonomous underwater vehicle with bow wings, Part I: ASE and CFD based estimations of hydrodynamic coefficients, Part II: Nonlinear dynamic modeling and experimental validations | |
CN102508819B (en) | Angular-speed-based quaternion Legendre approximate output method during extreme flying of aircraft | |
CN102495831B (en) | Quaternion Hermitian approximate output method based on angular velocities for aircraft during extreme flight | |
CN103336525B (en) | Stochastic system high weight convenient UKF filtering method | |
Suzuki et al. | Motion simulation of an underwater vehicle with mechanical pectoral fins using a CFD-based motion simulator | |
CN102495830B (en) | Quaternion Hartley approximate output method based on angular velocities for aircraft during extreme flight | |
Dantas et al. | AUV control in the diving plane subject to waves | |
Viswanath et al. | Separated trailing-edge flow at a transonic Mach number | |
CN102506866B (en) | Angle speed-based Chebyshev approximate output method of quaternion numbers in ultimate flight of aircraft |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |