Advanced
A study of the kinematic characteristic of a coupling device between the buffer system and the flexible pipe of a deep-seabed mining system
A study of the kinematic characteristic of a coupling device between the buffer system and the flexible pipe of a deep-seabed mining system
International Journal of Naval Architecture and Ocean Engineering. 2014. Sep, 6(3): 652-669
Copyright © 2014, The Society of Naval Architects of Korea
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Published : September 30, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Jae-Won Oh
Technology Center for Offshore Plant Industries, KRISO, Daejeon, Korea
Chang-Ho Lee
Technology Center for Offshore Plant Industries, KRISO, Daejeon, Korea
Sup Hong
Technology Center for Offshore Plant Industries, KRISO, Daejeon, Korea
Dae-Sung Bae
Mechanical Engineering Institute, Hanyang University, Seoul, Korea
Hui-Je Cho
Research & Development Institute, Virtual Motion, Inc, Seoul, Korea
Hyung-Woo Kim
Technology Center for Offshore Plant Industries, KRISO, Daejeon, Korea
Abstract
This paper concerns the kinematic characteristics of a coupling device in a deep-seabed mining system. This coupling device connects the buffer system and the flexible pipe. The motion of the buffer system, flexible pipe and mining robot are affected by the coupling device. So the coupling device should be considered as a major factor when this device is designed. Therefore, we find a stable kinematic device, and apply it to the design coupling device through this study. The kinematic characteristics of the coupling device are analyzed by multi-body dynamics simulation method, and finite element method. The dynamic analysis model was built in the commercial software DAFUL. The Fluid Structure Interaction (FSI) method is applied to build the deep-seabed environment. Hydrodynamic force and moment are applied in the dynamic model for the FSI method. The loads and deformation of flexible pipe are estimated for analysis results of the kinematic characteristics .
Keywords
INTRODUCTION
Many concepts for the commercial production of deep-seabed manganese nodules have been studied from the 1970s ( Brink and Chung, 1982 ; Chung, 1996 ; Herrouin et al., 1989 ; Amann et al., 1991 ; Liu and Yang, 1999 ; Hong and Kim, 1999 ; Deepak et al., 2001 ; Handschuh et al., 2001 ). The accumulate ground of the deep-seabed has a problem, in that the bearing capacity of the ground is not strong, because the accumulate ground is formed by fine particles with high moisture content. So it is impossible to carry manganese nodules in the collection system. Therefore, the validity of continuous mining by lifting pipe from ground to vessel is highly appreciated.
A continuous mining system is composed of mining vessel, lifting pipe, buffer system, transfer tube (flexible pipe), and self-propelled mining robot. The shape of the transfer tube between the buffer system and mining robot has a big impact on the driving efficiency of the mining robot. Also, the relative position between the buffer system and mining robot influences the efficiency of the mining robot. So, dynamic analysis of the integration of the mining system (mining vessel-lifting pipe-buffer system-flexible pipe-mining robot) is a very important technique in building a deep-sea mining system.
Currently, the lifting pipe and the flexible pipe are actively studied, according to the growth of offshore plant and the ocean floor industry. But studies on the buffer system are at an early stage. Just several functions to temporarily store nodules were mentioned by researchers ( Chung, 2003 ; Kotlinski et al., 2008 ).
The coupling device of the buffer system is the main subject of this study. The buffer system is currently developing. This coupling device makes a big impact on the dynamic movement of a mining robot. And it is a very important mechanic device in a deep-seabed mining system. The specs of the coupling device are determined by various efficiency tests and changes of design. But the production and experiment of a buffer system is costly and time-consuming. So the design specs of this device must be found by simulation. In general, this design method is called the simulation-based design.
Dynamic analysis of mechanical systems using computers is rapidly performed through the growth of computing power. The method of simulation-based design is a useful technique in otherwise impossible cases, which is verified by using an experiment with a model as an integrated deep-seabed mining system. The simulation technique is an excellent means of understanding qualitative (or quantitative) optimization design, and can skip the process of test model production, which is costly and time-consuming.
The mining robot and vehicle model used in this study is an MBD model. This MBD model was developed by Kim ( Kim et al., 2010 ). The integrated simulation model is developed by DAFUL (2012) .
In this study, the used equations are the joint constraints, a beam elastic equation and a multi-body dynamic solution. The verifications of the joint constraints and the beam elastic equation are written in DAFUL verification manual (2012) .
PPT Slide
Lager Image
Conceptual diagram of the deep-seabed mining system.
GOVERNING EQUATION
- Joint constraint
In this study, joint constraints ( Haug, 1989 ) were applied to find optimum kinematic characteristics of the coupling device. The joint constraints are as follows.
- Fixed constraint
A fixed constraint can remove motion for all degrees of freedom of rigid bodies or flexible bodies. Eqs. (1) and (2) show the fixed constraint equations ( Haug, 1989 ).
PPT Slide
Lager Image
PPT Slide
Lager Image
where, i is the base body’s number, j is the action body’s number, A is the orientation matrix of bodies, and f , g, h are axis components of matrix A .
Eq. (1) means a positional constraint, and Eq. (2) means a rotational constraint between the base body and action body. Eq. (1) is equal to a spherical constraint.
- Revolute constraint
A revolute constraint is able to allow rotation between two bodies about one rotation axis. Eqs. (1) and (3) show the revolute constraint equations ( Haug, 1989 ).
PPT Slide
Lager Image
where, i is the base body’s number, j is the action body’s number, A is the orientation matrix of bodies, and f , g , h are axis components of matrix A .
- Beam elastic
Beam elastic theory is used to build flexible pipe and transfer tube. Euler and Timoshenko had verified the beam method using beam stiffness. The equation of beam stiffness is as in Eq. (4).
PPT Slide
Lager Image
The beam stiffness matrix is a symmetric matrix. The components are as follows.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
where, E is Young’s modulus, G is the shear modulus, A is the cross-sectional area, L is the length, I is the area moment of inertia, and P is the coefficient by shear force. Commonly, P ’s value is zero when the beam stiffness is calculated.
The calculated stiffness is applied to the matrix force in DAFUL (2012) . This matrix force can apply forces of 3 translational and 3 rotational directions in a rigid body. These forces are dependent on the displacements and velocities of element bodies.
- Multi-body dynamic
All mechanical systems are composed as an assembly of many objects. This assembly has been called a multi-body. Multibody solutions are different from basic dynamic solutions. Multi-body solutions must find dynamic responses of each object, and exact solutions must be obtained by applying components, which are in correlation. Each body is defined to calculate the dynamic responses of a multi-body. The method to define an object is in two ways. First, an orthogonal coordinate system is attached to each body, and a motion of the body coordinate system is expressed by a generalized coordinate system. Second, connected bodies in series are numbered in couplers of the multi-body.
The body’s equation of motion is defined as a differential equation. The multi-body system’s equation of motion has nonlinearity of relative coordinate, velocity and acceleration.
An eidetic method to change from nonlinearity to linearity is that the equation is expressed by only the relative coordinate, velocity and acceleration. And relative coordinates must be expressed as equations of independent coordinates, velocities and accelerations. An interaction formula between an independent coordinate and dependent coordinate is able to lead from a constraint equation of position, velocity and acceleration.
The equation of motion is defined as follows ( Haug, 1989 ). A multi-body system can be modeled as a generalized coordinate of ‘ngc’ dimensions. These generalized coordinates are expressed as Eq. (13). Also, the velocity and acceleration coordinate are defined as Eqs. (14) and (15) respectively.
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
If a system has independent constraint equations of ‘ ncn ’ dimensions, this is expressed as Eq. (16). An acceleration equation (Eq. (17)) for the constraint conditional function is defined by a twice differential of Eq. (16) about the time domain.
PPT Slide
Lager Image
PPT Slide
Lager Image
where a is the acceleration vector of rigid bodies, and Φ q is the Jacobi matrix of the constraint conditional function about a generalized coordinate.
The equation of motion for a multi-body system having constraint is as in Eq. (18).
PPT Slide
Lager Image
where M is the mass matrix, λ is the Lagrange multipliers, and Q is the general force.
The equation of motion about the system is changed as a matrix, using Eqs. (17) and (18).
PPT Slide
Lager Image
The position, velocity and acceleration of a rigid body can be solved by these.
DYNAMIC MODEL OF DEEP-SEABED MINING SYSTEM
- Buffer system
A buffer system is a mechanical structure to save collected mineral lumps. A mass of the buffer system is 8 ton , a height is 7 m . The buffer system’s structural components are as in Table 1 .
Components of buffer system.
PPT Slide
Lager Image
Components of buffer system.
A flexible pipe in a buffer system provides passage to transfer minerals. The pipe that has generally been used is a PVC (polyvinyl chloride) pipe. The characteristics of the PVC pipe are given in Table 2 .
Characteristics of PVC pipe.
PPT Slide
Lager Image
Characteristics of PVC pipe.
The buffer and coupling device are built as a rigid body, to neglect the material deformation.
1) A swivel joint is built between the buffer and ground. The swivel joint is used to constrain the buffer’s rotation. Components of the swivel joint are the revolute joint and stopper. In this study, the stopper is built using a rotational spring damper and stiffness coefficient by the curved radius of the buffer. Also, the effect of the vessel is applied to the buffer.
2) Three constraints are used to research the kinematic characteristic of the coupling device between the flexible pipe in the buffer and the coupling device. They have 6 D.O.F, 5 D.O.F, and 3 D.O.F. Fig. 2 shows their constraints.
PPT Slide
Lager Image
Constraint between the coupling device and flexible pipe.
A buffer system model to perform numerical analysis is developed as in Fig. 3 .
PPT Slide
Lager Image
Buffer system model.
- Integration simulation model of mining system
An integration model is built with the buffer system model, mining robot and transfer tube. They are shown in Table 3 . The mining robot model is developed as in the Kim’s last study ( Kim et al., 2010 ). The transfer tube is built by beam elastic theory, to advance solving time. In common, dynamic responses by the beam model are equal to the results by the finite element method. But the solving time is very fast. The characteristics of the transfer tube are shown in Table 4 .
Components of the integration simulation model.
PPT Slide
Lager Image
Components of the integration simulation model.
Characteristics of the transfer tube.
PPT Slide
Lager Image
Characteristics of the transfer tube.
This transfer tube is affected by the fluid drag force, owing to its long length. So the drag force must be applied to the tube. Buoyancy modules are provided to prevent twisting and bending at particular positions of the transfer tube in the seabed environment. The characteristics of the buoyancy module are shown in Table 5 , and the positions of the buoyancy module are shown in Fig. 4 .
Characteristics of the buoyancy module.
PPT Slide
Lager Image
Characteristics of the buoyancy module.
PPT Slide
Lager Image
Positions of the buoyancy module.
Also, the inner flows generate in the pipe. The inner flows are a very complex chaotic phenomenon. Many studies are required to analyze the impact between nodules, the momentum change, the contact phenomenon between nodules and wall of the inner pipe and the expansion of compressed air. But the forces by inner flows are much smaller than the external forces. When we design the structures the effects of inner flows are not considered due to this reason. So in this study, inner flows are not applied.
The equation of the fluid drag force is as follows.
PPT Slide
Lager Image
where CD is the drag coefficient, ρw is the water density, A is the area of pipe, VB is the pipe velocity and VC is the current velocity. But the current velocity ( VC ) is close to zero in the deep-seabed environment.
The buoyancy equation is as follows.
PPT Slide
Lager Image
where ρ is the water density, V is the body volume and g is gravity.
Fig. 5 shows the developed integration simulation model.
PPT Slide
Lager Image
The integration simulation model.
SIMULATION CONDITION
- Forced motion
Forced motion is the basic test condition to estimate loads. The test conditions for movements of the coupling device are as follows.
Forced motion condition.
PPT Slide
Lager Image
Forced motion condition.
A maximum rotation angle of the coupling device for design is 20 deg . So a relative rotating angle of the coupling device is rotated from 0 degree to 20 degree . Drivers ( Haug, 1989 ) to rotate the device are applied at rotation 1 and 2 constraints.
PPT Slide
Lager Image
Driver function for rotation of the coupling device.
The rotation axes of conditions 1 and 2 are shown in Fig. 7 .
PPT Slide
Lager Image
Rotation axes of conditions 1 and 2 : (a) Rotation axis of condition 1, (b) Rotation axis of condition 2.
Total simulation time is 1.2 sec . and step size is 0.01.
- Driving motion
When the vessel and mining robot are driven, motions of the coupling device may be different from the results of the forced motion condition. Load estimation about the combined model must be performed. But the driving effect of the vessel is applied to the buffer system. Its effect is irrelevant to the movement of the coupling device. And the solving time is advanced, due to it being a simple model.
In this study, the driving conditions of the vessel and mining robot are determined by the results of Brink (1981), and Hong and Kim (2008) . Detailed driving conditions for the load estimation are shown in Table 7 .
Driving conditions for load estimation.
PPT Slide
Lager Image
ISP : Initial stationary state motion CTP : Translational motion by constant velocity of vessel and miner SP : Steering motion of miner DTP : Translational motion by different velocity of vessel and miner
NUMERICAL RESULT AND LOAD ESTIMATION
- Load estimation by forced motion
- Condition 1
Loads are estimated at the position of maximum stress and load in the flexible pipe. In figures, the color means size of stress. The red color is maximum stress of the pipe and the blue color is minimum stress. The position is shown in Fig. 8 .
PPT Slide
Lager Image
Measurement position for load estimation in condition 1.
Measured loads of the flexible pipe are different according to the kinematic constraints in condition 1. The loads are as follows.
PPT Slide
Lager Image
Load estimation of the flexible pipe according to constraints of the coupling device in condition 1.
Maximum load of flexible pipe in condition 1.
PPT Slide
Lager Image
Maximum load of flexible pipe in condition 1.
In fixed constraint, stress of pipe has come close to the yield strength. But other cases have not come close to the yield strength. When the coupling device has a rotational degree of freedom, the stability of the flexible pipe is better than with no rotation.
- Condition 2
Loads are estimated at the position of maximum stress and load in the flexible pipe. The position is shown in Fig. 10 .
PPT Slide
Lager Image
Measurement position for load estimation in condition 2.
The measured loads of the flexible pipe are different according to the kinematic constraints in condition 2. The loads are as follows. The loads in fixed & revolute constraints are exceeded the yield strength. These facts mean the breakdown of flexible pipe. When the kinematic constraint is fixed constraint, the load of the pipe has bigger than revolute constraint. And the tendency of load has been dramatic changes. Also, the revolute constraint has bigger than spherical constraint. This means that rotational degree of freedom is very important for safe design.
PPT Slide
Lager Image
Load estimation of the flexible pipe according to constraints of the coupling device in condition 2.
Maximum load of flexible pipe in condition 2.
PPT Slide
Lager Image
Maximum load of flexible pipe in condition 2.
We knew that deformations of the flexible pipe differ according to kinematic constraints of the coupling device. Fig. 12 shows their deformations. When the coupling device has rotational degrees of freedom, the load of the flexible pipe is smaller than with no rotation. Snapping of the pipe did not occur, so the mineral can be transferred easily. Also, when this device has 3 D.O.F for rotation, the stress of the flexible pipe keeps within the yield strength.
PPT Slide
Lager Image
Deformation of the flexible pipe by kinematic characteristic of the coupling device : (a) Fixed and revolute joint, (b) Spherical joint.
- Condition 3
In this case, measurement positions for the load estimation are not the same, because the maximum load and stress are applied by constraints at different points. The positions are shown in Fig. 13 .
PPT Slide
Lager Image
Measurement positions for load estimation of the flexible pipe.
The measured loads of the flexible pipe differ according to the kinematic constraints in condition 3. The loads are as follows. The tendency of results is similar to the condition 2. This means that rotational degree of freedom is very important for safe design in condition 3.
PPT Slide
Lager Image
Load estimation of the flexible pipe according to constraints of the coupling device in condition 3.
Maximum load of flexible pipe in condition 3.
PPT Slide
Lager Image
Maximum load of flexible pipe in condition 3.
The result by fixed joint is much larger than the result by spherical joint. This result shows that the structural stability differs according to the kinematic characteristics. And the bending phenomenon is expressed in the fixed joint and revolute joint. But it is not expressed in the spherical joint. Therefore, the fixed joint of kinematic characteristics is a very bad constraint. So result analysis by fixed joint is excluded in the integration analysis.
- Condition 4
Loads are estimated at the position of maximum stress and load in the flexible pipe. The position is shown in Fig. 15 .
PPT Slide
Lager Image
Measurement position for load estimation of the flexible pipe.
The measured loads of the flexible pipe differ according to the kinematic constraints in condition 4. The loads are as follows. The tendency of results is not similar to the condition 2. This means that the motion of rotation axis 2 had given a significant effect on the load.
PPT Slide
Lager Image
Load estimation of the flexible pipe according to constraints of the coupling device in condition 4.
Maximum load of flexible pipe in condition 4.
PPT Slide
Lager Image
Maximum load of flexible pipe in condition 4.
In fixed and revolute constraints, stress of the pipe has come close to the yield strength. But spherical constraint has not come close to the yield strength. When the coupling device has rotational degrees of freedom, the stability of the flexible pipe is better than with no rotation.
- Condition 5
Loads are estimated at the position of maximum stress and load in the flexible pipe. The position is shown in Fig. 17 .
PPT Slide
Lager Image
Measurement position for load estimation of the flexible pipe.
The measured loads of the flexible pipe differ according to the kinematic constraints in condition 5. The loads are as follows.
PPT Slide
Lager Image
Load estimation of the flexible pipe according to constraints of the coupling device in condition 5.
Maximum load of flexible pipe in condition 5.
PPT Slide
Lager Image
Maximum load of flexible pipe in condition 5.
In fixed, stress of the pipe has come close to the yield strength. But other constraints have not come close to the yield strength. When the coupling device has rotational degrees of freedom, the stability of the flexible pipe is better than with no rotation. Also, we can know that the rotational degree of freedom by rotation 2 is very critical parameter to prevent damage.
- Load estimation by driving motion
The load results for each driving condition are as follows.
Load estimation by driving condition.
PPT Slide
Lager Image
Load estimation by driving condition.
In all conditions, the spherical joint is better than the revolute joint. And when kinematic characteristic is spherical constraint, the stress of pipe doesn`t exceed the yield strength. This result shows that the rotational degree of freedom is a very important parameter in the design of a coupling device.
In DTP1 and DTP2 driving conditions, the loads of the flexible pipe are distinguished according to their kinematic characteristics. It is known that the driving velocities of the vessel and mining robot are important. Research into this behavior will be conducted in the future.
Load estimation in DTP1 of driving condition.
PPT Slide
Lager Image
Load estimation in DTP1 of driving condition.
Load estimation in DTP2 of driving condition.
PPT Slide
Lager Image
Load estimation in DTP2 of driving condition.
CONCLUSION
Through this study, a concept of the simulation-based design by using multi-body dynamics is introduced through the kinematic design of the buffer system. The results of this study are as follows.
1) Kinematic design of the coupling device of the buffer system
The loads of the flexible pipe depend on the kinematic characteristics of the coupling device. The size of the loads as follows.
Fixed constraint >> Revolute constraint >> Spherical constraint
When the coupling device has many the rotational degrees of freedom, the load of the pipe is very small. This result shows that the rotation of the coupling device is dominated with the stability of buffer system. Thus, the rotational degree of freedom is very important design parameter for the buffer system design.
The bending and twisting of the flexible pipe do not appear when there is the rotational degree of freedom. This enables the stable transfer of the minerals. And only, the result of the spherical constraint is not exceeded the yield strength in all conditions. The coupling device that has the rotational degrees of freedom must be produced for stable transfer of mineral and the stability of the buffer system. Also, we can be known that the rotational degree of freedom is a factor affecting equipment stability in operation of the vessel and the mining robot.
2) Simulation-base design
We can be known that the simulation-based design is the useful design process for a producing of deep-seabed equipment. The simulation-based design method is based on the multi-body dynamics, the nonlinearity is considered. For this reason, this design method considers influences by motions and forces of the actual equipment.
Yet, the theory on the simulation-based design is not formulated perfectly.
Finally, the potential of the simulation-based design is confirmed in this paper. In the future, we will formulate the theory of the simulation-based design method and design other mechanic systems of deep-seabed mining system by using the simulation-based design.
Acknowledgements
This study was initiated by an R&D Project, “Technology Development of Deep-Seabed Mining System for Manganese Nodules”, sponsored by the Ministry of Oceans and Fisheries of Korea. The authors are grateful for the full support shown for this research work.
References
Amann H. , Oebius H.U. , Gehbauer F. , Schwarz W. , Weber R. 1991 Soft ocean mining Proceeding Offshore Technology Conference Houston 5 June 1991 1 - 12
Brink A.W. , Chung J.S. 1982 Automatic position control of a 300,000-ton ship ocean mining system Journal of Energy Resource Technology 104 (4) 285 - 293
Chung J.S. 1996 Deep-ocean mining. technologies for manganese nodules and crusts International Journal of Offshore and Polar Engineering 6 (4) 244 - 254
Chung J.S. 2003 Deep-ocean mining technology: learning curve I. Proceedings of Fifth ISOPE Ocean Mining Symposium (ISOPE OMS-2003) Tsukuba, Japan 15-19 September 2003 1 - 5
DAFUL 4.2, User's manual Virtual Motion, Inc. Seoul
DAFUL 4.2, 2012. Verification manual Virtual Motion, Inc. Seoul
Deepak C.R. , Shajahan M.A. , Atmanand M.A. , Annamalai K. , Jeyamani R. , Ravindran M. 2001 Development tests on the underwater mining system using flexible riser concept Proceedings Fourth ISOPE Ocean Mining Symposium Szczecin, Poland 23-27 September 2001 94 - 98
Handschuh R. , Grebe H. , Panthel J. , Schulte E. , Wenzlawski B. , Schwarz W. , Atmanand M.A. , Jeyamani R. , Shajahan M. , Deepak C. , Ravindran M. 2001 Innovative deep ocean mining concept based on flexible riser and self-propelled mining machines Proceedings Fourth ISOPE Ocean Mining Symposium Szczecin, Poland 23-27 September 2001 99 - 107
Haug E.J. 1989 Computer-aided kinematics and dynamics of mechanical systems Allyn and Bacon Massachusetts basic methods. Needham Heights
Herrouin G. , Lenoble J.P. , Charles C. , Mauviel F. , Bernard J. , Taine B. 1989 A manganese nodule industrial venture would be profitable. summary of a 4-year study in France Proceedings Offshore Technology Conference Houston 5 January 1 - 12
Hong S. , Kim K. 1999 Proposed technologies for mining deep-seabed polymetallic nodules -chap 12 research and development of deep seabed mining technologies for polymetallic nodules in Korea Proceedings International Seabed Authority’s Workshop Kingston, Jamaica 3-6 August 1999 261 - 283
Hong S. , Kim H.W. 2008 Coupled dynamic analysis of underwater tracked vehicle and long flexible pipe Journal of the Korean Society of Oceanography 13 (3) 237 - 245
Kim H.W. , Hong S. , Lee C.H. , Choi J.S. , Yeu T.K. 2010 A study on steering characteristics of four-row tracked vehicle on extremely cohesive soft soil Proceedings of the 9th Asia-Pacific ISTVS Conference Sapporo, Japan 27-30 September 2010 1 - 4
Kotlinski R. , Stoyanova V. , Hamrak H. , Avramov A. 2008 An overview of the interoceanmetal deep-sea technology development (Mining and Processing) programme Proceeding of a workshop held by ISA Chennai, India 168 - 184
Liu F. , Yang N. 1999 Proposed technologies for mining deep-seabed polymetallic nodules - chap 9 environmentally friendly deep seabed mining system Proceedings International Seabed Authority’s Workshop Kingston, Jamaica 3-6 August 187 - 211