Next Article in Journal
Diffusivity Measurement by Single-Molecule Recycling in a Capillary Microchannel
Previous Article in Journal
Challenges and Perspectives for Biosensing of Bioaerosol Containing Pathogenic Microorganisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Boundary Proportion Differential Control Method of Micro-Deformable Manipulator with Compensator Based on Partial Differential Equation Dynamic Model

School of Mechanical Engineering, Hebei University of Technology, Tianjin 300000, China
*
Author to whom correspondence should be addressed.
Micromachines 2021, 12(7), 799; https://doi.org/10.3390/mi12070799
Submission received: 31 May 2021 / Revised: 29 June 2021 / Accepted: 2 July 2021 / Published: 5 July 2021

Abstract

:
It is challenging to accurately judge the actual end position of the manipulator—regarded as a rigid body—due to the influence of micro-deformation. Its precise and efficient control is a crucial problem. To solve the problem, the Hamilton principle was used to establish the partial differential equation (PDE) dynamic model of the manipulator system based on the infinite dimension of the working environment interference and the manipulator space. Hence, it resolves the common overflow instability problem in the micro-deformable manipulator system modeling. Furthermore, an infinite-dimensional radial basis function neural network compensator suitable for the dynamic model was proposed to compensate for boundary and uncertain external interference. Based on this compensation method, a distributed boundary proportional differential control method was designed to improve control accuracy and speed. The effectiveness of the proposed model and method was verified by theoretical analysis, numerical simulation, and experimental verification. The results show that the proposed method can effectively improve the response speed while ensuring accuracy.

1. Introduction

With the rapid development of modern manipulators, the traditional rigid manipulators with large mass and margin cannot satisfy fast response and accurate positioning requirements. However, micro-deformable manipulators have lighter weight, lower energy consumption, and less inertia than traditional manipulators; moreover, they have high precision, high efficiency, high speed, high flexibility, high adaptability, and intelligence [1]. Their application range is wider than that of rigid manipulators. Therefore, the lightweight and dexterous micro-deformable manipulator dynamic model and its precise control have gradually become a hot research topic [2].
Researchers have done much work on the dynamic modeling and precise control of micro-deformable manipulators. In [3], to control a single-link flexible manipulator, a hybrid method combining sliding mode and H theory was proposed. Furthermore, a linear optimal damping controller was used in [4] to adjust the flexible boom vibration mode to a lower level. Based on the partial differential equation (PDE) dynamic model of the flexible manipulator system, some researchers used the adaptive boundary control method [5,6,7] to control the manipulator, while others used the neural network control method [8]. Several researchers have combined the adaptive boundary control method with the approximation or compensation results of radial basis function (RBF) neural network [9,10,11,12,13] to optimize the control performance of flexible systems. Additionally, the RBF neural network proportional differential (PD) control method of the flexible manipulator was studied in [14]. Moreover, in [15], based on a wavelet neural network, a dynamic surface control method was proposed. Furthermore, the iterative learning-based adaptive control methods were designed in [16,17] to obtain trajectory tracking and vibration reduction. The dynamic model of the flexible manipulator system was established, and the dynamic characteristics were analyzed in addition to different methods in [18,19,20]. In [21,22], the boundary control of flexible Timoshenko arm was studied and analyzed. Furthermore, the position control method and residual vibration of the flexible system were analyzed in [23,24,25]. Moreover, some optimized fuzzy control methods were proposed in [26,27,28]. However, most studies mentioned above simplified the micro-deformation manipulator into a relatively simple ideal model; these idealized models were not consistent with the actual situation [6,8,17,23]. Additionally, most studies ignored the infinite dimensionality of the micro-deformation manipulator. These lead to overflow and instability; moreover, change of the system space state cannot be described accurately [18,19,26]. The overflow problems were considered in a few PDE modeling papers [5,7,9,12] and RBF neural network was used to estimate interference of one end, but the other end was ignored [9]. Due to the partial interference being neglected, the established model was not accurate enough. This study completely considered the interference at both ends of the boundary while establishing the PDE dynamic model. The estimated interference results were added to the control law using two outputs RBF neural network for improvement and optimization purposes. Thus, the accuracy of the model was improved, and the model fitted the actual situation better.
Furthermore, this study fully considered the vibration and deformation of the micro-deformation manipulator and the joint micro-deformation. Moreover, we studied the distributed parameter boundary PD control method of RBF neural network compensator for the micro-deformation manipulator. First, the infinite-dimensional PDE dynamic model of the micro-deformable manipulator system was established. After that, the PD control method and stability analysis of the micro-deformable manipulator based on the RBF neural network compensator were introduced. Finally, numerical simulation and analysis of the control method were performed. Additionally, the response speed of the proposed method was increased by at least 30% compared with the adaptive boundary control method and RBF neural network method. Experimental verification was carried out for the proposed method at last.

2. Dynamic Modeling of Micro-Deformable Manipulator

An accurate dynamic model is a basis for achieving high-performance control. The control methods of micro-deformation manipulators were developed based on the ordinary differential equation (ODE) dynamic model in most studies. However, it cannot accurately describe the distributed parameter characteristics of the micro-deformable structure and may cause overflow instability problems. The PDE dynamic model can reflect the dynamic characteristics of the micro-deformed structure more accurately than the ODE dynamic model [5,9]. The Hamilton method was used to derive the PDE dynamic equation of the micro-deformable manipulator system [5,7,9,12]; moreover, the corresponding boundary conditions of the system were obtained. This process did not require complex force analysis of the micro-deformable manipulator system. The dynamic model can be derived directly by mathematical methods.
Remark 1.
The deformation of the micro-deformable manipulator is spatiotemporal, but the time variable t does not affect the calculation and derivation of the dynamic model. For simplicity, we omitted the time variable, t, in the function variable. For example,  y ( x , t ) is expressed as   y ( x ) ,   l ( x , t ) is expressed as   l ( x ) , and   θ ( t ) is expressed as   θ .
Remark 2. 
The usage and explanation of subscripts in this paper are as follows:
( * ) x = ( * ) x , ( * ) x x = 2 ( * ) x 2 , ( * ) x x x = 3 ( * ) x 3 , ( * ) x x x x = 4 ( * ) x 4 , ( * ) t = ( * ) t , ( * ) t t = 2 ( * ) t 2 , ( * ) x t t = ( * ) x t 2 .
We considered the plane motion single-rod micro-deformable manipulator as the research object (Figure 1); its cross-section was circular, and its radius, R, was 0.01 m. The relevant symbols of the manipulator are shown in Nomenclature. According to the theory and formulas of structural dynamics, mean normal stress ratio to shear stress can be calculated, as shown in Equation (1). Therefore, σ a v e τ s . In other words, bending deformation is the main deformation; therefore, the shear deformation can be neglected. Hence, we developed the dynamic model based on the Euler–Bernoulli beam theory.
σ a v e τ s = σ max ( 2 / 3 ) τ s = ( M L y max / I z ) ( 2 / 3 ) M / A = 8 M / 3 π R 3 M / π R 2 = 8 3 R 266.67
The offset of any point, [ x , y ( x ) ] , in the follow-up coordinate system, x O y , from the micro-deformation manipulator in the inertial coordinate system, X O Y , is approximately expressed as l ( x ) .
l ( x ) = x θ + y ( x )
The manipulator is regarded as an Euler–Bernoulli beam here and clamped to a motor at x . The natural boundary conditions can be expressed as:
y ( 0 ) = y x ( 0 ) = 0
Using Equations (2) and (3), we get:
l ( 0 ) = 0 , l x ( 0 ) = θ , m l ( x ) x m = m y ( x ) x m , m 2
According to Hamilton principle [5,6,29,30], the PDE dynamic equations of the micro-deformable manipulator can be developed for every 0 t 1 < t < t 2 , as shown in Equation (5):
t 1 t 2 ( δ W k δ W p + δ W n c ) d t = 0
where δ W k , δ W p and δ W n c represent the kinetic energy, potential energy, and the variation of work done by non-conservative force, respectively; t i denotes a moment.
The total kinetic energy of the system can be obtained considering the rotational kinetic energy of the micro-deformation manipulator joint, the kinetic energy of the manipulator, and the kinetic energy of the load, as shown in Equation (6).
W k = 1 2 I m θ t 2 + 1 2 0 L ρ l t 2 ( x ) d x + 1 2 m l t 2 ( L )
The potential energy of the system is expressed as:
W p = 1 2 0 L E I y x x 2 ( x ) d x
The work done by the non-conservative force of the system is expressed as:
W n c = τ θ + p 1 θ + M l ( L ) + p 2 l ( L )
Substituting Equations (4), (6)–(8) into Equation (5), we get:
t 1 t 2 ( δ W k δ W p + δ W n c ) = t 1 t 2 0 L K 1 δ l ( x ) d x d t t 1 t 2 K 2 δ l x ( 0 ) d t t 1 t 2 K 3 δ l ( L ) d t t 1 t 2 K 4 δ l x ( L ) d t = 0
where:
{ K 1 = ρ l t t ( x ) + E I l x x x x ( x ) K 2 = I m l x t t ( 0 ) E I l x x ( 0 ) ( τ + p 1 ) K 3 = m l t t ( L ) E I l x x x ( L ) ( M + p 2 ) K 4 = E I l x x ( L )
Each monomial in Equation (9) is linearly independent because δ l ( x ) , δ l x ( 0 ) , δ l ( L ) , δ l x ( L ) are independent variables. Hence, K 1 = K 2 = K 3 = K 4 = 0 . Therefore, the PDE dynamic model of the micro-deformation manipulator system is obtained as Equation (10).
{ ρ l t t ( x ) = E I l x x x x ( x ) τ = I m l x t t ( 0 ) E I l x x ( 0 ) p 1 M = m l t t ( L ) E I l x x x ( L ) p 2 l x x ( L ) = 0
where: l t t ( x ) = x θ t t + y t t ( x ) .
The dynamic model of the manipulator system was developed from the mathematical model perspective. The motion characteristics are related to time and position; therefore, the micro-deformable manipulator is essentially a distributed parameter system. Hence, the distributed parameter model was established based on the PDE equation. Furthermore, the corresponding control method adopted the distributed parameter boundary PD control, which can effectively obtain the micro-deformable system control. The boundary control only needs a small number of thrusters to achieve a better control effect than the discrete distributed control.

3. RBF Neural Network Distributed Boundary PD Control Method

The paper [5] proposed an adaptive boundary control method suitable for the PDE dynamic model. The method was simple and accurate, but the response time remained long. Compared with other machine learning algorithms, RBF neural network has the ability of parallel information processing, stronger computing power, and faster running speed. Moreover, RBF neural network can avoid the local minimum problem. Therefore, the control scheme based on RBF neural network is more suitable for the requirement of real-time control. Thus, RBF neural network compensator was used to improve the adaptive boundary control method described in [5]. The joint position was adjusted, and the vibration was weakened by designing the Lyapunov function and boundary PD control law (Figure 2).

3.1. Design of Distributed Boundary PD Control Law Based on RBF Neural Network Compensator

The error function is defined as e = θ θ d , where θ d denotes the expectation angle. Consider that the uncertainty interference p 1 and p 2 of the actual model is unknown. The RBF neural network is used to estimate the uncertainty interference p 1 and p 2 . The estimated values are p ^ 1 and p ^ 2 . Furthermore, the input vector of the neural network is considered as x = [ e T e ˙ T θ d T θ ˙ d T θ ¨ d T ] . Moreover, the ideal RBF neural network algorithm is shown in Equation (11).
{ φ i = exp ( x c i 2 2 b i 2 ) , i = 1 , 2 , , m p 1 = W 1 * φ 1 ( x ) + γ 1 , p 2 = W 2 * φ 2 ( x ) + γ 2
where φ i = [ φ 1 , φ 2 , , φ m ] T denotes the output vector of hidden layer obtained by Gaussian function, m represents the number of neurons in the hidden layer, c i is the coordinate vector value of the Gaussian function’s center point, b i represents the width of Gaussian function, W i * denotes the ideal weight matrix of the neural network, and γ i is the network estimation error.
The output of the RBF neural network is shown in Equation (12).
{ p ^ 1 = W ^ 1 * φ 1 ( x ) p ^ 2 = W ^ 2 * φ 2 ( x )
where W ^ i * represents the estimated weight of the neural network for unknown parameter estimation.
The errors of the above estimation results are defined in Equation (13).
{ W ˜ 1 * = W 1 * W ^ 1 * , W ˜ 2 * = W 2 * W ^ 2 * p ˜ 1 = p 1 p ^ 1 , p ˜ 2 = p 2 p ^ 2
Equations (11)–(13) suggest that
{ p ˜ 1 = W ˜ 1 * φ 1 ( x ) + γ 1 p ˜ 2 = W ˜ 2 * φ 2 ( x ) + γ 2
Furthermore, the estimated results, p ^ 1 and p ^ 2 , of the RBF neural network were added into the control law as a compensator to compensate for the external interference, in order to obtain the angular response of the micro-deformation manipulator and suppress its deformation and vibration effectively. The RBF-based boundary PD control law is presented as:
{ τ = k p e k d e t + p ^ 1 M = k u α + m l x x x t ( L ) + p ^ 2
where: k p , k d , and k are control gain values and k p > 0 , k d > 0 , k > 0 , u α = l t ( L ) l x x x ( L ) , e t = θ t , e t t = θ t t .
Remark 3. 
All parameters mentioned in the designed control law (Equation (15)) are measurable or computable. We can use the position sensor, tachometer, laser displacement sensor, and shear force sensor to get the measurements,  θ , θ t , l ( x ) , and l x x x ( x ) , respectively. e , e t , e t t , l t ( x ) , and l x x x t ( x ) can be calculated from the measurements.

3.2. Stability Analysis Based on Lyapunov Function

Lemma 1. 
For ϕ 1 ( x , t ) , ϕ 2 ( x , t ) R , x [ 0 , L ] , t [ 0 , ) , the following inequality holds:
ϕ 1 ( x , t ) ϕ 2 ( x , t ) | ϕ 1 ( x , t ) ϕ 2 ( x , t ) | ϕ 1 2 ( x , t ) + ϕ 2 2 ( x , t ) | ϕ 1 ( x , t ) ϕ 2 ( x , t ) | 1 λ ϕ 1 2 ( x , t ) + λ ϕ 2 2 ( x , t ) , ( λ > 0 )
Lemma 2. 
For  p ( x , t ) R , x [ 0 , L ] , t [ 0 , ) , if p ( 0 , t ) = 0 , t [ 0 , ) , then p 2 ( x , t ) L 0 L p x 2 ( x , t ) d x , x [ 0 , L ] . Similarly, if p x ( 0 , t ) = 0 , t [ 0 , ) , then p x 2 ( x , t ) L 0 L p x x 2 ( x , t ) d x , x [ 0 , L ] .
Lemma 3. 
For  V : [ 0 , ) R , t t 0 0 , if V ˙ η V + g , then:
V ( t ) e η ( t t 0 ) V ( t 0 ) + t 0 t e η ( t s ) g ( s ) d s , η > 0 .
Theorem 1.
The closed-loop system described by Equation (10) is asymptotically stable, with the proposed RBF neural network compensator (Equation (12)) and control law (Equation (15)), namely,  θ θ d , θ t 0 , y ( x ) 0 , y t ( x ) 0 , as x [ 0 , L ] , t .
Proof of Theorem 1.
The Lyapunov function is defined to prove the stability of the controller, as in Equation (16).
H ( t ) = W 1 + W 2 + W α
□ 
where:
{ W 1 = 1 2 0 L ρ l t 2 ( x ) d x + 1 2 E I 0 L y x x 2 ( x ) W 2 = 1 2 I m e t 2 + 1 2 k p e 2 + 1 2 m u α 2 p 1 e p 2 t 1 t 2 u α d t W α = α ρ 0 L x l t ( x ) l e x ( x ) d x + α I m e e t
In Equation (17), W 1 represents the inhibition index bending deformation of the micro-deformation manipulator. The first two items in W 2 represent the control error-index, and the rest of the items are the auxiliary terms. Moreover, W α is the auxiliary part and α denotes a small positive constant. Furthermore,
l e ( x ) = x e + y ( x ) ,   l e x ( x ) = e + y x ( x ) ,   l e x x ( x ) = y x x ( x ) = l x x ( x )
According to Lemmas 1 and 2, we obtain:
| W α | α ρ L 0 L l t 2 ( x ) d x + 2 α ρ L 2 e 2 + 2 α ρ L 3 0 L l x x 2 ( x ) d x + α I m ( e 2 + e t 2 ) α 1 ( W 1 + W 2 )
where: α 1 = max ( 2 α L , 2 α ρ L 3 / E I , 2 ( α I m + 2 α ρ L 2 ) / k p , 2 α ) .
Therefore,
α 1 ( W 1 + W 2 ) W α α 1 ( W 1 + W 2 )
Take 0 < α1 < 1, i.e., 0 < max ( 2 α L , 2 α ρ L 3 / E I , 2 ( α I m + 2 α ρ L 2 ) / k p , 2 α ) < 1 , then α can be designed as 0 < α < 1 ( 2 L , 2 ρ L 3 / E I , 2 ( I m + 2 ρ L 2 ) / k p , 2 ) . We define 1 > α 2 = 1 α 1 > 0 2 > α 3 = 1 + α 1 > 1 ; hence,
0 α 2 ( W 1 + W 2 ) H ( t ) α 3 ( W 1 + W 2 )
According to inequality in Equation (21), the Lyapunov function is positively definite. Hence,
H t ( t ) = ( W 1 ) t + ( W 2 ) t + ( W α ) t
In Equation (22),
{ ( W 1 ) t = 0 L ρ l t ( x ) l t t ( x ) d x + E I 0 L y x x ( x ) y x x t ( x ) ( W 2 ) t = I m e t e t t + k p e e t + m u α ( u α ) t p 1 e t p 2 u α ( W α ) t = ( W α 1 ) t + ( W α 2 ) t + ( W α 3 ) t
In Equation (23),
{ ( W α 1 ) t = α ρ 0 L x l t t ( x ) l e x ( x ) d x ( W α 2 ) t = α ρ 0 L x l t ( x ) l t e x ( x ) d x ( W α 3 ) t = α I m ( e t 2 + e e t t )
Substituting Equation (10) into Equation (23) and combining Equation (4), ( W 1 ) t can be rewritten as:
( W 1 ) t = E I y x x x ( L ) l t ( L ) E I y x x ( 0 ) θ t = E I l x x ( 0 ) e t E I l x x x 2 ( L ) E I l x x x ( L ) u α
Using Equations (10) and (15) and combining Equations (23) and (25), we obtain:
( W 1 ) t + ( W 2 ) t = E I l x x ( 0 ) e t E I l x x x 2 ( L ) E I l x x x ( L ) u α + I m e t e t t + k p e e t + m u α ( u α ) t p 1 e t p 2 u α = e t ( E I l x x ( 0 ) + I m e t t + k p e p 1 ) + u α ( E I l x x x ( L ) + m ( u α ) t p 2 ) E I l x x x 2 ( L ) = e t ( τ + k p e p 1 ) + u α ( M m l x x x t ( L ) p 2 ) E I l x x x 2 ( L ) = k d e t 2 k u α 2 E I l x x x 2 ( L )
According to Lemmas 1–3, combining Equations (10) and (18), and using the method of integral by parts, we deduce:
( W α 1 ) t = α E I L l e x ( L ) l x x x ( L ) 3 2 α E I 0 L l x x 2 ( x ) d x α E I e l x x ( 0 ) ( 3 2 α 2 α L 2 2 α L 3 E I ) 0 L E I l x x 2 ( x ) d x + α E I L l x x x 2 ( L ) α E I e l x x ( 0 ) + ( 2 α E I L + 2 α L 2 ) e 2
( W α 2 ) t = 1 2 α ρ L l t 2 ( L ) 1 2 α ρ 0 L l t 2 ( x ) d x
( W α 3 ) t = α I m e t 2 + α e E I l x x ( 0 ) α k p e 2 k d α e e t ( α I m + k d α ) e t 2 ( α k p k d α ) e 2 + α e E I l x x ( 0 )
We obtain Equation (30) using Equations (23), (27)–(29).
( W α ) t ( 3 2 α 2 α L 2 2 α L 3 E I ) 0 L E I l x x 2 ( x ) d x + α E I L l x x x 2 ( L ) + ( 2 α E I L + 2 α L 2 ) e 2 + 1 2 α ρ L l t 2 ( L ) 1 2 α ρ 0 L l t 2 ( x ) d x + ( α I m + k d α ) e t 2 ( α k p k d α ) e 2
Substituting Equations (26) and (30) into Equation (22), we get Equation (31).
H t ( t ) = ( W 1 ) t + ( W 2 ) t + ( W α ) t k d e t 2 k u α 2 E I l x x x 2 ( L ) ( 3 2 α 2 α L 2 2 α L 3 E I ) 0 L E I l x x 2 ( x ) d x + α E I L l x x x 2 ( L ) + ( 2 α E I L + 2 α L 2 ) e 2 + 1 2 α ρ L l t 2 ( L ) 1 2 α ρ 0 L l t 2 ( x ) d x + ( α I m + k d α ) e t 2 ( α k p k d α ) e 2 = ( 3 2 α 2 α L 2 2 α L 3 E I ) 0 L E I l x x 2 ( x ) d x 1 2 α 0 L ρ l t 2 ( x ) d x ( k d α I m k d α ) e t 2 ( α k p k d α 2 α E I L 2 α L 2 ) e 2 k u α 2 + 1 2 α ρ L l t 2 ( L ) ( E I α E I L ) l x x x 2 ( L )
The inequality E I α E I L > 1 2 α ρ L can be established by choosing α, which guarantees that:
1 2 α ρ L l t 2 ( L ) ( E I α E I L ) l x x x 2 ( L ) η 0 ( l t ( L ) l x x x ( L ) ) 2 = η 0 u α 2
where: η 0 > max ( η 1 , η 1 η 2 η 2 η 1 ) , η 1 = 1 2 α ρ L , η 2 = E I α E I L .
According to Equations (31) and (32), we obtain Equation (33).
H t ( t ) ( 3 2 α 2 α L 2 2 α L 3 E I ) 0 L E I l x x 2 ( x ) d x 1 2 α 0 L ρ l t 2 ( x ) d x ( k d α I m k d α ) e t 2 ( α k p k d α 2 α E I L 2 α L 2 ) e 2 ( k η 0 ) u α 2 λ 0 ( W 1 + W 2 ) λ 0 H ( t ) α 3 = λ H ( t )
The following constraints must be met to ensure the validity of Equation (33).
ε 1 = 3 2 α 2 α L 2 2 α L 3 E I > 0 , ε 2 = 1 2 α > 0 , ε 3 = k d α I m k d α > 0 , ε 4 = α k p k d α 2 α E I L 2 α L 2 > 0 , ε 5 = k η 0 > 0 , min ( 2 ε 1 , 2 ε 2 , 2 ε 3 I m , 2 ε 4 k p , 2 ε 5 m ) λ 0 > 0 , λ = λ 0 α 3
The solution for the inequality in Equation (33) is H ( t ) H ( 0 ) e λ t . Therefore, the Lyapunov function H ( t ) is close to 0 if the condition H ( 0 ) is bounded.
According to Equation (21), W 1 + W 2 0 , therefore, e 0 and e t 0 , meaning θ θ d and θ t 0 . Additionally, l t ( x ) 0 , y t ( x ) 0 . Moreover, according to Lemma 2, y x 2 ( x ) L 0 L y x x 2 ( x ) d x , y 2 ( x ) L 0 L y x 2 ( x ) d x . Therefore, we can further deduce that y ( x ) 0 . Moreover, the closed-loop system tends to be stable.

4. Numerical Simulation Analysis

Three simulation tests were performed to verify the performance and effect of the proposed control method. During simulation, the discrete time was Δ t = 5 × 10 4 and discrete distance was Δ x = 0.01 m . The system parameters and controller coefficients were set as shown in Table 1. The initial state of the system and the initial interference compensation were set to 0. The number of neurons in the input layer, hidden layer, and output layer of the RBF neural network was 5, 5, and 2, respectively. According to the actual range of the input layer x , the parameters of the Gaussian basis function of the hidden layer were set to:  c j = [ 2.5 1.5 0 1.5 2.5 2 1 0 1 2 1.5 0.75 0 0.75 1.5 1 0.5 0 0.5 1 0.5 0.25 0 0.25 0.5 ] , b j = 0.5 . In order to ensure reasonable results, the parameters and coefficients in the simulation scenarios 1–3 should be consistent.
Scenario 1: with the adaptive boundary control method [5].
Scenario 2: with the RBF neural network control method [8].
Scenario 3: with the proposed control method.
Remark 4
. The design description of coefficients k p , k d , and in control law is shown in Appendix A.
Remark 5. 
An appropriately large control gain k p  can ensure the good tracking performance of the system, but if too large, the system will become invalid. Therefore, in practical operation, the system performance and the actuator saturation should be considered simultaneously when we design the coefficient.
The simulation results of the micro-deformable manipulator with the proposed control method are shown in Figure 3c, Figure 4, and Figure 5c, Figure 6c and Figure 7c; results with the adaptive boundary control method are shown in Figure 3a and Figure 5a, Figure 6a and Figure 7a; Figure 3b and Figure 5b, Figure 6b and Figure 7b show the simulation results with the RBF neural network control method.
Figure 3 shows the angle tracking and angular velocity response results of micro-deformable manipulator with different methods. As shown in Figure 3a,b, the adaptive boundary control method and RBF neural network control method affect the control performance due to the large error and long response time. On the contrary, the proposed control method can accurately adjust the joint angle to the expected value within 6 s, as shown in Figure 3c. The response time comparisons of three methods are shown in Table 2. As shown in Table 2, the response speed is improved by at least 30% than that of the other two methods. Therefore, the compensation results of the RBF neural network are satisfactory, as shown in Figure 4.
Figure 5 and Figure 6 show the elastic deformation and deformation rate of manipulator successively. The statistical results of deformation rate are shown in Table 2. As shown in Figure 5a,b, the elastic deformation of the micro-deformable manipulator is remarkable, and it is not suppressed obviously. Conversely, in Figure 5c, the elastic deformation reaches the peak in 5 s and is eliminated obviously in 7 s. Furthermore, compared with Scenarios 1 and 2, the deformation rate of the manipulator is reduced by an order of magnitude with the proposed method (Table 2).
Figure 7 depicts the corresponding control input with different methods. Furthermore, the related statistical results are shown in Table 2. In Figure 7a,b, the control input of the adaptive boundary control method and RBF neural network control method is unacceptable because of sizeable fluctuation and overshoot. Conversely, in Figure 7c, the fluctuation and overshoot are nearly weakened within 2 s.
The three simulations confirmed that the proposed method performs excellently in trajectory tracking and vibration and deformation suppression of the micro-deformable manipulator under the premise of considering the double boundary interference. Through the analysis of the above simulation results, the proposed method is superior to the adaptive boundary control method and RBF neural network control method. Compared with the other two methods, the response time of the proposed method is reduced by at least 30%. The deformation of the manipulator is restrained to a great extent.

5. Experimental Tests

A simple slender single link micro-deformable manipulator (L = 1 m, R = 0.01 m which are consistent with the numerical analysis) was used for the experimental test in this paper. Air experiment was carried out in the laboratory. The experimental apparatus is shown in Figure 8a. In each experiment, the real-time position, current, speed, and other information of the micro-deformable manipulator were sent to the computer through the serial port, and then the computer, as the host computer, sent the order to the driver which drives the thruster to control the motion of the manipulator. The schematic diagram of the experiment is shown in Figure 8b. The thruster here was a MOTEC DC servo motor of DSEM-V241230E60LN type. It was driven by MOTEC DC servo driver of ARES-80-15-E-A0 type. The host computer here was a notebook computer with Intel Core i7-5500u 2.5 GHz CPU and 12 GB of RAM, running under Windows 10 operating system. The control software was developed by Visual Studio 2017. The manipulator system was powered by 24 V DC power supply.
In each air experiment, the parameters of the motor were set to pulse/2π rad in advance. The manipulator was controlled to rotate 1592 pulses per time, which was about 0.5 rad, and swing back and forth every time. In order to ensure the accuracy of data analysis, the angle set in the experiment was consistent with that in the simulation. The data of one cycle were collected for analysis (Figure 9).
Figure 9 depicts the change of the torque values of those collected from the experiment and simulated with different methods. The errors of experiment and three simulations are recorded in Table 3. In Figure 9b,c, the other two methods are not satisfactory due to their sizeable errors. As shown in Figure 9d, the variation of the numerical torque value of the proposed method is similar to that of the experimental torque value. The maximum error did not exceed 0.7 N·m as shown in Table 3, which verified the correctness of the numerical results of the proposed method.

6. Conclusions

This paper studied the PDE dynamic model; based on this model, the distributed boundary PD controller with infinite-dimensional RBF neural network compensator was proposed. The RBF neural network was used to compensate for double boundary interference and uncertain interference. The boundary control can be used to obtain the angle tracking of the micro-deformable manipulator. At the same time, it effectively weakened the vibration generated by the micro-deformation manipulator during the movement. The response time of the proposed method was reduced by at least 30% compared with the adaptive boundary controller and the RBF neural network controller. The error between the simulated results and experimental results was no more than 0.7 N·m, which verified the accuracy of numerical results and the feasibility of the proposed method. Since the unknown interference distribution in the external space is random and unpredictable, it is essential to subdivide the spatially distributed interference further and study the corresponding compensation control method in the future.

Author Contributions

Conceptualization, X.P. and Y.T.; methodology, X.P.; software, X.P. and R.S.; validation, X.P., Y.T. and R.S.; formal analysis, X.P.; investigation, X.P. and M.Z.; resources, X.P. and M.Z.; data curation, R.S.; writing—original draft preparation, X.P.; writing—review and editing, Y.T.; visualization, X.P.; supervision, M.Z.; project administration, M.Z.; funding acquisition, M.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the National key R&D project of China (grant number: 2018YFB1309400 and 2018YFB1309401), the Science and Technology Research Project of Hebei Province (grant number QN2018090) and the Natural Science Foundation of Hebei Province (grant number E2021202032).

Acknowledgments

The computational work and experimental test for this research were performed at the Special Robot Laboratory of Hebei University of Technology, China.

Conflicts of Interest

We state that we do not use this work commercially and have not published this article. The authors declare no conflict of interest.

Nomenclature

L Length of micro-deformable manipulator, m
R Radius of micro-deformable manipulator, m
ρ Linear density of micro-deformable manipulator, kg/m
I m Center moment of inertia, kg·m2
E I Bending stiffness of uniform beam, N·m2
τ Motor control input torque at initial end point, N·m
M Control input force of end load, N
m Terminal load mass, kg
θ Joint turning angle (excluding deformation),
y(x)Elastic deformation of manipulator at point x, m
l(x)Offset of micro-deformable manipulator in the inertial coordinate system, m
F End force in the free-body diagram, N
T Moment obtained by force analysis, N·m
σ a v e Mean normal stress, N/m2
τ s Shear stress, N/m2
A Cross section area of micro-deformable manipulator, m2
p 1 Boundary interference of one side
p 2 Boundary interference of another side
p ^ 1 The estimated value of p 1
p ^ 2 The estimated value of p 2
k p The coefficient of control law
k d The coefficient of control law
k The coefficient of control law
λ Flexibility of micro-deformable manipulator

Appendix A

Design explanation of k p , k d , and k .
1.
According to the engineering trial method, the coefficient k p is taken as 60.
2.
According to Equations (19)–(21), 0 < α < 1 ( 2 L , 2 ρ L 3 / E I , 2 ( I m + 2 ρ L 2 ) / k p , 2 ) .
3.
From η 1 = 1 2 α ρ L , η 2 = E I α E I L , we obtain that η 0 > max ( η 1 , η 1 η 2 η 2 η 1 ) .
4.
Take k d = 30 according to the condition: ε 3 = k d α I m k d α > 0 ; take k = 20 according to the condition: ε 5 = k η 0 > 0 .
5.
Verifying other conditions in Equation (34):
ε 1 = 3 2 α 2 α L 2 2 α L 3 E I > 0 ε 2 = 1 2 α > 0 ε 4 = α k p k d α 2 α E I L 2 α L 2 > 0 ε 5 = k η 0 > 0 .
6.
Take α 3 = 1.5 according to the condition: 2 > α 3 > 1 .
7.
From min ( 2 ε 1 , 2 ε 2 , 2 ε 3 I m , 2 ε 4 k p , 2 ε 5 m ) λ 0 > 0 , we obtain λ = λ 0 α 3 .

References

  1. Alandoli, E.A.; Lee, T.S.; Lin, Y.J.; Vijayakumar, V. Dynamic Model and Intelligent Optimal Controller of Flexible Link Manipulator System with Payload Uncertainty. Arab. J. Sci. Eng. 2021. [Google Scholar] [CrossRef]
  2. Jing, Z.L.; Xu, Q.M.; Huang, J.Z. A review on kinematic analysis and dynamic stable control of space flexible manipulators. Aerosp. Syst. 2019, 2, 1–14. [Google Scholar] [CrossRef] [Green Version]
  3. Kharabian, B.; Mirinejad, H. Hybrid sliding mode/H-infinity control approach for uncertain flexible manipulators. IEEE Access 2020, 8, 170452–170460. [Google Scholar] [CrossRef]
  4. Lei, R.H.; Chen, L. Adaptative fault-tolerant and vibration-suppressed hybrid control for flexible space manipulator. J. Astronaut. 2020, 41, 472–482. [Google Scholar]
  5. Zhang, L.J.; Liu, J.K. Adaptive boundary control for flexible two-link manipulator based on partial differential equation dynamic model. IET Control Theory Appl. 2013, 7, 43–51. [Google Scholar] [CrossRef]
  6. He, W.; Ge, S.S.; How, B.V.E.; Choo, Y.S.; Hong, K.-S. Robust adaptive boundary control of a flexible marine riser with vessel dynamics. Automatica 2011, 47, 722–732. [Google Scholar] [CrossRef]
  7. He, W.; He, X.; Zou, M.; Li, H. PDE model-based boundary control design for a flexible robotic manipulator with input backlash. IEEE Trans. Control Syst. Technol. 2019, 27, 790–797. [Google Scholar] [CrossRef]
  8. Gao, H.; He, W.; Zhou, C.; Sun, C. Neural network control of a two-link flexible robotic manipulator using assumed mode method. IEEE Trans. Ind. Inform. 2019, 15, 755–765. [Google Scholar] [CrossRef]
  9. Li, L.; Liu, J.K. Neural-network-based adaptive fault-tolerant vibration control of single-link flexible manipulator. Trans. Inst. Meas. Control 2020, 42, 430–438. [Google Scholar] [CrossRef]
  10. Yang, H.J.; Liu, J.K. An adaptive RBF neural network control method for a class of nonlinear systems. IEEE-CAA J. Autom. Sin. 2018, 5, 457–462. [Google Scholar] [CrossRef]
  11. Zhang, S.; Dong, Y.; Ouyang, Y.; Yin, Z.; Peng, K. Adaptive Neural Control for Robotic Manipulators with Output Constraints and Uncertainties. IEEE Trans. Neural Netw. Learn. Syst. 2018, 29, 5554–5564. [Google Scholar] [CrossRef] [PubMed]
  12. He, C.; Zhang, F.; Jiang, J.H. Adaptive Boundary Control of Flexible Manipulators with Parameter Uncertainty Based on RBF Neural Network. Shock Vib. 2020, 2020. [Google Scholar] [CrossRef]
  13. Ren, Y.; Zhao, Z.; Zhang, C.; Yang, Q.; Hong, K.-S. Adaptive Neural-Network Boundary Control for a Flexible Manipulator with Input Constraints and Model Uncertainties. IEEE Trans. Cybern. 2020. [Google Scholar] [CrossRef]
  14. Zhang, Y.Y.; Kim, D.G.; Zhao, Y.D.; Lee, J.M. PD Control of a Manipulator with Gravity and Inertia Compensation Using an RBF Neural Network. Int. J. Control Autom. Syst. 2020, 18, 3083–3092. [Google Scholar] [CrossRef]
  15. Gu, Y.K.; Liu, H. Adaptative control of flexible-joint manipulator based on dynamics surface and neural networks. J. Huazhong Univ. Sci. Tech. (Nat. Sci. Ed.) 2018, 46, 64–69. [Google Scholar]
  16. Yang, Y.; Liu, Z.; Ma, G. Adaptive distributed control of a flexible manipulator using an iterative learning scheme. IEEE Access 2019, 7, 145934–145943. [Google Scholar] [CrossRef]
  17. Zhou, X.Y.; Wang, H.P.; Tian, Y.; Zheng, G. Disturbance observer-based adaptive boundary iterative learning control for a rigid-flexible manipulator with input backlash and endpoint constraint. Int. J. Adapt. Control Signal Process. 2020, 34, 1220–1241. [Google Scholar] [CrossRef]
  18. Wang, H.; Fu, B.C.; Xue, B.; Yang, C.L. Dynamics analysis of 6-DOF manipulators with flexible joints. China Mech. Eng. 2016, 27, 1096–1101. [Google Scholar]
  19. Wei, Y.K.; Jiao, G.T. Flexible robot arm dynamics analysis based on MATLAB. Mech. Eng. Autom. 2016, 1, 33–36. [Google Scholar]
  20. Liu, X.L.; Huang, Y.P.; Cui, P.J.; Xu, Z.X. Modeling and dynamic characteristic analysis of flexible robotic arm. Noise Vib. Control 2014, 34, 7–11. [Google Scholar]
  21. Endo, T.; Sasaki, M.; Matsuno, F. Contact-force control of a flexible Timoshenko arm. IEEE Trans. Autom. Control 2017, 62, 1004–1009. [Google Scholar] [CrossRef]
  22. Han, F.; Jia, Y. Boundary control and exponential stability of a flexible Timoshenko beam manipulator with measurement delays. IET Control Theory Appl. 2020, 14, 499–510. [Google Scholar] [CrossRef]
  23. Meng, Q.X.; Lai, X.Z.; Yan, Z.; Wang, Y.W.; Wu, M. Position control without residual vibration for a two-link rigid-flexible manipulator. Control Theory Appl. 2020, 37, 620–628. [Google Scholar]
  24. Sun, L.; Yin, W.; Wang, M.; Liu, J. Position Control for Flexible Joint Robot Based on Online Gravity Compensation with Vibration Suppression. IEEE Trans. Ind. Electron. 2018, 65, 4840–4848. [Google Scholar] [CrossRef]
  25. Pradhan, S.K.; Subudhi, B. Position control of a flexible manipulator using a new nonlinear self-tuning PID controller. IEEE/CAA J. Autom. Sin. 2020, 7, 136–149. [Google Scholar] [CrossRef]
  26. Sun, C.; Gao, H.; He, W.; Yu, Y. Fuzzy neural network control of a flexible robotic manipulator using assumed mode method. IEEE Trans. Neural Netw. Learn. Syst. 2018, 29, 5214–5227. [Google Scholar] [CrossRef]
  27. Dong, F.; Zhao, X.; Han, J.; Chen, Y. Optimal fuzzy adaptive control for uncertain flexible joint manipulator based on D-operation. IET Control Theory Appl. 2018, 12, 1286–1298. [Google Scholar] [CrossRef]
  28. Zhang, C.L.; Yang, T.; Sun, N.; Fang, Y.C. An Adaptive Fuzzy Control Method of Single-Link Flexible Manipulators with Input Dead-Zones. Int. J. Fuzzy Syst. 2020, 22, 2521–2533. [Google Scholar] [CrossRef]
  29. Meirovitch, L. Fundamentals of Vibrations; McGrow Hill: New York, NY, USA, 2002. [Google Scholar]
  30. Rahn, C.D. Mechatronic Control of Distributed Noise and Vibration; Springer: New York, NY, USA, 2001. [Google Scholar]
Figure 1. A typical single link micro-deformable manipulator: (a) coordinate diagram; (b) free-body diagram of the system.
Figure 1. A typical single link micro-deformable manipulator: (a) coordinate diagram; (b) free-body diagram of the system.
Micromachines 12 00799 g001
Figure 2. Boundary PD control principle with RBF neural network compensator.
Figure 2. Boundary PD control principle with RBF neural network compensator.
Micromachines 12 00799 g002
Figure 3. The angle tracking and response results of different methods: (a) adaptive boundary method; (b) RBF neural network method; (c) the proposed method.
Figure 3. The angle tracking and response results of different methods: (a) adaptive boundary method; (b) RBF neural network method; (c) the proposed method.
Micromachines 12 00799 g003
Figure 4. RBF neural network compensation results: (a) compensation results of p1; (b) compensation results of p2.
Figure 4. RBF neural network compensation results: (a) compensation results of p1; (b) compensation results of p2.
Micromachines 12 00799 g004
Figure 5. The deformation results of different methods: (a) adaptive boundary method; (b) RBF neural network method; (c) the proposed method.
Figure 5. The deformation results of different methods: (a) adaptive boundary method; (b) RBF neural network method; (c) the proposed method.
Micromachines 12 00799 g005
Figure 6. The deformation rate results of different methods: (a) adaptive boundary method; (b) RBF neural network method; (c) the proposed method.
Figure 6. The deformation rate results of different methods: (a) adaptive boundary method; (b) RBF neural network method; (c) the proposed method.
Micromachines 12 00799 g006
Figure 7. The control input τ and M results of different methods: (a) adaptive boundary method; (b) RBF neural network method; (c) the proposed method.
Figure 7. The control input τ and M results of different methods: (a) adaptive boundary method; (b) RBF neural network method; (c) the proposed method.
Micromachines 12 00799 g007
Figure 8. Experimental apparatus and schematic diagram: (a) experimental apparatus of the micro-deformable manipulator; (b) the schematic diagram of the experiment.
Figure 8. Experimental apparatus and schematic diagram: (a) experimental apparatus of the micro-deformable manipulator; (b) the schematic diagram of the experiment.
Micromachines 12 00799 g008
Figure 9. Comparison of τ  between experiment and simulations: (a) comparison of overall results; (b) comparison of partial enlarged results of adaptive boundary method; (c) comparison of partial enlarged results of the RBF neural network method; (d) comparison of partial enlarged results of the proposed method.
Figure 9. Comparison of τ  between experiment and simulations: (a) comparison of overall results; (b) comparison of partial enlarged results of adaptive boundary method; (c) comparison of partial enlarged results of the RBF neural network method; (d) comparison of partial enlarged results of the proposed method.
Micromachines 12 00799 g009
Table 1. Numerical simulation marks of micro-deformable manipulator system.
Table 1. Numerical simulation marks of micro-deformable manipulator system.
MarkEImImLρkpkdkλ
Value30.150.11.00.21603020200
UnitN·m2kg kg·m2m kg/m ////
Table 2. The statistical results of the simulation data of the proposed method and the comparison methods.
Table 2. The statistical results of the simulation data of the proposed method and the comparison methods.
ComparisonsAdaptive Boundary MethodRBF Neural Network MethodThe Proposed Method
Response time (s)796
Maximum of
deformation rate (m/s)
−0.2132−0.1396−0.0269
Mean value of M (N·m)0.29640.72330.3333
Mean value of τ (N·m)−0.3334−0.8040−0.3750
Table 3. The errors between experiment and three simulations.
Table 3. The errors between experiment and three simulations.
Errors (N·m) Adaptive Boundary MethodRBF Neural Network MethodThe Proposed Method
Mean value1.31906.67500.3008
Standard deviation0.92038.50600.1604
Maximum1.762020.9300.6644
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pei, X.; Tian, Y.; Zhang, M.; Shi, R. The Boundary Proportion Differential Control Method of Micro-Deformable Manipulator with Compensator Based on Partial Differential Equation Dynamic Model. Micromachines 2021, 12, 799. https://doi.org/10.3390/mi12070799

AMA Style

Pei X, Tian Y, Zhang M, Shi R. The Boundary Proportion Differential Control Method of Micro-Deformable Manipulator with Compensator Based on Partial Differential Equation Dynamic Model. Micromachines. 2021; 12(7):799. https://doi.org/10.3390/mi12070799

Chicago/Turabian Style

Pei, Xiangli, Ying Tian, Minglu Zhang, and Ruizhuo Shi. 2021. "The Boundary Proportion Differential Control Method of Micro-Deformable Manipulator with Compensator Based on Partial Differential Equation Dynamic Model" Micromachines 12, no. 7: 799. https://doi.org/10.3390/mi12070799

APA Style

Pei, X., Tian, Y., Zhang, M., & Shi, R. (2021). The Boundary Proportion Differential Control Method of Micro-Deformable Manipulator with Compensator Based on Partial Differential Equation Dynamic Model. Micromachines, 12(7), 799. https://doi.org/10.3390/mi12070799

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