CN114152245B - Multi-dimensional motion pose measurement system and calculation method for underwater suspension tunnel test - Google Patents

Multi-dimensional motion pose measurement system and calculation method for underwater suspension tunnel test Download PDF

Info

Publication number
CN114152245B
CN114152245B CN202210125047.5A CN202210125047A CN114152245B CN 114152245 B CN114152245 B CN 114152245B CN 202210125047 A CN202210125047 A CN 202210125047A CN 114152245 B CN114152245 B CN 114152245B
Authority
CN
China
Prior art keywords
tunnel
laser sensor
displacement
displacement value
underwater
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210125047.5A
Other languages
Chinese (zh)
Other versions
CN114152245A (en
Inventor
阳志文
崔永刚
季新然
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin Research Institute for Water Transport Engineering MOT
Original Assignee
Tianjin Research Institute for Water Transport Engineering MOT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin Research Institute for Water Transport Engineering MOT filed Critical Tianjin Research Institute for Water Transport Engineering MOT
Priority to CN202210125047.5A priority Critical patent/CN114152245B/en
Publication of CN114152245A publication Critical patent/CN114152245A/en
Application granted granted Critical
Publication of CN114152245B publication Critical patent/CN114152245B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C15/00Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
    • G01C15/002Active optical surveying means
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention relates to the technical field of ocean engineering, in particular to a multidimensional movement pose measuring system and a calculating method for an underwater suspension tunnel test, and solves the problem of measuring multidimensional and multimodal flexible movement deformation of an underwater suspension tunnel. The measurement and calculation method and the system can measure the bending, torsion or combined dynamic deformation of the suspension tunnel, obtain the displacement or posture of any section of the tunnel pipe body and realize the multi-dimensional and multi-mode flexible motion response measurement of the full-span structure of the suspension tunnel; the measuring method and the system provided by the invention have the advantages of clear principle, simple system composition, simple and convenient algorithm, high operation efficiency and convenience for installation and use, can be used for the research of the suspension tunnel technology, can also be used for the test measurement of other underwater long-span suspension and floating structures, and have wide application prospect.

Description

Multi-dimensional motion pose measurement system and calculation method for underwater suspension tunnel test
Technical Field
The invention relates to the technical field of ocean engineering, in particular to a multi-dimensional motion pose measuring system and a calculating method for an underwater suspension tunnel test.
Background
The floating tunnel is a large-scale cross-sea tunnel structure constructed in water and is a strategic and subversive technology for solving the crossing of deep water fjords in the future, but at present, the floating tunnel is not constructed in the world and is still in the research stage. Due to the large structural span and complex composition, the hydrodynamic response mechanism of the suspended tunnel structure system is very complex under the transient and variable marine power environment, and the problem to be broken through by researchers in various countries at present is urgently solved. The method directly simulates the effect of complex wave water flow load on the suspension tunnel in a test pool, obtains data such as the motion displacement and the attitude of the structure, and is the most direct means for researching hydrodynamic response of the suspension tunnel. The motion form of the long-span suspension tunnel structure is very complex, combined flexible motion deformation such as bending and torsion can be presented on the same tunnel section, more complex multi-dimensional and multi-modal flexible motion deformation can be presented on the whole span, and the measurement difficulty in a model test is greatly increased.
The existing non-contact floating body structure measuring system can only obtain six-dimensional motion data of a single target point, cannot measure multi-dimensional flexible motion of a full-span structure of a suspension tunnel, and cannot obtain displacement or posture of any section of a tunnel pipe body. Meanwhile, because the tunnel model is submerged underwater, the traditional contact-type measuring instrument can generate large interference on a wave water flow dynamic field, and the precision of test measurement is further influenced. In addition, the existing water measurement method needs to erect a higher and heavier instrument mounting bracket on the test model, and the additional mass and the vibration of the bracket caused by the instrument mounting bracket also have great influence on the measurement result.
Disclosure of Invention
Aiming at the problems, the invention aims to overcome the defects of the prior art, provides a multi-dimensional motion pose measurement system and a calculation method for an underwater suspension tunnel test, and solves the problem of multi-dimensional and multi-modal flexible motion deformation measurement of the underwater suspension tunnel. The invention adopts the underwater measurement method, thereby effectively avoiding the influence of the additional mass caused by the instrument mounting bracket and the self vibration of the bracket on the measurement result in the existing water measurement method; compared with the traditional single-target point six-dimensional floating body measurement system, the measurement method and the system provided by the invention can measure the bending, torsion or combined dynamic deformation of the suspended tunnel, obtain the displacement or posture of any section of the tunnel pipe body and realize the multi-dimensional and multi-modal flexible motion response measurement of the full-span structure of the suspended tunnel; compared with the traditional contact type measuring method, the measuring method and the system provided by the invention do not require the contact of the measuring instrument and the measured object, and the interference of the instrument on the wave water flow dynamic field is reduced to the maximum extent. Meanwhile, the measuring method and the system provided by the invention have the advantages of clear principle, simple system composition, simple algorithm, high operation efficiency and convenience in installation and use, can be used for research of the suspension tunnel technology, can also be used for test measurement of other underwater long-span suspension and floating structures, and have wide application prospects.
The technical scheme adopted by the invention for solving the technical problems is as follows: experimental multidimension degree motion position appearance measurement system of suspension tunnel under water, including suspension tunnel body, it is submarine, suspension tunnel surface is equipped with multiunit reflex reflector, reflex reflector is fixed connection with the suspension tunnel, the reflex reflector lower extreme is equipped with the displacement sensor device, displacement sensor device and submarine fixed connection, the displacement sensor device is connected with pencil C electricity, be equipped with turbidity appearance on the pencil C, turbidity appearance and pencil C electricity are connected, pencil C one end is equipped with multichannel data acquisition transceiver, multichannel data acquisition transceiver one end is equipped with the host computer, the host computer is connected with multichannel data acquisition transceiver electricity, the multichannel data acquisition transceiver other end is connected with turbidity appearance and pencil C electricity.
Further, the displacement sensor device comprises a laser sensor A, a laser sensor B, a laser sensor C and a wiring harness A, wherein the laser sensor A, the laser sensor B and the laser sensor C are electrically connected with each other, the laser sensor A, the laser sensor B and the laser sensor C are fixedly connected with the water bottom, and the laser sensor A, the laser sensor B and the laser sensor C are electrically connected with the wiring harness C through the wiring harness A.
The calculation method of the multi-dimensional motion pose measurement system for the underwater suspension tunnel test comprises the following steps:
s1: firstly, establishing a full-span coordinate system of an underwater suspension tunnel model;
s2: the coordinates y of each set of displacement sensor devices are then determinediX-direction distance DiAnd the radius R of the suspension tunnel;
s3: then obtaining the displacement value measured by each group of displacement sensor devices at the initial moment and the displacement value of the laser sensor A
Figure 351379DEST_PATH_IMAGE001
Displacement value of laser sensor B
Figure 176115DEST_PATH_IMAGE002
Displacement value of laser sensor C
Figure 905168DEST_PATH_IMAGE003
S4: then obtaining the displacement value measured by each group of displacement sensor devices at any other time and the displacement value of the laser sensor A
Figure 891579DEST_PATH_IMAGE004
Displacement value of laser sensor B
Figure 588139DEST_PATH_IMAGE005
Displacement value of laser sensor C
Figure 203185DEST_PATH_IMAGE006
S5: then calculating tube body at section position of suspension tunnel measured by each group of displacement sensor devices at any momentxTo the relative displacement value
Figure 36011DEST_PATH_IMAGE007
B, andzrelative displacement value of
Figure 193323DEST_PATH_IMAGE008
S6: and finally, calculating to obtain the displacement and attitude value of any section of the underwater suspended tunnel.
Further, the origin of coordinates of the full-span coordinate system of the suspended tunnel model in S1 is located at the center of one end of the tunnel, the y-axis of the coordinate system is located in the horizontal plane and points to the axial direction of the tunnel, the x-axis of the coordinate system is located in the horizontal plane and points to the downstream direction of the wave water flow perpendicular to the axial direction of the tunnel, and the z-axis of the coordinate system is vertically upward.
Further, y in S2iIs a projection point of the coordinate origin on the water bottom to the secondiDistance of the lines of the group displacement sensor devices, DiAnd R is the radius of the suspended tunnel pipe body and is the minimum value of the distance between every two adjacent laser sensors.
Further, among the S3 and S4
Figure 862333DEST_PATH_IMAGE009
And
Figure 28872DEST_PATH_IMAGE010
defined as the displacement value measured by the laser sensor a,
Figure 716206DEST_PATH_IMAGE003
and
Figure 293686DEST_PATH_IMAGE011
defined as the displacement value measured by the laser sensor C,
Figure 699260DEST_PATH_IMAGE002
and
Figure 669490DEST_PATH_IMAGE012
the displacement value measured by the laser sensor B is defined, the 0 time is indicated by the superscript "0", and the other arbitrary time is indicated by the superscript "n".
Further, in said S5
Figure 962062DEST_PATH_IMAGE007
,
Figure 726756DEST_PATH_IMAGE008
The calculation method of (2) is as follows: the vertical upward and horizontal downstream movement of the tube body is defined as positive.
When in use
Figure 619625DEST_PATH_IMAGE013
Figure 380164DEST_PATH_IMAGE014
Then, the following equation is satisfied:
Figure 776511DEST_PATH_IMAGE016
when in use
Figure 712106DEST_PATH_IMAGE017
Figure 843004DEST_PATH_IMAGE018
Then, the following equation is satisfied:
Figure 889457DEST_PATH_IMAGE020
in the formula (I), the compound is shown in the specification,D i R
Figure 920736DEST_PATH_IMAGE021
Figure 27232DEST_PATH_IMAGE022
Figure 894694DEST_PATH_IMAGE023
Figure 495571DEST_PATH_IMAGE024
Figure 866509DEST_PATH_IMAGE025
Figure 143907DEST_PATH_IMAGE026
all the measured quantities are known quantities, and the displacement value of the pipe body at each measured cross section at any n moments is calculated by simultaneous equation sets and utilizing the Newton iteration method and the least square method principle
Figure 497002DEST_PATH_IMAGE027
,
Figure 150837DEST_PATH_IMAGE028
Further, the calculation method in S6 is that, after obtaining the displacements at several known cross sections, the displacement and attitude data at any cross section along the tunnel at any time can be obtained by interpolation,
Figure 376282DEST_PATH_IMAGE029
Figure 575313DEST_PATH_IMAGE030
in the formula (I), the compound is shown in the specification,
Figure 417367DEST_PATH_IMAGE031
and
Figure 609314DEST_PATH_IMAGE032
respectively the x-direction and z-direction relative displacement values at the y-coordinate position along the tunnel at any moment,α k andβ k respectively at any time along the tunnely k The rotation angle of the tube body around the z-axis and the x-axis is at the coordinate.
The invention has the advantages that: the invention provides a multi-dimensional motion pose measuring system and a calculating method for an underwater suspended tunnel test, and the underwater measurement method is adopted, so that the influence of additional mass caused by an instrument mounting bracket and the vibration of the bracket on a measuring result in the existing water measurement method is effectively avoided; compared with the traditional single-target point six-dimensional floating body measurement system, the measurement method and the system provided by the invention can measure the bending, torsion or combined dynamic deformation of the suspended tunnel, obtain the displacement or posture of any section of the tunnel pipe body and realize the multi-dimensional and multi-modal flexible motion response measurement of the full-span structure of the suspended tunnel; compared with the traditional contact type measuring method, the measuring method and the system provided by the invention do not require the contact of the measuring instrument and the measured object, and the interference of the instrument on the wave water flow dynamic field is reduced to the maximum extent. Meanwhile, the measuring method and the system provided by the invention have the advantages of clear principle, simple system composition, simple algorithm, high operation efficiency and convenience in installation and use, can be used for research of the suspension tunnel technology, can also be used for test measurement of other underwater long-span suspension and floating structures, and have wide application prospects.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a schematic perspective view of the present invention;
FIG. 2 is a schematic plan view, partially in section, of the present invention;
FIG. 3 is a schematic diagram of the coordinates of a laser sensor according to the present invention;
FIG. 4 is a schematic plan view of a partial cross-sectional coordinate system of the present invention;
FIG. 5 is a schematic diagram of the displacement coordinates of the suspension conduit of the present invention;
FIG. 6 shows the measurement using the laser sensor A of the present invention
Figure 938533DEST_PATH_IMAGE027
,
Figure 823313DEST_PATH_IMAGE028
Reference is made to the figure;
FIG. 7 shows the measurement using the laser sensor B of the present invention
Figure 887084DEST_PATH_IMAGE027
,
Figure 367875DEST_PATH_IMAGE028
Reference is made to the figure;
FIG. 8 shows the measurement using the laser sensor C of the present invention
Figure 567912DEST_PATH_IMAGE027
,
Figure 623592DEST_PATH_IMAGE028
Reference is made to the figure;
FIG. 9 is a schematic diagram of the present invention for full span displacement interpolation using displacement at a known cross-section.
Wherein:
1. a suspension tunnel; 2. a light reflecting means; 3. a displacement sensor device; 31. a laser sensor A; 32. a laser sensor B; 33. a laser sensor C; 34. a wire harness A; 4. a wire harness C; 5. a turbidity meter; 6. an upper computer; 7. a multi-channel data acquisition transceiver; 8. the surface of the water.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the accompanying drawings, and it should be understood that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the description of the present invention, it should be noted that unless otherwise explicitly specified or limited, the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", and the like indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of description and simplicity of description, and do not indicate or imply that the referred device or element must have a specific orientation, be constructed in a specific orientation, and be operated, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance, and the terms "mounted," "connected," and "connected" are to be construed broadly and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
Example 1:
fig. 1 is a schematic view of a three-dimensional structure of the invention, fig. 2 is a schematic view of a plane partial section of the invention, a plurality of groups of light reflecting devices 2 are arranged on the outer surface of a suspension tunnel 1, the light reflecting devices 2 are fixedly connected with the suspension tunnel 1, a displacement sensor device 3 is arranged at the lower end of the light reflecting devices 2, the displacement sensor device 3 is fixedly connected with the water bottom, the displacement sensor device 3 is electrically connected with a wiring harness C4, a turbidity meter 5 is arranged on the wiring harness C4, the turbidity meter 5 is electrically connected with the wiring harness C4, a multi-channel data acquisition transceiver 7 is arranged at one end of the wiring harness C4, an upper computer 6 is arranged at one end of the multi-channel data acquisition transceiver 7, the upper computer 6 is electrically connected with the multi-channel data acquisition transceiver 7, and the other end of the multi-channel data acquisition transceiver 7 is electrically connected with the turbidity meter 5 and the wiring harness C4. The displacement sensor device 3 comprises a laser sensor A31, a laser sensor B32, a laser sensor C33 and a wiring harness A34, wherein the laser sensor A31, the laser sensor B32 and the laser sensor C33 are all sensors of L4-40 types, the laser sensor A31, the laser sensor B32 and the laser sensor C33 are electrically connected with each other, the laser sensor A31, the laser sensor B32 and the laser sensor C33 are fixedly connected with the water bottom, and the laser sensor A31, the laser sensor B32 and the laser sensor C33 are electrically connected with the wiring harness C4 through the wiring harness A34. The underwater displacement sensor device comprises a plurality of groups of underwater displacement sensor devices 3 arranged along the axis direction of a tunnel, at least 3 sensors in each group are arranged under a tube body side by side at equal intervals, the connecting line of the sensors is vertical to the axis direction of the suspended tunnel, one of the sensors is arranged on the axis of the tube body passing through the center of a circle, and the other two sensors are arranged on two sides of the tube body at equal intervals. The reflecting device 2 can be designed into a flexible white material with a rough surface, is attached to the lower surface of the suspension tunnel 1 test model, can form diffuse reflection on laser, and the area of the reflecting device 2 is enough to cover the irradiation range of a laser spot. Laser emitted by the laser sensor is reflected back by the reflecting device below the test model of the suspension tunnel 1 and then returns back to the laser sensor again, so that distance data between the laser sensor and the corresponding reflecting point is obtained. The turbidity meter 5 is arranged near any one group of displacement sensor devices 3 under water and is used for measuring the turbidity of water in real time so as to correct and compensate the distance measurement value of the underwater displacement sensor and ensure the measurement accuracy. The multichannel data acquisition transceiver 7 is connected with all the underwater displacement sensor devices 3 and the turbidity meter 5, can synchronously acquire distance data and underwater turbidity data of the underwater laser sensor, and converts the distance data and the underwater turbidity data into digital signals through an internal data conversion module to be transmitted to an operating system of the upper computer 6. The upper computer 6 operating system further integrates and analyzes the collected digital signals through a matched software analysis system, so that the multi-dimensional and multi-modal motion displacement and attitude data of the whole underwater suspension tunnel 1 are calculated.
Example 2:
the calculation method is as follows
FIG. 3 is a schematic diagram of the coordinates of a laser sensor according to the present invention, and FIG. 4 is a schematic diagram of the coordinates of a planar partial cross-section according to the present invention; fig. 5 is a schematic diagram of displacement coordinates of a suspension pipeline, and the calculation method of the underwater suspension tunnel test multi-dimensional motion pose measurement system shown in fig. 3 includes the following steps:
s1: firstly, establishing a full-span coordinate system of an underwater suspension tunnel 1 model;
s2: the coordinates y of each set of displacement sensor devices 3 are then determinediX-direction distance DiAnd the radius R of the suspension tunnel 1;
s3: then, the displacement value measured by each set of displacement sensor device 3 at the initial time, the displacement value of the laser sensor a31 is obtained
Figure 426857DEST_PATH_IMAGE033
Displacement value of laser sensor B32
Figure 695027DEST_PATH_IMAGE022
Displacement value of the laser sensor C33
Figure 15150DEST_PATH_IMAGE023
S4: then obtaining the displacement value measured by each group of displacement sensor devices at any other time and the displacement value of the laser sensor A
Figure 726885DEST_PATH_IMAGE034
Displacement value of laser sensor B
Figure 765248DEST_PATH_IMAGE025
Displacement value of laser sensor C
Figure 837110DEST_PATH_IMAGE026
S5: then calculating the relative displacement value of the pipe body at the section of the suspension tunnel 1 measured by each group of displacement sensor devices 3 at any moment
Figure 526586DEST_PATH_IMAGE035
,
Figure 658490DEST_PATH_IMAGE028
S6: and finally, calculating to obtain the displacement and attitude value of any section of the underwater suspended tunnel 1.
The origin of coordinates of the full-span coordinate system of the S1 suspension tunnel 1 model is arranged at the center of one end of the tunnel, the y axis of the coordinate system is arranged in the horizontal plane and points to the axial direction of the tunnel, the x axis of the coordinate system is arranged in the horizontal plane and is perpendicular to the axial direction of the tunnel and points to the downstream direction of the wave water flow, and the z axis of the coordinate system is vertically upward.
Y in S2iIs the distance from the projected point of the origin of coordinates on the water bottom to the connecting line of the i-th group of displacement sensor devices 3, DiFor the minimum of every group adjacent two laser sensor distance between, because for tunnel length, the bending deformation volume in tunnel is very little, so this application is looked displacement sensor device 3 and is measured the cross-section and be the circle cross-section, and R is the radius of suspension tunnel 1 body.
In S3 and S4
Figure 934882DEST_PATH_IMAGE021
And
Figure 810434DEST_PATH_IMAGE036
defined as the displacement value measured by the laser sensor a31,
Figure 839570DEST_PATH_IMAGE023
and
Figure 394572DEST_PATH_IMAGE037
defined as the displacement value measured by the laser sensor C33,
Figure 407528DEST_PATH_IMAGE022
and
Figure 86771DEST_PATH_IMAGE038
the displacement value measured by the laser sensor B32 is defined as the time 0 is indicated by the superscript "0", and the other arbitrary time is indicated by the superscript "n".
In S5
Figure 986725DEST_PATH_IMAGE035
,
Figure 460431DEST_PATH_IMAGE028
The calculation method of (2) is as follows: the vertical upward movement and the horizontal downstream movement of the pipe body are defined as positive respectively,
as shown in fig. 6, when
Figure 695103DEST_PATH_IMAGE039
Figure 161726DEST_PATH_IMAGE040
When the temperature of the water is higher than the set temperature,
Figure 165454DEST_PATH_IMAGE041
Figure 810062DEST_PATH_IMAGE042
Figure 282763DEST_PATH_IMAGE043
Figure 303808DEST_PATH_IMAGE044
Figure 162043DEST_PATH_IMAGE045
Figure 495329DEST_PATH_IMAGE046
Figure 704593DEST_PATH_IMAGE047
because:
Figure 263750DEST_PATH_IMAGE048
therefore:
Figure 992803DEST_PATH_IMAGE049
will be
Figure 244793DEST_PATH_IMAGE050
Figure 675774DEST_PATH_IMAGE051
Figure 287890DEST_PATH_IMAGE052
Figure 855138DEST_PATH_IMAGE053
Bringing in, can obtain:
Figure 294340DEST_PATH_IMAGE054
it is possible to obtain,
Figure 212618DEST_PATH_IMAGE055
as shown in fig. 7, when
Figure 113578DEST_PATH_IMAGE056
Figure 800911DEST_PATH_IMAGE057
When the temperature of the water is higher than the set temperature,
Figure 635182DEST_PATH_IMAGE041
Figure 40755DEST_PATH_IMAGE042
Figure 745406DEST_PATH_IMAGE058
Figure 303558DEST_PATH_IMAGE059
Figure 68251DEST_PATH_IMAGE060
Figure 695542DEST_PATH_IMAGE061
because of the fact that
Figure 187572DEST_PATH_IMAGE048
Therefore, it is not only easy to use
Figure 115076DEST_PATH_IMAGE062
Will be
Figure 801404DEST_PATH_IMAGE050
Figure 915990DEST_PATH_IMAGE051
Figure 962444DEST_PATH_IMAGE052
Figure 262232DEST_PATH_IMAGE053
Bringing in, can obtain:
Figure 368728DEST_PATH_IMAGE063
it is possible to obtain,
Figure 236190DEST_PATH_IMAGE064
as shown in fig. 8, when
Figure 571487DEST_PATH_IMAGE039
Figure 208005DEST_PATH_IMAGE040
When the temperature of the water is higher than the set temperature,
Figure 485402DEST_PATH_IMAGE041
Figure 558270DEST_PATH_IMAGE065
Figure 212105DEST_PATH_IMAGE066
Figure 703129DEST_PATH_IMAGE067
Figure 636581DEST_PATH_IMAGE068
Figure 478635DEST_PATH_IMAGE069
Figure 670582DEST_PATH_IMAGE070
Figure 268310DEST_PATH_IMAGE071
=
Figure 887510DEST_PATH_IMAGE072
Figure 951281DEST_PATH_IMAGE073
Figure 432072DEST_PATH_IMAGE074
Figure 632109DEST_PATH_IMAGE075
(ii) a Will be provided with
Figure 687790DEST_PATH_IMAGE076
Figure 222545DEST_PATH_IMAGE077
Figure 21874DEST_PATH_IMAGE071
And substituting the obtained result to obtain the final product,
Figure 92730DEST_PATH_IMAGE078
it is possible to obtain,
Figure 53732DEST_PATH_IMAGE079
in summary, when
Figure 92095DEST_PATH_IMAGE017
Figure 681733DEST_PATH_IMAGE040
Then, the following equation is satisfied:
Figure 590783DEST_PATH_IMAGE016
in the same way, when
Figure 988267DEST_PATH_IMAGE017
Figure 999079DEST_PATH_IMAGE040
Then, the following equation is satisfied:
Figure 140210DEST_PATH_IMAGE020
in the formula (I), the compound is shown in the specification,D i R
Figure 169346DEST_PATH_IMAGE021
Figure 721419DEST_PATH_IMAGE022
Figure 468795DEST_PATH_IMAGE023
Figure 413618DEST_PATH_IMAGE024
Figure 47992DEST_PATH_IMAGE025
Figure 521699DEST_PATH_IMAGE026
all the measured quantities are known quantities, and the displacement value of the pipe body at each measured cross section at any n moments is calculated by simultaneous equation sets and utilizing the Newton iteration method and the least square method principle
Figure 756371DEST_PATH_IMAGE027
,
Figure 237642DEST_PATH_IMAGE028
The calculation method in S6 is as follows, and after obtaining displacements at a plurality of known cross sections, displacement and attitude data at any cross section along the tunnel at any time can be obtained by using a polynomial interpolation method. In this embodiment, the data interpolation is performed by using the conventional commonly used lagrangian interpolation formula (Lagrange interpolation format).
As shown in FIG. 9, the number of the known cross sections of the suspension tunnel 1 is set asMThen is aty0xAt any n time in the planeMRelative displacement coordinates of the cross section of each suspension tunnel: (y 1
Figure 241370DEST_PATH_IMAGE080
)、(y 2
Figure 151557DEST_PATH_IMAGE081
)…(y i
Figure 624258DEST_PATH_IMAGE027
)…(y M
Figure 645304DEST_PATH_IMAGE082
) (ii) a In the same way, it isy0zAnd M relative displacement coordinates of the cross sections of the suspension tunnels exist at any n moments in the plane: (y 1
Figure 503538DEST_PATH_IMAGE083
)、(y 2
Figure 833894DEST_PATH_IMAGE084
)…(y i
Figure 43159DEST_PATH_IMAGE028
)…(y M
Figure 618628DEST_PATH_IMAGE085
)。
Directly applying Lagrange interpolation formula (Lagrange interpolation formula) to obtain the line of the suspended tunnel at any n momentyAt coordinatesxTo andzthe relative displacement value is as follows:
Figure 331369DEST_PATH_IMAGE029
then, by utilizing the relation between the slope of the curve and the derivative of the function and the principle of solving the angle by the inverse trigonometric function, the line of the suspended tunnel at any moment can be solvedy k Rotation angle of tube body around z axis at coordinateα k And angle of rotation about the x-axisβ k
Figure 583359DEST_PATH_IMAGE086
The displacement and attitude data of any section along the tunnel at any time are brought intoα k Andβ k the formula of the inverse trigonometric function is obtained:
Figure 532116DEST_PATH_IMAGE030
in the formula (I), the compound is shown in the specification,
Figure 894965DEST_PATH_IMAGE031
and
Figure 462212DEST_PATH_IMAGE032
respectively the x-direction and z-direction relative displacement values at the y-coordinate position along the tunnel at any moment,α k andβ k respectively at any time along the tunnely k The rotation angle of the tube body around the z-axis and the x-axis is at the coordinate.
The measuring method and the system provided by the invention have the advantages of clear principle, simple system composition, simple and convenient algorithm, high operation efficiency and convenience for installation and use, can be used for the research of the suspension tunnel technology, can also be used for the test measurement of other underwater long-span suspension and floating structures, and have wide application prospect.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.

Claims (7)

1. Experimental multidimension degree motion position appearance measurement system of suspension tunnel under water, including suspension tunnel (1) body, its characterized in that: the outer surface of the suspension tunnel (1) is provided with a plurality of groups of light reflecting devices (2), the light reflecting devices (2) are fixedly connected with the suspension tunnel (1), the lower end of each light reflecting device (2) is provided with a displacement sensor device (3), the displacement sensor device (3) is fixedly connected with the bottom of the water, the displacement sensor device (3) is electrically connected with a wiring harness C (4), a turbidity meter (5) is arranged on the wiring harness C (4), the turbidity meter (5) is electrically connected with the wiring harness C (4), one end of the wiring harness C (4) is provided with a multi-channel data acquisition transceiver (7), one end of the multi-channel data acquisition transceiver (7) is provided with an upper computer (6), the upper computer (6) is electrically connected with the multi-channel data acquisition transceiver (7), and the other end of the multi-channel data acquisition transceiver (7) is electrically connected with the turbidity meter (5) and the wiring harness C (4); the calculation method of the multi-dimensional motion pose measurement system for the underwater suspension tunnel test comprises the following steps:
s1: firstly, establishing a full-span coordinate system of an underwater suspension tunnel (1) model;
s2: then determining the coordinates y of each set of displacement sensor devices (3)iX-direction distance DiAnd the radius R of the suspension tunnel (1);
s3: then, the displacement value measured by each group of displacement sensor devices (3) at the initial moment and the displacement value of the laser sensor A (31) are obtained
Figure 679318DEST_PATH_IMAGE001
And the displacement value of the laser sensor B (32)
Figure 651297DEST_PATH_IMAGE002
And the displacement value of the laser sensor C (33)
Figure 742881DEST_PATH_IMAGE003
S4: then, the displacement value measured by each group of displacement sensor devices at any other time and the displacement value of the laser sensor A (31) are obtained
Figure 156545DEST_PATH_IMAGE004
And the displacement value of the laser sensor B (32)
Figure 915554DEST_PATH_IMAGE005
And the displacement value of the laser sensor C (33)
Figure 61364DEST_PATH_IMAGE006
S5: then theCalculating the x-direction relative displacement value of the pipe body at the section of the suspension tunnel (1) measured by each group of displacement sensor devices (3) at any moment
Figure 30457DEST_PATH_IMAGE007
And z-direction relative displacement value
Figure 592020DEST_PATH_IMAGE008
S6: and finally, calculating to obtain the displacement and attitude values of any section of the underwater suspended tunnel (1).
2. The underwater suspended tunnel test multi-dimensional motion pose measurement system according to claim 1, characterized in that: the displacement sensor device (3) comprises a laser sensor A (31), a laser sensor B (32), a laser sensor C (33) and a wiring harness A (34), the laser sensor A (31), the laser sensor B (32) and the laser sensor C (33) are electrically connected with each other, the laser sensor A (31), the laser sensor B (32) and the laser sensor C (33) are fixedly connected with the water bottom, and the laser sensor A (31), the laser sensor B (32) and the laser sensor C (33) are electrically connected with the wiring harness C (4) through the wiring harness A (34).
3. The calculation method of the underwater suspended tunnel test multi-dimensional motion pose measurement system according to claim 1, characterized in that: the origin of coordinates of a full-span coordinate system of the suspended tunnel (1) model in the S1 is arranged at the center of one end of the tunnel, the y axis of the coordinate system is arranged in the horizontal plane and points to the axial direction of the tunnel, the x axis of the coordinate system is arranged in the horizontal plane and points to the downstream direction of the wave water flow perpendicular to the axial direction of the tunnel, and the z axis of the coordinate system is vertically upward.
4. The calculation method of the underwater suspended tunnel test multi-dimensional motion pose measurement system according to claim 1, characterized in that: y in S2iIs the distance from the projection point of the coordinate origin point on the water bottom to the connecting line of the ith group of displacement sensor devices (3),Diand R is the radius of the pipe body of the suspension tunnel (1) and is the minimum value of the distance between every two adjacent laser sensors.
5. The calculation method of the underwater suspended tunnel test multi-dimensional motion pose measurement system according to claim 1, characterized in that: of the S3 and S4
Figure 595748DEST_PATH_IMAGE009
And
Figure 912460DEST_PATH_IMAGE010
defined as the displacement value measured by the laser sensor A (31),
Figure 244215DEST_PATH_IMAGE003
and
Figure 734102DEST_PATH_IMAGE011
defined as the displacement value measured by the laser sensor C (33),
Figure 202124DEST_PATH_IMAGE002
and
Figure 17633DEST_PATH_IMAGE012
the displacement value measured by the laser sensor B (32) is defined, and the time 0 is indicated by a superscript "0", and the other arbitrary time is indicated by a superscript "n".
6. The calculation method of the underwater suspended tunnel test multi-dimensional motion pose measurement system according to claim 1, characterized in that: in said S5
Figure 624633DEST_PATH_IMAGE007
,
Figure 59157DEST_PATH_IMAGE008
The calculation method of (2) is as follows: defining the wave water of the pipe body respectively vertically upwards and horizontallyThe movement of the downstream is positive,
when in use
Figure 506319DEST_PATH_IMAGE013
Figure 102516DEST_PATH_IMAGE014
Then, the following equation is satisfied:
Figure 533497DEST_PATH_IMAGE015
when in use
Figure 302870DEST_PATH_IMAGE013
Figure 11063DEST_PATH_IMAGE014
Then, the following equation is satisfied:
Figure 902796DEST_PATH_IMAGE016
in the formula, Di,R,
Figure 696440DEST_PATH_IMAGE009
Figure 331820DEST_PATH_IMAGE002
Figure 425678DEST_PATH_IMAGE003
Figure 160416DEST_PATH_IMAGE010
Figure 907268DEST_PATH_IMAGE017
Figure 611918DEST_PATH_IMAGE018
All the measured values are known quantities, and the x-direction relative displacement value of the pipe body at each measured cross section at any n moments is calculated and obtained by simultaneous equation sets and the Newton iteration method and the least square method principle
Figure 825862DEST_PATH_IMAGE007
And z-direction relative displacement value
Figure 200343DEST_PATH_IMAGE008
7. The calculation method of the underwater suspended tunnel test multi-dimensional motion pose measurement system according to claim 1, characterized in that: the calculation method in S6 is as follows, when the displacement at several known cross sections is obtained, the displacement and attitude data at any cross section along the tunnel at any time can be obtained by interpolation method,
Figure 562054DEST_PATH_IMAGE019
Figure 476920DEST_PATH_IMAGE020
in the formula (I), the compound is shown in the specification,
Figure 748633DEST_PATH_IMAGE021
and
Figure 418649DEST_PATH_IMAGE022
respectively is the x-direction and z-direction relative displacement value alpha of the y coordinate position along the tunnel at any timekAnd betakRespectively at any time along the tunnelkThe rotation angle of the tube body around the z-axis and the x-axis is at the coordinate.
CN202210125047.5A 2022-02-10 2022-02-10 Multi-dimensional motion pose measurement system and calculation method for underwater suspension tunnel test Active CN114152245B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210125047.5A CN114152245B (en) 2022-02-10 2022-02-10 Multi-dimensional motion pose measurement system and calculation method for underwater suspension tunnel test

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210125047.5A CN114152245B (en) 2022-02-10 2022-02-10 Multi-dimensional motion pose measurement system and calculation method for underwater suspension tunnel test

Publications (2)

Publication Number Publication Date
CN114152245A CN114152245A (en) 2022-03-08
CN114152245B true CN114152245B (en) 2022-04-29

Family

ID=80450352

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210125047.5A Active CN114152245B (en) 2022-02-10 2022-02-10 Multi-dimensional motion pose measurement system and calculation method for underwater suspension tunnel test

Country Status (1)

Country Link
CN (1) CN114152245B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969068A (en) * 2014-04-11 2014-08-06 中国科学院力学研究所 Method and device for simulating axial interaction between undersea pipe system structure and seabed soil body
CN105091858A (en) * 2015-08-02 2015-11-25 上海砺晟光电技术有限公司 Two-dimension inclination angle non-contact measurement method and system based on absolute distance measurement
CN109297427A (en) * 2018-10-17 2019-02-01 北京市政建设集团有限责任公司 A kind of portable laser monitoring instrument applied to tunnel convergence Deformation control

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN203069788U (en) * 2012-12-10 2013-07-17 中山大学 Orientation and three-dimensional attitude monitoring display system for floating and semisubmersible moving platform
CN108827448B (en) * 2018-07-28 2023-09-26 天津大学 Vibration and inclination angle measurement system and method based on plane reflector and photoelectric array
CN110044268B (en) * 2019-03-19 2021-08-31 天津大学前沿技术研究院有限公司 Shield tunnel joint opening and slab staggering monitoring system based on optical fiber reflection principle
CN110108469B (en) * 2019-05-28 2022-02-15 中交公路规划设计院有限公司 Suspended tunnel pipe section attitude measurement device, test system and test method
CN209927400U (en) * 2019-06-26 2020-01-10 广西大学 Suspension tunnel dynamic response test device under earthquake and wave current coupling action
KR20210026033A (en) * 2019-08-29 2021-03-10 한국과학기술원 The external collision monitoring system of a submerged floating tunnel
CN113218621A (en) * 2021-06-09 2021-08-06 招商局重庆交通科研设计院有限公司 Suspension tunnel dynamic response test device and method under solid migration and wave flow coupling
CN113589227B (en) * 2021-07-14 2022-03-29 飞马滨(青岛)智能科技有限公司 Underwater robot positioning system and method
CN114088344B (en) * 2022-01-19 2022-04-01 交通运输部天津水运工程科学研究所 Test device and method for measuring dynamic hydrodynamic load parameters of suspension tunnel

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103969068A (en) * 2014-04-11 2014-08-06 中国科学院力学研究所 Method and device for simulating axial interaction between undersea pipe system structure and seabed soil body
CN105091858A (en) * 2015-08-02 2015-11-25 上海砺晟光电技术有限公司 Two-dimension inclination angle non-contact measurement method and system based on absolute distance measurement
CN109297427A (en) * 2018-10-17 2019-02-01 北京市政建设集团有限责任公司 A kind of portable laser monitoring instrument applied to tunnel convergence Deformation control

Also Published As

Publication number Publication date
CN114152245A (en) 2022-03-08

Similar Documents

Publication Publication Date Title
US10697670B2 (en) Heliostat surface shape detection system and method based on multi-view image recognition
CN106886024B (en) Deep-sea multi-beam sound ray precise tracking method
KR100421428B1 (en) Measurement of 6-degree-of-freedom motions of rigid bodies using a 3-facet mirror
CN109436197B (en) Coupling motion and force measurement simulation experiment system of ocean floating structure under action of internal solitary wave
CN106197312B (en) A kind of settled date mirror surface-shaped rapid detection system and its method
CN110127078B (en) Optical fiber monitoring method for strain-deflection-bending moment state of helicopter blade structure
CN110412313A (en) A kind of scaling method of ship true wind measuring device
CN101706275A (en) Measuring device and method of dynamic ocean waves
Winkelmann et al. Wind actions on large-aperture parabolic trough solar collectors: Wind tunnel tests and structural analysis
CN110108469B (en) Suspended tunnel pipe section attitude measurement device, test system and test method
CN103197096A (en) Wind power plant ultrasonic wind speed detection method and device
CN114357571A (en) Inversion method and system for atmospheric boundary layer wind field characteristics in constructed building environment
CN114541264B (en) Cable-stayed bridge cable guide pipe integrated measurement positioning method
CN105043333A (en) Miniaturized underwater manipulator position angle measuring method
CN107145647A (en) A kind of ship ocean surface wind speed, wind direction measurement data deviation correction method
CN111288912A (en) Fiber bragg grating deformation measurement method for airborne distributed POS
CN114152245B (en) Multi-dimensional motion pose measurement system and calculation method for underwater suspension tunnel test
CN105758602A (en) Truss girder bridge section buffeting force synchronous measurement method
CN110470860B (en) Time difference method based ultrasonic anemometer calibration method
CN109543227A (en) Charming appearance and behaviour error in pointing method for real-time measurement, the radar antenna of large-scale reflector antenna
CN109344537A (en) A kind of dynamic displacement acquisition methods based on Kalman Filtering Data Fusion
CN217787131U (en) Differential fiber grating vector flow velocity sensor
CN109596295B (en) Laser absolute calibration device for linear vibration of multi-axial vibration table
CN106405147B (en) A kind of ultrasonic transducer surveys wind array and its wind detection method
CN212255061U (en) Ocean floating platform carbon dioxide flux measuring device

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