#### Abstract

A study is presented in this paper that uses a simple mechanical analogy to analytically tune the PD (proportional-derivative) controller of a linear flexible manipulator system. More specifically, the aim is to give simple closed-form solutions of the optimal P and D gains to yield the maximum bandwidth under a given damping requirement *or* conversely the maximum damping under a given bandwidth requirement. The idea of this study is based on the observation that the performance of the complete manipulator system is largely determined by the operational dynamics of the fundamental vibration mode. A lumped element method is thus applied to model this dynamics in terms of simple lumped mechanical elements. It subsequently turns out that the original servo control problem is analogous to a conventional Zener mount design problem, that is, mathematically, to optimize a third-order dynamic system consisting of the Zener model of a viscoelastic mount and an inertial object upon it. A design methodology is finally established to analytically determine the optimal elements of the mount, corresponding to the optimal control gains. Simulations and experiments were also conducted with a single-link flexible beam to support the model and the design methodology developed.

#### 1. Introduction

There are various methods to tackle vibration problems of a flexible manipulator system [1], which is a servo machine having a flexible joint or link. They include tuning methods of PD (proportional-derivative) gains [2] and PID (proportional-integral-derivative) gains [3]; feedforward control methods using various filters including input shaping [4], notch [5], and inverse [6] filters; and feedback control methods using additional vibration sensors such as an accelerometer [7]. Amongst, the PD-tuning method is the simplest and is of interest throughout this paper. It is particularly suitable for linear [8] and linearized [9, 10] manipulator systems, where no nonlinear dynamic effects (e.g., gravity and friction effects) are present. An understanding of this is important for control engineers as it is also a prerequisite to an understanding of the PID-tuning method that is then for general nonlinear manipulator systems.

Irrespective of the methods listed above, the fundamental vibration mode of the complete servo machine plays a key role in determining the system performance [1, 6]. The PD-tuning method is basically to tune the vibration frequency and the damping property of this mode, which are closely related to the bandwidth (i.e., speed) and damping (i.e., settling time) requirements of the system, respectively. The PD-tuning method is thus to give an optimal servo machine in terms of both speed and settling time. If an excessive speed had to be practiced due to some other requirements, this mode could then be suppressed by an additional means such as a control method. The control methods listed earlier are thus generally for rapid servo machines, which are nonoptimal in the sense of gain-tuning. The study presented in this paper is specifically concerned with the PD-tuning method for linear servo machines that are further nonrapid. Unlike the situation for rigid manipulators [2, 11], analytical studies towards the optimal PD gains for flexible manipulators are scarce in the literature. Instead, numerical and empirical tuning methods have been exclusively practiced by inspecting, for example, step responses. However, these are trial-and-error practices whose repetitions would not necessarily help engineers to build up an understanding of the dynamics and control mechanisms involved.

An analytical study is thus conducted in this paper that uses a simple mechanical analogy to determine the optimal gains. More specifically, the aim of this study is to give simple closed-form solutions of the optimal P and D gains to yield the maximum bandwidth under a given damping requirement *or* conversely the maximum damping under a given bandwidth requirement. To model the operational dynamics of the fundamental vibration mode, the lumped element method (LEM) [12–14] is applied to transform the mechanical manipulator into a system consisting of a base mass, an appendage mass, and a resilient member connecting the two masses [15]. The P and D gains are then, respectively, transformed to a spring and a damper underneath the base mass, thus yielding a simple positive-definite two-body lumped element system (LES) representing the complete flexible manipulator system. Because the dynamic effect of the base mass can be ignored by implementing a large P gain, which is an interesting finding, the original servo control problem is eventually transformed into a conventional viscoelastic mount design problem dealing with a one-body LES. More precisely it is, in fact, a Zener mount design problem [16, 17], where a given inertial object (the manipulator inertia) is supported by a Zener mount that is comprised of a given spring (the manipulator flexibility) in series with a parallel connection of another spring (P gain) and a damper (D gain). A design methodology is finally established by adapting the conventional design rule of the Zener mount system [18, 19], to analytically determine the optimal P and D gains.

The paper is organized as follows. Section 2 describes the modeling process to simplify the servo control problem to a Zener mount design problem. Section 3 describes the design (tuning) process to determine the optimal elements of the Zener mount, corresponding to the optimal P and D gains. New indexes are also proposed to quantitatively specify the system bandwidth and damping requirements. Simulation and experimental work are described in Section 4, followed by some conclusions in Section 5. There are also two appendices in this paper to summarize some important frequencies used in the modeling and design processes and also to review and compare two other alternative Zener mount design methods based, respectively, on convex minimization and frequency tuning.

#### 2. Modeling

##### 2.1. Lumped Element Modeling

Consider a servo positioning system driving a flexible manipulator in the horizontal plane as shown in Figure 1(a), where the beam represents the flexible manipulator. The motor rotor is PD-controlled using the P gain and the D gain to locate the tip at the position set by the command input . The P and D gains are generally related to the speed and the settling time of the servo machine, respectively. The Laplace operator is denoted by . A command force input is additionally shown for completeness. It is assumed throughout this paper that the dynamics of the mechanical manipulator are linear (*Condition* I: linearity).

**(a)**

**(b)**

Applying the novel LEM [12–14], the low-frequency operational dynamics of the mechanical manipulator can be closely approximated, as shown in Figure 1(b), by using a semidefinite system consisting of the base mass , the appendage mass , and a resilient member (a parallel connection of a spring and a damper ) connecting the two masses. A series connection may alternatively be used. The oscillator (, , and ) represents the first vibration mode of the beam in the fixed-free boundary condition. The physical variables are used in a generalized sense so that the force () can mean a force or a torque; and the motions ( and ) can be translational or rotational. The frequency dependence of these variables has been suppressed for simplicity, where , is the angular frequency, and is the time. Similarly, physical quantities are used in a generalized sense so that a mass can mean a moment of inertia; and the translational springs and dampers indicated can mean rotational springs and dampers, respectively.

Since the PD controller considered is a passivity-based controller [13], the control gains and can be, respectively, transformed into a spring and a damper underneath the base mass so as to complete the model of the positive-definite two-body LES representing the complete manipulator system including the servo actuation part. A hypothesis implicitly used here is that the servo actuator is nonrapid (not a very large ) such that it cannot excite any higher order modes of the mechanical manipulator but only the first one (*Condition* II: nonrapid operation).

Applying now the impedance approach [20–22], it is convenient to regard the two-body LES as a coupled system between the base and the appendage as indicated in Figure 1(b). The base impedance is given bywhere if uncontrolled. Thus, is variant depending on the control gains used. In contrast, the appendage impedance is fixed and can be written aswhere the vibration transmissibility is given bywhere is the angular natural frequency and is the damping ratio of the uncoupled appendage. Similarly, those of the uncoupled base can be written as and , respectively. The transmissibility defined above is a useful term, which is also widely used in the field of vibration isolation, to assess the vibration transmission characteristic between two points of a vibration system [23]. The impedance of the combined system seen by the source can then be written as

The combined system without control can be identified by curve-fitting this impedance expression to the measured impedance response [12, 13], as experimentally demonstrated in Section 4. The combined system with control can then be easily predicted by specifying the gains and used. As such, a complex-shaped manipulator system even with multiple links can be similarly identified. In fact, the LEM, more fundamentally the impedance approach, is a general approach that can be applied to any linear, passive operational dynamic systems that the conventional Newtonian method (i.e., the mobility approach) is also applicable [14].

If the system is now operated by the motion source only, the vibration transmissibility from to (i.e., the hub position) can be written asand that to (i.e., the tip position) is given bywhere , , and . Overall, the modeling process can be summarized as follows: given the (measured or simulated) driving point frequency response data of at the hub of a manipulator, the LEM is applied to identify the elements (, , , and ) of the LES; and then the complete LES further incorporating the control gains ( and ) is used to predict the hub and tip responses.

##### 2.2. Coupling Analysis

The dynamic behavior of the two-body LES is greatly dependent on the strength of coupling between the two subsystems. As an extreme case, the criterion for weak coupling can be written as [21]. For the specific base system of the variant spring with the fixed mass , this criterion can be more conveniently rewritten asThe former is the case when the servo actuator has a very small operation bandwidth (a very small ) far below the appendage natural frequency . The latter is the case when it has a very large operation bandwidth (a very large ) far above . In both cases, the base damper is little helpful for dissipating the appendage vibration. If the coupling is now made to be strong by appropriately adjusting (consequently ), which means that Equation (7) no longer holds, it is then possible to dissipate the appendage vibration by using . This is because the mode being formed is now a global mode. This is the characteristic that is exploited throughout this paper. It is also fortunate that the range of suggested here (i.e., neither very small nor very large) does not disobey Condition II.

To support the discussion above, Figures 2(a), 2(c), and 2(e) and Figures 2(b), 2(d), and 2(f) show the magnitudes of the transmissibility functions and for the system in Figure 1(b), respectively. Three values of the base damping ratio were tested for three values of the coupling term . The *resonance frequency * (i.e., the first peak frequency) and the damping property of the fundamental vibration mode of are closely related to the bandwidth and the damping requirement of the complete operational manipulator system, respectively. It can be seen in Figures 2(a) and 2(b) that if the strength of coupling is weak (e.g., ), the base damper has very little effect on the fundamental mode. It can also be seen in Figures 2(c) and 2(d) that, for a stronger coupling (e.g., ), the system damping improves but the bandwidth diminishes. The maximum coupling case () shown in Figures 2(e) and 2(f) is a special case, where the base system acts as “a *suspending* dynamic vibration absorber” to maximally damp the tip vibration of the appendage [24]. It can though be seen that the bandwidth greatly diminishes. Thus, there is a trade-off between the system damping and the system bandwidth; the damping can only be increased by compromising the bandwidth, and vice versa. These observations are theoretically verified later in this paper.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

##### 2.3. Model Reduction by Eigenvalue Analysis

From the discussion in the previous subsection, it is clear that the proximity between the two natural frequencies, and , can mean the strength of coupling. There are two values of to give the same proximity: a large is desirable such that (Figure 2); a small would result in a small bandwidth as it yields (not shown in Figure 2). This suggests an important design guideline as follows (*Condition* III: nonslow operation):Equivalently, it can be rewritten as or in which and . Note that this guideline sets the lower bound of the P gain. Its upper bound is then limited by Condition II as well as the strength of coupling.

Under the condition given by Equation (8), it is *remarkable* that the two-body LES in Figure 1(b) can be further reduced to the one-body LES shown in Figure 3(a) for the fundamental vibration mode and then to that in Figure 3(b) for the second vibration mode. The appendage damping has been ignored for simplicity. To justify these model reductions, Figure 4 compares the natural frequencies of the two one-body LESs with the exact ones obtained by solving the eigenvalue problem of the original two-body LES for a spring ratio of (a constant), according to the mass ratio (a variable). The region is shaded. It can be seen that the exact (solid lines) and approximate (dashed) eigenvalues agree well within the shaded region; the agreement gets better as increases (i.e., as increases). Thus, it can be stated that the two one-body LESs shown in Figures 3(a) and 3(b) are reasonable models if Equation (8) holds. A similar discussion can be made when the mass ratio is fixed: (a constant). The natural frequency of the first LES in Figure 3(a) is henceforth called the *nominal natural frequency* of the fundamental vibration mode and is denoted by , where is the static stiffness.

**(a)**

**(b)**

Overall, it has been shown that the design guideline given by Condition III has naturally acted as the model reduction condition to the one-body LES. The first one-body LES shown in Figure 3(a) is a simple mechanical analogy to the operational dynamics of the fundamental vibration mode of a complete flexible manipulator system. Here, the first LES model together with Conditions I–III is an important contribution of this paper. A merit of this is that it is now simple enough to offer a rigorous definition and a subsequent solution of the relevant optimal design problem of the system, as presented in the following section.

#### 3. Optimal Design

##### 3.1. Background

Of interest in this section is the simple one-body LES shown in Figure 3(a), where and mean the inertia and the flexibility of a given mechanical manipulator, respectively. It is *striking* to note that the mount () is, in fact, the series-type Zener model often used for describing the dynamic behavior of a viscoelastic material [16], where − and , respectively, indicate a series and a parallel connection of impedances. Such a mount that can be modeled by the Zener model may be called a *Zener mount*, whereas the complete vibration system further containing the inertial object is called in this paper a *Zener mount system*. The control problem here is also to somehow optimally tune the two elements ( and ) in consideration of the residual transient vibration of . Thus, it can be seen that the original servo control problem is analogous to a conventional Zener mount design problem. An advantage of this analogy is that the linear vibration theory established in the field of vibration isolation can also be applied to study the dynamics of the servo control system in the field of motion control. More specifically, the conventional design rule of a Zener mount system [18, 19] can now be applied to design the optimal PD controller. This task is undertaken in this section for the particular form of the Zener mount system shown in Figure 3(a), where and are fixed and known while and are tunable and unknown. As can be seen in what follows, it is a nontrivial but challenging task because the rule appears yet rather a loosely-defined concept on how to solve than any closed-form solution ready to use.

Prior to optimal design, we thus first review the dynamic behavior of the Zener mount in Figure 3(a) whose impedance can be conveniently written aswhere , , , and again . A parallel-type Zener model can also be similarly represented. The complex modulus of the material (i.e., mount) is , whose real and imaginary parts indicate the dynamic (frequency-variant) stiffness and loss functions, respectively. The stiffness function is related to the static and dynamic load-carrying capabilities while the loss function is to the damping capability. It is interesting that the material behavior is completely determined by the lead compensator whose characteristics are well known in the control society [2].

As an illustration, Figure 5 displays the polar, phase, and amplitude responses of a lead compensator for with . The loss angle and factor functions of the material can be, respectively, defined by and , where the latter is also . Since the loss factor is more fundamental and general than the damping ratio, it has been widely used in evaluating a viscoelastic material [25]. The transition frequency is critically important in mount design since it is the frequency offering the maximum use of the dynamic loss capability with respect to the dynamic stiffness capability. The maximum loss angle and factor at this frequency are then, respectively, given by and [2]. It is noteworthy that either one of and is uniquely determined by the other, which fact will later be used to give a practical definition of the optimization problem. There is also a frequency giving the maximum loss itself, denoted by in Figure 5(a). A mount is seldom tuned to this frequency as it induces a large increase in the dynamic stiffness but a small increase in the dynamic loss, compared with those at .

**(a)**

**(b)**

##### 3.2. Problem Definition and Solution

We then consider the optimal design of the Zener mount system that is a combined system between the inertial object and the Zener mount , as shown in Figure 3(a). From the dynamic equation with being the velocity of the inertial object, the characteristic equation can be written in the frequency domain asand more explicitly in terms of asIn general, this cubic equation has one real root and two complex conjugate roots. By adapting the *conventional* design rule for random vibration control: “” [18, 19], which is to tune the transition frequency of the mount to the resonance frequency of the fundamental vibration mode, we may now define the *general* design rule for residual vibration control considered in this paper as follows:where is the *fundamental damped natural frequency* of the complete system in free oscillation that can be given by the imaginary part of the complex conjugate roots. Mathematically, it can be seen that this is to find the solution of the third-order dynamic system in Equation (10b) satisfying the rule given by Equation (11). Here, one of the two unknowns, and , should be further (either explicitly or implicitly) given to yield a unique solution. As the rule literally suggests, this is also an optimization problem to maximize the loss factor of the system in free oscillation, to reach , at (*Problem definition*).

Exact analytical expressions of are though very lengthy to handle [17]. Use of some approximates would be more intuitive and practical. Two simple approximates considered in this paper are the nominal natural frequency defined earlier and the *standard natural frequency* defined now as in which is the dynamic stiffness at (*Approximations* I and II: fundamental vibration frequency). Using first for in Equation (11) such that , we haveUsing then such that , which is especially called the *standard* design rule in this paper, we haveEquations (12) and (13) are approximate solutions for whenever is explicitly given. Conversely, it is also possible to obtain from a given . Accuracy of each of the two approximate solutions will be assessed later in this section.

However, a problem that still remains but is little addressed in the literature is that design tasks are rarely specified in terms of and but in terms of the bandwidth and damping requirements. We thus require some indexes to quantitatively specify these. Noting that and are readily available unlike and , we firstly propose that the system bandwidth requirement be specified by one of these (*Convention* I: bandwidth index). Consider then a truly optimal Zener mount system perfectly satisfying Equation (11). As it is now a vibration system containing , the material property expression given earlier can be rewritten asAs it is also a vibration system fundamentally vibrating at in free oscillation because of Equation (11), defined at can now mean the system damping property. This is called here the *effective loss factor*, denoted by for distinction. We thus secondly propose that the system damping requirement be specified by (*Convention* II: damping index). These two conventions are valid throughout this paper.

Rewriting Equation (14) can give the generic trade-off relationship between the bandwidth and the damping requirement:which also theoretically explains the trend of the simulation results in Figure 2. Referring also to Figure 5(a), it can be seen that each of these gives the maximum value (equivalently ) to achieve a given level of (determined by ) or conversely the maximum value under a given level of . Therefore, it can be finally stated that the general design rule given by Equation (11) is *practically* to yield the maximum bandwidth (the optimization goal) under a given damping requirement (the constraint) *or* conversely the maximum damping under a given bandwidth requirement of the system in free oscillation (*Problem redefinition*). It is interesting to note that the bandwidth and damping requirements of the mount are analogous to the handling and ride comfort requirements of a vehicle [26], respectively. In contrast, a suspending dynamic vibration absorber can be more easily designed without any consideration of the bandwidth (or handling) requirement [24].

Rewriting Equation (14) further gives the optimal , whenever is specified, such thatIt should be noted that Equations (15)–(17) are generic relations valid for all truly optimal Zener mount systems perfectly satisfying Equation (11). A specific design case for a given manipulator of mass ratio then imposes some limitations to the values of Equations (15)–(17) such thatall because of the condition given by Equation (8).

Figure 6 displays a curve (thick line) representing Equation (16) and an area (shaded area) confined by Equation (18) for a manipulator of . Equation (15) may alternatively be used to draw the curve. The curve is called here the *optimal design curve*, while the area is the *design region*. The curve is invariant while the area enlarges as increases. The bandwidth and damping requirements can then be specified at a point along the optimal design curve within the design region. It subsequently gives the exact solution for as explicitly given by Equation (17). The corresponding D gain can then be approximately determined by Equations (12) or (13). It is desirable that the manipulator is of a sufficiently large mass ratio so that its design region includes a critical value of (e.g., ) giving the critical damping. This is often satisfied in many practical flexible manipulators (e.g., ).

##### 3.3. Simulations

Simulations were performed to compare Equations (12) and (13) by using the two-body LES in Figure 1(b). The data used and the three design cases () considered are given in Tables 1 and 2, respectively. The three effective loss factors specified are all within the valid design region since as indicated in Table 2. The optimal design curve, equivalently Equations (15)–(17), yielded and for each case. The gain was then obtained by Equations (12) and (13).

Figure 7(a) shows the resulting tip responses in frequency, calculated from Equations (5) and (6). It can be seen that (denoted by “ο”) estimates (denoted by “∇”) better than (denoted by “+”) does, where each indicated was numerically obtained from the roots of Equation (10) (“roots” in Matlab®). The deviation gets larger as increases. It is clear from the simulations that Equation (13) is more accurate than Equation (12). Consequently, the responses obtained from Equation (13) should more closely resemble the true optimal performances intended (Appendix B). Figure 7(b) then shows the corresponding tip responses in time excited by a step motion input, illustrating now the relationship between and the settling time. It can be seen that the responses for (lower pair), which is , are already quite highly damped.

**(a)**

**(b)**

Overall, a Zener mount design methodology based on the performance indexes of and has been established to analytically determine the optimal P and D gains: Equations (17) and (13) confined by Equation (18), respectively. Here, the design methodology also represented by Figure 6 together with Approximation II and Conventions I and II is an important contribution of this paper. A merit of this is that the optimal gains are determined without actually solving the cubic equation in Equation (10b) subject to Equation (11). The P gain obtained is the exact solution (satisfying ), while the D gain is a suboptimal solution (satisfying ). The D gain can be very close to its exact solution if , which is often satisfied when the damping requirement is not highly demanded; for example, (Table 2 and Appendix B). It should be finally emphasized that the design methodology presented in this paper is also applicable to various forms of Zener mount design problems in the field of vibration isolation. Its application can further extend to various forms of vehicle suspension design problems because the model in Figure 3(a) is dynamically the same as the quarter-car model having a negligible unsprung mass [27]. Both models are indeed governed by the same characteristic equation given by Equation (10).

#### 4. Experiments

##### 4.1. Experimental Setup and Modeling

The experimental setup driving a single-link flexible beam in the horizontal plane is shown in Figure 8. An AC servo motor (Yaskawa® SGMJV-02ADA2C, 200 Watts) was directly connected to the beam (1 mm × 25 mm × 270 mm) of spring steel, without the use of a reduction gear. A built-in 20-bit incremental encoder equipped within the motor was used to measure the rotation angle. An accelerometer was additionally attached to the end tip of the beam to monitor the vibration. A servo drive (Yaskawa® SGDV 1R6A01A) was used in “torque control mode” so as to act as a current amplifier, producing a mechanical motor torque output in proportional to the electrical drive voltage input. Consequently, the back electromotive force effect of the motor could have been excluded in the servo drive mechanism. The diagram in Figure 8 is thus exactly analogous to that in Figure 1(a) but merely described in terms of voltages and the dimensionless gains, and . The control circuit was implemented within a real-time prototyping machine (Matlab/Simulink xPC target module) running at a sampling frequency of 16 kHz.

To model the mechanical manipulator, a random signal was supplied with the controller being deactivated (i.e., ). The time derivative (i.e., velocity) of the measured encoder angle response to this excitation, which is termed the driving point mobility, is compared in Figure 9(a) with its identification using Equation (4) with setting . Although the mobility is displayed by convention, the impedance (i.e., the inverse of the mobility data) is of importance for the curve-fitting in the LEM [12–14]. Thus, it can be seen that the first antipeak of the mobility has been closely fitted. The lumped elements identified are already tabulated in Table 1. The measured transmissibility is then compared in Figure 9(b) with the simulation using Equation (3). The measured response was obtained by dividing the tip angle motion by the hub angle motion, where the tip angle motion was obtained by dividing the tip translational displacement (double integral of the acceleration measured by the accelerometer) by the arm length. Here, the double integration means dividing the acceleration response by in the frequency domain. The angle motions can be regarded as the generalized displacements indicated in Figure 1(b). Note that the use of the accelerometer was not for the modeling but merely for monitoring the tip motion. It can be seen in Figures 9(a) and 9(b) that the simulated responses from the semidefinite LES agree well with the measured responses at very low and low frequencies. The peaks higher than the first one in Figure 9(b) are the higher modes of the beam [13], which were deliberately unmodeled in this paper as they are unimportant under Condition II.

**(a)**

**(b)**

##### 4.2. Experimental Performance

The three control cases listed in Table 2 were then experimentally tested. The values of “ from Equation (12)” in the table were solely used for the D gains throughout the control experiments. (This was because Equation (13) has been only recently developed while finalizing this paper.) The experimental results obtained are still useful to validate the complete LES model (now containing and ) and further to investigate any unforeseen dynamics and instability issues in real applications. The control behavior was first examined in the frequency domain by supplying a random signal to , while setting in Figure 8. To validate the LES model, the generalized hub and tip displacements (solid lines) for the case of are compared in Figures 10(a) and 10(b), respectively, with their predictions (dashed) calculated from Equations (5) and (6). They agree well at low frequencies particularly for the fundamental vibration mode at around 7 Hz. Some slight deviations at very low frequencies are thought to be due to friction in the motor that was neglected in the model but becomes influential for slow motion [28].

**(a)**

**(b)**

The measured hub and tip responses for all three cases are then compared in Figures 11(a) and 11(b) to clearly demonstrate the damping effect . An additional case was also tested using an arbitrarily large (N·m/rad) and a small (N·m·s/rad), and the results (thick dotted) are also overlaid to illustrate how large vibrations could be induced when the gains were inappropriately chosen [7, 13].

**(a)**

**(b)**

The performance of the control system was finally examined in the time domain but now with a step input of 20° to the reference input while setting . The hub displacements and the tip accelerations for the three test cases are compared in Figures 12(a) and 12(b), respectively. The additional case considered in Figure 11 could not be tested because of the limitation of the output voltage of the controller used. The trend in the step responses is similar to that in the frequency response functions in Figure 11. A small offset error is also evident when the P gain was small (i.e., a soft spring), for example, .

**(a)**

**(b)**

Throughout the control experiments described above, no instability was encountered. This was not because the experimental control system was genuinely passive like the theoretical LES model presented but because only moderate gains were implemented for the optimal PD controllers designed under Condition II. If very large gains had to be used for some other requirements, the system could indeed go unstable because of the time delay (approx. 0.5 ms: measured) inherent within the servo drive used. This time delay is also that can cause a rapid rigid manipulator system unstable [11].

Overall, the experimental results obtained have been presented to verify the LES model developed. It has also demonstrated that some friction-induced offset error can be persistent if a soft (a small ) manipulator system is operating slowly. Such error could be effectively reduced by using a large * or* completely eliminated by introducing an integral action within the controller, becoming then a PID controller. Although this paper is specifically concerned with the PD controller of a linear flexible manipulator system, it should be emphasized that the optimal P and D gains presented can still be used, at least, as initial values of an intelligent PID controller for a flexible manipulator of unknown friction dynamics [29]. More fundamentally, the mathematical logic and physical insights presented in this paper provide a basic foundation upon which an analytical PID-tuning method can be established in the future for general nonlinear manipulator systems.

#### 5. Conclusions

A seemingly complex but fundamental problem in the field of motion control has been tackled, which is to analytically tune the PD controller of a linear flexible manipulator system. More specifically, the problem is to give simple closed-form solutions of the optimal P and D gains to yield the maximum bandwidth under a given damping requirement *or* conversely the maximum damping under a given bandwidth requirement of the system in free oscillation. It has turned out that the original servo control problem is analogous to a conventional Zener mount design problem. A design methodology has been finally established to analytically determine the optimal elements of the mount, corresponding to the optimal control gains. The P gain obtained is the exact optimal solution while the D gain is a suboptimal solution that can also be very close to its exact solution if the damping requirement is not highly demanded. The D gain can further be nearly exactly tuned, regardless of the damping requirement, if a numerical frequency tuning method is allowed. Mathematically, it should be emphasized that the methodology established in this paper offers an optimal design solution to a third-order dynamic system. Practically, it can thus be applied to various dynamic systems including Zener mounts and vehicle suspensions in the field of vibration isolation.

#### Appendix

#### A. Some Important Frequencies

With reference to the system in Figure 1(b), the appendage and base angular natural frequencies are, respectively,With reference to the system in Figure 3(a), the nominal and standard angular natural frequencies of the Zener mount system are, respectively,where is the static stiffness while is the dynamic stiffness at the transition frequency of the Zener mount given bywhere , , and . With reference to the frequency responses in Figures 2, 7(a), 10, and 11, the first peak frequency of each response is defined as the resonance frequency of the fundamental vibration mode of the corresponding servo system in operation under a certain *steady-state* excitation. With reference to the time responses in Figures 7(b) and 12, the fundamental residual vibration frequency is defined as the fundamental damped natural frequency of the servo system in *free* oscillation.

#### B. Alternative Approaches Using Analytical Convex Minimization and Numerical Frequency Tuning

The two design cases studied in this Appendix are depicted in Figures 13(a) and 13(b), where the system in Figure 13(a) is subjected to a motion input on both and , while that in Figure 13(b) is on only. The system in Figure 13(b) is the same as that considered in the main text. The dynamic equations of the first case can be written in the Laplace domain asThe vibration transmissibility function from to is then given bywhere , , and . Given the value of , the design problem here is to find and that minimize the resonant vibration amplitude of such thatwhere in which and . Under the assumption that the function has the global minimum, we can apply the convex minimization scheme [30] to get and consequently . It finally turns out that the resonance frequency solution is , which is *surprisingly* the same as the standard natural frequency defined in the main text. We can then apply to get and consequently that is the same as Equation (13). The global minimum for each is then given bywhere it can be seen that it is always for all positive .

**(a)**

**(b)**

It is further interesting that the vibration transmission problem here through a series-type Zener mount is exactly analogous to the force transmission problem through a parallel-type Zener mount widely studied in a textbook [30] and also in Reference [17]. In these standard cases to which the convex minimization scheme is applicable, the conventional design rule () turns into the *standard* design rule defined in Section 3.2 as

Under the optimal design framework defined mathematically by Equation (B.3) or intuitively by Equation (B.5), not only the P gain in Equation (17) but the D gain in Equation (13) is also the exact solution for the particular system shown in Figure 13(a). Strictly speaking, this framework is concerned with the steady-state vibration under a certain excitation condition while the *general* design rule given by Equation (11) is with any residual transient vibration. However, these two rules often yield similar results in many practical applications. This might be a reason why they were not clearly distinguished from one another in the literature [18, 19].

Similarly, the transmissibility function of the second case in Figure 13(b) can be written asApplying again now gives a trivial solution, which means that the convex minimization scheme is no longer applicable. In this nonstandard case, some approximate solutions would be more appropriate as illustrated below with simulations.

The standard and general design rules, respectively, given by Equations (B.5) and (11), were applied to the two systems in Figure 13 for the three design cases listed in Table 2. The P gain was analytically determined by Equation (17) for both rules. The D gain of the standard rule was then analytically determined by Equation (13), which is also the solution of the convex optimization method earlier. However, that of the general rule was numerically (i.e., iteratively) tuned until it precisely met the frequency tuning condition in Equation (11), which is a numerical frequency tuning method. The optimal D gains finally tuned (i.e., “ satisfying Equation (11)”) are also added to Table 2 for comparison. It can be seen that if either of the two rules is applied and if the damping requirement is not highly demanded, we haveand consequently the values of from Equations (13) and (11) are also very similar to each other, justifying Approximation II of the main text.

Using the values of and obtained as tabulated in Table 2, the generalized tip displacements of the two systems in Figures 13(a) and 13(b) could finally be calculated as shown in Figure 14. Those fulfilling the general (solid lines) and standard (dashed) rules are displayed together for comparison. Some symbols are additionally indicated to confirm that the rules have been correctly applied. It can be seen that (indicated by “∇” in the graphs of the -only excitation case, for convenience) and (indicated by “+” in the graphs of the excitation case) coincide with their corresponding values of (indicated by “ο” in the graphs of both cases), respectively. It can also be seen that little is discernable between each pair (solid and dashed lines) of the responses, particularly when is small. This supports that the simple standard rule () can be used instead of the general rule (), if the damping requirement is not highly demanded. Thus, the analytical design methodology presented in this paper can be applied to many practical residual as well as random vibration control problems, where so that .

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

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

#### Acknowledgments

This research was commenced in 2011 and was partly supported by the Leading Human Resource Training Program of Regional Neo Industry through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT, and Future Planning (No. 2016H1D5A1910490).