Fractional Modeling and Control of Lightweight 1 DOF Flexible Robots Robust to Sensor Disturbances and Payload Changes

: Model design and motion control are considered the cornerstones of the robotic ﬁeld that allow for achieving performance tasks. This article proposes a new dynamic modeling and control approach for very lightweight mechanical systems carrying payloads. The selection of the model and the design of the control are elaborated on using a fractional order framework under different conditions. The use of fractional order calculus is justiﬁed by the better performance that reveals a fractional order model compared to an integer order model of similar complexity. The mechanical structure of very lightweight manipulators has vibrations that impede the accurate positioning of their end effector. Moreover, they have actuators with high friction and sensors to measure the vibrations, which often are strain gauges, that have offset and high-frequency noise. All these mentioned problems might degrade the mechanical system’s performance. Hence, to overcome these inconveniences, two nested-loop controls are examined: an inner loop that controls the motor dynamics and removes the friction effects and an outer loop implemented to eliminate the beam vibrations by adapting the input-state feedback linearization technique. Then, we propose a new fractional order control scheme that (1) removes the strain gauge offset disturbances, (2) reduces the risk of the actuator’s saturation caused by the high-frequency noise of strain gauges and (3) reduces the dynamic effects of huge payload changes. We prove that our fractional controller has enhanced robustness with respect to the above-mentioned problems. Finally, the investigated approach is validated experimentally by applying it to a lightweight robot mounted on an air table.


Introduction
Dynamics modeling and control are keys to competence for industrial robot manufacturers.Current developments concerning them focus on improving the performance of robots, reducing costs related to robots, improving safety and introducing new functionalities as described in [1].Industrial robots are manufactured using heavy and rigid materials in order to impede the emergence of vibrations and static deflection in their mechanical structure and, hence, ensure accurate performance and avoid damage.However, rigid manipulators are also a source of problems as regards: 1.
More energy consumption due to the bulky structure.2.
Higher weight that often impedes their translation if mounted on a mobile platform.

3.
Higher design cost and operational risk.
Therefore, using robots whose links are very slender and are made of lighter materials is an excellent alternative to traditional robots [2].Unfortunately, these links show a flexibility phenomenon, which is a major cause of serious vibration problems that significantly decrease the precision of the positioning of the robot end effector and can engender damage to the mechanical system as well as to the environment if they are not properly modeled and controlled.Therefore, there is an important need to improve the models of these robots and their control methods in order to satisfy conflicting requirements such as increasing the performance of the robot by reducing its weight, which implies lower mechanical stiffness and larger oscillations, while trying to reduce the vibration modes (unwanted oscillations of the flexible link and the end effector) by implementing anti-vibration controllers.In order to develop suitable model-based control algorithms for flexible link robots, hereafter denoted FLR, it is necessary to construct computationally affordable dynamic models of these systems that incorporate such flexible behavior.In [3,4], literature surveys were reported on four commonly used dynamic modeling methods that depend on different assumptions: • The lumped parameter model: This method approximates the mass of the link through a finite set of lumped masses along its length that does not produce any kind of torque in the structure.The number of masses considered corresponds to the number of vibration modes in the model [5,6].
The assumed model: This represents flexibility by a truncated series of vibration modes.It consists of the linear combination of a set of space-dependent functions, called mode shapes, multiplied by time-dependent general coordinates [7,8].

•
The finite element method: The flexible link is modeled as a combination of a finite number of elements and the deflections are analyzed from the movement of small rigid bodies which leads to an important number of differential equations [9].

•
The transfer matrix method: This was originally used in optics and acoustics to analyze the propagation of electromagnetic or acoustic waves.The transfer matrix can represent each element of the flexible link system by transferring a state vector from one end of the element to the other.The whole system transfer matrix is obtained by multiplying the element transfer matrices together [10].
Dynamic models of FLR involve an infinite number of vibration modes.However, the methods mentioned above end up truncating models.Basically, the number of vibrational modes to be taken into account depends on the ratio between the links and the payload masses: the lower this ratio is (i.e, the lighter the links and the heavier the payload are), the smaller the number of significant vibration modes.
Links made of composite materials (e.g., the graphite-epoxy of fiberglass) are able to carry heavy payloads being much lighter than links made of metals such as aluminum.Then, robots with large and thin links made of these materials often exhibit very small linkpayload mass ratios and can be modeled for control purposes by only the first vibration mode of each link.Then, the lumped masses model of [5] can be applied assuming that all the moving mass is lumped at the end effector.Since this kind of link allows us to manufacture much lighter robots, this paper will be focused on their modeling and control.
The vibrations of the link are sensed using strain gauges.These sensors are cheap and simple to install.However, they present a relatively variable offset that depends on the temperature, and exhibit high-frequency noise [11].These are the major disturbances to be eliminated in this work.
In the context of modeling, Fractional-Order (FO) models have already been proposed for FLR and FJR in [12,13], where the authors proposed an algorithm to select the most suitable fractional model to describe the robot dynamics.They proposed an integer order transfer function which was then converted to an equivalent fractional order transfer function by tuning a fractional order λ ∈ (0, 1] in such a way that an error with the response of the system controlled by a PID was minimized.This approach is complicated and requires the assumption of the values of some parameters such as the damping ratio and the natural frequency of the mechanism and yields unnecessarily complicated fractional order models.In the context of FLR, fractionality is commonly used to model damping.Several of these models have been proposed in [14], but in many cases without experimental evidence.The model that has been more accepted and has been experimentally justified modifies the internal structural viscoelastic damping term of the model by changing the time derivative of this term by a fractional order time derivative between 0 and 2, e.g., ref. [15].
In the last few decades, several techniques have been reported for controlling FLR.Among the control strategies most used to dampen link vibrations, we mention Port-Hamiltonian modeling-based control [16], sliding mode control [17], adaptive control [18], optimal control [19], reinforcement learning control [20] and neural network control [8].
A paramount problem of FLR control systems is that they become untuned when the payload changes.Carrying different payloads from the nominal one-which is usually required in different tasks-produces changes in the vibration frequencies of the link, which yields degraded performance and even instability in the closed-loop system.This problem has been tackled using sliding controllers, e.g., [21].They yield accurate trajectory tracking even in the case of large payload variations (up to 200% of the nominal payload) but they need the fulfillment of the so-called matching condition; they provide a noisy control signal and their performance degrades if the measured signals that are feedback contain disturbances (like strain gauges).Robust H ∞ controller [22] has also been used to solve this problem.These controllers guarantee the stability of the the closed-loop system under payload variations or non-negligible spillover effects, but the quality of the reference tracking in the nominal case is not as good as that which is achieved with a controller specifically designed for the nominal payload, i.e., the stability margin is increased at the expense of losing system response performance (the system usually becomes slower).Moreover, these controllers have three more drawbacks: they allow only limited variations in the tip load, the quality of the trajectory tracking deteriorates considerably when the payload differs more from the nominal value and disturbances in the strain gauge sensor produce tracking errors and loss of precision in the end effector positioning.Another approach is using adaptive control.In [18], an adaptive controller based on an algebraic estimator of the payload and a generalized Proportional Integral control was designed to cancel the end effector vibration robustly with respect to payload variations.However, it did not take into consideration the presence of sensor offset disturbances nor the highfrequency noise, which degrades the payload estimate and the control performance.
The application of fractional calculus to control robotic systems has gained interest in the last few years because it yields controllers with enhanced robustness to external disturbances and parametric uncertainties.FO controllers have been applied to ensure control performances as tip position accuracy and suppression of the residual vibrations: FO adaptive back-stepping control [23], FO PID controller based on Kalman filter estimation [24], FO PID control [25], FO modeling and PI control of soft robotic arms [26], FO optimal control [27] and variable structure system fractional sliding mode controller [28], and a comparison was realized between the sliding mode controller and the fractional sliding mode controller in a non-commensurate order in [29].The latest review about fractional order modeling and control of flexible and rigid robotic manipulators is presented in [30].However, in the modeling part this paper focuses more on rigid than flexible structures, which are reported and illustrated with examples of rigid planar robots with different degrees of freedom and flexible joints.Finally, ref. [31] implements an FO controller combined with a standard PD controller [32] robust to strain gauge disturbances and payload variations.The design of this controller was based on the closed-loop pole allocation technique using an integer order model.This robust controller reduced the high-frequency noise introduced by the strain gauge sensor, eliminated the sensor offset effect and maintained the stability of the system during a large interval of payload variations.This article is used as the reference and the starting point of the present work.
The general contribution of the present work is introducing a new fractional order model of FLR and proving that the controller designed based on that model improves the experimental results attained by [31].Specific contributions are as follows: (1) proposing a new model of the dynamics of FLR that introduces a fractional order derivative in the inertia term of the Euler-Bernoulli equation instead of in the damping term, as is usually done; (2) proving experimentally the superior performance of this model; (3) modifying the control system [31]-which is based on a combination of the singular perturbation theory and the input-state linearization technique for an integer order model-to cope with the new FO model; and (4) proving experimentally that the resulting new control system outperforms the system presented in ref. [31].
This paper is organized as follows.Section 2 presents the experimental setup and its integer order dynamic model.Section 3 describes the new fractional order model and its experimental validation.Section 4 presents our control scheme.Sections 5 and 6 develop a control system robust to strain gauge disturbances and payload changes, respectively, based on the previous model.Section 7 describes the experimental platform and illustrates the results obtained experimentally with the new controller scheme compared to the one proposed in [31].Finally, Section 8 presents some conclusions.

Experimental Setup 2.1. Flexible Robot Setup Description
The contributions of this paper are experimentally validated in a setup located in the Robotics Laboratory at the Industrial Engineering School of the University of Castilla-La-Mancha.The setup consists of a lightweight single-link FLR mounted on an air table using a Harmonic Drive mini servo DC motor.The link is attached to the motor's output, while a disk rotates freely at the other end of the link on an air table which helps to cancel out the gravity and friction effects.The disk serves as the payload.Both the link and the payload move in a horizontal plane, supported by the three-legged metallic structure.The FLR is equipped with an encoder in the motor for precise measurement of the motor's angular position and a pair of strain gauges at the base of the link.Figure 1 shows this experimental platform and more details can be found in [31].The physical parameters of the motor and of the flexible mechanism are shown in Table 1 with numerical values.The actuator of the FLR is composed of three elements: a DC motor, a reduction gear and a current amplifier.The actuator generates the movement of the robot following the orders received from a control unit and, hence, it actuates the rotational joint of the robot.

Robot Dynamics
Under the assumption that the distributed mass of the link is neglected and the tip mass does not induce any rotational moment at the FLR tip because it can freely rotate around a vertical axis, a lumped-mass model of the robot with a single vibration mode is defined.Figure 2 shows a scheme of this FLR, in which θ m represents the joint angle, θ t is the angle of the tip, i.e., the angular position of the payload, m is the payload mass placed at the end of the arm and L is the length of the link.The stiffness of the flexible link is c =  The dynamic model of our FLR can be divided into two submodels: the first one represents the motor dynamics and the second one describes the behavior of the flexible link.These two submodels are linked by the coupling torque Γ (N•m) which is exerted between the motor and the flexible link and measured by the strain gauges placed at the base of the link.

Actuator Model
The actuator model connects the armature voltage to the motor angle at the gear output.It considers Γ c as the Coulomb friction torque, which is a piecewise constant perturbation dependent on the sign of the motor angular velocity (explained in references [31,33]).To address the static non-linearities of the DC motor, the Coulomb friction is compensated.Additionally, the coupling torque Γ is compensated to eliminate the impact of link dynamics on the motor.These compensations are achieved by defining a fictitious input u for the armature voltage, which incorporates the control signal, coupling torque and Coulomb friction.The actuator transfer function that relates the motor angle to the armature voltage fictitious input is defined by: where A = K/J and B = ν/J, K(N • m/v) is the electromechanical constant of the motor servo-amplifier system, J (kg•m 2 ) is the inertia of the motor and ν (N•m•s) is the viscous friction coefficient.More details about the actuator modeling are provided in [31].

Flexible Link Integer Order Model
The link is modeled as an Euler-Bernoulli beam with negligible distributed mass.Axial deformation is disregarded due to the presence of an air table that counteracts the effects of gravity and friction between the disk and the table surface.Since the torque at the base is a function of the second derivative of the tip position, second-order models are often used to describe the vibrations at the tip of this FLR.This model is derived based on the assumption of small deflections in the link, which allows us to approximate the dynamics of the FLR as linear.The input to this model is the motor angle θ m .The dynamics of the proposed lumped-parameter model are obtained by applying Newton's second law: where θt is the angular acceleration (rad.s−2 ) of the tip and the link deflection equation: Combining ( 2) and ( 3) and taking Laplace transforms on the resulting equation, the transfer function between the tip and motor angles is defined by: where is the frequency of the single vibration mode of our system, which is given by: Based on Equation (4), we obtain the transfer function between the coupling torque and the motor angle assuming zero initial conditions:

Fractional Order Model of the Flexible Link
Some systems have energy diffusion phenomena that can be better modeled using fractional order differentiation.Theoretical and experimental evidence shows that this phenomenon appears in FLR and is mainly associated with internal structural viscoelastic damping.As stated in the Introduction, this effect has been usually modeled by adding a damping term with a fractional order time derivative between 0 and 2 to the Euler-Bernoulli equation.

The Fractional Order Model
One main contribution of this paper is postulating a model for our FLR based on the fractional generalization of Newton's second law of motion.Then, we generalize the second-order ordinary differential equation of Newtonian mechanics by changing the order of the time derivative in the left-hand side of the Equation ( 2) to an arbitrary real number λ defined in the interval [1,2].This yields fractional differential equations such as the ones described in [14,34,35].This hypothesis is supported by recent studies of physics in which second-order classical wave equations, the Schrodinger equation, the Dirac equation and several others were generalized, changing some integer order derivatives with fractional order derivatives.Then, Equation (2) becomes which, combined with (3) gives Taking Laplace transforms on this equation assuming zero initial conditions, we obtain and, based on this equation, Note that Equations ( 9) and ( 10) become, respectively, ( 4) and ( 6) if λ = 2.

The Resonant Frequency
The resonant frequency ω r is the frequency at which G λ (j • ω) reaches its maximum value.Then, the frequency response of ( 9) is obtained taking into account that (j Let us represent by f (ω) the denominator of this equation.Then, the module squared of this function is: from which we obtain This expression shows that if 0 < λ ≤ 1 then the resonant frequency ω r does not exist.

Experimental Validation of the Dynamic Model
In this subsection, experiments justify the interest of the proposed model and the use of a constant fractional order λ that does not depend on payload changes.This is carried out by comparing the fitting to experimental data obtained from the motor encoder and the strain gauges of (1) the undamped integer order model (UD-IOM) given in ( 6), (2) the damped integer order model (D-IOM) given by (3) the integer order model with fractional order damping (D-FOM) given by and ( 4) the fractional order model (UD-FOM) given in (10).
As a first step, the parameters of these four models are identified from collected data of the motor angle θ m and the coupling torque Γ.The model is identified using the previous experimental measures, the Matlab software and the Mean Square Error (MSE) as the fitting function to be minimized.The identification method is presented in Algorithm 1: Algorithm 1 Model-based identification system algorithm 1. Collect experimental input-output data that represent the system's behavior to be modeled.2. Define transfer function structure.3. Create an initial model using the transfer function with arbitrary coefficients.4. Estimate the model parameters using optimization in Matlab to estimate the model parameters based on the data collected in step 1. 5. Evaluate the model obtained by computing the MSE which is the mean square of the difference between the predicted output of the system defined by ŷ and the experimental output defined by y.

if MSE > m ( ŷ is different from y) then
Repeat steps 4 and 5 until an acceptable level of accuracy is obtained (the value of the MSE is under the threshold m ).

end if
The trajectory of θ m (t) used to identify the dynamic models is a pulse signal 7 s wide with an amplitude of 1 rad.The pulse signal has a wide frequency content that excites the system at different frequencies and its impulse-like shape allows for a straightforward analysis of the system's response.Hence, for this trajectory, θ m (t), the MSE existing between the experimental torque Γ(t) and the coupling torque responses predicted via the four fitted models are 0.058 with UD-IOM, 0.00531 with D-IOM, 0.00524 with D-FOM and 0.00524 with UD-FOM.This means that introducing a fractional order derivative on the model of our FLR to account for the damping reduces the MSE of the standard second-order undamped model in 91.21%.Figure 3 presents the motor angle trajectory, the measured coupling torque and the coupling torque responses predicted via these four identified models.The proposed model is clearly able to introduce a damping behavior in the FLR.This model is more valuable and less complicated than an integer order model with a damping factor.Table 2 highlights that there is a significant difference in the MSE assessed from the response of the UD-IOM model and the rest: the error obtained from the first model is 10 times greater than the error of the rest of the models.The MSE has the same value for the UD-FOM and D-FOM.At this level, we introduce another criterion in order to evaluate the model quality.Accordingly, the AIC, which stands for Akaike Information Criterion, is needed to choose between the UD-FOM and D-FOM models as an additional selective indicator.It measures the relative quality of a model for a given set of data.It balances the goodness of fit of a model with the number of parameters used to define the model, penalizing models that are too complex and giving preference to models that fit the data well with fewer parameters.In the case of a huge number of observations (samples) n, its expression is where k is the number of parameters in the model and Lik is the likelihood function.Expression ( 17) can be written in terms of the MSE as: In general, a lower AIC score is better.Therefore, when comparing different models, the model with the lowest AIC score is usually considered the best among the candidate models.The estimated parameters obtained from the identification process of the four dynamic models presented, respectively, in ( 6), ( 10), ( 15) and ( 16) are summarized in Table 2, which also includes the performance indicators.Parameter c = 1.56 N•m was obtained via direct experimentation and we assume that it is constant in all present and future experiments.It can be concluded that the D-FOM and the UD-FOM have the lowest MSE but this last one has only two parameters to be tuned (λ, ), while the previous one has three (α, µ, ).The D-IOM has two parameters to be tuned (µ, ) like the UD-FOM, but it has an MSE slightly higher than the UD-FOM.Moreover, this last model has a lower AIC score compared to the D-FOM model, so it achieves a similar fit while using fewer parameters.All this is translated to the fact that the UD-FOM proposed in this paper has the lowest AIC value and, then, it is the most accurate model that, moreover, has the fewest parameters, i.e., is the simplest model.As a second step, we will prove that the fractional parameter λ of the identified UD-FOM model is invariant with respect to payload changes.Figures 4-6 compare the real and the predicted coupling torques for several values of the tip mass when integer and fractional order models are used.The system is excited by a Pseudo-Random Binary Sequence (PRBS) signal.This kind of signal ensures that the system is excited at a broad range of frequencies and amplitudes, allowing a more comprehensive analysis of the system's behavior over time.The real torque in Figure 4 presents a resonant frequency ω r,nom = 11.86 rad•s −1 for the nominal payload (m nom = 30 g).The real torque used in Figure 5 presents a resonant frequency ω r,min = 10.68 rad•s −1 for the maximal payload (m max = 40 g).The real torque in Figure 6 presents a resonant frequency ω r,max = 14.93 rad•s −1 for the minimal payload (m min = 20 g).The parameters of the damped systems (D-FOM and D-IOM) α and µ are considered constant for the three payload cases and are identified jointly for the three experiments in such a way that they give the minimum value of the sum of the MSE in the three payload cases.Then, α = 0.954 and µ = 1.37 were obtained.
The features of the three estimated torques are shown in Table 3 and with that table we can see that: 1.
The FO of the identified UD-FOM is kept constant during the payload changes, with a value λ = 1.92.

2.
The frequency changes significantly with the payload, in accordance with expression (5).

3.
The UD-FOM models have the lowest MSE compared to the rest of the models.
Finally, we remark that the identified natural frequencies are close but not equal to the resonant frequencies ω r obtained experimentally, because the damping coefficient or the fractional orders influence their relationship (see, e.g., ref. (14) for the case of UD-FOM).

General Control Scheme
The control system presented in this work incorporates a feedback mechanism using the angle of the actuator, denoted θ m , and the moment at the base of the link, denoted Γ.This control scheme, outlined and described in detail in the references [33,36], involves a nested structure consisting of two control loops based on the motor angle and the coupling torque: 1.
The inner loop employs the measured motor angle, θ m , to regulate the motor position.It effectively gets out the effects of non-linear Coulomb friction and time-varying viscous friction.Moreover, in order to ensure a faster dynamic response of the servocontrolled motor compared to the flexible link, this inner loop was implemented with a high-gain controller.

2.
The outer loop focuses on controlling the position of the payload at the tip to mitigate the occurrence of mechanical vibrations and ensure good trajectory tracking.
The control scheme presented here entails the use of the outer loop's control law output as a reference for the inner loop, splitting the control system design into two stages.The motor is then controlled in a manner that maximizes its response rapidity within the limits of saturation levels.To address the design of the outer loop, the singular perturbation theory is employed, which allows us to neglect the inner loop dynamics in the case that it is sufficiently fast.Subsequently, the outer loop controller is tuned to achieve the wanted performance related to the positioning of the tip.Both control loops are then merged into a singularly perturbed model and a stability analysis of the closed-loop plant is conducted.

The Inner Loop
The control scheme described in [36] is currently employed in this work.It is a robust controller based on two-degree-of-freedom PIDs.It incorporates a feedback term in order to compensate for the coupling torque-which makes the controlled motor dynamics less affected by the motion of the link-and mitigate the impact of Coulomb friction.This feedback term simplifies the motor model used in the controller design.
Hence, based on this simplified model, a two-degree-of-freedom PID control scheme is designed, as depicted in Figure 7.The control scheme consists of two PID controllers that incorporate a low-pass filter: According to the method outlined in [32], a pole-placement technique is implemented.The four poles of the closed-loop system are located at the same position, denoted z.Additionally, the two zeros are also positioned at z.As a result, two out of the four poles are eliminated by the two zeros, which gives the simplified closed-loop transfer function: By employing this control mechanism, precise trajectory tracking and rapid motor positioning are assured through effective compensation for disturbances such as unmodeled friction components and robustness against motor parameter uncertainties.To achieve very fast motor movements, a high value for |z| is selected, i.e., a small is looked for.

The Outer Loop
In order to remove vibrations in the flexible link, an outer control loop is designed by combining the input-state linearization technique and singular perturbation theory.The nonlinear controller proposed in [36] for a haptic device was adapted to suit our 1-DOF linear flexible link robot.By applying the principles of singular perturbation theory, the controller was initially designed assuming M(s) = 1 and, then, the stability of the closed-loop system including (21) was assessed.
The system to be controlled is characterized by the transfer function G(s) as given in Equation (9).By considering Figure 8 and employing the input-state linearization technique to G(s), the control law obtained is: In this context, θ t (s) corresponds to the intended trajectory of the manipulator, θ m (s) denotes the control signal produced by the outer loop that acts as the reference for the enclosed loop and C(s) represents the controller used in the outer loop.
In Figure 8, θ t is estimated using measurements of θ m and Γ.Specifically, θ t is calculated as mentioned in Equation (3), which gives θ t (t) = θ m (t) − Γ(t) c .If the measurements of the coupling torque were affected by a disturbance d (t), this would yield a modified estimated moment value Γ e (t) = Γ(t) + d (t).Consequently, the estimated value of θ t would be modified, becoming: Representing the disturbance in the estimation of the tip position as d(t) = − d (t) c , Figure 8 depicts the structure of our control system with the presence of this sensor disturbance.
After performing certain operations, the transfer matrix that establishes the relationship between θ t (t) and θ * m (t) with θ * t (t) and d(t) is determined: and substituting G λ (s) by its expression defined in (9) into this, The sensors used to measure the vibrations in the experimental platform are strain gauges [37].They are often used in robotic systems due to several advantages that make them highly recommended in applications.Their main advantages are recapitulated in the following: (1) they enable the measurement of vibrational signals as well as deflections, (2) they are easily implemented and used and can be reusable and (3) they are relatively low cost.However, strain gauges have several drawbacks [11,36], such as: (1) huge temperature variations, (2) high-frequency noise capable of saturating the actuators and (3) offset, which leads to a steady-state error that affects the end position of the link [38].These drawbacks will be addressed in the next section.

Control Robust to Strain Gauge Disturbances
Good trajectory tracking is a pledge of an efficient FLR controller.Therefore, some external and internal conditions should be taken into consideration when designing this control such as its ability to remove sensor disturbances such as the offset and the highfrequency noise [36] and its robustness to payload changes.

Strain Gauge Offset Disturbance Removal
Take into account the transfer function presented in Equation (25), which establishes the relationship between θ t (s) and d(s): As the offset can be represented as a step disturbance d(s) = d s , the steady state error e θ t induced in θ t can be assessed by applying the final value theorem: Considering that M(0) = 1, the aforementioned equation indicates that C(0) must be adjusted to a value of 1 to ensure that e θ t becomes zero.Therefore, the controller structure should take the following form: This implies that Ĉ(s) should have a zero at the origin to effectively eliminate the steady-state error caused by the offset of the strain gauges.

High-Frequency Noise Reduction
We represent the sensor noise by an unbiased high-frequency noise with an unknown statistical description, as was described in [39] from a relatively recent and non-standard analysis viewpoint.Then, this noise has a spectrum located at the highest frequency zone.
The control signal θ * m (t) is susceptible to amplified noise due to high-frequency noise in the sensor, which can result in actuator saturation, which ultimately leads to the rise of insufficient dynamic performance and potential instability.Let us examine the transfer function given in Equation (25), which establishes the relationship between θ * m (s) and d(s): By considering a controller in the form of Equation ( 28), we can deduce that: Assuming that the high-frequency noise in the strain gauge signal is represented by a sinusoidal signal with a significantly high frequency, the frequency response of P 2,2 as ω → ∞ demonstrates the impact of the amplification of this noise on the control signal θ * m of the outer loop.To further analyze this behavior, we divide the numerator and denominator of P 2,2 (s) by s λ resulting in: Assuming that lim it is easy to prove that expression (31) becomes lim ω→∞ |P 2,2 (j • ω)| = 0.Then, the impact of a high-frequency noise component of the signal θ * m (t) can be significantly diminished if condition (32) is verified.This condition involves that Ĉ(s) must be strictly proper.However, since high gains are used to close the inner loop, the influence of θ * m (t) highfrequency noise on the motor input is significantly amplified.As a result, the motor is prone to saturation.To avoid this undesirable effect, it is crucial to guarantee that condition (32) is satisfied.

Control Robust to Payload Changes
Variations in the payload mass m directly impact the natural vibration frequency , as expressed in Equation ( 9).These changes in frequency can potentially lead to underdamped or unstable closed-loop systems.To address this issue, this section focuses on the development of a controller C(s) to be incorporated in the control system shown in Figure 8 which enhances the robustness to payload changes.Additionally, the controller should exhibit robustness to the two disturbances studied in the previous section.Then, conditions ( 28) and ( 32) must be verified.
Based on the FO model defined in Equation ( 9), the first path of the Nyquist diagram (indicated in Figure 9 in brown color) corresponds to the frequency response of G λ (s), where ω is defined in the interval 0 < ω < ∞.This plot shows that the system exhibits a small positive phase margin.Therefore, though the closed-loop system is stable, its phase margin has to be increased in order to improve the damping.It can be achieved by designing a phase lead controller Ĉ(s).Let us assume that such a controller is an FO derivative controller: The FO operator s β is often approximated numerically.We use in this paper the Grünwald-Letnikov definition based on the backward difference.
The evident simplest choice for (33) is the classical integer order derivative controller, which is obtained by making β = 1 in this expression.It verifies that Ĉ(0) = 0 and can improve the damping of the closed-loop system ζ by increasing its phase margin φ m in π 2 rads.However, we cannot achieve desired closed-loop settling time t s because we can only tune one parameter, k d , in this structure.Then, we tune the two parameters β > 0 and k d of controller (33) by using a frequency domain technique and the following approximate relations: φ m ≈ 100ζ-which is an expression originally developed for second order systems (e.g., ref. [40])-and t s ≈ 2π/ω c .
The blue line of Figure 9 represents the first path of the Nyquist diagram of the open-loop transfer function of the controlled system: in which the dependence on has been made explicit in L() and G λ ().The gain crossover frequency ω c and the phase margin φ m are highlighted in this plot.According to this diagram, the control system is designed such that ω c > ω r .Next, we develop some results useful for designing a control system robust to payload changes.
Lemma 1.Consider the frequency response L(j • ω, ) given by (34).Its phase is a decreasing function with ω that decreases from π 2 • β rads to π 2 • (β − λ) rads.Moreover, changes in the value of only produce a homothety of the Nyquist plot of L(j • ω, ) with the center in the origin.
Moreover, expression (34) can be written for which illustrates that the shape of L(j • ω, ) is the same for any and, in this expression, variations of only produce changes in the scaling Proof.First of all, note that condition (39) is the same as and, taking into account that λ > 1, this inequality implies that Next, assume the open-loop frequency response (34).The relationship among its phase margin φ m , its gain crossover frequency ω c and its natural frequency , is given by, e.g., ref. [41]: Equating magnitudes in this expression and squaring yields and equating phases Differentiating expression (44) with respect to (note that and ω are linked by (43)), we obtain that Differentiating expression (43) with respect to and ω c , we obtain that Substituting (43) in the first and second terms of the left side member of (46) and operating yields and then Substituting ( 48) in ( 45) and operating gives In this expression: 1.

2.
Since 1 < λ ≤ 2, the numerator is less than or equal to zero.
Then, the right side member of (49) is less than or equal to zero and the theorem is proven.
Remark 1.In expression (48), both the numerator and denominator are positive, because ( 41) is verified and the denominator is the same as the one of (49).Then, under the conditions of this theorem, ω c is increasing with .
Next, we propose a theorem that facilitates the verification of condition (39).

Theorem 2.
Consider the open-loop transfer function (34) whose gain crossover frequency is represented by ω c .Assume that 1 < λ ≤ 2 and 0 < β < λ.Moreover, assume a value of the natural frequency of the link − that is greater than one and verifies condition (39).Then, function is less than zero in the interval − ≤ < ∞.
Proof.First, note that ω c is a function of , according to (43).For this reason, h has a single argument .Differentiating expression (50), we obtain that Substituting (48) in (51) and operating yields Since (39) is verified, inequalities (40) and ( 41) are also verified.Then, the left addend of the numerator is negative.

2.
Since inequality (37) is verified, the right addend of the numerator is also negative.

3.
The denominator is equal to the denominator of (49).Theorem 1 proved that it is positive.
The consequence of these three items is that dh( ) d < 0 at − .Subsequently, we will prove by contradiction that h( ) is negative in the interval − ≤ < ∞.Assume a value + at which h( + ) > 0. Expression (38) shows that ω c varies continuously with .Then, h( ) is a continuous function and, consequently, if h( − ) < 0 and h( + ) > 0, a value 0 exists such that h( 0 ) = 0 and h( ) < 0 in the interval ∈ [ − , 0 ).The following relation can be established: and, since it has been proved above that dh d < 0 if h( ) < 0, the previous integral in the interval [ − , 0 ] is negative.Consequently, h( 0 ) < h( − ) < 0, which is in contradiction with the previous condition that h( 0 ) = 0.
Assume that 1 < λ ≤ 2 and the payload m varies in the range [m, m].Then, the natural frequency varies in the range [ , ] obtained from the previous interval through (5).In the following, we propose a procedure to design controllers robust to payload changes, i.e., to changes, that verify the specifications of a minimum phase margin φ m (at m) and a minimum gain crossover frequency ω c (at m) in all the range [m, m]: 1.
Choose the gain crossover frequency ω c corresponding to .It must verify condition (39) and, according to Remark 1, it is the minimum ω c achieved in the interval [ , ].

2.
Apply Equation (43) to the pair , ω c .This gives the relation Then, to calculate k d , we need to obtain the value of β, which will be done in the next step.

3.
In accordance with Theorems 1 and 2, the minimum phase margin in the range [ , ] is located at .This value is represented by φ m .Then, applying condition (42) to and substituting k d by (54) yields This complex equation can be transformed into two real equations that allow us to obtain the two values ω c and β.

4.
Obtain k d from the left side equation of (54) using the already calculated χ and the fractional order β calculated in the previous step.

5.
Then, the designed controller is which can be modified as where n = 1 if β < 1 and n = 2 if 1 ≤ β < 2 in order to fulfill the high-frequency noise condition (32).In this expression, δ is chosen small enough to not modify the frequency characteristic of the open-loop transfer function L(j • ω, ).Since this modification produces a phase lag in L(j • ω, ), a criterion to design δ is to choose a value such that the minimum phase margin φ m is reduced by up to 10% of its initially designed value (or, alternatively, increase φ m in a quantity that accounts for the phase margin reduction subsequently caused by introducing (1 This algorithm is not iterative: once φ m and ω c have been defined, the algorithm gives or does not give the solution in a single realization.If the algorithm fails to give a solution, specifications have to be changed.Typically, the problem can be solved by reducing ω c .

Application to a Single Lumped Mass Flexible Link Manipulator
In this Section, FO controllers are designed for the robotic setup described in Section 2. First, a controller C M (s) is designed to apply the procedure described in Sections 5 and 6 to UD-FOM (9), (10).Then, a controller C M (s) is designed using UD-IOM (4), ( 6) and the procedure described in Section 5, in such a way that it has the same specifications as C M (s) for the nominal payload, in order to make a fair comparison between controllers and assess the effect of using the fractional order model instead of the integer one in the controller design.

Design of the Controllers
The application of the procedure of Section 6 to design C M (s) is as follows: 1.
The gain crossover frequency corresponding to the minimum natural frequency = 9.42 rad•s −1 is chosen as ω c = 16 rad•s −1 .Note that the condition > 1 of Theorem 2 is verified.Then, it is the minimum ω c achieved in the interval [ , ].

4.
Obtain k d from the left side equation of (54) and the fractional order β calculated in the previous step.We obtain k d = 0.148.

5.
The proper controller is designed by choosing δ = 0.16 and n = 2 instead of n = 1 (as stated in the procedure) in order to enhance the filtering of high-frequency noise.
Then, the designed controller is The gain crossover frequency and the phase margin of the UD-FOM (9), (10) with controller (58) at the nominal payload m nom = 30 g are ω c,nom = 19.73rad•s −1 and φ nom = 70.087• , respectively.
Subsequently, a control system with the same structure as the previous controller is designed, but assuming the standard undamped integer order model UD-IOM, i.e, all the parameters of the control system are tuned to this model making λ = 2 in all the expressions of Section 5.In order to make a fair comparison with the previous design, this controller is tuned to have the same frequency specifications ω c,nom and φ nom as the previous controller at the nominal payload.
Imposing the complex condition (42) and making λ = 2, φ = φ nom and ω c = ω c,nom in this equation, we obtain k d = 0.241 and β = 0.65.Then, the controller designed using the integer order model is Figure 10 illustrates the variations of the gain crossover frequency ω c and the margin phase φ m of the open-loop transfer function in function of the payload, when using the two controllers C M and C M and including or not including the inner loop motor dynamics, i.e., in the cases: In these expressions, M(s) has its poles in z = −80 (see expression (21)).We notice in this figure that 1.
The plots of ω c and φ m intersect at the nominal payload in the cases of L M (j ω c is decreasing with m (increasing with ) as stated in Remark 1.

3.
φ m is increasing with m (decreasing with ) as stated in Theorem 1.

Experimental Validation of the Controller
Figures 11-13 show the experimental trajectory tracking of the robot achieved using controllers C M (s) (58) and C M (s) (59), in the cases of the nominal (m = 30 g) payload, and the maximum (m = 40 g) and the minimum (m = 20 g) payloads that have been used in the design of the robust controller.These experiments were carried out using strain gauges that have a small offset, as Figure 14 shows in the steady state of an example of trajectory tracking.These results show that (1) both C M (s) and C M (s) controllers remove the steady-state error in the tip position caused by the strain gauge offset; (2) the responses using C M (s) are slightly more oscillatory than when using C M (s), especially in the case of the minimum payload; and (3) the robot's responses when tracking a trajectory in the cases of the nominal and the maximum payload are approximately the same using both controllers.To sum up, both controllers provide good trajectory tracking in the cases of the nominal and minimum resonant frequencies with steady-state error elimination and zero overshoot.In the case of the maximum resonant frequency-associated with the minimum mass m = 20 g-the tip oscillations are completely damped using the C M (s) controller, while small oscillations remain if C M (s) is used.Finally, we show that the closed-loop behavior is appropriate in the most challenging case of carrying very low payloads that are outside the payload interval used for the design of controller (58).Figures 15 and 16 show the responses of both controllers for m = 16 g and m = 8 g.Large oscillations can be seen in both figures in the responses obtained using C M (s), while controller C M (s) still damps the oscillations in the case of m = 16 g and yields an oscillation of much lower amplitude than the one of C M (s) in the case of m = 8 g.

Discussion
At first sight, one might think that using a fractional model of order 1.92 instead of the well-known integer model of order 2 in the design of the control system of our FLR would be irrelevant since these are very close orders.However, our analysis and experimental results show that it makes a significant difference in the closed-loop dynamics when low payloads have to be translated.Next, we analyze the reason for these differences.
We observe in Figure 10 that: 1.At payloads lower than the nominal payload, the gain crossover frequency of L M is higher than that of L M and the gain crossover frequency of LM M is higher than that of L M .2.
At payloads lower than the nominal payload, the gain crossover frequency of LM M is lower than that of L M and the gain crossover frequency of LM M is lower than that of L M .This can be easily understood because M(s) is a low-pass filter that reduces the amplitudes of L M and L M , causing a reduction in the gain crossover frequencies.

3.
The effect of adding M(s) to the open loop is a strong reduction in the phase margin at all the payloads, as a comparison of the middle and lower subplots shows.4.
At payloads lower than the nominal payload, the phase margin of L M is higher than that of L M . 5.
However, the lower subplot of Figure 10 shows that, at low payloads, the phase margin of LM M is higher than that of LM M .This is caused by the fact that, according to the higher subplot of Figure 10, the gain crossover frequency of LM M is lower than that of LM M .Then, M(s) produces less phase lag on L M than on L M .This compensates for the positive difference in the phase margin between L M and L M , yielding higher phase margins of LM M at low payloads.
This last issue has been corroborated via the experiments carried out with small and very small payloads.

Conclusions
The robust control of very lightweight flexible link robots has been addressed in this paper.In these robots, it is often assumed that all the mass is lumped in the tip payload.Then, their dynamics have a single vibration mode with little damping.A control strategy for this robot needs to guarantee robustness to (1) the nonlinear motor frictions, (2) the offset and the high-frequency noises of the strain gauge sensor implemented in the robot platform and (3) the variations of the carried payload.We achieved these objectives in previous work [31].
Usually, robot damping is modeled by adding to the Euler-Bernoulli beam equation a term proportional to the angular velocity of the robot (see (15)).However, we have proven experimentally that, for our specific robot, a better fitting of the damping phenomenon is achieved modifying Newton's second law of motion by using an FO derivative of the angle instead of the standard second order derivative (see (10)).Then, the input-state feedback linearization method based on the singular perturbation technique and an FO PD controller-which was used in [31]-was modified to account for such a fractional order model.
The contributions of this paper are therefore as follows: (1) experimental justification of the fact that, in some cases, the damping phenomena of the dynamics of flexible link robots are better fitted using models based on the generalization of Newton's second law of motion by using FO differentiation; (2) a new fractional input-state feedback linearization scheme; (3) a lemma and two theorems that allowed us to define a procedure to tune the FO controller embedded in this control scheme; and (4) proving analytically and experimentally the advantage of this new control system based on the FOM over the previous control based on the input-state feedback linearization method that, though yielding an FO PD controller, its design was based on the standard integer order model.These advantages are noticeable at low and very low payloads, cases in which the natural vibration frequency of the robot becomes very high and the closed-loop control is prone to instability.
Our future work will be focused on extending this FOM and controller to multilink robots.

Figure 2 .
Figure 2. Schematic Representation of a Single Link Flexible manipulator.

Figure 3 .
Figure 3.Comparison of the Different Predicted Torque Models.

Figure 4 .Figure 5 .Figure 6 .
Figure 4. Comparison of the Different Predicted Torques for the Nominal Payload.

Table 1 .
Values and units of the robot properties.
3•E•I L , where E • I (N•m 2 ) is the flexural rigidity of the link.

Table 2 .
Parameters of the Predicted Models.

Table 3 .
Features of the Predicted Models in the Different Payload Changes.
(33)a given controller Ĉ(j • ω) of the form(33), the phase margin of the open-loop transfer function (34) is non increasing with , i.e.,