Next Article in Journal
Nonlinear Finite Element Modelling of Thermo-Visco-Plastic Styrene and Polyurethane Shape Memory Polymer Foams
Previous Article in Journal
Vibration Control of a High-Rise Slender Structure with a Spring Pendulum Pounding Tuned Mass Damper
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Novel Singular-Perturbation-Based Adaptive Control with σ-Modification for Cable Driven System

1
School of Automation Science and Electrical Engineering, Beihang University, Beijing 100191, China
2
Ningbo Institute of Technology, Beihang University, Ningbo 315800, China
3
Department of Electrical and Electronic Engineering, University of Nottingham, Nottingham NG7 2RD, UK
*
Author to whom correspondence should be addressed.
Submission received: 1 February 2021 / Revised: 9 February 2021 / Accepted: 11 February 2021 / Published: 28 February 2021
(This article belongs to the Section Precision Actuators)

Abstract

:
Due to their large working space and fast response, cable driven systems have been widely applied in manufacturing, robotics and motion simulators, etc. However, the cable is flexible and tends to resonate at high frequencies, which raises challenges for the motion control of the cable driven system. To solve this problem, this paper proposes a singular-perturbation-based adaptive control method with σ -modification. Taking advantage of the multi-time scale characteristics, the flexible system is approximately decomposed into two subsystems, and then the damping compensation is designed in the boundary layer subsystem to enhance the tension stability. In addition, estimated parameters drift may occur for the reduced-order system. Thus, the σ -modification is proposed to ensure that the tracking and estimation errors converge to a bounded residual set. The Lyapunov stability theorem proves that the closed-loop system is stable and errors are ultimately uniformly bounded. A research prototype of a twin-motor cable driven system is developed, and experimental investigation is conducted on it. The experimental results show that the proposed control method can effectively suppress cable resonance at high frequencies. Compared with the conventional adaptive control method, it can significantly increase the system bandwidth.

1. Introduction

Cable driven systems have been widely used in lifting [1,2], telescopes [3], robotics [4,5], motion simulators [6,7] and other fields [8,9]. Different from a rigid transmission system, it has two inherent properties [5,10]: (1) the cable can only be pulled but not pressed, so redundant actuation must be adopted to make the end effector fully controllable; (2) the cable is elastic. At high frequencies, resonance often occurs, which severely affects the performance and stability of the system and brings great challenges to the system control. In existing literature, internal tension compensation is used to solve the problem of unidirectional force, and the motion control is achieved by treating the cable driven system as a rigid system [11,12,13]. Typical methods include the computed torque control method [14,15,16] and robust adaptive control method [17,18,19,20,21,22,23].
Inspired by the feedback linearization theory, Alp et al. [24] proposed a computed torque method based on the nominal model for the position control of a suspended robot, and proved the stability of the system. However, due to the parameter uncertainty in the system, the control performance is greatly compromised. To solve this problem, robust adaptive control is one good solution. To reduce the negative effect of dynamic parameters uncertainty, Li et al. [18] proposed a robust adaptive control, by adding online estimation of dynamic parameters. Wang et al. [19] proposed a time-delay estimation method, and combined it with terminal sliding mode control to achieve an adaptive integral terminal sliding mode control for cable-driven manipulators. Moreover, Wang et al. [20] introduced the super-twisting algorithm into the adaptive control scheme to ensure that the system can reach the sliding mode surface in limited time. Shang et al. [17] also conducted an in-depth study on the sliding mode control of the cable driven system. They combined the second-order sliding mode control with the synchronization control, and proposed a compound robust control with high precision and robustness. To further improve the control accuracy, Shang et al. [25] introduced the cross-coupled control and parameters adaptive compensation to the synchronization control, and then proposed an adaptive cross-coupled control scheme. Lamaury et al. [26] further extended the adaptive control to the workspace and joint space, and proposed a dual-space adaptive control strategy. The above studies are mainly on the motion control of systems in low dynamics. Because the high frequency characteristics of the cable are filtered, these control methods based on the rigid model can also perform well for the flexible systems.
However, for systems in high dynamic motion, the filtering effect weakens significantly. Thus, it is not equivalent to a rigid system anymore. Specifically, by considering the cable flexibility, the system is converted to a series-connected dual-mass system from a single-mass system. The controller based on a rigid model cannot guarantee the full state stability of a flexible system, especially the cable tension that easily becomes unstable due to insufficient cable damping. This problem has aroused the attention of researchers [27]. Khosravi et al. proposed the cascade control method to keep the tension stable by employing the internal tension control of the cable to avoid tension drift [28,29,30]. Theoretically, this method is feasible. However, it is difficult to collect high-quality tension signals, which affects the tension control. Furthermore, the dual closed loops in the cascade control is quite challenging. This method is poor to implement, especially for high-dynamics motion, a simpler controller is better, and cascade control contradicts this requirement. To solve this problem, the singular perturbation control was proposed [31,32]. The multi-time scales of the flexible system are employed to decompose the system into two subsystems, i.e., the boundary layer and the reduced-order system, and the damping compensation is designed in the boundary layer to keep the tension stable. Compared with the cascade controller, this method is simpler and more practical, and force sensors are not required in the system design.
However, so far the control algorithms of cable driven systems have not considered the drift of estimation parameters, which unavoidably compromises the motion control performance especially in high dynamics. Therefore, this paper proposed a singular-perturbation-based adaptive control with σ -modification to avoid the parameter drift instability, and to give the sufficient conditions for the system stability.
The main contributions of the paper are as follows.
(1)
Based on the singular perturbation theory, an adaptive control method with σ -modified adaptive law is proposed.
(2)
The stability of the flexible system is proved theoretically, and sufficient conditions for system stability are given.
(3)
This method is the first time to be used for the control of a flexible system with high-frequency motion.
(4)
The effectiveness of the proposed method is validated by experiments.
The rest is organized as follows. In Section 2, the schematic structure of the twin-motor cable driven system, a typical cable driven parallel mechanism, is introduced, and kinematics and dynamics models are formulated analytically. In Section 3, based on the singular perturbation theory, an adaptive composite control with σ -modification is proposed. In Section 4, the stability of the closed-loop system is proved, and sufficient conditions for the flexible system are given. In Section 5, the developed twin-motor cable driven system is presented, and an experiments are conducted. The experimental results are utilized to validate the proposed control method. Finally, the research work is concluded in Section 5.

2. Analysis of Flexible System Kinematics and Dynamics

2.1. System Architecture

The studied twin-motor cable driven system is shown in Figure 1. It is composed of motors, an arc guide rail, flexible cable, pulleys and the moving platform. Affected by the unidirectional force property of the flexible cable, the moving platform can only reciprocate along the arc guide under the cooperative traction of the motors on both sides. This system is a typical redundant actuation system. To facilitate subsequent discussion, the fixed coordinate system E: ( x e , y e ) on the base frame and the reference coordinate system R: ( x r , y r ) on the moving platform are set, and their origins O e and O r are coincident with the arc center and the centroid of the moving platform, respectively.

2.2. Analysis of System Kinematics

The vector of the ith cable is represented as
l i E = l i n i E = B i E X E R R E A i R , i = 1 , 2 ,
where l i denotes the cable length between points A i and B i , n i E the unit vector of cable tension, A i R the position vector of point A i , B i E the position vector of point B i , X E the position vector of the origin O r and R R E the rotation matrix of frame R with respect to frame E .
From (1), the velocity mapping between the cable length and the end effector is
L ˙ = J φ T φ ˙ ,
where L ˙ = [ l ˙ 1 , l ˙ 2 ] T , φ denotes the angle between O e O r and the positive x e , and
J φ = r a sin ( α 1 φ β a ) r a sin ( α 2 φ + β a ) , J φ R 1 × 2 ,
r a denotes the radius of the guide rail, α i the angle between the cable tension and the positive x e and β a the structure parameter as shown in Figure 1.
It is worth noting that, affected by the elasticity of the cable, the speeds of the two ends of the same cable are different, and they meet the following conditions
L ˙ a = L ˙ = J φ T φ ˙ = J q T q ˙ ,
L ˙ m = θ ˙ r m ,
where L ˙ a denotes the cable velocity at points A i , L ˙ m the cable velocity at motor connection, q = r a φ , θ the motor rotation angle and r m the radius of the winch,
J q = sin ( α 1 φ β a ) sin ( α 2 φ + β a ) , J q R 1 × 2 .

2.3. Analysis of System Dynamics

2.3.1. Dynamics Model

The force is shown in Figure 2. According to the theorem of momentum, the dynamics of the moving platform is
r c × m a c + I c · ω ˙ = r a , i × T i + r a , i × f i + r c × m g ,
where r c denotes the vector of the centroid of platform, m the platform mass, a c the acceleration vector of the center of mass, I c the moment of inertia of the platform about the centroid of mass, ω ˙ the rotational acceleration, r a , i the position vector of point A i , T i the tension vector, g the gravity acceleration, and f i the friction vector
f i = b v , i φ ˙ + f c , i sign ( φ ˙ ) ,
b v , i is the viscosity coefficient and f c , i is the coulomb friction.
In (7), the cable tensions are
K ( L m L a ) = T ,
where K denotes the cable stiffness, L m and L a the movement displacement of cables at motor and points A i , respectively. The dynamics of motors is
I m θ ¨ + T r m = τ ,
where I m R 2 × 2 denotes the inertia matrix, r m the radius of the winch, τ the motor torque.
Combining (7)–(10), the dynamics of the flexible system are
M q ¨ + B p q ˙ + F sign ( q ˙ ) + G ( q ) = J q K ( L m L a ) ,
I m θ ¨ + K ( L m L a ) r m = τ ,
where
M = m r c 2 + I c / r a 2 , B v = b v , 1 + b v , 2 / r a , F = ( f c , 1 + f c , 2 ) / r a , G ( q ) = G cos ( q / r a ) , G = r c m g / r a .

2.3.2. Dynamics Model Analysis

In (11) and (12), the stiffness K is much larger than other parameters. To facilitate analysis, define
K = k 1 ε 2 E 2 ,
where ε is a small positive parameter, k 1 is on the order of 1 and E 2 is the unit matrix. Substituting (9) and (14) into (12), the dynamics of the motor can be rewritten as
ε 2 r m 1 I m T ¨ + k 1 T r m = k 1 ( τ r m 1 I m L ¨ a ) .
Since ε is much smaller than other parameters, the changing rate of tension T in (15) far exceeds that of q in (11). On the time scale, the variables T and q can be divided into two types: fast and slow. When the fast variable T changes, the slow variable q can be considered frozen. Therefore, in (15), the transfer function from the torque to the cable tension can be approximated as a second-order system. At low frequencies, due to the filtering of the system itself, the high-frequency part of the system response is filtered, and the flexible system is equivalent to a rigid one. As the frequency increases, the system’s ability to filter high-frequency response decreases, and the impact of the flexible cable on the system gradually increases. When a certain frequency is reached, the impact of the flexible cable cannot be ignored, and even damage the system stability. To solve this problem, a composite adaptive control scheme based on singular perturbation is proposed.

3. Composite Control Scheme

Since the system has multiple time scales, the system can be divided into fast (boundary layer) and slow (reduced-order) subsystems, and then the control laws are designed for each subsystem to form the composite control law, i.e.,
τ = τ s + τ f ,
τ f = r m K v ( L ˙ a L ˙ m ) ,
where τ s denotes the slow control law, τ f the fast control law and K v the feedback gains. Define
K v = k 2 ε E 2 ,
where k 2 is on the order of 1. Substituting (16) into (15), the system can be written as
ε 2 r m 1 I m T ¨ + ε k 2 r m T ˙ + k 1 T r m = k 1 ( τ s r m 1 I m L ¨ a ) .
Moreover, from (19), it can be seen that the stability of the fast variable T can be improved by the compensation of fast control law without any impact on the reduced-order system. When ε = 0 , both (15) and (19) can be degenerated into the same reduced-order model as the rigid model.
M e ( q ) q ¨ + N e ( q , q ˙ ) q ˙ + B e q + F e sign ( q ˙ ) + G e ( q ) = J q τ s ,
where
M e ( q ) = r m M + r m 1 J q I m J q T , N e ( q , q ˙ ) q ˙ = r m 1 J q I m J ˙ q T , F e = r m F , B e = r m B , G e ( q ) = r m G ( q ) .
Therefore, the control schemes for the rigid system still work for the reduced-order system. Considering the uncertainty of dynamic parameters, the slow control law is designed as
τ s = τ e x + τ i n ,
τ e x = J q + [ M ^ e ( q ) q a + N ^ e ( q , q ˙ ) q v + B ^ e q ˙ + F ^ e sign ( q ˙ ) + G ^ e ( q ) k b q r ] ,
τ i n = ( E 2 J q + J q ) ξ ,
where τ e x and τ i n denote the effective force and internal force, respectively, ( ^ ) the estimated parameter, J q + the pseudo-inverse and ξ any 2-dimensional column vector,
z = q q d , q v = q ˙ d k a z , q a = q ˙ v , q r = q ˙ q v = z ˙ + k a z ,
q d is the desired trajectory, k a and k b are the adjustable gains. Then combining (20) and (22)–(25) gives the closed-loop dynamics of the reduced-order system
M e ( q ) q ˙ r + [ N e ( q , q ˙ ) + k b ] q r = Φ T ( q , q ˙ , q ˙ d , q ¨ d ) θ ˜ ,
where
Φ ( q , q ˙ , q ˙ d , q ¨ d ) = [ r m q ¨ d k a ( q ˙ q ˙ d ) , r m q ˙ , r m sign ( q ˙ ) , r m cos ( q / r a ) ] T ,
θ ˜ = [ M ˜ , B ˜ p , F ˜ , G ˜ ] T .
To prevent the estimated parameters from drifting, the adaptive law is designed as
θ ˜ ˙ = Γ Φ T ( q , q ˙ , q ˙ d , q ¨ d ) q r σ Γ θ ^ ,
where Γ is the gain matrix, σ a positive constant. According to the above, the designed composite adaptive controller is shown in Figure 3.

4. Stability Analysis

Based on the singular perturbation theory and Lyapunov’s stability theorem, the stability of a close-loop system is analyzed. It includes four parts: singular perturbation standardization, boundary layer stability analysis, reduced-order system stability analysis, and overall system stability analysis.

4.1. Singular Perturbation Standardization

Define the slow, fast and intermediate variables as
x = [ z , q r ] T ,
w = [ r m T T , ε r m T ˙ T ] T ,
τ * = τ s r m 1 I m L ¨ a .
Then from (11) and (19), the standard singular perturbation model of the closed loop is
x ˙ = f ( t , x , w , ε ) = A 1 x + Ψ θ ˜ + A 3 w + B 1 τ * ,
ε w ˙ = g ( t , x , w , ε ) = A 2 w + B 2 τ * ,
where
A 1 = k a 1 0 M e ( q ) 1 [ N e ( q , q ˙ ) + k b ] , Ψ = 0 1 × 4 M e ( q ) 1 Φ T , A 2 = 0 2 E 2 k 1 r m 2 I m 1 k 2 r m 2 I m 1 , B 1 = 0 1 × 2 M e ( q ) 1 J q , A 3 = 0 1 × 2 0 1 × 2 M e ( q ) 1 J q 0 1 × 2 , B 2 = 0 2 k 1 r m 2 I m 1 .
See Appendix A for the detailed derivation. When ε = 0 , the isolated root of w is
w ¯ = h ( t , x ¯ ) = A 2 1 B 2 τ * ,
where w ¯ is the quasi-steady state of w . However, due to the difference of the initial conditions and perturbation of ε , w ¯ deviates from the actual w . The deviation is defined as
y = w h ( t , x ) = w + A 2 1 B 2 τ * .
Substituting (37) into (33) and (34), and combining with (29), the standard singular perturbation model is rewritten as
S : x ˙ = f ( t , x , y + h ( t , x ) , ε ) = A 1 x + Ψ θ ˜ + A 3 y , θ ˜ ˙ = Γ Ξ x σ Γ θ ^ , ε y ˙ = g ( t , x , y + h ( t , x ) , ε ) ε h ˙ ( t , x ) = A 2 y + ε A 2 1 B 2 τ ˙ * ,
where Ξ = [ 0 , Φ T ] .

4.2. Boundary Layer and Stability Analysis

Define the fast-time scale η = ( t t 0 ) / ε , then (38)-3 can be represented as
d y d η = g ( t , x , y + h ( t , x ) , ε ) ε h ˙ ( t , x ) .
From the above formula, the boundary layer is
S b : d y d η = g ( t , x , y + h ( t , x ) , ε = 0 ) = A 2 y ,
where A 2 is a Hurwitz matrix that can be achieved by adjusting the gains k 1 and k 2 . Then for any given positive definite matrix Q , there exists a unique positive definite matrix P that satisfies
A 2 T P + P A 2 = Q .
Then for the boundary layer, choose the Lyapunov function candidate as
V b ( y ) = y T P y .
Along the boundary layer system S b , the fast-time derivative of V b is
V b = V b ( y ) y T g ( t , x , y + h ( t , x ) , ε = 0 ) = ( A 2 y ) T P y + y T P A 2 y = y T Q y λ min ( Q ) y 2 = γ 1 y 2 ,
where λ min ( Q ) = γ 1 denotes the smallest eigenvalue of matrix Q . Hence, as presented in (42) and (43), for all y , V b > 0 and V ˙ b < 0 , i.e., when η , y 0 , and the boundary layer is asymptotically stable.

4.3. Reduced-Order System and Stability Analysis

From (38), the reduced-order model of flexible system is obtained as
S r : x ˙ = f ( t , x , y + h ( t , x ) , y = 0 ) = A 1 x + Ψ θ ˜ , θ ˜ ˙ = Γ Ξ x σ Γ θ ^ .
Then for the reduced-order system, choose the Lyapunov function candidate as
V r ( x , θ ˜ ) = 1 2 x T H x + 1 2 θ ˜ T Γ 1 θ ˜ ,
where H = diag ( k a , M e ( q ) ) . Along the reduced-order system S r , the time derivative of V r is
V ˙ r = V r x T f ( t , x , y + h ( t , x ) , y = 0 ) + V r H T H ˙ + V r θ ˜ T θ ˜ ˙ = x T Λ x σ 2 θ ˜ T θ ˜ σ 2 θ ˜ + θ T θ ˜ + θ + σ 2 θ T θ γ 2 x 2 σ 2 θ ˜ 2 + σ 2 θ T θ x , θ ˜ Q r x θ ˜ + σ 2 θ ^ max T θ ^ max ,
where γ 2 = λ min ( Λ ) ,
Λ = k a 2 k a / 2 k a / 2 k b , Q r = γ 2 0 0 σ / 2 ,
Λ can be kept positive definite by adjustment of k a and k b . Define the sets:
Ω V r = Δ ( x , θ ˜ ) : x T Λ x + σ 2 θ ˜ T θ ˜ σ 2 θ ^ max T θ ^ max ,
Ω V r C = Δ ( x , θ ˜ ) : x T Λ x + σ 2 θ ˜ T θ ˜ > σ 2 θ ^ max T θ ^ max .
Hence, as presented in (45) and (46), for all ( x , θ ˜ ) Ω V r C , V r > 0 and V r ˙ < 0 , i.e., ( x , θ ˜ ) will convergence to the residual set Ω V r , and the reduced-order system is stable.

4.4. Stability Analysis of the Flexible System

According to the above two Lyapunov functions, the following composite Lyapunov function candidate for the flexible system is chosen as
V = ( 1 d ) V r ( x , θ ˜ ) + d V b ( y ) , 0 < d < 1 ,
where d is constant parameter. Then along the original flexible system S, the time derivate of V is
V ˙ = ( 1 d ) V r x T f ( t , x , y + h ( t , x ) ) + V r H T H ˙ + V r θ ˜ T θ ˜ ˙ + d ε V b y T g ( t , x , y + h ( t , x ) , ε ) ε h ˙ ( t , x ) = ( 1 d ) V r x T f ( t , x , y + h ( t , x ) , y = 0 ) + V r H T H ˙ + V r θ ˜ T θ ˜ ˙ + d ε V b y T g ( t , x , y + h ( t , x ) , ε ) , ε = 0 + ( 1 d ) V r x T f ( t , x , y + h ( t , x ) ) f ( t , x , y + h ( t , x ) , y = 0 ) + d ε V b y T g ( t , x , y + h ( t , x ) , ε ) ε h ˙ ( t , x ) g ( t , x , y + h ( t , x ) , ε = 0 ) .
Note that the first two terms in (51) are the time derivatives of V r and V b with respect to the reduced-order system and the boundary layer, respectively. From (43) and (46) we know that they are stable. The remaining two items are the interconnection between the slow/fast subsystem and the original system. Combining (38), (44) and (45) gives
V r x T f ( t , x , y + h ( t , x ) ) f ( t , x , y + h ( t , x ) , y = 0 ) = x T H A 3 y = x T 0 0 J q 0 y x J q y μ 0 x y ,
where μ 0 = sup J q . Combining (38), (40) and (42) gives
V b y T g ( t , x , y + h ( t , x ) , ε ) ε h ˙ ( t , x ) g ( t , x , y + h ( t , x ) , ε = 0 ) = 2 y T P ε A 2 1 B 2 τ ˙ * = 2 ε y T P E 2 0 τ ˙ * 2 ε y P τ ˙ * ,
where P is defined in (41).
Theorem 1.
For continuous and bounded q d , q ˙ d , q ¨ d , q d , there exists positive constants κ 1 , κ 2 , κ 3 and positive function κ 4 ( t ) satisfying
τ ˙ * ( κ 3 + 1 ε κ 2 ) y + κ 1 x + κ 4 ( t ) .
Proof of Theorem 1.
See Appendix B. □
Substituting (54) into (53), we get
2 ε y P τ ˙ * 2 ε y P ( κ 3 + 1 ε κ 2 ) y + κ 1 x + κ 4 ( t ) = ε ( μ 3 + 1 μ 2 ) y 2 + ε μ 1 x y + ε μ 4 ( t ) y ε ( μ 3 + 1 μ 2 ) y 2 + ε μ 1 x y + μ 4 2 ( t ) + 1 4 y 2 ,
where, μ 3 = 2 κ 3 P , μ 2 = 2 κ 2 P , μ 1 = 2 κ 1 P , μ 4 ( t ) = 2 κ 4 ( t ) P . Hence, substituting (43), (46), (52) and (55) into (51), we get
V ˙ x , y Q x y ( 1 d ) σ 2 θ ˜ 2 + d μ 4 2 ( t ) + ( 1 d ) σ 2 θ ^ max T θ ^ max ,
where,
Q = ( 1 d ) γ 2 1 2 ( 1 d ) μ 0 + d μ 1 1 2 ( 1 d ) μ 0 + d μ 1 d ε γ 1 d ( μ 3 + 1 ε μ 2 ) 1 4 d .
When
ε < d ( γ 1 μ 2 ) 1 4 ( 1 d ) μ 0 + d μ 1 2 ( 1 d ) γ 2 + d μ 3 + 1 4 d = ε u ,
Q is a positive definite matrix, which can be easily achieved by adjusting d in this system. Then define the following sets
Ω V = Δ Θ : x , y Q x y + ( 1 d ) σ 2 θ ˜ 2 d μ ¯ 4 2 ( t ) + ( 1 d ) σ 2 θ ^ max T θ ^ max ,
Ω V C = Δ Θ : x , y Q x y + ( 1 d ) σ 2 θ ˜ 2 > d μ ¯ 4 2 ( t ) + ( 1 d ) σ 2 θ ^ max T θ ^ max ,
where Θ = ( x , y , θ ˜ ) , μ ¯ 4 ( t ) = sup μ 4 ( t ) .
Therefore, according to (56), for all ( x , y , θ ˜ ) Ω V , V ˙ can be positive or negative, but for all ( x , y , θ ˜ ) Ω V C , V ˙ < 0 definitely. It means that for all ( x , y , θ ˜ ) Ω V , they either converge or diverge, but once they exceed the set Ω V , ( x , y , θ ˜ ) will remove towards the origin and converge to the set Ω V again. The singular-perturbation-based adaptive control with σ -modification method avoids the drift of estimated parameters, and enables all errors to converge to a residual set. The closed-loop of the flexible system is stable.

5. Experimental Validation

5.1. Experimental Setup

One research prototype of the twin-motor cable driven system is developed as shown in Figure 4. It is mainly composed of two parts: the transmission system and the control system. The transmission system includes motors, cables, an arranging mechanism, a guide rail and sensors. Cables are connected between the motors and the moving platform. To prevent cables from winding irregularly, the cable arrangement mechanism is also designed and installed on each motor, so that the cables can be arranged on the rotating shaft regularly.
The measurement and control system is shown in Figure 5, the hardware includes the host and slave computers, A/D and D/A boards, drivers, motors and sensors. Based on the Labview software, the control program is designed with functions of command input, status detection, data processing and storage. All data in the slave computer will be processed within 1 ms. Therefore, the applicability of the control method should be considered.
Experiments are conducted on the test bed to validate the proposed composite adaptive control with σ -modification ( σ -CAC), and compare with the traditional adaptive control algorithm (TAC). Specifically, it is to verify that the proposed method helps to improve the bandwidth of the cable driven system. The parameters of the two controllers are listed in Table 1.

5.2. Experimental Results and Discussion

The tracking command is
q d = A d [ 1 exp ( 0.1 t 3 ) ] sin ( 2 π f d t ) ,
where A d denotes the amplitude, and f d the frequency. The amplitude is fixed at 10 mm, and experiments with frequencies of f d = 2∼9 Hz have been done, respectively. As shown in Figure 6, both TAC and σ -CAC perform well at 3 Hz, except for the start-up phase. They both can keep the tracking error within ±1mm. However, because of the compensation of the fast control law in the σ -TCT, it has a stronger control effect and further reduces the error within ±0.5 mm. The high-frequency characteristics of the flexible cable show little effect at 3 Hz.
The frequency is further increased to 5 Hz, as shown in Figure 7. Compared with 3 Hz, the performance of the TAC deteriorates. Specifically, the tracking error increases from ±1 to ±1.5 mm, and the phase lags obviously. While the σ -CAC continues to maintain high accuracy, its error is kept within ±0.5 mm, and phase follows closely. The σ -CAC shows advantage at 5 Hz.
The frequency is increased to 7 Hz again, as shown in Figure 8. Compared with 3 and 5 Hz, the performance of the TAC gets worse as the frequency increases. The tracking error increases from the initial ±1 to ±3 mm, and there is a serious phase lag. The σ -CAC showed consistency in 3, 5, and 7 Hz experiments. This consistency is the evidence that the σ -CAC can make the flexible system more stable. Compared with TAC, the σ -CAC significantly improves the performance at 7 Hz.
The detailed comparison results of the two controllers in the frequency domain are recorded in Table 2. It shows that as the frequency increases, their performance deteriorates. However, the deterioration rate of the σ -CAC is much slower than the TAC. At 9 Hz, the σ -CAC can still ensure that the phase lag of system is less than 10 degrees and the amplitude error does not exceed 10%, while, the TAC is already worse than this performance at 6 Hz. The σ -CAC successfully increases the system bandwidth from 5 to 9 Hz. As shown in Figure 9, at 9 Hz, the σ -CAC can still work stably, and the phase lag is about −9.2 degrees.
The above comparative results consistently show that with the increase of the frequency, the influence of the flexible characteristics of the cable on the system control increases. At low frequencies that far away from the cable characteristics, TAC and σ -CAC both perform well. However, at high frequencies, the TAC performs poorly for the cable flexibility. While, the σ -CAC, due to the fast control law compensation, can effectively weaken the flexible cable influence. Hence, the σ -CAC shows obvious advantages in the high-frequency experiments. Figure 6, Figure 7 and Figure 8 and Table 2 directly show that the σ -CAC is significantly effective for increasing the bandwidth of the flexible system.

6. Conclusions

Based on the singular perturbation theory, a σ -modified adaptive control algorithm has been proposed to solve the control of a twin-motor cable driven system. The special feature of this method is that it uses the multi-time scale characteristics of the system to design a composite control law, which can be divided into fast and slow parts. In addition, to prevent instability of estimated parameters drift, a σ -modification was applied to the adaptive law. Then, stability was proved by the Lyapunov stability theorem. Results show that the closed-loop of the flexible system is stable, and all errors converge into a residual set. One research prototype of the twin-motor cable driven system has been developed, and experimental investigation was conducted on both the proposed σ -modified adaptive control method and conventional adaptive control method. The comparison of experimental results show that the proposed σ -modified adaptive control method can improve the tracking performance well, especially for high-frequency motion control, the fast control law compensation effectively reduces the risk of internal tension instability and enables the system to be stable with large control gains. The proposed σ -modified adaptive control algorithm successfully increases the system bandwidth from 5 Hz in the traditional adaptive control algorithm to 9 Hz.

Author Contributions

Conceptualization, B.L. and L.Y.; methodology, B.L. and L.Y.; software, B.L.; validation, B.L. and L.Y.; formal analysis, B.L.; investigation, B.L.; resources, B.L.; data curation, B.L.; writing—original draft preparation, B.L.; writing—review and editing, B.L., L.Y. and C.G.; visualization, B.L.; supervision, L.Y.; project administration, L.Y.; funding acquisition, L.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (NSFC) under grant 51875013, the Major Project of Ningbo Science and Technology Innovation 2025 Program under grant 2019B10071, the Major Project of the New Generation of Artificial Intelligence under grant 2018AAA0102900, the National Key R&D Program of China under grant 2017YFB1300400, and Fundamental Research Funds for the Central Universities.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request due to restrictions eg privacy or ethical.

Acknowledgments

The authors are thankful for the support of the Science and Technology on Aircraft Control Laboratory and Research Institute for Frontier Science, Beihang University, Beijing 100191, China.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
σ -CACcomposite adaptive control with σ -modification
TACtraditional adaptive control

Appendix A

Recall (11), the equation can be rewritten as
r m M q ¨ + r m B p q ˙ + r m F sign ( q ˙ ) + r m G ( q ) = r m J q T + J q τ * J q τ * .
Then combining (4), (32) and (A1), we get
M e ( q ) q ¨ + N e ( q , q ˙ ) q ˙ + B e q + F e sign ( q ˙ ) + G e ( q ) = r m J q T + J q τ s J q τ * .
Substituting the slow control law (22) into (A2) gives
M e ( q ) q ˙ r + [ N e ( q , q ˙ ) + k d ] q r = Φ T ( q , q ˙ , q ˙ d , q ¨ d ) θ ˜ + r m J q T J q τ * .
Then replace the variables in (A3) with the fast and slow variables defined in (30) and (31) to get
z ˙ q ˙ r = k a 1 0 M e ( q ) 1 [ N e ( q , q ˙ ) + k b ] z q r + 0 1 × 4 M e ( q ) 1 Φ T θ ˜
+ 0 0 M e ( q ) 1 J q 0 r m T ε r m T ˙ + 0 M e ( q ) 1 J q τ * ,
x ˙ = A 1 x + Ψ θ ˜ + A 3 w + B 1 τ * .
Substituting (30) and (31) into (19) gives
ε r m T ˙ ε r m T ¨ = 0 E k 1 r m 2 I m 1 k 2 r m 2 I m 1 r m T ε r m T ˙ + 0 k 1 r m 2 I m 1 τ * ,
ε w ˙ = A 2 w + B 2 τ * .

Appendix B

Recall the τ * defined in (32) and expand it as
τ * = τ s r m 1 I m ( J ˙ q T q ˙ + J q T q ¨ ) .
Combining (11) and (A8) gives
τ * + [ r m 2 M 1 I m J q T J q , 0 ] w = τ * + [ r m 2 M 1 I m J q T J q , 0 ] ( y A 2 1 B 2 τ * ) = τ s r m 1 I m J ˙ q T q ˙ + r m 1 I m J q T M 1 B p q ˙ + F sign ( q ˙ ) + G ( q ) .
Then define
Π ( q ) = E 2 + r m 2 M 1 I m J q T J q ,
Equation (A9) can be rewritten as
Π ( q ) τ * = τ s r m 1 I m J ˙ q T q ˙ + r m 1 I m J q T M 1 B p q ˙ + F sign ( q ˙ ) + G ( q ) [ r m 2 M 1 I m J q T J q , 0 ] y .
Substituting the slow control law (22) into (A11) gives
Π ( q ) τ * = J q + [ M ^ e ( q ) q a + N ^ e ( q , q ˙ ) q v + B ^ e q ˙ + F ^ e sign ( q ˙ ) + G ^ e ( q ) k b q r ] + ( E 2 J q + J q ) ξ r m 1 I m J ˙ q T q ˙ + r m 1 I m J q T M 1 B p q ˙ + F sign ( q ˙ ) + G ( q ) [ r m 2 M 1 I m J q T J q , 0 ] y ,
that is,
τ * = τ * ( x , y , θ ˜ , q d , q ˙ d , q ¨ d ) .
Hence the time derivative of τ * is
τ ˙ * = τ * x x ˙ + τ * y y ˙ + τ * θ ˜ θ ˜ ˙ + ρ ( t ) ,
where
ρ ( t ) = τ * q d q d ˙ + τ * q d ˙ q d ¨ + τ * q d ¨ q d .
Combining (A14) with (38) gives
τ ˙ * = Π τ * x ( A 1 x + Ψ θ ˜ + A 3 y ) + Π 1 ε τ * y A 2 y + Π τ * θ ˜ ( Γ Ξ x σ Γ θ ^ ) + Π ρ ( t ) .
From (A16), there exist corollaries satisfying
  • Corollary 1
Π τ * x A 3 y + Π 1 ε τ * y A 2 y κ 3 + 1 ε κ 2 y ,
where
κ 2 = sup Π τ * x A 2 = sup [ r m 2 M 1 I m J q T J q , 0 ] A 2 , κ 3 = sup Π τ * x A 3 .
  • Corollary 2
( Π τ * x A 1 Π τ * θ ˜ Γ Ξ ) x Π τ * θ ˜ σ Γ θ ^ + Π τ * x Ψ θ ˜ κ 1 x + σ κ 6 ( t ) + κ 7 ( t ) + κ 8 ( t ) ,
Π τ * x Ψ θ ˜ = Π τ * x 0 1 × r M e ( q ) 1 Φ θ ˜ κ 1 , 3 x + κ 7 ( t ) + κ 8 ( t ) ,
where
κ 1 , 1 = sup Π τ * x A 1 , κ 1 , 2 = sup Π τ * θ ˜ Γ Ξ , κ 1 , 3 = sup Π τ * x 0 1 × r M e ( q ) 1 Y 1 ( q , q ˙ ) , κ 1 = κ 1 , 1 + κ 1 , 2 + κ 1 , 3 , κ 6 ( t ) = sup Π τ * θ ˜ Γ θ ^ , κ 7 ( t ) = sup Π τ * x 0 1 × r M e ( q ) 1 Y 2 ( q , q ˙ , q ˙ d , q ¨ d ) , κ 8 ( t ) = sup Π τ * x 0 1 × r M e ( q ) 1 Y 3 ( q , q ˙ ) , Y 1 ( q , q ˙ ) = [ M ˜ e ( q ) k a 2 N ˜ e ( q , q ˙ ) k a + B ˜ e , M ˜ e ( q ) k a ] , Y 2 ( q , q ˙ , q ˙ d , q ¨ d ) = M ˜ e ( q ) q ¨ d + N ˜ e ( q , q ˙ ) q ˙ d + B ˜ e q ˙ d , Y 3 ( q , q ˙ ) = F ˜ e sign ( q ˙ ) + G ˜ e ( q ) .
  • Corollary 3
κ 9 ( t ) = sup Π ρ ( t ) .
Hence, combining (A17), (A19) and (A22) gives
τ ˙ * κ 3 + 1 ε κ 2 y + κ 1 x + κ 4 ( t ) ,
where,
κ 4 ( t ) = σ κ 6 ( t ) + κ 7 ( t ) + κ 8 ( t ) + κ 9 ( t ) .

References

  1. Shiang, W.J.; Cannon, D.; Gorman, J. Dynamic analysis of the cable array robotic crane. In Proceedings of the 1999 IEEE International Conference on Robotics and Automation, Detroit, MI, USA, 10–15 May 1999; pp. 2495–2500. [Google Scholar] [CrossRef]
  2. Vrabel, R. Design of the state feedback-based feed-forward controller asymptotically stabilizing the double-pendulum-type overhead cranes with time-varying hoisting rope length. Int. J. Nonlinear Sci. Numer. Simul. 2020, 1–12. [Google Scholar] [CrossRef]
  3. Nan, R.; Li, D.; Jin, C.; Wang, Q.; Zhu, L.; Zhu, W.; Zhang, H.; Yue, Y.; Qian, L. The five-hundred-meter aperture spherical radio telescope (FAST) project. Int. J. Mod. Phys. D 2011, 20, 989–1024. [Google Scholar] [CrossRef] [Green Version]
  4. Pott, A.; Mütherich, H.; Kraus, W.; Schmidt, V.; Miermeister, P.; Verl, A. IPAnema: A family of Cable-Driven Parallel Robots for Industrial Applications. In Cable-Driven Parallel Robots, 1st ed.; Bruckmann, T., Pott, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2013; Volume 12, pp. 119–134. [Google Scholar]
  5. Wang, Y.; Song, C.; Zheng, T.; Lau, D.; Yang, K.; Yang, G. Cable Routing Design and Performance Evaluation for Multi-link Cable-Driven Robots with Minimal Number of Actuating Cables. IEEE Access 2019, 7, 135790–135800. [Google Scholar] [CrossRef]
  6. Lachele, P.M.; Boss, R. The cablerobot simulator large scale motion platform based on cable robot technology. In Proceedings of the 2016 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Daejeon, Korea, 9–14 October 2016; pp. 3024–3029. [Google Scholar] [CrossRef]
  7. Tadokoro, S.; Murao, Y.; Hiller, M.; Murata, R.; Matsushima, T. A motion base with 6-DOF by parallel cable drive architecture. IEEE-ASME Trans. Mechatron. 2002, 7, 115–123. [Google Scholar] [CrossRef]
  8. Lee, T.; Kim, I.; Baek, Y.S. Design of a 2DoF Ankle Exoskeleton with a Polycentric Structure and a Bi-Directional Tendon-Driven Actuator Controlled Using a PID Neural Network. Actuators 2021, 10, 9. [Google Scholar] [CrossRef]
  9. Xiao, Y.; Lin, Q.; Zheng, Y.; Liang, B. Model Aerodynamic Tests with a Wire-driven Parallel Suspension System in Low-speed Wind Tunnel. Chin. J. Aeronaut. 2010, 23, 393–400. [Google Scholar] [CrossRef] [Green Version]
  10. Qian, S.; Zi, B.; Wang, D.; Li, Y. Development of Modular Cable-Driven Parallel Robotic Systems. IEEE Access 2018, 7, 5541–5553. [Google Scholar] [CrossRef]
  11. Xiong, H.; Diao, X. Cable tension control of cable-driven parallel manipulators with position-controlling actuators. In Proceedings of the 2017 IEEE International Conference on Robotics and Biomimetics (ROBIO), Macau, China, 5–8 December 2017; pp. 1763–1768. [Google Scholar] [CrossRef]
  12. Caverly, R.J.; Forbes, J.R. Flexible Cable-Driven Parallel Manipulator Control: Maintaining Positive Cable Tensions. IEEE Trans. Control Syst. Technol. 2018, 26, 1874–1883. [Google Scholar] [CrossRef]
  13. Oh, S.R.; Agrawal, S.K. Cable suspended planar robots with redundant cables: Controllers with positive tensions. IEEE Trans. Robot. 2005, 21, 457–465. [Google Scholar] [CrossRef]
  14. Shang, W.; Cong, S. Nonlinear computed torque control for a high-speed planar parallel manipulator. Mechatronics 2009, 19, 987–992. [Google Scholar] [CrossRef]
  15. Shang, W.; Zhang, B.; Fei, Z.; Shuang, C. Synchronization Control in the Cable Space for Cable-Driven Parallel Robots. IEEE Trans. Ind. Electron. 2019, 66, 4544–4554. [Google Scholar] [CrossRef]
  16. Najafi, F.; Bakhshizadeh, M. Development a fuzzy PID controller for a parallel cable robot with flexible cables. In Proceedings of the 2016 4th International Conference on Robotics and Mechatronics (ICROM), Tehran, Iran, 26–28 October 2016; pp. 90–97. [Google Scholar] [CrossRef]
  17. Jia, H.; Shang, W.; Xie, F.; Zhang, B.; Cong, S. Second-Order Sliding-Mode-Based Synchronization Control of Cable-Driven Parallel Robots. IEEE-ASME Trans. Mechatron. 2019, 25, 383–394. [Google Scholar] [CrossRef]
  18. Li, B.; Yan, L.; Zhang, L.; Gerada, C. Adaptive Robust Control of the Cable Driven System for Position Tracking. In Proceedings of the 15th IEEE Conference on Industrial Electronics and Applications (ICIEA), Kristiansand, Norway, 9–13 November 2020; pp. 253–257. [Google Scholar] [CrossRef]
  19. Wang, Y.; Yan, F.; Zhu, K.; Chen, B.; Wu, H. Practical Adaptive Integral Terminal Sliding Mode Control for Cable-Driven Manipulators. IEEE Access 2020, 6, 78575–78586. [Google Scholar] [CrossRef]
  20. Wang, Y.; Yan, F.; Chen, J.; Chen, B. Continuous Nonsingular Fast Terminal Sliding Mode Control of Cable-Driven Manipulators with Super-Twisting Algorithm. IEEE Access 2018, 6, 49626–49636. [Google Scholar] [CrossRef]
  21. El-Ghazaly, G.; Gouttefarde, M.; Creuze, V. Adaptive Terminal Sliding Mode Control of a Redundantly-Actuated Cable-Driven Parallel Manipulator: CoGiRo. In Cable-Driven Parallel Robots, 1st ed.; Bruckmann, T., Pott, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2015; Volume 32, pp. 179–200. [Google Scholar]
  22. Kino, H.; Yahiro, T.; Takemura, F.; Morizono, T. Robust PD Control Using Adaptive Compensation for Completely Restrained Parallel-Wire Driven Robots: Translational Systems Using the Minimum Number of Wires Under Zero-Gravity Condition. IEEE Trans. Robot. 2007, 23, 803–812. [Google Scholar] [CrossRef]
  23. Jabbari Asl, H.; Yoon, J. Robust trajectory tracking control of cable-driven parallel robots. Nonlinear Dyn. 2017, 89, 2769–2784. [Google Scholar] [CrossRef]
  24. Alp, A.B.; Alp, S.K. Cable suspended robots: Feedback controllers with positive inputs. In Proceedings of the 2002 American Control Conference, Anchorage, AK, USA, 8–10 May 2002; pp. 815–820. [Google Scholar] [CrossRef]
  25. Shang, W.; Xie, F.; Zhang, H.; Cong, S.; Li, Z. Adaptive Cross-Coupled Control of Cable-Driven Parallel Robots with Model Uncertainties. IEEE Robot. Autom. Lett. 2020, 5, 4110–4117. [Google Scholar] [CrossRef]
  26. Lamaury, J.; Gouttefarde, M.; Chemori, A.; Herve, P.-E. Dual-space adaptive control of redundantly actuated cable-driven parallel robots. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Tokyo, Japan, 3–7 November 2013; pp. 4879–4886. [Google Scholar] [CrossRef] [Green Version]
  27. Caverly, R.J.; Forbes, J.R. Dynamic Modeling and Noncollocated Control of a Flexible Planar Cable-Driven Manipulator. IEEE Trans. Robot. 2014, 30, 1386–1397. [Google Scholar] [CrossRef]
  28. Khosravi, M.A.; Taghirad, H.D. Robust PID control of fully-constrained cable driven parallel robots. Mechatronics 2013, 25, 87–97. [Google Scholar] [CrossRef]
  29. Khosravi, M.A.; Taghirad, H.D. Experimental performance of robust PID controller on a planar cable robot. In Cable-Driven Parallel Robots, 1st ed.; Bruckmann, T., Pott, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2013; Volume 12, pp. 337–352. [Google Scholar]
  30. Babaghasabha, R.; Khosravi, M.A.; Taghirad, H.D. Adaptive robust control of fully-constrained cable driven parallel robots. Mechatronics 2015, 25, 27–36. [Google Scholar] [CrossRef]
  31. Khosravi, M.A.; Taghirad, H.D. Dynamic modeling and control of parallel robots with elastic cables: Singular perturbation approach. IEEE Trans. Robot. 2014, 30, 694–704. [Google Scholar] [CrossRef]
  32. Vafaei, A.; Khosravi, M.A.; Taghirad, H.D. Modeling and Control of Cable Driven Parallel Manipulators with Elastic Cables: Singular Perturbation Theory. In Proceedings of the 4th International Intelligent Robotics and Applications, Aachen, Germany, 6–8 December 2011; pp. 455–464. [Google Scholar]
Figure 1. The schematics of the twin-motor cable driven system.
Figure 1. The schematics of the twin-motor cable driven system.
Actuators 10 00045 g001
Figure 2. Force exerted on the moving platform of the twin-motor cable driven system.
Figure 2. Force exerted on the moving platform of the twin-motor cable driven system.
Actuators 10 00045 g002
Figure 3. Control diagram of the composite adaptive controller.
Figure 3. Control diagram of the composite adaptive controller.
Actuators 10 00045 g003
Figure 4. Component of the flexible cable driven system.
Figure 4. Component of the flexible cable driven system.
Actuators 10 00045 g004
Figure 5. Framework of the control system.
Figure 5. Framework of the control system.
Actuators 10 00045 g005
Figure 6. Comparison results at 3 Hz. (a) Trajectory tracking of σ -CAC; (b) position error of σ -CAC; (c) trajectory tracking of TAC; (d) position error of traditional adaptive control (TAC).
Figure 6. Comparison results at 3 Hz. (a) Trajectory tracking of σ -CAC; (b) position error of σ -CAC; (c) trajectory tracking of TAC; (d) position error of traditional adaptive control (TAC).
Actuators 10 00045 g006
Figure 7. Comparison results at 5 Hz. (a) Trajectory tracking of σ -CAC; (b) position error of σ -CAC; (c) trajectory tracking of TAC; (d) position error of TAC.
Figure 7. Comparison results at 5 Hz. (a) Trajectory tracking of σ -CAC; (b) position error of σ -CAC; (c) trajectory tracking of TAC; (d) position error of TAC.
Actuators 10 00045 g007
Figure 8. Comparison results at 7 Hz. (a) Trajectory tracking of σ -CAC; (b) position error of σ -CAC; (c) trajectory tracking of TAC; (d) position error of TAC.
Figure 8. Comparison results at 7 Hz. (a) Trajectory tracking of σ -CAC; (b) position error of σ -CAC; (c) trajectory tracking of TAC; (d) position error of TAC.
Actuators 10 00045 g008
Figure 9. Results of σ -CAC at 9 Hz. (a) Trajectory tracking of σ -CAC; (b) position error of σ -CAC.
Figure 9. Results of σ -CAC at 9 Hz. (a) Trajectory tracking of σ -CAC; (b) position error of σ -CAC.
Actuators 10 00045 g009
Table 1. Parameters of two controllers.
Table 1. Parameters of two controllers.
SymbolSpecificationσ-CACTAC
M ^ 0 Initial Value of M ^ 20 kg20 kg
B ^ p 0 Initial Value of B ^ p 10 N/m/s10 N/m/s
F ^ 0 Initial Value of F ^ 100 N100 N
G ^ 0 Initial Value of G ^ 200 N200 N
k Γ 1 Gain of M ^ ˙ 0.20.2
k Γ 2 Gain of B ^ ˙ p 11
k Γ 3 Gain of F ^ ˙ 0.40.4
k Γ 4 Gain of G ^ ˙ 0.20.2
k a Feedback gain1515
k b Feedback gain100100
k 1 Coefficient10
k 2 Coefficient10
σ Constant positive value0.1
Table 2. Comparison results of the two control methods.
Table 2. Comparison results of the two control methods.
Frequency (Hz) σ -CACTAC
MagnitudePhase (deg)MagnitudePhase (deg)
23.65%−2.4−0.72%3.6
30.39%−3.24−1.38%−4.32
4−2.09%−1.441.09%−6.91
5−1.39%−2.16−6.54%−9.1
6−0.99%0.864−4.64%−12.96
7−1.17%−1.8−2.55%−16.13
8−2.58%−5.78--
9−2.01%−9.2--
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, B.; Yan, L.; Gerada, C. The Novel Singular-Perturbation-Based Adaptive Control with σ-Modification for Cable Driven System. Actuators 2021, 10, 45. https://0-doi-org.brum.beds.ac.uk/10.3390/act10030045

AMA Style

Li B, Yan L, Gerada C. The Novel Singular-Perturbation-Based Adaptive Control with σ-Modification for Cable Driven System. Actuators. 2021; 10(3):45. https://0-doi-org.brum.beds.ac.uk/10.3390/act10030045

Chicago/Turabian Style

Li, Bin, Liang Yan, and Chris Gerada. 2021. "The Novel Singular-Perturbation-Based Adaptive Control with σ-Modification for Cable Driven System" Actuators 10, no. 3: 45. https://0-doi-org.brum.beds.ac.uk/10.3390/act10030045

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop