I a ‘, A . _ . Fa...II.ul._..III_II.-.IIIIIIIIIIlg.w ,. qagwo'x :> GAN STATE UNIVERSITY LIBRARIES IIIII III III II II III III III 3 1293 006294 IL l I ' LIBRARY Michigan State I, University This is to certify that the dissertation entitled CONTROL OF SINGLE-LINK FLEXIBLE MANIPULA’I‘ORS FABRICATED FROM ADVANCED COMPOSITE IMNATES AND SMART MATERIALS INCORPORATING ELECI'RO-RHEOWICAL FLUIDS presented by SEUNG-BOK CHOI has been accepted towards fulfillment of the requirements for Ph-D. degree in M-E- CI \/ GonoIW/JFQ/E Major professor Date /fl% /7¢0 MS U i: an Affirmative Action/Equal Opportunity Institution 0- 12771 PLACE IN RETURN BOX to remove this checkout from your record. TO AVOID FINES return on or bdore dete due. DATE DUE DATE DUE DATE DUE MSU Is An Affirmative ActIorVEquel Opportunity Institution CONTROL OF SINGLE-LINK FLEXIBLE HANIPUIATORS FABRICATED I-‘ROH ADVANCED COHPOSITE IAHINATES AND SMART MATERIALS INCORPORATING ELECTRO-RHEOLOGICAL FLUIDS by Seung-Bok Choi A DISSERTATION Submitted to Hichigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHIIDSOPHY Department of Hechanical Engineering 1990 boSdQQc ABSTRACT CONTROL OF SINGLE-LINK FLEXIBLE MANIPUIATORS FABRICATED FROM ADVANC- COMPOSITE LAMINATES AND SMART MATERIALS INCORPORATING ELECTRO-RHEOLOGICAL FLUIDS 3? Seung-Bok Choi Three different control schemes for a single-link flexible manipulator fabricated from composite laminates or smart materials incorporating electro-rheological fluids are developed in this study ; an output feedback controller, a hybrid controller, and a nonlinear state feedback controller. The output feedback controller featuring two colocated sensors is designed by utilizing the root-locus technique and the Lyapunov method after a description of the dynamic modeling of the flexible manipulator. It is shown in the analytical and experimental investigation that the arm fabricated from composite laminates has superior system performances such as faster settling-time, smaller input torque, smaller overshoot, and superior stability characteristics relative to the arm fabricated from aluminum. The hybrid controller which consists of an output feedback controller and a pseudo-state feedback controller is proposed for the flexible manipulator fabricated from smart materials incorporating electro-rheological fluids. A phenomenological equation ofrmnflon is developed from experimental observations of the modal characteristics. .And an empirical model of a distributed electro-rheological fluid actuator which characterizes the pseudo-state feedback controller is obtained as an explicit function of the electrical intensities. It is shown in the experimental implementation that the robustness of the output feedback controller subjected to a disturbance which has a forcing frequency close to the clamped-mode natural frequency of the arm is impressively enhanced by employing the proposed pseudo-state feedback controller. I The nonlinear state feedback controller, which is robust to system parameter variations caused by uncertain parameters such as varying 'payload and mismatched lamina angle, is developed for the flexible manipulator fabricated from composite laminates. The effects of each uncertain parameter as well as the interaction of all parameters are considered without assuming matching conditions. The properties of the uniform and uniform ultimate boundedness of the solution of system state equations are employed to formulate the controller. It is demonstrated through the simulation that the proposed nonlinear feedback controller has much better capability of arresting the unwanted vibratflxnn than the linear controller in the presence of the uncertain parameters. iii This thesis is dedicated to the memory of my father, Kyung-Hyun Choi (1909-1988) and father-in-law, Bong-Nam Kim (1931-1989), who passed away during the course of this research program. Their immense love and sacrifices remain forever. iv ACKNOWLEDGMENTS I would like to express my deepest gratitude and appreciation to my co-advisors, Dr. Brian S. Thompson and Dr. Mukesh V. Gandhi, for their expert guidance, friendship and continuous encouragement throughout the course of this research program, and for providing the unique research environment. Sincere appreciation is extended to the other guidance committee, Dr. Steven.Wfl Shaw and Dr. David H.Y. Yen, for their valuable comments and suggestions. I also thank Dr. Clark Radcliffe for his constructive comments on this research. Particular gratitude goes to my lovely wife Yeon-Ook, my smart son Min-Kyu and my pretty daughter Bun-Young for being a wonderful family with patience and understanding. I, of course, wish to express sincere gratitude to my mother, Bo-Bae Cho, mother-in-law, Soon-Sun Jang, eldest brother, Seung-Bae and his wife Soon-Lim Lee, elder brother, Seung-Keun and his wife Keun-Soon Kim, younger sister, Seung-Yi and her husband Young-Hyun Kim, brother-in-law, Ha Kim and his wife Soon-Lim Bee, and sister-in-law, Hyo-Sook Kim, for their providing moral support and encouragement. Special thanks go to my friends at Michigan State University for making my graduate career a successful experience. Also, I would like to thank all the students in Machinery Elastodynamics Laboratory and Intelligent Materials and Structures Laboratory for their technical ass is tance . Finally, I wish to gratefully acknowledge the financial support of this research program by the Defence Advanced Research Projects Agency, under contract DAAL03-87-K-0018, the U.S. Army Research Office, under contract DAAL03-K-0022, and also the State of Michigan, Department of Commerce, Research Excellence and Economic Development Fund. vi LIST OF LIST OF TABLE OF CONTENTS TABLES ........................................... FIGURES .......................................... CHAPTER I. INTRODUCTION ................................. 1.1 Motivation ....................................... 1.2.1 Dynamic Modeling ............................. 1.2.2 Controller Design ............................ 1.2.3 Related Works on Modeling and Active Control of Large Flexible Structures ................. 1.2.4 Related Works on Smart Materials and Structures 1.3 1.4 CHAPTER 2.1 2.2 2.3 2.4 2.5 CHAPTER Thesis Organization ..I ............................ Principal Contributions .......................... II. MODELING OF A SINGLE-LINK FLEXIBLE MANIPULATOR FABRICATED FROM ADVANCED COMPOSITE LAMINATES Introduction ..................................... Derivation of Equations of Motion ................ Finite Element Formulation ....................... State Space Representation of the Model .......... Summary of the Chapter ........................... III. EXPERIMENTAL APPARATUS AND IDENTIFICATION OF SYSTEM MODEL PARAMETERS .................... Introduction ..................................... Experimental Apparatus ........................... 3.2.1 Flexible Arm ................................. .2.2 DC Motor and Servo-Drive 3.2.1.1 Aluminum Arm ............................. 3.2.1.2 Composite Arm ............................ vii eeeeeeeee 000000000 ......... 000000000 000000000 ......... 000000000 19 26 36 38 41 41 42 50 S3 S6 57 57 57 59 59 61 66 3.2.3 Servo-Positioning Assembly ............................ 67 3.2.4 Sensors ............................................... 68 3.2.4.1 Angular Position Sensor ........................... 68 3.2.4.2 Angular Velocity Sensor ........................... 68 3.2.4.3 Tip Deflection Sensor ............................. 69 3.2.5 Microprocessor ........................................ 70 3.3 Identification of System Model Parameters ................. 71 3.3.1 Predicted ............................................. 72 3.3.2 Measured .............................................. 81 3.4 Summary of the Chapter .................................... 87 CHAPTER IV. OUTPUT FEEDBACK CONTROL ASSOCIATED WITH TWO COLOCATED OUTPUT SENSORS ............................. 88 .1 Introduction .............................................. 88 .2 Controller Design ......................................... 89 .3 Implementation Results and Discussions .................... 99 4.3.1 Step Responses without Disturbances ................... 101 4.3.2 Step Responses with Disturbances ...................... 117 4.3.3 Instability of the System ............................. 128 4.4 Summary of the Chapter .................................... 136 CHAPTER V. HYBRID CONTROL FEATURING A DISTRIBUTED ELECTRO- RHEOLOGICAL FLUID ACTUATOR ............................ 138 5.1 Introduction .............................................. 138 5.2 Background on Electra-Rheological Fluids .................. 139 5.3 Phenomenological Equation of Motion ....................... 144 5.4 Controller Formulation .................................... 151 5.5 Experimental Implementation ............................... 162 5.5.1 Apparatus and Instrumentation ......................... 162 5.5.2 Results and Discussions ............................... 170 5.6 Summary of the Chapter .................................... 181 (HflAPTER VI. NONLINEAR STATE FEEDBACK CONTROL ROBUST TO SYSTEM PARAMETER VARIATIONS .......................... 183 6.1 Introduction .............................................. 183 6.2 Problem Formulation ....................................... 185 viii 6.3 Controller Design ......................................... 190 6.4 Illustrative Example ...................................... 209 6.5 Summary of the Chapter .................................... 222 CHAPTER VII. CONCLUSIONS AND RECOMMENDATIONS ..................... 224 7.1 Conclusions of the Thesis ................................. 224 7.2 Recommendations for Future Work ........................... 227 APPENDIX .......................................................... 231 List of Publications ...................................... 231 BIBLIOGRAPHY ...................................................... 236 ix LIST OF TABLES Page Table 3.1 Material and Geometrical Properties of the Aluminum Arm .......................................... 60 Table 3.2 Material and Geometrical Properties of the Composite Arm ......................................... 64 Table 3.3 Comparison of the System Model Parameters of the Aluminum and Composite Arm without Payload ............ 73 Table 3.4 Comparison of the System Model Parameters of the Aluminum and Composite Arm with Payload of 0.061 kg .... 73 Table 3.5 Comparison of the Calculated and Measured Natural Frequency without Payload ............................. 85 Table 3.6 Comparison of the Calculated and Measured Natural Frequency with Payload of 0.061 kg .................... 85 Table 4.1 Closed-Loop Poles with the Compensator Zero Zco- -2 .... 102 Table 4.2 Closed-Loop Poles with the Compensator Zero 200- -5 .... 107 Table 4.3 Closed-Loop Poles with the Compensator Zero Zco- -1 .... 111 Table 4.4 Closed-Loop Poles with Kp-0.406, Kv-0.254 for the Aluminum Arm, and Kp- 0.403, Kv- 0.252 for the Composite Arm ...................................... 114 Table 5.1 Coherence Factor of the First and Second Mode at Different Voltage ..................................... 149 Table 5.2 Measured System Model Parameters of the Smart ER Arm in the Absence of the Electric Field .................. 169 Table 6.1 Geometrical and Material Properties of the Nominal Composite Arm ................................. 210 Table 6.2 System Model Parameters of the Nominal System ......... 212 Table 6.3 Comparison of the Tip Displacement Response (Settling Time and Maximum Overshoot) ........................... 220 Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure .10 .11 .12 .13 .14 .15 LIST OF FIGURES Schematic Diagram of a Single-Link Flexible Manipulator Composite Laminate Configuration Two-Node Beam Element Photograph of Overall Robotic System Schematic Diagram of Overall Robotic System Photograph of the Experimental Aluminum Arm Design Methodology of the Composite Arm Curing Cycle for AS4/3501 Graphite-Epoxy Photograph of the Experimental Composite Arm Load-Deflection Curve for the Aluminum and Composite Arm Relationship between the Velocity and the Measured Voltage of the Tachometer Relationship between the Amplified Voltage and the Tip (End-Point) Deflection Frequency Variation with respect to the Ratio of the Hub Moment of Inertia to the Beam Moment of Inertia .... Frequency Variation with respect to the Ratio of the Payload to the Beam Mass First Three Pinned-Mode Shape for the Aluminum and Composite Arm Calculated Open-Loop Transfer Function for the Angular Velocity without Payload Calculated Open-Loop Transfer Function for the Angular Velocity with Payload of 0.061 kg Calculated and Measured Frequency Response of the Aluminum Arm without Payload xi 00000000000000000 OOOOOOOOOOOOOOOOOOOOOOO 51 58 58 59 62 63 65 66 69 7O 76 77 78 79 80 83 Figure 3.16 Calculated and Measured Frequency Response of the Composite Arm without Payload ........................ 84 Figure 3.17 Calculated Open-Loop Transfer Function for the Angular Velocity with é(s)-(s/(s/a+l))9(s) ................... 86 Figure 4.1 Block-Diagram of the Output Feedback Control System ....................................... 89 Figure 4.2 Open-Loop Pole-Zero Location with the Compensator Zero zco- -1 ........................................ 91 Figure 4.3 Root-Locus with the Compensator Zero Zoo. -1 ......... 92 Figure 4.4 Open-Loop Pole Zero Location with the Compensator Zero 2 - -15 ....................................... 94 co Figure 4.5 Root-Locus with the Compensator Zero ZCO- -15 ........ 95 Figure 4.6 Experimental Set-Up for Output Measurements ........... 99 Figure 4.7 Circuit Diagram of Two Channel Operational Amplifier ............................................ 100 Figure 4.8 Measured Step Responses without Payload. Kp-0.534, KVP0-257 for the Aluminum Arm, and Kp-0.292, KV-0.l46 for the Composite Arm ............. 103 Figure 4.9 Simulated Step Responses without Payload. Kp-0.534, Kv-0.267 for the Aluminum Arm, and Kp-0.292, Kv-0.146 for the Composite Arm ............. 104 Figure 4.10 Simulated Open-Loop Step Responses of the Aluminum Arm without Payload .................................. 105 Figure 4.11 Measured Step Responses without Payload. Kp-1.0 and Kv-0.2 for both Aluminum and Composite Arms ....... 108 Figure 4.12 Simulated Step Responses without Payload. Kp-1.0 and Kv-0.2 for both Aluminum and Composite Arms ....... 109 Figure 4.13 Measured Step Responses with Payload of 0.061 kg. Kp-Kv70.48 for the Aluminum Arm, and Kp-Kv-0.24 for the Composite Arm .................................... 112 xii Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure .14 .15 .16 .17 .18 .19 .20 .21 .22 .23 .24 Simulated Step Responses with Payload of 0.061 kg. Kp-Kv-0.48 for the Aluminum Arm, and Kp-Kv=0.24 for the Composite Arm OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO Measured Step Responses without Payload. Kp-0.406, Kv-0.254 for the Aluminum Arm, and Kp-0.403, er0.252 for the Composite Arm 0000000000000 Simulated Step Responses without Payload. Kp-0.406, Kv90.254 for the Aluminum Arm, and Kp-0.403, Kv-0.252 for the Composite Arm Block—Diagram of the PD Output Feedback Control System with the Disturbance of d(t) Measured Step Responses of the Aluminum Arm without Payload. Kp- 0.534, Kv- 0.267, and 3(t)-0.08 sin(3.25 Hz)t Comparison of the Measured Step Responses of the Aluminum Arm with and without the Disturbance. xp-o.534, Kv-0.267 and d(t)-0.08 sin(3.25 Hz)t ....... 121 Comparison of the Simulated and Measured Step Responses of the Aluminum Arm with the Disturbance. Kp-0.534, Kv-O.267 and d(t)-0.08 sin(3.25 Hz)t ....... 122 Measured Step Responses of the Composite Arm without Payload. Kp-0.292, Kv-0.146 and 3(c)-0.08 sin(4.75 Hz)t Comparison of the Measured Step Responses of the Composite Arm with and without the Disturbance. Kp-0.292, xveo.145, and 3(t)-0.08 sin(4.75 Hz)t ...... 125 Comparison of the Simulated and Measured Step Responses of the Composite Arm with the Disturbance. xp-o.292, Kv-O.146 and d(t)-0.08 sin(4.75 Hz)t ....... 126 Measured Step Responses of the Composite Arm without Payload. Kp- 0.292, Kv- 0.146, and 3(t)-o.4 sin(34.5 Hz)t ................................ 127 xiii Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure .25 .26 .27 .28 .29 .30 .10 .11 .12 .13 .14 In-Position Band of the Servo-Drive System ........... 129 Measured Step Responses of the Aluminum Arm without Payload. Kp-6.2 and Kv-0.4 ................... 131 Simulated Step Responses of the Aluminum Arm without Payload. Kp-6.2 and Kve0.4 .................. 132 Measured Step Responses with Feedback Gains of Kp-6.2 and KV-0.4 for both Aluminum and Composite Arms ...... 133 Measured Step Responses of the Composite Arm without Payload. Kp- 11.8 and Kv- 0.4 ............... 134 Simulated Step Responses of the Composite Arm without Payload. Kp- 11.8 and Kv- 0.4 ............... 135 Electra-Rheological Fluid in Liquid and Solid State ...................................... 140 Interaction between Charges on Electrodes and Those in Electro-Rheological Fluid ......................... 141 Photomicrograph of Electra-Rheological Fluid .......... 142 Smart Cantilevered Beam Structure .................... 146 Mode Shapes of the Smart Cantilevered Beam ........... 147 Measured Frequency Response of the Smart ER Arm without Payload ...................................... 155 Governing Characteristics of the Employed Electro- Rheological Fluid Actuator ........................... 157 Block-Diagram of the Hybrid Control System ........... 159 Bang-off-Bang Type Input Voltage ..................... 161 Experimental Set-Up for the Frequency Response Measurement of the Smart Cantilevered Beam ........... 163 Frequency Response of the Smart Cantilevered Beam ..... 165 Schematic Diagram of the Smart ER Arm ................ 166 Flow Diagram for Fabrication Procedures of the Smart ER Arm ......................................... 167 Photograph of the Smart ER Arm ....................... 168 xiv Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure Figure .15 .16 .17 .18 .19 .20 .21 .22 Measured Step Responses with Feedback Gains of K - 0.5 and K - 0.25 .............................. 172 p v Simulated Step Responses with Feedback Gains of K - 0.5 and K - 0.25 .............................. 173 p v Measured Step Responses with Feedback Gains of K - 1.4 and K - 0.25 .............................. 174 p v Simulated Step Responses with Feedback Gains of K - 1.4 and K - 0.25 .............................. 175 p v Measured Step Responses with Bang-off-Bang Type Input Voltage. Kp-1.4 and Kv-0.25 ......................... 176 Measured Step Responses with the Disturbance. Kp- 0.5, Kve 0.25 and a(:)-o.08 sin(4.1 Hz): ......... 178 Simulated Step Responses with the Disturbance. Kp- 0.5, Kv- 0.25 and a(:)-o.08 sin(4.1 Hz): ......... 179 Measured Step Response with the Disturbance. xp- 0.5, Kv- 0.25 and a(:)-o.a sin(16.8 Hz): ......... 180 Parameter Variation of the Composite Arm with respect to the Payload .......................... 211 Parameter Variation of the Composite Arm with respect to the Lamina Angle ..................... 211 Open-Loop Tip Displacement of the Nominal System ...... 216 Comparison of Step Responses with the Linear and the Nonlinear Controller ........... .............. 218 Comparison of Step Responses with (-0.1 and e-0.003 ... 221 XV CHAPTERI INTRODUCTION 1 . 1 Motivation Most of current day industrial robots depend on bulky design and can handle objects no heavier than five percent of their own weight. The bulkiness is mainly due to an inherent design requirement to minimize structural vibration by increasing the mechanical stiffness of each component. This massive structural design makes the robot slow and heavy; hence, requires large actuators and high mounting strength and rigidity, which consequently leads to high energy consumption and high overall operating cost [26]. The insatiable demand in the international marketplace for high performance articulating robotic systems has triggered a vigorous research thrust in various multi-disciplinary areas such as design and control. The high performance of robots may be quantified by high-speed of operation, high end-point positional accuracy, lower or elimination of friction and backlash incorporating indirect drive mechanism, and lower overall cost. In order to achieve these high performances, light-weight, flexible robotic arms have been proposed and developed. One can simply reduce the weight of the links of the robot in order to achieve these goals. However, it is not possible to move lightweight arms quickly without the onset of structural vibration which is a manifestation of inadequate structural damping as well as lower 1 2 structural strength and stiffness. Attempts have been and are being made to overcome these deficiencies through optimization of the stiffness-to- weight ratio by varying the length and cross section of each metallic link [67]. A more innovative approach employs advanced composite materials in the construction of high-speed manipulators to achieve larger structural damping and higher stiffness—to-weight ratio [38, 118, 145, 146, 147]. Recent studies [38, 146] have demonstrated both analytically and experimentally that transverse deflection can be reduced substantially when links are made of advanced composite materials. Moreover, to increase structural damping, optimization approaches have been developed [82] . This is usually called a passive damping method or passive control. However, present studies indicate that the increased damping due to passive means is not sufficient by itself to eliminate structural deflection [24]. Therefore, active control method, i.e. , the design of control algorithms that take into consideration of the flexible dynamics is needed. A variety of dynamic modeling and control strategies for flexible robotic manipulators, fabricated from conventional isotropic materials such as steel and aluminum have been developed, and a few of them have been experimentally implemented as described in the literature review. However, up to now, the research areas on dynamic modeling and control of flexible robotic manipulators fabricated from advanced composite materials are considerably rare. This study addresses the dynamic modeling and feedback control of a single-link flexible manipulator fabricated from advanced composite materials by undertaking an analytical and experimental investigation. The advantages of employing composite materials relative to the 3 conventional isotropic materials can be divided into three categories; design, control, and economics and safety. The advantages of the design area are well presented in [145, 146]. The advantages of control area may be described as a larger bandwidth due to the high stiffness-to-weight ratio of the material, an increase of the stability margin due to larger structural damping, a decrease of the control spillover effects due to the inherent high natural frequency of the unmodelled dynamics. Furthermore, more appropriate sampling time associated with observed and controlled modes can be achieved by adjusting the high frequency and large damping properties during the implementation [24] . The advantages of the last category will include lower energy consumption, lower overall cost due to high productivity on rate, and safer operation due to reduced inertia forces. Some of these advantages relating control and economic area will be observed both analytically and experimentally in this study. It is also addressed in this study that the closed-loop system with the conventional output feedback controller associated with two colocated output sensors is very sensitive to an external disturbance which.has a. forcing frequency close to the cantilevered-mode natural frequency of the flexible arm. This attempt is made here to introduce a single-link flexible manipulator fabricated from smart materials incorporating electro-rheological fluids. The. robustness of the conventional output feedback controller to the disturbance is to be expected by employing the distributed parameter actuator which is characterized by embedded electro-rheological fluids in the structure. It has been verified [41,56] that the vibrational characteristics such as natural frequencies and damping capacities of the flexible structures incorporating electro- rheological fluids can be controlled in a very short time by imposing 4 electric potentials, hence providing the possibility of the elimination of resonance phenomenon of the structures or systems. Another important issue to be considered for controlling of flexible robotic manipulators is system parameter variations such as natural frequencies which may always arise in practice due to a wide spectrum of various conditions associated with manufacturing of flexible arms, dynamic modeling, and operating environments. The parameter variations are more serious in the robotic systems featuring composite-based flexible arms than those in the robotic systems featuring metallic-based flexible arms. The reason is that the composite-based arms involve more design variables, more complicated manufacturing process, and also higher sensitivities to environmental conditions such as temperature and moisture. In order to respond this problem, a nonlinear state feedback controller robust to parameter variations caused by the varying payload which may be encountered in the operation, and the mismatched lamina angle which may be encountered in the manufacturing process, is developed in this study. These proposed observations which may provide significant impact on the diverse marketplace such as automobile industry and aerospace industry, coupled with the dearth of relevant research areas on the control of flexible robotic manipulators fabricated from advanced composite materials and smart materials have motivated the comprehensive research program reported in this study. 1 . 2 Literature Review Although this research program discusses on the dynamic modeling and control of a single-link flexible manipulator, some of the major 5 contributions on the dynamic modeling and control of multi-link flexible manipulators will be also reviewed here. FUrthermore, some research works on the modeling and active vibration control of the large flexible structures will be reviewed in order to provide a flavor of applicability of the proposed research to the relevant structural systems. Finally, the literature on the smart materials and smart structures will be reviewed in order to introduce the proposed single-link flexible manipulator fabricated from smart materials incorporating electro-rheological fluids. It is tuited that the materials employed for the fabrication of flexible arms in the following references are assumed to be conventional isotropic materials such as aluminum and steel unless specific statement is made. 1 . 2 . 1 Dynamic Modeling Schmitz [127] derived governing equations of motion and associated boundary conditions for a single-link flexible manipulator by employing the Hamilton's principle. The flexible link assumed to be governed by Euler-Bernoulli beam theory. The characteristic equation for the determination of natural frequencies and the eigenfunction for the determination of mode shapes are obtained for both cantilevered-mode and pinned-mode. Two different reduced-order models; reduced-order model with system modes and cantilevered modes are derived from finite modal expansion. The system model parameters such as natural frequencies and modal slope coefficient are verified through experimental identification, and shows good agreement between two results. Low [88,89,90] applied the Hamilton's principle in order to derive the coupled, nonlinear integro- t1), amnd the virtual displacements 6w(x,t) and 66(t) are independent and taken equal to zero at x-O and x=L. 2.3 Finite Element Formulation There are several ways to solve the governing equation (2.18) in order to obtain system model parameters such as natural frequencies and mode shapes which will be used for designing of feedback controllers. Herein, a finite element formulation incorporating Euler-Bernoulli beam theory is developed to obtain system model parameters. 12ma T ~ —1— 1 . . . n B I o 1 0 dx 0 dx (2.34) T 1 o d¢1(0) 0 . . . d¢n(0) 0 dx dx 6 _ o 1 o d¢1(o). . . o d¢n(o) dx dx L L o ¢1(L) o - - - ¢n(L) o . 55 In equation (2.34), mi and Ci are the natural frequency and damping ratio of the ith mode respectively, I is the total moment of inertia (hub, T beam and payload), d¢i(o)/dx is the modal slope coefficient and ¢i(L) is the tip sensor gain. It is noted that the output matrix ‘6' includes the gains associated with two colocated sensors - hub angle (first row) and hub angular velocity (second row), and one noncolocated sensor - tip position (third row). Since the dynamical model given by the equation (2.33) contains many modes in general, all of modes cannot be controlled in practice. Hence, the equation (2.33) should be decomposed into two subdynamic equations by considering the rigid-body mode and the first few critical flexible modes as primary modes, and the remaining modes as residual modes. Two subdynamic equations are Si-KSE+EG P P P P ~ ~ ~ (2.35) - C x yP P P and 31-3.; +§E r r r r ~ ~ ~ (2.36) yr - Crxr where subscript (p) denotes the primary modes and (r) represents the residual modes. In general, there is no certain rule differentiating the primary and residual modes. The differentiation heavily depends upon the geometrical and material properties of the flexible arm and the dynamic characteristics of the employed experimental apparatus such as actuators and sensors [74]. 56 2.5 Summary of the Chapter The governing partial differential equation and associated boundary conditions for the single-link flexible manipulator fabricated fnmn composite laminates were derived from the Hamilton's principle, and the equation was discretized using the finite element method incorporating the Euler-Bernoulli beam theory. Then, the decoupled linear dynamic model was (flytained.through the modal analysis and presented in the state space representation, which is convenient form for controller design. It is ‘noted that the relative simplicity of the dynamic model for the proposed composite arm disappears when considering multi-link flexible manipulators. However, the approach employed in this study could be extended to the case of multi-link flexible manipulators or other complicate flexible structures such as helicopter blades and appendages of the satellites fabricated from advanced composite materials. (HHHHER.III EXPEKHflflHml.APHMMUflEIANDIUMQHHITCAIION OF SiEflfllHDDELIHMMflflflEBS 3.1 Introduction A detailed description of the experimental hardware employed for the proposed single-link flexible robotic system featuring the aluminum or the composite arm will be presented in this chapter. Also, an experimental identification of the system model parameters such as natural frequencies and damping ratios will be undertaken. The measured parameters will be compared with the predicted ones. This chapter establishes a point of reference and sets a physical scale of the experiment for the subsequent chapter, which presents the design and implementation of the output feedback controller associated with two colocated sensors for the proposed robotic system. 3.2 Experimental Apparatus The experimental hardware of the proposed single-link flexible robotic system can be separated into five functional groups; flexible arms (aluminum or composite arm), a DC motor and a servo-drive, a servo- positioning assembly, sensors, and a microprocessor for controller implementation. Figure 3.1 presents a photograph showing these functional groups, and Figure 3.2 shows a schematic diagram of overall robotic 57 58 Figure 3.1 W Photograph of Overall Robotic System Flexible Arm Payload Strain Gage + Servo-Drive f Tachometer Status Servo-Positioning Assembly l Figure 3.2 Decoder x 4 Commanded Microprocessor Position Incremental Encoder Schematic Diagram of Overall Robotic System 59 system consisting of these functional groups. In this section, The detailed description of each functional group will be presented. 3 . 2 . 1 Flexible Arm 3.2.1.1 Aluminum Am The aluminum arm is a 0.812 m long, flexible beam which can be easily bended in the horizontal plane, but it is designed to be hard to bend in the vertical and torsional direction to prevent or reduce the torsional effects which were neglected in the modeling process. The arm is clamped on a rigid-body hub which has a radius of 4.3 cm and a moment 2 of inertia of 0.00496 kg-m, and is mounted directly on the vertical Figure 3.3 Photograph of the Experimental Aluminum Arm 6O shaft of the DC servo-motor as shown in Figure 3.3. Therefore the rotatfixnl of the arm in the horizontal plane is activated by the torque applied to the motor without speed reduction factor. 'It is seen from Figure 3.3 that the two extreme travel limit switches and one home position limit switch are mounted on the plate of the fixture supporting; the motor. The purpose of the installation of the limit switches is prevent overtravel in the emergency case such as instability of the system. When the arm touches one of the limit switches, the main power for the whole system will be turned off automatically. The home position switch is installed.in.order'to obtain high repeatability of the experimental results. With this home position detective switch tfiua'robot arm can always start to move from the same position like zero-position in typically industrial robotic manipulators. The payload made of aluminum is attached to the free end of the arm. The 30 mm (length) x 40 mm Table 3.1 Material and Geometrical Properties of the Aluminum Arm Properties Specifications Beam Material A2 606-T6511-QQA200/8 Young's Modulus 64.3 GPa Beam Length 0.812 m Beam Width 19.02 mm Beam Thickness 3.14 mm Beam Mass 0.132 kg 61 (width) x 8.56 mm (thickness) size of payload is used and its mass is 0.061 kg. Table 3.1 presents the material and geometrical properties of the experimental aluminum arm. 3.2.1.2 Composite Arm The design and fabrication of the composite arm is not as simple as the case of the aluminum arm. Once the geometrical and material properties of the aluminum arm has been determined, the design of the composite arm can be started based on some design constraints and design criteria which depend upon the objective of experimental investigation and the adaptability of the experimental apparatus. Herein, since the objective of the experimental investigation is to compare quantitatively the control performances between the aluminum and composite arm, following design constraints and a criterion are imposed. Design Constraints: beam length(L); (L)al - (L)com (3.1) beam w1dth(b) ; (b)a1 - (b)com Design Criterion; flexural rigidity(FR); (FR)al - (FR)com (3.2) where subscripts (a1) denotes the aluminum arm and (com) represents the composite arm. The hub which was employed for the aluminum arm can be easily used for the composite arm with the constraint of same width. The selection of the design criterion, i.e., flexural rigidity (bending stiffness) is somewhat rational way rather than logical manner. Someone 62 evaluation of design criterion for the aluminum arm ] selection of composite materials V imposition of design constraints 7 determination of design variables: _‘ fiber orientation and no. of ply V evaluation of design criterion for the composite arm comparison of design criterion unsatisfactory satisfactory r- -------- 1 ..... .. p: fabrication ‘ |—-________ Figure 3.4 Design Methodology of the Composite Arm 63 can choose other design criteria such as beam geometry. However, it seems to be reasonable start to choose such design constraints and a criterion, since it is well-known that the composite materials has higher strength- and stiffness-weight ratio than that of the aluminum, hence resultilu; in significant reduced inertia forces in the dynamic motion of the system. Once the design constraints and design criterion have been decided, the design procedure can be accomplished based on a design methodology presented in Figure 3.4. The design variables are restricted to the fiber orientatdxnl (lamina angle) and the number of the prepreg ply, since the manufacturing company generally provides the prepreg which has a certain. fixed fiber volume fraction. The flexural rigidity of the composite arm 300 — PRESSUR=85 PSI, VACUUM=28 IN HG * PRESSURE=1OO PSI, VACUUM=VENT (Do 200— if - ................................. 0: D l— 4 <1 0: “J 1 “2- L13 100-] O ‘Y I I I Y Y j r Y r Y 1 f l I I T Y I 0 60 120 180 240 300 TIME(sec) Figure 3.5 Curing Cycle for AS4/3501 Graphite-Epoxy 64 is calculated from the finite element analysis by changing the design variables to satisfy the imposed design criterion. Several combinations of the design variables were determined through simulation. However, in o of lamina angle was chosen for simplicity of this study, the 0 fabrication. Magnamite graphite prepreg tape, type AS4/3501-5A, manufactured by Hercules Inc., was employed to fabricate the composite arm. The arm was fabricated and cured by Auto Air Corporation, Lansing, Michigan, using the hand lay-up technique, vacuuming and autoclaving based on the curing cycle provided from Hercules Inc. as shown in Figure 3.5. The cure cycle depends upon the resin used, and the material properties are dependent Table 3.2 Material and Geometrical Properties of the Composite Arm Properties Specifications Beam Material AS4/3501 - 5A Young's Modulus 115.0 GPa Beam Length . 0.812 m Beam Width 19.02 mm Beam Thickness 2.5 mm Beam Mass 0.06 kg 0 Fiber Orientation 0 Ply Thickness 0.12 - 0.13 mm Number of Ply 20 65 Figure 3.6 Photograph of the Experimental Composite Arm annn the curing process provided from the manufacturer. The cure process also suggests a post-cure to allow for final cross linking of the polymer chains. The post-cure process has been done at 300°F for three hours with ambient pressure and no vacuum before the cutting laminate with the constrained geometry. Table 3.2 presents the geometrical and material properties of the composite arm, and Figure 3.6 shows the photograph of the composite arm mounted on the rigid hub. In order to experimentally evaluate the imposed design criterion, the load-deflection characteristics of the both aluminum and composite arms were obtained through the static testing. One end of the arm was clamped on the working table and monotonically increasing load was imposed on the other end of the arm, and the corresponding end-point (tip) deflection was measured subsequently. 66 4o . AHA ALUMINUM , o—e COMPOSITE E . g 30~ Z .J O .: . 0 LIJ . _J [:3 20— D d [— g 2 O ,L . C', 10— z 4 LL] 0 r T I I I T I I I f T 1 I I I T I j I 0 1o 20 3o 40 50 LOAD(grom) Figure 3.7 Load-Deflection Curve for the Aluminum and Composite Arm Figure 3.7 presents the load-deflection curve of the both arms. It is seen from this figure that the flexural rigidity (bending stiffness) of the composite arm is little smaller than that of the aluminum arm. The flexural rigidity of the composite arm is approximately 92 percent of that of the aluminum arm. It is recognized that the error is attributed to the manufacturing process such as curing cycle. From this figure, the Young's modulus of the aluminum was distilled by 64.3 GPa, and the Young's modulus in the fiber direction of the composite material was distilled by 115.0 GPa. 3.2.2 DC Motor and Servo-Drive The DC servo-motor shown in Figure 3.1 is a permanent magnet field, step-free torque motor, model 1338B-AV20, manufactured by Allen-Bradley. 67 The motor has some standard features such as the capability of 100 percent rated torque output at stall, the minimization of torque ripple and motor cogging by employing a skewed armature. The peak stall torque of this motor is 19.7 N-m and the peak stall current is 50 amperes at the temperature of 40°C. It is noted that the fundamental flexural bending and torsional natural frequencies of the motor shaft are 4.9 kHz and 2.1 kHz respectively, which are of course far beyond the natural frequencies of the aluminum and composite arm. The servo-motor is driven by a servo- drive, model l388-series B, manufactured by Allen-Bradley. This servo- drive is a pulse width modulated (PWM) driving system which consists of a 150V DC power supply, a transistorized power amplifier, and isolated current sensing and performance testing points such as current feedback signal, current error signal to the power amplifier, and tachometer signal. 3.2.3 Servo-Positioning Assembly The servo-positioning assembly consists of a servo-controller module, model l773-M3 and a servo-expander module, model l77l-ES, both manufactured by Allen-Bradley. The commanded position data from the microprocessor is sent to the servo-controller module, and the update status from the servo-expander is sent to the microprocessor. The servo- controller sends axis motion commands to the servo-expander which commands the axis motion by sending an analog signal to the servo-drive through D/A converter which has a 12 bits resolution. The servo sampling period is 2.4 milliseconds and the servo output voltage is ilOV DC maximum . 68 3.2.4 Sensors 3.2.4.1 Angular Position Sensor An optical incremental encoder was mounted on the shaft of motor using a flexible encoder coupling adapter,model l326DP-MOD-Ml-C2, nunnifactured by Allen—Bradley, in order to measure hub angular position. The incremental digital encoder, model 845N-SJH-N4—CRY2, manufactured by Allen-Bradley, has 2500 lines and provides feedback signals that indicate the magnitude and direction of any change of axis position. Also the phase relationship of the encoder output signals allows the decoding circuit to count either 1,2 of 4 feedback pulses for each line of the encoder, hence increasing the resolution of the feedback signal. 3.2.4.2 Angular velocity Sensor The tachometer mounted on the shaft of the motor was used to measure hub angular velocity feedback signal. The tachometer has voltage gradient of 0.3 vwsec/rad and 2.0 percent of ripple. The tachometer has 20 turn potentiometer that scales the tachometer feedback signal, and the gain of potentiometer is adjustable through the pot which has ten graduated increments. In order to verify the voltage gradient provided from the manufacturer, the joint was rotated at constant velocity, and the output voltage of the tachometer was measured, then the corresponding velocity was computed. Figure 3.8 presents the relationship between the velocity and the measured output voltage. It is distilled from this figure that the voltage gradient is 0.298 v-sec/rad which is almost same as the value provided by the manufacturer. 69 15 ’o‘ .. Q) U) \ ‘55 10a E U o -4 _l LIJ > § 5— D 0 Z < .1 0 T I r r I T r I I O 1 2 3 4 5 TACHOMETER OUTPUT(volt) Figure 3.8 Relationship between the Velocity and the Measured Voltage of the Tachometer 3.2.4.3 Tip Deflection Sensor The strain gages, type WD-DY-2SOBG-350, manufactured by Micromeasurement Group Inc., were mounted on the root of the arm in order to measure tip deflection, even though the signal was not used as a feedback signal in this study. The half-bridge circuit was arranged to measure the strain. The measured strain was conditioned and amplified through Wheatstone bridge conditioner/amplifier, model 2100, manufactured by Micromeasurement Group Inc. In order to obtain calibration factor between the amplified voltage (strain) and tip deflection, the series of load was imposed monotonically on the free end of arm and the corresponding deflection and amplified voltage were measured. Figure 3.9 70 shows the test results of this calibration for the both aluminum and composite arms. 3 . 2 . 5 Microprocessor The microprocessor, model PLC-2/l6, manufactured by Allen-Bradley was employed for the implementation of the output feedback controller associated with two colocated output sensors - angular position and velocity.ffluaudcroprocessor has features such as 3K words of memory size, computation capability up to 12 digits and conversion binary to binary coded decimal (BCD) which is a numbering system used to express individual decimal digits (0-9) in four bit binary notation. The personal 40 , A--A ALUMlNUM _ o—e COMPOSITE E . \E’ 30- z . Q 5 ~ Lu .. .1 L5 20—« O u l- ; .4 o l -4 CL 10— z LIJ '1 O I I I I I I F 0.0 0.5 1.0 1.5 2.0 AMPLIFIER 0UTPUT(voIt) Figure 3.9 Relationship between the Amplified Voltage and the Tip (End-Point) Deflection 71 computer, model ZVC-lS9-2, manufactured by Zenith Data Systems, was used as a program operating system. The program was written as relay-like language developed by ICOM Inc. [152]. 3.3 Identification of System.Hbde1 Parameters In this study, the rigid-body mode and the first two flexible modes are retained as the primary modes. Hence, the equation (2.35) becomes ~=K~+E~ xp pxp pu (3.3) ”=6; yP P P where P 0 0 - l 0 K- p -w1 ~2§1w1 0 O 0 0 0 l _ 0 0 -w2 -2§2w2‘ - 1 d¢1(o> d¢2(o) T Bp'T ° 1 0 'Ti" 0 —E_ T ' d¢ (o) d¢ (o) ' __1___ __2___ 1 0 dx 0 dx 0 Op ‘ d¢1 d¢2l XI + an El “transfer function from the input torque to angular position is given by YI(S) A _1 = 60(5) = c1 (51 - AP ) B (4.2) 5(5) p where s is the Laplace variable. The control torque 3(5) is given by u(S) - Kp( r(5) - y1(S)) - Kvy2(8) (4.3) Assuming the angular velocity sensor output is the exact derivative of the angular position sensor, the equation (4.3) becomes u(S) - KP ( r(S) - y1(3)) - sKvy1(S) (4.4) Now, the closed-loop transfer function of the system is obtained by §1(s) é o ( ) _ KD 60(5) c2 3 1 + Kv(s + Kp/Kv) 60(5) (4.5) Ecs) Therefore, the closed-loop poles of the system are roots of the following characteristic equation given by 1 + xv (s + KP / xv) 60(5) - o (4.6) The shape of the loot-locus of the closed-loop system, that is, the closed-loop system performance depends on the location of the compensator zero, defined by zcoé -Kp/KV. With the certain compensator zero, the rigid-body closed-loop poles may move from the origin of the s-plane (two open-loop poles are zeros) toward the open-loop zeros of the system which are characterized by cantilevered-modes, or toward the real axis of the s-plane resulting in sluggish motion. IMAGINARY IMAGINARY 91 390 ~ A open—loop poles - o open—loop zeroes ALUMINUM ARM " I compensator zero 260— - A . l30~ “ A 0 —I s A —130~ o . A —260 I I I r I I I r l I I I I l I I I —4.0 —3.0 —2.0 —-1.0 0.0 REAL 2590 ~ A open—loop poles - o Open—loop zeroes COMPOSITE ARM ‘ - compensator zero 260- I ‘ - 130- - A .. O 0 I - o .. A —130- - A . “-260 I T I I I I I I l I I I m I I I I -4.0 —3.0 —2.0 —1.0 0.0 REAL Figure 4.2 Open-Loop Pole-Zero Location with the Compensator Zero 2 - -1 co 92 Figure 4.2 shows the open-loop pole-zero location with the compensator zero Zoo: -1 for the experimental aluminum and composite arm whose material and geometrical properties are presented in Tables 3.1 and 3.2 respectively. It is seen from this figure that all non-zero poles are located left to the given compensator zero. It is noted that open-loop zeros and poles physically represent the clamped-mode and pinned-mode natural frequencies of the arm respectively. Figure 4.3 presents the root-locus of the closed-loop system obtained by increasing the feedback gain RV (or Kp, since the ratio is fixed) with the compensator zero Zco=-l. Since the root-locus is symmetric with respect to real axi5,only 250 : --- ALUMINUM - -— COMPOSITE < : 2004 W : mm mm ’3 q , 4 ' ‘ poms“): 901.239) 3 (D 150'“ “---0 Q _ 0.0, 0 0 0.0, 0.0 'o . 0.0, 0.0 0.0, 0.0 o . -1.288, 57 47 4.519. 66.07 L q -1.288,-57 47 -1.519,-66.07 E.” 100- —2.363, 165.23 -3.256, 223.03 a: _ -2.363,-165.23 -3.256,-223.03 2 : mosh): zsnosu» o . —————— 1- -0.525, 23.4.3 -0.777 32 79 g 50: ’ -0.525,-23.43 —o.777,-32.79 _ q ', -2.102, 147.01 -3.036, 207.95 . . __________ 4 -2.102,-147.01 -3.036.-207.95 0.: ..- fl IIOOMPENSATOR mam-1.0.0.0) “'50 I I I I l I I I I I -120 -80 -40 O REAL(rad/sec) Figure 4.3 Root-Locus with the Compensator Zero Z - -1 co 93 upper plane of the locus is presented. It is interesting to observe that the open-loop poles located at the origin (rigid-body mode) start to move toward to have imaginary part of the poles, and quickly move toward to have real part only as the feedback gain increases for the both aluminum and composite arms. The one rigid-body pole goes to the negative infinite on the real axis, while the other rigid-body pole ends at the given compensator zero as the feedback gain goes to the infinity. This phenomenon will results in the sluggish response of the closed-loop system without generating large overshoot to the step commanded input. As the feedback gain increases, the other open-loop poles move toward close open-loop zeros which correspond to the finite two cantilevered-mode natural frequencies. The same type of root-locus shape is obtained with the compensator zero Zco=-2 which is located between two non-zero poles on the real axis. Now, for the same dynamic model, consider the pole-zero plot on the s-plane with the compensator zero Zco--15 which is positioned far left to the all open-loop poles and zeros as shown in Figure 4.4. Figure 4.5 presents the root-locus as a function of the feedback gain Kv with the compensator zero Zco--15. Two poles located at the origin (rigid-body mode) move toward close two open-loop zeros which correspond to the fundamental cantilevered-mode natural frequencies. Therefore, it may be said that the bandwidth of the rigid-body mode is limited by the fundamental cantilevered-mode natural frequency {ll-3.73 Hz and I'll-5.22 Hz for the aluminum and composite arm, respectively. The maximum bandwidth of the composite arm is larger than that of the aluminum arm by 28.5 percent. This larger bandwidth results in faster response to the 94 390 .1 A open—loop poles ~ 0 open—loop zeroes ALUMINUM ARM ‘ I compensator zero 260~ . A. g 130— < -I Z; . (D A < ... Q3 0 -~I : 1 —130-— , a A “—260 I I I I I I I I I I I I I I I I ~18 ~15 ~12 ~9 ~6 ~3 REAL 390 ~ A open~loop poles - o open—loop zeroes COMPOSITE ARM ‘ I compensator zero 260~ j ‘6 E 130— E u __ ~ A ((9 . 0 2 O a 3 O 1 A ~130~ 1 - A? “-260 I I I I I I I I I T I I I I I I ~18 ~15 ~12 ~9 ~6 ~3 REAL Figure 4.4 Open-Loop Pole Zero Location with the Compensator Zero zco- -15 95 250 : --- ALUMINUM . - COMPOSWE 200— W 1 mm EMS—HP 1T3 I; ‘ POLES(A): POLES(A): 0 150~ m _ O£.0£ mo,mo \\ . 0b,0b mo,mo 13 . -1.288, 57.47 -1.519, 64.07 L_ q -;.§g§,-57.47 -1.519,-64.07 \_, _d - . , 165.23 -3.256. 223.03 E 100_ -2.363,-165.23 -3.256,-223.03 :21:- : zanosu ): zaaosu ): (p - -0.525, 23.43 -o.777, 32.79 <2: 50— -o.525,-23.t.3 —o.777,-32.79 __ . -2.102, 147.01 -3.036, 207.95 . —2.102.-147.01 -3.036,-207.95 O_ __ _ 1 = I :COHPENSATOR ZERO(-15.0,0.0) -50 I I I I I I I I I I I I ~120 ~80 ~40 O REAL(rod/sec) Figure 4.5 Root-Locus with the Compensator Zero Z ~ -15 co commanded step input. From these characteristics of the closed-loop system depending upon the certain fixed compensator zero, appropriate output feedback gains have to be determined to evaluate the control performances for the aluminum and composite arm. It is seen from Figures 4.3 and 4.5 that any choice of the feedback gains KP and Kv assure the stability of the closed-loop system. This can be interpreted that the control law is equivalent to a mechanism consisting of the torsional spring and dashpot mounted at the pinned end of the arm. Therefore this dashpot always provides a damping to the system ensuring the stability. However, this guaranteed stability is no longer valid in practical system, since unmodelled dynamics of the 96 flexible link, actuators and sensors, delays caused by measuring and sampling feedback signal and the inherent limitation of servo-drive are always present. This instability will be investigated through the experimental test. Another useful way to determine feedback gains which guarantee theoretically the stability of the system is to use the Lyapunov method. Consider a system described as x + fii + 6x = 0 (4.7) ~ According to the Lyapunov theorem for the stability, if both matrices P and a are symmetric and positive-definite, the system (4.7) is asymptotically stable. This theorem can be easily proved by choosing a Lyapunov function candidate V by l v = 5 0 and V < 0 for all x #0. The dynamic equation (3.3) can be reformulated in the form of the equation (4.7) using the modal coordinate given in the equation (2.34). Hence, the equation (3.3) becomes 0 + 25 + on - Bpfi (4.10) 97 2 2 where 0 (= diag (O,w1,w2)) is the modal frequency matrix, 2 (= diag (O, 2§1w1,2§2w2)) is the modal damping matrix and Bp is the control influence matrix which is given by A 1 d¢1(0) d¢2(0) T p I dx dx (4.11) Also the output matrix given by equation (4.1) can be rewritten as y1=clfl A . (4.12) y2=C2W where c — c - I B T (4.13) Since the system described by equation (4.10) is controllable, the PD output feedback controller is given by u - -pr1 - KVy2 (4.14) where KP is the constant position feedback gain and RV is the constant velocity feedback gain. The substitution of the equations (4.12) and (4.14) into the equation (4.10) yields 0 + (z + BpKv02)n + (a + BPKPC1)n = 0 (4.15) It is noted that the matrices Z and 0 are symmetric and positive semi- definite due to the rigid-body mode which is characterized by zero rmtural frequencies. Thus, the feedback gains KP and Kv in equation A A (4.15) should be carefully chosen such that the matrices (Z + BpKvCZ) and 98 A A (O + BprCl) are symmetric and positive definite. Then the resulting closed-loop system (4.15) will be asymptotically stable. The following proposithniis stated without a proof before determining the feedback gains. Proposition 4.1 : Let P and Q be any two symmetric matrices with the same dimension. If P is the positive definite and Q is the positive semi- definite, then (P+Q) is the positive definite. Now, under the condition stated in the Proposition 4.1, the feedback gains KP and RV can be chosen such that the closed-loop system be asymptotically stable. For example, let A (4.16) where P and Q are arbitrary symmetric and positive definite matrices. Then the feedback gains KP and Kv are obtained as follows. A T A -1 A T A AT A AT -1 K = B B B P c c c p 1 p pJ p 1 [ 1 1] (4.17) A T A -1 A T A AT A AT -1 RV - [Bp Bp] Bp Q 02 [02 c2] It is noted that the minimum norm solution(1east effort solution) and least error solution method [28] were employed to obtain the feedback gains. It is also noted that both matrices Z and 0 in equation (4.15) remain as symmetric and positive semi-definite, even though the model 99 parameters such as natural frequencies and damping ratios are perturbed in a certain threshold. Therefore, the perturbed system will be still asymptotically stable so long as the output feedback gains are chosen based on the above manner. Thus, it can be stated that the proposed output feedback controller assures the robustness of the system stability. The feedback gains determined from the afbrementioned two approaches will be used in the following experimental implementation. 4.3 Implementation Results and Discussions Figure 4.6 presents the experimental set-up for the output measurement. All experimental results presented in this section were obtained based on the following measurement procedures. The angular position feedback velocity feedback r 1 j strain gage servo positioning assembly " servo drive " motor ' flexible arm] commanded input torque velocity position status microprocessor analog filter Str232 Egggiéigger 1 1 A/D data D/A c~ acquisition - analog filter system Figure 4.6 Experimental Set-Up for Output Measurements 3:“ out 100 + channel 1 F——0 + channel 2 IO‘ out V + {IL _J_ channel 1 1"} in _01 V+ channel 2 +0— 1 in ~ 0— Circuit Diagram of Two Channel Operational Figure 4.7 Amplifier position data stored in the microprocessor as digital values were microprocessor to the data acquisition system transferred from the through D/A converter, model l771-OFE1, manufactured by Allen-Bradley. The D/A converter has output range of O to 10 volts DC. Since the data stored in the microprocessor was isolated from the main computer network which is employed for the data postprocessing and plotting the results, this transfer procedure was necessary. The employed data acquisition system, model MicroVMS II/GPX 21, manufactured by Digital Equipment Corporation, has 15 channels of A/D converter. The data sampling rate of 150 Hz, which is far beyond the considered natural frequencies of the arm, was employed for the measurement of all output signals except the second mode exciting response where 800 Hz of sampling rate was used. The input torque and 101 angular velocity from the servo-drive were amplified through the operational amplifier, since the actual output signals were relatively small resulting in significant effect of the computer noise of the data acquisition system to the actual output signals. Two channel operational amplifier producing the amplification factor of 100 was assembled as shown in Figure 4.7. The amplified signals then were fed to a dual analog filter, model 432, manufactured by Wavetek Inc., in order to eliminate high frequency components such as noises. The 100 Hz of cutoff frequency was employed before the filtered analog signals were digitized by A/D data acquisition system. The output signal from the strain gage mounted on the root of the flexible arm was conditioned and amplified through Wheatstone bridge amplifier, model 2100, manufactured by Measurement Groups Inc.,before the conditioned and amplified signal was fed to the low-pass analog filter with the cutoff frequency of 100 Hz. 4.3.1 Step Responses without Disturbances For the first comparison of the closed-loop performances between the aluminum and composite arm, the small compensator zero Zco--2 was chosen. With this compensator zero, the feedback gains KI) and RV were decided such a manner that the damping ratio of the closed-loop poles of the rigid-body mode is equal to 1 (critical damping). This means that the closed-loop poles of the rigid-body mode has no imaginary part for the first time as the feedback gain increases. From the root-locus analysis, the feedback gains Kp- 0.534 N-m, Kv- 0.267 N-m-sec and Kp- 0.292 N-m, Kv - 0.146 N-m-sec were obtained for the aluminum and composite arm 102 respectively. Table 4.1 the presents closed-loop poles of the system with these feedback gains. Table 4.1 Closed-Loop Poles with the Compensator Zero Zco= -2 Aluminum Arm (Kp — 0.534) Composite Arm (Kp - 0.292) Real Imaginary Real Imaginary -3.820 0 -3.638 0 -4.425 0 -4.538 O -l9.736 52.19 -lO.948 62.26 -l9.736 -52.19 -10.948 -62.26 -6.119 164.03 -4.279 222.89 -6.119 -l64.03 -4.279 -222.89 Figure 4.8 shows the measured step responses to a commanded step angular position of 7.2 degree. The hub angular position reaches its commanded position within 2.0 second for the both aluminum and composite arms. No overshoot in the angular position response is observed as expected by the employed feedback gains. The steady-state error between the commanded and actual position was 0.1 degree for the aluminum arm and 0.08 degree for the composite arm. This steady-state error is attributed to the joint friction and can be eliminated by applying the integral control action or manual calibration method. The significant different response between two arms is observed from the input torque response. The applied input torque for the composite arm is only 53 percent of the torque required for the aluminum arm. This amount of reduction of energy consumption will provide tremendous benefit of the economic area on the 103 .Eu< ouwmomEoo on» How o¢H.OI>x .NNN.0I a x was .5... 5:552 2t > a new mwN.0n.Vw fqmm.ol x .pwoaxmm udosuaz momcoamom moum pousmmoz w.¢ ouswwm Aommvmic. Aoomvmzz. n v n N F o m w m m __. o r s p t _ » _ L w p . — . b p p h . No.0' L h *1? p P b p h h b b p P ~ r p L Flu. .NI .0 n l l . O . u 8 D H n .. B ..m. s .N .H 9 - 7w. 1¢O.O . IN mtmonmzoo l m- Prion—200 II I 23233< 23233.2 8.0 n Aoomvmzfi noomvwzz. n v m N a o 0 ¢ n N F o p s p F » r P w . h t _ _ _ L _ p » .VO.OI r L p p _ p h b L p P _ — P _ _ p H . p p O w IN 9 . n _I w . no A 1¢ 3 1 1 O . D . m -m J D _l p I // It ................. J. S I 9 bx um 5509200 I utmonzoo II. . 23233.2 w 232.22.: de (UJUJ)NO|1331:130 dLL (fiapmomsoa awnONv 104 .Ep< ouwmomaou ago now qu.OI>M .NmN.ol a x can .52 .5552 05 > . 9 HOW NwN . OI M *wmm .OI. M .UNOHuhwm U50£UH3 mOwGOQmwm Qwum UmUQHUEHm 0 ..Q mufiwwh Acumvmzz. Aommvuzc. n w m m o m e n N F o h p . s p L p p L . _ p p . b I» p » No.0! » . p b _ . h L p . h h . . p ~ _ p b p _ . Fl! .N. .d n . 1 u. m .. no 0 0 H . n H u \G/ .. u. I .N ...” . we w m. $40.0 IN memouzoo I m. memoazoo I. . 23233.2. 1. 232.2312 00.0 m. Aoomvmfiz. Aoomvwic. w v m w n N F o p p c h t p p p b p . _ p — p p h p _ P p b L .— » ~ - p O W IN 9 I n I w . ud A 1+ 3 .. 7 O I D I m um W 1 D. .. / a I 0 ..m mtwonioo 1| utmoazoo ll . 23233.2 fl 232.23.? . ¢N.o o— (ww)Nou:331330 dLL (fiapmoulsoa awnONv 105 pmoahmm usozuws Eu< ESCHESH< man no noncommom doom mooq-como poumasaam oa.¢ ouswwm noonvwzc. . Aoomvmzrr 3.. n N F o a a. n N I o _ . h b b — . s p u b h h p » 00.0.] p b b — p b p p _ n b p . h . . P p _ b p - . CIVI- . f I 5.0: I 10d: . w 1 V n I r l . . I. . O . 100 O W 1 IO 0 1- n U . .B . . .N . .. (wk .. I 6.0 no u 1 . r . no.0 or. “003sz Aoomvwzc. #4 n N F O n .v — — b p _ L b . p r p _ . b n f OF.OI . . p _ . u . p m . p p . W p . p . W p p p t O 0.0 wood v w . . m . o N ... m I H.02. . . . no ... . m . T O .. I D I [one m was n W . . / I iii n . 3 III . Won 0 bx Tod 3.0 . 0. o (ww)Nou02maa dll (bep)N0Lusocl awnow 106 real production line. It is also observed from the tip deflection response that the magnitude of the vibration amplitude of the aluminum arm was reduced up to 30 percent by employing the composite arm. The second mode effect to the total tip deflection is very small for both aluminum and composite arms. The response decayed to zero within 2 second for the aluminum arm and 1.8 second for the composite arm respectively. It is noted from the angular velocity response that the ratio of the maximum angular velocity to the first pinned-mode natural frequency, i.e., 0 is 0.0033 and 0.0029 for the aluminum and composite arm max /wl’ respectively, hence eliminating the possibility of the resonance phenomenon due to the input torque. Figure 4.9 shows the simulated step responses corresponding to the measured response presented in Figure 4.8. The agreement between two results is excellent, even though there exit small mismatch in the response of tip tip deflection. For the comparison of the step responses between the open-loop and closed-loop system, the open-loop step input torque was introduced to the model. Figure 4.10 presents the simulated open-loop responses of the aluminum arm for the commanded angular position of 7.2 degree. The magnitude and period of the step input torque was determined to have almost same delay time (0.42 second) as the case of the closed-loop system. It is clear from this figure that the angular position and velocity responses exhibit oscillations as expected, and the settling-time of the tip deflection is increased by 45 percent comparing with the closed-loop response shown in Figure 4.9. For the second comparison of the control performances between the aluminum and composite arm, the compensator zero Zco- - 5 was chosen and 107 same feedback gains were employed to the both aluminum and composite arms. Table 4.2 presents the closed-loop poles of the system with feedback gains of Kp- 1.0 N-m and KV- 0.2 N-m-sec for both arms. From this table, the damping ratio of rigid-body mode was obtained by 0.633 and 0.855 for the aluminum and composite arm respectively. Table 4.2 Closed-Loop Poles with the Compensator Zero zco- -5 Aluminum Arm (Kp = 1.0) Composite Arm (Kp - 1.0) Real Imaginary Real Imaginary -2.831 4.645 -5.526 5.159 -2.831 —4.645 -5.586 -5.159 -15.233 55.47 -14.523 61.18 -15.233 -55.47 -l4.523 -6l.l8 -5.307 164.58 -4.630 222.80 -5.307 -164.58 -4.630 -222.80 Figure 4.11 presents the measured step responses to a command step angular position of 7.2 degree with above feedback gains. From the angular position response, the delay time was obtained by 0.26 second for the aluminum arm and 0.21 second for the composite arm. It is also observed from the angular position response that the maximum overshoot of the aluminum arm is 1.02 degree, while that of the composite arm is 0.3 degree. From the input torque response, it is also seen that the actual input torque required for the composite arm during the commanded motion is only 60 percent of the input torque required for the aluminum arm. The tip deflection of the aluminum arm induced from the step motion was decreased up to 50 percent by employing the composite arm. From these 108 65.4. oufimomaoo Foam §GHESH< Soon (59P)N0lll50d awnsw > a no.“ N.ol M can 04.: x .pmoahmm 8.65:3 momcoamom doom pondmmoz :4» ouswfim Foams: Foams: ¢ n N F o m b h b I- d b P b b m h p p h m p p F W h b F p O o 1...F.F..FF..F_.Fr._ Fmodl onl .: Iood l N d n l l.— o no.0 w n B .N . M . 0.20 .........H ..m 10.». E6828 .1 memoazoo I. . 23233.2 23233.2 . 2.0 n... Aoomvmzfi Aoomvmzfi .v n N F o m F. m N F o F F F F F F F F F F F F F F P F F F F _ N.OII r F F F — L F F F — F F F F — F F F _ F F F F O I W 1 ..... o 1n n W . no A I 3 m .m D m J IIII<4II ........ I m . D. V a 1m P . memoazoo I. . 85828 I . 23223... F 232.23% 0.0 (UJUJ)N011331:13C1 dli 109 > How N.0I x tam olfil Foomvmzc. a .msu< ouamomEoo paw ESCHESH< fiuon utmoaioo ..l F23z§33< F0352: utmoaioo .II 23253.? (w'N)3noam inle (Gas/PDJ)AIIOO'13A awnoNv M .pmoammm unozufia momcoamom moum nouMHSEHm NH.¢ ouswfim F832: I. N F o b h P P — b F p b h F L p h b h p b F p Oonl u d 0 I L 3 3 1 . 3 u o 0 u . O ... . w, ..m. w ”.... m. .....m -o.n 85828 .I . 232.23... OF. 83sz m F. N F o p u r h k u n h P u H h h n H h P n b O In tttttttt I In (53P)NOIllSOd WHONV wtmoazoo II 232.23.? 110 results, it may be summarized that with same feedback gains, the system performances of the composite arm have been remarkably enhanced over those of the aluminum arm, exhibiting the smaller delay time by 20 percent, the smaller maximum overshoot by 70 percent, the less energy consumption by 40 percent, and the smaller tip deflection up to 50 percent. It is observed from the tip deflection response that the tip deflection is not zero even after the angular position damped out to reach its commanded position. This phenomenon arises from the increased hysteretic joint friction due to the relatively high employed position gain. Figure 4.12 shows the simulated step responses corresponding to the measured results shown in Figure 4.11. It is observed from Figures 4.11 and 4.12 that two corresponding responses agree well except for the first 0.2 second of the motion. Next comparison of the control performances was undertaken with the payload of 0.061 kg which is 46.2 percent and 101.6 percent of the mass of the aluminum and composite arm respectively. The compensator zero Zoo- -1 was chosen with the feedback gains of Kp-0.48 N-m for the aluminum arm and Kp-0.24 N-m for the composite arm based on the same manner as in the first case. Table 4.3 presents the closed-loop poles of the system with these feedback gains, and from this table, it can be seen that the closed-loop poles of the rigid-body mode has no imaginary part for both arms, hence expecting the sluggish motion and no overshoot. Figure 4.13 presents the measured step responses to a commanded angular position of 7.2 degree for both arms with the payload of 0.061 kg. As expected, it took 3.5 second to reach its commanded position without occurring overshoot. The required energy of the composite arm for 111 this movement was only 48 percent that of the aluminum arm, and tip deflection of the aluminum arm was reduced by up to 60 percent by employing the composite arm. With smaller energy consumption, the more accurate end-position was obtained. It is expected that with the same ratio of the payload to the arm mass, the better control performances of the composite arm can be achieved with appropriate feedback gains. Figure 4.14 shows the simulated step responses corresponding to the measured responses presented in Figure 4.13. Again the correlation between the measured and simulated responses is impressive. The employment of the Lyapunov method is also another useful technique to determine the feedback gains which guarantee asymptotical stability of the closed-loop system as discussed earlier. Herein, the A A following symmetric and positive definite matrix P and Q in the equation (4.16) were chosen for the both aluminum and composite arms. Table 4.3 Closed-Loop Poles with the Compensator Zero Zco- -1 Aluminum Arm (Kp - 0.48) Composite Arm (Kp - 0.24) Real Imaginary Real Imaginary -1.234 0 -1.679 0 -6.418 0 -2.578 0 -36.197 28.93 -20.125 44.14 -36.197 -28.93 —20.125 -44.14 -9.656 126.82 -S.l7l 166.96 -9.656 -126.82 -5.l7l —166.96 .Eu< ouHmOQEoo o5» mom qN.0I>¥I a F as .52 5552 23 now wqdl MI M .mx .306 mo pmoahmm 4.53 monsoamom moum poHSmmoz max» ouswwm > a Foams: n F4 n N F O L F F r — P F L F p F F F F b H No.0! . 18.0 mtmOQIOU I r 232:23-2 2 no.0 1 1 F0852: m .v m. N F O F F H F F F F F F F F F F P ”0.0:! utmoazoo ...l 232.233< (w'N)3n0301 111le (ass/poJ)woo13/\ awnONv €93sz utmoaioo .1. 23233-2 FommvaF utmoazoo I 232.233< NF (59P)NOLUSOd awnst (ww)NOL103“L-130 dli 113 > a a .Su< ouHmomfioo ozu How sN.0I>MI M can .EH< ESGNESH< Gnu Foo ms.ou an x .mx Foo.o mo umoFFmF cuss mmmcoammm qum emuquaFm FF.3 mustF “93sz utmoa2oo .II 232.233< 1111:3‘. Iwod Foams: mod m._._moa2oo II 23223.? €3.sz (w'N)3noam 1(1le utmoa2oo II 232.23.? .. AoomvaF :0 Hr (ass/pmmIoma/x mnsw (bapmomsoa awnorw utm0d2oo I v 23223.? NF (ww)~ouoauaa dll . 30 o P=-[0 80 o o . so 0 Q= 0 so 0 0 With above matrices, system input matrix Bp 0 50 114 A (4.18) A , and the output matrices C1 and C without the payload, the feedback gains of KP - 0.406 N-m, K‘,- 2 0.254 N-m-sec for the aluminum arm, and K1) =0.403 N-m, Kv = 0.252 N-m-sec for the composite arm were obtained from the equation (4.17). Table 4.4 presents the closed-loop poles of the system with above feedback gains. As expected, all closedeloop poles fOr both arms have negative real part assuring the stability of the system. Table 4.4 Closed-Loop Poles with Kp-0.406, Kv-0.254 for the Aluminum Arm, and Kp- 0.403, Kv- 0.252 for the Composite Arm Aluminum Arm Composite Arm Real Imaginary Real Imaginary -2.303 0 -l.843 0 -5.560 0 -13.430 0 -18.806 52.60 -17.268 58.17 -18.806 -52.60 -17.268 -58.17 -5.960 164.13 -4.937 222.64 ~5.960 ~164.13 -4.937 -222.64 FF .5... 83858 an F8 NmN.o..>VF .83.? F 28 .E< 53.52 of you smN.oI>M .moqélaM .pmoammm ”50:33 noncommoM noum ooh—ammo: was 93w: 115 Foams: Foams: v n N F o n W h h P P m P h h. P m h P F b W - h b .F O o F F F — F F F F — F F F F P F F F F p F L F F NOIOI L 0 Fl - m d {7.11 n I. l 0 no 0 n \B .N H /w\ H memoazoo I . memoazoo .l .u I 23223.? . 232.23.? mod o.N F0862: Foommz: F. a F. n N F o r F F F _ L1 F L F F F — F F F P I— L F F r— L F F F _ F L F F O I. v T N . IN 9 I V n N m . o I n m .I. m w I m I S u I w \z :0 F . m V «41114 414.4... .....I. I w- e F .7wa (G rm I N 5501200 I I MEmOdzou I I 232_233< I 232.23% I 0N.o (ww)Nouo:rL-Iao dll a .eS «389:8 23 F8 «35:3. .23.? FF 28 6.2 5552 2t > . a now «mm.ol M 00.30! M .vmoahmm usonuHB mmmfioammm mmum Umudadawm 0H3 05%.: 116 AommeZ; Aoomvmfiz. Lr F F F W L F F F m F W W O m ¢ n N F O by P h b h h b h H .— F NOIOI h by n — n n h b —, h by L. P — F P b b P h h b h OIF| 1 r M d - - n l l O HO nu n \B .N .H W m mtwomfioo III F. utmOQIOO ll . 23232312 j 232—23.? w . 00.0 ON Aommvmzc. Aommvmzz. v 0 ¢ N N _. O P F F F — p P h P P h b b P P F b P h h \P h h b h b b h L h b b \P O w 1m 9 I n, I W i no A [v 3 .. «I O 1 3 w T ) IO m 1 p ...... r / . S 1 a («V 1w utmomioo ..lu .. wtmOaEOU II F. 232.23.? v 23233.? r ON.O OF (ww)Nou331330 le. (59P)NOILISOd awnow 117 Figure 4.15 presents the measured step responses to a commanded angular position of 7.2 degree with above feedback gains. It is observed form the angular position response that there is no overshoot and it takes about 2.5 second to reach its commanded position for the both aluminum and composite arms. From the input torque response, it is distilled that the control input torque required for the composite arm is approximately 60 percent of that of the aluminum arm. It is also seen from the tip deflection response that the maximum tip deflection of the aluminum arm was reduced up to 35 percent by employing the composite arm. A It should be noted that the different choice of the matrices P and £2 will produce different control performances of the system. Figure 4.16 shows the simulated step responses corresponding to the measured results shown in Figure 4.15. The agreement between two results is again excellent. It is noted from these results that no controller/observer spillover was occurred. This is attributed to the relatively high unmodelled frequencies (truncation effect is very small), large inherent damping of the system, and no state observer in the control system (no state estimation error, hence fast servo-sampling time). 4.3.2 Step Responses with Disturbances It is generally known that the closed-loop system associated with the PD type of the output feedback controller is very sensitive to uncertainties such as parameter variation, internal and external disturbances. Herein, the investigation of this sensitiveness was accomplished by considering external inputs, which have forcing 118 frequencies close to the cantilevered-mode natural frequencies of the arm, as the disturbances. This attempt is made here to introduce a single-link flexible manipulator fabricated from smart materials incorporating electro-rheological fluids which will be discussed in the subsequent chapter. The robustness of the PD controller to the disturbances is to be expected with the help of the distributed parameter actuator which is characterized by embedded electro-rheological fluids in the structure. With the introduced disturbance, the dynamic model given by the equation (3.3) becomes x-K§+§~+i§5t p pp p‘1 p” (4.19) ”-5; yP PP ;+ K L~~~~~~ yp p x- x+Bu+B d(t) 4» — PPPP P 3'2 Pi Kv-t Figure 4.17 Block-Diagram of the PD Output Feedback Control System with the Disturbance of d(t) 119 where d(t) is the imposed disturbance. Since the disturbance is applied to the servo-drive, the disturbance influence matrix is same as the input influence matrix F . Figure 4.17 presents a block-diagram of the closed- loop system with the PD output feedback controller subjected to the disturbance d(t). This distmubance is supplied from the function generator, model 275, manufactured by Wavetek, to servo-drive. Figure 4.18 shows the measured step responses of the aluminum arm without the payload to a commanded step angular position of 7.2 degree with the imposed disturbance of 0.08 sin flt N-m where f1= 3.25 Hz. The feedback gains of Kp- 0.534 and Kv= 0.267 used for Figure 4.8 were employed, and the disturbance was applied to the servo-drive from the starting of the motion. It is noted that the forcing frequency f1=3.25 Hz is very close to the clamped-mode fundamental natural frequency 01- 3.73 Hz in Table 3.5. As expected, the response of hub angular position and tip deflection are oscillated resulting near resonance phenomenon of the system. If the amplitude of the disturbance increase, the system response will be grown up resulting in the instability of the system. Figure 4.19 presents a comparison of step responses of the aluminum arm with and without the disturbance. It is interesting to observe from the angular position response that the reached angle after 2 second is 7.2 degree with the disturbance, while it is 7.1 degree without the disturbance. The reason is that the disturbance helps to overcome joint friction so that the arm can reach to the commanded position with the tradeoff the induced oscillation. Figure 4.20 shows a comparison of step responses between the ANGULAR POSITION(deg) TIP DEFLECTION(cm) 9i) 120 613* 31)— OI) TlME(sec) 1.2 llllLllelllelll l .0 0) 1 1 l J l l I —8 N O TlME(sec) Figure 4.18 Measured Step Responses of the Aluminum Arm without Payload. Kp- 0.534, Kv- 0.267, and d(t)-0.08 sin(3.25 Hz)t ANGULAR POSITION(deg) TIP DEFLECTION(cm) Figure 4.19 Comparison of the Measured Step Responses of the 121 12.0 with disturbance - —— without disturbance 8.0— 4.0-1 000 ' r* I 1 I T i I I I I T I T r I I I I 1' 0 1 2 3 TlME(sec) 1.5 : with disturbance . — without disturbance 1.0- o.5~ 5‘ 3. .55": J L J 0.. $‘fioo TlME(sec) Aluminum Arm with and without the Disturbance. xp-o.53a, Kv-O.267 and d(t)-0.08 sin(3.25 Hz)t 122 simulated rneasured TlME(sec) -- rneasured simulated 121) ’8 - O) 3 Z 9 t: (f) O 0. ‘5 D (D Z <( 1.5 1.0— A q E 0 V Z 9. 1.. 0 LL] ..l Ll. LLJ Q E i,— Figure 4.20 TlME(sec) Comparison of the Simulated and Measured Step Responses of the Aluminum Arm with the Disturbance. Kp-0.534, vao.267, and d(t)-0.08 sin(3.25 Hz)t 123 simulation and measurement with the imposed disturbance. The agreement between two results is satisfactory except for the beginning the motion. The reason of this mismatch at the beginning of the motion is that the motion is delayed due to the applied disturbance which induces the motion in the opposite direction to the commanded position. Figure 4.21 presents the measured step responses of the composite arm to a commanded angular position of 7.2 degree with feedback gains of Kp-0.29l6 N-m, Kv= 0.1458 N-m-sec and the disturbance of d(t) = 0.08 sin flt N-m where fl- 4.75 Hz, which is close to the fundamental clamped- mode natural frequency of 01= 5.22 Hz of the composite arm. As in the aluminum arm, near resonance phenomenon has been occurred even though the small oscillation amplitude is observed from the angular position response. Figure 4.22 shows a comparison of the measured step responses of the composite arm with and without the disturbance. Similar characteristics to the aluminum arm are observed. Figure 4.23 presents a comparison of the measured and simulated step response in the presence of the disturbances. The agreement of two results is excellent except first portion of the angular position response. Again, the delay of the response arises from the imposed disturbance which hinders the motion in opposite direction to the commanded step position. Figure 4.24 shows the measured step responses of the composite arm with the same feedback gains used for the Figure 4.21 and the disturbance d(t)=0.4 sin f2t where f2 = 34.5 Hz, which is very close to the second clamped-mode mode natural frequency 0 - 33.09 Hz of the composite arm. 2 ANGULAR POSITION(deg) TIP DEFLECTION(cm) 124 9.0 TlME(sec) -1.2 I l I r TlME(sec) Figure 4.21 Measured Step Responses of the Composite Arm without Payload. Kp-O.292, Kv-O.l46 and d(t)-0.08 sin(4.75 Hz)t ANGULAR POSITION(deg) TIP DEFLECTION(cm) Figure 4.22 125 121) with disturbance " — without disturbance 81)— 4a0‘ . 000 PI I I I I I I T I l I I I I —‘I I I fl I 0 1 2 3 4 TlME(sec) 1.5 2 with disturbance J — without disturbance 1 0—1 0.5— —o.54 ’ ‘ 1 —1.0 I I I I I I I I I I I I I I I I I I O 1 2 3 4 TlME(sec) Comparison of the Measured Step Responses of the Composite Arm with and without the Disturbance. xp-o.292, Kv-O.146, and d(t)-0.08 sin(4.75 Hz)t ANGULAR POSITION(deg) TIP DEFLECTlON(cm) 126 12X) simulated ‘ ——— rneasured 81)— 413— J: 000 I I I i I I I I I r I fl I 1 I I I I 0 ‘l 2 3 4 TlME(sec) 1.5 ‘ - simulated 1 ——- nneasured 1.0-3 1 .1 0.5-4 0.05 1 —o.5-3 4 1 —1.0 I 1 I fl 1 I I I I I I F fl I I i I fi I 0 1 2 3 4 TlME(sec) Figure 4.23 Comparison of the Simulated and Measured Step Responses of the Composite Arm with the Disturbance. Kp-0.292, KV-O.l46, and 3(:)-o.03 sin(4.75 Hz): ANGULAR POSITION(deg) HP DEFLECflON(cn0 127 9.0 6.0— 3.0— / 0.0 I I I l I I I l I I l I I I I l I O 1 2 3 4 TlME(sec) 1.2 0.6— O 0.: .,‘nfl11$IISI$Lw‘wtvxlI-v III £1H I:1"113W111MWINUIHI1M1attiimihy WU” —O.6- _1.2 l l I I I I I I 1 I I I I I I I I I I O 1 2 3 4 TlME(sec) Figure 4.24 Measured Step Responses of the Composite Arm without Payload. Kp- 0.292, Kv- 0.146 , and E(:)-o.a sin(34.5 Hz): 128 As observed from the step responses of the closed-loop system in the presence of the disturbance which has a forcing frequency close to the clamped-mode natural frequency of the arm, the responses are very sensitive to the imposed disturbance regardless of either the aluminum or the composite arm. Even though significant advantages of the control performances with the composite arm in the absence of disturbance were accomplished, there is a limit to improve the overall control system especially in the presence of uncertainties such as introduced disturbances. In order to circumvent this problem, a single-link flexible manipulator fabricated from smart materials incorporating electro- rheological fluids is proposed and investigated in the subsequent chapter. Since the natural frequencies of the arm is controllable in real time by employing an external voltage to the embedded electro-rheological fluid domains, the robustness of the PD output feedback control system is to be expected. 4.3.3 Instability of the System As discussed in Section 4.2, the closed-loop system with the PD type of the output feedback controller associated with two colocated angular position and velocity sensors is theoretically stable for any choice of feedback gains. However, this is not true in practice since the unmodeled dynamics such as actuator and sensors, the neglected higher flexible modes, the presence of time delay in the controller implementation, and the limited inherent servo-characteristics of the closed-loop system may cause the system instability. In this section, the instability of the closed-loop system caused by the limited inherent servo-characteristics 129 will be discussed through the experimental investigation. The servo-drive system employed in this experimental investigation has an in-position band which defines the maximum allowable position from the commanded angular position (in other words, the allowable positional accuracy) as shown in Figure 4.25. In Figure 4.25, the T1 and T2 represent reacting torques to force the actual position to the in-position band. The magnitude of these torques depend on the employed feedback gains, and the geometrical and material properties of the arm and the hub. The actual position in general remains within the in-position band with the small feedback gains (small T1 and T2 due the small gains, small spring and inertia force). However, as the feedback gains increase, the reacting Commanded Position \ / In-Position Band Hub Figure 4.25 In-Position Band of the Servo-Drive System 130 torques T1 and T2 become larger in order to force the actual position to the in—position band. Hence, with certain critical feedback gains, the reacting torques continuously increase resulting in the system instability. The experimental test was undertaken in a such a manner that the position feedback gain was increased until the instability of the system was occurred, while the velocity feedback gain was fixed. Since the instability may cause serious injury to persons or damage to adjacent hardwares of the system, the emergency switch has been activated after just 4 seconds of the commanded motion. Figure 4.26 presents the measured step responses of the aluminum arm to the commanded angular position of 7.2 degree with feedback gains of Kp- 6.2 N-m and Kv- 0.4 N-m-sec. With these feedback gains, the system became unstable resulting in the growing responses as time goes on. It is interesting to observe from the tip deflection response that the magnitude of oscillation does not grow up noticeably until 2 second, however after then, the response grows up rapidly. Figure 4.27 presents the simulated step responses corresponding to the measured responses given in Figure 4.26. As expected, the system is stable with these feedback gains. These feedback gains may be defined as maximum allowable feedback gains of the practical system featuring the aluminum arm. It is noted that the maximum allowable feedback gains have different combination of KP and Kv’ In other words, if the different velocity gain is chosen, the different maximum position feedback gain, which causes the instability of the closed-loop system, is determined. From these results, it may be recognized that the controller design of the system should be based on both analytical analysis and experimental ANGULAR POSITION(deg) TIP DEFLECTION(cm) l3l TlME(sec) I I ' ' 1 I 'T ' ' ' I O 1 2 3 4 TlME(sec) Figure 4.26 Measured Step Responses of the Aluminum Arm without Payload. Kp-6.2 and Kv-0.4 132 ANGULAR POSITION(deg) CD I O I I ' I I I I ‘ T l I U I I l I l I I O 1 2 3 4 TlME(sec) 3.0 A 1.5-1 E 8 . z 9 . B m 004 d w 4 o E -1 '- —1.5— _30 I T I I I I r I m I I I I I I I r I I O 1 2 3 TlME(sec) Figure 4.27 Simulated Step Responses of the Aluminum Arm without Payload. Kp-6.2, Kv-O.4 133 12 q - ALUMINUM —— COMPOSITE ”a. m -1 3 .' z . ' 0 —~ w~fi~ ;. {a e As— +~ ~ A E . (f) O 0. Cl: 5 I) (.9 Z < I I I f I I I I I l. I I I I 2 3 4 TlME(sec) 12 ~ ALUMINUM '-—- COMPOSWE 8.. E c1 8 . .. z _ '. 2 4 = '- ' : o . ,. . LIJ e'. . _l . . . u. ' 3 LLJ 4' D . '. : o. . -' : E i z. -4... o l —8 I I I I I I I I I I I I I I I I I I I O 1 2 3 4 TlME(sec) Figure 4.28 Measured Step Responses with Feedback Gains of Kp-6.2 and Kv-O.4 for both Aluminum and Composite Arms 134 ANGULAR POSITION(deg) CD 1 TlME(sec) TIP DEFLECTION(cm) O l TlME(sec) Figure 4.29 Measured Step Responses of the Composite Arm without Payload. Kp- 11.8, Ky? 0.4 TIP DEFLECTION(cm) 135 ANGULAR POSITION(deg) O) J O —I I I I —T n n I m I I I I I I I I I I O 1 2 3 4 TlME(sec) 3.0 1.5- 0.0-I v -1.5- —30 I I I I T I I I I I I I 1 I I If I I I O 1 2 3 4 TlME(sec) Figure 4.30 Simulated Step Responses of the Composite Arm without Payload. Kp- 11.8, Kv- 0.4 136 identification of the employed system. Figure 4.28 shows a comparison of the measured step responses of the aluminum and composite arm with feedback gains of Kp- 6.2 N-m and KV- 0.4 N-m-sec. With these feedback gains, the composite arm exhibits stable closed-loop responses of the system, hence shows significant improvement of the stability of the employed practical system. It is recognized that the reduced spring and inertial force of the composite arm is contributed to this result. The position feedback gain with the same fixed velocity feedback gain of Kv= 0.4 N-m-sec was increased until the instability of the composite arm was occurred. Figure 4.29 shows the measured step responses of the composite arm to a commanded angular position of 7.2 degree with the position gain of KP- 11.8 N-m. From the response of the tip deflection, it is also observed that the response quickly grows up after 2 second as the case of the aluminum arm. Figure 4.30 presents the simulated step responses corresponding to the measured responses shown in Figure 4.29. The closed-loop system is stable with these feedback gains as expected. From this experimental investigation, it may be summarized that the stability margin of this employed particular robotic system featuring the aluminum arm can be enhanced up to 45 percent by employing the composite arm. This is also another significant advantage of the composite arm over the aluminum arm. 4.4 Sumary of the Chapter The output feedback controller associated with two colocated sensors was designed utilizing the root-locus technique and the Lyapunov method, and experimentally implemented for the aluminum and the composite arm. 137 The comparative experimental investigation was undertaken to evaluate the characteristics of the control performances of the aluminum and composite arm. It has been demonstrated from the analytical and experimental investigation that the control performances of the overall closed-loop system for the both aluminum and composite arms were satisfactory, and the composite arm exhibited some superior control performances relative to the aluminum arm such as faster settling-time, smaller input torque, smaller tip deflection (better positional accuracy), auui superior increment of the stability of the practical system. These experimental results clearly indicate the significant payoffs associated with employing the composite materials featuring superior strength- and stiffness-to weight ratios to the commercial metals on the flexible manipulators. It has been also demonstrated that the aluminum and composite arm exhibited the inevitable sensitivity to the imposed. disturbance which has a forcing frequency close to the cantilevered-mode natural frequency of the arms. This attempt was undertaken here in order to introduce a single-link flexible manipulator fabricated from smart materials incorporating electro-rheological fluids for the purpose of enhancing the robustness of the output feedback controller to the disturbance. Subsequent chapter addresses this new class of flexible manipulator fabricated from smart materials. CHAPTERV HYBRID CONTROL FEATURING A DISTRIBUTED ELECTRO-RHEOIDGICAL FLUID ACTUATOR 5 . 1 Introduction As observed in the previous chapter, the closed-loop system associated with PD output feedback controller is very sensitive to the disturbance, which has a forcing frequency close to the clamped-mode natural frequency of the arm, regardless of the aluminum or the composite arm. In order to circumvent this problem, a single-link flexible manipulator fabricated from smart materials incorporating electro- rheological (ER) fluids (smart ER arm in short) will be investigated by proposing a hybrid control methodology. The hybrid controller consists of the output feedback controller employed for the aluminum and composite arm in the previous chapter, and a pseudo-state feedback controller characterized by a distributed ER fluid actuator [40]. A crucial ingredient in the development of an appropriate control algorithm for this class of smart materials is the development of an appropriate system (plant) model. Since the phenomenological behavior associated with the interaction of the structure and the embedded ER fluid is fairly complex, simplifying assumptions have been incorporated in the modeling process in order to furnish equations of motion which are analogous in mathematical structure to those encountered in beams fabricated from homogeneous monolithic materials. A phenomenological equation of motion, based on 138 139 these simplifying assumptions and experimental observations, admits a pseudo-state feedback controller to be implemented. The pseudo-control force (external input voltage) associated with the pseudo—state feedback controller accounts for the change in the global damping and stiffness properties of the structure due to the electrical potential applied to ER fluid domains. Since the natural frequencies and damping ratios of the flexible smart structures are tunable on real time by employing the external input voltage to the embedded ER fluid domain [41,56], the resonance phenomenon due to the imposed disturbance can be avoided by applying the pseudo- control forces. Hence, the robustness of the employed PD type output feedback controller is to be expected. In this chapter, a brief background on the ER fluids will be introduced prior to developing a phenomenological equation of motion of the proposed robotic system. Then a hybrid controller will be formulated based on the empirical model of the distributed ER fluid actuator. Finally, an experimental implementation of the proposed hybrid controller will be presented by investigating step responses of the system subjected to no disturbances and disturbances . 5 . 2 Background on Electra-Rheological Fluids Electra-Rheological (ER) fluids are typically suspensions of micron- sized hydrophilic particles suspended in the suitable hydrophobic carrier liquids, which undergo significant instantaneous reversible changes in material characteristics when subjected to electrostatic potentials. References [20,41,55,56,57,l36,137,138,140,151,160] provide a flavor of the research activities in this field. The most significant change in the 140 material characteristics of ER fluids is associated with the bulk viscosity of the suspension, which varies dramatically upon applying an electrical field to the fluid domain. This dramatic reversible change from a liquid-like state, as shown in Figure 5.1 (a) to a solid-like state as shown in Figure 5.1 (b) by the imposition of an appropriate electrical potential can be usefully exploited in vibration-control applications. I , . _ ’ w: _ a. ,» . ‘ ., ./ I. '. . (a) liquid state (b) solid state Figure 5.1 Electra-Rheological Fluid in Liquid and Solid State Figure 5.2 presents a mechanism showing that how this reversible change works in terms of interaction between charges placed on electrodes and those in the particles in the ER fluid. The charges in the fluid may 141 electrode ER particle electrode (a) without voltage (b) with voltage Figure 5.2 Interaction between Charges on Electrodes and Those in Electra-Rheological Fluid be either positive or negative, and free to move in the material in the absence of charge on the electrodes as shown in Figure 5.2 (a). When a voltage is applied to the electrodes, the charges in the fluid form a chain-like formulation in milliseconds as shown in Figure 5.2 (b). Figure 5.3 presents photomicrographs of an ER fluid subjected to electrical field intensities of 0 kV/mm and 2 kV/mm respectively. The photomicrographs were taken in the Biothermal Sciences Laboratory at Michigan State University using a Zeiss universal phase-contrast microscope with X40 magnification and a Chinon camera. The black regions at the top and bottom of the photographs are images of the electrodes employed to generate the electrical field in the ER fluid. Figure 5.3 (a) clearly illustrates the random structure of the suspension which imparts nominally isotropic global mechanical properties to the fluid mixture 142 when a potential difference is not generated between the electrodes. Figure 5.3 (b) clearly shows the truly dramatic change in the structure of the suspension upon developing a potential difference between the electrodes of magnitude 2 kV/mm. Under this condition, the particles in the suspension orientate themselves in relatively regular chain-like patterns to form a mixture with globally anisotropic mechanical properties. These columnar structures of particles restrict the fluid motion, thereby increasing the energy-dissipation or viscous characteristics of the suspension, in addition to increasing the stiffness characteristics of the fluid mixture. Thus, by imposing an (a) 0 kV/mm (b) 2 kV/mm Figure 5.3 Photomicrograph of Electro-Rheological Fluid 143 electric field upon an ER fluid, the stiffness and energy dissipation, or damping, characteristics of the electro-viscous suspension are changed. When the field returns to a state of zero potential, upon switching off the electrical energy supplied to the electrodes, the particles return to a state of random orientation in the carrier fluid as presented in Figure 5.3 (a). The voltages required to activate the phase-change in ER fluids are typically in the order of 4 kV per millimeter of fluid thickness, but since current densities are in the order 10 pA/cm2 or less, the total power required to trigger this phenomenon is quite low. Furthermore, the response of these ER fluids to the excitation voltage is almost instantaneous since the fluids can respond to voltage pulses with frequency up to approximately 11 kHz, which is, of course far beyond the requirements of most practical devices. ER fluid suspensions typically have viscosities of the order of 50 cP in the uncharged state, and some fluids are able to withstand a nonvanishing yield stress up to 60 kPa and shear rates up to 400/3. When an ER fluid suspension is stressed below this threshold, it behaves as a solid, but as the stresses are increased, then fluid flow occurs but the yield stress remains nominally constant. Since the characteristic aspects of ER fluids have been recognized by Winslow [160] in 1949 for the first time, numerous applications of ER fluids in engineering practice have been and being proposed and developed. Possible ER fluids devices and mechanisms proposed and developed so far are well documented in [128]. Typical applications include ER actuators, ER active damper systems, ER clutches, ER hydraulic valves, ER brakes, ER suspension systems, ER shock absorbers, rotary ER 144 viscometers, tunable ER isolators, and wide-band-high power ER vibrators. 5.3 Phenomenological Equation of Motion Phenomenologically, the global finite element equation governing the motion of the proposed smart ER robotic system has the following functional form [40,41,56,58] IMIII‘II + [C(V(t))]n‘n + [K(V(t))]{U} - {Q(t)l (5.1) where the global mass,damping and stiffness matrices are denoted by [M], [C(V(t))] and [K(V(t))] respectively. Clearly the stiffness and damping A matrices are functions of the external voltage V(t) applied to the ER fluid domains within the robotic arm. It must be emphasized that the damping matrix [C(Il(t))] in the equation (5.1) is, in general, not proportional to theistiffness and/or mass matrix. The nodal displacement is represented by the column vector {U} and the generalized force is denoted by the column vector {Q(t)}. Without loss of generality, the equation (5.1) may be reformulated as follows: [Mllfil + ([C]+[AC(V(t))]){fJ} + ([K]+[AK(V(t))l)lUl - {Q(t)l (5.2) where [K] and [C] are the global stiffness and damping matrices associated with zero electric intensity, and [AK(V(t))] and [AC(V(t))] are the changed global stiffness and damping matrices due to the change A in voltage from zero volts to V(t). 145 Consider now the modal analysis as did in Section 2.4, and in addition, assume that the ER fluid domains within the arm are not subjected to an external voltage. Then the equation (5.2) becomes {mm} + [CIII'II + [K]{U} = {Q(t)} (5.3) Now, by introducing modal coordinates {I7} and the modal matrix [] , i.e. , {Ul-[‘P]lnl (5.4) the equation (5.3) becomes {II} + IZIIéI + [mm = IBIS (5.5) 2 2 where [0] (- diag (0,w1, ...,wn)) is the modal frequency matrix with the pinned-mode natural frequencies wi, [Z] ( - [TC <1>]) is the transformed modal damping matrix, [B](- % [1 d¢1(o)/dx , . . . . , d¢n(o)/dx]T) T is the input influence matrix with the modal slope coefficient d¢i(o)/dx, and ii is the input torque. The transformed modal damping matrix [2] is coupled, since the damping matrix [C] is not assumed to be proportional to the stiffness and/or mass matrix. This coupling manifests itself in the non-zero values of the off-diagonal terms of the matrix [2]. In order to investigate the magnitude and hence significance of these terms, an experimental program was undertaken by considering a smart beam structure which has a similar modal characteristics of the proposed smart robotic arm structure. Figure 5.4 presents the schematic diagram of the employed smart cantilevered beam whose fabrication procedure is presented in details in Section 5.5. 146 Figure 5.5 presents the measured first and second mode shapes (magnitude and phase) of the proposed cantilevered beam. In Figure 5.5, the symbols and lines are the measured values and the curve-fitted interpolation respectively, and the magnitude of mode shapes is normalized as well as the beam length. These experimental values were obtained from the frequency response data utilizing the procedures described in Section 5.5. It is clearly evident from the phase of the mode shapes that the phases are all near 00 or 1800 at both 0 kV and 3 kV. Thus, the mode shapes are close to real [52], except in the region close to the nodal points where the phase transition occur. This is A / 223.5 m T / ./“’"““ 190 mm 1/ 27.5 i L. Silicone Rubber Aluminum ER Fluid 1.36 mm Section A—A Figure 5.4 Smart Cantilevered Beam Structure rnognfiude ofrnode shope phose ofrnode shape (degree) 147 00 0 WV " 0.- 3.0 kv ,' I 1 O- ”‘ " '7 I \ I I \ z ‘ I I ‘\ I I I .. i ‘- I I \ "‘ mode 2 I a I ‘ o ’ \ I ‘ / l \ ’ I \ ”- I '1 / V l .4 I ’z \ I l ' ‘ ’ A I I’ \ ’ (IS I , ‘ " r , ,I . \ 1 I ” I I .. I 1’ - I I I I " / ,’ I I O I I I 0 4—— , I’ ’ , I I O I -_‘ I I I I z’ I // ” l - 7 I I . .. I - I I I I ‘ I O 2— I, mode 1 I ,’ — i \ -". \ l _ ,- ‘I ’d’ d I OJ) CI2 0J4 C16 (18 1.0 27o.-..'...m.. ' 0,00k 0,-3.0kv - Figure 5.5 Mode Shapes of the Smart Cantilevered Beam 1148 particularly true for the employed cantilevered beam structure or the proposed smart robotic arm structure, since the distribution of all the eigenvalues of these structures are well separated. However, this fact may not be true for the complicated structures such as truss-like structures which have clustered eigenvalues. It is, therefore, reasonable to assume that each mode is uncoupled from the others. Hence, the off- diagonal terms in the transformed modal damping matrix [Z] can be ignored without causing serious errors. Assuming that the mode shapes are real, the coherence factor introduced in [75] is adopted here to investigate the viability of employing the modal matrix in the absence of externally-imposed voltage on the ER fluid domains, and also the situation where a non-zero voltage is imposed on these domains. The coherence factor is defined by ({¢§}T{¢§}>2 ID coh (5.6) i3 ({¢§}T{¢§})- ({¢§}T{¢§}) where coherence of 1 means two modes of i and j are identical, while coherence of 0 represents an orthogonal set of two modes. In equation (5.6), c and C represent the ith and jth cantilevered-mode modal i .1 vector respectively. Table 5.1 presents the experimentally-determined coherence factors for the first two flexible modes of the cantilevered beam structure investigated herein at four discrete state of voltage. Upon reviewing these data it is clearly evident that there is an insignificant difference between the coherence factor in the 0 kV state and the other states. Clearly, therefore the mode shapes are identical for the range of voltages considered in this study. 149 Table 5.1 Coherence Factor of the First and Second Mode at Different Voltage Mode 0 kV 1 kV 2 kV 3 kV First 1 0.993 0.993 0.993 Second 1 0.987 0.984 0.983 Ihmethese experimental observations, the modified equations of motion for the proposed robotic system in the presence of an electrical field intensity can be expressed as {n} + [c1{a} + [K]{n} - [BIG (5.7) where A T A . [Cl-[ZIH‘I’] [AC(V(t))][‘P]. [Z]=dlag(0. 2§1w1"”' 2§nwn) A T A [Kl-[OHM] [AK(V(t))][] Furthermore, from the experimental results presented in Table 5.1 [¢]T[A6(V(t))][¢]-[Z(§(t))] T A A (5.8) [é] [AK(V(t))][]=[0(V(t))] where [2(V(c))1 (- diag (o, 2c1)) is A the changed modal damping matrix due to the electric field and [0(V(t))] 2 2 ( - diag (0, w1(V(t)), ., wn(V(t))) is the changed undamped natural frequency matrix; namely 150 §.(V(t))w.(V(t)) - (§.w-)A - (§.w.)A 1 1 1 1 V(t)¢0 1 1 V(t)=0 2 " 2 2 (5.9) wi(V(t)) = (wi)x -(w.)x V(t)#0 1 V(t)=0 Equation (5.7) can be expressed in a state-space representation as follows. = A; + K(V(t))§ + ES X10 (5.10) ~ y =6; where state variable x, plant matrix K in the absence of an external ~ inaltage, input matrix 8, and the output matrix C is same as that defined in equation (2.33), respectively. And the changed plant matrix X(V(t)) due to the electrical potentials is given by O 0 0 0 o o o o O o o o o . - - o o 0 O 0 O o o o O 0 K(V(c))- 0 0 '121 '811 ° ' ° 0 O (5.11) 0 O O O o o o 0 O _0 0 0 0 . - . -a2n 181“; where 2" 2" 821- w1(V(t)). azn= wn(V(t)) 311' 2:1(V>w1(v). aln- 2cn> It is noted that only flexible modes are voltage dependent. 151 S.h Controller Formulation As did in Section 2.4, the equation (5.10) can be decomposed into two subdynamic equations by considering the rigid—body mode and the first few critical flexible modes as primary modes and remaining modes as residual modes. The two subdynamic equations are §=K§ +XVt§+Eu P P P P( ( )) P P ~ ~ ~ (5.12) = C x yP P P and XI. = rxr + A (V(t))x + Bru ~ ~ ~ (5.13) yr = Crxr where subscript (p) denotes the primary modes and (r) represents the residual modes. Before proceeding further, consider the stability problem of the system (5.12). Suppose the information about the changed plant matrix A KP(V(t)), which will be later termed as a pseudo-control force after reformulating it based on experimental observations, is not.awailable at A this moment. Then, the matrix KP(V(t)) can be considered as a time- A varying linear perturbation matrix due to the voltage V(t) with respect to the system in the absence of electrical potentials, that is, nominal system. Then following theorem for the stability can be established based on the Lyapunov function approach. 152 Round: 5.1 : In the following, the vector norm is taken to be Euclidean and the matrix norm is the corresponding induced norm, hence 1 2 ”A” - “max (AT A)) / where A (.) represents the operation of max(min) taking the maximum(minimum) eigenvalue. Lem 5.1 : Let Q be a symmetric and positive definite matrix. Then, .. 2 ~T .. ~ 2 Amin(Q)||xp|| s xp Q xp 5. Amax(Q)| prII (5.14) for any 32p (see [35] for the proof). Theorem 5.1 : Suppose the system (5.12) in the absence of electrical potentials is stable with a linear output feedback controller 5. Then, the system is stable if A (Q) .. A 1 "Ap(v(t))ll < 21m n(P) max (5.15) where P is a symmetric and positive definite matrix which is the solution of the Lyapunov equation ~ ~T PAc+AcP+Q-O (5.16) where Q is a given positive definite matrix and K -— K + E i2 ‘6 (5.17) C P P P where I? is output feedback gains associated with the controller 6. 153 Proof : With the equation (5.17), the state equation in (5.12) can be rewritten as N10 =3§+x Gt; 518 p Cp ¢<>>p < > Now, choose Lyapunov function candidate as ~ ~T ~ V x — x P x 5.19 ( p) p p ( ) It is noted that V(xp)> 0 for all xp# 0 and V(xp) 4 w as IprII 4 m. Now, 9(xp) < 0 is required to guarantee the stability of the system (5.12). w;)=§p; +flp; P P P P ~T ~T ~ ~ ~T~TA ~ — x A P + P A x + x A V t P x p(C p p p ¢<)) p +flpx(&a); P P P ~T ~ ~T~TA ~ = - x x + 2 x A V t P x 5.20 p Q p p p( ( )) p ( ) Hence, a; <01£ ( p) ~T ~ ~T ~T A ~ x x > 2 x A V t P x 5.21 p Q p p p( ( )) p ( ) Therefore, from the Lemma 5.1 and also ~T30%))?” IIKdMDIIHPIIHWI2 (5m) x t x S t x . P P P P P it follows that the inequality (5.21) is satisfied if 15h A .(Q) ||Kp(V(t))|| < 2Am1“(P> max Note that IIPII = AmaX(P). Physically, this theorem states that if the bound of the unknown A perturbed plant matrix KP(V(t)) due to the electrical potentials remains within the certain threshold value characterizing the closed-loop system: in the absence of the voltage, then the system (5.12) is stable. However, this theorem can be modified once the perturbed plant matrix is defined as a distributed ER fluid actuator, and the characteristics of the ER fluid actuator are identified through the experimental observation. Later'cn1 in this section, a modified theorem will be established for the stability of the system featuring the proposed hybrid controller which consists of the output feedback controller employed in the previous chapter and a pseudo-state feedback controller associated with the ER fluid actuator. The philosophy of developing the proposed hybrid control methodology for the system given by the equation (5.12) isrmnflwated by the complementary experimental program which is focussed on the development of an empirical actuator dynamics associated with the uniformly distributed ER fluid actuator. The empirical actuator dynamics should be A expressed as an explicit function of the applied voltage V(t) in order to formulate an appropriate control form. This can be accomplished by establishing the relationship between the stiffness and energy dissipation change and the applied electric field through the 155 experimental test. In order to obtain this relationship, frequency response was measured at various voltage levels. Figure 5.6 presents the measured frequency response of the smart ER arm without payload at two different discrete voltage levels. This frequency response was obtained by applying random torques (white noise) input from the dynamic analyzer to the servo-drive, and measuring output signals from the strain gage mounted at the root of the smart ER arm. In Figure 5.6, the peaks correspond to the clamped-mode natural frequency having the maximum strain energy at the root of the arm, and the valleys represent the pinned-mode natural frequency having the minimum strain energy. It is noted that this phenomenon is opposite to the case shown in Figure 3.15. It is clearly evident from Figure 5.6 that the natural frequencies of the arm, and also the arm's energy dissipation characteristics are both voltage dependent. ‘ '“ 0.0 kV ‘ '—- 5£)kV MAGNITUDE(dB) FREQUENCY(Hz) Figure 5.6 Measured Frequency Response of the Smart ER Arm without Payload 156 The increment of the natural frequency and the damping ratio of each mode at prescribed voltages relative to the no voltage state was obtained from Figure 5.6. And the distilled results suggest that the empirical actuator dynamics is governed by the following linear relationships of the form ciwi(V . w - aliv(t) ‘1 i 2 A (5.23) 2 = a21V(t) (.0. 1 where 0:11 and (121 are experimentally determined constants which characterize the increment of the damping ratio and natural frequency of A each mode as an explicit function of applied voltage V(t). Figure 5.7 presents this governing relationship for the first two flexible modes. In Figure 5.7, symbols are the measured values and the lines are the least- square error curves. From the slope of these curves,a -0.418, a -0.656, 11 12 a21 - 0.303 and a22- 0.666 were distilled. It should be noted that the governing relationship described by the equation (5.23) may be modified once a viable microstructural-level model of the system is developed. This microstructural-level model would include the different chemical composition of the ER fluids, the interface conditions between.t1u3 solid and fluid layers, the conductivity of the face materials (electrodes), and also the dielectric properties of the insulation between the electrodes. The integration of all this information may provide a viable model for the distributed "ER fluid actuator dynamics, which is a function of the electric field intensity. And also, the precise investigation on 157 1.8 0 FIRST MODE 1 A SECOND MODE A.- N... 1.2-~ .4! S 3‘ _ (E . NéH . O.6~ 5A1 0 .1 .A " o 0.0 . I I l I I l l 0.0 1.0 2.0 3.0 ELECTRIC FIELD (RV/mm) 1.8 0 FIRST MODE 1 A SECOND MODE A, .3: 1.2—4 k <2: 9:4 .1 < E 0.6“ .‘A. u .1 9., 000 I I r I I I l I 0.0 1.0 2.0 3.0 ELECTRIC FIELD (kV/mm) Figure 5.7 Governing Characteristics of the Employed Electro- Rheological Fluid Actuator 158 the ER effect to the different excitation amplitude will improve the accuracy of the dynamics. Now with the equation (5.23), the system equation (5.12) may be reformulated as follows. § - K § + E u + E V(t) p p p p p ~ ~ ~ (5.24) y - C x p p p where r0 0 0 o - . . o o ‘ o 0 o o - . - o o O O O O o o o 0 0 E _ o o -k21 -k11 - - - o o (5.25) p 0 O O O O O O O O 0 O O 0 o o o O o _0 o o o - - . -k2p ~k1p‘ k11‘ 2"‘11§1“’1' klp- 2°1p§pwp 2 2 k21' “zlwl' kzp‘ a2pwp V(t) - § V(t) P Clearly the equation (5.24) admits a hybrid control featuring the output feedback control {1 and the form of state feedback control V(t) . Even though external input voltage is imposed to an ER fluid domain in a fashion of the open-loop control, the distributed ER fluid actuator possesses a characteristic aspect showing the pseudo-state feedback controller which produces the associated pseudo-control force in a fashion of the closed-loop control. In fact, this pseudo-control force accounts for the increased stiffness and damping properties of the structure due to the imposed external voltage. Figure 5.8 presents a 159 71- H {>1 X! + U) C + 7: < A f? V @ ER Beam l I I B 8 K 8200 Ling D.S. Force Transducer V411 Shaker I » B & K 2635 ,7" Charge Amplifier HP 35660A Hafler P500 5‘ Dynamic Signal Amplifier Signal Analyzer Results Figure 5.10 Experimental Set-Up for the Frequency Response Measurement of the Smart Cantilevered Beam 164 forcing function (white noise) was applied from a dynamic signal analyzer, model 35660A, manufactured by Hewlett-Packard, to the structure. This forcing signal was passed through an amplifier, model P500, manufactured by Hafler and then employed to energize a shaker, model V411, manufactured by Ling Dynamic Systems Limited. The amplitude of this excitation was kept below a certain level to reduce the influence of nonlinearity of the ER fluids and other parasitic effects such as air damping, cable drag,etc. A force transducer, type 8100, manufactured by Bruel and Kjaer was mounted on the head of the shaker in order to determine the input force to the beam, whereas the displacement probe was positioned at several locations along the beam in order to measure the output signal which was monitored by a non-contacting displacement transducer, model KD2400, manufactured by Kaman Corporation. Then these two signals were processed by the signal analyzer to obtain the frequency response (magnitude and phase) corresponding to each location. In order to reduce the influence of the noise and get more consistent results ,eulaverage of 40 times measurements was used. Figure 5.11 presents one of the frequency responses of the cantilevered beam at two discrete voltage levels using the displacement probe mounted on the tip of the beam. The phase-change in the ER fluid was initiated by energizing the upper and lower faces of the beam specimen. The various voltages were generated by a Hewlett-Packard dc power supply, model 6255A, and amplified through a dc-to-dc converter, model SC30, manufactured by Standard Energy. The desired input voltage was measured using a dc high voltage probe, model 3411A and digital multimeter, model 3566A which are both manufactured by Hewlett-Packard. 165 40 :18 20.5 Hz _ 0 RV --"‘ 3.42 kV 1 10 dB - Q/div _J 3 E _ E 176.0 52 L- .......................... -4O 1 1 1 1 1 1 1 l I Start: 15 Hz Stop: 215 HZ FREQUENCY (Hz) RMS‘ 4° Figure 5.11 Frequency Response of the Smart Cantilevered Beam The detail-design features of the smart ER arm are schematically presented in Figure 5.12, and the salient features of the fabrication process for the smart robot arm incorporating the ER fluid are highlighted in the flow diagram presented in Figure 5.13. RTV silicone rubber adhesive was selected as an insulator in order to provide both excellent bonding between the electrodes and to also furnish a high dielectric constant. The arm was fabricated by employing an appropriate die to generate the hollow section and also to control both the thickness of the ER fluid layer and the insulator . The die was subsequently removed after partial curing of the sealant, and the arm was maintained at room temperature for 24 hours in order to achieve the full strength of the silicone rubber adhesive. During this time an air blower, 166 manufactured by Air Control Installations, Ltd. , U.K. , was employed to remove the acid from the silicone rubber by exposing the arm to an air stream for two hours. After curing at room temperature, the arm subsequently cured at 65°C for 2 hours in a hygrothermal chamber, model TH27S, manufactured by BMA Inc. The arm was then exposed to room temperature for 24 hours before checking the dielectric behavior. After checking the dielectricity of the arm, an acrylic sheet was employed at one end where it needs to be inserted into the hub. The hollow arm, which featured a cavity of uniform rectangular cross-section, was filled with A __ / 72.4 cm __ ll ——— .//——O'5 cm . 70 cm 1/ l ,l.4_cm| Silicone Rubber Aluminum ER Fluid 2,2um1 Section A-A Figure 5.12 Schematic Diagram of the Smart ER Arm 167 I Preparation of Structural Components 6 Bonding I l ItfitialPartialCuring l - 20 minutes at 245°C Runovfng Die 6 4 Air Drying Intermediate :aItial Curing - 24 hours at 245°C - 2 hours at 65°C - 24 hours at 245°C Not Acceptable Check Dielectric Filling Beam with ER Fluid and Scaling of Entry Port - Final Curing ’ - 30 hours at 245°C Figure 5.13 Flow Diagram for Fabrication Procedures of the Smart ER Arm “.~. \A‘ , § . I“. n ~ Figure 5.14 Photograph of the Smart ER Arm the silicone-oil-based ER fluid. The ER fluid volume fraction to the total volume was 36 percent. Finally a Measurements Group strain gages, type WD-DY-ZSOBG-350 were mounted at the root of the arm. It is noted that the fabrication procedure of the cantilevered beam shown in Figure 5.4 is same as that for the smart ER arm. Figure 5.14 presents a photograph of the smart ER arm clamped on the rigid hub. From the measured frequency response shown in Figure 5.6, system model parameters such as natural frequencies and damping ratios were identified. The modal slope coefficient d¢i(o)/dx was determined as followings. The open-loop transfer function G°(s) in the equation (4.2) is also given by [127] 169 l0 2 1 p (s /0. + 2§is/Oi + 1) NH 2 (5.28) I 5 i=1 (5 /w. + ZCiS/wi + l) H 1" where p is the number of primary mode, 0i is the ith clamped-mode natural frequency and wi is the ith pinned-mode natural frequency. (hi the other hand, the 60(3) in equation (4.2) can be expressed as 1 1 P 2 2 2 G (s) = ——*—3‘ + —_— 2 (d¢.(o)/dx) (l/(s +2§.w.s+w.)) (5.29) o I 5 IT i-l 1 i i 1 T Therefore, the modal slope coefficient d¢(o)/dx can be determined by equating equations (5.28) and (5.29). Table 5.2 presents the measured Table 5.2 Measured System Model Parameters of the Smart Robotic Arm Parameters Values wl (Hz) 15.2 ”2 (Hz) 30.5 01 (Hz) 4.3 02 (Hz) 17.1 {1 0.0526 {2 0.0426 d¢1(0) dx 3.12 d¢2(0) dx 5.38 170 system model parameters of the smart ER arm in the absence of the electrical field. In this investigation, rigid-body mode and first two flexible modes were considered as the primary modes. 5.5.2 Results and Discussions Figure 5.15 presents the measured step responses to a commanded step angular position of 7.2 degree with the compensator zero Zc--2. The feedback gains of Kp=0.5 and Kv=0'25 were employed. The external input voltage (pseudo—control force) of 2.0 kV/mm was employed continuously throughout the motion to the ER fluid domain. It is clear from this figure that the responses with and without input voltage are basically same except the tip deflection response of the arm. This aspect may be expected since the mass of the robotic arm remains same regardless of the input voltage, and the pseudo-control force does not affect the rigid- body motion. The amplitude of the tip deflection of the robot arm in the absence of the electric field was reduced up to 35 percent by employing the input voltage of 2.0 kV/mm. It is also observed from this response that the settling-time (decaying time) is about 1.3 second with the input voltage, while it is about 1.6 second without the input voltage. This can be interpreted as the enhancement of the positional accuracy of the robotic system. No second flexible mode effect to the total tip deflection is observed for both cases. The corresponding simulated results are presented in Figure 5.16 which compare very favorably with the measured results shown in Figure 5.15. For the second comparison of control performances with and without the input voltage, the compensator zero ZCO=-5.6 was chosen. Figure 5.17 171 shows the measured step responses with feedback gains of Kp-l.4 and KV-0.25. From the angular position response, the rise time was obtained by 0.42 second for both cases, and the maximum overshoot was about 0.89 degree without the input voltage and 0.93 degree with the input voltage. This is basically same for both cases, and the input torque is alt“) same for both cases as expected. The tip deflection without the input voltage was reduced up to 20 percent by employing the input voltage for whole duration of the response. With these feedback gains, it is also observed jittering phenomenon as in the case of the aluminum or composite arm (see Figure 4.11) (hue to the hysteretic joint friction of the motor. Figure 5.18 presents the simulated step responses corresponding to the measured. results shown in Figure 5.17. In order to implement the bang—off-bang type input voltage shown int Figure 5.9, on-off switched response was obtained by imposing a non-zero input voltage at certain time. Figure 5.19 presents the measured step responses with the same output feedback gains used for the results in Figure 5.17. The input voltage was imposed from the beginning up to the peak time (0.61 second) of the angular position response in Figure 5.17, and was released throughout the response. It is observed that the tip deflection without the input voltage was reduced up to same amount as achieved in Figure 5.17. This experimental result illustrates the speed with which the smart ER.arm can.respond to the imposed input voltage (pseudo-control force) in order to synthesize the desired vibrational response of the flexible robot arm. It is noted that the determination of the best switching time depends on various parameters such as rise time, peak time and settling-time, and also depends upon the system performance requirements such as minimum time or minimum effort problem. .mN.0 .I. M van m.0 I x mo mcwmc songbook £33 momaoemom moum pouammmz mad mummfim > a $80sz 0 s. n N _. 0 a b h — s a p u — L p h h b h . p n F . . . . Fl ... 1N 000:? SQE 53 II. . 002.9 Sac. 505$ . 2 n 7 1 Seams: n v n N F o p p p — _ b b L w p P b p h h p u p _ b n h F O ..N 1* 10 1m 00029, 59.; 5? II 003.? “:9: 3053 . 0— (ww)NO|iO_Tl:130 dLL (fiap)NOLusoa HV'IOONV AoomVMZF m e n N F o h b b h F a p p h — F p L n ? b b n P — p P p - No.0I J , i. rfl Hood Tuod L H _. .....od 1 700.0 000:9 39: £3 ll . omoxo> SeE 5055 1 no.0 Aoomvmifi m e m. N P 0 . . p p _ - p h p P b p p p p L . u p — p p u h O.F‘ .................................................................................. 10.0 to; ON 106 EE\>x o.~ ..I . EE\>x 0.0 . OJ. (w'N)3nOElO.L lflle (mm/maovnm 10le 173 .mN.0 x .«o mcwmu xomnpmmm £33 momsoemom moum pouoasaam win muswwm > . a n M 28 m o u 08m}: e n N _. 0 L L p — p p h p w L P p . w p L . p b p p h . Flu 1N 000:? 39: £3 I . 003.3 59: 5055 n A8302: v p L .— — L . L p m . n r L W L . L . ~— L L L u 00 IN ...v 10 10 omozo> uses 53 I: s 603.9, «:9: 3053 . 0_. (ww)Nou:>31~130 an (fiapmowsoci awnorw Scams: ¢ n N F 0 L L L L — p r L L p _ u L L ~ P L L . — . P F L No.0-l 100.0 [No.0 w $40.0 100.0 omogo> 305 53 I . ooogo> 59: 5055 0. 00.0 A8302: .v n N _. 0 L p u p b . p u n b L Fr _ — u p . L — p p - r onFIl. j .................................................................................. IO.O 10; ON. 10.n EE\>x 0.N II . EEx>x 0.0 e (w'N)3noao.L inle (ww/M)39Vl‘lO/\ 1!)le 174 a Aoomvmzc. omo:o> 59: 5:» ll 003.9 SaE L305? m 330%: év L L h L L L L m L L L L W L L L L —F L L L L O o T 1N 1 ...v W 1m w 10 10L 003.9 59: 53 l ooozo> Sec. 3053 . (59P)NOLLISOd whom (ww)~ouom3c1 dLl .mN.0 > 1 x can «A I M mo mchU xomnpoom LENS momcoamom moum pousmmox mad ouswwm “03va0. a w m N L o L L L p p L p L L— b b P — L L L n — L b L L h @000! 1 00.0 00.0 0L0 000:? 39: 53 I. ooogo> Sec. «305:; . .vN.0 A0352: 0 ¢ n N L 0 L L L L b L L Li P — L L L L — L L LL L L L L L O F'- ................................................................................... 10.0 10; ON. 10.0 EE\>V_ o.~ I . EE\>.. o.o T (w'N)3noaOi .10le (ww/MBQVDOA inle 175 .mN.0 «A 1 M mo mcwmo xomnpoom €53 newcommmm moum pmumaaawm 0H6 muswwm > a ... x can 063%: n N L 0 L L L L — L L L L b L L L L .— L L L L .V1 .. 10 003.9, 59: £25 I 003.9 59: 5053 1 m 380%: .v m. L L L L — L P L L W L L L L w L L L L o O 1N [v 1 1.0 10 10L 002.3 59.; £5 I 00029 «:9: Song; 1 NL (UJUJ)NOl103'L:lEJO dli (Bapmomsoa awnsnv 330%: m L. m. N L 0 L L L L F L L L L b L L L L b L L L L h L r L L mO.OI 100.0 100.0 10L.0 omo..._o> SQE :23 I 1 002.? 39: .305} VNd Loomvmz: m e n N L 0 L L L L _ fILI L L _ L L L L h L L L L b L L L L O.Fll ................................................................................... 10.0 10.L 0N 10.0 EE\>x ON I 1 EE\>V_ 0.0 T 0% (11110300301 .10le (wW/M)39\/i10/\ 111le 176 3i) .1 ’é‘ 2.0 E .4 \\ >. . i‘, 15 . g 1.0— .—J .1 O 2> 4 i... :3 L a. Z 0.0... ......... .—1 O I I I I l I I I I I I I I I I I I Ifi l fi I I I 0 1 2 3 4 5 TlME(sec) 8 without input voltage — with input voltage A 5-1 ‘z’ 9 I... S uJ c3 0. p: -I L “4 T 1 1 1 I r 1 If r I I 1 1 1 I r 1 1 1 I 1 1 ft 0 1 2 3 4 5 TlME(sec) Figure 5.19 Measured Step Responses with Bang-off-Bang Type Input Voltage. Kp-l.4, Kv-0.25 177 'The proposed hybrid controller gets its distinguished credit when the system is subjected to a disturbance which has a forcing frequency close to the clamped-mode natural frequency of the arm. In other words, the robustness of the employed PD type output feedback controller to the imposed disturbance can be improved by employing the pseudo-control force which accounts for the increased stiffness and energy dissipation of the robot arm. With the introduced disturbance, the system given by the equation (5.24) becomes > (13>)p + (80+ AB(p) + AB(£))u (6.6) xp(to) - xp0 where To 1 o o - - . o o ' 0 0 0 0 . - . 0 0 0 0 0 1 o o o 0 0 2 0 O -w01 -2§1001- - 0 0 0 A... o o o o o O 0 0 (6.73) o C O O O O C O O C o O 0 o 1 2 2 o O 0 'w ' w . 0q (q 0q. 189 To 0 o o - - - o 0 o o o o - . - o o 0 O O O o o g 0 O O 0 p11 -p21 - . . O 0 AA(p) E O O O O O O O O O (6.7b) o o o o . - - o o .0 o o o - - - -plq -p2q - 2 p p _ p p11 wOlwr1(wrl+ 2) ' p21 2i1“’01“’r1 2 1 p p _ p plq ququ(qu+ 2) , p2q zngqurq ”o o o o . - . o o I O 0 O 0 o o o 0 O O O O 0 o o o O 0 o o -211 -221 . . . o o AA(£) = o o o o o o o o o (6.7C) O O O O o o o o o o o o o . . . - - L 1Q 2Q - 2 - 02 2 ( 2 + 2) 2 - 2 w mg 11 01°r1 ”:1 ' 21 ‘1 01 r1 2 2 2 2 - w w w + , 2 - 2 w w lq Oq rq( rq ) 2Q {q 0q rq ’o o o o . . - o o ‘ o o o 0 . . - o o o o o o - . . o o o 0 -p21 0 - - - o o AA(p’2)- o o o o o o o o o (6_7d) o o o o - . - o o o o o o - - - -p2 o . q . 2 2 1 p 3 a 9 p21 2w01wrlwrl , pflq 2w0qququ B - [o b o b - - - o b ]T (6 8a) 0 oo 01 Oq ' - P P . . . P T AB(p) [o boobr0 o b01br1 o boqbrq] (6.8b) 2 2 2 T _ . . . .8 AB(2) [O bOObrO O bOlbrl O qubrq] (6 c) It is clear from the equations (6.7) and (6.8) that the uncertainty Inatrices AA(-) and AB(-) can not be decomposed into a manner such that the matching conditions are met, i.e., AA<-> - Bo D(.) AB(°) " BO E(°) Where D(-) and E(-) will be defined in the subsequent section. Therefore, without consideration of the matching conditions, a nonlinear robust state feedback controller which guarantees the stability of the system (6.6) for arbitrary initial conditions x will be designed in the next p0 section. 6.3 Controller Design For generality, consider following dynamical system described by the state equation i - (A0 + AA(01(t)) + AA(02(t)) + AA x(t) + (B0 + AB(ol(t)) + AB(02(t))) u(t) (6.9) x(to) - x for t e [t0,t o' 1] where x(t) 6 RF is the state, u(t) E R? is the input, and A0, B0 are the prescribed nominal system matrices of appropriate dimensions. In equation (6.9), AA(al(t)) and AA(02(t)) are plant uncertain matrices associated with uncertain parameters 01(t) and 02(t) respectively, and AA(al(t), 191 02(t)) is the plant uncertain matrix associated with the both parameters 01(t) and 02(t). AB(al(t)) and AB(a2(t)) are input uncertain matrices depending on parameters 01(t) and 02(t), respectively. Now following assumptions and definitions are made before proceeding further. Assumption 6.1 : The uncertain matrices AA(-) and AB(-) are prescribed functions which are continuous on RP. Assumption 6.2 : Uncertain parameters .14 p 01( ) . R. F1 C R .1_, p 02().R 1‘ch (6.10) are Lebesgue measurable, where 1‘1 and I‘2 are prescribed compact subsets of appropriate spaces. Assumption 6.3 : The pair (AO,BO) is controllable. This assumption in fact can be released for.the system (6.6). Assumption 6.4 : All state variable x(t) are available from direct measurements or state estimators. Remark.6.1.: In this chapter, the vector norm is taken to be Euclidean and the matrix norm is the corresponding induced norms, hence 1 2 [IAII - (A )(ATA)) / where A )(.) denotes the operation of max(min max(min taking the maximum (minimum) eigenvalue. 192 Definition 6.1 : Given a solution x(o): [t0,t1] n RP, x(to) ==xo for the equation (6.9), it is said that the solution is uniformly bounded if' there is a positive constant 6(x0) < 00, possibly dependent on initial but not on initial time t such that value x0, 0, ||x(:)|| s 6(x0) for all : e 1:0,:1] Definition 6.2 : Given a solution x(o): [t0,w) 4 RP, x(to) = x0 for time equation (6.9), it is said that the solution is uniformly ultimately bounded with respect to a set S if there is a non-negative constant and S, but not on t such that 1(xo S) < m, possibly dependent on x 0, 0 x(t) 6 S for all t 2 t0+ r(xo,S) Definition 6.3 : The equation (6.9) is said to be decomposable if there exist continuous matrix functions D(al(t)): Planpx“ D(02(t)): r21RPX“ D(al(t),02(t)): (r1,r2)422x“ (6.11) E(al(t)): rlenpx“ E(02(t)): r2422x“ and AA(°) and AE(-) having following properties. AA - Bon(a1) 193 AA(02(t)) = BOD(02(t)) + AK(02> (6.12) AB(al(t)) = BOE(01(:)) + A8(al(t)) AB(02(t)) — BOE(02(t)) + AE(a2(t)) Remark 6.2 : The argument t is often omitted throughout this section when no confusion can arise. It is observed from the equation (6.12) that the uncertain matrices AA(-) and AB(-) are decomposed into two parts ; the matched part BOD(-) and BOE(-), and the mismatched part AA(-) and AE(-). It is noted that this decomposition is highly nonunique. In some sense, a 'good' decomposition may be achieved by lumping uncertainties as much as into the matched part, hence minimizing the portion of the mismatched part. One way to obtain a 'good' decomposition is to apply the least-square error method [28]. For example, by letting AA<~> = BOD(-) AB(o) =BOE(°) (6.13) following decomposition can be accomplished in the least-square error sense. T -l T D(01) - (B0 BO) BO AA(01) T -1 T D(a2) - (B0 BO) BO AA(02) 194 T -1 T D(01,a O 2) 2) T -l T E(01) - (BO BO) BO AB(ol) 3(02) - (33 BO)'1 Bo AB(oz) (6.14) AA(01) - AA(01) - BO D(al) AA(02) - AA(02) - BO D(02) AA(01,02) - AA(a B D(a 1’02) ' o 1'02) AE(01) - AB(al) - BO E(al) AB(OZ) - AB(02) - B0 E(02) The purpose of this problem is to design a state feedback controller which assures that every solution of the equation (6.9) is uniformly bounded and uniformly ultimately bounded with respect to a set S to be specified, no matter what the uncertainties and initial conditions are. For the sake of generality, suppose that the nominal plant matrix A0 is not stable. However, there exist a constant feedback gain matrix K such that AC-A +B K (6.15) is stable, since the pair (A0, B0) is assumed to be controllable. This implies that the real parts of all eigenvalues of Ac are negative. Now, consider a state feedback controller which has the form of u - Kx + p(x) for all x 6 H2 (6.16) 195 where the function p(-): 22422 is such that I 'r B0 PX . T - [IBT lel (p1(x) + p2(X) + p3(X)). 1f IIBO PXII>€ 0 p(x)-< (6.17) B3 Px T - -———:—-—- (p1(x) + p2(x) + p3(x)), if [[30 PxIISe where e is a prescribed positive constant, P is the solution of the Lyapunov equation AT 9 + P A = -Q (6 18) c c ' for a given constant positive definite (nxn) matrix Q, and pi(-): RF 4 R1 will be specified subsequently. Substituting equations (6.12) and (6.16) into the equation (6.9) yields x - Acx + Bop(x) + Bo(em1(x,t) + em2(x,t) + em3(x,t)) + 31(x,:) + 52(x,:) + 33(x,:) (6.19) where em1(x,t) - D(al)x + E(01)Kx + E(al)p(x) em2(x,t) - D(02)x + E(02)Kx + E(02)p(x) em3(x,t) - D(al,az)x E1(X,t) " (AZ(01) + A§(OI)K)X + AB(gl)p(x) (620) 32(x,:) - (AA(02) + AE(02)K)x + AE(02)p(x) 196 E3 = (AK>x for all (x,t) 6 RP x R1. In equation (6.19), BOemi(x’t) and ei(x,t) are called as the matched error vector and mismatched error vector respectively. From the equation (6.20) Ilem1(x,t)|| 5 max ||D(al)x|| + max IIE(01)KxII aleFl alEFI + max IIE(al)IIp1(x) 9 p1(x) aleFl Ilem2(x,t)ll s max IID(02)xII + max IIE(02)KxII (6.21) 0 EF 0 EP 2 2 2 2 + max ||E(02)||p2(X) g p2(X) 026F2 A llem3(x,t)ll s max IID(01,02)x|| = p3(x) a GP lerl ”22 The defined parameter p (x) in equation (6.21) can be solved provided max ||E(al)|| < 1 OIEFI (6.22) max ||E(02)|| < l 026F2 Thus, with conditions (6.22), pi(x) are obtained as follows. p1(x) - [1 - max IIE(01)I|]-1[ max ||D(al)xll + max I|E(ol)Kx|| alePl OIEFI alefl p2(x) - [1 - max ||E(02)II]-1[ max ||D(02)xll + max ||E(02)lel (6.23) 02€P2 02€P2 026P2 197 p3(x) - max I|D(al,02)xll a GP 1 1 02€F2 Now define 31(x): RP* R: by ~ A P1(x) - 01 IIXII ~ A P2(x) - 02 llxll ~ A p3(x) 3 03 lell where a1 - [1 - max ||E(al)||]'1[ max ||D(al)|| + max I|E(01)KII 6P 0 61‘ a 61‘ ”1 1 1 1 1 1 a2 - [1 - max IIE(02)II]-1[ max ||D(02)|| + max IIE(02)KII azerz 02€P2 026F2 a3 - max I|D(al,02)ll a GP OIEPI 2 2 Thus, from the equations (6.23) and (6.24) p1(X) S 31(X) IA p2(X) 32(X) p3(X) S 33(X) for all x 6 RP. Furthermore, from the equations (6.21) and (6.26) llem1(X.t) + em2(X.t) + em3(X.t)I| s Ilem1(x,t)|| + Ilem2(x,t)|| + llem3(x,t)ll S p1(X) + p2(X) + p3(X) (6.24) (6.25) (6.26) 198 s Zl + E3 = (a1+ a2+ a3)||x|| <6-27> Now in order to handle the mismatched error vector in equation (6.20), following definitions are introduced. Definition 6.4 : Let the measure of mismatch be defined as follows. :2 I max IIAK(01)|| + max ||A§(01)K|| + a max IIA§(01)|I 1 1 aleFl alefl aleFl m2 é max IIAK(02)II + max ||A§(02)K|| + a2 max ||a§(a2)|| (6 28) 6P a 6P 0 GP ”2 2 2 2 2 2 A ~ M3 - max IIAA(01,02)I| a GP lerl ”2 2 Definition 6.5 : Define the mismatch threshold as A (Q) * A min M - 2A (P) (6.29) max Definition 6.6 : Define the mismatch index as A (”1 + M2 fl - x M + M3) (6.30) It is noted in Definition 6.4 that M1 and M2 account for the mismatch associated with uncertain parameters 01(0) and 02(-) respectively, and M3 accounts for the mismatch associated with both parameters. From the equations (6.20) and (6.28) 199 ||51(x.t>|| 5 M1 ||X|| ||32(x.t)|| 5 M2 leII (6.31) ||53(X.t)|| 5 M3 IIXII Thus , ||31ll + l|53(x,t)ll 5 (M1 + M2 + M3) ||x|| (6.32) x The mismatch threshold M defined in (6.29) implies a limit on the amount of the mismatched uncertainties, defined as the measure of mismatch Mi’ which can be tolerated. The mismatch index [3 in (6.30) indicates the normalized total measure of mismatch (M1+ M2+ M3) ‘flith respect to the mismatch threshold. Now consider the uniform and uniform ultimately bounded sets which satisfy the conditions given in Definitions 6.2 and 6.3. With a hypothesis that the total measure of mismatch is not two large, i.e., fi<1 or (M1 + M + M3) < u* (6.33) 2 form the scalar function h(-): [0, w)+R given by A 2 h(n) - (Amin(Q) - ZAmax(P)(M1+ M2+ M3))n -((al+ a2+ a3)6/2)n (6.34) Then, letting A 5 max ( n : h(n) s o } (6.35) 200 yields (a1+ a + a3)e 2 fl ' (6.36) 2(Amin(Q) - 2Amax(P)(M1+ M2+ M3)) It is noted that the E in (6.36) implies the radius of the largest closed ball which is wholly contained within the bounded set (see Theorem 6.1). Remark 6.3 : It is observed from the condition (6.33) that the * mismatch threshold M depends on eigenvalues of the matrices Q and P. The * * larger'li tolerates the larger uncertainties. The maximization of the M will be discussed later in this section. Since Ac in equation (6.15) is stable and the matrix Q in equation (6.18) is positive definite, and hence the matrix P is positive definite, following Lyapunov ellipsoids can be constructed. X(k) - { x e R? | xTP x s k ) (6.37) where k is positive constant. Now define the smallest ellipsoid containing ball 8(5) as 2 km 9 min ( k | 3(5) g X(k) ) = Amax(P) ‘5 (6.38) And also, for iven an E > km,and an initial condition x X(E),define 8 Y Y 0 the constant P x ' (6.39) and 201 ... ~ 1 2 1 2 c(x0,k) 9 min ( h(n) : (k/Amax(P)) / s n s (kO/Amin(P)) / } (6.40) Finally, following theorem regarding to the boundedness of the solution of the equation (6.9) can be established. Theore-.6.1 : Consider the dynamical system given by equation (6.9). And suppose that the conditions given by the equations (6.22) and (6.33) are satisfied. Then, the controller given by the equation (6.16) assures the uniform boundedness of every solution x(-):[t0, t1] 4 Rh, and uniform ultimate boundedness of every solution x(-) extended over [t0,m). ~ m . . . . Furthermore, given any k > k , every solution corresponding to initial condition (x0,to) E Rpx R is uniformly bounded with F 1/2 ... (Amax(P)/Amin(P)) leoll for xo ¢ X(k) 6(xo) - 1 (6.41) .. 1/2 .. L(k/Amin(P)) for x0 6 X(k) and is uniformly ultimately bounded with respect to a set S e X(E) with r (k0- E)/c(xo, E) for xo ¢ 3 ?(xo,3) - < (6.42) 0 for x0 6 S Proof : Denote the right-hand side of the equation (6.19) by A f(x,t) - Acx + B0p(x) + Bo(em1(x,t) + em2(x,t) + em3(x,t)) + 31(x,t) + 52(x,c) + 33(x,c) (6.43) 202 where f(x,t) is Caratheodory functhn1; f(x,t) is continuous in x and x - Lebesgue measurable in t. Now let X(k ) be a superset of X(k) containing all initial conditions. Then for given initial conditions, f(x,t) can be * uniformly bounded over the domain D e X(k ) x R from the observation IA llf|| IlAcll lell + (.1. a2+ a3>|lBoll llxll + IIBO(em1(x,t) + em2(x,t) + em3(x,t))|| + ||51 + 52(x,t> + E3cx.c>|| IA IIACII llxll + 2||Boll llxll (6.44) + (Ml+ M2+ M3)||x|| lll> g(X) * Since g(-) is continuous [19] on X(k ),there exist a finite maximum value A k - max g(x) (6.45) xeX(k*) A Thus, it is clear that k plays a role of a uniform bound for.f(x,t) over domain D. Therefore,it follows from [43] that there exist a solution x(-) :[to,t1] 4 RP. The extension of this solution will be shown later of this proof. Now, introduce a Lyapunov function candidate given by V(x) - xTP x for all x 6 RP (6.46) where V(o):RP4R is continuously differentiable. To show that it is indeed a Lyapunov function for the system (6.9) with controller (6.16), consider the Lyapunov derivative as follows. 203 T 9(x,t) - x P x + xTP x - - xTQ x + 2(33? x)T (p(x) + em1(x,t) + em2(x,t) + em3(x,t)) + 2 xTP (31(x,t) + 22(x,c) + 33(x,t)) - xTQ x + 2(33? x)T (p(x) + ((ng x/||ng x||) (6.47) IA (p1 + p2>> + 2 xTP (51(x.t> + 32(x,c) + E3(x,c)) It is observed from the equation (6.17) that, if IIBg P xII > e, the term of p(x) + (BE P x / IIBg P x||) (pl(x) + p2(x) + p3(x)) vanishes, and if [[83 P xll < 6, then it has maximum value of e (p1(x) + p2(x) + p3(x))/2. Hence, from the equations (6.27) and (6.32), the equation (6.47) becomes V(x,t) s - xTQ x + ((a1+ a2+ a3)e/2)I|x|| 2 + 2Amax(P)(M1+ M + M3)I|xll (6.48) 2 Thus v(x,t) < 0 (6.49) for all t E R; and all x such that TQ x - ((a + a + a )6/2) x x 1 2 3 2 - 2Amax(P)(M1+ M + M3)||x|| > 0 (6.50) 2 However, from the Lemma 4.1 in the previous chapter 204 T 2 2 Amin(Q)||x|| s x Q x s Amax(Q)||x|| (6.51) where Amin(Q) > 0, since Q >0. Thus,the condition (6.49) is satisfied for all t e R; and all x such that 2 (Amin(Q) - 2Amax(P)(M1+ M2+ M3))||xll - ((a1+ 62+ a3)e/2)llx|| > 0 (6.52) Hence, it can be said that the inequality given by the equation (6.49) is guaranteed for all te R;and all x ¢'B(3),where E is given by the equation (6.36), provided.that the condition (6.33) is satisfied. It should be noted, however, that the condition (6.52) is only a sufficient condition. Now consider the smallest Lyapunov ellipsoid that contains the ball 8(3), that is, X(km), where km is given by the equation (6.38). From the Lemma 4.1, 2 T 2 Amin(P)||x|| s x P x s Amax(P)||x|| (6.53) where Amin(P) >0, since P >0. The condition (6.49) implies boundedness of all solution. That is, if x(o):[t0,tl]4RP with initial condition x(tO)-xo is a solution of the equation (6.9) with the controller given by the equation (6.16) and allowable uncertainties, then there are following two cases to be considered: if x0 Q'X(E), that is, x(t) e X(ko) for all t e [to,t1],then following inequality can be obtained in view of the equation (6.53). 205 2 T T 2 0 < Amin(P)||x(t)|| s x (t) P x(t) 5 x0 P x0 5 Amax(P)||xO|| (6.54) If xo 6 X(E), that is x(t) e X(E) for all t e [t0,tl], then following < < Therefore, 6(x given by the equation (6.41) can be achieved from the 0) equations (6.54) and (6.55) for x0 ¢ X(E) and x e X(E) respectively. And 0 from this turLform boundedness and existence of the solution, it follows that a solution x(o): [t0,t1] 4 RP can be extended over time interval [to.m>. The proof of ultimate boundedness follows again from the condition (6.49). If x0 6 X(E), it is an immediate consequence of the boundedness result, hence it suffices to take ?(x0,3) - 0. If xO ¢ X(E), the V(x(t)) decreases as long as x(t) Q X(km) and x(t) reach the boundary 6X(E) in a finite time. From the equations (6.54) and (6.55) ~ 12 12 (k/Amax> / s llx(t>|| s (ko/Amin

) / (6.56) whenever x(t) 6 X(kO)/X(E)>. Hence, from the equations (6.48) and (6.52) v(x,t) IA 'h(llx(t)||) . ~ . ~ 1/2 ~ 0 1/2 S -m1n {h(n). (k/Amax(P)) S n S (k /Xmin(P)) } --e(xo,E) < 0 (6.57) 206 where 2 h(||x(t)|| - (2(M1+ M2+ M3)||P|| - Amin(Q))||x|I + ((al+ 02+ a3)e/2)||x|| (6.58) ~ * * - ~ Suppose x(t) g X(k) for all t e [t0,t ], where t - (kO- k)/c(x0,k). Then from the equation (6.57) ' E)t* s - (k0- E) (6.59) V(x(to+ t*)) V(xo) s - c(x o, k. Then this contradicts the supposition. Thus IA so that V(x(to- t*)) A x A ~ there exist a time t S t such that x(t) E X(k) = S. Therefore, from the x * ~ boundedness and x(t) e S for all t 2 t, t - r(xo,S). This concludes the proof of the Theorem 6.1. Now, consider the maximization of the mismatch threshold M* given by the equation (6.29) as noted in Remark 6.3. It is clear from this equation that the threshold depends on the eigenvalues of the positive definite matrices Q, and P which is a solution of the Lyapunov equation associated with closed-loop eigenvalues of the nominal system. Hence, the * maximization of M can be accomplished by considering either given matrix Q or the linear feedback gains in equation (6.16). Herein, the * maximization of the M is achieved by considering the matrix Q with fixed feedback gains. Before establishing a theorem regarding to the * maximization of M , following lemma is introduced. 207 lemma 6.1 (see [113] ) : Let P be the unique positive definite solution of the Lyapunov equation T ~ ~ - Ac P + P Ac = -Q = - 1Q (6.60) where 1 is positive scalar constant, and Ac and Q are in equation (6.18). Then * .. .. M = Amin(Q)/24max(P) - Amin(Q)/2Am (P) (6.61) 3X where P is the solution of the Lyapunov equation (6.18). Proof : Rewrite the equation (6.60) as T ~ ~ AC (P/v) + (P/1)Ac - - Q (6 62) Since P >0 and P >0, and also there is an unique solution of the equation (6.18) and (6.60) respectively, it follows from the equations (6.18) and (6.62), P - P/V . Hence Amax(P) = Amax(P)/7 (6.63) Thus the equation (6.61) is satisfied Theorem 6.2 : Suppose the linear feedback gains in equation (6.16) * are fixed. Then the mismatch threshold M is maximum when the positive definite matrix Q - I (identity matrix). 208 Proof : Let the matrix Q in equation (6.60) be arbitrary positive 0 O O - ~ 3 ~ 0 0 th definite matrix and let 1 l/Amin(Q). Then Amin(Q) 1, where Q is 111 e equation (6.60). The solution of the equation (6.60) is given [18] by ~ AZC ~ Act P = I” e Q e dt (6.64) 0 And also consider T A A AP+PA=-I (6.65) c c with the solution A AZt Act P - Im e I e dt (6.66) 0 Now define ~ A ~ M - l/Amax(P) " A " (6.67) M l/Amax(P) Then from the equations (6.64) and (6.66) ~ a Act ~ Act P - P - [m e (Q - I) e dt (6.68) 0 However Q-IZO since all eigenvalues of the Q greater or equal to 1. Hence P - P 2 0, A (P) 2 A (P) and M 2 M from the equation (6.67). Since M max max is invariant under any choice of the positive scalar 1 by the Lemma 6.1” M 2 M holds for all Q. This completes the proof of the Theorem 6.2 209 Therefore, the mismatch threshold defined by the equation (6.29) can be modified as * A 1 M - 2A (P) (6.69) max by choosing Q - I. Other quantities involving Amin(Q) are modified by replacing Amin(Q) by 1. For example,the radius of the largest closed ball 5 defined by the equation (6.36) becomes (a1+ a + a e 2 3) 5 = - 2 4Amax(P) (M1+ M + M3) 2 6 . 4 Illustrative Example In order to demonstrate the robustness of the developed nonlinear state feedback controller, the composite arm, which has geometrical and material properties given in Table 6.1, is considered. It is clear from this table that the nominal system has no payload and 200 of lamina angle (fiber orientation) for each lamina. These two parameters will be considered as uncertain parameters in this example. The parameter variations of the proposed robotic system associated with the varying payload and the mismatched lamina angle are presented in Figures 6.1 and 6.2 respectively. It is observed from Figure 6.1 that the modal slope coefficient d¢1(o)/dx is very sensitive to a whole range of the payload, while the first mode natural frequency is sensitive to only a small range of the payload. And it is clear from Figure 6.2 that both the first and second mode natural frequencies are very sensitive to the lamina angle. Table 6.1 210 Geometrical and Material Properties of the Nominal Composite Arm Parameters and Properties Specifications (Graphite-Epoxy) Beam Length (L) 1.5 m Beam Width (b) 3.0 cm Beam Thickness (h) 1.8 mm Payload (mt) 0.0 kg Hub Moment of Inertia (1“) 0.00496 kg-m2 Number of Lamina 6 di - 0.3 mm Lamina Thickness (1 - 1. .6) 01 - 20° Fiber Orientation (1 - 1. .6) v; - 0.60 Fiber Volume Fraction (i - l, .,6) Efiber - 220 GPa Young's Modulus E . -4.27 GPa matrix -1750 k /m3 pfiber 5 Density 3 ”matrix.1200 kg/m v - 0.20 Poisson's Ratio fiber u - 0.34 matrix w1(rad/sec) w1(rad/sec) 25 13 01) I I I r I I I I I (15 Fi I I I I 1 .0 PAYLOAD(kg) gure 6.1 I I I 1.5 I I I 21) N O l U] l IIIII I I I I 60 r I I I LAMINA ANGLE(deg) Figure 6.2 90 211 20 mi 3; a 5... 0 0 w2(rad/sec) IFIIIIIIITrI .0 (l5 1 .0 PAYLOAD(kg) [W I 1.5 Parameter Variation of the Composite Arm with respect to the Payload I I 213 70 30~ 10 IIIII LAMINA ANGLE(deg) with respect to the Lamina Angle Parameter Variation of the Composite Arm 212 The system model parameters of the nominal system are presented in Table 6.2. For the determination of these parameters, two-node element was employed with two degrees of freedom at each node, and.the number of elements was ten. In this simulation, the rigid-body mode and the first flexible mode are considered as the primary modes. Hence the system matrices in equation (6.6) become (6.70a) (6.70b) _ l 0 0 0 0 0 A0- 0 1 _ -539.9 -O.3253‘ . 0 - 10.2135 B0- 0 _30.7323_ Table 6.2 System Model Parameters of the Nominal System d¢i Mode wi(Hz) E);— (0) ¢i(L) 1 3.698 3.009 -0.960 2 8.603 2.798 0.993 3 19.800 1.139 -1.156 4 37.931 0.566 1.180 5 62.484 0.338 -1.188 213 -0 0 - AA(p) - 0 0 (6.70c) 0 0 -0 0 'p11 'p213 a p p = 9 p11 539.9wr1(wr1+ 2), p21 0.32536r1 [0 0 0 ' 0 AA(£) - (6.70d) 0 _0 0 -211 -221_ 2 2 2 £11- 539.9wr1(wr1+ 2), £21= 0.325360r1 [0 0 0‘ AA(p,2) - 0 0 0 (6.70e) 0 0 L 0 -p21 0, _ 9 p21 1079.9wr1wr1 - 0 - p . AB(p) - 10°2135br0 (6.70f) 0 30.7323bp - r1- AB(2) - 0 (6.70g) It is noted that the damping ratio of the first flexible mode was chosen by 0.007. And the input uncertain matrix AB(2) was chosen by 0, since the variation of the lamina angle does not affect the variation of the control influence coefficient. 214 Now the decomposition of the system uncertain matrices can be accomplished in view of the equation (6.14) as follows. - _ p P _ p D(p) [ 0 0 15.82wr1(wr1+2) 0.00953er] 2 2 2 0(2) - [ 0 0 -15.82wr1(wr1+2) -0.00953er] (6.71) D(p 2) — [ 0 0 -31 64wp wi 0] ’ ' r1 r1 _ p p E(p) 0.0994615r0 + 0.9br1 The feedback gains K in equation (6.15) were chosen by K - [ -35 -5 ~50 -4 ] (6.72) Hence with the positive definite matrix Q-I, the solution of the Lyapunov equation (6.18) was obtained as ' 8.683 0.730 1.295 -0.242 ‘ P _ 0 730 0.189 1.397 -0.0548 (6.73) 1.295 1.397 27.133 -0.343 _-O.242 ~0.0548 -0.343 0.0194‘ The variation of the control influence coefficient was assumed to be bounded with -0.06 s bios 0, and -0.05 S bgls 0 so that the condition (6.22) is satisfied. Then the p(x), which is nonlinear portion of the controller, is obtained as follows in view of the equations (6.23) and (6.26) p(x) - - + p2(x) + p3(x)) sgn(BgPX). (6.74a) if |0.00143x + 0.2461x + 3.7167x + 0.0376x4I > e l 2 3 or 215 p(x) - -(p1(x) + p2(x) + p3(x))(0.00143x1+ 0.2461x2 + 3.7167X3+0.0376X4)/6, (6.74b) if I0.00143x1+ 0.2461x2+ 3.7167x3 + 0.0376x4I s e where p p 2 pl(x) - msx (1.05374 (15.82(2.0|wr1| + (wrl) )Ix3| wrl + 0.051(35|x1| + 5Ix2| + 50Ix3| + 4|x4l))} (x) - max {15 82(2 0|m£ | + (m2 )2)|x | (6 75) 92 ' ° r1 r1 3 ' 2 r1 w + 0 00953lw£ | |x I) ' r1 4 2 p3(x) - mgx {31.64(|w:1| Iwr1|)|x3l} wil r1 0: and sgn(Bng) is defined by [ 1, if Bng > 0 T T sgn(8o P x) - < 0, if BOPx - 0 (6.76) L-1, if Bng < 0 Figure 6.3 presents the open-loop tip displacement of the nominal system. The vibration of the first flexible mode was invoked by taking the initial conditions as T x(O) - { 0 0 0.06 0} (6.77) It is clear from the figure that the amplitude decays with only assumed damping ratio of the first mode. 216 TIP DISPLACEMENT(cm) 53 O L _8.0 r I I f I I I T 0.0 1.0 2.0 3.0 4.0 TlME(sec) I T I i I r I Figure 6.3 Open-Loop Tip displacement of the Nominal System Figure 6.4 presents the simulated step responses for the commanded tip position of 15 cm with the linear and nonlinear controller. Herein the linear controller is given by u - Kx, and the nonlinear controller is defined by equation (6.16). The typical value of e - 0.003 was chosen, and the actual bound on the variation of the uncertain parameters were chosen by -0.2 5 6:1 5 0 (6.78) -0.4 s w! .<_ 0 P The variation of wrl corresponds to the variation of payload from 0 kg to about 0.17 kg. It is noted that the mass of the flexible arm is 0.124 kg. 217 . 2 . . 0n the other hand, the variation of wrl corresponds to the variation of o o the lamina angle from 20 to approximately 48 . For the larger variation . 2 . of the uncertain parameters, both wgl and wrl were taken as negative. The situation of the positive mil and wfil, which correspond to an increment of the fundamental natural frequency is not considered in this study. In fact, an increment of the natural frequency due to the payload can not be happened in practice. Furthermore, choosing the positive of w makes the system be undergone with the smaller perturbation due to rl P 2 r1 and positive wrl' the compensation between the negative w It is observed from the tip displacement response in Figure 6.4 that the response with the nonlinear controller settles in the commanded position in 5.8 seconds (settling-time) inducing maximum overshoot of 9 cm, while the response with the linear controller settles in the commanded position in 12 seconds (settling-time) inducing maximum overshoot of 13.8 cm. This shows that the nonlinear controller have a much better capability for regulation of the response than the linear controller in the presence of system uncertainties. This can be also interpreted that the nonlinear controller can tolerates larger system parameter variations than the linear controller. It is also interesting to observe from the input torque response that the required input torque of the nonlinear controller is about 60 percent that of the linear controller for the simulated response time. Table 6.3 presents a comparison of the tip displacement response (settling-time and maximum overshoot) between the linear and nonlinear INPUT TORQUE(N.m) TIP DlSPLACEMENT(cm) 218 —- NONUNEAR UNEAR N .p. l l L l l (I) l .1 .0. I I _6 I I I r I I 0.0 3.0 6.0 TlME(sec) 12.0 —- NONUNEAR UNEAR I— I "0.3 I I r I T 1* 0.0 3.0 6.0 TlME(sec) 9.0 12.0 Figure 6.4 Comparison of Step Responses with the Linear and the Nonlinear Controller ANGULAR VELOCITY(rad/sec) 219 ANGULAR POSITION(deg) IJNEAR -— NONUNEM? 0.. o 0 on o ...-... ‘00....- 3.0 I I I l I . I ' 6.0 9.0 ' 12.0 TlME(sec) .0 O l l l l l l l l l l .0 01 IJNEAR -— NONUNEAR .0 0 Figure 6.4 I I I I I I 61) 9&) 121) TlME(sec) (Continued) 220 Table 6.3 Comparison of the Tip Displacement Response (Settling Time and Maximum Overshoot) Variation Settling Time (sec) Maximum Overshoot (cm) Bound Nonlinear Linear Nonlinear Linear 451- 0 2.9 5.8 7.1 11.3 -o.isw1r’lso 4.5 8.6 8.2 12.6 025651; 5.8 12.0 9.4 13.8 03545150 11.6 22.1 10.8 15.5 $456550 24.3 44.5 12.4 17.0 controller under various bounds of oil. The variation of the lamina angle was fixed as -0.4 S mils 0. It is clear from the table that the settling- time with the nonlinear controller is half of that with the linear controller in all variation bounds. The bound -0.3 S ”:15 0 corresponds to the varying payload from 0 kg to about 2.0 kg which is about 16 times of the arm mass. Similar results are anticipated from tflua'various . . 2 . . . . ‘variations of the wrl With the fixed variation of the wgl. Therefore, it is possible to adjust the variation bound of each uncertain parameter for achieving the required system performances. Next, in order to observe the sensitiveness of the controller parameter' c in the equation (6.17), the different value of e - 0.1 was 221 36 - c=0.1 ‘ — e=0.003 30“ E - 3 I.— 2 Lu 2 Lu Q 1’) O. Q Q E }_ '—'6 I I I I I I I I I I 0.0 3.0 6.0 9.0 12.0 TlME(sec) 0.6 =04 - -— e=0.003 ’23 _ E; 0.3 LL] . :3 .. 0 C! O I— f— D 0. ...... 11 .2. "0.3 I l I I l I I I I I 0.0 3.0 6.0 9.0 12.0 TlME(sec) Figure 6.5 Comparison of Step Responses with 6-0.1 and 6-0.003 222 chosen. Figure 6.5 presents the comparison of the step responses with two different values of c. It is clear from this figure that the smaller 6 produces the better system responses. This fact is obvious from the largest ball I; defined by the equation (6.36). For stabilization problem, it is desirable to have I). as small as possible, hence increasing the robustness of the system stability. As 6 4 O, n 4 0. However, the smaller 6 may produce the larger unwanted vibration in practice due to the larger chattered input torque as shown in Figure 6.5. It is finally noted that even though two uncertain parameters;the varying payload and the mismatched lamina angle were considered in this study, the approach developed herein can be extended to the case which accounts for many uncertain parameters. The manner of choosing uncertain parameters should be dependent on the physical environments conditions. For example, if the robotic system is supposed to operate under high temperature or moisture environment, these uncertain parameters should be taken into consideration. 6.5 Sui-nary of the Chapter The nonlinear robust controller for the composite arm was designed by utilizing the properties of the uniform and uniform ultimate boundedness of the solution. The varying payload and the lamina angle were chosen as uncertain parameters. It has been demonstrated that the nonlinear controller had much better capability for regulation of the response than the linear controller in the presence of the uncertain parameters. This implies that the larger tolerance of the fabrication error (lamina angle) of the composite arm can be allowed for the nonlinear control system, and 223 also the larger bound of the unknown payload is allowed with guaranteed system stability. Even though the proposed nonlinear controller was not experimentally implemented in this study, the realization of this controller would accelerate the relevant researdh area in the future. The extension of this controller to flexible multi-link robotic systems, and other flexible structure systems fabricated from composite materials are remained as further research areas as well as the experimental realization of this controller. CHAPTER'VII CONCLUSIONS AND RECOMMENDATIONS 7.1 Conclusions of the Thesis Three different control schemes for a single-link flexibLe manipulator fabricated from composite laminates or smart materials incorporating electro-rheological (ER) fluids were developed in this thesis ; the output feedback control, the hybrid control, and the nonlinear state feedback control robust to system parameter variations. The output feedback controller featuring two colocated angular position and velocity sensors was designed and experimentally implemented in order to investigate and compare the control performance of the flexible manipulator fabricated from aluminum (aluminum arm) and composite laminates (composite arm). The controller design was accomplished by utilizing the root-locus technique and the Lyapunov method. In the root-locus analysis, the effect of the compensator zero was investigated by considering various location of the compensator zero. When aLlIunhzero open-loop poles and zeros were located left to the ,given compensatory the system exhibited a sluggish motion without inducing overshoot. However, when the compensator zero was positioned far left to all open-loop poles and zeros, the system exhibited a fast motion with tradeoff the overshoot. It has been demonstrated through the analytical and experimental 224 225 investigation that the composite arm had superior system performances such as faster settling-time, smaller input torque, smaller overshoot and a superior increment of stability relative to the aluminum arm. It was obtained from the measured step responses that with certain position and velocity feedback gains the delay time was 0.26 second for the aluminum arm and 0.21 second for the composite arm, and the maximum overshoot was 1.02 degree for the aluminum arm and 0.3 degree for the composite arm. It was also observed that the reduction of energy consumption by 40-50 percent for the commanded motion, and the reduction of the tip deflection up to 50 percent was accomplished by employing the composite arm. These experimental results clearly indicate the significant payoffs associated with implementing composite materials featuring superior strength- and stiffness-to-weight ratios to the commercial metals on the flexible manipulators. It has been also shown that the implemented output feedback controller was very sensitive to the external disturbance which has a forcing frequency close to the cantilevered-mode natural frequency of the arm, regardless of the aluminum and composite arm. This attempt was undertaken herein to introduce a single-link flexible manipulator fabricated from smart materials incorporating electro-rheological fluids (smart ER arm) to obtain the robustness of the output feedback controller to the imposed disturbance with the help of distributed parameter actuator, which is characterized by embedded electro-rheological (ER) fluids in the structure. As a sequel, the hybrid controller featuring the output feedback controller employed for the aluminum and composite arm, and the pseudo- state feedback controller associated with the distributed ER fluid 226 actuator was proposed and implemented on the smart ER arm in the Chapter V. Before developing the hybrid control scheme, the phenomenological equation of motion of the smart robotic system has been developed based on the experimental observations of the modal characteristics of the smart beam structure. Then, the hybrid controller was formulated based on the phenomenological equation of motion and the empirical dynamic model of the ER fluid actuator. The pseudo-control force (input voltage) associated with the pseudo-state feedback controller accounts for the change in the global damping and stiffness properties of the smart robotic arm due to electrical potentials imposed on the ER fluid domains. The stability criterion has been established for the robotic system by considering the changed properties as the plant perturbation due to electrical potentials. Moreover, this stability criterion has been modified by incorporating the empirical dynamic model of the ER fluid actuator. The modified criterion says that the stable closed-loop system with the output feedback controller remains as stable for positive input voltage considered in this investigation. It has been shown through the experimental implementation that the robustness of the employed conventional output feedback controller was significantly enhanced by employing the proposed pseudo-state feedback controller in the presence of the disturbance. The nonlinear robust state feedback controller for the composite arm was developed for the first time to overcome uncertain system parameters associated with the varying payload encountered in the operation and the mismatched lamina angle encountered in the manufacturing process. The controller design has been accomplished by utilizing the properties of uniform and uniform ultimate boundedness of the solution. The major 227 difference of this approach from the currently existing approaches was that the effect of each uncertain parameter as well as the effect of the interaction of the uncertain parameters to the total system performance was considered in this study, hence providing a certain measure of the magnitude of each uncertain parameter. It has been demonstrated from the simulation that the proposed nonlinear controller could tolerate larger system parameter variations than the conventional linear controller, showing much faster settling-time and smaller overshoot with the nonlinear controller. It has been also shown that the system performance was very sensitive to control parameters such as e. The author believes that the research prosecuted in this thesis will provide significant impact on the diverse marketplace such as automobile industry and aerospace industry. And also, this complementary theoretical, computational and experimental investigation on the control methodologies for the flexible manipulators fabricated from composite laminates or smart materials incorporating ER fluids will furnish fundamental and important guidelines, and precious information for the development of the next generation robotic system characterized by high- speed of operation, high end-point positional accuracy, and also high adaptability to unstructured environments. 7 . 2 Reco-endations for Future Work As observed in Chapter V of this thesis, the proposed hybrid controller was implemented on the smart robot arm by applying the pseudo- control force (input voltage) in the fashion of the open-loop control. Hence, the information of the disturbance such as a forcing frequency should be identified prior to applying the proper input voltage. However, 228 it is possible to construct a closed-loop control system of the pseudo- state controller by employing appropriate output sensors to obtain the information of the imposed disturbance, and by identifying the system parameters as a function of the input voltage in the frequency domain. From the disturbance information and identified system parameters, the appropriate pseudo-control force can be obtained and supplied from the microprocessor to the embedded ER domains in real time. Thisiclosed-loop system will highlight the distinct merit of employing the smart robotic system which can respond to unstructured environmental conditions. Another potential future work is to develop a flexible robotic system featuring a high level of hybrid controller associated with conventional actuators and various type of actuators incorporating smart materials. The high level of controller consists of more than three actuators which provide own individual actuating force for the different purpose and also in a different manner. Since each smart material such as ER fluids, piezoelectric materials and shape memory alloys has its own distinct characteristic as observed in the literature review in Chapter i1, it is possible to construct the high level of controller which can adaptively respond to highly unstructured environments. As an example, consider a flexible robotic system featuring a conventional actuator, an ER fluid actuator and a piezoelectric actuator. When the system needs to be settled very fast, then the control force associated with the conventional actuator and the piezoelectric actuator are supposed to be activated to increase damping properties. However, if the system is subjected to a disturbance which has a forcing frequency close to the natural frequency of the robot arm, the control force associated with the ER actuator and the conventional actuator are supposed to be activated in 229 order to avoid the resonance phenomenon of the system. This high level of hybrid controller will provide tremendous advantages over the conventional controller in real world. The potential advantages include the self-controllability of the plant, the robustness of the controller to the uncertainties such as external and internal disturbances, the flexibility of the trajectory or work space, and the flexibility of all control parameters. In parallel with above mentioned future work, it may be useful to develop a hierarchical control methodology for the smart flexible structures. The principal idea behind the hierarchical control methodology is that the overall control system is divided into more than two sub-control system ; for example, global controller, intermediate controllers and regional controllers in the flexible structure. By considering each finite segment actuator incorporating ER fluids or piezoelectric materials as a regional controller, sum of several segment actuators as an intermediate controller, and sum of all intermediate controllers as a global controller, the total system can be controlled for the desired system performances. This control scheme will be very useful especially for the very large flexible structures such as space station. Since there always exist wave propagation delay during the sensing of the outputs in the large flexible structures, the system instability may be arisen.due to the observation spillover caused by the delay. This observationxspillover problem can be eliminated by employing the hierarchical control methodology. The optimal placement of the appropriate actuator associated with each regional controller is one of the problem to be solved in the process of formulating the hierarchical controller. 230 Finally, the extension of the nonlinear robust state feedback controller developed in Chapter VI to the multi-links flexible manipulators, and experimental implementation of this controller are remained as a future work. APPENDIX APPENDIX List of Publications Thesis 1. S.B. Choi, " On Robust Tracking in Uncertain Systems - A Variable Structure Approach," Master Thesis, Department of Mechanical Engineering, Michigan State University, June 1986. S.B. Choi, " Control of Single-Link Flexible Manipulators Fabricated from Advanced Composite Laminates and Smart Materials Incorporating Electra-Rheological Fluids," Ph D Thesis, MEL/IMSL, Department of Mechanical Engineering, Michigan State University, June 1990. Refereed Journal Elications 1. M.V. Gandhi, B.S. Thompson, S.B. Choi and S. Shakir, " Electro- Rheological-Fluid-Based Articulating Robotic Systems, " ASME Jougnal 9f Meghanisns, Transmissions, and Automation in Design, Vol.111, No. 3, 1989, pp. 328-336. S.B. Choi, M.V. Gandhi, B.S. Thompson and C.Y. Lee, " An Experimental Investigation of an Articulating Robotic Manipulator with a Graphite- Epoxy Composite Arm," ournal of Robotic Systems, Vol.5, No.1, 1988, pp. 73-79. M.V. Gandhi, B.S. Thompson and S.B. Choi, " A New Generation of Innovative Ultra-Advanced Intelligent Composite Materials Featuring Electra-Rheological Fluids: An Experimental Investigation," Jouznal 9f Qonngsige Materials, Vol.23, 1989, pp. 1232-1255. M.V. Gandhi, B.S. Thompson and S.B. Choi, " A Proof-of-Concept Experimental Investigation of a Slider-Crank Mechanism Featuring a Smart Dynamically-Tunable Connecting-Rod Incorporating Embedded Electra-Rheological Fluid Domains," Journal of Sound and Vibration, Vol. 135, No.3, 1989, pp. 511-515. S.B. Choi, C.Y. Lee, B.S. Thompson and M.V. Gandhi, " Elastodynamic Characteristics of Smart Beams Incorporating Electra-Rheological F1uids.," AIAA Journal of Guidance, Control, and Dynamics (under review). 231 10. ll. 12. 13. 14. 232 S.B. Choi, A. Magolan, B.S. Thompson and M.V. Gandhi, " An Experimental and ,Theoretical Investigation of the Static and Elastodynamic Responses of an Industrial Robotic Manipulator Featuring a Graphite-Epoxy Composite Arm," ournal of Robotic Systems (under review). S.B. Choi and R.L. Tummala, " Modeling and Control of a Flexible Single-Link Robotic Arm Fabricated from Composite Materials," lEEE Iransgggions on Robotics and Automation (under review). S.B. Choi, B.S. Thompson and M.V. Gandhi, " Smart Structures Incorporating Electra-Rheological Fluids for Vibration-Control and Active-Damping Applications: An Experimental Investigation," M Joninnl 9f Vibration, Acousticfis, Stress and Reliability in Design (under review) . S.B. Choi, B.S. Thompson and M.V. Gandhi, " Modeling and Position Control of a Single-Link Flexible Robotic Manipulator Fabricated from Advanced Composite Materials," Journal of Robotic Systems (under review). S.B. Choi, B.S. Thompson and M.V. Gandhi, " Modeling and Output Feedback Control of a Single-Link Flexible Robotic Manipulator Featuring a Graphite-Epoxy Composite Arm,” W Systens, Measnrenent, and Control (under review). " A Nonlinear Control Robust to System Parameter Variations with Application to a Single-Link Flexible Manipulator Fabricated from Composite Laminates," (in preparation). " Modeling and Control of Smart Flexible Structures Incorporating Piezoelectric Materials," (in preparation). " Active-Vibration-Tuning of a Simply-Supported Beam Featuring Electra-Rheological Fluids,“ (in preparation). " Modeling and Control of a Single-Link Flexible Robotic Manipulator Fabricated from Smart Materials Incorporating Electra-Rheological Fluids, " (in preparation) . Confergngg Elications 1. S.B. Choi and S. Jayasuriya, " A Sliding Mode Controller Incorporating Matching Conditions Applied to Manipulators," Pinceedings oi the 10th IFAC World Congress, Munich, West Germany, July 1987. S. Jayasuriya and S.B. Choi, " On the Sufficiency Condition for Existence of a Sliding Mode," 0 eed o Ame a o Coniegence, Minneapolis, Minn., 1987, Vol.1, pp. 84-89. 10. 11. 233 S.B. Choi, F. Sun, B.S. Thompson and M.V. Gandhi, " An Experimental Investigation of the Static and Dynamic Responses of an Industrial Robotic Manipulator," Proceedings of 1982 ASME Design Automation Conferance, Boston, Mass., Sept. 1987, ASME Paper No. 87-DAC-50, DE- Vol.10-2, pp. 421-427. M.V. Gandhi, B.S. Thompson, S.B. Choi and S. Shakir, " Electro- Rheological-Fluid-Based Articulating Robotic System,” Proceedings of l982 ASME Design Automation Conference, Boston, Mass. , Sept. 1987, ASME Paper No. 87-DAC-55, DE-Vol.10-2, pp. 1-10. M.V. Gandhi, B.S. Thompson and S.B. Choi, " Smart Ultra-Advanced Composite Materials Incorporating Electra-Rheological Fluids for Military, Aerospace, and Advanced Manufacturing Applications," Proceedings 9f ghe Bid Annual ASM/ESD Advanced Composiga Confeijence and Expositign, Detroit, MI, Sept. 1987, pp. 23-30. S.B. Choi, B.S. Thompson and M.V. Gandhi, ” Electra-Rheological Fluids Technology Stimulates a New Generation of Robotic and Machine Systems," Pioceedings of the 10th Annlied Mechanism Conference, New Orleans, La., December 1987, Vol.1, pp. 3B.1-381.8. S.B. Choi, M.V. Gandhi, B.S. Thompson and C.Y. Lee, " An Experimental Investigation of an Articulating Robotic Manipulator with a Graphite- Epoxy Composite Arm," Pioceedings 9f the l0§h Annliea Mechanism Conteranae, New Orleans, La., December 1987, Vol.1, pp. 3B2.1-3B2.8. M.V. Gandhi, B.S. Thompson and S.B. Choi, " Ultra-Advanced Composite Materials Incorporating Electra-Rheological Fluids," Prageedings 9i ghe 4th Japan-US Composite Conference, Washington D.C., June 1988, pp. 875-884. M.V. Gandhi, B.S. Thompson and S.B. Choi, " Smart Materials Incorporating Electra-Rheological Fluids: A Theoretical and Experimental Investigation," Proceedings of the 18th Annual Meeting of ghe Eina Particle Society, 1987 Particle/Power Technology, Microcontamination and Biotechnology Exhibition and Technology Forum. S.B. Choi, B.S. Thompson and M.V. Gandhi, " An Experimental Investigation on the Active-Damping Characteristics of a Class of Ultra-Advanced Intelligent Composite Materials Featuring Electro- Rheological Fluids," Proceedings of the Damping '89 Conferenge, West Palm Beach, Fl. , February 1989, Organized by the Flight Dynamics Laboratory of the Air Force Wright Aeronautical Laboratories, Wright Patterson Air Bases, Ohio. S.B. Choi, B.S. Thompson and M.V. Gandhi, " An Active Vibration- Tuning Methodology for Smart Flexible Structures Incorporating Electra-Rheological Fluids: A Proof-of-Concept Investigation,” Proceedings of 1989 American Control Confegenae, Pittsburgh, Pa. , June 1989. 12. l3. 14. 15. 16. 17. 18. 234 S.B. Choi, B.S. Thompson and M.V. Gandhi, " Smart Structures Incorporating Electra-Rheological Fluids for Vibration-Control and Active-Damping Applications: An Experimental Investigation," roceedi s of the t B ennial ASME Con ere ce 0 echa ical Vibration and Noise, Montreal, Canada, Sept. 1989. S.B. Choi, A. Magolan, B.S. Thompson and M.V. Gandhi, ” An Experimental and Theoretical Investigation of the Static and Elastodynamic Responses of an Industrial Robotic Manipulator Featuring a Graphite-Epoxy Composite Arm,“ P oceed n s o the lst Natipnal Qpnference on Applied Mechanisms ana Robotics, Cincinnati, Ohio, Nov. 1989, Vol.2, pp. 8B5.1-8B5.11. S.B. Choi, L.P. Chao, B.S. Thompson and M.V. Gandhi, " An Integrated Design Strategy for the Optimal Structural Fabrication and the Optimal Control of a Single-Link Flexible Manipulator Fabricated from Advanced Composite Laminates," Proceedings of the 1st National Conference an Applied Mechanisms and Robotics, Cincinnati, Ohio, Nov. 1989, Vol.2, pp. 8B4.1-8B4.10. S.B. Choi, M.V. Gandhi and B.S. Thompson, " An Experimental Investigation of the Elastodynamic Responses of a Slider-Crank Mechanism Featuring a Smart Connecting-Rod Incorporating Embedded Electra-Rheological Fluid Domains," Proceedings of the lst National Conference on Applied Mechanism and Robotica, Cincinnati, Ohio, Nov. 1989, Vol.2, pp. 9B3.1-9B3.7. C.Y. Lee, S.B. Choi, B.S. Thompson and M.V. Gandhi, " A Variational Formulation for the Finite Element Analysis and Control of Linkage and Robotics Systems Featuring Smart Materials Incorporating Piezoelectric Materials," Proceedings of the lst National Conference an Applieg Mechanisms and Robotics, Cincinnati, Ohio, Nov. 1989, Vol.1, pp. 2C1.1-2C1.10. S.B. Choi, B.S. Thompson and M.V. Gandhi, " A Theoretical and Experimental Investigation on the Control of a Single-Link Flexible Robotic Manipulator Fabricated from Composite Materials," to ha ppesanpaa a; tha 1990 American Control Conference, San Diego, CA, May 1990. S.B. Choi, B.S. Thompson and M.V. Gandhi, " Modeling and Feedback Control of a Single-Link Flexible Robotic Manipulator Featuring a Graphite-Epoxy Composite Arm," to be presented at the l990 lEEE lntepnapional Conference on Robotics and Autonation, Cincinnati, Ohio , May 1990. Technical Mr}; 1. S. Shakir, L.P. Chao, S.B. Choi, M.V. Gandhi and B.S. Thompson, " Smart Materials Incorporating Electra-Rheological Fluids Technology for Articulating Systems and Advanced Structures," MSU Machinery Elaapodynanias Lab, Report, Report No. 870108801, Feb. 1987. 2. 235 S.B. Choi, B.S. Thompson and M.V. Gandhi, " The Synthesis of Smart Materials Incorporating Electra-Rheological Fluids for Military, Aerospace and Advanced Manufacturing Applications: An Interim Report," MSQ Maahinegy Elaagpaynaniaa Lab, gepopt t9 Laid Winn, NC, Report No. 87082501, July 1987. S.B. Choi, B.S. Thompson and M.V. Gandhi, " The Synthesis of Smart Materials Incorporating Electra-Rheological Fluids for Military, Aerospace and Advanced Manufacturing Applications: An Interim Report," U Ma iner Elastod n m cs ab e art to Lo (1 Winn, NC, Report No. 87082502, December 1987. S.B. Choi, B.S. Thompson and M.V. Gandhi, " The Synthesis of Smart Materials Incorporating Electra-Rheological Fluids for Military, Aerospace and Advanced Manufacturing Applications: A Final Report," SU a e E a tod am cs a e o to Lord Cor oratio , NC, Report No. 87082503, March 1988. BIBLIOGRAPHY 10. 11. BIBLIOGRAPHY Abou-Hanna,J.J. and Evces,C.R. , " Dynamics and Control of Flexible Manipulators," Proceedinga pi ASME Winger Annual Meeting, Boston, Mass., DSC-Vol.6, 1987, pp. 269-276. Adriani, P.M. and Cast, A.P., " A Microscopic Model of Electro- Rheology," Janinal 9f Phyaics ',Plnida, Vol.31, No.10, 1988, pp. 2757-2768. Arbe1,A. and Gupta,N.K., " Robust Colocated Control for Large Flexible Space Structures," J urna 0 Cu da e Con d Mina, Vol.4, No.5, 1981, pp. 480-486. Archer, J.S., " Consistent Mass Matrix for Distributed Mass Systems," t ca e of v e Vol.89, 1963, pp.161-171. Asada,H. , Ma,Z.D. and Tokumaru,H., " Inverse Dynamics of Flexible Robot Arms for Trajectory Control," Proceedings of ASME Winter Annnal Meeping, Boston, Mass., DSC-Vol.6, 1987, pp. 329-336. Bailey, T. and Hubbard, J.E., " Distributed Piezoelectric-Polymer Active Vibration Control of a Cantilever Beam," 1AA Jougnal pf Guidanae and Control, Sept. 1985, Vol. 8, No. 5, pp. 605-611. Balas,M.J. , " Modal Control of Certain Flexible Dynamic Systems," WWW. Vol-16. No.13. 1978. pp- 450-462. Balas,M.J., " Feedback Control of Flexible Systems," lEEE Iransaationa on Augomagic Control, Vol.AC-23, No.4, 1978, pp. 674- 679. Balas,M.J., " Direct Velocity Feedback Control of Large Space Structures," I Jo al Guida ce Co trol and D amic , Vol.2, No.3, 1979, pp. 252-253. Balas,M.J., " Trends in Large Space Structure Control Theory: Fondest Hopes, Wildest Dreams," lEEE Tranaactiona on Automatic mainlflol. AC-27, No.3, 1982, pp. 522-535. Balas, M.J., Kaufman, H. and Wen, J., " Stable Direct Adaptive Control of Linear Infinite-Dimensional Systems Using a Command Generator Tracker Approach." WW WW. San Diego. CA. June 1984. 236 12 13. 14. 15. 16. 17. 18. 19. 20. 21. 22. 23. 24. 25. 237 Bar-kana,1. and Kaufman, H. " Some Applications of Direct Adaptive Control to large Structural Systems," MAW W. Vol 7 No 6 1984 PP 717- 724 Barmish,B.R. and Leitmann,G., " On Ultimate Boundedness Control of Uncertain Systems in the Absence of Matching Assumptions," IEEE WWW. Vol. AC-27. No.1. 1982. PP- 153- 158. Barraco,A., Cuny,B. and Ishiomin,G., " Dynamic Models for Flexible Robots: Different Approaches," Proceedings of international Qpnference an Robotics and Automation, 1986, pp. 1038-1043. Baruh,H. and Silverberg,L. , " Robust Natural Control of Distributed Systems," AIAA Jpnmal oi Guidance. .Control. and Dynamics, Vol.8, No.6, 1985, pp. 717-724. Bayo,E. , " A Finite-Element Approach to Control the End-Point Motion of a Single-Link Flexible Robot," goninal oi Robotic Syaiema, Vol.4, No.1, 1987, pp. 63-75. Baz,A. and Poh, S. , " Performance of an Active Control System with Piezoelectric Actuators, " na 0 and Vibration, Vol . 126 , No. 2, 1988, pp.327-343. Bellman,R., Siabiliiy Iheog oi Qiifeieniial Pguapipns, Dover, New York, 1969. ' Berge,C., Wei, Oliver and Boyd, London, England, 1963. Block, H. and Kelly, J.P., " Electro-Rheology," Journal oi Physics D; Appliad Phyaica, Vol.21, 1988, pp. 1661-1677. Bodden,D.S. and Junkins,J.L. , " Eigenvalue Optimaization Algorithms for Structure/Controller Design Iterations,” AIAA gpurnal of WW. Vol.8. No.6. 1985. PP-697-706- Book,W.J. , " Analysis of M‘assless Elastic Chains with Servo Controlled Joints," D c s ems ea e t Maggi, Vol.101, 1979, pp. 187-192. Book,W.J., " Recursive Lagrangian Dynamics of Flexible Manipulator Arms," ur Robo cs Resea ch, Vol.3, No.3, 1984, pp. 87-101. Book,W.J. , Alberts,T.E. and Hastings,G.G. , " Design Strategies for High Speed Lightweight Robots," m e chan ca En e n Vol.26, 1986, pp. 26-33. Book,W.J. and Cetinkunt,S., " Near Optimum Control of Flexible Robot Arms on Fixed Paths," c s f 4 C fe ence 0 De 0 d Qpnpzol, Ft.Lauderdale, FL, 1985, pp. 1522-1528. 26. 27. 28. 29. 30. 31. 32. 33. 34. 35. 36. 37. 38. 238 Book,W.J. and Majette,M. , " Controller Design for Flexible, Distributed Parameter Mechanical Arms via Combined State Space and Frequency Domain Techniques," ASME Jouinal of Dynamic Systems, Maasuiemeni, and Qonipgl, Vol.105, 1983, pp. 245-254. Book,W.J. , Maizza-Neto,0. and Whitney,D.E. ," Feedback Control of Two Beam, Two Joint Systems with Distributed Flexibility," ASME Janina], pi Dynania Syaiens, Maasuieneni, ana Control, Vol.97, No.4, 1975, pp. 424- 431. Brogan, W.L. , Madein Qantrol Ihaory, Quantum Press, New York, 1974. Bullough, W.A. and Foxon, M.B., " A Proportionate Coulomb and Viscously Damped Isolation System," 0 a ound a d V bratio Vol.56, No.1, 1978, pp. 35-44. Burke, S.E. and Hubbard, J.E., ” Active Vibration Control of a Simply Supported Beam Using a Spatially Distributed Actuator," lEEE Qpnpiol Syaiem Magazine, Aug. 1987, pp. 25-30. Cannon,R.H. and Schmitz,E., " Initial Experiments on the End-Point Control of a Flexible One-Link Robot," The International Journal of protica Researah, Vol.3, No.3, 1984, pp. 62-75 Chalhoub,N.G. and Ulsoy,A.G. , " Control of a Flexible Robot Arm: Experimental and Theoretical Results,” SM ou na na W. Vol.109. 1987. PP- 299-309- Chang,L.W. , " Dynamic Analysis of Robotic Manipulators with Flexible Links," W, Purdue University, 1984. Chassiakos,A.G. and Bekey,G.A., " On the Modeling and Control of a Flexible Manipulator Arm by Point Actuators," MM aniaianca an Degiaion and Contipl, Athens, Greece, December 1986, pp. 1145-1150. Chen, C.T., Linea; Sysian Ihaoiy ang Dasign, Holt, Reinehart and Winston Co., New York, 1984. Chen, Y.H., " On the Robustness of Mismatched Uncertain Dynamical Systems," am c te s easu emen Vol.109, 1987, pp. 29-35. Choi, S.B., Chao, L.P., Thompson, B.S. and Gandhi, M.V., " An Integrated Design Strategy for the Optimal Structural Fabrication and the Optimal Control of a Single-Link Flexible Manipulator Fabricated from Advanced Composite Laminates," Pioceedinga oi the e ed ha ms 0 c Cincinnati, Ohio, Nov. 1989, Vol.2, pp. 8B4.1-8B4.10. Choi, S.B. , Gandhi, M.V. , Thompson, B.S. and Lee, C.Y. , "An Experimental Investigation of an Articulating Robotic Manipulator with e Graphite-EPOXY Arm." Wm. Vol. 5. 1988, pp.73-79. 39. 40. 41. 42. 43. 44. 45. 46. 47. 48. 49. 239 Choi, S.B. , Thompson, B.S. and Gandhi, M.V., " Electra-Rheological Fluids Technology Stimulates a New Generation of Robotic and Machine Systems," e A ec m C ere ce, New Orleans, La., Dec. 1987, Vol.1, pp. 38.1-3Bl.8. Choi, S.B. , Thompson, B.S. and Gandhi, M.V., An Active Vibration- Tuning Methodology for Smart Flexible Structures Incorporating Electro-Rheological Fluids: A Proof-of-Concept Investigation," Pippappinga pi 1,982 MEIIEQB Conirol Qpniaience, Pittsburgh, Pa., June 1989. Choi, S.B., Thompson, B.S. and Gandhi, M.V., " Smart Structures Incorporating Electra-Rheological Fluids for Vibration-Control and Active-Damping Applications: An Experimental Investigation, " , . ----I: . me a ;-e _a a M Conference ._ ec W, Montreal, Canada, Sept. 1989. Claus, R.O., Mckeeman, J.C., May, R.G. and Bennett, K.D., " Optimal Fiber Sensors and Signal Processing for Smart Materials and Structures Applications," In esearch Of i W s Snap; Maiepiala, Stznppniea and Maihematical lssnea, Blacksburg, VA, 1988. Coddington, E.A. and Levinson, N., Theor of Ordinar Dif e e a1 Minna, McGraw-Hill, New York, 1955. Corless, M.J. and Leitmann, G., " Continuous State Feedback Guaranteeing Uniform Ultimate Boundedness for Uncertain Dynamic Systems , " lEPE Iransactiona on Automatic Control , Vol . AC- 26 , No . 5 , 1981, pp. 1139-1144. Crawley, E.F. and de Luis, J., " Use of Piezoelectric Actuators as Elements of Intelligent Structures," AIAA Jouinal, Vol. 25, No.10, 1987, pp. 1373-1385. Dado,M. and Soni,A.H. , " A Generalized Approach for Forward and Inverse Dynamics of Elastic Manipulators," Piopeedings oi lPEE WW. 1986. PP- 359- 364. Davis,J.H. and Hirschorn,R.M., " Tracking Control of a Flexible Robot Link," lEEE Iianaaptipns on Auppmatic Conprol, Vol. 33 , No. 3 , 1988, pp. 238-248. Davis,R., Henshe11,R.D. and Warburton,G.B. , " A Timoshenko Beam Element," Jppinal pf Sound anp Vibration, Vol.22, No.4, 1972, pp. 475-487. Denavit, J. and Hartenberg, R.S. , " A Kinematic Notation for Lower- Pair Mechanisms Based on Matrices," ASME Jpninal pi Appliap Mpphanipa, No.22, June 1955. 50. 51. 52. 53. 54. 55. 56. 57. 58. 59. 60. 61. 62. 240 Elliott,L.E., Mingori,D.L. and Iwens,R.P., " Performance of Robust Output Feedback Controller for Flexible Spacecraft," ee 3 of o o tooLae xibeSacecaf, Blacksburg, Virginia, 1979, pp. 409-420. Everett,L.J., " Modeling and Robust Control of Compliant Robotic Manipulator via the Finite Element Method,” Ph,D, Ihesis, Texas A&M University, 1983. Ewins,D.J. , Mppal leaping ', Iheoiy and Practice, Research Studies Press LTD, Bruel 6: Kjaer, Letchworth, England, 1984. Fukuda,T. , " Flexibility Control of Elastic Robotic Arms," Journal oi Ppbopic Systems, Vol.2, No.1, 1985, pp. 73-88. Gamache,D. and Thompson,B.S., " The Finite Element Design of Linkages-A Comparison of Timoshenko and Euler-Bernoulli Elements," Pioceedinga pi ASME Qpnfeipnce on Design Engineeiing, Hartford, Connecticut, 1981, ASME Paper No. 81-DET-109. Gandhi, M.V. and Thompson, B.S. , ma er s and u es ev lut o , Technomic Publishing Company, Inc., Lancaster, PA, 1989 (in press). Gandhi, M.V., Thompson, B.S. and Choi, S.B., " A New Generation of Innovative Ultra-Advanced Intelligent Composite Materials Featuring Electro-Rheological Fluids: An Experimental Investigation, " Journal of Composite Mateiials, Vol.23, 1989, pp. 1232-1255. Gandhi, M.V., Thompson, B.S. and Choi, S.B., " A Proof-of-Concept Experimental Investigation of a Slider-Crank Mechanism Featuring a Smart Dynamically-Tunable Connecting-Rod Incorporating Embedded Electra-Rheological Fluid Domains," gournal oi Sound anp Vibgation, Vol. 135, No.3, 1989, pp. 511-515. Gandhi, M.V., Thompson, B.S., Choi, S.B. and Shakir, S., " Electro- Rheological-Fluid-Based Articulating Robotic Systems," ASME Journal pf Mephaniama, Iianamissions and Autpmation in Design, Vol.111, No.3, 1989, pp. 328-336. Gevarter,W.B. , " Basic Relations for Control of Flexible Vehicles," AW, Vol.8, 1970, pp. 666-672. Goldenberg,A.A. and Rakhsha,F., " Feedforward Control of a Single- Link Flexible Robot," Mechanism and Maphine Iheoiy, Vol.21, No.4, 1986, pp. 325-335. Gutman, S., ” Uncertain Dynamical Systems - A Lyapunov Min-Max Approach," IEEE Iiansappipna on Antpmapip Qpntipl, Vol. AC-24, No.3, 1979, pp. 437-443. Hale,A.L., Lisowski,R.J. and Dah1,W.E., " Optimal Simultaneous Structural and Control Design of Maneuvering Flexible Spacecraft," AIAA Jpninal pi Dnipance, gpntzol, and Dmamics, Vol.8, No.1, 1985, pp. 86-93. 63. 64. 65. 66. 67. 68. 69. 70. 71. 72. 73. 74. 75. 76. 241 Hanagud, S., Won, C.C. and Obal, M.V., " Optimal Placement of Piezoceramic Sensors and Actuators," Pgoceedings of American Control anfigzgnng, Aug. 1988, pp. 1884-1889. Hastingg, G.G. and Book, W.J., " Experiments in Optimal Control of a Flexible Arm," d th E e ence o eci 0 nd 99351.11, Athens, Greece, Dec. 1986, pp. 86-90. Hastings, G.G. , " Controlling Flexible Manipulators, An Experimental Investigation," W3, Department of Mechanical Engineering, Georgia Institute of Technology, 1986. Hughes,P.C., " Space Structure Vibration Modes: How Many Exit ? Which Ones are Important 7," IEEE Contrgi Systems Magazine, Feb. 1987, pp. 22-28. Iman, I. and Sander, G.N., " High-Speed Mechanism Design- A General Analytical Approach," ASME gnnznai oi Engineering for Industgy, Vol. 97, No. 2, 1975, pp. 609-238. Johnson,T.L., " Progress in Modeling and Control of Flexible Spacecraft," Jguzna]. of the Franklin institute, Vol.315, No.5/6, 1983, pp. 495-520. Jones, R.M., Mechanigs of Composite Materials, McGraw-Hill, Inc., New York, 1975. Joshi,S.M. , " Robustness Properties of Collocated Controllers for Flexible Spacecraft." MW M, Vol.9, No.1, 1986, pp. 85-91. Kapur,l(.l(. , " Vibrations of a Timoshenko Beam, Using Finite-Element Approach," Jgngnnl nfi the Acongticnl figgiegy 9i Amegiga, Vol.40, No.5, 1966, pp. 1058-1063. Khot,N.S. , Eastep,F.E. and Venkayya,V.B., " Optimal Structure Modifications to Enhance the Optimal Active Vibration Control of Large Flexible Structures," ed t t W. 1985. PP- 134-142- Korolov,V.V. and Chen,Y.H. , " Controller Design Robust to Frequency Variation in a One-Link Flexible Robot Arm," ASME Jgugnnl 9f Dynamic WM. Vol.119. 1989. pp- 9-14- Krishnan,H. and Vidyasagar,M., " Control of a Single-Link Flexible Beam using a Hankel-Norm-Based Reduced Order Model," 0 Inggznngignal anigrenge 9n Rohogics gnd Automation, 1988, pp. 9-14. Kuang,J.H. and Tsuei,Y.G. , " A More General Method of Substructure Mode Synthesis for Dynamic Analysis, " W, Vol.23, No.4, 1985, pp. 618-623. Laurenson,R.M., " Modal Analysis of Rotating Flexible Structure," 51W, Vol.14, No.10, 1976, pp. 1444-1450. 77. 78. 79. 80. 81. 82. 83. 84. 85. 86. 87. 88. 242 Lawson, A.W., " The Vibration of Piezoeletric plates," Dhygigal Ming, Vol.62, July 1942, pp. 71-76. Lee,H. and Castelazo,A. , " Nonlinear Feedback Control of a Flexible Robot Arms," W te n u eet , Boston, Mass., DSC-Vol.6, 1987, pp. 307-314. Leitman, G., " Guaranteed Asymptotic Stability for Some Linear Systems with Bounded Uncertainties," ASME Journai oi Dyngmic Sygtgng, Mgaguzemgnt, ang Control, Vol.101, 1979, pp. 212-216. Leitman, G., " On the Efficacy of Nonlinear Control in Uncertain Linear Systems," SME Jou o amic S stems Measuremen and Contrgi, Vol.102, 1981, pp. 95-102. Leitman,G, Ryan,E.P. and Steinberg, A., " Feedback Control of Uncertain Systems : Robustness with respect to Neglected Actuator and Sensor Dynamics," inteinaiional Journal of Control, Vol.43, No.4, 1986, pp. 1243-1256. - Liao,D.X. , Sung,C.K. and Thompson,B.S. , " The Design of Flexible Robotic Manipulators with Optimal Arm Geometries Fabricated from Composite Laminates with Optimal Material Properties," The Intemational ggumni oi Bobotigg Research, Vol . 6 , No . 3 , 1987 , pp. 116-130. Liegeois,A., Dombre,Z. and'Borre1,P., " Learning and Control for a Compliant Computer-Controlled Manipulator," IEEE Transactigng 9n WM, Vol.AC-25, No.6, 1980, pp. 1097-1102. Lin, D.X., Ni. R.G. and Adams, R.D., " Prediction and Measurement of the Vibration Damping Parameters of Carbon and Glass Fiber- Reinforced Plastics Plates," Qurnai of Comnosite hateiials, Vol. 18, March 1984, pp. 132-152. Lin,J.G. , Hegg,D.R., Lin,Y.H. and Keat, J.E. , " Output Feedback Control of Large Space Structures: An Investigation of Four Design Methods," e o 2 as on amics and Co 01 of Largs_Elsxihle_finesssra£s. Blacksburg. Virginia. 1979. pp. 1-18- Lin,Y.H. and Lin,J.G. , " An Improved Output Feedback Control of Flexible Large Space Structures," Diogeedingg 9i 12gb IEEE gnnigzgnge 9n Dggision nngi Continl, Albuquerque, NM, 1980, pp. 1248- 1250. Lindberg,R.E. and Longman,R.W. , " On the Number and Placement of Actuators for Independent Modal Space Control," AIAA Jguing]. 9i Q2iQen2si_ggn£rgli_end_nxnamiss. Vol.7. No 2. 1984. pp- 215-221- Low,K.H. , " A Systematic Formulation of Dynamic Equations for Robot Manipulators with Elastic Links," Jouinai oi ggbotic §y§tem§, Vol.4, No.3, 1987, pp. 435-456. 89. 90. 91. 92. 93. 94. 95. 96. 97. 98. 99. 100 . 101 . 243 Low,l(.H. and Vidyasagar,M., " Dynamic Study of Flexible Manipulators with Open and Closed Chain Mechanism," WWW mugging, Boston, Mass., DSC-Vol.6, 1987, pp. 277-286. Low,l(.H. and Vidyasagar,M. , " A Lagrangian Formulation of the Dynamic Model for Flexible Manipulator Systems," SM u a 0 MW. Vol.110. 1988. pp. 175- 181. Matsuno,F., Fukushima,S., Ohsawa,Y., Kiyohara,M. and Sakawa,Y. , " Feedback Control of a Flexible Manipulator with a Parallel Drive Mechanism," Ihg Ingennaiinnn]. Jouinai 9i Egbotics Resenich, Vol.6, No.4, 1987, pp. 76-84. McLaren,M.D. and Slater,G.L., " Robust Multivariable Control of Large Space Structures using Positivity," AIAA Journal oi Guidancg, Conizoi, and Dynamics, Vol.10, No.4, 1987, pp. 393-400. Meck1,P. and Seering,W. , " Active Damping in a Three-Axis Robotic Manipulator , " eed o S C re 0 Willi-£11318. Dearborn, MI, Sept. 1983, ASME Paper No. 83- DET-80. Meirovitch, L., Analytical Methgd in Vibraiions, MacMillan Publishing Co., New York, 1967. Meirovitch,L. , Baruh,H. , Montgomery,R.C. and Williams,J.P. , " Nonlinear Natural Control of an Experimental Beam," 51W Qnidance, Qonizni, nnd Dynamics, Vol.7, No.4, 1984, pp. 437-442. Meirovitch,L. and Silverberg,L.M., " Globally Optimal Control of Self-Adjoint Distributed Systems," Eiggeedingg 9f 4th Synnngim 9n 0 t o of ar e S c c aft, Blacksburg, Virginia, 1981, pp. Meirovitch,L. and Silverberg,L.M. , " Active Vibration Suppression of a Cantilever Wing," Journal, of Sound and Vibration, Vol.97, No.3, 1984, pp. 489-498. Meldrum,D.R. and Balas,M.J., " Application of Modal Reference Adaptive Control to a Flexible Remote Manipulator Arm," W 2W. Seattle. Washington. 1986. pp- 825- 832. Mesquita,L. and Kamat,M.P. , " Structural Optimization for Control of Stiffened Laminated Composite Structures," u S un d W, Vol.116, No.1, 1987, pp.33-48. Miller, S.B. and Hubbard, J.E., " Smart Components for Structural Vibration Control,” e e can Con ro Co fer nce, 1988, pp. 1897-1902. Mindlin, R.D. , " High Frequency Vibrations of Piezoelectric Plates," WWW. Vol.8. 1972. pp-895-906- 102 . 103 . 104 . 105 . 106 . 107 . 108 . 109 . 110. 111. 112. 113. 114. 244 Moerder, D.D. and Calise, A.J., " Convergence of a Numerical Algorithm for Calculating Optimal Output Feedback Gains," iEED Iznngngiinng on nuigmatig Cnntioi, Vol. AC-30, No.9, 1985, pp. 900- 903. Morimoto,Y. , Inamura,J. and Hizoguchi,l(. , " Dynamic Control of a Flexible Robot Arm by using Experimental Modal Analysis," Einggegings 9i ASME Wintei Annual Meeting, Boston, Mass., DSC-Vol.6, 1987, pp. 337-343. Mufti,I.H. , " Modal Reference Adaptive Control for Large Structural Systems,” AIAA Jouinai 9i Guidance, Coniroi, and Dynamics, Vol.10, No.5, 1987, pp. 507-509. Naganathan,G. and Soni,A.H., " Coupling Effects of Kinematics and Flexibility in Manipulators," e ernationa Journa 0 ob 5 31m, Vol.6, No.1, 1987, pp.75-84. Nelson,W.L. , Mitra,D. and Boie,R.A. , " End-Point Sensing and Load- Adaptive Control of a Flexible Robot Arm," oceedi 4th Qonieience on Decisign and Contiol, Ft.Lauderda1e, FL, December 1985, pp. 1410-1415. Nemir,D., Koivo,A.J. and Kashyap,R.L., ” Control of Gripper Position of a Compliant Link using Strain Gage Measurements," Erggeedings 9i Zfiih Cnnigigncg on Decigion and aniroi, Athens, Greece, December,l986, pp. 1140-1144. Neto,O.M., " Modal Analysis and Control of Flexible Manipulator Arms," W, Massachusetts Institute of Technology, 1974. Ni, R.G. and Adams, R.D. , " A Rational Method for Obtaining the Dynamic Mechanical Properties of Lamina for Predicting Stiffness and Damping of Laminated Plates and Beams," Composites, Vol. 15, No. 3, 1984, pp. 193-199. Nicosia,S., Tomei,P. and Tornambe,A., " Dynamic Modeling of Flexible Robot Manipulators," o eed n s f ter ationa Con e e ce Wm. 1986. pp. 365-372. Omatu,S., Koide,S. and Soeda,T., " Optimal Sensor Location Problem for a Linear Distributed Parameter Systems," EEE ransa o s on Aniomaiin antioi, Vol.AC-23, No.4, 1978, pp. 665-673. Onoda,J. and Haftka,R.T., " An Approach to Structure/Control Simultaneous Optimization for Large Flexible Spacecraft," ALLA mg, Vol.25, No.8, 1987, pp. 1133-1138. Patel, R.V. and Toda, M. , " Quantitative Measures of Robustness for Multivariable Systems," Eigggeding§ oi JACQ, San Francisco, CA, 1980, TP8-A. Pestel, E.C. and Leckie, F.A., Mat h d las ome ani s, McGraw-Hill Co., 1963. 115. 116. 117 . 118. 119. 120. 121. 122 . 123. 124. 125. 126. 127. 245 Plump, J.M., Hubbard, J.E. and Bailey, T., " Nonlinear Control of a Distributed System : Simulation and Experimental Results," ASME Journal of Dynamic Systems, Measuremeni, and Control, Vol. 109, June 1987, pp. 133-139. Preumont,A., " Spillover Alleviation for Nonlinear Active Control of Vibration," IAA ournal of Guidance Co and amics, Vol.11, No.2, 1988, pp. 124-130. Rakhsha,F. and Goldenberg,A.A. , " Dynamics Modeling of a Single-Link Flexible Robot," Droceedings oi inteinationai Conference on Roboiics and Automation, 1985, pp. 984-989. Ramirez,R.B. , " Design of a High Speed Graphite Composite Robot Arm," ASME Paper No. 84-DET-fiL Reddy, J.N. , An Introduction to the Finite Element Method, McGraw- Hill, Inc., New York, 1986. Rogers, C.A., Liang, C. and Barker, K.D., " Dynamic Control Concepts Using Shape Memory Alloy Reinforced Plates," U,S, Army Reseaich Office Worksho on Smart Materials St uctures and Mathematica Issues, Blacksburg, VA, 1988. Rogowski, R.S., Heyman, J.S. and Claus, R.O., " The Evolution of Smart Composite Materials," AS e h e 3, Vol. 12, No. 10, 1988, pp. 20-22. Rovner,R.M. and Cannon, R.H. , " Experiments Toward On-Line Identification and Control of a Very Flexible One-Link Manipulators,” e e atio a1 ourna o bo ic ese c , Vol.6, No.4, 1987, pp. 3-19. Sakawa,Y., Matsuno,F. and Fukushima,S. , " Modeling and Feedback Control of a Flexible Arm," Jnurnal of Robotic Systems, Vol.2, No.4, 1985, pp. 453-472. Schaechter,D.B. and Eldred,D.B., " Experimental Demonstration of the Control of Flexible Structures," AIM Journai of Guidance, Contiol, ang Dynamics, Vol.7, No.5, 1984, pp. 527-534. Schafer,B.E. and Holzach,H. , " Experimental Research on Flexible Beam Modal Control," I A J ur a e t o W, Vol.8, No.5, 1985, pp. 597-604. Schmitz,E. , " Dynamics and Control of a Planar Manipulator with Elastic Links,” Diogeedings of 2§th Conferenge on Deniginn ans; Qontroi, Athens, Greece, 1986, pp. 1135-1139. Schmitz, E., " Experiments on the End-Point Position Control of a Very Flexible One-Link Manipulator,“ Ph,D, Ihesis, Department of Aeronautics and Astronautics, Stanford University, 1985. 128. 129. 130. 131. 132. 133. 134. 135. 136. 137. 138. 139. 140 . 246 Shulman, Z.P., Gorodkin, R.G., Korobko, E.V. and Gleb, V.I(., " The Electro- Rheological Effect and its Possible Uses," Journal of Non- Newtoninn Fluid Mechanics, Vol. 8, 1981, pp. 29-41. Siciliano,B. and Book,W.J., " A Singular Perturbation Approach to Control of Lightweight Flexible Manipulators," Ihe International Journai oi Robotics Research, Vol.7, No.4, 1988, pp. 79-90. Siciliano,B. , Book,W.J. and Maria,G., " An Integral Manifold Approach to Control of a One Link Flexible Arm,” Dioceedings of Zfiih Qpnierence on Decision and Contipi, Athens, Greece, December 1986, pp. 1131-1134. Siciliano,B., Calise,A.J. and Jonnalagadda,V.R.P., " Optimal Output Fast Feedback in Two-Time Scale Control of Flexible Arms," oceed n of 25t Co erenc Decis and Control, Athens, Greece, 1986, pp. 1400-1403. Siciliano,B. , Yuan,B.S. and Book,W.J., " Model Reference Adaptive Control of a One Link Flexible Arm,” Dioceesiingg oi 25th Conierence on Decision nnd Control, Athens, Greece, 1986, pp. 91-95. Singh,S.N. and Coelho, A.A.R. , " Nonlinear Control of Mismatched Uncertain Linear Systems and Application to Control of Aircraft,“ ASME Jpnpnni pf Mamip Systems, Measuieneni, nng Controi, Vol.106, 1984, 00. 203-210. Singh,S.N. and Schy,A.A., " Control of Elastic Robotic Systems by Nonlinear Inversion and Modal Damping," MW §y§ien§, uen§nrement, nng Qontpol, Vol.108, 1986, pp. 180-189. Skaar,S.B., Michel,A.N. and Miller,R.K., " Stability of Viscoelastic Control Systems," EE Transactions u t on , Vol.33, No.4, 1988, pp. 348-357. Stanway, R., Sproston, J. and Firoozian, R., " Identification of the Damping Law of an Electra-Rheological Fluid : A Sequential Filtering Approach," M ou a o D nami ste a re t Coniioi, Vol. 111, 1989, pp. 91-96. Stevens, N.G., Sproston, J.L. and Stanway, R. , " An Experimental Study of Electra-Rheological Torque Transmission," ASME Jon,r_nai pi sm ss 0 n utoma ion in s , Vol. 110, 1988, pp. 182-188. Stevens, N.G., Sproston, J.L. and Stanway, R., " The Influence of Pulsed D.C. Input Signals on Electra-Rheological Fluids," quina], pi Electiostatics, Vol.17, 1985, pp. 181-191. Stevens, N.G. , Sproston, J.L. and Stanway, R., " On the Mechanical Properties of Electra-Rheological Fluids," ASME Jouinai pi Applied Mechanipg, Vol.54, 1987, pp. 456-458. Stangroom, J.E., " Electrorheological Fluids," Journai pf Dhygics Technology, Vol. 14, 1983, pp. 290-296. 141 . 142 . 143 . 144 . 145 . 146 . 147 . 148 . 149 . 150. 151. 152. 153. 247 Streit,D.A., Krousgrill,C.M. and Bajaj,A.K., " Dynamic Stability of Flexible Manipulators Performing Repetitive Tasks," Proceedings oi Bobotipg and Manufacturing Automation, 1985, pp. 121-136. Sunada,W.H. and Dubowsky,S. , " The Application of Finite Element Methods to the Dynamic Analysis of Flexible Spatial and Co-Planar Linkage Systems," SME Journal oi Mechanipai Design, Vol.103, 1981, pp. 643-651. Sunada,W.H. and Dubowsky,S. , " On the Dynamic Analysis and Behavior of Industrial Robotic Manipulators with Elastic Members," ASME urna of Mechanisms T ansmissions and Automation in Desi , Vol.1035, 1983, pp. 43-51. Thomas,D.L. , Wilson,J.M. and Wilson,R.R., " Timoshenko Beam Finite Elements , " ournal of Sound and Vibiaiipn, Vol . 31 , No . 3 , 1973 , pp . 315-330 . Thompson,B.S., " Composite Laminate Components for Robotic and Machine Systems: Research Issues in Design,” ASME Applied Mechanics Review, Vol.40, No.11, 1987, pp. 1545-1552. Thompson,B.S. and Sung,K.S. , " A Variational Formulation for the Dynamic Viscoelastic Finite Element Analysis of Robotic Manipulators Constructed from Composite Materials," ASME gpurnai of Mechanismg, Tiangniggion, and Antonntipn in Design, Vol.106, 1984, pp. 183-190. Thompson,B.S. and Sunk,C.K. , " The Design of Robots and Intelligent Manipulators using Modern composite Materials," ec d Machinp Ihgoiy, Vol.20, No.6, 1985, pp. 471-482. Thomson, W.T. , Iheog of Vibration with Applicaiiong, Prentice-Hall Inc., Englewood Cliffs, NJ, 1091. Tracy, J.J. and Pardoen, C.C., " Effect of Delamination on the Natural Frequency of Composite Laminates," MW Materialg, Vol.23, No.12, 1989, pp. 1200-1215. Truckenbrodt,A. , " Modeling and Control of Flexible Manipulator Structures," roceedin s 0 4th ISM- FToMM m os um Theor and Practice of Robots and Manipulators, Warsaw, 1981. Uejima, H., " Dielectric Mechanism and Rheological Properties of Electra-Fluids," gapnnese gouinal of Applied Dhysics, Vol. 11, No. 3, 1972, pp. 319-326. User's ua o LC-2 a d PLC 2 1 icro ro esso , Allen- Bradley, 1987. Usoro,P.B. , Nadira,R. and Mahil,S.S. , " A Finite/Lagrange Approach to Modeling Lightweight Flexible Manipulators," ASME Journal of Dynamic §y§pens, Measuiement, and Control, Vol.108, 1986, pp. 199- 205 . 154. 155. 156. 157. 158. 159. 160. 161. 162. 163. 248 Waite, S.R., Tatan, R.P. and Jackson, A., " Use of Optical Fibre for Damage and Strain Detection in Composite Materials," W, Vol. 19, No.6, 1988, pp. 435-442. Wang,P.K.C. and Wei,J.D., " Vibrations in a Moving Flexible Robot .Arm,” quinal pf Sound.and.Vibpption, Vol.116, No.1, 1987, pp.149- 160. Wang,S.H. , Hsia,T.C. and Wiederrich,J.L., " Open-Loop Control of a Flexible Robotic Manipulator," international qurnal of Robotics and Automation, Vol.1, No.2, 1986, pp. 54-57. Wang,W.J., Lu,S.S. and Hsu,C.F., 1" Output Feedback Control of a Flexible Robot Arm," Droceedings of 25th Conference on Deci§ion pnp Control, Athens, Greece, 1986, pp. 86-90. Wang,W.J., Lu,S.S. and Hsu,C.F., " Experiments on the Position Control of a One-Link Flexible Robot Arm," IEEE Iransactions on prppics and Automation, Vol.5, No.3, 1989, pp. 373-377. West-Vukovich,G.S., Davison,E.J. and Hughes,IK(3., " The Decentralized Control of Large Flexible Space Structures," lEEE Transactions on Automatic Control, Vol.AC-29, No.10, 1984, pp. 866- 879. Winslow, W.M., ” Induced Fibration Suspensions, " Journal of Applied Dhysicp, Vol. 20, 1949, pp. 1137-1140. Yaeger, J.R., " A Practical Shape-Memory Electromechanical Actuator,” Eipppedings pi Intpinational fiynpopinn on Autpmopive Ipchnology pnd Automation, Milan, Italy, 1984, pp. 633-642. Yuan,B.S. and Book,W.J., " A Robust Scheme for Direct Adaptive Control of Flexible Arms," roceedin s o S W1 te A a1 Meeting, Boston, Mass., DSC-Vol.6, 1987, pp. 261-268. Zalucky,A. and Hardt,D.E., ” Active Control of Robot Structure Deflection," u f na stems easu eme t d Conpiol, Vol.106, 1984, pp.63-69.