Mathematical Modelling of Bonded Marine Hoses for Single Point Mooring (SPM) Systems, with Catenary Anchor Leg Mooring (CALM) Buoy Application—A Review

The application of mathematical analysis has been an essential tool applied on Catenary Anchor Leg Mooring (CALM) buoys, Wave Energy Converters (WEC), point absorber buoys, and various single point mooring (SPM) systems. This enables having mathematical models for bonded marine hoses on SPM systems with application with CALM buoys, which are obviously a requisite for the techno-economic design and operation of these floating structures. Hose models (HM) and mooring models (MM) are utilized on a variety of applications such as SPARs, Semisubmersibles, WECs and CALM buoys. CALM buoys are an application of SPM systems. The goal of this review is to address the subject of marine hoses from mathematical modeling and operational views. To correctly reproduce the behavior of bonded marine hoses, including nonlinear dynamics, and to study their performance, accurate mathematical models are required. The paper gives an overview of the statics and dynamics of offshore/marine hoses. The reviews on marine hose behavior are conducted based on theoretical, numerical, and experimental investigations. The review also covers challenges encountered in hose installation, connection, and hang-off operations. State-of-the-art, developments and recent innovations in mooring applications for SURP (subsea umbilicals, risers, and pipelines) are presented. Finally, this study details the relevant materials that are utilized in hoses and mooring implementations. Some conclusions and recommendations are presented based on this review.


Introduction
The oil and gas sector is currently a challenging and dynamic environment, necessitating the exploration of additional subsea natural resource deposits in order to extract and meet global energy demand. Operators must develop and implement innovative methods of natural resource extraction, which offer more efficient and cost-effective extraction solutions, in order to stay ahead in this challenging industry. Furthermore, composite material utilization in bonded flexible hoses is still relatively new. There is currently minimal literature on the subject, necessitating the need for further investigations and examinations of the potential for general acceptability within the industry [1][2][3][4][5][6][7]. The application of mathematical analysis has been an essential tool applied in Catenary Anchor Leg Mooring (CALM) buoys, Single Anchor Leg Mooring (SALM) buoys, Wave Energy Converters (WEC), Paired Column Semisubmersibles (PCSemi), point absorber buoys, and single point mooring (SPM) systems. Thus, having mathematical models for bonded marine hoses on CALM buoys and SPM systems is obviously a requisite for the techno-economic design fundamental cause of this failure was discovered to be out-of-plane (OPB) bending-induced fatigue, which lowered the fatigue life of the chain links for the moorings by 95%. Due to the multiple parameters required in the formulations and the unpredictability of mooring setups, the methodology to integrate OPB/IPB fatigue for failure assessments is a difficult process. Thus, the need for Offshore Monitoring Systems (OMS). Notable HMs have been considered in recent studies as reported in Amaechi [110][111][112][113][114], and Tonatto [115][116][117][118][119] which considered the statics and dynamics of hoses in separate presentations. In the study by Szczotka [120], an HM was presented with a numerical model developed using a rigid finite element model (RFEM) for the static and dynamic reeling operation in pipe-laying. The study showed that by using this model, active and passive tensioner reels could be studied and the reeling installation of pipelines in sea conditions could be simulated well, as validated via a pipe laying simulation with an active reel drive. Analytical models of marine risers [121][122][123][124][125][126][127][128] used in developing the statics and dynamics of marine hoses, except that material behaviour also induces hose motion [129][130][131]. Conversely, there are different approaches on these mathematical presentations on MMs and HMs. The present study portrays different technical reports on the complexity of various HMs and MMs, hence the need for discussions on some governing mathematical equations, as presented in this review. A typical sketch on the mathematical model for a CALM buoy hose is shown in Figure 1.
Images for jmse-1295533 -Proof Correction sent by Amaechi, Chiemela Victor, 29/10/2021.  The goal of this work is to fill this gap by analyzing mathematical modeling of bonded marine hoses on SPM mooring systems, with the application using CALM buoys. To achieve this goal, this paper will identify various mathematical modeling methods for bonded marine hoses on SPM mooring systems with CALM buoy application, and review their use in the literature relating to its design configurations and analysis. The study aims also to mathematically determine the benefits and challenges of the modeling methods for different applications in SPM analysis, such as on CALM. In this review, Section 2 covers an overview of SPM moorings, Section 3 presents model methodologies and software tools, Section 4 presents mathematical models and other model types, Section 5 presents governing equations and motion characteristics, and Section 6 presents concluding remarks.

Single Point Mooring (SPM): An Overview
A single point mooring buoy is a buoy that is securely anchored to the seabed by several mooring lines/chains/anchors in a permanent position, permitting a liquid petroleum product cargo to be transferred. A bearing system on the buoy allows a portion of it to   Harnois [160,161] provided a comparison of various mooring line materials. The inertia, elastic stiffness, and damping of a mooring line are affected by the material used. An anchor's purpose is to secure a mooring line to a fixed place on the bottom. The ability to resist high, horizontal, and in some cases vertical loads in a specific seabed type (soft to hard), cost-effectiveness, and ease of installation are the major requirements for an anchor.
There are several types of anchors available, including dead weight, drag embedment, pile anchor, and plate anchor. It is noteworthy to add that the use of hawser is dependent on the size of the vessel to be anchored to the buoy, as hawser systems can use one or two ropes, as depicted in Figures 1-3.
A mooring system is made up of many materials and components that are organized in a specific way, as shown in Figure 4. The other SPM components are as follows: • The access to the buoy deck is provided by a boat landing; • The buoy is protected by fenders; • The material handling equipment includes lifting and handling equipment; • Maritime visibility aids and a fog siren are used to keep moving vessels alert and attentive; • The navigation aids or other equipment are powered by the electrical subsystem; • The sources of power systems are batteries and solar systems. While the batteries are replenished on a regular basis, the solar power systems employ sun-sourced renewable energy and maintain the charge in the battery packs, for electrical power; • A hydraulic system can be added for remote operation with PLEM valves, if needed.
CALM buoy bridle, anchor, and anchor chain or mooring chain. Harnois [160,161] provided a comparison of various mooring line materials. The inertia, elastic stiffness, and damping of a mooring line are affected by the material used. An anchor's purpose is to secure a mooring line to a fixed place on the bottom. The ability to resist high, horizontal, and in some cases vertical loads in a specific seabed type (soft to hard), cost-effectiveness, and ease of installation are the major requirements for an anchor. There are several types of anchors available, including dead weight, drag embedment, pile anchor, and plate anchor. It is noteworthy to add that the use of hawser is dependent on the size of the vessel to be anchored to the buoy, as hawser systems can use one or two ropes, as depicted in Figures 1-3.
A mooring system is made up of many materials and components that are organized in a specific way, as shown in Figure 4. The other SPM components are as follows: • The access to the buoy deck is provided by a boat landing; • The buoy is protected by fenders; • The material handling equipment includes lifting and handling equipment; • Maritime visibility aids and a fog siren are used to keep moving vessels alert and attentive; • The navigation aids or other equipment are powered by the electrical subsystem; • The sources of power systems are batteries and solar systems. While the batteries are replenished on a regular basis, the solar power systems employ sun-sourced renewable energy and maintain the charge in the battery packs, for electrical power; • A hydraulic system can be added for remote operation with PLEM valves, if needed.

Components of CALM Buoy System
The Catenary Anchor Leg Mooring (CALM) buoy system has a buoy with a pivot, called the turntable. This rotates around the vertical axis of the pivot, as the tanker is

Components of CALM Buoy System
The Catenary Anchor Leg Mooring (CALM) buoy system has a buoy with a pivot, called the turntable. This rotates around the vertical axis of the pivot, as the tanker is moored to it. The floating hose is also connected to the turntable, at an angle through the hose manifold. The elastically moored buoy of radius a is acted upon by a wave train of irregular waves and wave height H progressing in x-direction, as illustrated in Figure 1. The turntable on the mooring buoy can spin around its vertical axis. The tanker is moored to the turntable and is connected to the floating hose strings that are also attached to the turntable. Due to the forces imposed by the currents and waves, the entire system can freely rotate, which is termed weather-vaning. Figures 1 and 3 show Catenary Anchor Leg Mooring (CALM) buoy systems. Basically, there are three CALM system mooring components, namely, the anchors, the chain anchors and the chain stoppers. The anchors are used to hold things together, including the piles or gravity anchors for connecting the seabed with the mooring chain. The most common chain anchors are systems with either six or eight anchor chains. The third component is the stoppers for chains, which are for connecting the buoy with the mooring chains. The anchor chains help to keep the buoy in place. The fluid is transferred to the submarine hose strings via a swivel, which links to the undersea pipeline via the pipeline end manifold.

Different Mooring Configuration
There are other types of offshore mooring systems, aside SPM, as summarized in Table 1. Based on the mathematical modeling, HMs and MMs, considering SPMs for bonded marine hoses, have been developed over 45 years based on earlier works on point moorings and simple floating buoys. The application of offshore hoses has also led to advances in different mooring systems used in fluid transfer, as seen in Figure 5.   Moorings are also applied on shipping vessels for other oil field operations like CO 2 oil recovery [162] as seen in Table 1. Several studies assessed mooring statics and dynamics for CALM buoy, as well as with attached hoses, which were considered as a single point mooring (SPM) terminal [163][164][165][166][167]. The design of each hose-mooring system considers different loadings, predictive motion responses with structural statics/dynamics [168][169][170][171][172][173][174][175][176][177][178], and governing theories on the hydrodynamics of floating structures [179][180][181][182][183][184][185][186][187][188]. In addition, the design of FOS is based on different industry standards [189][190][191][192][193]. The application of a mooring configuration is based on the application requirement, the type of (un)loading operation, and the environmental conditions. Some of these mooring applications require floating, catenary, and reeling hoses, while others require submarine hoses.

Review on Physical Models on Hoses and SPMs
The selection of hose systems for single point mooring (SPM) systems has been described by Ziccardi and Robbins [194]. Setting up buoys in low-tide areas was discussed as the authors also wanted to stimulate more hose and flexible rubber pipeline designs and applications. They studied the SPM deployments at Tokyo Bay's Hakozaki and Koshiba terminals. They also included a timeline of hose design and trends. They claimed that the basic designs of under-buoy hoses and floating hoses are comparable. The strong crush resistance of sub-surface hoses, on the other hand, was shown to be dependent on the water depth. This was accomplished by increasing either the wire's area or the diameter of the helical wire, or both. The rated operating pressure was found to be 5 to 6 times the design burst pressure. They highlighted the abrasion and abuse that the floating hoses attached to the tanker from the buoy were subjected to. They came to the conclusion that developing flexible rubber lines that could sustain high operating pressures and external crush, particularly in severe environments, was critical. The hose system for an SPM terminal was also reliant on both the operational and environmental conditions, according to the report. Physical tests are also used to develop environmental wave spectra, such as the Joint North Sea Wave Project (JONSWAP) wave spectrum and regular wave types like Airy waves [195][196][197][198][199]. Typical recent numerical model of CALM buoy model conducted in Orcina's Orcaflex by the authors can be seen in Figure 6.
Images for jmse-1295533 -Proof Correction sent by Amaechi, Chiemela Victor, 29/10/2021.   The operational requirements, such as the system's working pressure, necessitated the transportation of well-specified products with an adequately defined nature. They listed several factors that must be considered when determining the length of hose strings in hose designs, including mean water depth, tide depth (low/high), maximum wave height, buoy position relative to pipeline header, maximum mooring distance, rated working pressure, desired throughput, and product(s) to be transported. They advised that the ultimate design of the under-buoy system should ensure that the hose does not come into direct contact with the seabed of the moored ship under high tide conditions. They looked at the primary design criteria for SPM hoses, underbuoy system, floating hose systems, float sinks, hose designs, and hose diameters, and encouraged greater research from hose manufacturers, using the two case studies that were employed by the US military on unloading from SPM tankers.
Earlier investigations on marine hoses depended on some lengthy calculations and experiments. Brady et al. [200] with the help of Shell B.P Petroleum Company of Nigeria Limited, built a test apparatus that was connected to 60.96 cm (24 inches) hoses attached to a CALM buoy off the coast of Nigeria. A Medilog 4-24 small four-channel cassette recorder, a 4-366 pressure transducer, and a Beaulieu S.P 16 mm Cine camera with a lens width of 5.9 mm were also included in the setup. To measure the strains on the hose of a monobuoy, a strain-gauge measuring spool was installed between the buoy manifold and the first-off buoy hose. They claimed that the hoses closest to the buoy have a lower life expectancy because they carry the majority of the hose stresses. Correlation of the measured loads was achievable using the statistical method described for calculating the 60 s recordings and visual records of the sea conditions. However, this was limited due to a lack of environmental data. Rather than the trial-and-error method employed previously, this technology enabled the investigation of the forces on buoy hoses. They came to the conclusion that the hose problem was primarily caused by fatigue rather than high loads. As a result, increasing the hoses' strength will improve their performance. SPM terminals were subjected to model testing by Pinkster and Remery [201]. The test results were also used to describe SPM terminal features and hose phenomena found in CALM and SALM mooring systems. They also stressed the need for selecting the appropriate scale for model tests. They cited water depth, the accuracy of the results, and the capacity to generate the needed wave height and period at a certain scale in the basin as critical variables. Water depth, current, wave generators, and wind are some of the variables that can be modified to affect environmental conditions. They also went over the model testing technique, measurements, and analysis in detail. They concluded that nonlinearities were exploited in the construction of the equation of motion, which was then integrated in small time increments step by step. Additionally, based on uncertainties in the prototype's drag coefficients, the inaccuracies in the estimates of the findings obtained from the model tests due to scale effects should be applied without modification. There was additional discussion of the Pierson-Moskowitz spectrum, wind forces, current forces, first-order wave forces, second-order wave forces, and drag wave drifting equations. However, the methodologies for calculating these forces were not sufficiently developed for design consistency. An industry collaboration with academia was conducted on the feasibility of using geodesic IGW designs for offloading hoses, as reported in Nooij [202]. Another important study that was carried out on the load response of offshore hoses by Lassen T. et al. [203] involved finite element models and full-scale testing for a 20inches-bonded hose with steel end fittings. The study presents limits based on API 17K [204] criteria for the extreme load capacity assessments. The study also included a methodology for predicting the fatigue life of bonded loading hoses' response to applied bending, tension and pressure using a catenary configuration, with reeling loadings repeated and significantly tensioned. The study emphasized the fatigue life prediction methods, as well as the load impacts on the hose during reeling operations, for both rubber and steel parts.

FPSOs for Marine Hose Operations
There are different types of FPSOs, as summarized (with details on their characteristics) in Table 2, that are used in SPMs for transfer, loading and offloading operations. The turret systems are the most common because of their freedom of movement, ease of anchorage, and accessibility during mooring and deployment. A typical turret FPSO in catenary mooring is shown in Figure 7a, and an Offloading FPSO attached to an SPM's CALM buoy is shown in Figure 7b. A variety of numerical models on other mooring systems can be seen in the literature and existing industry projects on marine hose, as earlier discussed. Suitable for large riser numbers with capability of additional tie-ins.

Station Keeping Performance
Lower number of anchor legs, offset is minimized.
Large number of anchor legs, offset is variable.

Vessel Motions
Motions are reduced as the vessel orients itself into the most suitable environmental direction.
Varies from small to large depending upon the relative direction of vessel and environment.
Riser Connection Turret provides the connection point for the risers.
Risers hang from the porch on the port/starboard side of FPSO.
Offloading Performance Better as the FPSO is aligned with the mean environment.
Depends on vessel/environment orientation.   Large number of anchor legs, offset is variable.

Vessel Motions
Motions are reduced as the vessel orients itself into the most suitable envi-Varies from small to large depending upon the relative direction of vessel

Model Methodologies and Software Tools
The designs of hose models are considered using some design parameters, as illustrated in Figure 8. Details of these parameters for the mathematical model are given in Sections 4 and 5. However, these mathematical models are also designed numerically using both commercial and in-house software packages. They are also verified using experimental models. The availability of model tests and numerical simulation tools determines the design of marine hoses, mooring lines, marine risers and other floating structures for deep and ultra-deep waters. The latter has been the subject of various research investigations in the past, leading to the classification of analysis approaches as uncoupled, completely coupled, or hybrid, as in the following sections.

Uncoupled Methodology
This was the traditional method for analyzing mooring lines and risers, which involved numerical tools based on uncoupled formulations. In general, this approach starts with motion analyses, where mooring lines and risers are represented by simplified models (e.g., scalar coefficients), and vessel motions are calculated in terms of static offsets, as well as wave frequency (WF) and low frequency (LF) components. In a second stage, those offsets and motions are used as input for structural studies of the lines modelled by finite element (FE) meshes. Typically, such studies are carried out for each individual riser. The fundamental benefit of the uncoupled techniques is that each simulation has a cheap processing cost. However, they have significant, well-known flaws that could result in significant inaccuracies in deep-water circumstances and Floating Production systems (FPSs) with a large number of risers. Such flaws were identified by [205][206][207][208][209], addressing the following topics: (1) The uncoupled approach does not normally account for mean current loads on moorings and risers, but in those scenarios, the interactive effects of current forces on the Submarine elements and the mean offset and LF motions of the floater are significant; (2) A significant damping effect of the moorings and risers on LF motions must be included in a simplified way, usually as linear damping forces. Because multiple characteristics are involved, creating simplified models of this event is difficult. It should be noted that this classification of "uncoupled techniques" could include a variety of analysis methodologies.
Based on the classification, the fundamental or "classic" uncoupled technique can be described in two ways: It is vital to distinguish between these two features because, even when only uncoupled programs are available, some amount of integration for both hose risers and moorings is achievable using variations on the "traditional" uncoupled methodology. These enhancements, which have been applied for over 25-35 years, include the use of

Uncoupled Methodology
This was the traditional method for analyzing mooring lines and risers, which involved numerical tools based on uncoupled formulations. In general, this approach starts with motion analyses, where mooring lines and risers are represented by simplified models (e.g., scalar coefficients), and vessel motions are calculated in terms of static offsets, as well as wave frequency (WF) and low frequency (LF) components. In a second stage, those offsets and motions are used as input for structural studies of the lines modelled by finite element (FE) meshes. Typically, such studies are carried out for each individual riser. The fundamental benefit of the uncoupled techniques is that each simulation has a cheap processing cost. However, they have significant, well-known flaws that could result in significant inaccuracies in deep-water circumstances and Floating Production systems (FPSs) with a large number of risers. Such flaws were identified by [205][206][207][208][209], addressing the following topics: (1) The uncoupled approach does not normally account for mean current loads on moorings and risers, but in those scenarios, the interactive effects of current forces on the Submarine elements and the mean offset and LF motions of the floater are significant; (2) A significant damping effect of the moorings and risers on LF motions must be included in a simplified way, usually as linear damping forces. Because multiple characteristics are involved, creating simplified models of this event is difficult. It should be noted that this classification of "uncoupled techniques" could include a variety of analysis methodologies.
Based on the classification, the fundamental or "classic" uncoupled technique can be described in two ways: (a) Uncoupled formulations utilized in the numerical tools deployed in the analysis of the floating structures, in which the vessel's hydrodynamic response is unaffected by the lines' nonlinear dynamic behavior; (b) The analytical technique takes into account the low or negligible integration for both the hose risers and mooring systems.
It is vital to distinguish between these two features because, even when only uncoupled programs are available, some amount of integration for both hose risers and moorings is achievable using variations on the "traditional" uncoupled methodology. These enhancements, which have been applied for over 25-35 years, include the use of improved algorithms for determining the scalar coefficients used in the vessel's motion equations to describe the behavior of these lines. Even these refinements contain a number of limitations. They do not capture the non-linear dynamics and their interplay in both the lines and vessel completely. This includes the nonlinearity obtained in the response of the damping of the sections that control the vessel's low-frequency (LF) motions. When the water depth increases, this fluid-structure interaction (FSI) increases its impact and has additional relevant effects, thus it could yield unreliable results and erroneous outputs for these deep-water scenarios.

Coupled Methodology
Using coupled analytical tools is the most accurate approach for considering the interaction between the components of an FPS and properly predicting the individual responses of the vessel, mooring lines and hose risers. In general, such tools use a timedomain solution method and a rigorous representation of the lines using FE models, taking into account any nonlinearities in the system's dynamic response. In the literature [210][211][212], descriptions of different coupling formulations investigated for the fully coupled analysis of FPS (referred to as "weak coupling" WkC and "strong coupling" StC) have been presented. The WkC formulation focuses on the hull equations of motion: forces operating on the right-hand side of the hull equations perform the coupling between the hydrodynamic model of the hull and the hydrodynamic/structural model of the lines. The finite element mesh of all moorings and hose risers is integrated into a single set of equations in the StC formulation, and the hull is regarded as a "node" of this model. Only coupled analytical tools can simultaneously produce vessel motions and detailed structural responses of the mooring lines and risers, resulting in a "completely coupled technique." The mooring lines and hose risers should be represented with a finite element mesh fine enough to precisely establish all key features of the structural response for this purpose (not only their top tensions, but also, for instance, tensions and moments near the touch-down zone). However, this may result in high computational costs for fully coupled analysis tools; in any case, emphasizes the importance of performing fully coupled analyses for at least some critical cases, or for specific studies wherein small-scale tests are not possible due to the limited depth in wave tanks.

Hybrid Coupled Methodology
Although devising alternate solution methods and algorithms to minimize the computer costs of coupled analyses (e.g., optimized time-domain methods with adaptive time-step variation, or frequency-domain methods) is an important line of research, the goal of achieving a balance between accuracy and computational efficiency has also been sought by proposing "hybrid" analgesics. This classification can include a variety of procedures, such as "coupled motion analyses" and "semi-coupled analyses", which will be discussed next. Since these hybrid processes may be distinguished by the use of coupled models that introduce simplifications or approximations, they shall be referred to as "simple coupled analyses" from here on.

Coupled Motion Analysis
This is a two-step technique for analysis, as shown in Figure 8. The first phase (the "coupled motion analysis" itself) uses a coupled model with FE meshes that are coarser than those used in fully coupled analyses. The meshes are created to accurately represent the hose risers' overall contribution in terms of top tensions, mass, and damping (therefore considering the non-linear dynamic interaction between vessel and lines). As a result, their vessel motions are more precise than those obtained from uncoupled motion analyses. With lower processing costs, the accuracy is comparable to that of fully coupled analyses. The detailed structural response of the hose risers (e.g., tensions and moments near the touchdown zone) is obtained in the second step of the procedure, which consists of prescribing the motions at the tops of the uncoupled FE models of each individual riser, now modeled with more refined meshes to obtain these results. Reviews of marine risers provide a more extensive overview of the approximations and benefits of this technique [123][124][125][126][127].

Semi-Coupled (S-C) Motion Analysis
This is a three-step technique for analysis. In the first phase, a coupled model (similar to the one described in Section 3.3.1) is used to perform a nonlinear static analysis with the static components of the loadings, yielding the system's mean equilibrium position. Other static studies begin with the static equilibrium configuration in the second step, using modest values of prescribed displacements at the platform's center of gravity (CG) for each of the six rigid-body DoFs. An equivalent 6DoF global stiffness matrix is calculated from the acquired comparable forces, representing the global contribution of all mooring lines and hose risers. This comprises a linearized tangent stiffness matrix around the static equilibrium configuration. Finally, an uncoupled dynamic analysis of the platform is performed, in which only the 6DoF vessel equations of motion are solved, and the equivalent global matrix is added to the platform's hydrostatic matrix. From the perspective of the lines, this technique is quasi-static, yielding conclusions that are not as stringent as those obtained from a completely coupled analysis, but they are more accurate than those obtained from a standard uncoupled formulation (in which the lines are represented by scalar models or catenary equations). The equivalent stiffness matrix of the lines is calculated using a full FE model, allowing it to be evaluated for a wide range of complex line configurations, including effects that are typically overlooked by traditional quasi-static scalar formulations, such as the influence of current loads along the lines, the coupling of two or more lines, and bending moments when the lines bend. The inertial and damping contributions of the lines can be integrated to increase the precision of the platform motions supplied by this S-C method by computing equivalent global six-DOF matrix coefficients for its mass and damping. Numerical decay tests can be used to achieve these calculations, when the primary goal of analysis is to observe the motion response.

Software Packages
Based on the mathematical model of marine hoses, marine risers and mooring lines, there are different types of software packages that can be applied, as summarized in Tables 3 and 4. There are different types of riser analysis software that can be applied in both the riser-specific and general-purpose software packages currently available. Table 3 discusses the most widely used riser software and its availability for this research at Lancaster University (denoted as LancsUni). Table 4 is also discusses different software that are applicable for the design and analysis of various mooring systems. This is detailed in the body of the paper, with more discussion of their formulation, governing equations and bounding theories. This section presents different software packages and their applications. For instance, Orcaflex ( [213,214]) has been applied on different offshore applications, from pipelaying, reeling, loading hoses to subsea lifting operations.   -FLEXCOM [223] x Open-source: In-house:

Mathematical Model and Other Model Types
Different theories related to the hydrodynamics and statics of bonded marine hoses attached to a CALM buoy have been proposed, as presented in this section. Since these hoses are high-pressure, high-temperature (HPHT) structures, it is important to also discuss some of the mathematical background of the test methods for the hoses.

Mathematical Formulation of Hose Model
Based on the given problem, the mathematical formulation of the hose model was developed with some simplifications, based on the experimental findings of the author and some experiences reported in offshore engineering practice. These applications can be seen on moored WECs [231], steel catenary risers [232][233][234][235][236][237][238], flexible risers [239,240] and pipelaying operations [241][242][243][244][245]. The seabed considered is a horizontal and a rigid plane, and the material for the hose is isotropic and in an elastic state. The marine environment for this mathematical model is a stable environment. To simplify the model, the contact points of the pipeline to the pipeline end manifold (PLEM) are not considered; however, the tension at hose ends is considered. The frame of reference is the Cartesian coordinate system (CCS), with the origin shown in Figure 1. This is the locus at which the pipeline meets the seabed, called the touch-down point (TDP). The natural coordinate system for the model is established along the hose string. It is also noteworthy that the arc length of the hose string, s, determines the physical properties of the submarine hose string.
Considering the wave-structure interaction, the Laplace Equation is used to formulate the equation governing the motion. This is derived from the Continuity Equation for fluids, as shown in Equation (1).
∂u ∂x The motion of the system can be represented by Equation (2). This presents the Newtonian force, F, derived from the external load of the system as the sum of the inertia force of the system, the viscous damping load and the elastic force components (also called the stiffness load of the system).
For incompressible flow, as considered here, the Continuity Equation applies, wherein the components of the flow domain are denoted by u, v and w.
For an irrotational motion, all the vector components of the rotation are equal to zero. Thus, For simplicity, let us denote the vector components in the Cartesian coordinates as u = ∂φ ∂x , v = ∂φ ∂y , and w = ∂φ ∂z , where the scalar functions ϕ(x,y,z) are the relations used in Equation (4).
Introducing the function φ to the Continuity Equation gives the second-order linear differential equation, called the Laplace Transform, as given in Equation (5).
where ∇ is the Laplace grad operator, u is the velocity of the fluid and the motion of the fluid can be expressed as ∇u = 0 based on Laplace formation. Thus, in a plane coordinate system, the Laplace Equation is given by the relation: For the conventional cylindrical coordinate or polar coordinate (r, θ, z), the Laplace Equation is given as Equation (7),: For a buoy of radius a and height L, the bounding surface is r = a and the two flat parts are z = 0 and z = L. Thus, the boundary conditions can be depicted by the following for the given function, called the Velocity Potential, φ(r,θ):

Assumptions
The buoy system is considered to comprise the buoy, the mooring cables and the connected submarine hoses. For this study, the hawser lines and the floating hoses are not included, as the responses from other FPSO and transport vessels are not considered. The buoy is also considered as a single-system rigid body with 6DoFs, as shown in Figure 1. The following is assumed: The fluid is incompressible, irrotational, and bounded by the surface of the buoy, the rigid bottom and the free surface; 2.
The seabed is horizontal and on a rigid plane. For the diffraction analysis, the fluid motion is in a cylindrical coordinate system of form (r, θ, z); 3.
The submarine hose is considered as a beam undergoing pure bending;

4.
The internal and external forces will place longitudinal forces on the hose. However, the effects can be negligible at depths with small effects; 5.
The hose curvature is the inverse of the minimum bend radius (MBR), and the curvature calculation can be approximated using 1 r = ∂ 2 z ∂x 2 . The measurement of the bend radius of the hose is never less than the MBR; 6.
The influences of both the horizontal forces and the shear forces on the curvature are negligible, depending on the bending moment; 7.
Due to some nonlinearities within the hose geometry, there will be some nonlinearities in the motions of fluids within the hose; 8.
The hose is considered to possess a solidly rigid body with constant bending stiffness for all given cross-sections transverse to the axis of the hose. The hose also transports (or carries) fluid under high pressure, and the fluid can be oil or water; 9.
The hose can be made of different sections, flanges, reinforced ends, floated sections and unfloated sections, and can have different section radii. A uniform density of the hose is assumed for both the rubber and steel sections.

Boundary Condition Formulation
The boundary conditions for the system are formulated with some assumptions in Section 4.1.2, as discussed herein. Illustrations of one (1) floating hose connected between the floating buoy and the FPSO, and two (2) submarine hoses connected underneath the buoy in a Chinese lantern configuration (with environmental forces on the system, such as waves), are shown in Figure 1. The buoy is cylindrical in shape and the submarine hoses are connected to the seabed by a fixed connection to the pipeline end manifold (PLEM). The seabed is assumed to be a horizontal plane. For simplification in modeling this system, a free-floating buoy was utilized. As such, in the initial boundary condition formulation, the mooring lines are not considered. Potential theory is used in the formulation of the boundary conditions, as given herein. The motion of the system is based on Equation (1), which is the Laplace Transform derived from the Continuity Equation of fluids. The following boundary conditions (represented as a-f) are related to Equation (2): (a) Dynamic boundary conditions: (c) Free surface boundary conditions: z = η(x, y, t) represents the free surface [246]. The free surface boundary conditions are given in Equations (6) and (7), which are to be satisfied by both the wave elevation η and velocity potential φ.
(d) Body surface boundary conditions: (e) Seabed (or bottom) boundary conditions: For an impermeable seabed of depth h(x,y), carrying a floating buoy, the seabed boundary condition is given by Equation (13); (f) Radiation boundary conditions: With the assumption of infinity, the radiation boundary condition is given by Equation (9), where φ S denotes the scattered wave potential for this condition and φ I denotes the incident wave potential.
Thus, Equations (15) and (16) will satisfy the radiation potential at infinity;

Boundary Layer
The boundary layer for the flow around the hose line attached to a CALM buoy can be idealized using a beam with a uniform cross-section, acted upon by velocities U 1 and U 2 , as shown in Figure 9. At initial conditions, the hose is stationary and the flow around it is assumed to be stationary. To formulate the boundary value problem (BVP), the wave potential theory was considered with grouped set equations for the boundary conditions presented in Section 4.1.3, following Lighthill's approach ( [247,248]), which was based on the Laplace Equation. The coordinate systems are those used in Figure 1, where XYZ is the CCS with a z-axis facing upward, and X-Y coincides with the axis acting across the free surface that is unperturbed. The fluid domain is not limited in the horizontal position. The seabed is also rigid and horizontal, and Equation (13) is valid for the z-axis. The total velocity potential for the system is given by Equation (18). This includes the incoming wave on both the hoses and the CALM buoy, as considered in the wave mechanics of HM-MM systems.   We used these following definitions: where g is acceleration due to gravity, h is the water depth below a mean water level, φ is the total velocity potential, φ s is the scattered wave potential, φ I is the incident wave potential, η is the height of the free surface, a is the radius of the buoy, k is the wave number, T is the wave period, C is the wave celerity and λ is the wavelength.
The wavelength, λ, the wave number function (called the wave celerity, C) and the wave period, T, are set as in the equations C = λ/T and k = 2π/λ [249][250][251][252][253][254][255]. Slender bodies such as risers, hose strings, piles and mooring lines make second-order wave loading contributions to the flow domain. Discussions on the effect of the theory on large vertical cylinders are available in the literature [181][182][183]. In principle, these properties contribute to the behavior of the structure in water or an ocean body. Equations (15) and (16) are derived from the perturbation series expressions for higher-order corrections, with Taylor series expansion.
For submarine hoses, the effect of the wave's action on the buoy motion and the hoses is considered during the design of the manifold and the connections. This helps the hose string to withstand the minimal horizontal movements made by the buoy [256,257]. The initial displacement of the hose under static conditions is denoted by Y 0 , the displacement at equilibrium is denoted by Y e and the sloped angle the hoses make with the manifold of the buoy underneath is θ 0 . At equilibrium, the stresses experienced by the hoses of the buoy are not dependent on the stresses experienced at the PLEM connection. However, there must be a sufficient length of hose string from the buoy to the PLEM.
Thus, the boundary condition at the buoy is given by Equation (19): While the boundary condition at the hose end is given by Equation (20): Another condition that is considered in the BVP is the fluid domain, as in Equation (21): 4.1.5. Modeling the Submarine Hose For this system, there are nonlinearities in the buoy system. These include the points of connection on the buoy for both the riser hoses and the mooring lines, the different materials used in the hose sections, and the wave trains for the irregular waves acting on the system. There are also nonlinearities in the free surface under the boundary conditions, presented via Equations (22) and (23), which are a result of the wave form, the current and the propagating waves.
Consider the submarine hose string depicted in Figure 10, with hose radius r and length S. The different sections of the hose string are designed differently. The front sections of the hose, connected to the manifolds on the buoy and the PLEM, are called the first-off hoses and are designed higher than the rest of the hoses, depending on the environmental condition. This involves more reinforcements and no floats at the beginning section, or no floatation covers, which enable it to withstand heavy stresses. Thus, the complexity of the problem can be reduced with the following equations, which are used to describe the submarine hose radius, where r u is the radius of the unfloated hose section and r f is the radius of the floated hose section. Figure 9. Boundary layer developing around a hose beam. Spline is used in the static calculation of submarine hoses, because it handles floating lines subjected to buoyancy better than the catenary method, and also enables better shape parameterization. However, when in equilibrium, the Bezier spline is used to obtain the shape of the submarine hose as a Bezier curve. This involves some control points obtained by integration, and this determines the order of the spline. Newman and Lee [258] discussed the utility of B-splines for velocity potential approximations in wave-structure interactions and boundary element method (BEM) computation. An approximation based on B-splines in a given order k is applied to simplify the velocity potential of the problem. For the boundary where z = 0, the radiation condition is given in Equation (24). Spline calculations use the sequence of the points on the hose at a given axis, as shown in Figure 10, where P i (x ki ), k = 1, 2, 3, i = 1.

Governing Differential Equations
Let us consider a section of the offshore hose, using the assumptions given in Section 4.1.2. Considering the hose's motion along the direction that is transverse to mean water level (MWL), the equation of motion is given in Equation (25), as presented in the literature [259], where the hose loading, Q, depends on the hose's weight, w, alongside the radius of the hose, r, located on the sea, with a water depth h and a hose section bending stiffness of EI z .
A small hose string segment that is relatively short, as given in Figures 11-13, can be easily accessed via the governing differential equations [121,122]. In Figure 11, the resultant force, T 0 , is localized at point A on the arc length, s. The horizontal force, H 0 , is assessed at the CCS locus point, called the origin, O, while V 0 represents the vertical force. Due to the motion of the hose string varying with different durations, different angles will emerge between the axis of the hose string and the horizon, θ (i = 1,2,3,..n) . The times required for a full wave cycle is given by n. The angle between the resultant force's direction and the horizontal is given by θ 0 .     the force per unit length of the hose element with respect to the depth of the hose from the mean water level (MWL) or surface of the sea is taken to be relative to the existing structure and the global loadings ( [263][264][265][266][267]). Thus Equation (31)

Hose Bending and Lateral Deflection
Considering the end at point A as depicted in Figure 14, the distance is x = 0. As the submarine hose strings move due to waves and some vibration motions, an asymptotic relationship can be derived. This is expressed in Equation (33), which gives the movement at the end of the hose string.
where ℎ is the amplitude of the buoy motion, σ is the circular frequency of the wave, ω is the angular frequency, and k is the wave number. It is noteworthy to state that this equation is only accurate when there is steady harmonic motion relative to a straight and lateral configuration of the buoy, positioned backwards. However, for a simpler hose motion in phase, the perturbation from waves can be represented as: The submarine hose being modeled is considered as a beam under tension, as shown in Figure 12. Normally, two types of differential equations are applied in the analysis of this type of problem. These are the second-order and third-order differential equations; however, a third-order equation is more suitable at the start of the operation. Different approximations have been made from analytical models by some researchers, but these have mostly dealt with floating hoses. The more popular approximations carried out on hoses include the tanh approximation, elastic foundation approximation and buoyancy approximation. The methods applied in the integration of these equations for marine hose's static analysis include the Runge-Kutta integration method, the fixed mesh difference method and the generalized-α integration method [260][261][262].
Considering the assumptions in Section 4.1.2, the Navier-Stokes equation is applied for an incompressible fluid carrying out a nonrotational motion at sea depth, z, on a floating buoy of depth, d. The velocity potential may be expressed in terms of the fluid quantities, thus yielding: Considering diffraction theory, there will be no normal flux or normal velocity for impermeable cases, as given in Equation (27); the equation may be simplified and solved by considering the problem in two dimensions (2D), in terms of the velocity potential, denoted as φ(x,y), as follows: However, the force on the submarine hose element, F, can be deduced using Cartesian coordinates (x, y, z) or polar coordinates (r, θ, z). As shown in Figure 11, the resultant force is a function of the of the pressure of the fluid , the sea depth, z and the angle made by the hose element, θ. Similarly, in Figure 13, the total force for the hose segment will be a function of the pressure of the fluid, (P o and P i), and the sea depth, z. Thus, The length of the submarine hose is relative to the sea depth, z and the longitudinal depth profile of the hose, L as seen in Figure 13. A very small section of the submarine hose is considered as the hose element, as represented in Figure 13. The forces acting on it, include the inner pressure denoted as P i , the outer pressure denoted as P o , the moment denoted as M, the tension force denoted as T while ds is the small length of the hose element. Also, the supporting floating structure (the CALM buoy) with buoy depth, d is influenced by some dynamic states, which can be discretised. However, for a sea depth z, the force per unit length of the hose element with respect to the depth of the hose from the mean water level (MWL) or surface of the sea is taken to be relative to the existing structure and the global loadings ( [263][264][265][266][267]). Thus Equation (31) can be expressed as:

Hose Bending and Lateral Deflection
Considering the end at point A as depicted in Figure 14, the distance is x = 0. As the submarine hose strings move due to waves and some vibration motions, an asymptotic relationship can be derived. This is expressed in Equation (33), which gives the movement at the end of the hose string.
u(x = 0) = a h cos σt (33) where a h is the amplitude of the buoy motion, σ is the circular frequency of the wave, ω is the angular frequency, and k is the wave number. It is noteworthy to state that this equation is only accurate when there is steady harmonic motion relative to a straight and lateral configuration of the buoy, positioned backwards. However, for a simpler hose motion in phase, the perturbation from waves can be represented as: Considering the hose bending and deflection, a displaced beam element is illustrated in Figure 14. It shows the tension forces, bending moments and torsional moments for the displaced hose beam element. This was developed by considering a beam model element undergoing forces and moments with longitudinal displacement u and lateral displacement w, as depicted in Figure 15. The equations of motion in the x-direction and z-direction can be represented as in Equation (35). Considering the hose bending and deflection, a displaced beam element is illustrated in Figure 14. It shows the tension forces, bending moments and torsional moments for the displaced hose beam element. This was developed by considering a beam model element undergoing forces and moments with longitudinal displacement u and lateral displacement w, as depicted in Figure 15. The equations of motion in the x-direction and z-direction can be represented as in Equation (35).   In the x-direction, the equation of motion is: Diving both sides by dθ and simplifying gives: Diving both sides by dθ and simplifying gives: Taking moments about the edge of the elements in LHS: Diving both sides by dθ and simplifying gives: Limiting δS to 0 gives: Applying beam bending theory gives: Applying the hose modeling method by O'Donoghue [39][40][41][42], assuming that θ = ∂z ∂x , Rdθ = dx, and w = z, and incorporating these into Equations (41) and (42), the equation of motion in the longitudinal or x-direction gives rise to Equation (43), and that in the lateral or y-direction gives Equation (44), similar to Equation (25): Equation (44) gives the equation of motion for beam bending for the lateral deflection z(x,t) along the length of the beam section [268][269][270], while Equation (43) presents the longitudinal motion of the beam, with the point of force applied at only x = 0, which is the position used in the beam model for the free lateral vibration's inertial end condition.

Hydrodynamic Model
In this section, a simple treatise on hydrodynamic forces will be presented, applying the boundary element method (BEM), which is applied when carrying out the hydrodynamics study.

Hydrodynamic Forces
Hydrodynamics affect both the CALM buoy and the attached hoses. When the hose string experiences some hydrodynamic forces, the hose bends. If it is a floating hose, it will bend on the water surface, but if it is a submarine hose, it will bend inside the water body. The damping force is the part of the hydrodynamic force that is proportional to the hose velocity, whereas the added mass force is proportional to the acceleration. The Euler beam approximation of the damping force and added mass force components in the z-axis, can be used to model the dynamical state of the hose beam. The force components for both the damping and the added mass coefficient can be respectively represented as −Q ∂z ∂t and (−C a m ∂ 2 z ∂t 2 ), where the damping constant is denoted as Q while the coefficient of added mass is denoted as C a . By incorporating these terms into Equation (44), we obtain: Similar to the Euler beam approximation, the hydrodynamic damping and added mass force are seen as insignificant in the context of the inertial end condition.

Snaking Model of Hose
Following the snaking model given in Figure 15, the solution without damping can be considered by taking V = z(x,t) + iy(x,t); assuming that the Euler beam equation without damping is represented by V as K approaches 0 [274] gives and when the inertial end condition is satisfied by z(x,t), this gives Equation (49), thus: The solution is given in the form V = Ae −iωt when Equation (49) The terms of B and D tend towards zero in this mathematical formulation based on physical grounds, whereby B = D = 0, as they negate the terms for the waves in A and i, thus: where the first term of the RHS of Equation (53) is the traveling wave propagated away by the floating buoy, while the standing wave is represented by the second term, which becomes exponentially lower as x decreases.
By considering the inertial end condition, we obtain: An accurate depiction of the motion of the hose string without damping can be represented as Equation (57), with consideration of the uniqueness of the equation under real conditions. Thus: However, as the hose string's narrow end approaches x = 0, when x > x 0 , then x 0 > 2a s . On the other hand, the solution with damping can be considered by taking V = z(x,t) + iy(x,t); by assuming that the Euler beam equation with damping is represented by V as K approaches 0 under inertial end conditions ( [42]), yields: Based on the dynamical state of the hose, the form V = Ae −iωt in Equation (58)

Buoyancy Force
When the marine hose with hose radius r and length S is not submerged in water, it is not acted upon by any buoyancy force because the marine hose is out of the ocean (B f = 0). This considers the fully submerged hose string depicted in Figure 1, as is the case with submarine hoses, and the force of buoyancy per unit of hose length is denoted by B f , as shown in Figure 16 for the case of floating hoses. The buoyancy force does not depend on the sea depth when it is completely submerged, thus the buoyancy force is the maximum, B f (max) . However, if it is partially submerged, as is the case with floating hoses, the buoyancy force becomes equal to the weight of the water the hose displaces, and this also depends on the depth below the MWL of the axis of the hose. Thus, the force of buoyancy when the submarine hose is out of water is given by Equation (60), while the buoyancy force for the fully submerged submarine hose is given in Equation (61):

Hose Material Models
The mathematical design of hoses cannot be complete without presenting some hose design models, while discussing material models for the behavior of hoses. Chesterton [6] and Amaechi [19,275] presented a full model of reeling hose in finite element modelling with unique multi-layered scheme, for the hose body and the reinforcement helix, as typified respectively in Figure 17a,b. However, Gao et al. [276] presented a simplified hose model for a 3D hose, with steel helix wire, reinforcement layers and rubber, and addressed symmetry using a 3D quarter model, as shown in Figure 17c.

Hose Material Models
The mathematical design of hoses cannot be complete without presenting some hose design models, while discussing material models for the behavior of hoses. Chesterton [6] and Amaechi [19,275] presented a full model of reeling hose in finite element modelling with unique multi-layered scheme, for the hose body and the reinforcement helix, as typified respectively in Figure 17a,b. However, Gao et al. [276] presented a simplified hose model for a 3D hose, with steel helix wire, reinforcement layers and rubber, and addressed symmetry using a 3D quarter model, as shown in Figure 17c. This helped to save computational time in carrying out the analysis and solving the equations behind the model. However, the later (simple model could not consider reeled sections which the former (full scale model) covered. There are other interesting numerical hose models recently been published on the effect of internal pressure on bonded marine hoses, using methods similar to those of the earlier investigation on multi-layered tubular structures, bonded marine composite hoses and composite risers [277][278][279][280].
In another model by Gonzalez et al. [281], a new approach for modeling the rebars was proposed for marine hoses using polymeric reinforcements. As shown in Figure 18, rebar layers are used to represent polymeric reinforcements embedded in an elastomeric material. Each lamina of fiber reinforcement, referred to as a rebar layer, is suitably positioned within the elastomeric matrix, which is depicted by a continuum finite element similar to the one discussed previously. The relative height of the lamina in relation to the height of the host continuum element totally determines the positioning. This parameter also acts as a constraint between the rebar layer elements and the host element, preventing changes in the relative height value, even after nonlinear analysis deformation. The lay angle, which is given with respect to a chosen global direction, the cross-sectional area of one fiber cable, and the distance between two successive fiber cords are all aspects of the reinforcement that must be defined. This helped to save computational time in carrying out the analysis and solving the equations behind the model. However, the later (simple model could not consider reeled sections which the former (full scale model) covered. There are other interesting numerical hose models recently been published on the effect of internal pressure on bonded marine hoses, using methods similar to those of the earlier investigation on multi-layered tubular structures, bonded marine composite hoses and composite risers [277][278][279][280].
In another model by Gonzalez et al. [281], a new approach for modeling the rebars was proposed for marine hoses using polymeric reinforcements. As shown in Figure 18, rebar layers are used to represent polymeric reinforcements embedded in an elastomeric material. Each lamina of fiber reinforcement, referred to as a rebar layer, is suitably positioned within the elastomeric matrix, which is depicted by a continuum finite element similar to the one discussed previously. The relative height of the lamina in relation to the height of the host continuum element totally determines the positioning. This parameter also acts as a constraint between the rebar layer elements and the host element, preventing changes in the relative height value, even after nonlinear analysis deformation. The lay angle, which is given with respect to a chosen global direction, the cross-sectional area of one fiber cable, and the distance between two successive fiber cords are all aspects of the reinforcement that must be defined. An important feature of hose models is the layers, which are usually made of rubberized or elastomeric materials, while the reinforcements are usually made of steel materials. The end-fitting of the marine hose, as shown in Figure 19, is usually bonded to its layers. As seen in Figure 19a,b, the end fitting has a flange, retention ring or locking ring, nipple, holding plies or textile plies, lining, main plies, helical wires and binding wires or reinforcement wires. The cross-section of the marine hose depicted in Figure 19c shows a hose material with the following layers: the liner, cord-1, bend stiffener, cord-2, cover, cord-3, liner, cord-4, and liner. An important aspect of material assessment is presented in Table 5, regarding the use of a matrix with similar rubber properties for marine hoses, offshore fenders and offshore safety barriers. Rubber models are considered based on different elastomeric attributes when designing, modelling and manufacturing hoses. This is important because different grades of elastomers have different behaviours to certain chemicals, high temperature, corrosive fluids and sea water. Milad et al. [282] and Aboshio et al. [283] analytically investigated the hyperelastic material behavior of a PVC/nitrile elastomer with a woven continuous nylon reinforcement composite sheet based on experimental findings under loading cases of uniaxial extension and pure shear. This was achieved via wide strip tension testing using a novel advanced non-contact optical strain measurement technique, carried out on an Imetrum system, and validated using ABAQUS hyperelastic material models for the curve fittings.
Based on the fiber design, Nooij [202] presented a geodesic model for modeling offloading hoses by applying Integral Geodesic Winding (IGW) technology and Netting theory. In the model, the theory had two limiting assumptions: that the stiffness matrix contribution is negligible and that the wall thickness remains small. This implies that the composite structure will be classified as a thin-walled structure with a thin revolving shell, under external load (A) and internal pressure (Pr), as illustrated in Figure 20. However, when the laminate's in-plane shear stress equals zero, this design can optimally utilize the maximal strength of its fiber bundle. The study was able to use composites in winding marine hoses, which have better flexibility than helical paths and higher fatigue resistance. However, the trade-off with high pressure (burst load) reported on the IGW hoses [202] was improved by composites in the geodesic design, which helped it to maintain a level position (balance), stable position (stability), and improve resistance to high motion response (vibration). An important feature of hose models is the layers, which are usually made of rubberized or elastomeric materials, while the reinforcements are usually made of steel materials. The end-fitting of the marine hose, as shown in Figure 19, is usually bonded to its layers. As seen in Figure 19a,b, the end fitting has a flange, retention ring or locking ring, nipple, holding plies or textile plies, lining, main plies, helical wires and binding wires or reinforcement wires. The cross-section of the marine hose depicted in Figure 19c shows a hose material with the following layers: the liner, cord-1, bend stiffener, cord-2, cover, cord-3, liner, cord-4, and liner. An important aspect of material assessment is presented in Table 5, regarding the use of a matrix with similar rubber properties for marine hoses, offshore fenders and offshore safety barriers. Rubber models are considered based on different elastomeric attributes when designing, modelling and manufacturing hoses. This is important because different grades of elastomers have different behaviours to certain chemicals, high temperature, corrosive fluids and sea water. Milad et al. [282] and Aboshio et al. [283] analytically investigated the hyperelastic material behavior of a PVC/nitrile elastomer with a woven continuous nylon reinforcement composite sheet based on experimental findings under loading cases of uniaxial extension and pure shear. This was achieved via wide strip tension testing using a novel advanced non-contact optical strain measurement technique, carried out on an Imetrum system, and validated using ABAQUS hyperelastic material models for the curve fittings.
Based on the fiber design, Nooij [202] presented a geodesic model for modeling offloading hoses by applying Integral Geodesic Winding (IGW) technology and Netting theory. In the model, the theory had two limiting assumptions: that the stiffness matrix contribution is negligible and that the wall thickness remains small. This implies that the composite structure will be classified as a thin-walled structure with a thin revolving shell, under external load (A) and internal pressure (Pr), as illustrated in Figure 20. However, when the laminate's in-plane shear stress equals zero, this design can optimally utilize the maximal strength of its fiber bundle. The study was able to use composites in winding marine hoses, which have better flexibility than helical paths and higher fatigue resistance. However, the trade-off with high pressure (burst load) reported on the IGW hoses [202] was improved by composites in the geodesic design, which helped it to maintain a level position (balance), stable position (stability), and improve resistance to high motion response (vibration).  Figure 19(a-b) depicts two end fittings, showing the flange, retention ring or locking ring, nipple, holding plies or textile plies, lining, main plies, helical wires and binding wires or reinforcement wires. The cross-section of the marine hose depicted in Figure 19(c) shows a hose material with the following layers: the liner, cord-1, bend stiffener, cord-2, cover, cord-3, liner, cord-4, and liner.

Elastomers General Attributes
Natural rubber Excellent physical properties, high elasticity, flexibility, very good abrasion, limited resistance to

Hose Stability Models
Hoses, like submarine pipelines, can be subjected to lateral forces; thus, lateral stability, motion stability, and material stability models exist, as presented in this section.

Hose Coordinate Systems
In modeling marine hose systems, three coordinate systems or frames of reference are considered, as illustrated in Figure 21. They are as follows: (a) The global frame of reference (x, y, z); (b) The CALM buoy frame of reference (x w , y w , z w ); (c) The curvilinear distance along the hose line, s, also used for moorings.

Hose Stability Models
Hoses, like submarine pipelines, can be subjected to lateral forces; thus, lateral stability, motion stability, and material stability models exist, as presented in this section.

Hose Coordinate Systems
In modeling marine hose systems, three coordinate systems or frames of reference are considered, as illustrated in Figure 21. They are as follows: (a) The global frame of reference (x, y, z); (b) The CALM buoy frame of reference (xw, yw, zw); (c) The curvilinear distance along the hose line, s, also used for moorings. Both the global and CALM buoy frames of reference are three-dimensional, defining two horizontal and one vertical Cartesian coordinate. The CALM buoy frame of reference's origin, Ow, is specified at the CALM buoy's center of mass, while the global frame of reference's origin, O, can be defined at any appropriate fixed position. The curvilinear distance is a one-dimensional coordinate measured along the hose line (also used for moorings), with s = 0 corresponding to the bottom end linked to the PLEM for hoses and anchors for moorings. Additionally, s = L, corresponding to the top end connected to the CALM buoy (L being the length of the hose line, also used for moorings).
In Figure 21, some angles and vectors are denoted, defined as follows: Φ is a vector containing three angles-Φ1 corresponds to the angle between x and xw due to the CALM buoy pitch displacement; Φ2 corresponds to the angle between y and yw due to the CALM buoy roll displacement; Φ3 corresponds to the angle between z and zw, due to the CALM buoy yaw displacement.
θw is the angle between the mooring line's top end and the CALM buoy at the connection point; θt is the angle between the mooring line and the ground at the touch-down point on the sea floor; Both the global and CALM buoy frames of reference are three-dimensional, defining two horizontal and one vertical Cartesian coordinate. The CALM buoy frame of reference's origin, O w , is specified at the CALM buoy's center of mass, while the global frame of reference's origin, O, can be defined at any appropriate fixed position. The curvilinear distance is a one-dimensional coordinate measured along the hose line (also used for moorings), with s = 0 corresponding to the bottom end linked to the PLEM for hoses and anchors for moorings. Additionally, s = L, corresponding to the top end connected to the CALM buoy (L being the length of the hose line, also used for moorings).
In Figure 21, some angles and vectors are denoted, defined as follows: Φ is a vector containing three angles-Φ 1 corresponds to the angle between x and xw due to the CALM buoy pitch displacement; Φ 2 corresponds to the angle between y and yw due to the CALM buoy roll displacement; Φ 3 corresponds to the angle between z and zw, due to the CALM buoy yaw displacement.

Environmental Conditions
Waves, currents, wind, and variations in water depth owing to the tide are the environmental inputs that drive the mooring line dynamics, as shown in Figure 1. Due to the action of the waves and currents, external forces are directly imparted onto the mooring lines from the fluid. Furthermore, because to the CALM buoy's motion, driven by waves, currents and winds, which create some variations in the water depth, a force is applied to the end of the hose line and the mooring line. This also generates some top tension at the hose's top connection to the buoy manifold. As a result, most hose lines have heavily reinforced ends at the top end to control internal pressure loads, tension forces and wave exciting forces. The exciting forces are a combination of three excitation modes based on environmental inputs, are as follows: Firstly, a mean force, or static loading, is created by the steady current, mean wind, and mean wave drift forces. Secondly, high-frequency (HF) forces are considered for oceans with higher wave characteristics in typical range of 0-0.35Hz, while low-frequency (LF) forces are induced by slowly varying wave drift forces, unsteady wind forces, and slowly varying tidal forces in the general range of 0-0.02 Hz. Thirdly, the wave frequency (WF) forces are induced by first-order wave forces in the general range of 0.03-0.3 Hz. This is applicable in larger floating structures like semisubmersibles. This is also considered in the mathematical computation of the Response Amplitude Operators (RAO) for the motion behaviour of the floating structure.

Hose Slenderness Ratio
In principle, the stability of marine hoses is a function of their slenderness ratio, as they are slender structures. Due to the high span-to-diameter ratio of a marine hose, especially when connected as a hose string, as well as its flexible characteristics, the bending effects are not negligible, unlike flexible riser systems or models [208,284,285].

Hose Floats and Buoyancy Module Models
A buoyancy connection on the hoses was designed with a float incorporated as part of the hose line. The design of float materials is performed in accordance with the industry requirements for OCIMF GMPHOM hoses [3,21,204]. The buoyancy of the submarine hose line is obtained by designing a series of floats arranged together, as depicted in Figure 22. With the application of the equivalence principle of the hydrodynamic loads per unit length and buoyancy load for the buoyancy section, as presented in [286], the equivalent float weight w e , equivalent float outer diameter D e , and equivalent hydrodynamic coefficients C de and C τe for the buoyancy section can be presented as in Equations (62)-(65), where w is the weight per unit length of riser, ρf is the material density of the buoyancy block, mfh is the mass of attached rigging of the buoyancy float block, and Cτn is the tangential drag coefficient acting on the cross-section of the buoyancy float block. Note that Equation (62) is purely for the case in [286], but not for other hoselines with floats. It defers from other baseline mathematical models for hoselines in Orcaflex documentation [213,214]. However, the equivalent normal and tangential added mass coefficients for the buoyancy section in this case by Ruan et al. [286] can also refer to the equivalent process of drag force coefficients [287], given as: the hose. Figure 19(a-b) depicts two end fittings, showing the flange, retention ring or locking ring, nipple, holding plies or textile plies, lining, main plies, helical wires and binding wires or reinforcement wires. The cross-section of the marine hose depicted in Figure 19(c) shows a hose material with the following layers: the liner, cord-1, bend stiffener, cord-2, cover, cord-3, liner, cord-4, and liner.

Governing Equations and Motion Characteristics
In this section, modeling methods and software tools for marine hoses are presented.

Static Analysis of Marine Hoses and Risers Subjected to Submerged Self-Weight
The mathematical model of marine hoses that are in a Lazy-S configuration is considered to be similar to that of marine risers under submerged self-weight. To design this, different methods are considered for the elastic line analysis. Of note is the finite difference method as reviewed herein, using earlier considerations for marine hoses on articulated platforms [287,288]. The finite difference method gives a more robust and efficient solution for nonlinear problems that have large angles within the geometrical aspect of the hose's bending. However, the hose must be designed with limits (2D or 4D, where D is diameter) to its diameter, depending on the type of the hose, as specified in the OCIMF standards [289][290][291].
For marine hoses, most are designed to be buoyant. Marine risers are not designed to include buoyancy, but buoyancy loads are considered among the loads in its design. As such, configurations such as Lazy-S are considerable for marine hoses and marine risers by ensuring that the elastic line has minimal flexural stiffness and hangs between both points of the mid-arch buoy (or floating buoy) and the floater (or FPSO or shuttle vessel). It is this segment that is exposed to most of the environmental loadings, which is considered in the following mathematical modeling.
To develop this, some assumptions have been made for the mathematical model under static and quasi-static loads, as follows: i.
The loads acting on the marine hose or marine riser are defined in the analysis; ii.
The marine hose is considered to be a hose string, similar to a marine riser in a single line; iii.
The marine hose has its own buoyancy, which is the only buoyancy considered, and it is assumed that there is no additional buoyancy in the system; iv.
The marine hose exists in a plane that lies in two dimensions, as produced by both static and quasi-static forces; v.
Horizontal tension is applied at the hose end attached to the floater (called the floater end), which predominantly controls the marine hose profile; vi.
It is assumed that the applied horizontal tension is supported by the floater's anchoring system; vii.
The two ends of the marine hose are hinged, whereby the end attached to the mid-arch buoy (called the mid-arch buoy end) has negligible moment as a result of the minimal flexural stiffness of the marine hose.

Motion Behaviour of Marine Hoses and Risers
The equation of motion for marine hoses has been defined in an earlier submarine hose model by Amaechi et al. [50], as given in Equation (66). With particular details regarding the degrees of freedom (DoFs) along the translational direction, the general equation of motion of the body in the horizontal plane can be represented as Equation (66) ..
However, for an articulated hose model, the equation of motion includes the measured movements and the descriptive load effects, as given in Equation (67).
where the damping forces and restoring forces may theoretically be nonlinear.

Static Equilibrium of Marine Hoses and Risers
In addition to the mathematical illustrations shown in Figure 23, some equations from the mathematical model of static equilibrium in marine hoses are presented in this section. The portion of the hose line that is suspended between the hose end and the floater end is depicted in Figure 1. Equation (68) depicts the relationship relating to the static equilibrium of the hose element, with a given length of hose, represented as ∂s.
where the damping force, F d can be expressed as:

Static Equilibrium of Marine Hoses and Risers
In addition to the mathematical illustrations shown in Figure 23, some equations from the mathematical model of static equilibrium in marine hoses are presented in this section. The portion of the hose line that is suspended between the hose end and the floater end is depicted in Figure 1. Equation (68) depicts the relationship relating to the static equilibrium of the hose element, with a given length of hose, represented as .  Rewriting Equation (68), we get: Differentiating Equation (72) with respect to x yields Via the substitution of ∂V ∂x and ∂T ∂x into Equation (73), using the terms fromEquations (67) and (70), For simplicity, it can be convenient to rewrite Equation (74) as: Equation (75) is a fourth-order differential equation of equilibrium with four boundary conditions in terms of w and x.
In terms of displacement, Jain [288] expressed the equation as: where q p and F d are respectively expressed as: For the bending moment, Equation (79) holds: where At a static equilibrium position, the bending moments at the buoy end and the floater end are assumed to equal zero. Thus, the equilibrium relationship in Equation (75)

Equation of Motion of Marine Hoses and Risers
Here, we consider a short hose line assumed to be in a straight line, similar to marine risers, whereby the deformation is relative to its orientation. Bennett and Metcalf [292] presented equations of motion for a structure along a 2D plane, with the equilibrium equations of Equations = − + ( , ) where ̂ is applied moment/length, V is transverse shear, u is axial displacement, w is lateral displacement, P is tensional force, x is the axial axis, T is the concentrated applied

Equation of Motion of Marine Hoses and Risers
Here, we consider a short hose line assumed to be in a straight line, similar to marine risers, whereby the deformation is relative to its orientation. Bennett and Metcalf [292] presented equations of motion for a structure along a 2D plane, with the equilibrium equations of Equations (83)- (85), constitutive Equations (86) and (87) and kinematic Equations (88)- (90); wheret is applied moment/length, V is transverse shear, u is axial displacement, w is lateral displacement, P is tensional force, x is the axial axis, T is the concentrated applied tension, k is the change in curvature, and r is rotational spring stiffness/length. As regards the constitutive equations, we have: With the application of rotation theory under conditions of large deflection and moderate conditions of ( δw δx ) 2 1, there are limitation that restrict ( δw δx ). Thus, the curvature becomes:

Analysis of Lazy-Wave Configuration
Permissible bending at the two bends in a "near" position and allowable tension at the upper end in a "far" position are both characteristics of an acceptable design. Vertical equilibrium and compatibility equations can be created, allowing unknown design parameters to be calculated. Three diverse design situations will be addressed with solutions. Initial higher and lower bend positions, as well as the effective weights of segments, are provided; segment lengths can be computed.

First Option
Initial upper and lower bend positions, as well as effective weights of segments, are provided; segment lengths can be computed.
When the desired bend positions and effective weights of all segments are determined, the required segment lengths may be calculated. There are six unknown length parameters in Figure 25, but we also have six equations; four are based on compatibility, and two are based on vertical equilibrium. One assumption here is that the initial condition lies below the MWL. Segment 1 is the lowest segment that partially rests on the bottom. In the initial condition, the length L 1 is vertical and the effective tension T B equals zero at the touch-down point (TDP). Segment 2 is the segment with uniform buoyancy. Its total length is L 2 . This segment will, in the initial condition, have a horizontal tangent at distance Z u from the upper end (or surface on the figure). This point is referred to as the upper bend. It is noteworthy to state that a horizontal tangent means that tension in the riser must be zero in the initial condition, since the horizontal force equals zero and the hose has no shear force. Segment 3 is the segment that ends at the floater. Segment length is L 3 . This segment will in the initial condition have a horizontal tangent at a distance Z L from upper end. This point is referred to as lower bend. The effective axial force at the upper end will be T T . For this analysis, w i is the effective (or submerged) weight per unit length of segment i (N/m), while D is the vertical distance from the upper end to the bottom.

Second Option
The length of the buoyant zone and upper segment, as well as the effective weight of the segments, are specified, and the bend position are computed.
Using some known variables, the unknown parameters include the tension of the upper end, denoted as TT, the length of the suspended part of segment 1, denoted as L1, the vertical position of the upper bend, denoted as Zu, and the vertical position of the lower bend, denoted as ZL.
In this case, there are three equilibrium equations, and one that defines compatibility:

Second Option
The length of the buoyant zone and upper segment, as well as the effective weight of the segments, are specified, and the bend position are computed.
Using some known variables, the unknown parameters include the tension of the upper end, denoted as T T , the length of the suspended part of segment 1, denoted as L 1 , the vertical position of the upper bend, denoted as Z u , and the vertical position of the lower bend, denoted as Z L .
In this case, there are three equilibrium equations, and one that defines compatibility: Equilibrium top-lower bend: T T = w 3 Z L Equilibrium bottom-upper bend: {(D − L 1 − Z U )w 2 + L 1 w 1 } = 0 Global equilibrium: L 1 w 1 + L 2 w 2 + L 3 w 3 = T T (105) Global compatibility: Standard manipulation of these equations gives the following solutions: L 1 = (A 3 + w 3 Z L )/w 1 (109)

Third Option
Upper and lower bend initial positions, effective weight of upper and lower segments, and upper segment length are all specified. The buoyancy segment's length and effective weight must be computed.
If the desired bend positions are known, vertical equilibrium at the initial position can be used to determine the length and weight of the buoyancy segment. L 1 , L 4 , L 5 , L 6 , and w 2 are unknown parameters. The following relationships can be used to find them:

CALM Catenary Configuration Equations
The governing equation used in the calculation of the statics for the mooring lines is the catenary equation. It is also applied in other applications, such as steel catenary risers (SCR) and cable structures. In the case of mooring lines, which suspend from the CALM buoy to the anchor on the seabed, the catenary shape is approximate, as shown in Figure 26. The catenary line is defined by the following equation (118), where w denotes weight per unit length and H denotes the tension in the horizontal component.

CALM Catenary Configuration Equations
The governing equation used in the calculation of the statics for the mooring lines is the catenary equation. It is also applied in other applications, such as steel catenary risers (SCR) and cable structures. In the case of mooring lines, which suspend from the CALM buoy to the anchor on the seabed, the catenary shape is approximate, as shown in Figure  26. The catenary line is defined by the following equation (118), where w denotes weight per unit length and H denotes the tension in the horizontal component.
= . (120) The catenary equations for CALM systems are presented in Equations (118)-(124), where h (or z) denotes the height above seabed, s denotes the arc length, x denotes the section length of the mooring cable, ws denotes the submerged weight, TH denotes the tension along the horizontal component, Tv denotes the tension along the vertical The catenary equations for CALM systems are presented in Equations (118)-(124), where h (or z) denotes the height above seabed, s denotes the arc length, x denotes the section length of the mooring cable, w s denotes the submerged weight, T H denotes the tension along the horizontal component, T v denotes the tension along the vertical component, T denotes the tension component along the plane, V denotes the body's volume and W denotes the body's weight. The schematic of the catenary mooring system of a CALM buoy is illustrated in Figure 27. The suspended part of the line, s [172] is given by: ponent, T denotes the tension component along the plane, V denotes the body's volume and W denotes the body's weight. The schematic of the catenary mooring system of a CALM buoy is illustrated in Figure 27. The suspended part of the line, s [172] is given by: = ℎ + − + (124) Figure 27. Schematic of catenary mooring system of CALM buoy. Figure 28 depicts the geometry of a lazy-wave riser as well as all of the required geometry parameters. If the horizontal force H, all weights, and segment lengths are known, and the bottom segment is considered to be long enough to maintain zero vertical force at the bottom end, the unknown parameters are: L1 is the lower suspended segment length, and z2 is the upper bend vertical position.  Figure 28 depicts the geometry of a lazy-wave riser as well as all of the required geometry parameters. If the horizontal force H, all weights, and segment lengths are known, and the bottom segment is considered to be long enough to maintain zero vertical force at the bottom end, the unknown parameters are: Figure 27. Schematic of catenary mooring system of CALM buoy. Figure 28 depicts the geometry of a lazy-wave riser as well as all of the required geometry parameters. If the horizontal force H, all weights, and segment lengths are known, and the bottom segment is considered to be long enough to maintain zero vertical force at the bottom end, the unknown parameters are: L1 is the lower suspended segment length, and z2 is the upper bend vertical position. Figure 28. Definition of lazy-wave geometry parameters. Figure 28. Definition of lazy-wave geometry parameters.

Catenary Analysis of Lazy-Wave Configuration
L 1 is the lower suspended segment length, and z 2 is the upper bend vertical position. The riser geometry can be simply determined using classical catenary theory if these parameters can be found. For the two unknown parameters, a direct solution is not available. As a result, an iterative technique must be devised. First, the relationship between L 1 and z 2 must be discovered, and then equations describing compatibility for all segments must be developed.
The link between segments 1 and 2 (bottom and buoyancy segments) is z 2 in the vertical position. The following relationship between z 1 and L 1 is derived from the catenary theory in Equation (125) for a given horizontal force: Consider segment 1 or the portion of segment 2 between the connection point and the higher bend to find tension in the riser at the connection point between segments 1 and 2. Since we assume there is no shear force in the riser, this is possible. As a result, we have Equation (126): Since the buoyancy of segment 2 will obviously be greater than its weight, w 2 in Equation (125) is negative. If we assume that z 2 is known, this problem can be solved in Equation (127) with respect to z 1 : We may find the relationship between L 1 and z 2 without any other unknowns by entering Equation (127) into Equation (125), resulting in Equation (128) and (129): 0.5 (129) Compatibility can be determined by combining all of the segments' vertical projections and comparing the result to the known vertical distance between the upper and lower ends, using the parameters as defined in Figure 28.
Thus, Equations (130) and (131) gives the vertical distances z 3 and z 4 : When the segment lengths, horizontal force, and riser weights are known, the angles at segment boundaries and higher ends can also be determined, as shown in Equations (132)- (134). Finally, Equation (135) may be used to determine the vertical position of upper bend z 2 . See the literature [293] for more information on this calculation.
With the above, the vertical position of the upper bend z 2 can be obtained, using:

Fundamental Approaches to the Motion of Marine Hoses and Risers
There are three fundamental approaches that can be used in the mathematical modeling of marine hoses. The dynamic responses of the marine hoses can be modeled as follows: (a) Time history analysis by the direct integration approach; (b) Time history analysis by the mode superposition approach; (c) Frequency domain analysis by the steady-state approach.
Firstly, the direct integration method considers time, and involves step-wise numerical integration. The equations of motion with respect to time have already been defined in earlier sections. They show that the computation can be taken step-by-step, and could be time-consuming. For the second, the mode superposition approach, there is a transposition of the equations of motion into the modal space. The modal shapes are considered in the deformation and buckling of the marine hose, by using the dynamic responses of the marine hose system. This approach is utilised in the FEM of flexible risers [284,285,294,295]. In the third approach, it is assumed that the time-dependent terms are harmonic, with the form shown in Equation (136), with respect to time. With the elimination of time in this equation, the transience of the hose response can be neglected to improve the convergence time of the steady-state solution by making it faster.
where Y is the lateral response of the marine hose, w is the angular frequency and t is the time.
In the consideration of hydrodynamic loads, it is assumed that the current, waves and hose motions occur in the same plane along the beam-column axis, with the same small angle conditions and linear elastic behavior.

Marine Hose and Riser Response Equations
The hose curvature is given by the inverse of the minimum bending radius (MBR), as: The limit of permission for the bending radius is subject to the stiffness, EI, which thus becomes: Hose bending is a key parameter that requires some sensitive investigation, due to strains in the hose material and other forces on the hose. Figure 29 shows the equivalent force system of marine riser or hose segment for internal fluid and external fluid flows. When external pressure pe is also present, a similar approach can be used. Both the lateral pressure effects are removed by adding the force systems acting on the pipe section and the internal fluid, then subtracting the force system operating on the displaced fluid. When external pressure pe is also present, it is approached in a similar fashion, as depicted in Figures 29 and 30. Conversely, the schematic in Figure 30 presents the force loads and pressures acting along a hose string segments with float collars, for the loads on the (a) floating hose, and (b) submarine hose.
Based on the hose content or marine riser content, the pressure field acting on the internal fluid column is closed and in equilibrium with the weight of the internal fluid. The lateral pressures acting on the pipe wall are equal and opposite to those acting on the internal fluid. Hence, by the superposition and addition of the two force systems, those lateral pressures are eliminated. However, the supposedly axial tension in the fluid column, denoted as −piAi, remains, where pi is the internal pressure and Ai is the internal cross-sectional area of the pipe [121,122]. This leads to equations for the effective tension Te and apparent weight wa of the equivalent system.
The force systems operating on the pipe section and the internal fluid are added together, and then the force systems are subtracted. In Figures 29 and 30, wt denotes the equivalent system's weight, we denotes the weight of the displaced fluid column, wa denotes the weight of the internal fluid column, and wi denotes the weight per unit length of the tube.
The effective tension, , denotes the axial tension at any point of the riser calculated by considering only the top tension and the apparent weight of the intervening riser segment [49,50]. The equations of effective tension Te are expressed as in Equation (139) and (140). When external pressure p e is also present, a similar approach can be used. Both the lateral pressure effects are removed by adding the force systems acting on the pipe section and the internal fluid, then subtracting the force system operating on the displaced fluid. When external pressure p e is also present, it is approached in a similar fashion, as depicted in Figures 29 and 30. Conversely, the schematic in Figure 30 presents the force loads and pressures acting along a hose string segments with float collars, for the loads on the (a) floating hose, and (b) submarine hose.
Based on the hose content or marine riser content, the pressure field acting on the internal fluid column is closed and in equilibrium with the weight of the internal fluid. The lateral pressures acting on the pipe wall are equal and opposite to those acting on the internal fluid. Hence, by the superposition and addition of the two force systems, those lateral pressures are eliminated. However, the supposedly axial tension in the fluid column, denoted as −p i A i , remains, where p i is the internal pressure and A i is the internal cross-sectional area of the pipe [121,122]. This leads to equations for the effective tension T e and apparent weight w a of the equivalent system.
The force systems operating on the pipe section and the internal fluid are added together, and then the force systems are subtracted. In Figures 29 and 30, w t denotes the equivalent system's weight, w e denotes the weight of the displaced fluid column, w a denotes the weight of the internal fluid column, and w i denotes the weight per unit length of the tube.
The effective tension, T e , denotes the axial tension at any point of the riser calculated by considering only the top tension and the apparent weight of the intervening riser segment [49,50]. The equations of effective tension T e are expressed as in Equation (139) and (140).  The apparent weight, w a , can then be represented by Equation (141): The resolution of forces along the axial direction of an element of length ds can be represented as: dT e ds = w a cos ϕ Considering the resolution of forces along the vertical plane of an element of length ds with small vertical angles yields: dT e ds = dT e dx = w a (143)

Concluding Remarks
A review on the mathematical modeling of bonded marine hoses for single point mooring (SPM) systems, with the application of Catenary Anchor Leg Mooring (CALM) buoys, has been conducted. Due to the need to design, analyses, and optimize bonded marine hoses, mathematical modeling for SPM mooring systems involves many methodologies from other offshore engineering domains. Various modeling strategies of increasing complexity and fidelity are available, with the most realistic models being computationally expensive in comparison to their simpler equivalents. These models are based on wave loads on the offshore structures. The sort of model chosen will be determined by the CALM buoy's operating principles, the mooring system employed, and the analytic application.
This review has achieved the following: • A mathematical modeling review of marine hoses, SPM moorings and CALM buoys; • Assessment of marine hoses used for CALM buoys and single point moorings; • An overview of single point moorings with the contribution to other marine applications; • Assessment of marine industry application of mooring models (MMs) and hose models (HMs); • The mathematical modeling of hose behavior, the effect of waves and hydrodynamics.
Studies reporting on hose mooring array analysis, CALM buoy hose system simulations, extreme loads, CALM buoy hose system control, the use of HMs and MMs in design, hose strength analysis, and hose configuration, are compiled separately in Sections 2 and 3. As regards mooring, fatigue and heavy load studies suggest the need for the highest quality MMs capable of capturing dynamic tensions and effects across the mooring line. The effect of the mooring system on CALM buoy dynamics, i.e., the mechanical impedance of the mooring system, is essential to CALM buoy modeling and control applications, which generally allow for the use of more basic MMs, but this is rarely detailed in research reports. For the real-time computation of CALM buoy control, as well as long-term simulation, it is beneficial to consider time periods, weather conditions and sea states as case studies by using simpler models in fully developed sea state simulations of HM-MM models of CALM buoy systems. Simple models are used in CALM buoy hose systems, hose mooring array layouts, and the optimization of HM-MM design, with model fidelity increasing as the design and analysis become more advanced. The mathematical model has also been reported as applied in field experiences and safety of marine hoses.
Models of the hydrodynamics of marine hoses are presented in Section 4.2. It is noteworthy that the whole range of phenomena that HMs and MMs can capture covers multiple time scales, including ultra-high-frequency effects such as low-frequency (LF) drift motions, low-frequency (LF) motion variations, wave frequency (WF) oscillations, dynamic tensions, snap loads, load disconnection effects and ultra-low-frequency motion responses. To capture the HF effects, the simulations must use very small time-steps, but large simulation durations are necessary to catch the LF effects. Furthermore, the mooring system and the CALM buoy hose models should be subjected to a variety of sea states and environmental variables, which should be investigated using different simulations. These include utilising FPSOs, floating buoys and the attached marine hoses. To achieve the optimal balance between computational requirements and the accuracy of outcomes, pragmatic decisions about the sort of mathematical models to utilize and which inputs to examine must be made. The review also covers challenges encountered in hose installation, connection and hang-off operations. The state-of-the-art, developments, and recent innovations in mooring applications for SURP (subsea umbilicals, risers and pipelines) have been presented. Finally, this study details the materials required for model implementations in hoses and moorings. Funding: The Department of Engineering, Lancaster University, UK is highly appreciated. In addition, the funding of Overseas Scholarships by Nigeria's NDDC (Niger Delta Development Commission) is also appreciated, as well as the support of the Standards Organisation of Nigeria (SON), F.C.T Abuja, Nigeria. We also acknowledge National Natural Science Foundation of China (NSFC) for supporting the Projects 51922064 and 51879143, including this study. The article processing charges (APC) for this article was funded by Author 1, with support from MDPI's JMSE.
Institutional Review Board Statement: Not applicable. It was purely research on marine hoses, and no questionnaire was used or prepared for ethical committee review board. Also, there was no human trails or clinical tests involved as it is on marine hose review, thus it is presently not categorised under studies for institutional review. It was purely a review on mathematical modelling on hoses, so the research on marine hoses was based on hose modelling, author's experience on hose mechanics, and not from questionnaires or samples.

Informed Consent Statement: Not applicable.
Data Availability Statement: The data supporting the reported results cannot be shared at this time, as they have been used in producing more publications on this research.

Acknowledgments:
The author acknowledges the technical support from my research team in Lancaster University Engineering Department and Lancaster University Library staff, especially during the COVID-19 pandemic. We are also grateful to Abiodun K. Oyetunji of Lancaster University, Lancaster Environmental Centre (LEC), Lancaster, UK for reviewing this manuscript. The permissions granted to use images by: Bluewater for image used in Figure 2, Offspring International (OIL) for adapted image in Figure 3, PennWell Books for adapted image in Figure 11 and Elsevier Publisher for image used in Figure 17c are also duly appreciated and well acknowledged.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.