Research Article | Open Access

# Backstepping Control Using Barrier Lyapunov Function for Dynamic Positioning Control System with Passive Observer

**Academic Editor:**Jorge Rivera

#### Abstract

This paper presents a backstepping controller using barrier Lyapunov function (BLF) for dynamic positioning (DP) system. For safety reasons, the position and heading of DP ship are to be maintained in certain range. Thus, in this paper, a control law based on BLF and backstepping technique is proposed to limit the position and heading. The closed-loop system is proved stable in the sense of Lyapunov stability theories. In addition, since the velocities of ship are not measurable and the wave frequency (WF) motion is unavailable, a passive observer is adopted to estimate the velocities and the effect of WF motion. The simulation results show that the proposed controller can limit the position and heading of the vessel in a predefined range and verify the performance of the proposed controller and the passive observer.

#### 1. Introduction

With the exploration of ocean resources such as oil and gas, dynamic positioning (DP) control system is increasingly applied in deep sea drilling operations. The number of ships with DP control system has risen to thousands, and the trend continues. DP ship is defined as a ship that only depending on active thrusters maintains its position and heading (fixed location or predetermined track) (see [1] (Srensen A J, 2011)). DP control system is a complex nonlinear system with uncertain disturbances, parameter uncertainties, input saturation, time-delay, and so on. Thus, the study of DP control system focuses on nonlinearity, unmodeled dynamics, unknown environmental disturbances, input saturation, time-delay, and so on (see [1–4] (Srensen A J, 2005; 2011; Rabanal O M R, Brodtkorb A H, Breivik M, 2016; Xia G, Xue J, Jiao J et al., 2016)). Among those control strategies, backstepping method is a common technique to achieve a control law via defining error variables and a corresponding Lyapunov function of each subsystem to ensure system stability. The paper in [5] (Kim Y S, Kim J, Sung H G, 2016) applied nonlinear backstepping technique in weather-optimal controller designed for DP control system to maintain ship in a fixed position. The study in [6] (Xia, G., Xue, J., Jiao, J, 2018) combined backstepping technique with fuzzy approach to resolve the problem of time-delay and uncertainties in DP control system. The authors in [7] (Morishita H M, Souza C E S, 2014) developed an attenuating control law on basis of observer backstepping methodology for DP control system considering the saturation of actuators.

Moreover, in marine operation, the position of ship with DP control system shall be limited in certain range to ensure safe operation (see [2, 8] (Srensen A J, 2005; Liu X, Abrahamsen B, 2010)). And the study of state constraints has drawn much attention in area of control system (see, for instance, [9–11] (Liu Y J, Tong S, 2016; 2017; Wei H, Zhao Y, Sun, 2017)). Using diffeomorphisms to transform nonlinear constrained control problem into an unconstrained one, the authors in [12] (Doria-Cerezo A, Acosta J A, Castano A R et al., 2014) applied nonlinear state and input constraints control theories to the DP control system to ensure safe operation. To ensure state constraints, BLF has been widely accepted as an effective scheme since this method needs less restrictive initial conditions comparing with other methods, such as reference setting method, model predictive control, and control theories based on invariant set or admissible set. In order to tackle the nonlinearities of control signals and guarantee a safe operation, the study in [13] (Tu F, Ge S S, Choo Y S et al., 2017) presented an adaptive neural constrained control with nonlinear adaptive filter based on BLF for DP of an accommodation ship in deep water. To tackle the constraint problems, the study in [14] (Kong L, He W, Yang C et al., 2018) utilized an asymmetric time-varying BLF to the design of adaptive fuzzy neural network control, proved that closed-loop system was stable by Lyapunov stability theory, and verified the effectiveness of proposed control by comparative simulations. Combining the BLF with backstepping technique is one method of the application of BLF to tackle the problem of state constraint in control areas. For the purpose of solving the problem of the error constraint, the authors in [15] (Ghommam J, Ferik S E, Saad M, 2018) proposed a new robust adaptive controller incorporating tan-BLF with the backstepping control and guaranteed the proposed control law semi-global uniform ultimately bounded stability using Lyapunov stability theory. With the aim of solving the problem of output constraint, the study in [16] (Zheng Z, Huang Y, Xie L et al., 2018) combined an asymmetric time-varying BLF with backstepping technique, employed the hyperbolic tangent function to approximate the saturation function to overcome the input saturation, and verified the effectiveness of the proposed controller by simulation results.

In practice, the wave frequency (WF) motion of DP ship is the most interesting aspect of DP such as estimating the first order WF motions of the ship and excluding it from the feedback loop. In order to estimate the position, heading, and velocities of ship, the authors in [17] (Du J, Hu X, Liu H et al., 2015) proposed a high-gain observer and an adaptive robust output feedback control scheme to estimate velocities and unknown parameters. Using disturbance observer to estimate unknown time-varying disturbances, the study in [18] (Du J, Hu X, Krstić M et al., 2016) came up with a robust nonlinear control law for DP control system to handle input saturation. For the sake of improving performance of model-based observers for DP ship during transients, the authors in [19] (Svenn A. Værn, Brodtkorb A H, Skjetne R et al., 2017) developed a model-based observer with time-varying gains, and the performance of designed observer was verified on research ship Gunnerus.

The main contribution of this paper is the design of a backstepping control law using BLF for DP control system. With this controller, the close-loop system is proved stable in the sense of Lyapunov stability, and the position of the vessel is guaranteed to be limited in a predefined range. Besides, a traditional backstepping control law is compared to the newly designed control law in order to intuitively reveal the advantages of the two control laws. In addition, since the velocities of DP ship are not measurable and the WF motion is unavailable, a passive observer is adopted to estimate the low frequency (LF) motion of DP ship and the effect of WF motion. The rest of the paper is organized as follows. Problem formulation is presented in Section 2. Using BLF, a backstepping control law is developed based on Lyapunov stability theories, the traditional backstepping control law is designed, and a passive observer is introduced in Section 3. The stability of designed control law is proved in Section 4. MATLAB simulation results demonstrate the effectiveness of the designed controller and observer in Section 5. Finally, the conclusion is summarized in Section 6.

#### 2. Problem Formulation

This section describes the mathematical model of DP ship, and the control objective is introduced in this part.

##### 2.1. The Kinematics and Dynamics of DP Ship

In order to analyze the dynamics of DP ship, the earth-fixed reference frame denoted as and body-fixed reference frame denoted as are introduced and shown in Figure 1.

As shown in Figure 1, ship motion is described in 6-degree-of-freedom (DOF) model including three-translational motion and three-rotational motion along the , , and axes. DP model can be derived under the assumption of low speed to keep position; it is widely accepted in the research of DP control system that only three-horizontal motion is considered. The position vector in earth-fixed reference is defined as and the velocity in body-fixed reference is defined as . And the kinematics of DP ship is described aswhere is the rotation matrix transforming body-fixed velocities into earth-fixed velocities with the properties of , , represents identity matrix, and is expressed by

In DP control system, the model of ship is often separated into LW model and WF model (see [20] (Fossen, T. I, 2002)). The 3-DOF model of DP ship in LF motion and WF motion is written aswhere is the WF motion state vector; , and are the constant matrices; represents the zero-mean Gaussian white noise; is position and orientation vector due to one-order wave-induced motion; is the total matrix including inertia mass matrix and added mass matrix; is the hydrodynamic damping matrix; is control force caused by the thrusters of DP ship; is environmental disturbances such as wind, waves, and ocean currents and can be determined by formulas of environmental forces; is the measurement vector, and is constant matrix.

##### 2.2. Control Objective

For DP ship, the horizontal displacement of DP ship should be maintained in certain range in order to ensure safety operation. That is, the position in surge and sway should be limited such as , and the yaw angle is always in region of . One control objective is to design a control law that steers ships to the desired position and heading in the condition of constrained position and heading with the environmental disturbances. The other objective is the designed control law maintaining the position and head of ship at the desired position, heading in the condition of constrained position, and heading under the effect of environmental disturbances.

#### 3. Backstepping Control Using BLF

In this section, the backstepping control law is to be designed using BLF to steer DP ship to the desired position and heading under state constraints and disturbances. And the close-loop system of the proposed control scheme is proved stable in sense of Lyapunov stability theories. The traditional backstepping control law is compared to the proposed control law with BLF in this section. Moreover, a passive observer is adopted in this section to estimate the available and unavailable states of DP control system.

##### 3.1. Backstepping Control Law Using BLF

The basic knowledge of BLF is introduced before designing the control law.

*Definition 1 (see [21, 22] (Zhang T, Wang N, Xia M., 2017; Ren B, Ge S S, Tee K P et al., 2010)). *A BLF is a scalar function , defined with respect to the system on an open region containing the origin; BLF is continuous, positive definite and has continuous first-order partial derivatives at every point of ; the BLF has the property that as approaches the boundary of and satisfies , along the solution of with and some positive constant .

In this paper, the BLF is chosen aswhere is the error; is the positive constraint value on ; that is .

Lemma 2 (see [21, 22] (Zhang T, Wang N, Xia M., 2017; Ren B, Ge S S, Tee K P et al., 2010)). *For any positive constant , let and be open sets. Consider the systemwhere is the state and function is piecewise continuous in and locally Lipschitz in , uniformly in on . Supposing that there exist two continuously differentiable and positive definite functions and in their respective domain, so thatwhere are class functions. Let and . If the following inequality establisheswhere and are positive constants, then remains bounded and .*

Lemma 3 (see [22] (Ren B, Ge S S, Tee K P et al., 2010)). *For any constant , when , the following inequation is held:*

Based on the above introduction, expanding the one-dimensional theory of BLF into three-dimension, a recursive backstepping control law based on BLF is designed in the following steps.

*Step 1. *The desired position and heading vector is defined as . And the position and heading error variable is defined asCombining with (1) and (9), the time derivate of is derived asAnd the virtual control law is designed aswhere is positive constant vector and it is the constraint condition on , and is a positive definite and diagonal matrix.

*Step 2. *The velocity error variable is defined asSo (10) is rewritten asCombining with (3) and (12), the time derivate of is derived aswhere is the time derivate of and is computed asIn view of the BLF theory to solve the problem of state constraint, thus the BLF is utilised into the design of control law. Therefore, the control law is designed aswhere is the item based on BLF to handle with the problem of state constraint, and is a positive definite and diagonal matrix.

##### 3.2. Traditional Backstepping Control Law

The traditional backstepping control law [20, (Fossen, T. I, 2002)] is designed in the following steps.

*Step 1. *The position error is considered as , and the error dynamics is defined aswhere , are the position and velocity vectors, respectively.

Combining with (1) and (17), the time derivate of is determined as

Thus, the virtual control vector is chosen aswhere is a positive definite and diagonal matrix.

*Step 2. *The velocity error is defined as , and the error dynamics is defined asCombining with (3) and (20), the time derivate of is derived asCombining with (18) and (19), the time derivate of is determined asAnd the traditional backstepping control law is designed aswhere parameter matrix is a positive definite and diagonal matrix.

##### 3.3. Passive Observer

For DP ship, the velocities are not measurable and the WF motion is unavailable and this is necessary in controller designing. In order to obtain the changing tendency of the velocities and the effect of WF motion, a passive observer is adopted in this section [20, (Fossen, T. I, 2002)].

The passive observer is expressed aswhere are the estimation state vectors; is the observer matrix and , are the gain matrices; is the vector of estimation error; and is the element of vector . And the passive observer was guaranteed uniformly globally exponential stability detailed in [20, (Fossen, T. I, 2002)].

#### 4. Stability Analysis

##### 4.1. Stability Analysis for Backstepping Control Law Using BLF

Based on the above introduction of BLF, the first BLF candidate is chosen as

So, the time derivate of is derived by combining with (9), (12), and (25) as

Substituting (11) into (26), can be rewritten as

According to the positive definite property of , when is zero, thus, the variable of is stabilization.

The second Lyapunov function candidate is chosen as

So the time derivate of is determined as

Substituting (10), (27) into (29), is derived as

Due to the fact that mass matrix is skew symmetric matrix and , thus, ; therefore, is rewritten as

Substituting the designed control law (16) into (30), is determined as

So the variable of is stabilization when . Therefore, the close-loop system under control of proposed control law is stable in the condition of constrained position and heading.

##### 4.2. Stability Analysis for Traditional Backstepping Control Law

The Lyapunov function candidate is selected as

Combining with (17), (18), and (20), the time derivate of is computed as

Substituting (19) into (34), is determined as

Due to positive definite property of , thus when is zero. Therefore, the variable of is stabilization.

The Lyapunov function candidate is selected as

So the time derivate of is determined as

Substituting (18), (21), and (35) into (37), is derived as

Substituting (23) into (38), is determined as

Therefore, the variable of is stabilization when . So the control system under traditional control law is stability.

#### 5. Simulation Study

Based on LF motion of DP ship in Section 2, the mass matrix of ship introduced in (3) is given as

The damping matrix of ship introduced in (3) is given as

The desired position and heading are chosen as . The expected position is defined as and the expected velocity is defined as . The initial values of position and heading vector are zero, and the original values of velocity vector are zero.

Thus, to generate a smooth desired path, the 2-order filter [20, (Fossen, T. I, 2002)] is adopted aswhere is the natural frequency, is the relative damping ratio, and is the designed parameter. In simulation, , , and are chosen as . Thus the desired states are generated smoothly and the figure of the desired states is shown in Figure 2. The desired position, desired heading, and desired velocities in three directions are smooth and converged to the given value of within . Besides, the bias term introduced in (3) including wind, waves, and ocean currents disturbances is computed in [6, 20] (Xia, G., Xue, J., Jiao, J, 2018, Fossen, T. I, 2002).

What is more, the constant vector considered in Section 3.1 is chosen as , which means that the position error of ship is allowed to move within a square with a side length of and the heading error of the ship may vary within range. And the gain matrices of backstepping control using BLF mentioned in Section 3.1 are chosen as

The gain matrices of traditional backstepping control introduced in Section 3.2 are listed as

The parameters in passive observer introduced in Section 3.3 are selected as

There are four groups of simulation results in every element of position and velocity vectors: one group is the desired position, desired heading, and desired velocities with green solid line; the other one is under presented controller with red chain dotted line; the third one is under traditional controller with cyan dotted line; and last one is under passive observer with red dashed line shown in Figures 3, 4, 5, 6, 7, and 8.

Figures 3, 4, and 5 are the time history of position in three directions. From Figure 3, position smoothly moves along the desired position in the condition of constrained positon from and environmental disturbances. And the position maintains at the desired position under the effect of environmental disturbances. Thus, the two control objectives are implemented under the control law . Also position smoothly and precisely moves to the desired positon and maintains the desired position. And the estimated state precisely and smoothly approaches the position in -direction. But from the partial enlarged details in Figure 3, since the WF motion in -direction is nonlinear motion and affects the whole motion of ship, the states and estimated state oscillate slightly. However, is smoother and more precise than . In this view, the -direction position under presented control law is more precise than under traditional control law .

With the same situation of Figure 4, position smoothly moves along the desired position in the condition of constrained positon from and environmental disturbances and maintains the desired position under the effect of environmental disturbances. Thus, the two control objectives are implemented under the control law . Also position smoothly and precisely moves to the desired position and maintains the desired position. And the estimated state precisely and smoothly approaches the position in -direction. But from the partial enlarged details in Figure 4, since the WF motion in -direction is nonlinear motion and affects the whole motion of ship, the states and estimated state oscillate slightly. However, under the presented control law is smoother and more precise than under control law .

In Figure 5, the heading angle smoothly moves along the desired heading angle in the condition of constrained heading angle from and environmental disturbances and maintains the desired heading angle in the case of constrained heading angle and disturbances. Thus, the two control objectives are implemented under the control law . Also, the heading angle smoothly moves along the desired heading angle and maintains the desired heading angle. Furthermore, the rising time of is longer than the rising time of . And the estimated state smoothly approaches to the heading angle. But from the partial enlarged details in Figure 5, since the WF motion in -direction is nonlinear motion and affects the whole motion of ship, the states and estimated state oscillate slightly. In this view, under the presented control law is faster, smoother, and more precise than under control law .

Figures 6, 7, and 8 are velocities in three directions. From Figure 6, the velocities and are smoothly and precisely moving along the desired velocity. And the estimated velocity shows slight oscillation in partial enlarged details due to adding the estimation of nonlinear WF velocities in direction of . With the same condition in Figure 7, the velocities and are smoothly and precisely moving along the desired velocity. And the estimated velocity oscillates slightly in partial enlarged details due to adding the estimation of nonlinear WF velocities in direction of . From Figure 8, the velocities and are smoothly moving along the desired velocity . But the top value of velocity is not arriving at the desired top value of desired velocity. In this view, the velocity is more precise than . Besides, there are some oscillations in estimated velocity due to adding the estimation of nonlinear WF velocities in -direction.

In addition, the simulation results of control forces and moment in three directions under control laws and are shown in Figure 9. The control forces and moment are smoother under the control law than control law .

Therefore, the position, heading, and velocities under presented control law show better performance than traditional control law .

#### 6. Conclusions

In this paper, in order to limit the position and heading into certain range for safety reasons, a control law is proposed combining backstepping technique with BLF. Utilizing an error constrained vector in BLF, the virtual control law and backstepping control law are derived to limit the position and heading. And the designed control law is proved stable in the sense of Lyapunov stability theories. The simulation results show that the position and heading under proposed control law move along the desired position and heading in the constrained position, heading, and the effect of environmental disturbances, and maintain the desired position and heading under the constraint position and the effect of environmental disturbances. Therefore, the two control objectives are achieved under the designed control law. And a traditional backstepping control law is also introduced as comparison. Through the comparisons of the two methods, the position and heading appear faster, smoother, and more precise under the proposed control law than traditional control law. Therefore, the proposed control law shows better performance than traditional control law. Moreover, a passive observer is adopted to estimate position and velocities of ship. Considering the nonlinear WF motion of observer, there are small oscillations during control of position, heading angle, and estimation of position and velocities, and the future work is to research this phenomenon.

#### Data Availability

The authors declare that data sources are public and are referred to in [6, 20] (Xia, G., Xue, J., Jiao, J, 2018, Fossen, T. I, 2002).

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this study.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 51609046), the Fundamental Research Funds for the Central Universities (HEUCFM170403, HEUCFJ180404), and the 7th Generation Ultra Deep Water Drilling Unit Innovation Project.

#### References

- A. J. Sørensen, “A survey of dynamic positioning control systems,”
*Annual Reviews in Control*, vol. 35, no. 1, pp. 123–136, 2011. View at: Publisher Site | Google Scholar - A. J. Sørensen, “Structural Issues in The Design And Operation of Marine Control Systems,”
*Annual Reviews in Control*, vol. 29, no. 1, pp. 125–149, 2005. View at: Publisher Site | Google Scholar - O. M. R. Rabanal, A. H. Brodtkorb, and M. Breivik, “Comparing Controllers for Dynamic Positioning of Ships in Extreme Seas,”
*IFAC-PapersOnLine*, vol. 49, no. 23, pp. 258–264, 2016. View at: Google Scholar - G. Xia, J. Xue, J. Jiao et al., “Adaptive Fuzzy Control for Dynamic Positioning Ships with Time-Delay of Actuator,” in
*Proceedings of the 2016 OCEANS MTS/IEEE Monterey, OCE 2016*, pp. 1–6, California, Calif, USA, September 2016. View at: Google Scholar - Y.-S. Kim, J. Kim, and H.-G. Sung, “Weather-Optimal Control of A Dynamic Positioning Vessel Using Backstepping: Simulation And Model Experiment,”
*IFAC-PapersOnLine*, vol. 49, no. 23, pp. 232–238, 2016. View at: Google Scholar - G. Xia, J. Xue, and J. Jiao, “Dynamic Positioning Control System with Input Time-Delay Using Fuzzy Approximation Approach,”
*International Journal of Fuzzy Systems*, vol. 20, no. 2, pp. 630–639, 2018. View at: Publisher Site | Google Scholar | MathSciNet - H. M. Morishita and C. E. S. Souza, “Modified Observer Backstepping Controller for A Dynamic Positioning System,”
*Control Engineering Practice*, vol. 33, no. 33, pp. 105–114, 2014. View at: Publisher Site | Google Scholar - X. Liu and B. Abrahamsen, “Risk Assessment of Riser Operations on DP Drilling Units,” in
*Proceedings of the ASME 2010 29th International Conference on Ocean, Offshore and Arctic Engineering, OMAE2010*, pp. 823–829, China, June 2010. View at: Google Scholar - Y.-J. Liu and S. Tong, “Barrier Lyapunov Functions-Based Adaptive Control for A Class of Nonlinear Pure-Feedback Systems with Full State Constraints,”
*Automatica*, vol. 64, pp. 70–75, 2016. View at: Publisher Site | Google Scholar | MathSciNet - Y.-J. Liu and S. Tong, “Barrier Lyapunov Functions for Nussbaum Gain Adaptive Control of Full State Constrained Nonlinear Systems,”
*Automatica*, vol. 76, pp. 143–152, 2017. View at: Publisher Site | Google Scholar | MathSciNet - W. He, Z. Yin, and C. Sun, “Adaptive Neural Network Control of a Marine Vessel With Constraints Using the Asymmetric Barrier Lyapunov Function,”
*IEEE Transactions on Cybernetics*, vol. 47, no. 7, pp. 1641–1651, 2016. View at: Publisher Site | Google Scholar - A. Doria-Cerezo, J. A. Acosta, A. R. Castano et al., “Nonlinear State-Constrained Control. Application to The Dynamic Positioning of Ships,” in
*Proceedings of the 2014 IEEE Conference on Control Applications, CCA 2014*, pp. 911–916, France, October 2014. View at: Google Scholar - F. Tu, S. S. Ge, Y. S. Choo et al., “Adaptive dynamic positioning control for accommodation vessels with multiple constraints,”
*IET Control Theory & Applications*, vol. 11, no. 3, pp. 329–340, 2017. View at: Publisher Site | Google Scholar - L. Kong, W. He, C. Yang et al., “Adaptive Fuzzy Control for A Marine Vessel with Time-Varying Constraints,”
*IET Control Theory Applications*, vol. 12, no. 10, pp. 1448–1455, 2018. View at: Google Scholar - J. Ghommam, S. El Ferik, and M. Saad, “Robust Adaptive Path-Following Control of Underactuated Marine Vessel with Off-Track Error Constraint,”
*International Journal of Systems Science*, vol. 49, no. 7, pp. 1540–1558, 2018. View at: Publisher Site | Google Scholar | MathSciNet - Z. Zheng, Y. Huang, L. Xie, and B. Zhu, “Adaptive Trajectory Tracking Control of a Fully Actuated Surface Vessel With Asymmetrically Constrained Input and Output,”
*IEEE Transactions on Control Systems Technology*, vol. 26, no. 5, pp. 1851–1859, 2018. View at: Publisher Site | Google Scholar - J. Du, X. Hu, H. Liu, and C. L. Chen, “Adaptive Robust Output Feedback Control for A Marine Dynamic Positioning System Based on A High-Gain Observer,”
*IEEE Transactions on Neural Networks and Learning Systems*, vol. 26, no. 11, pp. 2775–2786, 2015. View at: Publisher Site | Google Scholar | MathSciNet - J. Du, X. Hu, M. Krstić et al., “Robust Dynamic Positioning of Ships with Disturbances under Input Saturation,”
*Automatica*, vol. 73, no. C, pp. 207–214, 2016. View at: Publisher Site | Google Scholar | MathSciNet - S. V. Svenn, A. H. Brodtkorb, R. Skjetne et al., “Time-Varying Model-Based Observer for Marine Surface Vessels in Dynamic Positioning,”
*IEEE Access*, vol. 5, pp. 14787–14796, 2017. View at: Publisher Site | Google Scholar - T. I. Fossen,
*Marine Control Systems-Guidance, Navigation, And Control of Ships, Rigs And Underwater Vehicles, Marine Cybernetics*, Norway, 2002. - T. Zhang, N. Wang, and M. Xia, “Adaptive Neural Output Feedback Control of Stochastic Nonlinear Systems with Unmodeled Dynamics,”
*Mathematical Problems in Engineering*, vol. 32, no. 1, Article ID 529862, 13 pages, 2015. View at: Publisher Site | Google Scholar | MathSciNet - B. Ren, S. S. Ge, K. P. Tee et al., “Adative Neural Control for Output Feedback Nonlinear Systems Using A Barrier Lyapunov Function,”
*IEEE Transactions on Neural Networks and Learning Systems*, vol. 21, no. 8, pp. 1339–1345, 2010. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2019 Guoqing Xia et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.