1. Introduction
UAVs have been widely utilized in resource exploitation, anti-terrorism surveillance, civil aerial photography, and other sectors due to static hovering and stability. However, UAVs struggle with endurance and the wirelessly transmitted data they relay are readily hacked. To make the UAV platform operate for a long time and assure the effectiveness and safety of information transmission, it can be attached to an optoelectronic composite tether. Tethered UAVs are this kind of UAV.
Tethered UAVs have a wider range of potential applications and the benefits of lengthy hovering durations and plentiful power over regular UAVs. These benefits have led to the widespread deployment of tethered UAVs in areas such as marine emergency communication [
1], forest firefighting [
2], and bridge safety detection [
3].
The photoelectric composite cable serves as both the power source and the information conduit for the tethered UAVs. Due to its interconnected nature, the UAV platform’s movement will have an impact on the internal strain of the tethered cable. When the strain is too high, the tethered cable will be harmed, which will also cause the tethered UAVs to malfunction. As a result, it is not sufficient to focus just on the UAV platform’s dynamic reaction when it is impacted; The way UAV platform affects the tethered cable must be taken into account. In order to guarantee the long-term functioning of tethered UAVs, it is necessary to explore their dynamic properties in a challenging environment. For tethered UAVs to operate steadily, a large number of professionals have worked on the dynamic model and control principles. Wu neglected the cable’s local bending in favor of simplifying the tethered cable to the spring-mass-damping model while creating the control algorithm for the ship-borne tethered UAVs [
4]. Nicotra focused on creating a stabilizing control rule for a cable-connected aerial vehicle attached to a base station [
5]. Using the UAV as a constraint at one end of the cable, Muttin simplified the ship-borne tethered UAVs system and built a nonlinear dynamic model [
6]. Before optimizing UAV pathways to maintain tight tethers, Ya looked at the dynamic available wrench set and maximum maneuver acceleration [
7].
Tethered cables often consist of flexible ropes with tiny diameters, making them easier to bend and having short portions. The tethered cables will swing as a result of the external wind load and UAV motion, cycling between a taut and slack condition. Because the cable’s optical fiber has inadequate bending resistance in this instance, the signal transmission will be hampered. According to certain research, the snap load in the cable is 5 times or even more than 10 times the usual tension when the mooring cable structure is under relaxation tension [
8]. According to Basile Audoly’s research, a localized increase in curvature caused by the abrupt relaxation of the curvature at this end causes a burst of flexural waves, which permits a cascade failure process [
9]. Tjavaras investigated the dynamic behavior of a tethered near-surface buoy subject to wave excitation and discovered that, at a threshold wave amplitude, the mooring system exhibits first zero tension, then a snapping reaction, and the buoy shows chaotic motion [
10]. In addition, he discovered that the curvature may significantly alter the magnitude of the shock wave’s stress during the research of synthetic curved cables [
11]. With the use of the Hamiltonian principle, Perkins was able to construct the three-dimensional dynamic equation for the moving cable and discovered that the equilibrium configuration’s nonzero curvature is what links the cable tension and velocity [
12]. Two qualitatively different in-plane waveforms were discovered by Behbahani Nejad while researching elastic wave free propagation in the cable. One wave type mostly generates longitudinal deformations and propagates with a speed that is similar to that of an elastic rod, whereas the other wave type primarily induces transverse deflections and propagates with a speed that is similar to that of a taut string [
13]. Additionally, he studied the forced response of an elastic cable that was infinitely long and slightly curved while being exposed to focused harmonic stimulation using the Green’s function approach. According to research, transverse excitation can produce a sizable, long-distance tension wave [
14]. Lin studied the structural wave generated and propagated in the interaction between the seabed and mooring structure through the numerical solution and found that the seabed has a significant influence on the propagation of transverse waves in the cable [
15].
There is currently no study on the dynamic response of ship-borne tethered UAVs under taut–slack conditions. This is due to the fact that the ship-borne tethered UAVs are impacted by ship motion at one end and attached to the UAV platform at the other end, and the boundary conditions are dynamic. In this study, a set of wave dynamics equations for the ship-borne tethered UAVs are created, and a numerical analysis of the tension and velocity of ship-borne tethered UAVs under the conditions of the ship’s heaving is carried out by the characteristic line technique. After that, the trial confirms the analyses’ findings, which form the basis for the ship-borne tethered UAVs’ design.
The remainder of the paper is organized as follows. In
Section 2, the dynamic equations of ship-borne tethered UAVs are derived and then solved by the stress wave method. In
Section 3, the numerical simulation is carried out to analyze the dynamic behaviors under different conditions. In
Section 4, the model experiment is established to verify the former results. Finally,
Section 5 concludes the paper.
2. Theoretical Analysis
2.1. Governing Equations of the Ship-Borne Tethered UAVs
A theoretical model is derived that describes the nonlinear, three-dimensional motion of ship-borne tethered UAVs. The tethered cable is considered to be a homogeneous, linearly elastic, one-dimensional continuum with negligible torsional, bending, and shear rigidities and is subject to gravity, tension, and wind forces.
Figure 1 illustrates the tethered UAVs’ equilibrium configuration (dashed curve)
and the dynamic configuration (solid curve)
following a disturbance. The position vectors
and
describe the location of a material point on the cable centerline in the equilibrium and dynamic configurations, respectively. The three-dimensional dynamic response of the cable about equilibrium
is projected onto the equilibrium Serret–Frenet triad composed of the unit tangent
, normal
, and binormal
directions. Hence,
, where
denotes the equilibrium arc length coordinate and
denotes time.
The strain energy of the tethered UAVs in the dynamic configuration is:
where
is the strain energy in the equilibrium configuration,
is the equilibrium cable length,
is the equilibrium cable tension distribution,
is the cable cross-sectional area, and
is Young’s modulus. Here,
denotes the dynamic component of the Lagrangian strain measure for centerline extension.
where
is the equilibrium curvature.
The gravitational potential energy of the system in the dynamic configuration is:
where
is the angle between the
direction and horizontal direction.
The kinetic energy of the system in the final configuration is
where
is the material density of the mooring cable in the ship-borne tethered UAVs, and
is the absolute velocity of each point on the tethered cable, and its expression is:
The virtual work caused by wind load is considered as
where
,
are the components of external force on tethered UAVs along
, and
directions, respectively.
Substituting those expressions into Hamilton’s principle:
The equations of motion for the ship-borne tethered UAVs are derived:
It can be obtained that the in-plane equations of motion are coupled with each other by the equilibrium curvature , and the out-of-plane equation is coupled, so only the in-plane equations of motion are discussed.
The tethered cable is fixed to the ship at its lower end and is attached to the UAV platform at its upper end. The lower end of tethered UAVs is determined by ship motion. The force diagram of the UAV platform is shown in
Figure 2. The boundary condition of the upper end is derived as follows:
Where
is the body coordinate system in the plane of the UAV,
is the horizontal wind speed,
is the total lift of the UAV,
is the inclination between
and the horizontal direction,
is the included angle between
and the horizontal direction,
is the gravity of the UAV,
is the total tension of the mooring cable to the tethered UAVs, and
and
are the aerodynamic force received by the UAV in
and
directions, respectively. The expressions are:
where
is the component of
in
and
is the component of
in
, and the expressions are:
Decomposing the force on the UAV along
directions, the boundary condition of the upper end is:
2.2. Stress Wave Solution of Ship-Borne Tethered UAVs’ Motion Equations
The tangential, normal, and sub-normal dynamic equations of ship-borne tethered UAVs are obtained by the Hamilton principle as shown in Equations (17)–(19):
The dynamic strain, equilibrium tension, and equilibrium curvature are as follows:
where
A is the cross-section area of the cable,
ρ is the material density of the cable, and
P0 is the horizontal component of the equilibrium tension.
In the small-curvature hypothesis
, Equations (21) and (22) can be changed into:
By substituting Equations (23) and (24) into Equations (17)–(19), the following results are obtained:
It can be obtained from Equations (25)–(27) that the dynamic equation of the auxiliary normal tethered UAVs is the tension string equation, so we mainly focus on the tangential and normal dynamic equations in the plane.
The characteristic line is an optimized finite difference method. When the characteristic line equation and the compatibility relationship are determined, the solution of the partial differential equations can be reduced to the corresponding ordinary differential equations, which are much easier to solve. The characteristic line method can give clear physical images of wave propagation, and can directly deal with wave images with discontinuities. The characteristic line method is widely used in the analysis of the KP equation [
16] and explosion equations [
17].
According to the linear elastic assumption, the tangential strain and normal strain are [
11]:
In the case of neglecting the external force of air acting on the cable, substituting Equations (28) and (29) into Equations (25) and (26), the dynamic equation is obtained:
Because
are derivatives of
about
, and
are monotonic and continuous, we have a compatible equation:
Two algebraic equations are obtained from Equations (30)–(33) by the indefinite line method:
As there is more than one characteristic line, the solutions of Equations (34) and (35) are not unique, and the determinants are zero:
According to Equations (36) and (38), the velocities of the longitudinal wave and transverse wave in tethered UAVs are as follows:
Substituting Equations (40) and (41) into Equations (37) and (39), respectively, the compatibility relationship along the characteristic line is obtained:
It can be obtained from Equations (42) and (43) that the compatibility relationship along the characteristic line contains and , so the phase velocity of the wavefront is not constant. It describes the dispersion characteristics of a transverse wave and a longitudinal wave in ship-borne tethered UAVs. The coefficients of these terms are related to the equilibrium curvature . It can be seen that even if the cable has a linear elastic constitutive relation, the internal stress wave will propagate. However, the phenomenon of geometric dispersion caused by the geometric size factor still appears.
Generally speaking, the velocity of a longitudinal wave is greater than that of a transverse wave. For each point on the tethered cable, when the longitudinal wave arrives first and the transverse wave does not arrive, it should only show the characteristics of a longitudinal wave, and when the transverse wavefront arrives, the position should show the coupling characteristics of a transverse wave and longitudinal wave.
Considering the influence of the ship’s heave on the tethered UAVs, the lower end of the ship-borne tethered UAVs moves in a simple harmonic direction, so the lower boundary conditions of the stress wave are shown in Equations (44) and (45):
Among them, Ae is the excitation amplitude of the ship to the tethered UAVs, and ω is the excitation circle frequency of the ship to the tethered UAVs.
The top of the ship-borne tethered UAV keeps its balance by adjusting its lift at all times, so the upper boundary conditions are as follows:
It can be obtained from the compatibility equation that it is challenging to solve the analytical solution, which is solved by the characteristic line numerical method. The difference method is used to establish the discrete scheme of the compatible equation.
The difference scheme of boundary points and interior points is shown in
Figure 3.
3. Numerical Simulation
From the differential format, it can be seen that the dynamic strain in the tethered cable is influenced by the material parameters of the cable, the initial pretension, and the bottom excitation frequency and amplitude. In order to investigate the effects of parameters such as the length, material parameters, excitation frequency, and excitation amplitude on the dynamic tension of the cable, the working conditions shown in
Table 1 are designed. According to studies [
18,
19], the ship’s heave motion can be approximated as a simple harmonic, so the excitation frequency of the bottom end of the tethered UAVs is set between 0.1 Hz and 0.3 Hz with an interval of 0.05 Hz. The excitation amplitude ranges from 1 m to 5 m with 1 m. The linear density of the cable
is 0.02 kg/m (not changing with the cross-sectional area of the cable).
Taking condition 1.1 as an example, the excitation amplitude
Ae = 2 m and the excitation frequency
= 0.1 Hz, 0.15 Hz, 0.2 Hz, 0.25 Hz, 0.3 Hz. The relationship between the maximum absolute value of the tangential, normal, and dynamic strain with frequency and the cable’s body coordinates is shown in
Figure 4.
From
Figure 4, it can be obtained that with constant pretension and excitation amplitude, the maximum absolute value of the top dynamic strain is less than that of the bottom dynamic strain at low excitation frequencies, and the maximum absolute value of the top dynamic strain is slightly more than that of the bottom dynamic strain at high excitation frequencies. However, at the excitation frequency of 0.25 Hz, an anomaly occurs. It can be seen that the maximum absolute value of tangential, normal, and dynamic strain is significantly larger than those at other excitation frequencies regardless of the cable’s body coordinates.
In order to analyze the anomalies appearing at the excitation frequency of 0.25 Hz, the excitation frequency change interval is shortened to 0.1 Hz, and the top dynamic strain time course curve is subjected to spectral analysis with a sampling frequency of 30 Hz and 2048 sampling points. The number and magnitude of the amplitudes obtained from the spectrum analysis are as shown in
Figure 5, and the frequency spectrum curves according to different excitation frequencies are as shown in
Figure 6.
When the tethered UAV hovers, the structure is similar to the cable-stayed bridge [
20]. It can be approximated that the natural frequency of the tethered UAVs in the hovering state is
. This equation is substituted by the parameters of working condition 1.1; the first natural frequency of the hovering tethered UAVs is predicted to be about 0.1235 Hz.
Figure 4 shows that a rich dynamics phenomenon emerges when the excitation frequency crosses the first and second natural frequencies.
Figure 5a shows that when the excitation frequency approaches the first natural frequency, the amplitude jumps from the single-periodic component to the multi-periodic component, and subharmonic responses, the multiply periodic response, and subharmonic response appear. Subsequently, the excitation frequency moves away from the intrinsic frequency, and in the response of
Figure 5b, it can be seen that in addition to the components of the excitation frequency, there are also components of the first and second natural frequencies, and the subharmonic response disappears. The amplitude jumps again when the excitation frequency is close to the second natural frequency, and the multiplicative periodic and subharmonic responses appear. When the excitation frequency is far from the second natural frequency, the first and second natural frequencies and the excitation frequency components exist in addition to the excitation frequency, and the multiply periodic response disappears. In order to avoid the amplitude jump phenomenon, the first natural frequency of the tethered UAVs must be made as large as possible by changing the relevant parameters, such as the length of the tethered cable and the tension provided at the tip, thus avoiding the excitation frequency close to the natural frequency.
Figure 7 shows the curve of dynamic tension of the top with time in the first few seconds. The time of 0 s is the first time that the bottom of the tethered UAV is excited. It can be seen that in the first few seconds, there are two sudden changes in the top dynamic tension. The first is because the longitudinal wave generated by the bottom excitation reaches the top, and the second is because the transverse wave reaches the top. Because of the existence of the transverse wave, the shape of the cable is changed, resulting in sudden changes in the dynamic tension. In the second change, the dynamic tension increases rapidly, then decreases, and then the top dynamic tension is affected by both transverse and longitudinal waves.
Figure 8 shows the distribution of the dynamic tension in the cable with time and space for different excitation frequencies. It can be found that the variation in the cable’s dynamic tension consists of both standing waves and traveling waves. When the excitation frequency is close to the natural frequency, the standing wave dominates and forms a large response at the end of the tethered UAV. Moreover, when the excitation frequency is far from the natural frequency, the traveling wave dominates, and the end response of the tethered UAV is smaller.
It can be seen from
Figure 8b,c that the maximum dynamic tension occurs not only at the upper and lower end of the cable but also in the middle of the cable at a certain time.
Then, the variation in dynamic cable tension with excitation amplitude is analyzed. Taking working condition 1.1 as an example, the excitation frequency
= 0.15 Hz and the excitation amplitude
Ae = 1 m, 2 m, 3 m, 4 m, and 5 m. The strain versus time curves at different locations are shown in
Figure 9.
It can be seen from
Figure 9 that under the same excitation frequency, with the increase in excitation amplitude, the maximum absolute value of dynamic strain at each point of the tethered cable gradually increases.
Then, the variation in dynamic tension in the cable with cable length and pretension is analyzed. Taking working conditions 1.2 and 1.3 as examples, under the excitation frequency
= 0.12 Hz and the excitation amplitude
Ae = 2 m, the distribution of dynamic tension with time and cable position under the two working conditions is shown in
Figure 10.
It can be seen from
Figure 10 that by reducing the length of the cable and increasing the pretension, the sudden increase in tension and end point aggregation caused by the excitation frequency close to the natural frequency can be effectively avoided. Because the tethered UAVs need to maintain a certain height when working, increasing the lift of the UAV is an effective way to keep the tethered UAVs operating for a long time.
Finally, the variation in dynamic tension with cross-sectional cable area and elastic modulus is analyzed. Taking working conditions 1.4 and 1.5 as examples, under the excitation frequency
= 0.19 Hz and the excitation amplitude
Ae = 2 m, the distribution of dynamic tension with time and cable position under these two working conditions is shown in
Figure 11.
When the linear density remains unchanged, increasing the elastic modulus and the cable’s cross-sectional area will increase the longitudinal wave’s velocity. It can be seen from
Figure 11 that the maximum amplitude tends to increase.
4. Experiment
In order to conduct experimental research on tethered UAVs, it is necessary to determine the experimental model and experimental devices. According to the
theorem, the relationship between the parameters of the model and the prototype is given by the implicit function in Equation (56):
The dimensional matrix is shown in
Table 2.
Five independent
criteria are identified as follows:
In addition,
is:
Because the forces on the tethered UAVs are mainly analyzed, the dimensionless equation is changed to:
The parameters with subscript
represent the model’s parameters, and the parameters with subscript
represent the parameters of the prototype. The prediction equation is:
Assuming that the scale ratio of length is
and the Young’s moduli of the prototype and model are the same, the scale ratio relationship between model and prototype parameters can be obtained, as shown in
Table 3.
As shown in
Figure 12, a simplified model is designed to study the influence of the frequency of periodic displacement excitation of the mooring point on the tension in the tethered UAVs. The frequency converter applies vertical harmonic displacement excitation at the lower end of the cable through the fixed pulley, the upper and lower tension sensors measure the tension at the end of the mooring line, and five FBG sensors are used to measure the strain at all parts of the cable. The UAV platform keeps hovering, and the cable keeps a small curvature. Because it is difficult to meet all the scale conditions, it is mainly to ensure that the material parameters of the prototype and model cable are the same. Model and prototype physical parameters of the tethered UAVs and cables are shown in
Table 4, and an experimental diagram of the tethered UAVs model is shown in
Figure 13.
When the amplitude of the crank slider mechanism is 0.1 m, the excitation frequency is changed to obtain the time-varying curves of the tension at the bottom and top of the tethered UAVs, as shown in
Figure 14.
The distribution of maximum dynamic tension of tethered UAVs along with different cable positions under different excitation frequencies is shown in
Figure 15.
The spectrum of the dynamic tension at the top under different excitation frequencies is shown in
Figure 16.
It can be seen that the bottom end excitation will lead to two and three times the fundamental frequency components of the top tension, and there is a subharmonic response of the top tension at the excitation frequencies of 0.1 Hz and 0.2 Hz. At the same time, it can be seen from
Figure 15 that when the excitation frequency remains unchanged, the maximum dynamic tension at different positions of the tethered UAVs system does not increase nonlinearly but shows a trend of increasing fluctuation, which is consistent with the theoretical results.