CN110834703A - Analysis method for simulating relationship between change of deck roll inclination angle along with external force - Google Patents
Analysis method for simulating relationship between change of deck roll inclination angle along with external force Download PDFInfo
- Publication number
- CN110834703A CN110834703A CN201911115932.XA CN201911115932A CN110834703A CN 110834703 A CN110834703 A CN 110834703A CN 201911115932 A CN201911115932 A CN 201911115932A CN 110834703 A CN110834703 A CN 110834703A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- external force
- ship
- inclination angle
- model
- 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.)
- Pending
Links
Images
Abstract
The invention relates to a method for analyzing the relation of a simulation deck roll inclination angle along with external force change. The method specifically comprises the following steps: the method comprises the following steps: a ship body model with reasonable design and an ideal experimental environment fix the ship bottom in a positive floating state, so that only rolling motion occurs. Step two: establishing a geodetic coordinate system OXYZ and a hull coordinate system OX1Y1Z1So as to more clearly describe the motion attitude of the ship and obtain a transformation matrix when the roll inclination angle phi is obtainedStep three: and (5) analyzing the stress of the model. Step four: by the formula
The moments of gravity G, buoyancy F and external force W acting on the ship model can be obtained. According to the condition that the resultant moment is 0 in the steady state process, the relation of the change of the simulated deck roll inclination angle phi with high accuracy along with the action of an external force can be obtained.
Description
The technical field is as follows:
the invention relates to an analysis method for ship swinging angle along with external force, and particularly discloses an analysis method for simulating the relation between the change of a deck rolling inclination angle along with external force.
Background art:
under the action of sea waves and sea winds, the ship can generate attitude motions such as rolling, pitching and lifting, and the like, so that the operation of the ship-based aircraft and other equipment and personnel can be seriously influenced. In order to reproduce the motion posture of a ship, the invention provides a 70m multiplied by 20m full-size deck which can float in a still water area, realize controllable swinging motion under the drive of a plurality of groups of power systems carried by the deck, and provide an experimental platform for training or practical installation and verification of carrier-based aircrafts or other equipment and personnel. The invention relates to a method for analyzing the change relation of the rolling angle of a ship bottom fixed simulation deck along with the action of an external force.
The invention content is as follows:
the invention aims to provide an analysis method for simulating the relationship between the change of a deck roll inclination angle along with external force.
The method is realized by the following steps:
the method comprises the following steps: a ship body model with reasonable design and an ideal experimental environment fix the ship bottom in a positive floating state, so that only rolling motion occurs.
When the shape design of a simulation deck is carried out, the following conditions are met:
(1) at the maximum design inclination angle phimaxWhen the deck is not submerged below the water line;
(2) at the maximum design inclination angle phimaxWhen the ship body does not overturn;
(3) the device can automatically return to the initial state without the aid of driving force.
From the above conditions, a hull model was designed as shown in fig. 1. The designed ship model is shown in figure 1, the designed experimental environment is shown in figure 2, and the experimental environment consists of a ship model (2), a driving device (1), a fixing device (3) and a water pool (4). The driving device (1) provides external acting moment for the ship model, and the fixing device (3) is installed at the bottom of the ship model, so that the ship model can only generate rolling motion.
Step two: establishing a geodetic coordinate system OXYZ and a hull coordinate system OX1Y1Z1So as to more clearly describe the motion attitude of the ship and obtain a transformation matrix when the roll inclination angle phi is obtained
To accurately traceAnd respectively establishing a ship body coordinate system and a geodetic coordinate system according to the swing attitude of the simulated deck. The origin of the geodetic coordinate system is established on the intersection points of the bottom plane, the middle transverse section and the middle longitudinal section of the hull model and is marked as a point O, the intersection line of the centerline plane and the base plane of the hull model is taken as an X axis, the stern points to the bow of the ship as a positive direction, the intersection line of the middle station plane and the base plane is taken as a Y axis, the stern points to the starboard as a positive direction, the intersection line of the middle station plane and the centerline plane is taken as a Z axis, and the vertical direction is taken as a positive direction, so that the geodetic coordinate system OXYZ is. Similarly, establishing a hull coordinate system OX1Y1Z1And in the initial state, the geodetic coordinate system is superposed with the ship body coordinate system.
In the formula (I), the compound is shown in the specification,
is the vector of a certain particle in the geodetic coordinate system OXYZ.
The radius is the origin of a ship coordinate system;
the vector diameter of a certain particle in a ship body coordinate system;
and
a unit vector of a geodetic coordinate system;
and
is the unit vector of the ship coordinate system.
Multiplying the expression (1) by a unit vector in order
And
then, obtaining:
the geodetic coordinate system coinciding with the hull coordinate system, i.e. X0＝0，Y0＝0，Z0＝0。
In the above formula, the unit vector products are respectively:
from the above formula, the nine cosine angles for determining the hull coordinate system and the earth coordinate system can be represented by orthogonal conditions except 3 ones. The angle representation of the physical model coordinate system in the geodetic coordinate system can be represented by 3 Euler angle cosine matrices.
Transformation matrix of the roll angle Φ:
step three: and (5) analyzing the stress of the model.
And applying external force to the designed model to enable the model to generate rolling motion, and analyzing the stress when the inclination angle is phi according to the fact that the motion process is a static process to obtain the stress of the hull model in the state, namely gravity G, buoyancy F and external force W. Because the dragging force provided by the fixing device passes through the origin of coordinates, the moment generated by the fixing device on the ship is 0, and the moment is not considered in calculation.
Step four: by the formula
The moments of action of gravity G, buoyancy F and external force W on the hull model can be obtained, and the change relation of the simulated deck roll inclination angle phi along with the external force action can be obtained according to the condition that the resultant moment is 0 in the steady-state process.
In the process of calculating the moment, the direction of gravity is always vertical downward relative to a fixed coordinate system, and the magnitude of gravity is G, and according to the conditions, a moment formula is obtained:
Z1G-gravity at Z1Coordinate position on axis, m.
The components of the moment generated by gravity on each coordinate axis are obtained by the formula:
MxG-the moment of gravity in the geodetic coordinate system, n.m.
Similarly, the position of the external force in the hull coordinate system is expressed in the geodetic coordinate system:
Wx-the component of the external force on the X-axis in the geodetic coordinate system, N;
Wy-the component of the external force on the Y axis in the geodetic coordinate system, N;
Wz-the component of the external force on the Z-axis in the geodetic coordinate system, N;
w- -the representation of the external force in the hull coordinate system, N.
By analogy, the relative position change of the external force action point in the ship coordinate system relative to the geodetic coordinate system is as follows:
XW-the component of the external force in the direction of the X-axis in the geodetic coordinate system, m;
YW-the component of the external force in the Y-axis direction in the geodetic coordinate system, m;
ZW-the component of the external force in the direction of the Z-axis in the geodetic coordinate system, m;
Z1W-the component of the external force in the Z-axis direction in the vessel coordinate system, m.
MWx-the component of the external moment in the X-axis direction in the geodetic coordinate system, n.m;
the invention takes the rolling research of the ship body as an example, namely the ship body model rotates around an X axis in the established geodetic coordinate system, so that the coordinate position of the acting point of the buoyancy on the X axis does not change in the ship body coordinate system and is (0, Y)1F，Z1F). Since the buoyancy magnitude and the center of buoyancy position change with the change in the volume of the hull model submerged in water, a plurality of variable calculations are introduced in the calculation of the buoyancy. Firstly, converting a buoyancy coordinate point from a ship body coordinate system to a geodetic coordinate system:
XF-the component of buoyancy in the direction of the X axis in the geodetic coordinate system, m;
YF-the component of buoyancy in the Y-axis direction in the geodetic coordinate system, m;
ZF-the component of buoyancy in the direction of the Z axis in the geodetic coordinate system, m;
Y1F-the component of buoyancy in the direction of the Y axis in the vessel coordinate system, m;
Z1F- - - -superficialThe component of the force in the Z-axis direction in the vessel coordinate system, m.
MxF-the component of the buoyancy moment in the X-axis direction in the geodetic coordinate system, n.m;
from the above calculation, taking the counterclockwise direction as positive, we can obtain the moment equation of each force to the physical model in the geodetic coordinate system:
Mx-the component of the total moment in the geodetic coordinate system in the direction of the X-axis, n.m;
and (3) considering that any state is a static process, the resultant moment is 0, and the expression of the change relation of the external force along with the inclination angle is obtained as follows:
description of the drawings:
fig. 1 is a model view of a ship hull according to the present invention.
FIG. 2 is a schematic diagram of an experimental environment of the present invention.
Fig. 3 is a schematic diagram of a coordinate system of the present invention.
FIG. 4 is a force analysis diagram of the present invention.
The specific implementation mode is as follows:
the method comprises the following steps: a ship body model with reasonable design and an ideal experimental environment fix the ship bottom in a positive floating state, so that only rolling motion occurs.
When the shape design of a simulation deck is carried out, the following conditions are met:
(1) at the maximum design inclination angle phimaxWhen in use, the deck does not submerge into the water lineThe following;
(2) at the maximum design inclination angle phimaxWhen the ship body does not overturn;
(3) the device can automatically return to the initial state without the aid of driving force.
From the above conditions, a hull model was designed as shown in fig. 1. The designed ship model is shown in figure 1, the designed experimental environment is shown in figure 2, and the experimental environment consists of a ship model (2), a driving device (1), a fixing device (3) and a water pool (4). The driving device (1) provides external acting moment for the ship model, and the fixing device (3) is installed at the bottom of the ship model, so that the ship model can only generate rolling motion.
Step two: establishing a geodetic coordinate system OXYZ and a hull coordinate system OX1Y1Z1So as to more clearly describe the motion attitude of the ship and obtain a transformation matrix when the roll inclination angle phi is obtained
In order to accurately describe the swing attitude of the simulation deck, a ship body coordinate system and a geodetic coordinate system are respectively established. The origin of the geodetic coordinate system is established on the intersection points of the bottom plane, the middle transverse section and the middle longitudinal section of the hull model and is marked as a point O, the intersection line of the centerline plane and the base plane of the hull model is taken as an X axis, the stern points to the bow of the ship as a positive direction, the intersection line of the middle station plane and the base plane is taken as a Y axis, the stern points to the starboard as a positive direction, the intersection line of the middle station plane and the centerline plane is taken as a Z axis, and the vertical direction is taken as a positive direction, so that the geodetic coordinate system OXYZ is. Similarly, establishing a hull coordinate system OX1Y1Z1And in the initial state, the geodetic coordinate system is superposed with the ship body coordinate system.
In the formula (I), the compound is shown in the specification,
is the vector of a certain particle in the geodetic coordinate system OXYZ.
The radius is the origin of a ship coordinate system;
the vector diameter of a certain particle in a ship body coordinate system;
and
a unit vector of a geodetic coordinate system;
and
is the unit vector of the ship coordinate system.
Multiplying the expression (1) by a unit vector in orderAnd
then, obtaining:
the geodetic coordinate system coinciding with the hull coordinate system, i.e. X0＝0，Y0＝0，Z0＝0。
In the above formula, the unit vector products are respectively:
from the above formula, the nine cosine angles for determining the hull coordinate system and the earth coordinate system can be represented by orthogonal conditions except 3 ones. The angle representation of the physical model coordinate system in the geodetic coordinate system can be represented by 3 Euler angle cosine matrices.
Transformation matrix of the roll angle Φ:
step three: and (5) analyzing the stress of the model.
And applying external force to the designed model to enable the model to generate rolling motion, and analyzing the stress when the inclination angle is phi according to the fact that the motion process is a static process to obtain the stress of the hull model in the state, namely gravity G, buoyancy F and external force W. Because the dragging force provided by the fixing device passes through the origin of coordinates, the moment generated by the fixing device on the ship is 0, and the moment is not considered in calculation.
Step four: by the formulaThe moments of action of gravity G, buoyancy F and external force W on the hull model can be obtained, and the change relation of the simulated deck roll inclination angle phi along with the external force action can be obtained according to the condition that the resultant moment is 0 in the steady-state process.
In the process of calculating the moment, the direction of gravity is always vertical downward relative to a fixed coordinate system, and the magnitude of gravity is G, and according to the conditions, a moment formula is obtained:
Z1G-gravity at Z1Coordinate position on axis, m.
The components of the moment generated by gravity on each coordinate axis are obtained by the formula:
MxG-the moment of gravity in the geodetic coordinate system, n.m.
Similarly, the position of the external force in the hull coordinate system is expressed in the geodetic coordinate system:
Wx-the component of the external force on the X-axis in the geodetic coordinate system, N;
Wy-the component of the external force on the Y axis in the geodetic coordinate system, N;
Wz-the component of the external force on the Z-axis in the geodetic coordinate system, N;
w- -the representation of the external force in the hull coordinate system, N.
By analogy, the relative position change of the external force action point in the ship coordinate system relative to the geodetic coordinate system is as follows:
XW-the component of the external force in the direction of the X-axis in the geodetic coordinate system, m;
YW-the component of the external force in the Y-axis direction in the geodetic coordinate system, m;
ZW-the component of the external force in the direction of the Z-axis in the geodetic coordinate system, m;
Z1W-the component of the external force in the Z-axis direction in the vessel coordinate system, m.
MWx-the component of the external moment in the X-axis direction in the geodetic coordinate system, n.m;
the invention takes the rolling research of the ship hull as an example, namely, the ship hull model rotates around an X axis in the established geodetic coordinate system, so that the ship hull model is arranged on the ship hullIn the coordinate system, the coordinate position of the acting point of the buoyancy on the X axis does not change and is (0, Y)1F，Z1F). Since the buoyancy magnitude and the center of buoyancy position change with the change in the volume of the hull model submerged in water, a plurality of variable calculations are introduced in the calculation of the buoyancy. Firstly, converting a buoyancy coordinate point from a ship body coordinate system to a geodetic coordinate system:
XF-the component of buoyancy in the direction of the X axis in the geodetic coordinate system, m;
YF-the component of buoyancy in the Y-axis direction in the geodetic coordinate system, m;
ZF-the component of buoyancy in the direction of the Z axis in the geodetic coordinate system, m;
Y1F-the component of buoyancy in the direction of the Y axis in the vessel coordinate system, m;
Z1F-the component of buoyancy in the Z-axis direction in the vessel coordinate system, m.
MxF-the component of the buoyancy moment in the X-axis direction in the geodetic coordinate system, n.m;
from the above calculation, taking the counterclockwise direction as positive, we can obtain the moment equation of each force to the physical model in the geodetic coordinate system:
Mx-the component of the total moment in the geodetic coordinate system in the direction of the X-axis, n.m;
and (3) considering that any state is a static process, the resultant moment is 0, and the expression of the change relation of the external force along with the inclination angle is obtained as follows:
Claims (1)
1. a method for analyzing the relation between the roll inclination angle of a simulation deck and the change of an external force comprises the following steps:
the method comprises the following steps: a ship body model with reasonable design and an ideal experimental environment fix the ship bottom in a positive floating state, so that only rolling motion occurs.
Step two: establishing a geodetic coordinate system OXYZ and a hull coordinate system OX1Y1Z1So as to more clearly describe the motion attitude of the ship and obtain a transformation matrix when the roll inclination angle phi is obtained
Step three: and (5) analyzing the stress of the model.
Step four: by the formulaThe moments of action of gravity G, buoyancy F and external force W on the hull model can be obtained, and the change relation of the simulated deck roll inclination angle phi along with the external force action can be obtained according to the condition that the resultant moment is 0 in the steady-state process.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911115932.XA CN110834703A (en) | 2019-11-15 | 2019-11-15 | Analysis method for simulating relationship between change of deck roll inclination angle along with external force |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911115932.XA CN110834703A (en) | 2019-11-15 | 2019-11-15 | Analysis method for simulating relationship between change of deck roll inclination angle along with external force |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110834703A true CN110834703A (en) | 2020-02-25 |
Family
ID=69575078
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911115932.XA Pending CN110834703A (en) | 2019-11-15 | 2019-11-15 | Analysis method for simulating relationship between change of deck roll inclination angle along with external force |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110834703A (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20050097014A (en) * | 2004-03-30 | 2005-10-07 | 삼성탈레스 주식회사 | Method for correcting sight error of aiming apparatus established in ship |
| CN101881620A (en) * | 2010-06-25 | 2010-11-10 | 哈尔滨工程大学 | Ship swaying and surging information measuring method based on optical fiber gyro inertial measuring system |
| CN103727937A (en) * | 2013-11-20 | 2014-04-16 | 中国人民解放军海军大连舰艇学院 | Star sensor based naval ship attitude determination method |
| JP2014165511A (en) * | 2013-02-21 | 2014-09-08 | Nippon Telegr & Teleph Corp <Ntt> | Satellite tracking antenna system, and satellite tracking antenna control method |
| CN104459608A (en) * | 2014-11-27 | 2015-03-25 | 中国船舶重工集团公司第七二四研究所 | Ship pitch and roll comprehensive vector coordinate transformation antenna beam pointing direction correcting method |
| CN105825061A (en) * | 2016-03-17 | 2016-08-03 | 大连海事大学 | Method for calculating random floating state of ship on basis of STL model |
-
2019
- 2019-11-15 CN CN201911115932.XA patent/CN110834703A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20050097014A (en) * | 2004-03-30 | 2005-10-07 | 삼성탈레스 주식회사 | Method for correcting sight error of aiming apparatus established in ship |
| CN101881620A (en) * | 2010-06-25 | 2010-11-10 | 哈尔滨工程大学 | Ship swaying and surging information measuring method based on optical fiber gyro inertial measuring system |
| JP2014165511A (en) * | 2013-02-21 | 2014-09-08 | Nippon Telegr & Teleph Corp <Ntt> | Satellite tracking antenna system, and satellite tracking antenna control method |
| CN103727937A (en) * | 2013-11-20 | 2014-04-16 | 中国人民解放军海军大连舰艇学院 | Star sensor based naval ship attitude determination method |
| CN104459608A (en) * | 2014-11-27 | 2015-03-25 | 中国船舶重工集团公司第七二四研究所 | Ship pitch and roll comprehensive vector coordinate transformation antenna beam pointing direction correcting method |
| CN105825061A (en) * | 2016-03-17 | 2016-08-03 | 大连海事大学 | Method for calculating random floating state of ship on basis of STL model |
Non-Patent Citations (3)
| Title |
|---|
| 刘胜 等: "《自动控制原理》", 30 September 2015, 哈尔滨工程大学出版社 * |
| 苏振超 等: "《理论力学》", 31 August 2018, 华中科技大学出版社 * |
| 蒋维清 等: "《船舶原理》", 30 November 1979, 人民交通出版社 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103942383B (en) | The dynamics and kinematics method of estimation of a kind of deep ocean work type underwater robot | |
| Xin et al. | A vectored water jet propulsion method for autonomous underwater vehicles | |
| Yue et al. | Mechatronic system and experiments of a spherical underwater robot: SUR-II | |
| Hai et al. | Investigation on the mechanical design and manipulation hydrodynamics for a small sized, single body and streamlined I-AUV | |
| Wang et al. | Dynamics modeling of an unmanned wave glider with flexible umbilical | |
| Shen et al. | Numerical simulations of ship self-propulsion and maneuvering using dynamic overset grids in OpenFOAM | |
| Jiao et al. | Experimental investigation of wave-induced ship hydroelastic vibrations by large-scale model measurement in coastal waves | |
| CN107037821A (en) | The estimation of underwater hiding-machine athletic posture and control method under repetitive shocks | |
| CN110321639B (en) | Quasi-physical numerical method for realizing underwater robot four-degree-of-freedom vertical plane self-navigation diving prediction | |
| Su et al. | Dynamics modeling and simulation of autonomous underwater vehicles with appendages | |
| CN111709151A (en) | Method for creating ship rapid mathematical model | |
| Akimoto et al. | Finite-volume simulation method to predict the performance of a sailing boat | |
| Qiu et al. | 3D motion model for the freefall lifeboat during its launching from a moving ship | |
| CN110834703A (en) | Analysis method for simulating relationship between change of deck roll inclination angle along with external force | |
| CN108460206B (en) | Wave glider motion prediction method | |
| Jain et al. | Design and Development of an Open-frame AUV: ANAHITA | |
| CN103400517A (en) | Semi-physical simulation system and manipulating method during lifesaving process of deep submergence rescue vessel | |
| Richardt et al. | Maneuvering simulations for ships and sailing yachts using FRIENDSHIP-Equilibrium as an open modular workbench | |
| Chen et al. | Experimental study on coupled motions of mother ship launching and recovering of human-occupied vehicle in regular waves | |
| CN110688733B (en) | Numerical simulation method for water surface separation of carrier type submarine-launched missile | |
| Guangyi et al. | Research on underwater safety inspection and operational robot motion control | |
| CN112307419A (en) | Rapid forecasting method for vertical motion response of ship | |
| CN107463754B (en) | Offshore boundary condition simulation method for FLNG device liquefaction process | |
| Gauvain | The un-restrained sailing yacht model tests–a new approach and technology appropriate to modern sailing yacht seakeeping | |
| Koskinen | Numerical simulation of ship motion due to waves and manoeuvring |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200225 |
|
| WD01 | Invention patent application deemed withdrawn after publication |