Next Article in Journal
Design of an Acoustic Synthetic Jet Actuator for Flow Control
Next Article in Special Issue
Target Tracking of Snake Robot with Double-Sine Serpentine Gait Based on Adaptive Sliding Mode Control
Previous Article in Journal
An Arbitrary Lagrangian Eulerian Approach for Estimating Energy Dissipation in Micromirrors
Previous Article in Special Issue
Implementation of Iterative Learning Control on a Pneumatic Actuator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive Dynamic Surface Technology-Based Electromechanical Actuator Fault-Tolerant Scheme for Blair Mine Hoist Wire Rope Tension Control System

1
State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mines, Anhui University of Science and Technology, Huainan 232001, China
2
Coal Mine Safety Mining Equipment Innovation Center of Anhui Province, Anhui University of Science and Technology, Huainan 232001, China
3
Post-Doctoral Research Station, Huaibei Mining Corporation Limited, Huaibei 235099, China
4
School of Mechatronic Engineering, China University of Mining and Technology, Xuzhou 221116, China
5
College of Mechanical and Electronic Engineering, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Actuators 2022, 11(10), 299; https://doi.org/10.3390/act11100299
Submission received: 22 September 2022 / Revised: 7 October 2022 / Accepted: 15 October 2022 / Published: 18 October 2022
(This article belongs to the Special Issue Applications of Intelligent Control in Actuators Systems)

Abstract

:
As an important equipment for deep well hoisting, the safe and stable operation of the Blair mine hoist is vital for the development and utilization of deep mineral resources. However, it is always a challenging task to keep consistent wire rope tension in the event of an actuator fault. In this study, an adaptive dynamic surface technology-based actuator fault-tolerant scheme is proposed. A fault observer with a neural network adaptation term is designed to estimate the loss of actuator efficiency caused by faults. Considering the redundant characteristic of the two actuators, a novel dynamic surface technology-based controller with a fuzzy assignment and state constraints is developed to eliminate the impact of fault. The stability of the closed-loop system under the proposed strategy is theoretically proved by rigorous Lyapunov analysis. Comparative experiments under various conditions are carried out on a xPC based mine hoist platform, and the results show the applicability together with the superiority of the proposed scheme.

1. Introduction

As the throat equipment connecting the underground and the ground of the ultra-thousand-meter-deep mine, the efficient and reliable operation of the Blair mine hoist is directly related to the life safety of the miners and the development of deep mineral resources. Although the safety level of the mine hoist has been continuously improved with the progress of research, safety risks still exist due to the harsh working environment of the equipment [1]. One of the important problems is that the wire rope tension of the mine hoist will fluctuate greatly under the action of multiple factors such as the partial load of the skip [2], the guide rail friction [3], and the inconsistent winding [4], which seriously affects the service life of the wire rope [5]. According to coal mine safety regulations in China, the difference between the tension of any hoisting wire rope and the average tension shall not exceed ±10%. Therefore, real-time measurement and control of the wire rope tension is necessary to ensure the safe operation of the hoist.
To prevent accidents, a wire rope tension control method has always received extensive attention from scholars. Wang et al. [6] presented a backstepping-based boundary control law to stabilize the mine hoist wire rope tension oscillation state. Zhu et al. [7] designed a nonlinear disturbance observer-based adaptive scheme to improve the performance of the wire rope tension coordination controller. In the literature [8], an intelligent learning control method is used to enhance the tension control effect. Li et al. [9] proposed an approach consisting of the disturbance observer and flatness-based controller to overcome the influence of system uncertainty on wire rope tension. By summing up the above research results, it can be found that nonlinear factors, including unknown disturbance [10], sensor delays [11], and model uncertainties [12], are the main consideration in the design process of the current tension coordination controllers. In addition, it should be pointed out that not only is the bounded interference quite important in engineering applications, but also handling the fault/failure in the system component is necessary for reliability [13]. Otherwise, the goal of tension coordination control cannot be achieved, and even a system crash can be caused. Fortunately, the fault-tolerant technique, which has been widely used in industrial equipment, such as the spacecraft [14], unmanned marine vehicles [15], battery assembly system [16], and robotic manipulators [17], provides a new direction for improving the reliability of the mine hoist closed-loop control system.
In general, fault-tolerant control systems can be classified into two main types: the passive fault-tolerant control system (PFTCS) and active fault-tolerant control system (AFTCS) [18,19]. For the PFTCS, controllers typically operate in offline mode under both normal and abnormal conditions, with predefined parameters to overcome the effects of faulty components. In [20], a novel passive fault-tolerant strategy made up of a Kalman filter and a sliding mode controller was applied to accomplish the path-following control of the electric vehicle under a steering system fault. A passive fault-tolerant switched state-feedback controller designed for linear multivariable systems with time-varying actuator faults was proposed in [21]. A weighted tube model predictive scheme-based passive fault-tolerant controller was developed in [22] to deal with an actuator efficiency loss fault of hypersonic vehicles. In [23], a hybrid fault-tolerant control strategy consisted of a Takagi-Sugeno fuzzy model, and a sliding mode control scheme was proposed to handle unmanned marine vehicles’ dynamic positioning problem in the case of a time-varying delay and thruster faults. An adaptive second-order fault-tolerant sliding mode controller with a proportional-integral-differential sliding surface was designed in [24] to achieve the speed control of an electromechanical system under a structural damage fault condition. As for the AFTCS, an online detection subsystem provides fault conditions constantly, which will guide the supervisory subsystem to adjust the structure and parameters of the controller to cut down faults’ effect on the whole system [25,26]. In [27], a DK-iteration approach-based multi-controller switching strategy was proposed to achieve an active fault-tolerant purpose for unmanned aerial vehicles. A novel AFTC consisting of an artificial neural-network-based fault detection and isolation subsystem and a nonlinear dynamic inversion controller with a new feedback structure for flight was presented in [28]. In [29], the design problem of the actuator fault-tolerant control for a class of two coupling permanent magnet synchronous motors system was discussed. An adaptive robust observer was built to estimate the actuator faults. An improved third-order sliding mode observer was developed in [30] for robot manipulator systems to heighten the effect of the nonsingular fast terminal sliding mode-based fault-tolerant controller. In [31], a dynamic proportional-integral observer-based nonlinear fault-tolerant controller was put forward to stabilize a nonlinear system with a partially unknown dynamic. A neuroadaptive observer and command filtered technique-based fault-tolerant controller was applied to overcome the position tracking problem for induction motors with possible actuator faults and load disturbances [32]. In [33], nonlinear disturbance observers, an asymmetric barrier Lyapunov function, and an adaptive compensation technique were combined to deal with the partial loss of effectiveness problem of the electro-hydraulic actuator system. An iterative learning algorithm-based fault-tolerant controller was established for repetitive systems with Brownian motion in [34]. To our best knowledge of relevant research, there are few existing works taking the fault-tolerant control problem of the mine hoist into account. In fact, electromechanical actuator failure will increase the risk of system oscillation as large flexible wire ropes are used in the mine hoist to carry the skip, which makes a fault-tolerant design of the tension coordination control system of critical importance [35]. On the other hand, the existing fault-tolerant schemes designed for other equipment, such as the H technique [36], sliding mode approach [37,38], adaptive method [39,40], model predictive control scheme [41], etc., are very difficult to be applied on the tension coordination control system without modifications owing to its special nonlinear structure.
Motivated by the above research, in this work, the wire rope tension coordination control problem for the mine hoist in the presence of actuator faults is considered. To achieve efficient actuator fault detection, a neural-network-based adaptive observer is designed and constructed. The controller thus can obtain more efficient fault information. In order to drive the wire rope tension difference into the allowable range under an actuator fault condition, an adaptive dynamic surface technology-based fault-tolerant control scheme is devised. In the control scheme, taking into account the redundant characteristic of the actuators, fuzzy rules are used to construct the control variable allocation factor. Furthermore, to improve control performance, the barrier function-based control law is proposed. The presented controller can not only realize coordination control of wire rope tension under normal working conditions but also maintain acceptable performance under an actuator fault condition. The main contributions of this paper can be summarized as: (1) Considering the Blair mine hoist wire rope tension control system with unknown disturbances, modelling errors and actuator faults, a novel adaptive dynamic surface technology-based fault-tolerant controller with a fuzzy assignment and state constraints, instead of the traditional dynamic surface controller design method, is developed, and the stability of the proposed strategy is mathematically proved; (2) The proposed design does not need to distinguish the type and specific location of actuator failures accurately, and determines the real-time reallocation of actuators with the fuzzy efficiency loss function. This helps to realize fault-tolerant control quickly, even if one side of the actuator is completely stuck which is a challenge for other controllers; and (3) The effectiveness of the proposed controller is verified by the experiments’ results. Unlike existing studies, relevant theories proposed in this paper are proved feasible in the experimental environment, which has strong engineering guiding significance.
This paper is organized as follows: Section 2 provides the nonlinear equation of the wire rope tension coordination control system with an actuator fault. In Section 3, the design of the proposed fault observer and fault-tolerant controller is illustrated. Besides, the stability proof of the overall closed-loop control system is proposed. In Section 4, the experimental design and results’ analysis are presented and discussed. Conclusions are given in Section 5.

2. Mathematical Model of Blair Mine Hoist

2.1. Wire Rope Tension Control Model

A typical Blair mine hoist with an electric cylinder tension adjustment device is shown in Figure 1. It can be seen that under the winding action of the drum, the skip is driven up or down by the vertical rope to realize the transportation of coal. Obviously, if the length difference appears between the left side of the wire rope and the right side, it is easy to cause overloading of one side of the wire rope as the skip withstands the constraint of the guide rail, which is the root cause of the requirement of tension coordination control. At this time, the active adjustment actuator is installed under the headgear sheave to compensate the rope length difference. For the design of the actuator, there are two forms, the hydraulic drive or electric drive. This paper selects the electric drive form for research.
As shown in Figure 2, the mine hoist dynamic model can be divided into three sub-models: the hoisting subsystem model, skip subsystem model, and electric drive actuator model.
In the actual operating conditions, the dynamic behaviors of the mine hoist caused by the interaction between the drum, the wire rope, and the headgear sheave are very complex [42]. Considering the complexity, accuracy, and calculation amount of the control model, the following assumptions are made: (1) The structure of the wire rope is simplified as an elastic wire with a distributed mass; (2) The elastic modulus of the catenaries and the vertical ropes is constant along its entire length; that is, its deformation characteristics conform to Hooke’s law; (3) Compared with the length of the vertical ropes, the lateral displacement of the ends of the vertical ropes on both sides is very small, and the lateral inclination angle of the vertical ropes on both sides can be ignored. Choose the state vector of the hoisting subsystem as
x h s = [ θ w d θ s ( 1 ) θ s ( 2 ) d s ( 1 ) d s ( 2 ) d c n ( 1 ) d c n ( 2 ) ] T
where θ w d is rotation angle of the drum; θ s ( 1 ) and θ s ( 2 ) represent the rotation angle of the left and right headgear sheave, respectively. d s ( 1 ) and d s ( 2 ) denote the vertical displacement of the left and right headgear sheave. d c n ( 1 ) and d c n ( 2 ) indicate displacement of the end of the left and right vertical rope.
According to the above assumption, the displacement of any point on the catenary can be expressed as
x y c ( i ) = θ s ( i ) r s ( i ) d s ( i ) sin θ q x ( i ) + ( l c ( i ) y c ( i ) ) ( l c ( i ) ) 1 ( θ w d r w d ( i ) θ s ( i ) r s ( i ) + d s ( i ) sin θ q x ( i ) )
where r w d ( 1 ) and r w d ( 2 ) denote the radius of the left and right part of the drum. The radii of the left and right headgear sheave are represented by r s ( 1 ) and r s ( 2 ) . The inclination angles of the left and right catenary are denoted by θ q x ( 1 ) and θ q x ( 2 ) . The real-time lengths of the left and right catenary are l c ( 1 ) and l c ( 2 ) . For the point on the left catenary, y c ( 1 ) indicates the displacement between the point and the left rope release point (C point). In addition, for the point on the right catenary, y c ( 2 ) denotes the displacement between the point and the right rope release point (D point).
The displacement of any point on the vertical rope can be represented as
x y v ( i ) = d c n ( i ) + ( l v ( i ) y v ( i ) ) ( l v ( i ) ) 1 ( θ s ( i ) r s ( i ) + d s ( i ) d c n ( i ) )
where the real-time length of the left and right vertical rope are l v ( 1 ) and l v ( 2 ) . For the point on the left vertical rope, y v ( 1 ) is the displacement between the point and the upper point (A point) of the vertical rope. In addition, for the point on the right vertical rope, y v ( 2 ) represents the displacement between the point and the right rope release point (B point).
Analyze the kinetic, potential, and dissipation energy of the hoisting subsystem; Appendix A presents the detailed calculation process, and then the nonlinear dynamic model of the hoisting subsystem can be represented as
ρ w r l c ( i ) r s ( i ) r w d ( i ) 6 θ ¨ w d + [ J s ( i ) + ρ w r ( r s ( i ) ) 2 ( l c ( i ) + l v ( i ) ) 3 ] θ ¨ s ( i ) k v ( i ) r s ( i ) d c n ( i ) + ρ w r l v ( i ) r s ( i ) 6 d ¨ c n ( i ) c c ( i ) r s ( i ) r w d ( i ) θ ˙ w d + ( c c ( i ) + c v ( i ) ) ( r s ( i ) ) 2 θ ˙ s ( i ) + ( c v ( i ) c c ( i ) sin θ q x ( i ) ) r s ( i ) d ˙ s ( i ) c v ( i ) r s ( i ) d ˙ c n ( i ) k c ( i ) r s ( i ) r w d ( i ) θ w d + ( k c ( i ) + k v ( i ) ) ( r s ( i ) ) 2 θ s ( i ) + ( k v ( i ) k c ( i ) sin θ q x ( i ) ) r s ( i ) d s ( i ) + ( ρ w r l v ( i ) r s ( i ) 3 ρ w r l c ( i ) r s ( i ) 3 sin θ q x ( i ) ) d ¨ s ( i ) = ρ w r g l c ( i ) r s ( i ) sin θ q x ( i ) 2 ρ w r g l v ( i ) r s ( i ) 2
ρ w r l c ( i ) r w d ( i ) sin θ q x ( i ) 6 θ ¨ w d + ( ρ w r l v ( i ) r s ( i ) 3 ρ w r l c ( i ) r s ( i ) 3 sin θ q x ( i ) ) θ ¨ s ( i ) + ρ w r l v ( i ) 6 d ¨ c n ( i ) + ( m s ( i ) + ρ w r l v ( i ) 3 + ρ w r l c ( i ) 3 sin 2 θ q x ( i ) ) d ¨ s ( i ) + c c ( i ) r w d ( i ) sin θ q x ( i ) θ ˙ w d + ( c v ( i ) c c ( i ) sin θ q x ( i ) ) r s ( i ) θ ˙ s ( i ) + ( c v ( i ) + c c ( i ) sin 2 θ q x ( i ) ) d ˙ s ( i ) c v ( i ) d ˙ c n ( i ) + k c ( i ) r w d ( i ) sin θ q x ( i ) θ w d + ( k v ( i ) k c ( i ) sin θ q x ( i ) ) r s ( i ) θ s ( i ) + ( k v ( i ) + k c ( i ) sin 2 θ q x ( i ) ) d s ( i ) k v ( i ) d c n ( i ) = F e m a ( i ) ρ w r g l v ( i ) 2 m s ( i ) g ρ w r g l c ( i ) r s ( i ) 2 sin θ q x ( i )
ρ w r l v ( i ) r s ( i ) 6 θ ¨ s ( i ) + ρ w r l v ( i ) 6 d ¨ s ( i ) + ( m c n ( i ) + ρ w r l v ( i ) 3 ) d ¨ c n ( i ) c v ( i ) r s ( i ) θ ˙ s ( i ) c v ( i ) d ˙ s ( i ) + c v ( i ) d ˙ c n ( i ) k v ( i ) r s ( i ) θ s ( i ) k v ( i ) d s ( i ) + k v ( i ) d c n ( i ) = F c n ( i ) m c n ( i ) g ρ w r g l v ( i ) 2
where ρ w r indicates the linear density of the wire ropes. J s ( 1 ) and J s ( 2 ) represent the equivalent moment of inertia of the left and right headgear sheave, respectively. k v ( 1 ) = E A / l v ( 1 ) and k v ( 2 ) = E A / l v ( 2 ) denote the stiffness coefficients of the left vertical rope and the right vertical rope which are a function of the rope length. EA is the stiffness constant. The stiffness coefficients of the left catenary and the right catenary are represented by k c ( 1 ) = E A / l c ( 1 ) and k c ( 2 ) = E A / l c ( 2 ) .
Permanent magnet synchronous motors are used to drive the actuator action; the simplified equation of the motor can be represented as
{ u m t d ( i ) = R m o ( i ) i m t d ( i ) + L m t d ( i ) d d t i m t d ( i ) ω e ( i ) L m t q ( i ) i m t q ( i ) u m t q ( i ) = R m o ( i ) i m t q ( i ) + L m t q ( i ) d d t i m t q ( i ) + ω e ( i ) ( L m t d ( i ) i m t d ( i ) + ψ m t f ( i ) ) J h s ( i ) θ ¨ h s ( i ) = T m t ( i ) B h s ( i ) θ ˙ h s ( i ) T h s ( i ) ( i = 1 , 2 )
where θ h s ( 1 ) and θ h s ( 2 ) denote the rotation angles of the left and right actuator-driven motor, respectively. J h s ( 1 ) and J h s ( 2 ) indicate the equivalent moments of inertia of the left and right actuator-driven motor. B h s ( 1 ) and B h s ( 2 ) represent the damping coefficients of the left and right actuator-driven motor. T m t ( 1 ) and T m t ( 2 ) are the electromagnetic torque values of the left and right actuator-driven motor and can be defined as
T m t ( i ) = 1.5 p m t ( i ) i m t q ( i ) [ i m t d ( i ) ( L m t d ( i ) L m t q ( i ) ) + ψ m t f ( i ) ] = 1.5 p m t ( i ) i m t q ( i ) ψ m t f ( i ) + χ m t ( i ) ( i = 1 , 2 )
where p m t ( 1 ) and p m t ( 2 ) indicate the numbers of the left and right actuator-driven motor pole pair, respectively. i m t q ( 1 ) and i m t q ( 2 ) represent the currents of the left and right actuator-driven motor along q axes. i m t d ( 1 ) and i m t d ( 2 ) denote the currents of the left and right actuator-driven motor along d axes. L m t d ( 1 ) and L m t d ( 2 ) are the inductances of the left and right actuator-driven motor along d axes. The inductances of the left and right actuator-driven motor along q axes are represented by L m t q ( 1 ) and L m t q ( 2 ) , respectively. ψ m t f ( 1 ) and ψ m t f ( 2 ) denote the flux linkage coefficients of the left and right actuator-driven motor. As the selected motors are generally the surface-mount type, thus, L m t d ( i ) L m t q ( i ) and the nonlinear term in Equation (8) satisfies χ m t ( i ) 0 .

2.2. Problem Formulation of Actuator Fault

During the long-term service of the mine hoist wire rope tension control system, from the structure of the electric actuator, it can be known that key components such as the motor, gearbox, and screw mechanism may turn out faults [43,44]. The fault tree of electromechanical actuators in the tension coordination control system is shown in Figure 3. The failure event and the corresponding ID are shown in Table 1. It is not difficult to know that each abnormal situation may correspond to a variety of component failures; for instance, electrical failure, structural deformation, and excessive frictional resistance may cause the actuator to be stuck. In addition, with the development of the same type of component failure to different extents, it will lead to different abnormal response situations. For example, in the early stage of screw deformation, the electromechanical actuator may show a response situation with reduced efficiency. When the bending deformation of the screw is large, the actuator will show a stuck response situation.
From a macro perspective, the abnormal response of the actuator subsystem of the tension coordination control system mainly includes several types, such as stuck situations, out-of-control situations, efficiency reduction, random fluctuation, and is reflected in the deviation between the actual output response and the expected response [45,46]. More specifically, the expected electromechanical actuator position refers to the action of the non-faulty actuator under a given input voltage, and the actual electromechanical actuator position refers to the action of the faulty actuator under the given input voltage. From the actual situation, for the electromechanical actuator, the movement resistance anomaly caused by mechanical deformation faults is defined as f a f 1 ( i ) , the unstable control voltage of the motor due to an electronic unit failure as f a f 2 ( i ) , and the transmission coefficient change as f a f 3 ( i ) , and then the mathematical model of the electromechanical actuator with faults can be expressed as
{ J h s ( i ) θ ¨ h s ( i ) = 1.5 p m t ( i ) ( i m t q ( i ) + f a f 2 ( i ) ) ψ m t f ( i ) B h s ( i ) θ ˙ h s ( i ) T h s ( i ) f a f 1 ( i ) d ˙ e m a ( i ) = f a f 3 ( i ) z b c ( i ) z g b ( i ) θ ˙ h s ( i ) F e m a ( i ) = T h s ( i ) / f a f 3 ( i ) z b c ( i ) z g b ( i ) ( i = 1 , 2 )
where z b c ( i ) and z g b ( i ) indicate the transmission ration of the screw i and the gearbox i, respectively.

3. Adaptive Actuator Fault-Tolerant Controller Design

In this section, the dynamic surface technique-based adaptive fault-tolerant scheme will be developed for the wire rope tension control system’s actuator fault. The overall architecture of the control system is shown in Figure 4. In this approach, in light of the fact that the faulty actuator is unknown, adaptive actuator fault observers consisting of the left fault observer and right fault observer are firstly designed for the left actuator and the right actuator, respectively, which are applied to provide actuators’ fault information timely. Through the received fault indicators Δ f a u l t ( i ) and sensor feedback information, utilizing the redundancy properties of the two actuators, the tension coordination control loop, which is the first design step of the controller for adaptive control distribution, and the actuator control loop, which is the second design step of the controller for fast reference tracking, are put forward. In this way, the electromechanical actuator can continue to drive the headgear sheave to run up and down under fault and disturbance conditions finally to achieve the control goal.

3.1. The Adaptive Actuator Fault Observer

For the actual fault detection system and control system, accurate estimation of the specific fault type and amplitude requires a large number of sensor signals, field data, and complex calculation processes, which is not conducive to the timely response and action of the control system. Therefore, for safety reasons, it is more meaningful to detect the abnormal response state and degree of the actuator in time for the actual system.
In this section, a fault observer is designed separately for each motor of electromechanical actuator. Combining the faults of the same channel in Equation (9), the state equation of i-th electromechanical actuator motor can be rewritten as
{ z ˙ f o ( i ) = m f o ( i ) z f o ( i ) + h f o ( i ) u f o ( i ) + G f o ( i ) + β f o ( i ) ( z f o ( i ) ) y f o ( i ) = C f o z f o ( i )
where z f o ( i ) = [ z f o 1 ( i ) z f o 2 ( i ) ] T = [ θ h s ( i ) θ ˙ h s ( i ) ] T (i = 1,2) represent the i-th electromechanical actuator state variables. T represents the standard matrix transposition operation. m f o ( i ) ( z f o ( i ) ) = [ 0 1 0 B h s ( i ) / J h s ( i ) ] . h f o ( i ) = [ 0 1.5 p m t ( i ) ψ m t f ( i ) / J h s ( i ) ] T denote the input coefficient vector. u f o ( i ) = i m t q ( i ) indicate input of the i-th motor. G f o ( i ) ( z f o ( i ) ) = [ 0 T h s ( i ) / J h s ( i ) ] T are load items of the actuator. β f o ( i ) ( z f o ( i ) ) = [ β f o 1 ( i ) β f o 2 ( i ) ] T = [ 0 ( 1.5 p m t ( i ) f a f 2 ( i ) ψ m t f ( i ) f a f 1 ( i ) + ( 1 f a f 3 ( i ) ) F e m a ( i ) z b c ( i ) z g b ( i ) ) / J h s ( i ) ] T represent the sum of the unmodeled states, uncertainties, and faults. C f o = I 2 represent the observation matrix.
Considering the universal approximation characteristics of neural networks and the real-time calculation requirements of the controller, a radial basis neural network structure is used to approximate β f o ( i ) ( z f o ( i ) ) . The structure of the employed neural network is shown in Figure 5.
Thus, we have
β f o ( i ) ( z f o ( i ) ) = ( ψ f o ( i ) ) T R f o ( i ) κ f o ( i )
where ψ f o ( i ) , R f o ( i ) , and κ f o ( i ) indicate the weights, hidden layer output, and approximation error of the radial basis neural network.
Assumption 1.
Comparing the calculated results of the established model with the actual actuator dynamic characteristics, the error is bounded. In addition, for the actual mine hoist, the parameters or coefficient changes caused by faults are bounded due to hardware limitations.
Assumption 2.
Parameters of the radial basis neural network which is used to equivalently indicate the model simplification error and unknown disturbance of electromechanical actuators are bounded.
Assumption 3.
The target failure has the characteristics of slow changes, which means that derivative of  ψ f o ( i ) can be regard as zero.
This study does not make a more detailed distinction between motor failures and electronic control unit failures, and deals with them in a unified manner. On the other hand, the velocity signal is generally obtained by differentiating the rotation angle signal, so the displacement and velocity signals thus can be obtained by observation and calculation.
For i-th actuator motor faults, an adaptive fault observer is designed in the form of
{ z ^ ˙ f o ( i ) = m f o ( i ) z ^ f o ( i ) + h f o ( i ) u f o ( i ) + G f o ( i ) α f o ( i ) e f o ( i ) δ f o ( i ) χ f o ( i ) β ^ f o ( i ) y ^ f o ( i ) = C f o z ^ f o ( i )
where z ^ f o ( i ) = [ z ^ f o 1 ( i ) z ^ f o 2 ( i ) ] T denote the estimated state. e f o ( i ) = [ e f o 1 ( i ) e f o 2 ( i ) ] T = y ^ f o ( i ) y f o ( i ) = C f o e z f o ( i ) represent the error between the observed state and the measured state. α f o ( i ) = diag ( α f o 1 ( i ) , α f o 2 ( i ) ) and δ f o ( i ) = diag ( δ f o 1 ( i ) , δ f o 2 ( i ) ) indicate the coefficient matrix; all elements inside are positive constants. χ f o ( i ) are functions of the observation error and can be designed as
χ f o ( i ) = [ χ f o 1 ( i ) χ f o 2 ( i ) ] T = [ ( | e f o 1 ( i ) | p s o f 1 ( i ) ( p x o f 1 ( i ) ) 1 + θ f o 1 ( i ) ) sign ( e f o 1 ( i ) ) ( | e f o 2 ( i ) | p s o f 2 ( i ) ( p x o f 2 ( i ) ) 1 + θ f o 2 ( i ) ) sign ( e f o 2 ( i ) ) ]
where θ f o 1 ( i ) , θ f o 2 ( i ) , p s o f 1 ( i ) , p x o f 1 ( i ) , p s o f 2 ( i ) , p x o f 2 ( i ) are positive constants.
β ^ f o ( i ) = [ β ^ f o 1 ( i ) β ^ f o 2 ( i ) ] T = ( ψ ^ f o ( i ) ) T R f o ( i ) are the adaptive compensation terms, and the change rate of which can be given as
ψ ^ ˙ f o ( i ) = γ f o ( i ) R f o ( i ) ( e f o ( i ) ) T W f o ( i )
where γ f o ( i ) represent the coefficient matrix.
Theorem 1.
For the electromechanical actuator with faults, design fault observers, as shown in Equation (12), if there exist positive scalars α f o ( i ) matrix W f o ( i ) such that the following inequalities hold:
[ ( m f o ( i ) α f o ( i ) ) T W f o ( i ) + W f o ( i ) ( m f o ( i ) α f o ( i ) ) + I 2 W f o ( i ) ( W f o ( i ) ) T ( σ f o ( i ) ) 2 I 2 ] < 0
with the adaptive compensation term given by (14), then the actuator state estimation error and fault estimation error are both ultimately uniformly bounded.
Proof of Theorem 1.
Following the Lyapunov function, the estimation error is defined as
V f o ( i ) ( t ) = ( e z f o ( i ) ( t ) ) T W f o ( i ) e z f o ( i ) ( t ) + t r [ ( ψ ^ f o ( i ) ψ f o ( i ) ) T ( γ f o ( i ) ) 1 ( ψ ^ f o ( i ) ψ f o ( i ) ) ]
As C f o is the identity matrix, we have e f o ( i ) = e z f o ( i ) , from Equations (14) and (16); according to Assumption 3, derivatives of Lyapunov functions can be expressed as
V ˙ f o ( i ) ( t ) = ( ( e ˙ z f o ( i ) ( t ) ) T W f o ( i ) e z f o ( i ) ( t ) + ( e z f o ( i ) ( t ) ) T W f o ( i ) ( e ˙ z f o ( i ) ( t ) ) ) + t r [ 2 ( ψ ^ f o ( i ) ψ f o ( i ) ) T ( γ f o ( i ) ) 1 ψ ^ ˙ f o ( i ) ] = ( e z f o ( i ) ( t ) ) T [ ( m f o ( i ) α f o ( i ) ) T W f o ( i ) + W f o ( i ) ( m f o ( i ) α f o ( i ) ) ] e z f o ( i ) ( t ) + ( κ f o ( i ) δ f o ( i ) χ f o ( i ) ) T W f o ( i ) e z f o ( i ) ( t ) + ( e z f o ( i ) ( t ) ) T W f o ( i ) ( κ f o ( i ) δ f o ( i ) χ f o ( i ) ) + tr [ 2 ( ψ ^ f o ( i ) ψ f o ( i ) ) T ( γ f o ( i ) ) 1 ( γ f o ( i ) ( R f o ( i ) ( e f o ( i ) ) T W f o ( i ) ) ) ] + 2 ( e z f o ( i ) ( t ) ) T W f o ( i ) ( ψ ^ f o ( i ) ψ f o ( i ) ) T R f o ( i )
As ( e z f o ( i ) ( t ) ) T W f o ( i ) ( ψ ^ f o ( i ) ψ f o ( i ) ) T R f o ( i ) are real numbers and equal to the trace of matrix ( ψ ^ f o ( i ) ψ f o ( i ) ) T ( R f o ( i ) ( e z f o ( i ) ( t ) ) T W f o ( i ) ) , the function can be written as
V ˙ f o ( i ) ( t ) = ( e z f o ( i ) ( t ) ) T O f o ( i ) e z f o ( i ) ( t ) + ( λ f o ( i ) ) T W f o ( i ) e z f o ( i ) ( t ) + ( e z f o ( i ) ( t ) ) T W f o ( i ) λ f o ( i )
where ( m f o ( i ) α f o ( i ) ) T W f o ( i ) + W f o ( i ) ( m f o ( i ) α f o ( i ) ) = O f o ( i ) , ( κ f o ( i ) δ f o ( i ) χ f o ( i ) ) = λ f o ( i ) .
Letting ξ f o ( i ) = [ ( e z f o ( i ) ( t ) ) T ( λ f o ( i ) ) T ] T , it is clear that
V ˙ f o ( i ) ( t ) = ( ξ f o ( i ) ) T [ O f o ( i ) W f o ( i ) ( W f o ( i ) ) T 0 ] ξ f o ( i )
In this case, the indicator H∞ is defined as
0 t ( e z f o ( i ) ( t ) ) T ( e z f o ( i ) ( t ) ) d t < ( σ f o ( i ) ) 2 0 t ( λ f o ( i ) ) T ( λ f o ( i ) ) d t + V ˙ f o ( i ) ( 0 )
where σ f o ( i ) are positive real numbers.
Equation (19) can be reorganized as
V ˙ f o ( i ) ( t ) + ( e z f o ( i ) ( t ) ) T ( e z f o ( i ) ( t ) ) ( σ f o ( i ) ) 2 ( λ f o ( i ) ) T ( λ f o ( i ) ) = ( ξ f o ( i ) ) T [ O f o ( i ) + I 2 W f o ( i ) ( W f o ( i ) ) T ( σ f o ( i ) ) 2 I 2 ] ξ f o ( i ) = ( ξ f o ( i ) ) T Ω f o ( i ) ξ f o ( i )
Letting Ω f o ( i ) < 0 , we have
V ˙ f o ( i ) ( t ) + ( e z f o ( i ) ( t ) ) T ( e z f o ( i ) ( t ) ) ( σ f o ( i ) ) 2 ( λ f o ( i ) ) T ( λ f o ( i ) ) 0
Performing an integration of both sides of the formula (22) and moving items on, we can obtain
V f o ( i ) ( t ) ( σ f o ( i ) ) 2 0 t ( λ f o ( i ) ) T ( λ f o ( i ) ) d t + V f o ( i ) ( 0 ) 0 t ( e z f o ( i ) ( t ) ) T ( e z f o ( i ) ( t ) ) d t ( σ f o ( i ) ) 2 0 t ( λ f o ( i ) ) T ( λ f o ( i ) ) d t + V f o ( i ) ( 0 )
According to Assumption 2, considering the boundary of the approximation error, setting up 0 t ( λ f o ( i ) ) T ( λ f o ( i ) ) d t λ ¯ f o ( i ) , Equation (23) can be written as
V f o ( i ) ( t ) ( σ f o ( i ) ) 2 λ ¯ f o ( i ) + V f o ( i ) ( 0 )
Then, it is clear that
( W f o ( i ) ) min e z f o ( i ) ( t ) 2 ( e z f o ( i ) ( t ) ) T W f o ( i ) e z f o ( i ) ( t ) ( σ f o ( i ) ) 2 λ ¯ f o ( i ) + V f o ( i ) ( 0 )
which means that the state estimation error and the fault estimation error of the motor of the electromechanical actuator are finally unanimous, and convergence results are e z f o ( i ) ( t ) 2 [ ( σ f o ( i ) ) 2 λ ¯ f o ( i ) + V f o ( i ) ( 0 ) ] / ( W f o ( i ) ) min .□

3.2. The Improved Dynamic Surface Technology-Based Fault-Tolerant Controller

Dynamic surface control technology is an important nonlinear control technology developed on the basis of backstepping control technology. Compared with the traditional backstepping control strategy, the dynamic surface control technology can better solve the “item expansion” problem in the backstepping controller design process, and has obvious advantages in the design of high-order system controllers. This section focuses on the fault-tolerance performance improvement of the controller by considering the actuator faults’ condition during the design process.
Define the state variables as z m h = [ z m h 1 z m h 2 ( i ) z m h 3 ( i ) ] T = [ Δ F c n F e m a ( i ) θ ˙ h s ( i ) ] T , ( i = 1 , 2 ) , from Equation (4)–(6). The tension coordination governing equations are defined as
{ z ˙ m h 1 = k e m a [ ( f a f 3 ( 1 ) z b c ( 1 ) z g b ( 1 ) z m h 3 ( 1 ) d ˙ s ( 1 ) ) ( f a f 3 ( 2 ) z b c ( 2 ) z g b ( 2 ) z m h 3 ( 2 ) d ˙ s ( 2 ) ) ] / ( 1 + sin θ q x ) ( 0.5 ρ w r l v ( 1 ) r s ( 1 ) ) θ s ( 1 ) + ( 0.5 ρ w r l v ( 2 ) r s ( 2 ) ) θ s ( 2 ) ( m s ( 1 ) + 0.5 ρ w r l v ( 1 ) ) d s ( 1 ) + ( m s ( 2 ) + 0.5 ρ w r l v ( 2 ) ) d s ( 2 ) ( J s ( 1 ) + 0.5 ρ w r r s ( 1 ) l v ( 1 ) ) θ s ( 1 ) sin θ q x ρ w r Δ l ˙ v g + ( J s ( 2 ) + 0.5 ρ w r r s ( 2 ) l v ( 2 ) ) θ s ( 2 ) sin θ q x + ν t d 1 z ˙ m h 2 ( i ) = k e m a ( f a f 3 ( i ) z b c z g b z m h 3 ( i ) d ˙ s ( i ) ) + ν t d 2 ( i ) z ˙ m h 3 ( i ) = β m h 1 ( i ) u m h ( i ) + ( 1.5 p m t ( i ) f a f 2 ( i ) ψ m t f ( i ) f a f 1 ( i ) + ( 1 f a f 3 ( i ) ) z m h 2 ( i ) z b c ( i ) z g b ( i ) ) / J h s ( i ) β m h 2 ( i ) z m h 3 ( i ) β m h 3 ( i ) z m h 2 ( i ) + ν t d 3 ( i )
where k e m a indicates equivalent stiffness between the electromechanical actuator and the headgear sheave. θ q x denotes the initial inclination angle of the catenary. ν t d 1 , ν t d 2 ( i ) , and ν t d 3 ( i ) represent the sum of the unmodeled states and disturbances in each equation, respectively. β m h 1 ( i ) = 1.5 p m t ( i ) ψ m t f ( i ) / J h s ( i ) , β m h 2 ( i ) = B h s ( i ) / J h s ( i ) , and β m h 3 ( i ) = z b c ( i ) z g b ( i ) signify coefficients of the actuator state equation.
For the subsequent control allocation, the fault indicator Δ f a u l t ( i ) is defined according to the results of the observation
Δ f a u l t ( i ) = 0 t ( α f o 2 ( i ) e f o 2 ( i ) + δ f o 2 ( i ) χ f o 2 ( i ) + β ^ f o 2 ( i ) ) d t
Considering the nonlinear characteristics of the tension coordination control system under a fault condition and the use of the diagnosis of the fault observer, the following introduces the design process of the wire rope tension coordination controller:
Step 1: First, define the tracking deviation of the tension difference as
E m h 1 = z m h 1 F m h r
According to Equation (26), the differential of Equation (28) can be calculated as
E ˙ m h 1 = z ˙ m h 1 F ˙ m h r = k e m a [ ( z b c z g b z m h 3 ( 1 ) d ˙ s ( 1 ) ) ( z b c z g b z m h 3 ( 2 ) d ˙ s ( 2 ) ) ] / ( 1 + sin θ q x ) + f m h 1 F ˙ m h r
where z b c and z g b indicate the design values of the transmission coefficient.
In addition, f m h 1 is expressed as
f m h 1 = k e m a [ f a f 3 ( 1 ) z b c ( 1 ) z g b ( 1 ) z m h 3 ( 1 ) z b c z g b z m h 3 ( 1 ) + z b c z g b z m h 3 ( 2 ) f a f 3 ( 2 ) z b c ( 2 ) z g b ( 2 ) z m h 3 ( 2 ) ] / ( 1 + sin θ q x ) ( 0.5 ρ w r l v ( 1 ) r s ( 1 ) ) θ s ( 1 ) + ( J s ( 2 ) + 0.5 ρ w r r s ( 2 ) l v ( 2 ) ) θ s ( 2 ) sin θ q x + ( 0.5 ρ w r l v ( 2 ) r s ( 2 ) ) θ s ( 2 ) ( m s ( 1 ) + 0.5 ρ w r l v ( 1 ) ) d s ( 1 ) ρ w r Δ l ˙ v g + ( m s ( 2 ) + 0.5 ρ w r l v ( 2 ) ) d s ( 2 ) ( J s ( 1 ) + 0.5 ρ w r r s ( 1 ) l v ( 1 ) ) θ s ( 1 ) sin θ q x + ν t d 1
Remark 1.
On the one hand, the structure of the hoisting system is symmetrical; the corresponding parameters on each sides can be considered consistent, which means that r s ( 1 ) r s ( 2 ) , m s ( 1 ) m s ( 2 ) , J s ( 1 ) J s ( 2 ) , z b c ( 1 ) z b c ( 2 ) , z g b ( 1 ) z g b ( 2 ) . On the other hand, due to the large quality or rotation inertia of mechanical devices, high-level items of displacement or the corner can be regarded as a boundary. The sum of transmission department faults, uncertainties, and model errors thus meet f m h 1 f ¯ m h 1 and can be expressed as
f m h 1 = ( ψ m h 1 ) T R m h 1 κ m h 1
Assumption 4.
Parameters of neural networks for an equivalent representation of modeling uncertainties and the unknown disturbance excitation of the hoisting subsystem are bounded, which means that ψ m h 1 ψ m h 1 M , κ m h 1 κ ¯ m h 1 . κ ¯ m h 1 represents the upper bound of the approximation error of the neural network adaptive term.
Take the first Lyapunov function as V m h 1 = ( E m h 1 ) 2 / 2 to make V ˙ m h 1 = E m h 1 ( z ˙ m h 1 F ˙ m h r ) < 0 ; set z m h 3 ( 1 ) and z m h 3 ( 2 ) as the virtual input, and, to reduce the effect of actuator faults, the fuzzy control allocation factor is constructed in two steps.
(1)
According to the failure status of the hoisting system and the fault indicator Δ f a u l t ( i ) (i = 1, 2) in Equation (27), define fuzzy sets Z n i ( i )
(2)
Utilize the following fuzzy rules’ construct control distribution factors μ m h .
Μ m h ( j ) : IF Δ f a u l t ( 1 ) is Z n 1 ( 1 ) , and Δ f a u l t ( 2 ) is Z n 2 ( 2 ) THEN μ m h is K n 1 , n 2
where K n 1 , n 2 is a conclusion fuzzy set; Μ m h ( j ) denotes the j-th rule, and the total number of rules is Π i = 1 2 ξ a f f p ( i ) ; ξ a f f p ( i ) represents the number of fuzzy sets designed for Δ f a u l t ( i ) .
According to the above rules, the specific calculation process of the distribution factor can be obtained as follows:
(1)
Use the product inference engine to solve the premise inference results of the rules Π i = 1 2 γ Z n i ( i ) ( Δ f a u l t ( i ) ) ;
(2)
Take the central value of K n 1 , n 2 as λ n 1 , n 2 ;
(3)
Let λ n 1 , n 2 be a free parameter and write to the set ψ a f f p using a product inference engine and a central average defuzzifier; the control allocation factor can finally be written as
μ m h = s a t [ ψ a f f p T R a f f p ( Δ f a u l t ( i ) ) ]
where R a f f p ( Δ f a u l t ( i ) ) is Π i = 1 2 ξ a f f p ( i ) , the dimensional vector, whose n 1 n 2 -th element can be expressed as
R a f f p , n 1 n 2 ( i ) ( Δ f a u l t ( i ) ) = Π i = 1 2 γ Z n i ( i ) ( Δ f a u l t ( i ) ) n 1 = 1 ξ a f f p ( 1 ) n 2 = 1 ξ a f f p ( 2 ) ( Π i = 1 2 γ Z n i ( i ) ( Δ f a u l t ( i ) ) )
Remark 2.
When the performance of the actuator is seriously ineffective or even completely locked in place, the fault magnitude may exceed the fixed control gain action in the conventional controller and lead to the effect of a decrease in tension coordination control. Designing the fuzzy factor to achieve adaptive allocation of actuator performance will enable the actuator, which is a low level of failure to take on more functions, and thus bring about a fault-tolerant control purpose.
Making use of distribution factor μ m h obtained from the real-time solution, the reference values of z m h 3 ( i ) are given as
{ z m h d 3 ( 1 ) = μ m h z b c z g b ( d ˙ s ( 1 ) d ˙ s ( 2 ) + 2 α m h 1 [ F ˙ m h r σ m h 1 E m h 1 ( ψ ^ m h 1 ) T R m h 1 η m h 1 sgn ( E m h 1 ) ] ) z m h d 3 ( 2 ) = ( μ m h 1 ) z b c z g b ( d ˙ s ( 1 ) d ˙ s ( 2 ) + 2 α m h 1 [ F ˙ m h r σ m h 1 E m h 1 ( ψ ^ m h 1 ) T R m h 1 η m h 1 sgn ( E m h 1 ) ] )
where σ m h 1 and η m h 1 are positive coefficients. α m h 1 = ( 1 + sin θ q x ) / 2 k e m a . ( ψ ^ m h 1 ) T R m h 1 represent the mechanical transmission fault compensation item, the rate changes of which are defined as follows
ψ ^ ˙ m h 1 = ϖ m h 1 [ E m h 1 R m h 1 ρ m h 1 ψ ^ m h 1 ]
Remark 3.
As the actuator speed reference signal, z m h d 3 ( i ) is brought forward considering the sum of errors due to modeling simplification, random disturbance because of asymmetric drum winding and friction, and inaccurate structure parameters owing to nonlinear characteristics. The traditional robust term needs to overcome the sum of the above factors, which is easy to cause overshoot and oscillation. The introduction of the adaptive term reduces the amplitude that the robust term needs to compensate, and the tension control effect can be improved to some extent.
Define low-pass filter outputs of z m h d 3 ( i ) as z ¯ m h d 3 ( i ) with the following form
{ ϑ m h 2 ( 1 ) z ¯ ˙ m h d 3 ( 1 ) + z ¯ m h d 3 ( 1 ) = z m h d 3 ( 1 ) , z m h d 3 ( 1 ) ( 0 ) = z ¯ m h d 3 ( 1 ) ( 0 ) ϑ m h 2 ( 2 ) z ¯ ˙ m h d 3 ( 2 ) + z ¯ m h d 3 ( 2 ) = z m h d 3 ( 2 ) , z m h d 3 ( 2 ) ( 0 ) = z ¯ m h d 3 ( 2 ) ( 0 )
where ϑ m h 2 ( 1 ) and ϑ m h 2 ( 2 ) are filter coefficients.
Remark 4.
The problem of the differential explosion of the traditional backstepping controller due to repeated differentiation is solved by introducing the filter, which is one of the characteristics of the dynamic surface controller.
Step 2: Taking the low-pass filtered output z ¯ m h d 3 ( 1 ) and z ¯ m h d 3 ( 2 ) as reference values, deviations of the tracking speed of the left and right electromechanical actuator are defined as
E m h 2 ( 1 ) = z m h 3 ( 1 ) z ¯ m h d 3 ( 1 ) , E m h 2 ( 2 ) = z m h 3 ( 2 ) z ¯ m h d 3 ( 2 )
According to (26), it is clear that
E ˙ m h 2 ( i ) = z ˙ m h 3 ( i ) z ¯ ˙ m h d 3 ( i ) = β m h 1 ( i ) u m h ( i ) β m h 2 ( i ) z m h 3 ( i ) β m h 3 ( i ) z m h 2 ( i ) + f m h 2 ( i ) z ¯ ˙ m h d 3 ( i )
where f m h 2 ( i ) is written as
f m h 2 ( i ) = ( 1.5 p m t ( i ) f a f 2 ( i ) ψ m t f ( i ) f a f 1 ( i ) + ( 1 f a f 3 ( i ) ) z m h 2 ( i ) z b c ( i ) z g b ( i ) ) / J h s ( i ) + ν t d 3 ( i )
Taking advantage of the observation results of the previous section, f m h 2 ( i ) can be estimated by the actuator fault observations. Define the estimation error as Δ f m h 2 ( i ) = f m h 2 ( i ) + α f o 2 ( i ) e f o 2 ( i ) + δ f o 2 ( i ) χ f o 2 ( i ) + β ^ f o 2 ( i ) , which can be expressed as
Δ f m h 2 ( i ) = ( ψ m h 2 ( i ) ) T R m h 2 ( i ) κ m h 2 ( i )
where ψ m h 2 ( i ) , R m h 2 ( i ) , and κ m h 2 ( i ) denote the approximation coefficient, ideal Gaussian output, and approximation error of the neural networks, respectively.
Remark 5.
According to Assumption 1 and Theorem 1, neural network parameters for an equivalent denoting fault observer error are bounded, which means that ψ m h 2 ( i ) ψ m h 2 M ( i ) , κ m h 2 ( i ) κ ¯ m h 2 ( i ) . The upper bound of the approximation error of the neural network adaptive term is denoted by κ ¯ m h 2 ( i ) .
To ensure the tracking performance under the fault condition and improve the control effect, design the second Lyapunov function as
V m h 2 ( i ) = 1 2 ln ( ( w b d ( i ) ) 2 [ ( w b d ( i ) ) 2 ( E m h 2 ( i ) ) 2 ] 1 ) + 1 2 ( E m h 2 ( i ) ) 2
where w b d ( i ) are positive constants as the upper bound value of the speed tracking error of a unilateral electromechanical actuator.
The differential of Equation (41) can be written as
V ˙ m h 2 ( i ) = [ ( w b d ( i ) ) 2 ( E m h 2 ( i ) ) 2 ] 1 E m h 2 ( i ) E ˙ m h 2 ( i ) + E m h 2 ( i ) E ˙ m h 2 ( i ) = ( γ b d ( i ) + 1 ) E m h 2 ( i ) E ˙ m h 2 ( i )
To make V ˙ m h 2 ( i ) < 0 , design the input value of the electromechanical actuator as
u m h ( i ) = 1 β m h 1 ( i ) ( β m h 2 ( i ) z m h 3 ( i ) + β m h 3 ( i ) z m h 2 ( i ) η m h 2 ( i ) sgn [ ( γ b d ( i ) + 1 ) E m h 2 ( i ) ] ( ψ ^ m h 2 ( i ) ) T R m h 2 ( i ) σ m h 2 ( i ) E m h 2 ( i ) + α f o 2 ( i ) e f o 2 ( i ) + δ f o 2 ( i ) χ f o 2 ( i ) + β ^ f o 2 ( i ) ( γ b d ( i ) + 1 ) 1 χ m h 1 ( i ) E m h 2 ( i ) γ b d ( i ) ( γ b d ( i ) + 1 ) 1 χ m h 2 ( i ) E m h 2 ( i ) ( γ b d ( i ) ) 2 ( γ b d ( i ) + 1 ) 1 χ m h 3 ( i ) E m h 2 ( i ) + z ¯ ˙ m h d 3 ( i ) )
where η m h 2 ( i ) , σ m h 2 ( i ) , χ m h 1 ( i ) , χ m h 2 ( i ) , and χ m h 3 ( i ) are positive coefficients. ( ψ ^ m h 2 ( i ) ) T R m h 2 ( i ) denote the disturbance and fault adaptive compensation item, the change rate of which can be defined as
ψ ^ ˙ m h 2 ( i ) = ϖ m h 2 ( i ) [ E m h 2 ( i ) R m h 2 ( i ) ( γ b d ( i ) + 1 ) ρ m h 2 ( i ) ψ ^ m h 2 ( i ) ]
where ϖ m h ( i ) and ρ m h ( i ) are positive constants.
Remark 6.
u m h ( i ) as the real control law will be converted into the control current of the motor driver of the electric actuator. The stability and performance of electromechanical actuator speed closed-loop tracking control can be ensured by choosing coefficients in (43) and (44). By introducing the adaptive term, the compensation of the fault and modeling error of actuators is realized. To handle the approximation errors which affect the control effect easily, the error-related barrier function and robust term are designed into the control law whose parameters are adjusted according to how well the designed neural network fits with uncertainties.
Theorem 2.
For the tension coordination control system with electromechanical actuator faults described by Equation (26), using control law (34) and (43), and the filter Equation (36), then all the signals in the closed-loop control system finally have a unanimous boundary in the end, and there are tracking errors to converge into a small boundary range.
Proof of Theorem 2.
First, define the Lyapunov function as shown below:
V m h w = E m h 1 2 2 + i = 1 2 1 2 ln ( ( w b d ( i ) ) 2 [ ( w b d ( i ) ) 2 ( E m h 2 ( i ) ) 2 ] 1 ) + i = 1 2 ( E m h 2 ( i ) ) 2 2 + i = 1 2 ( Γ m h 3 ( i ) ) 2 2 + ( ψ ˜ m h 1 ) T ψ ˜ m h 1 2 ϖ m h 1 + i = 1 2 ( ψ ˜ m h 2 ( i ) ) T ψ ˜ m h 2 ( i ) 2 ϖ m h 2 ( i )
where Γ m h 3 ( i ) = z ¯ m h d 3 ( i ) z m h d 3 ( i ) (i = 1, 2) indicate the filter error, ψ ˜ m h 1 = ψ ^ m h 1 ψ m h 1 , ψ ˜ m h 2 ( i ) = ψ ^ m h 2 ( i ) ψ m h 2 ( i ) .
Making the differential on both sides of Equation (45), we can obtain
V ˙ m h w = E m h 1 E ˙ m h 1 + i = 1 2 ( γ b d ( i ) + 1 ) E m h 2 ( i ) E ˙ m h 2 ( i ) + i = 1 2 Γ m h 3 ( i ) Γ ˙ m h 3 ( i ) + ( ϖ m h 1 ) 1 ( ψ ˜ m h 1 ) T ψ ˜ ˙ m h 1 + i = 1 2 ( ϖ m h 2 ( i ) ) 1 ( ψ ˜ m h 2 ( i ) ) T ψ ˜ ˙ m h 2 ( i )
Based on (36) and (37), the derivative of the filter error can be expressed as
Γ ˙ m h 3 ( i ) = ( ϑ m h 2 ( i ) ) 1 Γ m h 3 ( i ) + ζ m h 3 ( i )
where ζ m h 3 ( i ) = z ˙ m h d 3 ( i ) .
Combining (29), (34), and (37), the derivative of E m h 1 can be written as
E ˙ m h 1 = k e m a [ z b c z g b Γ m h 3 ( 1 ) z b c z g b Γ m h 3 ( 2 ) + z b c z g b E m h 2 ( 1 ) z b c z g b E m h 2 ( 2 ) ] / ( 1 + sin θ q x ) ( σ m h 1 E m h 1 + ( ψ ^ m h 1 ) T R m h 1 + η m h 1 sgn ( E m h 1 ) ) + f m h 1
From (38) and (43), the following inequality is yielded:
E ˙ m h 2 ( i ) = η m h 2 ( i ) sgn [ ( γ b d ( i ) + 1 ) E m h 2 ( i ) ] ( γ b d ( i ) + 1 ) 1 χ m h 1 ( i ) E m h 2 ( i ) γ b d ( i ) ( γ b d ( i ) + 1 ) 1 χ m h 2 ( i ) E m h 2 ( i ) ( γ b d ( i ) ) 2 ( γ b d ( i ) + 1 ) 1 χ m h 3 ( i ) E m h 2 ( i ) ( ψ ˜ m h 2 ( i ) ) T R m h 2 ( i ) κ m h 2 ( i )
Substituting (47), (48), (49) to (46), V ˙ m h w becomes
V ˙ m h w = E m h 1 ( [ z b c z g b Γ m h 3 ( 1 ) z b c z g b Γ m h 3 ( 2 ) + z b c z g b E m h 2 ( 1 ) z b c z g b E m h 2 ( 2 ) ] / 2 α m h 1 ( σ m h 1 E m h 1 + η m h 1 sgn ( E m h 1 ) ) ( ψ ˜ m h 1 ) T R m h 1 κ m h 1 ) + i = 1 2 ( γ b d ( i ) + 1 ) E m h 2 ( i ) ( η m h 2 ( i ) sgn [ ( γ b d ( i ) + 1 ) E m h 2 ( i ) ] ( γ b d ( i ) + 1 ) 1 χ m h 1 ( i ) E m h 2 ( i ) γ b d ( i ) ( γ b d ( i ) + 1 ) 1 χ m h 2 ( i ) E m h 2 ( i ) ( γ b d ( i ) ) 2 ( γ b d ( i ) + 1 ) 1 χ m h 3 ( i ) E m h 2 ( i ) ( ψ ˜ m h 2 ( i ) ) T R m h 2 ( i ) κ m h 2 ( i ) ) + i = 1 2 Γ m h 3 ( i ) ( ( ϑ m h 2 ( i ) ) 1 Γ m h 3 ( i ) + Ψ m h 3 ( i ) ) + ( ϖ m h 1 ) 1 ( ψ ˜ m h 1 ) T ψ ˜ ˙ m h 1 + i = 1 2 ( ϖ m h 2 ( i ) ) 1 ( ψ ˜ m h 2 ( i ) ) T ψ ˜ ˙ m h 2 ( i )
Select parameters that meet the following conditions:
η m h 1 > κ ¯ m h 1   , η m h 2 ( i ) > κ ¯ m h 2 ( i )
Remark 7.
It can be concluded from (48) and (49) that by selecting suitable η m h 1 and η m h 2 ( i ) , the adaptive error will be effectively restrained.
Substituting the adaptive laws (35) and (44), we can further obtain
V ˙ m h w E m h 1 z b c z g b [ Γ m h 3 ( 1 ) Γ m h 3 ( 2 ) + E m h 2 ( 1 ) E m h 2 ( 2 ) ] / 2 α m h 1 + i = 1 2 Γ m h 3 ( i ) ( ( ϑ m h 2 ( i ) ) 1 Γ m h 3 ( i ) + ζ m h 3 ( i ) ) σ m h 1 ( E m h 1 ) 2 i = 1 2 χ m h 1 ( i ) ( E m h 2 ( i ) ) 2 i = 1 2 γ b d ( i ) χ m h 2 ( i ) ( E m h 2 ( i ) ) 2 i = 1 2 χ m h 3 ( i ) ( γ b d ( i ) ) 2 ( E m h 2 ( i ) ) 2 ρ m h 1 ( ψ ˜ m h 1 ) T ψ ^ m h 1 i = 1 2 ρ m h 2 ( i ) ( ψ ˜ m h 2 ( i ) ) T ψ ^ m h 2 ( i )
It is clear that the following inequality is established [47]:
{ ( ψ ˜ m h 1 ) T ( ψ ^ m h 1 ) 1 2 ( ψ m h 1 ) T ( ψ m h 1 ) 1 2 ( ψ ˜ m h 1 ) T ψ ˜ m h 1 ( ψ ˜ m h 2 ( i ) ) T ( ψ ^ m h 2 ( i ) ) 1 2 ( ψ m h 2 ( i ) ) T ( ψ m h 2 ( i ) ) 1 2 ( ψ ˜ m h 2 ( i ) ) T ψ ˜ m h 2 ( i ) ln γ b d ( i ) ( w b d ( i ) ) 2 ( E m h 2 ( i ) ) 2 γ b d ( i )
Combining (52) and (53) yields
V ˙ m h w = [ z b c z g b α m h 1 σ m h 1 ] ( E m h 1 ) 2 + i = 1 2 ( z b c z g b 4 α m h 1 χ m h 1 ( i ) ) ( E m h 2 ( i ) ) 2 + i = 1 2 ( χ m h 2 ( i ) ln γ b d ( i ) ( w b d ( i ) ) 2 ) 0.5 ρ m h 1 ( ψ ˜ m h 1 ) T ( ψ ˜ m h 1 ) + i = 1 2 ( z b c z g b 4 α m h 1 + ( ζ m h 3 ( i ) ) 2 2 1 ϑ m h 2 ( i ) ) ( Γ m h 3 ( i ) ) 2 + ι m h w i = 1 2 ( 0.5 ρ m h 2 ( i ) ( ψ ˜ m h 2 ( i ) ) T ψ ˜ m h 2 ( i ) )
where ι m h w = 1 + 0.5 ρ m h 1 ( ψ m h 1 M ) 2 + i = 1 2 ( 0.5 ρ m h 2 ( i ) ( ψ m h 2 M ( i ) ) 2 ) . As indicated in Assumption 3 and Remark 4, ι m h w is bounded.
Let the coefficients satisfy the following conditions:
σ m h 1 z b c z g b ( α m h 1 ) 1 + ε m h , χ m h 1 ( i ) z b c z g b ( 4 α m h 1 ) 1 + ε m h , χ m h 2 ( i ) ε m h ρ m h 1 2 ε m h ϖ m h 1 , ρ m h 2 ( i ) 2 ε m h ϖ m h 2 ( i ) , 1 ϑ m h 2 ( i ) z b c z g b 4 α m h 1 + ( H m h 3 ( i ) ) 2 2 + ε m h
where ε m h is a positive scalar. H m h 3 ( i ) = ζ m h 3 ( i ) max represent the upper bound of the ζ m h 3 ( i ) .
Then, it can be derived that
V ˙ m h w ε m h ( E m h 1 ) 2 i = 1 2 ε m h ( E m h 2 ( i ) ) 2 ε m h χ m h 2 ( i ) ln γ b d ( i ) ( w b d ( i ) ) 2 i = 1 2 ε m h ( Γ m h 3 ( i ) ) 2 + ι m h w ( ϖ m h 1 ) 1 ε m h ( ψ ˜ m h 1 ) T ψ ˜ m h 1 i = 1 2 ( ϖ m h 2 ( i ) ) 1 ε m h ( ψ ˜ m h 2 ( i ) ) T ψ ˜ m h 2 ( i ) 2 ε m h V m h w + ι m h w
V m h w is thus bounded; in terms of the Lyapunov stability theory, all signals of the steel wire rope tension coordination control system represent a semi-global unity. Theorem 2 is proved. □

4. Experiment Result and Analysis

To verify the proposed electromechanical actuator fault-tolerant control method, a xpc technology-based experimental system involving the mine hoist mechanical structure, different kinds of sensors, and two computers was built as shown in Figure 6. The tension sensor connects the end of the vertical rope and skip which is used to measure the tension of the wire rope. Three encoders are fixed on the drum shaft end and two headgear sheaves’ shaft end, respectively. Each electric actuator consists of a 400 W motor, a gearbox whose reduction ratio is 1:10, and a screw structure. Sensor signals are acquired by PCI-6601 and PCI-6024E data collection cards which are installed in an Advantech target computer mainboard slot with sampling frequency 1 khz. According to the dimensions and specifications of the experimental system hardware, relevant parameters of the mathematical model can be set as follows: θqx = 45°, ρ w r = 0.004 kg/m, r s ( i ) = 0.03 m, ψ m t f ( i ) = 0.32 Wb, L m t d ( i ) = L m t q ( i ) = 2.8 mH , J s ( i ) = 0.005 kg·m2, p m t ( i ) = 4 , E A = 2 × 10 6 N.
Considering that the failure of the motor and the electronic unit during the actual operate process takes place inside the equipment, the fault type, time of occurrence, and fault amplitude are not easy to control. Meanwhile, as faults mainly affect the opening response of the electromechanical actuator, add additional disturbance signals to the motor drive input to simulate faults of the right headgear sheave actuator subsystem as follows:
{ f a f 2 ( 2 ) = z m u l f i m t q ( 2 ) f a f 1 ( 2 ) = z a d d f sin 5 π ( t 2.5 ) f a f 1 ( 2 ) = z a d d f sin 5 π ( t 8.2 ) 0 2.5   s < t < 4   s 2.5   s < t < 4   s 8.2   s < t < 10.1   s else and 8.2   s < t < 10.1   s
where z m u l f and z a d d f are fault coefficients and can be modified in the computer program. The reason for picking the birth time of faults is that the difference between the length of the two ropes occurs in the winding process, which can reflect the actuator performance.
To evaluate the performance of the developed actuator fault-tolerant control method for wire rope tension coordination control, two indicators of the root mean square error (RMSE) and maximum absolute error (MAE) are employed for the evaluation. The calculation method of the two indicators can be expressed as
RMSE = N t t = 1 N t s | F c n , N t t ( 1 ) F c n , N t t ( 2 ) | 2
MAE = max N t t = 1 N t s | F c n , N t t ( 1 ) F c n , N t t ( 2 ) |
where F c n , N t t ( i ) denotes the rope tension value at the Ntt-th sampling point. N t s represents the number of sampling points during the experimental process.
For illustrating the effectiveness of the above fault-tolerant controller design more clearly, the following two controllers are compared.
(1)
AFSMC: This is the adaptive fuzzy sliding mode controller with a fuzzy smooth switching item as discussed in [48]. The tension coordination control loop applies an improved sliding mode algorithm to overcome the interference, and the actuator control loop employs the traditional proportion-integration-differentiation (PID) algorithm. The input voltage amplitude of the left side actuator is equal to the right side actuator, and the symbols of the two side inputs are opposite. The gains of the controller are chosen as: K a f = 2 × 10 5 , K v f = 2.25 , K d f = 0.028 , σ = 2.25 , γ = 0.042 , η = 10 3 .
(2)
AFODST: This is the adaptive fault observer (12) and dynamic surface technology-based fault tolerant controller (43) proposed in this paper according to the theoretical design procedure. The parameters of the actuator fault observer are obtained by using the LMI toolbox to solve according to formula (15). The final gains are: α f o ( 1 ) = α f o ( 2 ) = diag ( 10 2 , 0.016 ) , δ f o ( 1 ) = δ f o ( 1 ) = diag ( 0.003 , 0.015 ) , γ f o ( 1 ) = γ f o ( 2 ) = diag ( 0.06 , 1.28 ) , θ f o 1 ( 1 ) = θ f o 1 ( 2 ) = 10 6 , θ f o 2 ( 1 ) = θ f o 2 ( 2 ) = 2 × 10 4 , p s o f 1 ( 1 ) = p s o f 2 ( 1 ) = p s o f 1 ( 2 ) = p s o f 2 ( 2 ) = 3, p x o f 1 ( 1 ) = p x o f 2 ( 1 ) = p x o f 1 ( 2 ) = p x o f 2 ( 2 ) = 5. In addition, the gains of the improved dynamic surface fault-tolerant controller are chosen as: ϑ m h 2 ( 1 ) = ϑ m h 2 ( 1 ) = 10 3 , σ m h 2 ( 1 ) = σ m h 2 ( 2 ) = 0.02 , w b d ( 1 ) = w b d ( 2 ) = 10 3 , σ m h 1 = 0.035 , ϖ m h 1 = 2.4 , η m h 1 = 2.1 × 10 5 , ϖ m h 2 ( 1 ) = ϖ m h 2 ( 2 ) = 0.035 , χ m h 1 ( 1 ) = χ m h 1 ( 2 ) = 0.028 , χ m h 2 ( 1 ) = χ m h 2 ( 2 ) = 0.019 , χ m h 3 ( 1 ) = χ m h 3 ( 2 ) = 0.026 , η m h 2 ( 1 ) = η m h 2 ( 2 ) = 3 × 10 4 , ρ m h 1 = 0.035 , ρ m h 2 ( 1 ) = ρ m h 2 ( 2 ) = 0.022 .
The rope tension coordination control effect performances under the two controllers in a normal case are shown in Figure 7. As can be seen, the proposed AFODST and AFSMC controllers have approximately the same effect and are better than the traditional PID controller which was used in the literature [49] dealing with random drum winding disturbance since the above two controllers both made use of an adaptive function to adjust the control input and compensate dynamic modeling uncertainties. As shown in Figure 7d, the maximum ratio of the tension difference to the average tension under the effect of the AFSMC is 5.22%, and the maximum ratio under the action of the AFODST is 4.74%; the fluctuation ranges are both within 10%.
The experimental results of two controllers under actuator fault conditions are provided in Figure 8, Figure 9 and Figure 10. In addition, a qualitative comparison of tension control performance indices between the two controllers under four conditions is given in Table 2. It can be seen that along with the value of the fault parameter zmulf, increasing from 0.3 to 1, fluctuation ranges for the RMSE and MAE under the action of the AFSMC are [20.062 N~66.527 N] and [0.927 N~2.802 N]. In the case of condition 2, the RMSE and MAE of the tension difference under AFSMC action grew by 13.94 percent and 7.55 percent, respectively, compared with condition 1. The RMSE and MAE rose by 36.45 percent and 27.29 percent with the effect of the AFSMC in the case of condition 3 compared with the no faults condition. With the action of the AFSMC, in the case of condition 4, the RMSE and MAE increased 231.61 percent and 202.27 percent compared with the case of the normal condition. Those indicate that the AFSMC can overcome the set fault f a f 1 ( 2 ) very well; the tension coordination control effect is basically the same as in the no faults condition when amplitudes of the set addition and multiplication electromechanical actuator faults are small (zaddf = 0.1, zmulf = 0.3). However, the effect is gradually getting worse with the coefficient increase of f a f 2 ( 2 ) , as the effectiveness of the single actuator is getting weaker, and the tension compensation effect cannot meet the requirements.
Figure 9 illustrates the outputs of fault observers in the above fault case. In Figure 9a,b, it can be seen that the actuators implement actions in real time to control rope tension, and the states of actuators are well tracked by the estimated state in all cases. Meanwhile, when the actuators are in a normal state, the running state of the actuators is basically maintained in many experiments, and the actuators on both sides carry out the control task together. During the period of actuators’ faults, the right actuator gradually loses control as the fault severity increases; by contrast, the left actuator takes on more of the control task under the action of the fuzzy assignment factor so that the controller can complete the control objectives. Figure 9c,d present the estimation of fault influence factors. It can be noticed that although there is no fault set for the left actuator, there is still a certain fault estimation value, which is caused by modeling error and uncertainties. In addition, when the fault parameter of the right actuator increases, the amplitude of the fault index is also becoming larger, which can be well tracked and estimated. By analyzing the fault observer results in Figure 9, it is obvious that the fault observers have a good fault detection effect for one thing and can effectively guide the control allocation of the controller for another.
From Figure 10, it can be shown that fluctuation ranges for the RMSE and MAE under the action of the proposed strategy are [18.171 N~21.399 N] and [0.803 N~1.015 N]. In the case of condition 2, the RMSE and MAE of the tension difference under AFODST action are reduced by 20.21 percent and 19.46 percent, respectively, compared with under AFSMC action. The RMSE and MAE lessen by 34.47 percent and 28.22 percent with the effect of the AFODST in the case of condition 3 compared with under the action of the AFSMC. With the action of the AFODST, in the case of condition 4, the RMSE and MAE lower 67.83 percent and 63.78 percent compared with under the effect of the AFSMC.
For the experimental results under the four conditions, the superiority of the proposed electromechanical actuator fault-tolerant controller is definitely proven. The proposed strategy can not only effectively overcome the disturbance effect as the conventional controller, it can also reduce the impact of the time-varying actuator fault. From the point of view of the rope tension control effect and controller robustness, the proposed method is more applicable than the comparative AFSMC method to the tension coordination control under variable operating conditions.

5. Conclusions

In this paper, an adaptive dynamic surface technology-based electromechanical actuator fault tolerant scheme is designed for the Blair mine hoist. Comparative experimental results of the proposed strategy with the AFSMC under drum winding disturbance and actuator fault conditions are obtained to illustrate the superiority of the designed hybrid controller under complex conditions in tension coordination controlling. The results of this work can be summarized and concluded as follows:
(1)
According to the actuator fault analysis and dynamic characteristic analysis of the mine hoist tension coordination control system, the overall mathematic model of the system includes the hoisting subsystem, skip subsystem, electric drive actuators, and the electromechanical actuator fault is built and expressed with a state space formulation.
(2)
In order to solve the tension coordination control problem under fault conditions, the RBF neural-network-based adaptive observer is applied to estimate the electromechanical actuator fault state. The hybrid method consists of the adaptive dynamic surface controller, fuzzy allocation factor, and the error barrier function is developed to achieve a better real-time control effect dealing with the sudden fault situation and unexpected disturbance. Each step of this algorithm is explicitly presented, and the stability of the tension coordination closed-loop control system is proved.
(3)
Through the simulated actuator faults’ experimental test, the proposed hybrid strategy shows stronger robustness and reacts faster when an electromechanical actuator performance loss fault occurs. Even when one of the actuators fails completely, the designed controller can achieve almost the same performance as the normal operating conditions. Compared with the AFSMC, the RMSE is reduced by 67.83 percent, and the MAE is reduced by 63.78 percent, in the worst setting situation. The validity of the theoretical analysis is thus verified.
However, since the designed controller contains various parameters that need to be set and adjusted, specific tune rules and optimization algorithms are urgently needed to be developed, which will be the focus in further research. Moreover, the design fault simulation strategy closer to the real situation is also the subsequent exploration direction.

Author Contributions

X.C. (Xiao Chen) designed the experiments and wrote the paper; Z.Z. contributed analysis tools and analyzed the data; T.M. performed the experiments and contributed materials; J.C. organized and analyzed experimental results; X.C. (Xiangdong Chang) reviewed and edited the paper; W.Z. performed formal analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Key Projects of Natural Science Research in Anhui Universities (KJ2021A0426), the University-level key projects of Anhui University of science and technology (xjzd2020-14), the Program for the University Synergy Innovation Program of Anhui Province (GXXT-2019-029), the start-up fund for introducing talents scientific research of Anhui University of Science and Technology. This work was also supported by the Institute of Energy, Hefei Comprehensive National Science Center under Grant No.21KZS215.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The kinetic energy of the hoisting subsystem can be denoted as
T h s = i = 1 2 1 2 m c n ( i ) d ˙ c n ( i ) + i = 1 2 1 2 J s ( i ) ( θ ˙ s ( i ) ) 2 + i = 1 2 1 2 m s ( i ) ( d ˙ s ( i ) ) 2 + 1 2 J w d θ ˙ w d 2 + i = 1 2 1 2 ρ w r + i = 1 2 1 2 ρ w r 0 l v ( i ) ( x ˙ y v ( i ) ) 2 dy v ( i ) + i = 1 2 1 2 ρ w r 0 l c ( i ) ( x ˙ y c ( i ) ) 2 dy c ( i ) = i = 1 2 1 6 ρ w r l c ( i ) [ ( θ ˙ w d r w d ( i ) ) 2 + ( θ ˙ s ( i ) r s ( i ) d ˙ s ( i ) sin θ q x ( i ) ) θ ˙ w d r w d ( i ) + ( θ ˙ s ( i ) r s ( i ) d ˙ s ( i ) sin θ q x ( i ) ) 2 ] + i = 1 2 1 6 ρ w r l v ( i ) [ ( θ ˙ s ( i ) r s ( i ) + d ˙ s ( i ) ) 2 + ( θ ˙ s ( i ) r s ( i ) + d ˙ s ( i ) ) d ˙ c n ( i ) + ( d ˙ c n ( i ) ) 2 ] + i = 1 2 1 2 m c n ( i ) d ˙ c n ( i ) + i = 1 2 1 2 J s ( i ) ( θ ˙ s ( i ) ) 2 + i = 1 2 1 2 m s ( i ) ( d ˙ s ( i ) ) 2 + 1 2 J w d θ ˙ w d 2
The potential energy of the hoisting subsystem can be represented as
U h s = i = 1 2 m c n ( i ) g d c n ( i ) + i = 1 2 m s ( i ) g d s ( i ) 1 2 i = 1 2 ρ w r l c ( i ) g [ ( θ w d r w d ( i ) + θ s ( i ) r s ( i ) ) sin θ q x ( i ) + d s ( i ) ] + 1 2 i = 1 2 ρ w r l v ( i ) g ( θ s ( i ) r s ( i ) + d s ( i ) + d c n ( i ) ) + 1 2 i = 1 2 k v ( i ) ( θ s ( i ) r s ( i ) + d s ( i ) d c n ( i ) ) 2 + 1 2 i = 1 2 k c ( i ) ( θ w d r w d ( i ) θ s ( i ) r s ( i ) + d s ( i ) sin θ q x ( i ) ) 2
The dissipation energy of the hoisting subsystem can be
D h s = 1 2 i = 1 2 c v ( i ) ( θ ˙ s ( i ) r s ( i ) d ˙ c n ( i ) + d ˙ s ( i ) ) 2 + 1 2 i = 1 2 c c ( i ) ( θ ˙ w d r w d ( i ) θ ˙ s ( i ) r s ( i ) + d ˙ s ( i ) sin θ q x ( i ) ) 2
The generalized force of the hoisting subsystem can be expressed as
F h s = [ F h s 1 F h s 2 F h s 3 F h s 4 F h s 5 F h s 6 F h s 7 ] T
where F h s 1 = τ d r u m + i = 1 2 0.5 ρ w r g l c ( i ) r w d ( i ) sin θ q x ( i ) , F h s 2 = 0.5 ρ w r g l c ( 1 ) r w d ( 1 ) sin θ q x ( 1 ) 0.5 ρ w r g l v ( 1 ) r s ( 1 ) , F h s 3 = 0.5 ρ w r g l c ( 2 ) r w d ( 2 ) sin θ q x ( 2 ) 0.5 ρ w r g l v ( 2 ) r s ( 2 ) , F h s 4 = F e m a ( 1 ) 0.5 ρ w r g l v ( 1 ) m s ( 1 ) g 0.5 ρ w r g l c ( 1 ) r w d ( 1 ) sin θ q x ( 1 ) , F h s 5 = F e m a ( 2 ) 0.5 ρ w r g l v ( 2 ) m s ( 2 ) g 0.5 ρ w r g l c ( 2 ) r w d ( 2 ) sin θ q x ( 2 ) , F h s 6 = F c n ( 1 ) m c n ( 1 ) g 0.5 ρ w r g l v ( 1 ) , F h s 7 = F c n ( 2 ) m c n ( 2 ) g 0.5 ρ w r g l v ( 2 ) . τ d r u m indicates the input torque of the drum.
The nonlinear dynamic equations of the hoisting subsystem (4)–(6) are then calculated by the following the Euler-Lagrange Equation [50].
d d t ( T h s x ˙ h s ) T h s x h s + U h s x h s + D h s x ˙ h s = F h s

References

  1. Giraud, L.; Galy, B. Fault tree analysis and risk mitigation strategies for mine hoists. Saf. Sci. 2018, 110, 222–234. [Google Scholar] [CrossRef]
  2. Pi, Y.; Zhang, J.; Tang, X.; Zhu, J. Three-dimensional dynamic modeling and simulation of a multi-cable winding hoister system considering bidirectional coupling between cage and flexible guides. J. Vib. Control 2022. [Google Scholar] [CrossRef]
  3. Yao, J.; Ma, Y.; Ma, C.; Xiao, X.; Xu, T. Effect of misalignment failures of steel guides on impact responses in friction mine hoisting systems. Eng. Fail. Anal. 2020, 118, 104841. [Google Scholar] [CrossRef]
  4. Peng, X.; Gong, X.S.; Liu, J.J. The study on crossover layouts of multi-layer winding grooves in deep mine hoists based on transverse vibration characteristics of catenary rope. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2019, 233, 118–132. [Google Scholar] [CrossRef]
  5. Chang, X.; Peng, Y.; Zhu, Z.; Gong, X.; Yu, Z.; Mi, Z.; Xu, C. Breaking failure analysis and finite element simulation of wear-out winding hoist wire rope. Eng. Fail. Anal. 2019, 95, 1–17. [Google Scholar] [CrossRef]
  6. Wang, J.; Pi, Y.; Krstic, M. Balancing and suppression of oscillations of tension and cage in dual-cable mining elevators. Automatica 2018, 98, 223–238. [Google Scholar] [CrossRef]
  7. Zhu, Z.C.; Li, X.; Shen, G.; Zhu, W.D. Wire rope tension control of hoisting systems using a robust nonlinear adaptive backstepping control scheme. ISA Trans. 2018, 72, 256–272. [Google Scholar] [CrossRef]
  8. Chen, X.; Zhu, Z.C.; Li, W.; Shen, G. Fault-tolerant control of double-rope winding hoists combining neural-based adaptive technique and iterative learning method. IEEE Access 2019, 7, 50476–50491. [Google Scholar] [CrossRef]
  9. Li, X.; Zhu, Z.; Shen, G.; Tang, Y. Wire Tension Coordination Control of Electro-Hydraulic Servo Driven Double-Rope Winding Hoisting Systems Using a Hybrid Controller Combining the Flatness-Based Control and a Disturbance Observer. Symmetry 2021, 13, 716. [Google Scholar] [CrossRef]
  10. Ding, X.; Shen, G.; Tang, Y.; Li, X. Axial vibration suppression of wire-ropes and container in double-rope mining hoists with adaptive robust boundary control. Mechatronics 2022, 85, 102817. [Google Scholar] [CrossRef]
  11. Zang, W.; Shen, G.; Rui, G.; Li, X.; Li, G.; Tang, Y. A flatness-based nonlinear control scheme for wire tension control of hoisting systems. IEEE Access 2019, 7, 146428–146442. [Google Scholar] [CrossRef]
  12. Wang, J.; Pi, Y.; Hu, Y.; Zhu, Z. State-observer design of a PDE-modeled mining cable elevator with time-varying sensor delays. IEEE Trans. Control Syst. Technol. 2019, 28, 1149–1157. [Google Scholar] [CrossRef]
  13. Ding, S.X.; Li, L. Control performance monitoring and degradation recovery in automatic control systems: A review, some new results, and future perspectives. Control Eng. Pract. 2021, 111, 104790. [Google Scholar] [CrossRef]
  14. Shen, Q.; Yue, C.; Goh, C.H.; Wang, D. Active fault-tolerant control system design for spacecraft attitude maneuvers with actuator saturation and faults. IEEE Trans. Ind. Electron. 2018, 66, 3763–3772. [Google Scholar] [CrossRef]
  15. Wang, Y.; Jiang, B.; Wu, Z.G.; Xie, S.; Peng, Y. Adaptive sliding mode fault-tolerant fuzzy tracking control with application to unmanned marine vehicles. IEEE Trans. Syst. Man Cybern. Syst. 2020, 51, 6691–6700. [Google Scholar] [CrossRef]
  16. Tian, G.; Yuan, G.; Aleksandrov, A.; Zhang, T.; Li, Z.; Fathollahi-Fard, A.M.; Ivanov, M. Recycling of spent Lithium-ion Batteries: A comprehensive review for identification of main challenges and future research trends. Sustain. Energy Technol. Assess. 2022, 53, 102447. [Google Scholar] [CrossRef]
  17. Chen, Z.; Yang, X.; Liu, X. RBFNN-based nonsingular fast terminal sliding mode control for robotic manipulators including actuator dynamics. Neurocomputing 2019, 362, 72–82. [Google Scholar] [CrossRef]
  18. Amin, A.A.; Hasan, K.M. A review of fault tolerant control systems: Advancements and applications. Measurement 2019, 143, 58–68. [Google Scholar] [CrossRef]
  19. Gao, Z.; Cecati, C.; Ding, S.X. A survey of fault diagnosis and fault-tolerant techniques—Part I: Fault diagnosis with model-based and signal-based approaches. IEEE Trans. Ind. Electron. 2015, 62, 3757–3767. [Google Scholar] [CrossRef] [Green Version]
  20. Chen, T.; Chen, L.; Xu, X.; Cai, Y.; Jiang, H.; Sun, X. Passive fault-tolerant path following control of autonomous distributed drive electric vehicle considering steering system fault. Mech. Syst. Signal Process. 2019, 123, 298–315. [Google Scholar] [CrossRef]
  21. Başak, H.; Kemer, E.; Prempain, E. A passive fault-tolerant switched control approach for linear multivariable systems: Application to a quadcopter unmanned aerial vehicle model. J. Dyn. Syst. Meas. Control 2020, 142, 031004. [Google Scholar] [CrossRef]
  22. Hu, C.; Wei, X.; Ren, Y. Passive fault-tolerant control based on weighted LPV tube-MPC for air-breathing hypersonic vehicles. Int. J. Control Autom. Syst. 2019, 17, 1957–1970. [Google Scholar] [CrossRef]
  23. Hao, L.Y.; Zhang, H.; Li, T.S.; Lin, B.; Chen, C.L.P. Fault tolerant control for dynamic positioning of unmanned marine vehicles based on TS fuzzy model with unknown membership functions. IEEE Trans. Veh. Technol. 2021, 70, 146–157. [Google Scholar] [CrossRef]
  24. Aydin, M.N.; Coban, R. PID sliding surface-based adaptive dynamic second-order fault-tolerant sliding mode control design and experimental application to an electromechanical system. Int. J. Control 2022, 95, 1767–1776. [Google Scholar] [CrossRef]
  25. Abbaspour, A.; Mokhtari, S.; Sargolzaei, A.; Yen, K.K. A survey on active fault-tolerant control systems. Electronics 2020, 9, 1513. [Google Scholar] [CrossRef]
  26. Xiao, B.; Yin, S. A deep learning based data-driven thruster fault diagnosis approach for satellite attitude control system. IEEE Trans. Ind. Electron. 2020, 68, 10162–10170. [Google Scholar] [CrossRef]
  27. Rudin, K.; Ducard GJ, J.; Siegwart, R.Y. Active fault-tolerant control with imperfect fault detection information: Applications to UAVs. IEEE Trans. Aerosp. Electron. Syst. 2019, 56, 2792–2805. [Google Scholar] [CrossRef]
  28. Abbaspour, A.; Yen, K.K.; Forouzannezhad, P.; Sargolzaei, A. A neural adaptive approach for active fault-tolerant control design in UAV. IEEE Trans. Syst. Man Cybern. Syst. 2018, 50, 3401–3411. [Google Scholar] [CrossRef]
  29. Xiong, H.; Liao, Y.; Chu, X.; Nian, X.; Wang, H. Observer based fault tolerant control for a class of Two-PMSMs systems. ISA Trans. 2018, 80, 99–110. [Google Scholar] [CrossRef] [PubMed]
  30. Nguyen, V.-C.; Tran, X.-T.; Kang, H.-J. A Novel High-Speed Third-Order Sliding Mode Observer for Fault-Tolerant Control Problem of Robot Manipulators. Actuators 2022, 11, 259. [Google Scholar] [CrossRef]
  31. Han, J.; Liu, X.; Wei, X.; Sun, S. A dynamic proportional-integral observer-based nonlinear fault-tolerant controller design for nonlinear system with partially unknown dynamic. IEEE Trans. Syst. Man Cybern. Syst. 2021, 52, 5092–5104. [Google Scholar] [CrossRef]
  32. Lei, Q.; Ma, Y.; Liu, J.; Yu, J. Neuroadaptive observer-based discrete-time command filtered fault-tolerant control for induction motors with load disturbances. Neurocomputing 2021, 423, 435–443. [Google Scholar] [CrossRef]
  33. Truong, H.V.A.; Ahn, K.K. Actuator failure compensation-based command filtered control of electro-hydraulic system with position constraint. ISA Trans. 2022. [Google Scholar] [CrossRef]
  34. Li, L.; Yao, L.; Wang, H.; Gao, Z. Iterative learning fault diagnosis and fault tolerant control for stochastic repetitive systems with Brownian motion. ISA Trans. 2022, 121, 171–179. [Google Scholar] [CrossRef]
  35. Chen, X.; Zhu, Z.-C.; Ma, T.-B.; Shen, G. Model-based sensor fault detection, isolation and tolerant control for a mine hoist. Meas. Control 2022, 55, 274–287. [Google Scholar] [CrossRef]
  36. Guzmán-Rabasa, J.A.; López-Estrada, F.R.; González-Contreras, B.M.; Valencia-Palomo, G.; Chadli, M.; Pérez-Patricio, M. Actuator fault detection and isolation on a quadrotor unmanned aerial vehicle modeled as a linear parameter-varying system. Meas. Control 2019, 52, 1228–1239. [Google Scholar] [CrossRef] [Green Version]
  37. Mao, Z.; Yan, X.G.; Jiang, B.; Chen, M. Adaptive fault-tolerant sliding-mode control for high-speed trains with actuator faults and uncertainties. IEEE Trans. Intell. Transp. Syst. 2019, 21, 2449–2460. [Google Scholar] [CrossRef] [Green Version]
  38. Shiravani, F.; Alkorta, P.; Cortajarena, J.A.; Barambones, O. An Enhanced Sliding Mode Speed Control for Induction Motor Drives. Actuators 2022, 11, 18. [Google Scholar] [CrossRef]
  39. Wang, Y.; Xu, N.; Liu, Y.; Zhao, X. Adaptive fault-tolerant control for switched nonlinear systems based on command filter technique. Appl. Math. Comput. 2021, 392, 125725. [Google Scholar] [CrossRef]
  40. Zhao, Z.; Ren, Y.; Mu, C.; Zou, T.; Hong, K.-S. Adaptive neural-network-based fault-tolerant control for a flexible string with composite disturbance observer and input constraints. IEEE Trans. Cybern. 2021. [Google Scholar] [CrossRef]
  41. Huang, J.; Cao, Y.; Xiong, C.; Zhang, H.-T. An echo state Gaussian process-based nonlinear model predictive control for pneumatic muscle actuators. IEEE Trans. Autom. Sci. Eng. 2018, 16, 1071–1084. [Google Scholar] [CrossRef]
  42. Guo, Y.; Zhang, D.; Zhang, X.; Wang, S.; Ma, W. Experimental study on the nonlinear dynamic characteristics of wire rope under periodic excitation in a friction hoist. Shock Vib. 2020, 2020, 8506016. [Google Scholar] [CrossRef] [Green Version]
  43. Chaoui, H.; Khayamy, M.; Okoye, O. Adaptive RBF network based direct voltage control for interior PMSM based vehicles. IEEE Trans. Veh. Technol. 2018, 67, 5740–5749. [Google Scholar] [CrossRef]
  44. Qiao, G.; Liu, G.; Shi, Z.; Wang, Y.; Ma, S.; Lim, T.C. A review of electromechanical actuators for More/All Electric aircraft systems. Proc. Inst. Mech. Eng. Part C: J. Mech. Eng. Sci. 2018, 232, 4128–4151. [Google Scholar] [CrossRef] [Green Version]
  45. Smaeilzadeh, S.M.; Golestani, M. A finite-time adaptive robust control for a spacecraft attitude control considering actuator fault and saturation with reduced steady-state error. Trans. Inst. Meas. Control 2019, 41, 1002–1009. [Google Scholar] [CrossRef]
  46. Costa, T.V.; Sencio, R.R.; Oliveira-Lopes, L.C.; Silva, F.V. Fault-tolerant control by means of moving horizon virtual actuators: Concepts and experimental investigation. Control Eng. Pract. 2021, 107, 104683. [Google Scholar] [CrossRef]
  47. Edalati, L.; Sedigh, A.K.; Shooredeli, M.A.; Moarefianpour, A. Adaptive fuzzy dynamic surface control of nonlinear systems with input saturation and time-varying output constraints. Mech. Syst. Signal Process. 2018, 100, 311–329. [Google Scholar] [CrossRef]
  48. Li, X.; Zhu, Z.C.; Shen, G. A switching-type controller for wire rope tension coordination of electro-hydraulic-controlled double-rope winding hoisting systems. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2016, 230, 1126–1144. [Google Scholar] [CrossRef]
  49. Chen, X.; Zhu, Z.; Shen, G.; Li, W. Tension coordination control of double-rope winding hoisting system using hybrid learning control scheme. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2019, 233, 1265–1281. [Google Scholar] [CrossRef]
  50. Huang, J.; Luo, C.; Yu, P.; Hao, H. A methodology for calculating limit deceleration of flexible hoisting system: A case study of mine hoist. Proc. Inst. Mech. Eng. Part E Ournal Process Mech. Eng. 2020, 234, 342–352. [Google Scholar] [CrossRef]
Figure 1. The structure of Blair mine hoist.
Figure 1. The structure of Blair mine hoist.
Actuators 11 00299 g001
Figure 2. Mine hoist dynamic model. (a) Hoisting subsystem model; (b) Skip subsystem model; (c) Electric drive actuator model.
Figure 2. Mine hoist dynamic model. (a) Hoisting subsystem model; (b) Skip subsystem model; (c) Electric drive actuator model.
Actuators 11 00299 g002
Figure 3. Fault Tree of Electromechanical Actuators in Tension Coordination Control System.
Figure 3. Fault Tree of Electromechanical Actuators in Tension Coordination Control System.
Actuators 11 00299 g003
Figure 4. Overall Architecture of Electromechanical Actuators Fault-Tolerant Control System.
Figure 4. Overall Architecture of Electromechanical Actuators Fault-Tolerant Control System.
Actuators 11 00299 g004
Figure 5. Radial Basis Neural Network Structure.
Figure 5. Radial Basis Neural Network Structure.
Actuators 11 00299 g005
Figure 6. xpc-based mine hoist platform.
Figure 6. xpc-based mine hoist platform.
Actuators 11 00299 g006
Figure 7. Measured rope tension and tension difference under normal condition: (a) Rope tension with AFSMC; (b) Rope tension with AFODST; (c) Tension difference; (d) Ratio of the tension difference to the average tension.
Figure 7. Measured rope tension and tension difference under normal condition: (a) Rope tension with AFSMC; (b) Rope tension with AFODST; (c) Tension difference; (d) Ratio of the tension difference to the average tension.
Actuators 11 00299 g007
Figure 8. Measured rope tension and tension difference under fault condition with AFSMC: (a) Zmulf = 0.3; (b) Zmulf = 0.6; (c) Zmulf = 1; (d) Tension difference.
Figure 8. Measured rope tension and tension difference under fault condition with AFSMC: (a) Zmulf = 0.3; (b) Zmulf = 0.6; (c) Zmulf = 1; (d) Tension difference.
Actuators 11 00299 g008
Figure 9. Fault observer output under fault conditions: (a) Estimated state of left actuator; (b) Estimated state of right actuator; (c) Fault estimation output of left actuator; (d) Fault estimation output of right actuator.
Figure 9. Fault observer output under fault conditions: (a) Estimated state of left actuator; (b) Estimated state of right actuator; (c) Fault estimation output of left actuator; (d) Fault estimation output of right actuator.
Actuators 11 00299 g009
Figure 10. Measured rope tension and tension difference under fault condition with AFODST: (a) Zmulf = 0.3; (b) Zmulf = 0.6; (c) Zmulf = 1; (d) Tension difference.
Figure 10. Measured rope tension and tension difference under fault condition with AFODST: (a) Zmulf = 0.3; (b) Zmulf = 0.6; (c) Zmulf = 1; (d) Tension difference.
Actuators 11 00299 g010
Table 1. Electromechanical Actuator Failure Event and the Corresponding ID.
Table 1. Electromechanical Actuator Failure Event and the Corresponding ID.
Event IDEventEvent IDEvent
S11Motor Winding FailureS222Board Circuit Failure
S111Winding Short CircuitS31Gear Failure
S112Winding Open CircuitS311Gear Surface Wear
S121Rotor EccentricityS312Gear Teeth Broken
S131Unstable Supply VoltageS313Gear Teeth Deformation
S14Motor Sensor FailureS32Bearing Failure
S141Encoder FailureS321Structure Plastic Deformation
S142Current Sensor FailureS322Surface Cracking and Spalling
S21Cable FailureS323Surface Corrosion
S211Cable Open CircuitS33Screw Mechanism Failure
S212Cable Short CircuitS331Guide Wear
S22Control Board FailureS332Screw Wear
S221Board LooseS333Screw Deformation
Table 2. Tension Control Performance indexes.
Table 2. Tension Control Performance indexes.
ControllerRMSEMAERMSEMAE
Condition 1—No faultsCondition 2—zaddf = 0.1, zmulf = 0.3
AFSMC20.0620.92722.8580.997
AFODST18.1710.83318.2390.803
Condition 3—zaddf = 0.1, zmulf = 0.6Condition 4—zaddf = 0.1, zmulf = 1
AFSMC27.3741.18066.5272.802
AFODST17.9380.84721.3991.015
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, X.; Zhu, Z.; Ma, T.; Chang, J.; Chang, X.; Zang, W. An Adaptive Dynamic Surface Technology-Based Electromechanical Actuator Fault-Tolerant Scheme for Blair Mine Hoist Wire Rope Tension Control System. Actuators 2022, 11, 299. https://doi.org/10.3390/act11100299

AMA Style

Chen X, Zhu Z, Ma T, Chang J, Chang X, Zang W. An Adaptive Dynamic Surface Technology-Based Electromechanical Actuator Fault-Tolerant Scheme for Blair Mine Hoist Wire Rope Tension Control System. Actuators. 2022; 11(10):299. https://doi.org/10.3390/act11100299

Chicago/Turabian Style

Chen, Xiao, Zhencai Zhu, Tianbing Ma, Jucai Chang, Xiangdong Chang, and Wanshun Zang. 2022. "An Adaptive Dynamic Surface Technology-Based Electromechanical Actuator Fault-Tolerant Scheme for Blair Mine Hoist Wire Rope Tension Control System" Actuators 11, no. 10: 299. https://doi.org/10.3390/act11100299

APA Style

Chen, X., Zhu, Z., Ma, T., Chang, J., Chang, X., & Zang, W. (2022). An Adaptive Dynamic Surface Technology-Based Electromechanical Actuator Fault-Tolerant Scheme for Blair Mine Hoist Wire Rope Tension Control System. Actuators, 11(10), 299. https://doi.org/10.3390/act11100299

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