Distributed Sensor Networks for Health Monitoring of Civil Infrastructures
View this Special IssueResearch Article  Open Access
Anchor Dragging Analysis of RockBerm Using Smoothed Particle Hydrodynamics Method
Abstract
This study presents dynamic responses of rockberm structural system under anchor dragging and accordingly provides the characteristics of the stresses and displacements obtained. For the purpose, first, a rockberm was modeled by the SPH (smoothed particle hydrodynamics) method and piecewise DruckerPrager material model by facilitating the associated software package—ANSYSAUTODYN. Second, 2ton stockless anchor was modeled as a rigid body and eventually dragging external force was obtained. Then, the dragging velocity (1 and 2 m/s) was considered as a parameter to investigate the effect of its variation on the responses. Finally, the dragging tensile forces of the anchor cable were obtained and compared according to the dragging velocities. It is shown that the fourlayer rockberm gives the safety margin to the submarine power cable according to the unaffected gauge points near the cable. This safety is accomplished by the four layers (related to rockberm height) and the number of rock particles at each layer (related to rockberm widths).
1. Introduction
Submarine power cables have been used to connect dots and accordingly transmit electric power from the main land to an island. Recent developments in ocean renewable energies make the conventional energy flows in another direction such as offshore wind farms to the main land. This innovation requires the frequent use of submarine power cables; hence, their safety assurance methods have been developed.
Among the methods, rockberms have been used in relatively deep waters because of their easy installation and inexpensive cost. A recent study showed that rockberm height and rock size are important parameters against an accidental load—anchor collision [1]. Along with anchor collision, anchor dragging is also important for the protecting structures [2]. Yoon and Na [2] investigated the dynamic responses of a flexible concrete mattress through field experiments. However, no one has investigated the behaviors of rockberms under anchor dragging so far.
This fact extends to numerical studies because the conventional Lagrange method is not applicable in the modeling of rocks [3] but applicable in tunneltype protective structures [4, 5]. Accordingly, it is demanded to develop a numerical procedure, to evaluate the dynamic behaviors of a rockberm under anchor dragging, and to acquire the safety assessment of submarine power cables, which are installed under rockberms.
This study presents dynamic responses of rockberm behaviors under anchor dragging and accordingly provides the characteristics of the stresses and displacements obtained. For the purpose, first, a rockberm was modeled by the SPH (smoothed particle hydrodynamics) method and piecewise DruckerPrager material model by facilitating the associated software package—ANSYSAUTODYN [6]. Second, 2ton stockless anchor was modeled as a rigid body and eventually dragging external force was obtained. Then, the dragging velocity (1 and 2 m/s) was considered as a parameter to investigate the effect of its variation on the responses. Finally, the dragging tensile forces of the anchor cable were obtained and compared according to the dragging velocities.
It should be noted here that this study is limited to investigating the dynamic responses of the target structure under anchor dragging. Other issues such as monitoring of the transient responses [7], multistage structural damage diagnosis [8], and more detailed signal process [9] should be considered in further research activities.
2. Materials and Methods
The SPH method was initially developed by Gingold and Monaghan [10] and Lucy [11] for astrophysical problems. It has been used in many research fields such as astrophysics, ballistics, volcanology, and oceanography because it is a meshfree Lagrange method and resolution can easily be adjusted according to variables such as the density. Most initial problems are fluidrelated by dividing the fluid into a set of discrete elements (particles). Libersky and Petschek [12] extended the SPH to solid mechanics. The main advantage of SPH in solids is the possibility of dealing with the larger local distortion than gridbased methods; hence, it has been exploited in many applications in solid mechanics such as metal forming, collision, crack growth, fracture, and fragmentation. The detailed description of the SPH method and associated references can be found in the literature [13]. Here, we introduce the brief concept of the SPH method.
The SPH method uses kernel approximation, which is based on randomly distributed interpolation points with no assumptions about which points are neighbors, to calculate spatial derivatives and particle density using Kernel estimation [14]. The particle density can be expressed using kernel approximation as follows: where is the density of th particle, is the mass of th particle, and is the weight function according to particle size . and are coordinates of th and th particle. In the SPH method, the particle density can be obtained as mass of neighbor particle multiplied by weight function. Figure 1 shows a schematic of the weight function in SPH method. In the study, the Kernel Bspline is used as the weight function. The basic form of weight function is as follows: where is a parametric variable defined by .
Two contact conditions were considered, as shown in Figure 2. First, the interface between rock and other materials (basement and anchor) is based on the contact algorithms used in Lagrangeoriented finite element codes. Second, the interface between the stones (particle to particle contact) is based on the meshless approach [15]. If other material intersects the circle of a particle, then contact is assumed to be occurred. When contact occurs, the restoring force is applied along the normal to the surface. In the particle to particle contact, restoring forces are applied between particles.
(a) Particle to surface contact
(b) Particle to particle contact
Rockberm is a typical example on which the SPH method is applicable because rockberm is a structure constructed by stacking rocks to protect lifelines (e.g., submarine power cables and pipelines) and to build breakwaters. Figure 3 shows a rockberm consisting of rock particles (the upper layers) modeled by the SPH method and lower seabed layers modeled by 8node solid elements. The linear elastic material properties were used for modeling of the seabed layers, as shown in Table 1. The piecewise DruckerPrager material model was used for modeling of the rock layers because the material model is frequently used for the discrete materials such as concrete, rocks, and soils [16, 17]. In the rock material model, we used the material properties of Table 2, which is proposed by Chen et al. [18] and the associated pressureyield stress curve, shown in Figure 4.


The rockberm size was determined according to Cberm, which was proposed by the Korean Electric Power Corporation [19]. The bottom width is 7 m, the top width is 2 m, the height is 1.2 m, and the slope angle is 25.64° as shown in Figure 3. The diameters of the rocks typically range from 100 to 300 mm in the practice; hence, the upper bound 300 mm was used in the study. In reality, the rockberm extends to the and directions for the continuous layout; hence, two symmetry planes (the boundary planes A and B in Figure 3) were assumed and constrained in the horizontal ( and ) directions. In addition, the seabed bottom plane was constrained in the vertical direction.
The target anchor (2ton stockless) is one of anchors frequently used in commercial and other vessels in the Korean waters. It was modeled according to the Korea Standard [20] and modeled by 4node shell elements for the dragging analysis, as shown in Figure 5. Moreover, the anchor is assumed to be a rigid body because the density (7200 kg/m^{3}) and Young’s modulus (170 GPa) are much larger than those of the rocks and seabed. Table 3 shows the anchor considered and its representative dimensions [20]. As shown in Figure 3, the 2ton stockless anchor was guided (attached) by the anchor cable, which is modeled by line (link) elements (density: 7850 kg/m^{3} and Young’s modulus: 210 GPa). The cable is always under tension; hence, the conventional cable element which can resist the tension only is not necessary. Besides, the geometrical nonlinearity of the anchor cable was not considered because of the external force is tensiononly. The dragging location was assumed to be the bottom of the 2nd layer (Layer B) at the boundary plane A, as shown in Figure 6. Two dragging scenarios were made by considering two dragging velocities (1 and 2 m/s).

To get the dynamic responses, 11 gauge points are specified as shown in Figure 7. In the figure, each layer has its own name, for example, Layer A, Layer B, Layer C, and Layer D; hence, the gauge points are named to indicate the layers. For example, A1 denotes the first gauge point locating at Layer A. The coordinates of the gauge points are given in Table 4. It should be noted here that the location of submarine power cables is just below A1 because the position of the number “1” is the center of the rockberm (axis).

3. Results and Discussions
The vonMises stress (equivalent stress) is often used in determining whether a material will yield when subjected to a complex loading condition. Figure 8 shows the vonMises stresses at the gauge points (A1–D2) when dragging velocity is 1 m/s (first scenario). Initially, the stresses at A3 and B3 are distinct because the anchor dragging starts at the right side between Layer A and Layer B. The stresses at B3 are bigger than those at A3 due to the dragging direction (left or direction) and the bottom boundary condition. After approximately 1.5 sec, the stresses at C2 and D2 exceed the stresses at A3, which indicates that the anchor is moving up. Then, the stresses at C1 and D1 exceed the stresses at A3 after around 2.0 sec, which denotes that the anchor reaches at the top layer (Layer D). It should be emphasized that the stresses at A1 are almost zero; hence, the submarine power does not tends to experience stresses. Over the period, the maximum vonMises stress is approximately 4.3 kPa, which occurs at C2 and D2, as shown in Table 5.

Figure 9 shows the horizontal displacements (UX). Initially, the front particles (A3, A4, and B3) move to the left ( direction) according to the anchor dragging and then followed by C2 and D2 and accordingly D1. The horizontal displacements at the gauge points A1, A2, and B1 are almost zero, which indicate that the particles do not receive main effect from the dragging movement. The maximum UX are about 855 mm at D1. Figure 10 shows the vertical displacements (UY). Initially, the front particles (A3 and B3) move upward ( direction) and then followed by C1 and D1. An interesting factor is that the D2 experiences the downward movement ( direction) after around 2.0 sec. This means the corresponding particle is disintegrated. It is also shown that the vertical displacements at A1, A2, A4, and B1 are almost zero over the time interval, standing for less effective zone of the anchor dragging. The maximum UY is approximately 283 mm at D1.
To pinpoint the particle movement, Figure 11 shows the movement sequence according to time. From the figure, it is clearly observed that the stockless anchor penetrated, dragged, and disintegrated the rocks as time goes. However, it does not cause a significant response at A1, which locates at the top of the submarine power cable. Thus, the dragging scenario does not threaten the safety of the power cable. Figure 12 shows the tensile force of the anchor cable. The initial tensile force is 50 ton, then the force tends to decrease to 10 to 25 ton when the anchor penetrated the rocks (between 0.7 and 2.0 sec), and then the force increases to 10 to 35 ton when the anchor gives a deeper penetration into the rocks (after 2 sec). This late increase of the tensile force is caused by the upper layer rocks, which do not roll down from the rockberm as shown in Figure 11.
(a) 0.5 sec
(b) 1.0 sec
(c) 2.0 sec
(d) 3.0 sec
Figure 13 shows the vonMises stresses at the gauge points (A1–D2) when dragging velocity is 2 m/s. Initially, the stresses at A3, A4, and B3 are distinct because the anchor dragging starts at the right side between Layer A and Layer B. The stresses at B3 are bigger than those at A3 and A4 due to the dragging direction ( direction) and the bottom boundary condition. After approximately 1.3 sec, the stresses at C2 and D2 exceed the stresses at A3, which indicates that the anchor is moving up. Then, the stresses at C1 and D1 exceed the stresses at A3 after around 1.5 sec, which denotes that the anchor reaches at the top layer (Layer D). It should be emphasized that the stresses at A1 are almost zero; hence, the submarine power does not tend to experience stresses. Over the period, the maximum vonMises stress is approximately 4.5 kPa, which occurs at D1, as shown in Table 6. In comparison with the first scenario, the second scenario shows the similar characteristics but little earlier occurrence than those of 1 m/s because of the faster dragging velocity (2 m/s).

Figure 14 shows the horizontal displacements (UX) in the case of 2 m/s of dragging velocity. Initially, the front particles (A3, A4, and B3) move to the left ( direction) according to the anchor dragging and then followed by C2 and D2 and accordingly D1. The horizontal displacements at the gauge points A1, A2, and B1 are almost zero, which indicate that the particles do not receive main effect from the dragging movement. The maximum UX are about 710 mm at D1. Figure 15 shows the vertical displacements (UY) in the second scenario. Initially, the front particles (A3 and B3) move upward ( direction) and then followed by C1 and D1. An interesting factor is that the D2 experiences the downward movement ( direction) after around 1.7 sec. This means the corresponding particle is totally disintegrated. It is also shown that the vertical displacements at A1, A2, A4, and B1 are almost zero over the time interval, standing for less effective zone of the anchor dragging. The maximum UY is approximately 365 mm at B3, unlike the first scenario showing that the maximum UY occurs at D1. This indicates that higher dragging velocity causes larger vertical movement at the initial period. In comparison with the first scenario, the second scenario shows the similar characteristics but little earlier movement than those of 1 m/s because of the faster dragging velocity (2 m/s).
To pinpoint the particle movement, Figure 16 shows the movement sequence according to time. From the figure, it is clearly observed that the stockless anchor penetrated, dragged, and disintegrated the rocks as time goes. However, it does not cause a significant response at A1, which locates at the top of the submarine power cable. Thus, the dragging scenario does not threaten the safety of the power cable. In comparison with the first scenario, as shown in Figures 8 and 13, the higher dragging velocity causes much faster anchor movement, and accordingly the stockless anchor is already separated from the rockberm at 3.0 sec. Thus, more disintegration or disturbance of the rocks happens in the second scenario. Figure 17 shows the tensile force of the anchor cable. The initial tensile force is 65 ton, and then the force tends to decrease as time goes. The tensile force of the second scenario has a more distinct decreasing pattern than that of the first scenario (see the trend lines and their characters of Figures 12 and 17). This can be explained by the early escape of the anchor from the rockberm in the second scenario.
(a) 0.5 sec
(b) 1.0 sec
(c) 2.0 sec
(d) 3.0 sec
4. Conclusion
By facilitating the SPH method, it is shown that the global behavior of the rockberm is not much sensitive to the dragging velocity. This indicates that the gauge points give similar responses although the responding moments are little different according to the velocity variation. However, the faster dragging force (2 m/s) gives more rapid disintegration of the rockberm, and accordingly the late responses are relatively distinguishable. Moreover, considering the tensile force of the anchor cable, the higher velocity (2 m/s) gives a more distinct decreasing pattern of the tensile force than that of the lower velocity (1 m/s) as time goes because of the early escape of the anchor from the rockberm. The fourlayer rockberm gives the safety margin to the submarine power cable according to the unaffected gauge points near the cable. This safety is accomplished by the four layers (rockberm height) and the number of rock particles at each layer (related to rockberm widths and slope angle). Therefore, the dimensions (height = 1.2 m, lower width = 7 m, upper width 2 m, and slope angle = 25.64°) provide the safety of submarine power cables under the anchor dragging scenarios.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This research was a part of the project titled “Development on the Guidelines of Safety Assessment for Submarine Cable Protection Facilities in Shallow Water” funded by the Ministry of Land, Transport and Maritime Affairs, Korea (KIMT201220120018).
References
 J. Woo and W. B. Na, “Safety assessment of rockberms under anchor collision for protecting submarine power cables,” Ocean Engineering. In press. View at: Google Scholar
 H.S. Yoon and W.B. Na, “Safety assessment of submarine power cable protectors by anchor dragging field tests,” Ocean Engineering, vol. 65, pp. 1–9, 2013. View at: Publisher Site  Google Scholar
 J. Woo, D. Kim, H.S. Yoon, and W.B. Na, “Characterizing Korean general artificial reefs by drag coefficients,” Ocean Engineering, vol. 82, pp. 105–114, 2014. View at: Publisher Site  Google Scholar
 J. Woo and W. B. Na, “Analyses of the maximum response of cylindersconnected protectors under anchor colliding and dragging,” Journal of Ocean Engineering and Technology, vol. 24, pp. 81–87, 2010 (Korean). View at: Google Scholar
 J. Woo, W. B. Na, and H. T. Kim, “Numerical simulation of archtype submarine cable protector under anchor collision,” Journal of Ocean Engineering and Technology, vol. 23, pp. 96–103, 2009 (Korean). View at: Google Scholar
 ANSYS, AUTODYN User Manual Version 12.0, ANSYS, Canonsburg, Pa, USA, 2009.
 T. H. Yi, H. N. Li, and M. Gu, “A new method for optimal selection of sensor location on a highrise building using simplified finite element model,” Structural Engineering and Mechanics, vol. 37, no. 6, pp. 671–684, 2011. View at: Publisher Site  Google Scholar
 T.H. Yi, H.N. Li, and H.M. Sun, “Multistage structural damage diagnosis method based on “energydamage” theory,” Smart Structures and Systems, vol. 12, no. 34, pp. 345–361, 2013. View at: Publisher Site  Google Scholar
 H. Li, T. Yi, M. Gu, and L. Huo, “Evaluation of earthquakeinduced structural damages by wavelet transform,” Progress in Natural Science, vol. 19, no. 4, pp. 461–470, 2009. View at: Publisher Site  Google Scholar
 R. A. Gingold and J. J. Monaghan, “Smoothed particle hydrodynamics: theory and application to nonspherical stars,” Monthly Notices of the Royal Astronomical Society, vol. 181, pp. 375–389, 1977. View at: Google Scholar
 L. B. Lucy, “A numerical approach to the testing of the fission hypothesis,” The Astronomical Journal, vol. 82, pp. 1013–1024, 1977. View at: Google Scholar
 L. D. Libersky and A. G. Petschek, “Smooth particle hydrodynamics with strength of materials,” in Advances in the FreeLagrange Method Including Contributions on Adaptive Gridding and the Smooth Particle Hydrodynamics Method, vol. 395 of Lecture Notes in Physics, pp. 248–257, Springer, Berlin, Germany, 1991. View at: Publisher Site  Google Scholar
 J. Lin, H. Naceur, D. Coutellier, and A. Laksimi, “Efficient meshless SPH method for the numerical modeling of thick shell structures undergoing large deformations,” International Journal of NonLinear Mechanics, vol. 65, pp. 1–13, 2014. View at: Publisher Site  Google Scholar
 C. J. Hayhurst, R. A. Clegg, I. H. Livingstone, and N. J. Francis, “The application of SPH techniques in AUTODYN2D to ballistic impact problems,” in Proceeding of the 16th International Symposium on Ballistics, San Francisco, Calif, USA, 1996. View at: Google Scholar
 S. Seo and O. Min, “Axisymmetric SPH simulation of elastoplastic contact in the low velocity impact,” Computer Physics Communications, vol. 175, no. 9, pp. 583–603, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 W. K. Liu, S. Jun, S. Li, J. Adee, and T. Belytschko, “Reproducing kernel particle methods for structural dynamics,” International Journal for Numerical Methods in Engineering, vol. 38, no. 10, pp. 1655–1679, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 X. H. Zhu and Y. J. Jia, “3D mechanical modeling of soil orthogonal cutting under a single reamer cutter based on DruckerPrager criterion,” Tunnelling and Underground Space Technology, vol. 41, no. 1, pp. 255–262, 2014. View at: Publisher Site  Google Scholar
 S. G. Chen, J. Zhao, A. Makurati, and C. Madshus, “Discrete element modeling of an explosion test in granite,” in Proceedimgs of GeoEng2000: An International Conference on Geotechnical & Geological Engineering, 19–24 November 2000, Melbourne Exhibition and Convention Centre, Melbourne, Australia, pp. 1474–1481, 2000. View at: Google Scholar
 KEPCO, Field Test Report for Installation of Rock Berm in the Shallow Water of JindoJeju Island HVDC Interconnection Project, Korea Electric Power Corporation, 2010 (Korean).
 KS V 3311, Anchors, Korean Standards Association, 2012.
Copyright
Copyright © 2015 Jinho Woo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.