Understanding and Modeling The Behavior of A Harmonic Drive Gear Transmission
Understanding and Modeling The Behavior of A Harmonic Drive Gear Transmission
Timothy D. Tuttli
DTIC
ELECTE
SJN2
21993
•TR~l,,,N TAEMEN,
1 98
21~~~ 038 ~ 93-01051
11NIng -"I38•R
93
Form Approved
REPORT DOCUMENTATION PAGE I OMB No. oo4-0o188
Pu•ic reorting burden for this coffeclion of information is esttmated to are"ae I hour per response. including the time or reviewing instructions.' searching existing data sources,
gatt•ring and maintaining the data needed, and completing and reviewing the collection ot information. Send comments regarding this burden esunate or any other .,oect of this
collection of information. including suggestions for reducing this burden, to Washington Hleadkuarsers Se•vKieS. Directorate for information Ooerations and Reports. 121S Jefferson
Davis Highway. Suite 1204. Arlington. VA 22202-4302. and to the Office of Management and Bludgt. Paperwork Reduction Project (0704-0 1BB), Wal..hngton. DC 20503
1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED
May 1992 technical report
4. TITLE AND SUBTITLE S. FUNDING NUMBERS
Timothy D. Tuttle
(continued on back)
14. SUBJECT TERMS(key words) IS. NUMBER OF PAGES
harmonic drive harmonic drive gear dynamic modeling 1P
harmonic gear harmonic transmission 16. PRICE CODE
17. SECURITY CLASSIFICATION 1. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 20. LIMITATION OF ABSTRACT
OF REPORT OF THIS PAGE OF ABSTRACT
UNCLASSIFIED UNCLASSIFIED UNCLASSIFIED UNCLASSIFIED
NSN 7540-01-280-5500 Standard Form 298 (Rev 2-89)
P2"3c1,02d b ANSI $I Z39. 0
298-102)
Block 13 continued:
escalate during system resonance. Due to the interaction of these transmission properties,
dynamic response measurements showed surprisingly agitated behavior, especially around
system resonance. This behavior was rationalized in terms of the observed properties and
used to develop dynamic models to describe harmonic-drive operation. From these
models, I realized that non-linear frictional effects cannot be ignored in any accurate
representation. Additionally, if behavior around system resonance must be replicated, both
kinematic error and compliant behavior are essential. By introducing the effects of friction,
compliance, and kinematic error into a simplified representation of the gear-tooth meshing
geometry in the actual transmission, model performance, especially around system
resonance, improved substantially. Unfortunately, since accurate measurements of
harmonic-drive properties were elusive and highly dependent on operating conditions, it
seems unlikely that any harmonic-drive representation will be able to deliver reliable
performance unless detailed experimental observations of the specific harmonic-drive
system are available.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY
ARTIFICIAL INTELLIGENCE LABORATORY
Byu...........
....... .
Distribution I
Availability Codes
Avail and / or
Dist Special
&t I
0 1992 MASSACHUSETTS INSTrUT• oFTECHNOLOGY
ALL RIGHTs RESERVED
BY
TIMOTHY DOUGLAS TUTTLE
ABSTRACT
It was fun.
List of Figures .............................................................................. 7
3
Table of Contents
4
Table of Contents
5
6
Figure 1.1: The three components of a harmonic drive ................................... 15
Figure 1.2: An assembled harmonic-drive transmission ..................................... 15
Figure 1.3: The harmonic-drive operating mechanism ................................... 17
Figure 1.4: Typical transmission error displayed by harmonic drives ..................... 21
Figure 1.5: Typical stiffness profile displayed by harmonic drives ...................... 23
Figure 2.1: Harmonic-drive testing apparatus ................................................. 31
Figure 2.2: Cross-section of joint 1......................................................... 32
Figure 2.3: Joint 1 sensor locations ........................................................ 33
Figure 2.4: Cross-section of joint 2 ....................................................... 34
Figure 2.5: Joint 2 sensor locations ........................................................ 35
Figure 2.6: Cross-section of joint 3....................................................... 36
Figure 2.7: Joint 3 sensor locations ........................................................ 37
Figure 2.8: Experimental data acquisition system and controller ......................... 38
Figure 2.9: Joint 1 position-error waveform at the input rotation frequency .......... 45
Figure 2.10: FFT of joint I position error in cycles per input revolution ................. 45
Figure 2.11: Joint 2 position-error waveform at input rotation frequency ................. 46
Figure 2.12: FFT of joint 2 position error in cycles per input revolution ................. 46
Figure 2.13: Joint 3 position-error waveform at input rotation frequency ................. 47
Figure 2.14: FFT of joint 3 position error in cycles per input revolution ................. 47
Figure 2.15: High resolution FFT of joint 1 primary position-error ........................ 49
Figure 2.16: FFT of joint I position error in cycles per output revolution ................. 51
Figure 2.17: FFT of joint 3 position error in cycles per output revolution ................. 51
Figure 2.18: FFT of joint 2 high-frequency position error .............................. 53
Figure 2.19: FFT of joint 3 high-frequency position error .............................. 53
Figure 2.20: Apparatus for harmonic-drive stiffness tests ................................ 56
Figure 2.21: Joint 1 harmonic-drive stiffness curve ....................................... 58
Figure 2.22: Joint 2 harmonic-drive stiffness curve ....................................... 58
Figure 2.23: Joint 3 harmonic-drive stiffness curve ...................................... 59
Figure 2.24: Joint 1 stiffness curves for different wave-generator angles ................ 60
Figure 2.25: Joint I stiffness curves for different wave-generator preloads .............. 60
Figure 2.26: Joint I linear approximations to stiffness data .............................. 62
Figure 2.27: Joint 2 linear approximations to stiffness data .............................. 62
Figure 2.28: Joint 3 linear approximations to stiffness data .............................. 63
Figure 2.29: Cubic-fit of joint 1 harmonic-drive stiffn, data ........................... 64
Figure 2.30: Cubic-fit of joint 2 harmonic-drive stiffn. - data ........................... 65
Figure 2.31: Cubic-fit of joint 3 harmonic-drive stiffness data ........................... 65
Figure 2.32: Joint 1 aggregate-friction and motor-damping curves ...................... 74
Figure 2.33: Damping curve for the joint 1 harmonic drive .............................. 74
Figure 2.34: Joint 2 aggregate-friction and motor-damping curves ...................... 75
Figure 2.35: Damping curve for the joint 2 harmonic drive .............................. 75
Figure 2.36: Joint 3 aggregate-friction and motor-damping curves ...................... 76
Figure 2.37: Damping curve for the joint 3 harmonic drive .............................. 76
Figure 2.38: Joint I harmonic drive efficiency ........................................... 78
Figure 2.39: Joint 2 harmonic drive efficiency ............................................... 78
Figure 2.40: Joint 3 harmonic drive efficiency ............................................... 79
Figure 2.41: Cubic fit of joint I non-resonance harmonic-drive damping ........ 80
Figure 2.42: Cubic fit of joint 2 non-resonance harmonic-drive damping ........ 81
7
List of Fieures
Figure 3.1: Schematic of apparatus design for joints 1 and 2 ................................ 101
Figure 3.2: Schematic of apparatus design for joint 3 ........................................ 101
Figure 33: Lumped-parameter model for joints 1 and 2 ..................................... 102
Figure 3.4" Lumped-parameter model for joint 3 .............................................. 103
Figure 3.5: Variable notation for the model of joints 1 and 2 ................................ 104
Figure 3.6: Variable notation for the model of joint 3......................105
Figure 3.7: An ideal harmonic-drive representation ........................................... 109
Figure 3.8: Notation for the ideal harmonic-drive model ..................................... 110
Figure 3.9: Joint I ideal-model input-velocity for a 1.6 amp step ........................... 110
Figure 3.10: Joint 2 ideal-model input-velocity for a 1.4 amp step ......................... 111
Figure 3.11: Joint 3 ideal-model input-velocity for a 0.28 amp step ........................ 111
Figure 3.12: A harmonic-drive model with friction ........................................... 113
Figure 3.13: Notation for the harmonic-drive model with friction .......................... 115
Figure 3.14: Joint 1 input-velocity response for a model with friction ..................... 118
Figure 3.15: Joint 2 input-velocity response for a model with friction ..................... 119
Figure 3.16: Joint 3 input-velocity response for a model with friction ..................... 120
Figure 3.17: A harmonic-drive model with friction and compliance ........................ 123
Figure 3.18: Notation for the model with friction and compliance .......................... 124
Figure 3.19: Harmonic-drive model 4 .......................................................... 128
Figure 3.20: Notation for harmonic-drive model 4 ........................................... 129
Figure 3.21: Joint 1 stiffness profile for harmonic-drive model 4 .......................... 134
Figure 3.22: Joint 2 stiffness profile for harmonic-drive model 4 .......................... 134
Figure 3.23: Joint 3 stiffness profile for harmonic-drive model 4 ....................... 135
Figure 3.24: Close-up of harmonic-drive gear-tooth action .................................. 141
Figure 3.25: Schematic of harmonic-drive model with gear-tooth geometry .............. 142
Figure 3.26: Harmonic-drive gear-tooth contact zones ....................................... 145
Figure 3.27: Revised transmission model for negative-torque loading ..................... 146
Figure 3.28: Kinematic notation and sign conventions for model 5 ........................ 148
Figure 3.29: Torque notation and sign conventions for model 5 ............................ 149
Figure 3.30: Joint 1 stiffness profile for harmonic-drive model 5 .......................... 155
Figure 3.31: Joint 2 stiffness profile for harmonic-drive model 5 .......................... 156
Figure 3.32: Joint 3 stiffness profile for harmonic-drive model 5 .......................... 156
8
List of Figures
9
List of Fieures
Figure F.1.8: Joint 1 simulated dynamic position error at 4.0 amps............ 298
Figure F.2. 1: Joint 2 input-position simulated step-responses ............................. 299
Figure F.2.2: Joint 2 input-velocity simulated step-responses ...................... 300
Figure F.2.3: Joint 2 input-current simulated step-responses .............. 301
Figure F.2.4: Joint 2 output-position simulated step-responses ............................ 302
Figure F.2.5: Joint 2 output-velocity simulated step-responses ............................ 303
Figure F.2.6: Joint 2 output-velocity simulated step-responses ............................ 304
Figure F.2.7: Joint 2 output-torque simulated step-responses .............................. 305
Figure F.2.8: Joint 2 simulated dynamic position-error at 2.6 amps ....................... 306
Figure F.3. 1: Joint 3 input-position simulated step-responses ............................. 307
Figure F.3.2: Joint 3 input-velocity simulated step-responses ............. 308
Figure F.3.3: Joint 3 input-current simulated step-responses .............. 309
Figure F.3.4: Joint 3 output-position simulated step-responses ............................ 310
Figure F.3.5: Joint 3 output-velocity simulated step-responses ............................ 311
Figure F.3.6: Joint 3 output-velocity simulated step-responses ............................ 312
Figure F.3.7: Joint 3 output-torque simulated step-responses .............................. 313
Figure F.3.8: Joint 3 simulated dynamic position error at 0.44 amps ...................... 314
10
T Le. ef Tab eso
Table C. 1: Resistance values for the motor amplifier model ................................. 234
11
12
A good mechanical designer holds a clear awareness of the capabilities and
limitations of his design. For most designers, this awareness rests heavily on a firm
understanding of the individual components that comprise the entire assembly. For
mechanical systems employing harmonic drives, a designer can gain confidence in the
performance of the overall mechanism if the operation of the harmonic drive is clearly
understood. Additionally, for applications, such as robot control, which require
mathematical models to describe system performance, a thorough understanding of the
behavior of harmonic-drive transmissions is essential. The purpose of my investigation is
to identify important attributes that influence the operation of harmonic-drive transmissions
and to relate these attributes to overall dynamic behavior. On the basis of insight gained
from this empirical analysis, I hope to develop mathematical models that can be used to
predict reliably the dynamic performance of harmonic-drive transmissions.
13
Chapter 1: Introduction
Results from these models were encouraging and informative, but more detailed
representations are required to capture fully the subtleties of harmonic-drive behavior.
Every harmonic drive is composed of the three distinct parts illustrated in figure
1.1: the wave-generator, the flexspline, and the circular spline. The wave-generator is a
ball-bearing assembly with a rigid, elliptical inner-race and a flexible, thin-walled outer
race. The flexspline is a thin-walled, flexible cup adorned with small, external gear teeth
around its rim. When assembled, the wave-generator is nested inside the flexspline,
causing the flexible gear-tooth circumference on the flexspline to adopt the elliptical profile
14
Chapter)1: Introduction
15
Chapter)I: Introduction
of the wave-generator. The circular spline is a rigid ring with internal teeth machined along
a slightly larger pitch diameter than those of the flexspline. When the flexspline-and-wave-
generator assembly is inserted into the circular spline, the external teeth on the flexspline
mesh with the internal teeth on the circular spline along the major axis of the wave-
generator ellipse. A correctly assembled harmonic drive with a cut-away section is shown
in figure 1.2.
Since the harmonic drive has three rotational ports, it boasts a versatility unavailable
to standard gear transmissions. Specifically, in its most popular configuration, the circular
spline is fixed to ground and a low-torque, high-speed motor driving the wave-generator
can produce high-torque, low-speed rotation on the flexspline. Similarly, with the
flexspline mounted to ground, the torque on the wave-generator can be magnified and
transmitted through the circular spline. If either the circular spline or flexspline is used as
an input, the wave-generator can be driven at high velocities and low torques. In general,
by using different combinations of rotations on the three harmonic-drive components,
numerous differential- gearing functions and reduction ratios can be achieved.
16
Chapter 1: Introduction
. . . .
.r . . ... ...
.fh . ........ ~ a k L r II
Intaly th0 lx ln n
efficiencies can approach 90 percent. The smallest transmissions can provide a respectable
maximum torque output of about 3.5 N-m (30.0 in-lbs), while the heavy-duty units boast
up to 10,000 N-m (89,000 in-lbs) of torque capacity.
Because of its unique operating principles, the harmonic drive displays performance
features both superior and inferior to conventional gear transmissions. These advantages
and disadvantages are itemized in Table 1.1 and Table 1.2.
The unique performance features of the harmonic drive have captured the attention
of designers in many fields. Specifically, due to their accuracy and simple construction,
High torque capacity Since torque is transmitted through muftiple-tooth contact, harmonic
drives can withstand high loads at small pitch-diameters.
Concentric geometry Since all three harmonic-drive components are concentric and coaxial,
designers can drastically reduce power-train size and complexity.
Lightweight and Requiring only three basic elements, the harmonic drive can deliver
compact design extremely high gear ratios in a small package.
Zero backlash Natural gear preloading and predominantly radial tooth-engagement
eliminate virtually all transmission backlash.
High Efficiency Ifproperly lubricated, typical efficiencies of harmonic-drive
transmissions can reach 80 to 90 percent.
Backdrivability Due to their high efficiency, wave-generator rotation can be driven
through the flexspline or circular spline.
High flexibility Due to the high loads seen by the wave-generator and gear teeth,
moderate operating torques can produce substantial transmission
torsion.
Kinematic error Due to manufacturing inaccuracies, harmonic drives exhibit a small but
ubiquitous position error across the transmission.
Resonance vibration Since torque fluctuations produced by kinematic error can interact with
the low stiffness of the transmission to excite resonances, high
vibration amplitudes may be generated in some operating ranges.
Non linearity Both the flexibility and frictional losses in the drive exhibit highly
nonlinear behavior.
Poorly understood Compared to conventional gear transmissions, relatively little
exploration has been done of the unusual operating mechanisms in
the harmonic drive.
18
Chapter 1: .- rroduction
harmonic drives have been used successfully in industrial robots, assembly equipment,
and measuring instruments. Heavy-duty applications such as machine tools and printing
presses have also utilized harmonic drives for their high-torque capabilities. Additionally,
space and aircraft systems often employ harmonic drives for their lightweight and compact
geometry.
Throughout its short existence, the harmonic drive has enjoyed continuously
increasing international attention from designers as well as researchers. In the Soviet
Union, substantial research, predominantly performed on a harmonic drive featuring a
wave-generator with two rollers instead of an elliptical bearing, has provided valuable
insight into the structural and dynamic properties of the drive as well as its suitability for
heavy-duty applications. More recently, Japanese researchers have explored such areas as
transmission stiffness, positioning accuracy, and tooth-meshing mechanisms in an effort to
develop new gear-tooth geometries and improve harmonic-drive performance.
Investigating harmonic drive performance in robotics and aerospace applications,
researchers in the United States have also contributed constructive insights about the limits
of harmonic-drive applicability.
The reference list at the end of this document contains citations for numerous papers
encompassing the scope harmonic-drive research worldwide. Primarily, this research has
been focused in six major areas: (1) vibration and dynamic behavior, (2) kinematic error,
(3) torsional stiffness, (4) structural analysis, (5) robotics applications, and (6) gear-tooth
geometry. Since my investigation is particularly focused on the dynamic and static
performance of harmonic drives, past research on dynamic behavior, positional error, and
torsional stiffness received most of my attention. In an effort to overview the current state
of harmonic drive technology, a description of the fundamental mechanisms of harmonic-
drive operation as well as relevant insights presented in existing literature will be
summarized in this section.
19
Chapter): Introduction
rotation on two of the three harmonic-drive ports as well as a value for N, the ideal rotation
of the third harmonic-drive port can be predicted by the equation
where Ows is the rotation of the wave-generator, 0cs is the rotation of the circular spline,
and Of, is the rotation of the flexspline. All three rotations in this equation are defined in
the same frame of reference. Similarly, given that N is constant under ideal assumptions,
the derivative of this relationship yields a similar velocity constraint:
,swg
= (N + 1) wcs - N (,fs, (1.2)
where cowg. (0cs, and cofs represent the angular velocities of the three harmonic-drive
components relative to the same velocity reference. From these equations, it can be seen
that, if the circular-spline velocity is zero, the wave-generator will rotate N-times faster than
the flexspline in the opposite direction. In a different configuration, a grounded flexspline
dictates that the wave-generator will spin (N+I) times faster than the circular spline in the
same direction. By applying the law of power conservation to the three ports of the
harmonic drive, the ideal torque behavior of the transmission can also be defined:
T (N+I1)N
Tw= I T.. =- - Tfs, (1.3)
where Twg, TCs, and Tfs are the torques, defined with identical sign conventions, seen by
the three harmonic-drive components. From this identity, it can be seen that the torque on
one port of the harmonic drive dictates the ideal torque seen by the other two ports. Notice
that the torque on the flexspline is nearly equivalent to the torque on the circular spline,
which is approximately N-times greater than the torque on the wave-generator. By
applying the torque, velocity, and position equations presented above, the dynamic
behavior of a harmonic drive transmission under ideal assumptions can be completely
defined.
20
Chapter 1: Introduction
the harmonic drive from the input rotation scaled by the ideal gear ratio for the given
transmission configuration:
0 err = 0i -
0 out
gear ratio (1.4)
21
Chapter 1: Introduction
Additional observations in [281 and [24] record the presence of the significant
amplitude modulation or beating in position error waveforms of most harmonic drives
22
Chapter 1: Introduction
while a few transmissions showed little or no amplitude modulation. Both authors concur
that the presence of amplitude modulation over several wave-generator cycles is caused by
the summation of the circular spline and flexspline gear-error components; since these two
error waveforms appear at slightly different frequencies, constructive and deconstructive
interference of the two error signals is manifested by low-frequency beating.
I
two characteristic properties of harmonic-drive flexibility: increasing stiffness with
23
Chapter 1: Introduction
24
Chapter 1: Introduction
a single drive can vary as much as ±30 percent. Nye, [65], confirms the magnitude of this
variability and observes its dependence on the integrity of internal fit-up. Specifically, he
notes that stiffness typically increases with tighter meshing of gear teeth and increased
internal loading. Experimental results from Nye and Hikada, et al., [ 11], indicate that,
since wave-generator orientation often influences gear-tooth loading, transmission stiffness
can vary significantly with wave-generator rotation. Further investigation by Nye
produced theoretical loading models which successfully described the influence of gear-
tooth preloading on stiffness and predicted preload- and stiffness-relaxation due to gear-
tooth wear.
25
Chapter]: Introduction
opted to capture harmonic-drive dynamics using linear stiffness and damping relationships
while others, [44, 51, 52], required non-linear representations to ensure model accuracy.
In order to capture the torque fluctuations induced by kinematic error, a few researchers,
such as [11, 18], supplement their models with a sinusoidal torque function derived from
position-error fluctuations or measured directly. For many applications in which non-ideal
transmission behavior becomes dominant, harmonic-drive manufacturer's frequently
recommend adjustment of the operating envelope to a region less plagued by transmission
dynamics. Some research, such as [13], has been focused on developing simple models to
predict operating regions in which harmonic drives are well-behaved.
Depending on the design and operating conditions, several researchers note that
strikingly non-ideal behavior can be exhibited in harmonic-drive systems. Specifically,
experimental measurements by Hikada, et al., [111, and Volkov, et al, [18], reveal that
harmonic drives typically exhibit vibrations at twice the frequency of wave-generator
rotation and several subsequent harmonics. It was also noticed that the amplitude of
vibration increased when one of these frequencies coincided with the natural frequency of
the harmonic-drive system. Furthermore, plots of resonance vibration presented by
Hikada, et al., illustrated the same twice-per-output-revolution amplitude modulation seen
in position-error signatures. By developing a simple non-linear vibration model which
captured the behavior of tooth-meshing errors, Hikada, et al., illustrated the direct influence
of position error on dynamic vibration. Further comparison between model predictions and
experimental results revealed the importance of representing the variation in stiffness with
wave-generator angle to ensure model accuracy. Another non-linear vibration model
developed by Volkov, et al., was able to predict observed system non-linearities such as
vibrations at subharmonics of the wave-generator excitation frequency due to the non-linear
harmonic-drive stiffness.
26
Chapter 1: Introduction
and dynamic friction demonstrated that energy dissipation in harmonic drives varies non-
linearly with velocity and periodically with rotatic and can escalate during system
resonance. By comparing these empirical observations to manufacturer's estimates, I
discovered that, due to the fickle nature of most harmonic-drive attributes, catalog
predictions often fail to describe the observed behavior. Following this investigation of
individual properties, overall dynamic response measurement showed surprisingly agitated
transmission behavior. Specifically, due to the problematic influence of system resonance
in typical harmonic-drive operating ranges, dynamic response data demonstrated violent
torque and velocity fluctuations as well as sudden unexpected increases in operating
velocity. Using the understanding of transmission properties gained from experimental
analysis, this unusual dynamic behavior was rationalized.
27
Chapter 1: Introduction
kinematic-error, (2) stiffness, (3) starting torque, and (4) dynamic friction. For each of
these topics, the testing procedure will be described, experimental results will be presented,
and conclusions will be made about the observed behavior. Using this empirical insight,
dynamic-response measurements will then be characterized and rationalized in terms of
important transmission properties. With experimental dynamic measurements to serve as a
metric for model performance, my theoretical exploration of harmonic-drive operation will
then be outlined. This presentation begins with a discussion of the dynamic models
developed to simulate the testing apparatus followed by a description of five different
harmonic-drive representations. For each of these transmission models, equations of
motion will be derived, accuracy and reliability will be evaluated, and conclusions will be
drawn about the potential usefulness and applicability of each representation. Based on
results from these theoretical models as well as the findings of the experimental inquiry,
conclusions and recommendations will be made about rewarding areas for future research.
28
Cha te 2: E p rm na l Anay i
29
Chapter 2: Description of Testing Apparatus
S..
. . .... ... .
1M 5C 1C
160:1 160:1 80:1
.... . '.............
...................
30
Chapter 2: Description of Testing Apparatus
31
Chapter 2: Description of Testin q Apparatus
32
Chapter 2: Description of Testing Apparatus
The design of the first harmonic-drive testing station, joint 1, is shown in figure
2.2. From this cross-section, it can be seen that the flexspline of the harmonic drive is
mounted to ground while the wave-generator is driven by a DC motor mounted to the
circular spline. The output link of this joint is supported by a pair of angular-contact
bearings and is driven by the circular spline of the harmonic drive. This joint also contains
two encoders, one torque sensor, and a tachometer whose locations are illustrated in figure
2.3. From this diagram, it can be seen that the input encoder and tachometer measure the
rotation of the motor shaft with respect to the motor housing which is mounted to the
circular spline. The output encoder is mounted to ground and measures the absolute
rotation of the circular spline through a gear pair with an 8:1 reduction ratio. I-astly, since
the only rigid connection between the drive train and ground is the torque sensor at the
flexspline output, this sensor can be used to measure directly the torque across the
transmission as seen by the flexspline.
33
Chapter 2: Description of TestinR. Apparatus
34
Chapter 2: Description of Testing Apparatus
0 0 mil 0
35
Chapter 2: Description of Testing Apparatus
The third and final testing station, as illustrated in figure 2.6, relies on a design
configuration substantially different from the two other joints. In particular, while the
wave-generator is still driven by a motor mounted to the circular spline, the output rotation
of the joint is transmitted through the flexspline and the circular spline is attached to
ground. This unique design not only allows for unrestrained rotation of the output link but
also provides a facility for examining harmonic-drive behavior under a different operating
configuration. From the sensor locations outlined in figure 2.7, it can be seen that, unlike
the other two joints, the input encoder directly measures absolute motor rotation and there
is no input tachometer. On the output side, a resolver measures the flexspline rotation with
respect to ground through a gear-pair with a 7.2:1 reduction ratio. Since the torque sensor
connects both the motor and circular spline to ground, it can measure directly the combined
torque of the wave-generator and circular spline, or equivalently, the torque seen at the
flexspline.
36
Chapter 2: Description of Testing Apparatus
.IIT
In order to monitor sensors and control motors on each joint, the data-acquisition
system and motor-controller illustrated in figure 2.8 was used. Using a software
environment called CONDOR, developed at the MIT Artificial Intelligence Laboratory,
computer code was compiled on a Sun 3/180 workstation and loaded onto a series of
Ironics processor boards. When executed, this code could then collect sensor information
from the data-acquisition boards and sensor-processing circuits shown in figure 2.8. To
control the motors, desired motor currents were specified on the D to A board and then
amplified through an Aerotech 4020-LS DC Servo Amplifier. By monitoring the amplifier
current with the A to D board, measurements of motor current and torque could be
gathered.
37
Chapter2: Description of Testing Apparatus
2.3, the input and output rotations with respect to ground are defined in table 2.4, and the
sensor-measurement relationships are defined in table 2.5. Second, as presented in table
2.6, sensor resolution and accuracy -must be characterized, and, lastly, sensors must be
calibrated. Experimental calibration procedures and results are reported in Appendix A.
From this catalog of sensor information, experimental data can now be used to accurately
measure the dynamic, static, and kinematic performance of the three harmonic-drive test
specimens.
On. Ofs, and 0=s Angular position of the wave-generator, flexspline, and circular spline,
respectively, as referenced from ground.
64g, o"1s, and (w Angular velocity of the wave-generator, flexspline, and circular spline,
respectively, as referenced from ground.
Twg, Tt,, and T. Torque experienced by the wave-generator, Ilexspline, and circular
spline, respectively.
.....
.... . A
38
Chapter 2: Description of Testin e A Pparatus
Tabe .5:Sesormesurmetsintem
of haronic-driv vaoiables
Output
~Vetpct............ .... o.eno
Output; noou
senornsot O
.. I.on
~h .... 2 t
.. ...
......
Chapter 2: Kinematic Error Measurement
In order to collect position-error data, the three harmonic-drive testing stations were
used to measure input and output rotation over a range of motion. By assuming quasi-
static conditions, kinematic-error data was collected while the joint was rotated at a small
but non-zero velocity. The typical procedure for data-acquisition and processing is listed
below.
n By sending a constant current to the joint motor that exceeds the transmission
starting torque, rotate the harmonic drive at the slowest velocity at which
resonance vibration is minimal.
F21 Collect input and output position information from the joint sensors at equally
spaced time intervals over a given motion range.
40
Chapter 2: Kinematic ErrorMeasurement
err = geagerratio
tio - out. (2.1)
Plot the resulting position-error signal versus the number of input, or wave-
generator, revolutions to analyze the shape of the error signature.
F5] Linearly interpolate the time-based position-error signal to produce data points
that are equally spaced by input rotation.
[•] Take an FFT of the interpolated data vector to identify the frequency
components of the error signature in terms of input revolutions.
By adjusting the range of transmission rotation and the number of data points
accumulated, FFT spectra spanning a wide range of frequencies were collected. Assuming
that at least two data points per cycle are required to measure position-error frequency
components, position-sensor resolution limits the cycles per angular position, or spatial
frequency, that can be measured on each harmonic drive. These spatial frequency upper-
bounds are listed in table 2.7.
outpt enoder250.0
output resolver 1365.0
41
Chapter 2: Kinematic Error Measurement
occurring once every output revolution, and harmonics, were contaminated by inaccuracies
in the output sensor gear which appeared at the same frequency. To distill the contribution
of the harmonic drive from this aggregate error waveform, a second set of position-error
data in which the harmonic drive was physically rotated by a fraction of a revolution
relative to the output sensor gear was required. In other words, after the first set of
kinematic-error data was obtained, the angular position of the harmonic-drive relative to the
output gear was shifted by a known amount before the second data set was collected. By
taking an FFT of these two signals, subtracting them, and scaling the result appropriately,
the harmonic-drive error-component could be accurately identified.
The theoretical basis for this experimental procedure can be clarified by assuming
that the position-error waveform collected from the harmonic drive testing apparatus has
three components: (1) hi(0), a first error signal at a spatial frequency of once-per-output
revolution , (2) h2 (0), a second signal at a spatial frequency of once-per-output revolution,
and (3) g(0), a third signal composed of various other components at different frequencies,
where 0 is a measurement of angular position. Therefore, the aggregate position-error
signal, h,. 1 (0), can be defined as
Now assume that a second signal is collected in which the first component was rotated by a
known phase shift, AO, in the positive 0-direction. This signal, hsu 2 (0), is described by
where co is the spatial frequency in radians per rotational unit. By subtracting these two
equations, it can be seem that the H2 (4o)and G(n) terms disappear leaving an equation in
terms of H1 (w) and the two aggregate error signals:
42
Chapter 2: Kinematic Error Measurement
The procedure used to collect and process two sets of position-error data according
to this theory is itemized below.
W Set the zero setpoint for input and output position sensors on the given testing
station.
n2 Move the joint at a slow, non-resonating velocity by sending a current to the
joint motor.
[•J Read input and output position data at equally spaced time intervals.
F51 Filter this position error data at every time step using a linear FIR digital low-
pass filter to prevent aliasing of high-frequency position error components
into sampled data.
n6 Return the joint to its original position after this data run is complete.
[F] Interpolate the first set of position error data to generate data points that are
equally spaced by input rotation rather than time.
r8] Introduce a phase shift into either the harmonic drive, the output gear, or both
by rotating that component a measured number of degrees.
n9 Collect the second set of data using the same procedure used to gather the first
set.
1-0 Interpolate this data set over the identical range of rotation used in the first set
of data.
1 Subtract the two resulting vectors of frequency domain data at all frequencies
and scale the difference according to equation (2.6) to yield the magnitude
spectrum of the phase-shifted position-error component.
F'3ý Take an inverse FFT of this spectrum to produce the waveform of the phase-
shifted position-error component if desired.
43
Chapter 2: Kinematic Error Measurement
Figures 2.9 through 2.14 illustrate the typical position-error waveforms observed
on the three harmonic drives as well as the spatial FFT results which indicate the frequency
content of these error signatures. For all of these plots, position-error magnitude is
measured in degrees of output (flexspline or circular spline) rotation and plotted in terms of
input (wave-generator) revolutions. In both the FFT and waveform plots, significant
energy content can be seen at frequencies lower than one cycle-per-input-revolution. These
low-frequency error components are due primarily to position inaccuracies in the output
sensor gear train and should not be mistaken for harmonic-drive error. Unfortunately, the
tooth-passing frequency of the output sensor gear-train appears on the joint 2 FFT
spectrum at 3.0 cycles-per-input-revolution. When viewing this plot, the error component
observed at this frequency should not be attributed to the harmonic-drive transmission.
Table 2.8 summarizes the peak-to-peak error magnitudes observed on the three harmonic
drives.
44
Chapter 2: Kinematic Error Measurement
0.040
0.030-
0.020 ...
- 0.010..........
0.000 .4 ..
.2-
-0.020 ............. .... ....
....
-0.030 - --
0 0.5 1 1.5 2
Input Revolutions
S0.010
FL. .. ... ... ... ... ... ... ... .... ... ... ... ... ... ... ... ....
... . . . . . . .
S 0 .0 0 6 .......................
..... ... .............. ........
.2 0.004
.g0.0020
45
Chapter 2: Kinematic Error Measurement
0.040-
0 .0 3 0 i...............
.... ..... ....
....
....
(A
0.020
............... ......... .... ......... .... ...
W 0.1 ......- . ..... .........
w -0.010 . ..
.0.020
-0.040 F - "
..... .... . . . . . . .
0 0.5 1 1.5 2
Input Revolutions
a)
P)
S0.0101
0
0 .00 8 -.............................
....... ......
Cu -
r 0 .004 .........-.................
-2 0.002
0.00
0 2 4 6 8 10 12 14 16 18 20
Frequency (cyclestinput revolution)
46
Chapter 2: Kinematic Error Measurement
0.040
0.030-
- 0.010 ....-.
.. ...... ........... .
wU -0.010.............
U)-0.020 .. ....
0
-0.030 ..........-..........
..... .... 7.......
....
-0.040-
0 0.5 1 1.5 2
Input Revolutions
S0.012
4)
'0 0.010
0.008
0...0
0.6...
. ........
................ ... ......
...
0. 0 .. ...... ............ ... ... ... ...
. .. . . ... ... ....... ....... ............;... ... ... ....
0.002...............
4 ..... ...-..
0.000
0 2 4 6 8 10 12 14 16 18 20
Frequency (cycles/input revolution)
47
Chapter 2: Kinematic Error Measurement
on
Ion ont 3
St es .u
... ••.....ut
....... .. 0.045 0.030 0.060
(utut*)0.042
QU~G~vlu 0.053 0.095
48
Chapter 2: Kinematic Error Measurement
particular, for a few tests, the largest position-error component appeared at four cycles-per-
wave-generator revolution rather than two, and harmonic content of the aggregate error
waveform increased substantially. Fortunately, no significant changes in the overall
position-error magnitude were observed in the testing apparatus due to joint reassembly.
However, as Nye and Kraml note in [24], since alignment can influence the magnitude of
kinematic error in harmonic drives, I would not have been surprised if less careful
reassembly resulted in different kinematic-error measurements.
) )-I
0.004
0.0031m
..
..
..
..
.. . ......
........ ... ... ..
..........
................. ...... ... ... i...... .........
............ .......... .............! ............. . .
E
c~0.002........
0.00
0..001............ . .............
:.......;. . . . ......... 7
.2............................... ..............
0.001
0.000
1.96 1.98 2 2.02 2.04
Frequency (cycles/input revolution)
49
Chapter 2: Kinematic Error Measurement
placement error and realizing that, during joint rotation, the flexspline moves relative to the
circular spline. Specifically, since the input position sensor on joint I measures wave-
generator rotation with respect to the circular spline, one-half revolution of the wave-
generator will rotate the tooth-meshing zone exactly 180 degrees around the circular spline
and complete one-cycle of the circular spline gear-error waveform. At the same time, since
the circular spline will rotate one tooth ahead of the flexspline after this 180 degree wave-
generator rotation, the flexspline gear-tooth error waveform will have corrn!eted one
complete cycle plus one gear tooth or (1 + (1/N)) cycles. Therefore, since the gear ratio,
N, on joint 1 is 160, for every complete rotation of the wave-generator the circular spline
gear-tooth error will have completed two complete cycles and the flexspline gear error will
have completed 2 (1 + (1/160)) or 2.0125 cycles. From this explanation, it can be seen that
the error component in figure 2.15 that appears at two cycles-per-input-revolution is due to
the circular spline gear-tooth-placement error and the error component that appears at
2.0125 cycles-per-input-revolution is due to the flexspline gear-tooth-placement error.
Figure 2.15 also illustrates an error occurring at 2.00625 cycles-per-input-revolution which
is probably due to the non-sinusoidal nature of the other two neighboring components.
50
Chapter 2: Kinematic Error Measurement
c 0.012 -
oa.
0.000
0 1 2 3 4 5 6 7 8 9 10
Frequency (cycles/output revolution)
0...L.........
....
... ... ... .......... ... ... ...
E. ... ..;..
........... .........
............ .......
.... ...
eu006-
0.0
0 1 2 3 4 5 6 7 8 9 10
Frequency (cycles/output revolution)
51
Chapter 2: Kinematic Error Measurement
to the frequency of output (circular-spline or flexspline), rather than input, rotation. From
figure 2.16, it can be seen that the joint 1 harmonic drive exhibits a low-frequency error
component occurring once per output revolution. The FFT peak that appears at 8.0 cycles-
per-output-revolution is caused by inaccuracy in the output sensor pinion and should be
ignored. In Figure 2.17, the error components for the joint 3 harmonic drive appear
primarily at 1.0 and 3.0 cycles-per-output-revolution. Again, the FFT spike occurring at
7.2 cycles-per-output-revolution is due to the output sensor pinion and should not be
mistaken for harmonic-drive error. In both of these plots, the aggregate. low-frequency
error component due to harmonic-drive transmission-error remains significantly less than
the harmonic-drive error observed at multiples of the input rotation frequency.
Unfortunately, due to the limitations of the testing apparatus, the low-frequency position
error of the joint 2 harmonic drive could not be measured.
52
Chapter 2: Kinematic Error Measurement
0.. ..........
...... . ....
.............. .......
S0.004 ..
............. ............ ...... .. . . .
........................... . . ......
=0 0.002 -... . . . . . . . . . . .
0.000 mm - . E-
0 100 200 300 400 500 600
Frequency (cycles/input revolution)
0
...................................
........ .....
CL
m 0.006 ..........
0.000 *- - iunm .
0 50 100 150 200 250 300 350 400
Frequency (cycles/input revolution)
53
Chapter 2: Kinematic Error Measurement
revolution. Similarly, since the joint 3 harmonic drive contains 160 gear teeth, the tooth-
passing frequency should be seen at 160 cycles-per-input-revolution in figure 2.19.
However, for both joints, no detectable error component was found at these frequencies.
This unexpected result can be justified by recognizing that, at any point in time, the
harmonic drive has several teeth in contact due to the concentric nature of its design and the
shape of the wave-generator. Close observation revealed that the harmonic drives tested
typically had up to twenty teeth in contact on both sides of the wave-generator. Given this
result, it is unlikely that the individual error of a single gear tooth could manifest itself in
the aggregate position-error spectrum. Unfortunately, the high-frequency error
components on the joint 1 harmonic drive could not be measured due to resolution
limitations of the output position sensor.
Oeffn = Al sin (cow01 + 01) + A2 sin (0, 202 + 02) + .... (2.7)
54
Chaoter 2: Kinematic ErrorMeasurement
modeling efforts to be discussed later, I relied mainly on the two largest error components,
which appeared at twice- and four-times every input revolution. The resulting harmonic-
series representation of this error is
where 0 inis the rotation of the harmonic-drive input. To be truly accurate, as indicated in
the high-resolution tests presented above, two distinct error components should be included
for each of the two terms in this series to account for the individual error contributions of
the flexspline and circular spline which occur at slightly different frequencies. However,
since it is unlikely that this small change will significantly influence performance of a
dynamic model, this effect will be ignored.
55
Chapter 2: Stiffness Tesiine
56
Chapter 2: Stiffness Tesline
drive which is typical of the apparatus for the other two joints. The testing procedure used
on each joint is as follows.
[•] Read and store output rotation data from the output position sensor and output
torque data from the torque sensor.
r4] Continue collecting data for a range of positive and negative applied loads
until the stiffness curve is complete.
This straightforward procedure was used to collect multiple sets of stiffness data on each
transmission and generate the results presented below.
Figures 2.21, 2.22. and 2.23 illustrate the measured stiffness curves on the three
harmonic drives. Due to torque-sensor loading limitations, the joint I harmonic drive was
loaded to 35 percent of its rated maximum torque, the joint 2 harmonic drive was loaded to
45 percent of its rated maximum torque, and the joint 3 harmonic drive was loaded to 90
percent of its rated maximum torque. These torque limitations, however, encompass the
entire range of dynamic behavior seen by the harmonic drive in experiments to be discussed
later. All three stiffness profiles exhibit noise due to limitations on position sensor
resolution.
57
Chapter 2: Stiffness TestinR
0.08
....... . .4................
.2 0
-0.04 ....
111
-0.08 . ..... I~ ' ......... ... .......
.......
035
0
0 .2... ........... . ..........
......... ............ ........
..... ...
-0.2
-0.31
-042-
-50....-40 -3.-0. 100.020304050607
-7-60...
......... utpu Torque .(.....N
).....
.. ..................
-0.3S
..
..
.. ...
...
5. ....
....
.......
....
....
.....
.....8..
...
Chapter 2: Stiffness Testine
0.3
.......... i........ i........ ....... ....... ........ ........ ........ ....... • ....... .. ....
• ........
i ............ ... .... ...
........ ........
S........
;...... •..... .... ...............
,.... ....... • ....... ...............
. . . ... ;.... ....... ....... ........ ........
•-• .......:... .... " .... ........ ."....... " ... .... . ..... ..$ ..... .. ........-.'........ i........ i....... -.
....... ..........- ....... .........
.......... ÷ ....... "............... . ........ ......... ....... i ....... .. ....... ..
....... ...............
........
...- .. ....... ....... .......
4.... ...........
-0 1 .... .. ......... .........i ....... . .i.. 1....
... ............................... ........ ........ •....... ........ ............
-0.3
-15 -10 -5 0 5 10 15
Output Torque (N'm)
Another set of stiffness tests that I performed were aimed at gaining a qualitative
understanding of how stiffness varied with environmental and operating conditions. In
particular, I used the joint 1 testing station to collect data for a range of different preloads
and wave-generator angles. Figure 2.24 shows three stiffness curves generated from data
59
Chapter 2: Stiffness Testing
0.1 I11V
de r e
0............
0.09...... .-..........
0.08 60 degrees
~, 0.06 .....
0 0.05 . ...
0.02
0
0 10 20 30 40 50 60 70 80 90 100
Output Torque (N~m)
0.1
0.09 low preload ....... .............. .... ....
008edium preload
0.0
0 .01 ....................................
0 0.0 . ..... * 0 ~ 0
06
Chapter 2: Stiffness Testing
collected with the wave-generator major-axis locked in different positions relative to the
circular spline. As the difference between these curves indicates, variations in stiffness of
up to 25 percent may be measured due to different orientations of the wave-generator.
Figure 2.25 illustrates three different stiffness curves collected when the preload between
the harmonic drive gear teeth was varied by changing the depth of the wave-generator in the
flexspline. While wave-generator depth is not a direct indicator of gear-tooth preload, it
was observed to influence the radial loading of the harmonic drive in many cases. The test
labeled "high preload" corresponds to the case where the wave-generator was pushed into
the flexspline to its maximum depth while the low-preload case saw the wave-generator at
minimum depth. From the stiffness curves in this plot, it can be seen that slight variations
in preload can also influence the value of harmonic-drive stiffness.
In general, from the range of behavior illustrated in the stiffness results, it can be
concluded that the variation in stiffness profiles between different harmonic drives can be
substantial. As seen on joints 1 and 2, one drive can behave relatively linearly in the low-
to-medium torque range while another can exhibit a strikingly non-linear profile over the
same loading range. Additionally, the presence of hysteresis due to phenomenon such as
friction and soft wind-up may be dominant in one harmonic drive while unnoticeable in
another. Furthermore, the substantial variations in harmonic drive flexibility due to such
factors as preloading and input orientation indicate the sensitivity of stiffness
measurements. This gloomy forecast for the predictability of harmonic-drive flexibility
grows more overcast when these results are compared to manufacturer's quoted values.
61
Chapter 2: Stiffness Testing'
0.35/
.......................... . . . .
.................. . .
0 10 20 30 40 50 60 70 80 90 100
Output Torque (N~m)
CP
0
~0 .2 ......... .......... ..............
.............
0 10 20 30 40 50 60 70
Output Torque (N*m)
62
Chapter 2: Stiffness Testing
0.35
..... ...........-..........
. .........I ......... . .........
L . ........
... ......... ............ :.........-......... ...........
.. .........
. .......... ........
" .......... "......... :"......... .. .
i.......... :......... [......... ........... •.......... ..".........i .......... ......
0. ...i...... ....
0.3
S.......-.......... .......... .......... i.........- .
:.........
-'i .,....... :... . L....... ......... ii........
.......... .... ........i
0•.~~
. .........
.......
.........
• •.........
-• •,......... .........
........... . ~...
M,
• ....
0 .25 ~~ .......... .........
.....
~ ~" ~ ~ ~ ~ ,...... ....
...................... ........ .......... ..........
. ....... : :.:... .:...... ......... i......... ......... .......... !......
.. .z.....
0 ........ i........... ... . ........
CL
20.25 o............
- ..........
... .. .........
....
..
.........
..........
.-.......
............ .........
.........
.........
.........
.........
i.....
at..........a...
S2 o
.0
05
S0.1
....... .........
........
i......... ...i.....i.....i.............
;' .....i.
i.......... ...... oexperimental stiffness data
0 5 10 15
Output Torque (N'm)
standard stiffness curve. As noted in past research, this unexpectedly high stiffness can
probably be explained by a high gear-tooth preload in the harmonic drive. The lack of
hysteresis loss observed on joint 1 also supports this postulation siice a higher preload
would eliminate any soft wind-up in the transmission. On joint 2, the measured stiffness
while loading the harmonic drive is slightly lower than the standard stiffness approximation
while the optimal quoted stiffness agrees closely with the unloading portion of the
measured stiffness curve. On joint 3, both the loading and unloading portions of the
measured stiffness curve have a slope that agrees with the standard stiffness approximation
within about 10 percent. In partial confirmation of the discrepancies between predicted and
measured stiffness values, manufacturer's catalog [3] warns that a 30 percent variation
between quoted and actual stiffness values is not uncommon. In terms of hysteresis,
catalog [31 states that hysteresis loss rarely exceeds 0.033 degrees of output torsion and
that soft wind-up is kept below 0.05 output degrees of rotation. This prediction holds true
for joints 1 and 3 but is greatly under'~stimated for the joint 2 harmonic drive.
This comparison between experimental data and catalog predictions exposes the
shortcomings of manufacturer's estimates. Multiple stiffness curves collected on each
harmonic-drive assembly illustrated that experimental results were fairly repeatable to
63
Chapter 2: Stiffness Testing
within about ten or fifteen percent. However, measured values of stiffness and hysteresis
loss differed unpredictably and significantly fr )m quoted values. In a few instances such
as the joint I stiftness profile and the joint 2 hysteresis loss, measured values even fell
outside of the manufacturer's predictability tolerances.
Figures 2.29, 2.30, and 2.31 show cubic approximations to the stiffness curves for
the three harmonic drives. In all cases, the cubic polynomial was fit to the portion of the
stiffness data collected during joint loading rather than unloading. From these results, it
can be seen that cubic approximations can accurately capture the shape of experimental
stiffness data. Based on this result, cubic representations of non-linear harmonic-drive
stiffness will be used in dynamic modeling efforts to be discussed later.
.....!.............
..........!.......
0.1 0.1 .......!.. . . . . . . . . .. .. .. .. .
,. . . . . . . . !.. ! !.. . . . . . . ......
0.090.0. ... . . . . . . . . . . . . ... . . . . . . ,. . . . . . . . . . . . . . . . L-
. . . . . . . . . . . . . . . ...............
....... ......
. •............ ...-..
S0 .0 8 ........... ...-............... •......... ...... ............... .. ................
•...............
.......... -'-- --•.".......... ....
S0
0 ... ................ .:............... ......................... .. .. .. .... ................... ............... ........... .
0.08 ........
............. ............... i............... ..
. ....... ........... .. .............. ............ .•. ......... . ........ . .
.........
-~0.06 .....
S. ............. • ............... ............... ; . . . . . .............. o!.. . ..........
. . . ......... ..... . . . ..i...............-. ..............
0 S.
.0 ............
.20 .05 ....... -- .......................
:f.......
. . .•. . . . ... .•
i.. ...... ............ .. . .........i.........
..... . .. ........ .. ...............
.....
... ........ ............... • ................ .. . . . .. .•. . . ... .. ...
........ •.... .
.
•..........
0 ........ . ............ ;...
0 .0 ............. -..
,. • ............... i.............. ............... i............... i..............
........... ............ ....
S0.04
00 ,..... ........... i...... -.........i............... ............... ............... i............... i................. ............... i..............
0 10 20 30 40 50 60 70 80 90 100
Output Torque (N'm)
64
Chapter 2: Stiffness Testing
. .. . . . . . . ... . .. . . . .. . . . . . .. . . . . .... . . . .. .. .
.90.25 ..........
.-....... .....
0.05................ ......
0 .2 5 ..... .....
..... ...........
4..... ........... .....
(Z
S0.2 :W-
....
Is 0.15
0.1.......
.. .......... ........... T..... ..... .......... .....
.....
0
0 10 3 4 100 715
Output Torque (N-m)
65a
Chapter 2: Stiffness TestinR
hysteresis loss, the measured stiffness profiles can be used to make estimates of the friction
present in harmonic drives. As discussed later in this document, dynamic models of
harmonic drives were developed which capitalized on the friction information stored in the
hysteresis loss of the stiffness measurements.
66
Chapter 2: Stiffness Testing
The most significant conclusion that can be derived from the results in this section
in that accurate measurement and prediction of harmonic-drive stiffness profiles can be very
difficult. Not only are catalog estimates likely to be significantly different from actual
values, but experimental measurement can also vary substantially as a function of
operational conditions such as wave-generator angle and preload. To foreshadow issues
which appear later in this document during dynamic modeling discussions, the static and
coulomb friction which are manifested as hysteresis loss in stiffness measurements
compound the problem of extracting accurate stiffness information. Since these friction
phenomena are unpredictable and difficult to separate from pure transmission stiffness, the
accuracy and usefulness of harmonic drive-stiffness measurements becomes even more
nebulous.
67
Measurement
C 2!ue
The procedure used for these experiments was designed to collect accurate
measurements of starting torque in order to understand the amount of static friction present
in a stationary harmonic drive. This procedure is listed below.
[•] Gradually apply a current to the input motor while collecting input torque data
from the current sensor, and position data from the input encoder.
[n] Calculate the harmonic-drive starting torque by subtracting the catalog value
for the motor starting torque from the maximum motor torque reached just
before the joint began to rotate.
The total starting torque in each harmonic-drive testing station is composed of the
individual starting torques due to (1) the input motor, (2) the harmonic drive, and (3) the
output bearings. Since the starting torque exerted by the output bearings, as measured
from the input, will be very small, accurate values for the harmonic-drive starting torque
can be found simply by subtracting the motor starting-torque from the measured value.
68
Chapter 2: StartinR-Toraue Measurement
Table 2.9 lists the measured starting torques for all three harmonic drives. These
starting-torque measurements, however, were observed to vary considerably with the
rotational position of the harmonic-drive output (flexspline or circular spline). These
variations, as listed in table 2.9, can be explained by noting that any misalignment of the
flexspline or circular spline, for example, is likely to influence gear-tooth meshing and,
consequently, gear-tooth friction over one output revolution. Since these starting-torque
values provide a direct measurement of static friction, the substantial variation in starting
torque with output position illustrates the sensitivity of static friction in harmonic drives.
To further support this unfortunate consequence, additional observations revealed a strong
dependency between harmonic-drive starting torque and operating temperature. In a few
cases, if testing equipment was not warmed-up before data collection, starting torques
could increased by up to 80 percent. This troublesome effect can be explained partially by
the increase in lubricant viscosity which increases interface friction at lower temperatures.
Additional variables which can constrain the accurate measurement and prediction of
harmonic-drive static friction are wave-generator orientation and gear-tooth preload. As
noted in previous research, since the orientation of the wave-generator as well as the
applied preload can influence the tightness-of-fit between the flexspline and circular spline,
static friction has been observed to waver considerably under the influence of these two
factors. For example, experimental tests presented in [65] showed that starting torque
could vary as much as 50 percent over one rotation of the wave-generator. Due to the
limitations of the testing equipment, I was not able to confirm this result, however
qualitative observations clearly demonstrated the importance of input angle and preload on
starting-torque measurement.
69
Chapter 2: Startin,-Torque Measurement
TablInuttartng-
2.: orqevluesfo th ithe hamoi doivtes
Since, by definition, static friction occurs only when velocity is zero, developing an
accurate model of this phenomenon is not critical for predicting harmonic-drive behavior at
non-zero velocities. However, by assuming that velocity-independent friction, such as
Coulomb friction, behaves similarly to static friction, conclusions can be drawn about
friction components at non-zero velocities. For example, due to the preload between the
harmonic-drive gear teeth, a non-zero contact and rubbing force will always exist. From
the definition of Coulomb friction, this non-zero normal force should always generate
proportional frictional losses at all speeds of harmonic-drive operation. A truly accurate
model of this coulomb friction behavior should use the harmonic-drive geometry to derive a
value for the gear-tooth preload which can then be used to calculate the resulting Coulomb
friction. Alternatively, for the sake of simplicity, these frictional losses due to preloading
can be approximated by either a constant frictional torque or a function which replicates the
typical behavior of the Coulomb friction. By assuming that the static friction results have a
magnitude and behavior similar to this velocity-independent friction function, estimates of
harmonic-drive frictional losses at non-zero operating velocities due to Coulomb friction
70
Chapter 2: Starting-Torgue Measurement
can be made. These insights about velocity-independent friction gained from starting-
torque measurements will also be supplemented by the dynamic damping measurements
presented below.
71
Chapter 2: Dynamic Friction Measurement
Assuming that the rotational velocity during data collection was truly constant, the input
motor torque should exactly balance the damping torque at the given velocity. Therefore,
by collecting input torque data over a range of different velocities, the damping profile for
each harmonic drive was constructed.
72
Chapter 2: Dynamic Friction Measurement
By using the procedure described above, the input motor torque required to operate
the harmonic drive at a constant velocity was used to determine the damping torque of the
transmission at different velocities. Due to the harmonic-drive kinematic error, the
equilibrium velocity reached by the harmonic drive in many trials experienced noticeable
fluctuation. In these cases, the average equilibrium velocity was used to determine the
damping torque. The resulting torque and velocity data that was collected on each testing
station directly characterized the total friction torque in the system. This total damping
torque included not only the friction present in the harmonic drive but also the damping
pr::;ent in the input motor and output bearings as well. By assuming that, since the
harmonic drive output rotates at relatively small velocities, the damping contribution from
the output bearings was negligible, the frictional torque due solely to the harmonic-drive
transmission was calculated by subtracting the motor damping from the measured aggregate
friction-torque. Figures 2.32, 2.34, and 2.36 illustrate the total measured friction-torque
plotted with the motor damping, while figures 2.33, 2.35, and 2.37 show the extracted
harmonic drive friction-torque component.
From the six plots shown below, several observations about dynamic friction
losses in harmonic drives can be made. First, as figures 2.32, 2.34, and 2.36 illustrate,
the friction torque generated in the harmonic drive is much larger than any other damping
losses in the testing equipment. This anticipated result pronounces the importance of
accounting for transmission losses in harmonic-drive models and calculations. Second, as
labeled in figures 2.33, 2.35, and 2.37, each drive passes through operating regions in
which frictional losses are amplified due to dynamic resonance vibration. More
specifically, at operating velocities where the position error of the harmonic drive excites
joint resonances, the resulting torque fluctuations can increase friction in the transmission
73
Chapter 2: Dvnamic Friction Measurement
0.3
0.25 .-...
totalfrictionto
....... ... .... ...-....
... .............. .... .......?......
1.......
.................... ...... ...
0.
.... .. . ... .. . . . .. .. . . . .. . . . .. . . . . ..-. . . . . . . .. . . . . . .. . . . . . . . . . . .. . . .
0.25
02..
......
..................... ......
.. ......
............ ......
......
0 .2 ... ...... ............
......
..... .....
......... .... ...
E0. ..... ...... ........I
... ................... ...... ...... .............. .....
74
Chapter 2: Dynamic Friction Measurement
.0.12
....
*..
....
..... total friction torquel
(D 0.06 ~ . . ......-
....
... ........ . ......
......
. ... ... ..
. . . . . -. . .-.
0.12
0.
-0.04(
0
0 10000 20000 30000 40000 50000
Input Velocity (deg/sec)
75
Chapter 2: Dynamic Friction Measurement
0.015
.......
............
... ... ...... .. torq....
..... ue..
....
:...
190.0151
* ~ ~. .
. .0 . *.
. . . . .
0 0
0-0
.. . .,. .. . . . ..t .. . ....
!... .. .... ... ... .. ..
0.005 tru
... L...
.......
L.. .... . ............. .J.......
. .... ]........J.....
. ..L
......... ......
....
10000
0...... 20000. 30000. 40000 50000
..
Inu reloiy desc)
.... ..........
............... ...
76.. .......
...
Chapter 2: Dynamic Friction Measurement
and dissipate extra energy. Depending on the location of these resonances as dictated 'y
the design of the harmonic-drive system, the frictional losses excited in these regimes are
likely to dominate harmonic-drive behavior for a significant range of operating velocities.
Third, the shape of the damping curve and the torque-offset experienced at low rotational
velocities can lend insight into the nature of the non-zero, velocity-independent friction-
torque resulting from the harmonic-drive gear-tooth preload. As discussed in the previous
section, Coulomb friction at the gear-tooth interface due to non-zero gear-tooth loading
influences the starting-torque of the harmonic drive as well as friction-torque offset
experienced at non-zero velocities. Rough estimates of this velocity-independent offset-
torque are derived from the damping-curve profiles in figures 2.33, 2.35, and 2.37 and are
labeled "constant friction torque". As explained previously, these frictional offsets may be
constant with respect to velocity but are likely to be influenced by other parameters such as,
input and output rotational angle, operating temperature, and preload.
Harmonic drive catalogs, [3, 4], characterize the dynamic losses experienced in the
transmission by providing efficiency ratings at different operating velocities. These
efficiency ratings are derived by comparing the friction torque in the transmission to the
rated torque of the harmonic drive at different operating velocities. Harmonic-drive
catalogs specify the rated torque for each drive by using bearing-life equations to calculate
the operating torques required to ensure a constant wave-generator bearing life at all
operating velocities. Using this definition of the rated torque and the experimental frictional
77
Chapter 2: Dynamic Friction Measurement
* 0.5 -....
2. ... .......
... 3 . . ..........
0 . 4.... . .....
.......................... ......
....... .......... ... ......... ..... ... ..
............... ...... .. .. .
S 0.3 -..... ........................ ...................-.......
cc
~~ ~ ~ ~ ~ ~ ~ a.- ........ ...... .....
....... eff....ic.....iency..
measured................
0.0 2 ................-.... ..............
. e..ic.enc
...... catalog..........
.......
....... ..................... ...... ............................... .......... ......
0
S 0.3 :...............
. ...... ....... ......
.........
MA ........ . - o.....-..
........... measured efficienc ........
0 .2 . ....
I................................ .....
S~. catalog efficiency~*
' 01 ...... i
0 100 800020 000 30000 40000 50000
Input Velocity (deg/sec)
1 78
Chapter 2: Dynamic Friction Measurement
1 r r F • + i r r i - - i- - r + •
............
......-......... ' .....................
"..."...."."".""
a,) 0.9 ..
2' 7 .... ~~~..
... ;"" ......
.................. " ' ; ...... . ..
0.8
,, 0.5
- +• -•......................... ....
... .- . + . .. • .i..+......-- -
....
05
~.. .....
uJ 03
"".3 --...... i.. .+..............+.........-• ...... ........-......
• ... ..- •
.....
-....i .. .....
.... .. + ] _....• _
~~~
.~
I0-~ i . .. . . . . . .. ..... ........ ............. . ... i. . . . ..... . ..... . . . .. . . . .. .. . .
: ......
i. ....... ....... ....... . ....
o.. ....... ............ +
... . . •..+•.••...•....•..].. .•; .. .•............
losses, the resulting transmission efficiency was calculated for each harmonic drive.
Figures 2.38, 2.39, and 2.40 illustrate these experimental efficiencies as well as the catalog
efficiency curves from [41. Given that the catalog values are provided as rough estimates,
these three plots demonstrate that the catalog efficiencies are relatively similar to the
experimental results. For all three experimental curves, the dramatic efficiency loss that
occurs during resonance is apparent, and, since the location of these system resonances is
dependent on harmonic-drive compliance as well as the inertias of the specific harmonic-
drive system, it is not expected that catalog efficiency curves should capture this behavior.
Additional discrepancies between experimental and predicted results can be accounted for
by a number of environmental factors such lubrication or operating te+mperature.
79
Chaoter 2: Dynamic Friction Measurement
gear-tooth preload. In its simplest form, this component can be represented by a constant
torque at all velocities, but, as discussed above, a truly accurate model should incorporate
influences from such factors as input and output orientation and gear-tooth preload. The
s5',cad friction component, namely non-linear velocity-dependent damping, can be better
understood by ignoring the resonance losses in the damping profiles of figures 2.33, 2.35,
and 2.37. The resulting damping curves, as shown in figures 2.41, 2.42, and 2.43, can be
approximated accurately by cubic functions. Given these two models of harmonic-drive
friction components, the representation will be complete if the behavior of resonance losses
can be captured. Unfortunately, this last friction component is somewhat more difficult to
model since it depends on the dynamic torques experienced by the harmonic drive.
Consequently, a representation which produces these effects must be developed in
conjunction with a complete dynamic model of the transmission as will be discussed later.
I
0.32 S............ .............. __
......... .............. ............. .............. .............. ............ ......................... . .
0 .2 8 -. ........ . .......... .............. i............ ............. ............. !............ ... ... ....... ..
I. .......... ............. t............. ............ ...... ....... ............. ............ ..... ... ........
0.24
E .....
..... ........ ............ 1............
..............
. . . . . . . .. • r . . . ............. ............ .......................... ......................
Z
0 • ... .......... ............ ........ ............................ ........... ..............
. .0Y............ ............. .......... ........... .......
S0.16
01 . .............. ..... ................ ............. ............. .......... ....... ............. ............. ......... ........
. . ................
E
0.1
008.. ...... ....... ....... .. .. .. ..-
....................................... . . ......
04
... .. ....... ... .... .... . . . . .
............. ........
...... ...
04. ...... ............. ............ ............ ......... ............ ......... ............. ............ ............ ............
. . .......................
0.0
80
Chaptecr 2: Dynamic Friction Measurement
o.12 -1-rfI!TIr
11 1 111 1 F
0.0
4M 0.06
00,4
..................... ...... ... ...... ......
............ ................ ......
...
0.0....1I.. 5........ ....
z.
..... .....
.... .... ...
............
........ ..... . ..
........ . . ....
...... .....
0.0............
........ . . ... . . . . .
.. ..
.. ...
...
.... .....
...
..
...... .... ... .
.........
Chapter 2: Dynamic Friction Measurement
torque, (2) a non-linear damping torque, and (3) an additional torque loss during system
resonance. Useful measurements for the first two of these frictional components can be
made by analyzing the shape of the damping profile of a given harmonic drive while
accurate representation of the third friction mechanism is more elusive. Due to the
dependence of transmission friction on environmental and operating conditions, catalog
efficiency curves can only be used to get a rough approximation of harmonic-drive losses at
different velocities.
82
Chapter 2: Dynamic Response Measurement
83
Chapter 2: Dynamic Response Measurement
Collect data from the position sensors, current sensor, and torque sensor for
several seconds of joint operation.
n5] Derive velocity and kinematic error response from the position-sensor data.
The motor-current commands used in each trial ranged from the lowest current that could
overcome transmission starting-torque to the current-saturation limit of the motor
amplifiers. The resulting velocity effected by these current commands spanned the typical
recommended operating range for the three harmonic drives. In all cases, current step-
84
Chapter 2: Dynamic Response Measurement
-*
- - --- -- I
Xo
....
.. i... n ..:::.Joint 3
3.8e-5 kg-m 2
......... 2.8e-5 kg-rn 2 2.0e-6 kg-m 2
tota~il
nput InetIa 2.2e-4 kg-rn 2 7.6e-5 kg-m 2 5.2e-6 kg-M 2
t .......... 2.828 kg 0.7027 kg 0.7027 kg
•' .. ..,.............. ...... . .. .
commands were incremented at regular intervals small enough to capture the variations in
dynamic behavior with increased current. From the resulting myriad of data sets, a handful
were extracted which effectively characterized the dynamic response over the entire
experimental range of operation.
85
Chapter 2: Dynamic Response Measurement
structured by first outlining typical behavioral patterns as observed in the velocity response,
and then discussing the influence of these patterns on torque, current, position and position
error measurements.
F21 Th1ese regions of resonance act as energy sinks which constrain the increase
of average rotational velocity for increasingly larger current commands.
[n] An unpredictable jump in velocity can frequently occur when the transmission
harnesses just enough energy to push through a system resonance.
[4] Due to variations in frictional resistance over one output revolution of the
transmission, significant variation in rotational velocity can result.
86
Chapter 2: Dynamic Response Measurement
These surprising features which are common to the measured response of all three
harmonic-drives will be discussed in order.
From the listed time-response illustrations, it can be seen that the input velocity of
the three harmonic drives passes through regions where significant resonance vibration can
occur. These resonant regions occur when the frequency of kinematic-error fluctuations,
which occurs at multiples of input, or wave-generator, rotational speeds, coincides with the
primary temporal resonance of the harmonic-drive system. Since the harmonic-drive
transmission flexibility behaves non-linearly, there is no distinct frequency of system
resonance. Instead, as the torque-fluctuation amplitude increases, the effective stiffness of
the transmission increases which results in an increase of the frequency of system
resonance. This variation in temporal resonant behavior is manifested in a range of
operating velocities which can excite system resonance. On joint 1, two separate
resonances can be observed at input speeds of around 4000 and 8000 degrees/second. The
vibration in the first of these two regions can be attributed to the kinematic-error excitation
which occurs four times every input revolution and the second region can be blamed on the
twice-per-input-revolution position-error component. As implied, the resulting frequency
of velocity fluctuation that occurs at each of these resonances is identical to the rotational
speed of the harmonic-drive input times the spatial frequency of the responsible kinematic-
error excitation. On joint 2, the amplitude fluctuations at resonance are smaller than
observed on joint 1, but regions of resonant behavior can still be found at about 10000
degrees/second and , to a lesser extent, 20000 degrees/second. These two regimes can be
attributed to the kinematic-error fluctuations occurring at 2.0 and 1.0 cycles-per-input-
revolution, respectively. On joint 3, three resonant areas can be observed at the
approximate input velocities of 5000, 10000, and 20000 degrees/second, and can be
directly associated with the kinemaLic error occurring at 4.0, 2.0, and 1.0 cycles-per-input-
revolution, respectively. The lack of any significant system resonance on joint 1
corresponding to a frequency of once-per-input-revolution is not surprising since, as
illustrated previously in figure 2.10, no appreciable kinematic error component was
measured at that frequency. Additionally, kinematic-error harmonics appearing at
frequencies higher than 4.0 cycles-per-input-revolution have no effect on any of the
observed dynamic responses since the resulting excitation frequencies exceeded the system
resonant regions at the slowest experimental operating velocities.
87
Chapter 2: Dynamic Response Measurement
roughly the same final rotational velocity at the first region of resonance, and current
commands ranging from 2.8 to 4.0 amps all produce the same final velocity in the second
resonance regime. Similarly, as illustrated in figures D.2.2 and D.3.2, identical behavior
can be observed around the resonance regimes on joints 2 and 3. The velocity limitation at
resonance can be attributed to the dramatic increase in Coulomb-like friction at system
resonance due to the large torque fluctuations experienced in the transmission. In other
words, as a step command in current brings the system closer to resonant velocity, the
torque fluctuations due to resonance begin to increase. With this increase in torque
amplitude, energy is increasingly dissipated by a Coulomb-like friction mechanism
presumably at the gear-teeth interface. Due to this increasing energy loss as the velocity
nears the resonant areas, fluctuations in velocity increase but the average velocity remains
roughly the same. Because of this self-reinforcing frictional mechanism, the harmonic-
drive velocity is prolonged at the frequency at which kinematic error fluctuations excite
system resonance. This problem is heightened by the non-linear nature of the stiffness
since increases in torque fluctuation as the system approaches resonance also increase the
natural frequency of vibration. This effect further prevents the input current from pushing
through the harmonic-drive resonances since the temporal resonant frequency increases
alongside the increase in the rotational velocity and resulting spatial excitation frequency, of
the harmonic drive. In summary, because of the escalating frictional losses and the
increase in natural frequency which occur as rotational velocities approach system
resonance, a wide range of input current steps can fall victim to these resonance traps. As
an explanatory note, the velocity truncation present in the response of the highest current
step on each velocity plot is due to the saturation limits of the current amplifiers and should
not be attributed to resonance losses.
88
Chapter 2: Dynamic Response Measurement
harmonic drive, a resonance singularity often appears when the frequency of resonance
excitation slightly exceeds the natural frequency of vibration. More specifically, as the
amplitude of torque vibration increases, the equivalent stiffness of the harmonic drive also
increases as does the natural frequency of system vibration resulting in a frequency
spectrum similar to the one show in figure 2.45. However, when the harmonic-drive
velocity breaks through this resonant region, resonance vibration drops immediately due to
the unstable decline in both torque-vibration amplitude and natural frequency of resonance.
This non-linear dynamic behavior is commonly referred to as a jump resonance. The
second accomplice in this jump-resonance behavior is the frictional loss at resonance. Due
to these losses, as discussed in the previous paragraph, velocity is restricted in resonant
regions. However, when the harmonic-drive jumps to non-resonating behavior, the
resonance losses disappear and the leftover energy increases rotational velocity
dramatically. If the non-linear stiffness and resonance losses were the only factors
influencing jump-resonance behavior, it is likely that this unpredictable response would
only appear a very small range of input current commands. However, as observed in the
experimental results, these jump-resonant responses occur for a surprisingly large range of
input currents and appear far too often during normal harmonic-drive operation. The origin
of this unfortunate consequence can be found in the large variation in frictional resistance
that can exist in harmonic drives over one output rotation. Due to this variation in energy
dissipation, all operating velocities somewhat near jump-resonance behavior can harness
enough energy to break through the resonance when frictional variations with output
89
Chapter 2: Dynamic Response Measurement
rotation decrease. Depending on the magnitude of these frictional variations, a wide range
of operating velocities can be kicked out of resonance due to this decrease in friction. In
conclusion, due to the interaction of non-linear compliance, resonance losses, and friction
variation, the stability of harmonic drive operation can be seriously undercut by drastic and
sudden jumps in operating velocity.
ai (N
i ) 0.6 1.3 0.03
90
Chapter 2: Dvnamic Response Measurement
behavior. Also note that the beating in output velocity oscillation amplitude at high
velocities seen on joints 1 and 2 is due to aliasing in the sampled data rather than an actual
physical mechanism.
The current response measured from the motor amplifiers is illustrated in figures
D.1.3, D.2.3, and D.3.3. Since the motor current is linearly proportional to the motor
torque, as dictated by the motor torque constant, the behavior of the harmonic drive input
torque can be understood. From the response curves for each joint, it can be seen that at
low current commands, the amplifiers can sustain a fairly constant input torque in the face
of harmonic-drive torque fluctuations. However, at higher current commands, the current
amps have difficulty keeping pace with the higher frequency torque fluctuations and
noticeable variations in current amplitude occur. The highest current-command response in
each figures illustrate the vulnerability of the current signal when the amplifier saturates
from the high motor velocity. Due to imperfections in the A to D board, significant
electrical noise appears in the current signals on joints 2 and 3. Fortunately, this noise had
no effect on the dynamic response of the harmonic-drive systems.
The output torque response measured from the torque sensor on each joint is
illustrated in figures D. 1.7, D.2.7, and D.3.7. As a reminder, because of the location of
the torque sensor on each joint, this measured response is a direct reading of the torque
experienced at the flexspline of the harmonic drives. This output torque has two
contributions: (1) the dynamic torques generated by the acceleration of the output inertia,
and (2) the damping torque produced by the output bearings. Since the damping torque in
the output bearings is very small compared to the inertial torques generated by the position-
error fluctuations, the measured torque response in all cases oscillated around zero at all
velocities. As shown on all three joints, torque fluctuations can become extremely large at
velocities which excite system resonance. For example, on the largest harmonic-drive
system, joint 1, torque readings reached a powerful 80 N-m during system resonance. As
expected, the frequency of torque fluctuation is identical to the frequency of velocity
fluctuation which is equal to the rotational speed of the input times the position-error
excitation frequency at the appropriate resonance.
The input and output position response measured on the three harmonic-drive
specimens is illustrated in figures D. 1.1, D.2.1, and D.3.1, and figures D. 1.4, D.2.4, and
D.3.4. From these plots, it can be seen that substantial fluctuation noticed on the velocity
plots is virtually unnoticeable on position measurements. However, due to increased
frictional losses at resonances, the increases in final position for constant increases in the
91
Chapter 2: Dynamic Response Measurement
current-ste, input is by no means constant. For all tests, the input and output position
responses .splay identical behavior except that they differ in magnitude by the given
harmonic-drive reduction ratio.
In general, changes in the dynamic response for identica' step-response trials was
minimal as long as harmonic-drive parameters remained unchanged. For example, since,
as discussed earlier, factors such as harmonic-drive preload, assembly, lubrication, input
and output orientation, and operating temperature, can greatly influence system parameters
such as stiffness and frictional losses, changes in any of these conditions between trials can
greatly influence dynamic response. Care was taken to make sure that these parameters
remained consistent across all step-response trials on each joint.
92
Chapter 2: Dynamic Response Measurement
measured relatively easily and modeled by mathematical functions which approximate the
observed behavior. The second item, t -)wever, relies on the interaction between torque
fluctuations in the harmonic-drive and rubbing at the gear-tooth interface. Therefore, any
model that is supposed to accurately describe frictional losses at resonance must include a
representation which predicts the dynamic forces acting on the gear-tooth geometry.
Additionally, as illustrated in the motor-current response curves, a dynamic model of the
current amplifiers might be necessary in some domains to capture the saturation behavior
and dynamic fluctuations.
93
94
The goal of modeling any physical mechanism is to discover the simplest
representation which can replicate system performance to the desired level of accuracy.
When modeling a mechanical assembly with a harmonic-drive transmission, several levels
of detail can be incorporated into the harmonic-drive model to improve the integrity of the
system representation. In particular, as seen in the experimental results, harmonic-drive
properties which can significantly influence system performance include
The goal of this section is to characterize the modeling accuracy that can be achieved by
incorporating different combinations of these transmission properties into a harmonic-drive
representation. From the complex behavior observed in the measured dynamic-response
data, finding a model which accurately characterizes harmonic-drive behavior proves to be
a formidable task. However, if such a representation can be created, design and control of
systems containing harmonic drives can be substantially improved. Additionally, the
knowledge embedded in such a dynamic model can become a powerful medium for
understanding and improving the design and operation of harmonic-drive transmissions. I
hope that the results presented in this section can equip and enlighten the expanding
community of harmonic-drive users and designers.
95
Chapter 3: Theoretical ModelinR
theoretical predictions. To address this issue, compliance and then kinematic error were
included into the transmission model. This improved representation produced simulated
results which more closely matched the experimental observations but failed to capture the
significant energy losses which occur at system resonance. To eliminate this discrepancy,
a final model which incorporated gear-tooth geometry and coulomb friction was developed
which improved predicted performance around resonance. Unfortunately, due to the
sensitivity of this model to harmonic-drive properties and the inherent difficulty in
accurately measuring these properties experimentally, further research must be done to
refine this representation into a viable prediction tool. Nevertheless, the novel framework
developed in this model provides a simple and powerful representation for hatching and
nurturing an intuitive understanding of harmonic-drive operation.
96
Chapter 3: Modeling Techniques and Notation Conventions
F5] All positions and velocities are defined to be positive in the clockwise
direction when viewed from the input-side of the harmonic drive.
F61 All torques are defined to be positive in the direction in which they will act
when a positive load is applied to the input inertia and the output inertia
remains stationary.
F] lThe positions, velocities, and torques seen by the three ports of the harmonic
drive will have the subscripts wg for the wave-generator, fs for the flexspline,
and cs for the circular spline.
Additional variables such as model parameters and constants that are used for the testing-
station models are defined in table 3.1. Similarly, table 3.2 lists descriptions of the
variables which will appear in the five different harmonic-drive models to be developed
later. Sporting this notation lifejacket, I will now plunge into the depths of equation
derivation and analysis for five different harmonic-drive models.
97
Chapter 3: ModelinR Techniques and Notation Conventions
station models, five different harmonic-drive models will be developed that can be inserted
easily into the overall system models. Since the inertias of the wave-generators are
included with the motor armature inertias, no inertia elements are needed in any of these
harmonic-drive representations. Consequently, the equations of motion for each will
consist solely of algebraic rather than differential equations. These algebraic equations will
98
Chapter 3: Modeline Techniques and Notation Conventions
outline the relationships between the torques and rotations on the three harmonic-drive
ports. In particular, since the inertias in the testing stations dictate the position and velocity
seen at the input and output of the harmonic-drive, equations will be developed to relate
these rotations to the torques experienced by the transmission. Additionally, because these
equations completely characterize the three-port behavior of the harmonic-drive, any two-
port configuration of the transmission can be represewted easily by setting the rotation of
the non-rotating harmonic-drive port to zero.
Due to the complexity of the differential equations generated when the harmonic-
drive models are inserted into the overall system models, analytical solutions describing
dynamic behavior were rarely available. However, using a Runga-Kutta numerical
integration scheme with adaptive step-size control, the dynamic behavior specified by a set
of differential equations could be simulated over a given time period for a specific set of
input conditions. The code which was developed to generate dynamic time-response data
for different harmonic-drive models is listed in Appendix B. In addition to calculating the
forces, positions, and velocities predicted by the differential equations, this simulation code
determines the energy distribution in the system and can be used generate static stiffness
data as well.
When executed on a Sun SPARCstation ELC, typical simulation times for a four-
second time-response using the more complex harmonic-drive models were much less than
one minute. Although the exact time for simulation is influenced not only by the
complexity of the dynamic equations but by the desired accuracy and nature of the system
parameters as well, in general, the small computational requirements of the harmonic-drive
models presented below suit them well for rapid numerical simulation. Additionally, it is
likely that simulation performance can be further improved through refinement of the code
which calculates values using the dynamic equations.
99
Chapter 3: Modeling the Experimental Apparatus
The relative location of important model parameters can be understood from these
figures. First, for both joint configurations, the DC motor applies a torque to the input
inertia which is composed of the motor armature and the harmonic-drive wave-generator.
Additionally, each joint sports an output inertia consisting of a large mass mounted to the
end of the output link on the output port of the transmission. Due to motor damping, a
friction torque can be seen between the motor shaft and the motor housing, and similarly,
the output bearings introduce a finite friction torque between output rotation and ground.
Since, the motor housing on joints 1 and 2 rotates with the circular spline, the absolute
rotation seen by the input inertia is equal to the sum of the input and output rotation. The
input inertia on joint 3, on the other hand, sees the direct rotation of the motor shaft relative
to ground.
100
Chapter 3: Modeling the Experimental Apparatus
input damping
input motor torque
input rotation
input inertia
[DF-harmonic drive
output rotation
output inertiLaL_ýýý
output damping
E3
output rotation
output inertia -7:
101
Chapter 3: ModelinR the Experimental Apparatus
-I-
KIM1,,. 0 - 2 *114
102
Chapter 3: ModelinR the Experimental Apparatus
-03
103
Chapter: Modeling the Experimental Apparatus
Using the lumped-parameter models illustrated and described above, the generation
of dynamic equations is straightforward. The variable names shown for both models in
figures 3.5 and 3.6 will be used throughout this derivation and governed by the notation
conventions itemized in section 3.1. For further clarity, important model parameters were
summarized previously in table 3.1. Using the given notation, naming conventions, and
parameter definitions, the equations of motion for joints 1 and 2 will now be derived
followed by a similar derivation for joint 3.
0o. T..
Tb-., = f(O)
104
Chapter 3: ModelinR the Experimental Apparatus
Tb Wt= f(w
T,
By applying Newton's Law to the input and output inertias in the model shown in
figure 3.3, the resulting equations of motion for joints 1 and 2 are
and
J Tc
t= - Tb_ + Tb in (3.2)
In order to solve these equations, the torques on the right-hand side of each equation must
be expressed in terms of the state variables, OI and Ornt. This can be done readily for the
damping torques by assuming that they are linear functions of velocity:
105
Chapter 3: Modeling the Experimental Apparatus
T ` bout 0
=.out out • (3.4)
Since the motor torque is linearly proportional to the current produced by the motor
amplifiers, Tm can also be determined:
Tm = Kt im, (3.5)
where the amplifier current, ir, can be found from the equations presented in Appendix C
which describe the saturation limits of the amplifiers. The remaining unknown torque
components, Twg and T,,, can be found by specifying the kinematic constraints imposed by
the system on the three harmonic-drive ports and calculating the resulting torques using one
of the harmonic-drive models presented later in this section. Since each port of the
harmonic drive is directly coupled to either ground or one of the inertias, these kinematic
constraints are simply
O= Oout - (3.8)
The equations of motion for the joint 3 testing station are very similar to those listed
above for joints 1 and 2. Using figure 3.4, they can be determined readily to be
The output damping torque for this model can still be calculated using equation (3.4)
above, but the input damping torque must be determined from
106
Chapter 3: Modeling the Experimental Apparatus
Additionally, the motor torque expression remains the same as equation (3.5) above, but
the kinematic constraints on the harmonic-drive ports are now
OCS= 0. (3.14)
Having derived the governing equations for the two different harmonic-drive testing
stations, dynamic response can be simulated once values have been determined for the
given model parameters. First, the input inertia can be found by summing the catalog
inertia values for the wave-generator and the motor armature, and the output inertia can be
calculated from the weight and location of the output mass on each joint. Since all other
output inertia components on each joint were negligible compared to the inertia of the
output mass, they were ignored in the output inertia calculations. All important inertia
values are summarized in table 2.10 and reproduced in table 3.3 below. Second, a value
for the input damping coefficient can be obtained from the DC motor catalog, and the
damping in the output bearings can be found from roller-bearing catalog specifications.
2.e47.6e-5 5.2e-6
107
Chapter 3: Modeline the Experimental Apparatus
These values are also listed in table 3.3. Lastly, as mentioned above, the torque produced
by the DC motor is linearly proportional to the current it receives from the amplifiers. This
current is specified by the desired input command and remains true to this value unless high
motu.t velocities cause the current amplifiers to saturate. Using the catalog motor
specifications listed in table 3.3 and the measured motor velocity, appendix C presents a
simple algorithm which can be used to capture saturation behavior. Having developed
representations which describe testing-station operation and determined values for
corresponding model parameters, attention can now be focused directly on the harmonic-
drive transmission.
108
Chapter 3: Harmonic-Drive Model 1
Note that a positive value of Tfs will cause a negative rotation of the flexspline since it is
defined to be positive in the direction opposite to Twg and T, As previously discussed, N
is the catalog gear ratio and the remaining variables are illustrated in figure 3.8. Using
these equations, the positions, velocities, and torques for any harmonic-drive configuration
are completely defined.
109
Chapter 3: Harmonic-Drive Model I
damping torque of the motor and output bearings. This behavior is illustrated in the
simulated input velocity responses for the three harmonic-drive testing stations shown in
figures 3.9, 3.10 and 3.11.
80000
0.....
0.,...
... .......
.....
•.. . . S... . . i.......... .,...
......
...
........ .. ....
. . . . i.. . . . I........
i........... ........
., ,............
....
.......
.........
.....
. . . . . . . . . . . . . . . . i..........
. . . . . . . . . . . . !.. . . . ........
. . . . ........
i........
S60000
. •. . . .......
i...
.....
.......
..i.... .i...
....
....
..i....
.......... . ........ ........ ..
S.... ....
........
........
/........
.....
............... :::
::::
........
:::
. . .. ........ -........ .......... ........
.. .. ... . ........ ........ ........ •....................; .... ...•........,
......
• ...-. .......... ........ ........ ......... ......... ......... , .... ........... ;.. . . ................. ............ ..... .............................. ......... . . . ... . . .
...................
. 20000
" --.-.. ..... ............... .. ........ •........ •.......
........ "........ -........ ........ ........ i........- ....... . . ....i ..... ...-.........i ........i ..... ... -........ •........ ........
2000 . 11
..
Chapter 3: Harmonic-Drive Mode! I
..........................
.0000................
20000
. . . . . .
........... . . . .
2*10000 -. . . . . . . . . . . . .
.. S
8 0 0 -... . ...............
.4 . ........................................... .........
4. ..........
600 .. .-
40000-.........
... ....
.... .........----..............
............. ....
.......... ........ .......
0
0 12 3 4
time (seconds)
120000
Chapter 3: Harmonic-Drive Model I
kinematic error in the transmission model. Given this justification for poor model
performance, it can be concluded that the only circumstances where this representation may
prove useful is when harmonic-drive frictional losses are small compared to other losses in
the mechanical system and when velocity and torque fluctuations are unimportant. For
example, under static loading, when dynamic friction is zero, this model can be used to
accurately predict torques when applied loads are significantly greater than the static friction
in the transmission. Additionally, in static positioning applications where the desired
accuracy is lower than the kinematic error in the transmission, this model may be useful to
calculate the rotational displacement of the harmonic-drive ports.
112
Chapter 3: Harmonic-DriveModel 2
113
Chapter 3: Harmonic-Drive Model 2
In addition to conceptually capturing the rubbing losses between the flexspline and
circular spine gear teeth, the location for the frictional losses shown in figure 3.12 was
selected for several additional reasons. First, as shown in figures D.1.7, D.2.7, and
D.3.7, since torque-sensor readings in experimental dynamic response measurements
display very large fluctuations but virtually no frictional bias torque, it can be concluded
that the torque sensors are blind to any frictional torque present in the harmonic drive. On
joint 1, for example, if the friction losses in the model were placed between the circular
spline and ground, the torque sensor mounted to the flexspline would not only see any
torque fluctuations on the output inertia but would also directly measure the frictional
torque from the transmission losses. By moving these frictional losses to act between the
flexspline and circular-spline, the torque sensors for both joint configurations ignore these
losses and return results that more closely resemble the experimental observations. As
another justification for the selected location of transmission losses, any static friction
included in this dissipation element will be manifested as hysteresis loss in a stiffness
profile collected from either output port of the transmission. On the other hand, if the
friction mechanism occurred between the wave-generator and either output port, hysteresis
loss would appear only in the stiffness profile collected from the output port connected to
that friction element. From this distinction, it seems highly probable that frictional torques
in the harmonic-drive occur primarily between the flexspline and circular spline gear-teeth
and can account for the hysteresis loss observed in experimental stiffness profiles. To
substantiate this claim empirically, results were gathered from models employing different
friction mechanisms and the integrity of the representation derived in this section was
verified.
114
Chapter 3: Harmonic-Drive Model 2
1.T.,T.0)..
0. eTnm i co ft
In these equations, Tb is the aggregate friction-torque function and, as implied by the sign
conventions stated previously, Tnfs is defined to be positive in the opposite direction of
Tn.Cs. The variables used in these equations and the rest of this derivation are illustrated in
figure 3.13. Using the ideal-transmission torque relationship in equation (3.17), the torque
on the wave-generator can be directly related to Tn..f and Tncs:
Given these three equations, the torque constraints on the model are complete,
however, the friction-torque function, Tb, must be defined in order to finish the
representation. As observed in experimental results, harmonic-drive friction typically has
three components: (1) a velocity-independent friction term, (2) a non-linear velocity-
dependent damping function, and (3) a periodic friction torque which varies with output
rotation. Experimental observations also suggest that frictional losses can be greatly
influenced by other parameters such as preload, wave-generator orientation, and operating
temperature. For the sake of model simplicity, these frictional influences will be ignored.
Additionally, experimental results show that harmonic-drives exhibit a substantial static-
friction torque which must be overcome before rotation can occur. Since this stiction
torque is also influenced by several parameters such as preload and input angle, an accurate
model of this property relies heavily on careful experimental characterization. Once this
experimental understanding is complete, a stiction model can be implemented by placing
appropriate constraints on the minimum motor torque required to effect rotation. Although
115
Chapter 3: Harmonic-Drive Model 2
this model was developed in the computer simulation, it will be ignored here since it only
influences dynamic behavior when velocities are close to zero and consequently has little
effect on the experimental and theoretical dynamic response. Given these assumptions, the
friction function which provided optimal results with minimal complexity is defined as
where
As equation (3.21) states, this model captures the three facets of friction behavior noticed in
the experimental results. First, the velocity-independent friction component, as described
in equation (3.22), is represented by a single value, bcontn. Second, since experimental
data indicated that velocity-dependent friction could be accurately captured by a cubic
relationship, the dynamic-friction torque defined in equation (3.23) is a cubic function of
the relative velocity between the flexspline and circular spline. In this equation b, and b2
represent the linear and cubic damping coefficients, respectively. Since all observed
experimental behavior indicates that dynamic damping is highly non-linear, a linear
damping representation is likely to cause inaccurate model performance. Third, periodic
variations in the frictional torque are captured in the sinusoidal function in equation (3.24).
As this relationship indicates, frictional torque fluctuations of half-amplitude Ab complete
one cycle every time the circular-spline makes one complete rotation relative to the
flexspline. To match this moel to experimental observations, a phase-shift of Ob is also
included in this equation. Using the experimental results from sections 2.4, 2.5, and 2.6
of this document, values for the friction parameters bconstant, bl,9b2, Ab, and Ob can be
determined readily.
With the given definition of the harmonic-drive friction, the transmission model is
virtually complete. The only remaining task is to realize that the variables 0 nfs, 0 n-cs,
wn_f, and (oncs are identical to the variables Os 0, cofs and wcs respectively, due to the
kinematic constraints in the model. Given this identity, the friction-torque function can be
calculated based on the rotations specified on the harmonic-drive ports and the remaining
torques in the transmission can be found.
116
Chapter 3: Harmonic-Drive Model 2
Using values for the harmonic-drive friction parameters taken directly from the
experimental observations, several numerical simulations were performed for each
harmonic-drive over a range of current commands. As a sole indicator of model
performance, the resulting input-velocity response for the three different harmonic drives is
illustrated in figures 3.14, 3.15, and 3.16. These three plots can be compared to the
experimental input-velocity response shown in figures D.1.2, D.2.2, and D.3.2 in
Appendix D to evaluate model integrity.
117
Chapter 3: Harmonic-Drive Model 2
25000 T TI
. .. .. .6.0 amps
5.0 amps4
20000-
0
0
4.0 amps
0.
-. 3.2 amps
5000
.. ...... .... .... ..... .... ..... ..... ..S.. .. ........... ...
118amp
Chapter 3: Harm onic.Drive Model 2
4.5 amps
3.8 amps
0
.............. .. ..... ... .... .. ............. .. .. .... ... ...
CD
!D
300003.4 amps
C.2
3.0 amps
.G
..... ... .... .... .... .. . ......... ...
20000
10000.................... .......................
119. ...
..4
...
.. ...
....
Chapter 3: Harmonic-Drive Model 2
50000 T-7-r rr 1
0.60 amps
........
.*.
..........
. . .........
4030000rr.N 02ap
.........
..............
120.... .......
Chaoter 3: Harmonic-Drive Model 2
results. However, at low velocities, these periodic variations are in closer agreement with
experimental response. This unfortunate behavior indicates that the amplitude of the cyclic
friction function is not constant as anticipated, but may indeed vary with such variables as
rotational velocity.
In conclusion, the results of this section illustrate that the harmonic-drive model
with friction performs significantly better than the ideal transmission model. Due to the
incorporation of observed energy dissipation into the model, average steady-state velocities
show fairly close agreement for operating regions where resonant vibration is minimal.
Unfortunately, since the model ignores kinematic position-error variations, theoretical
predictions are devoid of the velocity fluctuations which appear throughout the
experimental response. Additionally, because these velocity fluctuations become intense
during system resonance, theoretical results become very inaccurate over a large portion of
the observed operating range. Based on these shortcomings, it is clear that an accurate
model must capture velocity turbulence and resonance behavior. Since, as discussed
previously, this behavior is spawned from the interaction between kinematic error and
121
Chapter 3: Harmonic-Drive Model 2
122
Chapter 3: Harmonic-Drive Model 3
123
Chapter 3: Harmonic-Drive Model 3
reflected to either of the other two ports of the transmission with no effect on the dynamic
behavior of the system. Due to the location of the friction element with respect to the
compliance in this model, the stiffness profile measured from either output port of the
transmission will display the characteristic hysteresis loss due to the velocity-independent
friction component in the friction-torque function.
Given this simple lumped-parameter model and the notation definitions illustrated in
figure 3.18, the derivation of the equations of motion is straightforward. First, using
force-balance equations in (3.18) and (3.19) as well as the friction function in (3.21), the
torque on the output ports of the harmonic-drive, Tf&and Tcs, can be related to the variables
Tn~f5 and Tn_cs. Second, using the ideal gear-reduction relationship provided in equation
(3.17), these torques can then be equated to the torque across the spring and the torque at
the wave-generator:
Tk = f()
Tb Z=f( I
124
Chapter 3: Harmonic-Drive Model 3
Given this identity, the torques on all three harmonic-drive ports can be calculated once the
torque across the spring is known.
3
Tk = k, (Owg - en_wg) + k 2 (Owg - en-wg) (3.26)
where k, and k2 are the linear and cubic stiffness constants, respectively. To further refine
this stiffness characterization, a more complex equation could be used which includes the
influence of wave-generator orientation and gear-tooth preload on the stiffness profile.
Since, as experimental observations indicate, this influence can cause fluctuations in
stiffness profiles of over thirty percent, a truly accurate stiffness function cannot ignore
these effects. However, since the necessity for this level of detail is unclear, these
influences will be neglected and the simple cubic relationship in (3.26) will be featured
exclusively in the dynamic model.
The derivation of equations for this model can be completed by recognizing that
kinematic constraints dictate 0 n-f$ is equal to Of, and On-,s is equal to 0,,. From these
identities, 0nwg can now be found from the kinematic requirement of the ideal
transmission model in equation (3.15):
Given values for Owg, 0 fs, and Ocs, these equations can be used to calculate the resulting
torques on the three ports of the harmonic-drive transmission.
To evaluate the validity of the stiffness representation used in this model, two
hallmarks can be employed. First, the stiffness function should replicate the experimental
dynamic behavior observed at system resonance. In particular, resonant frequencies and
vibration amplitudes in the model should coincide with observed experimental behavior.
125
Chapter 3: Harmonic-Drive Model 3
Second, when static stiffness data is generated from the transmission model, it should
closely resemble the experimental stiffness profiles. Specifically, if the interaction between
the model friction and compliance is accurate, the hysteresis loss and curve shape in the
actual stiffness data should be captured. In order to observe the influence of resonant
behavior, both of these criteria will be discussed in light of the harmonic-drive model
which includes kinematic transmission error to be discussed later.
126
Chapter 3: Harmonic-Drive Model 4
127
Chapter 3: Harmonic-Drive Model 4
S * 0 0 Elm0
can be reduced with little penalty on model accuracy by combining these two components
into a single function.
128
Chapter 3: Harmonic-Drive Model 4
Second, using the non-linear damping function developed in model 2, the constitutive
relationship for the friction element in the model is
.T,.= f(0)
(•le,,=(e),OZ= _f(0) T
IT,= f (o))
129
Chapter 3: Harmonic-Drive Model 4
where the three independent friction-torque functions are defined in equations (3.22),
(3.23) and (3.24). Third, as discussed in the third model and defined in equation (3.26),
the torque function provided by the stiffness element can be described by
3
Tk = ki (Owg - Ok.out) + k 2 (ewg - Okout) (3.32)
Given these definitions for the gear-reduction, frictional, and stiffness properties of the
transmission, the constitutive constraints in the model will be complete once the torque and
rotation behavior of the kinematic-error element is characterized.
Notice that, in this equation, each term is a sinusoidal function of the wave-generator
rotation, Owg, with a specific amplitude, Ai, and phase, Oi, describing the magnitude and
relative orientation of each error component. Since these error components result from gear
error on the flexspline and circular spline, a strictly accurate representation would have
these functions vary with harmonics of the gear-tooth meshing frequencies. For example,
in joint configuration I where the flexspline is attached to ground, every rotation of the
wave-generator propagates the tooth-engagement zone completely around the flexspline
circumference but fails to reach two teeth on the circular spline. Consequently, the
frequency of the flexspline gear-tooth error components correspond exactly with the
rotational frequency of the wave-generator while the frequency of the circular-spline gear-
error components vary with the difference between the wave-generator and circular spline
130
Chapter 3: Harmonic-Drive Model 4
rotation. However, since the input rotation is much larger that the output rotation in the
transmission, this distinction can be neglected and the aggregate kinematic error can be
varied as a function of the wave-generator rotation alone. In models which desire to
capture the kinematic behavior uf the harmonic drive to a higher level of detail, this
clarification may be important.
Given this representation for the position error in harmonic-drives, the constitutive
relationships for the kinematic-error element in figure 3.20 can be defined. First, using the
error-function, Odn, given in equation (3.33) the position constraint on this element can be
specified as
0 0 (3.34)
n-wg = 0k_out + N erfn
Since the magnitude of the error function is measured in terms of degrees of output
rotation, it must be scaled by the transmission gear-ratio. The velocity constraint on this
element can be defined by taking the derivative of the position constraint:
Using this definition, the velocity constraint in equation (3.35) can now be rewritten as
By applying the law of power conservation, this velocity constraint can be used to derive
the torque constraint across the position-error element. First, the power at the input port of
the element must always equal the power at the output:
131
Chapter 3: Harmonic-DriveModel 4
By substituting the velocity constraint in (3.39) into this equation, the torque constraint
becomes
Since the velocity variations across the spring will be small compared to average velocity
values, cOk_,out can be assumed to be identical to wwg and the torque relationship reduces to:
Therefore, using this relationship and the equation for ,efn in (3.38) the torque constraint
required to ensure power conservation across the kinematic-error element is complete.
Having developed the position, velocity, and torque constraints for the position
error in the model, the constitutive relationships governing model performance are
complete, and the equations of motion can be finished by specifying the remaining
kinematic and torque-balance relationships. Specifically, as seen in equations (3.18) and
(3.19), a torque-balance across the output ports of the harmonic-drive yields the
relationships
Also, by performing a similar force balance across the spring on the input of the
transmission the following relationship results:
Lastly, by specifying the kinematic constraints at the output ports of the harmonic-drive,
namely
On fs = Ofs, (3.48)
the equations of motion are complete. Therefore, using the equations presented in this
section, the positions and velocities at the three harmonic-drive ports can be directly related
to the torques across the harmonic-drive.
132
Chapter 3: Harmonic-Drive Model 4
Estimating the values for the non-linear stiffness element in model 4 proved to be
more complex than capturing friction behavior. In particular, due to the precarious
influence of friction on the static stiffness measurements outlined in section 2.3 of this
document, the resulting stiffness values delivered poor dynamic performance when inserted
into the model. As a more direct method of determining the dynamic stiffness properties of
each transmission, the experimental dynamic-response measurements outlined in section
2.6 were used to identify the resonant frequencies of the three harmonic-drive systems.
From this known frequency and the given inertia values, the linear component of the ideal
stiffness element could be adjusted to match the observed resonance. Given this linear
component, the cubic stiffness coefficient required by equation (3.32) was selected so that
the static stiffness profile generated by the model most closely resembled the experimental
static stiffness behavior. Figures 3.21, 3.22, and 3.23 illustrate the resulting stiffness
profiles used in the three joint-models as compared to experimental data and catalog
predictions. In each of these plots, it can be seen that the theoretical stiffness profiles show
significant deviation for the static stiffness measurements. In particular, the hysteresis loss
in each curve, as dictated by the constant friction component used in the model, Tb_constt,
while producing good dynamic response, is unreliable as a predictor for static hysteresis
loss. Additionally, the shape of the theoretical stiffness profiles required to ensure good
dynamic response in the three harmonic-drive models provides a poor estimate of the
measured static-stiffness profiles. As this comparison illustrates, for all three joints, the
experimental measurements typically yield static stiffness values that are higher than the
actual dynamic stiffness in the drive. This discrepancy is not surprising since Coulomb-
like friction in the transmission can increase the measured torque values during static
stiffness measurement. To address this problem, a new model will be developed later
which predicts these Coulomb-friction effects.
133
Chapter 3: Harmonic-Drive Model 4
0.35 1 1 1 1[ T . i
................ ................
S............... ....... . ............... i.......... . .
I............... ............... "............... .........
~............... .... ...
i...............................
S........ !...............
1
0.3 .-...
............................................................................................... experim ental stiff ness data
.c .......0.5
~j - /- .............. - - i... curve
s0.25 stiffness
standard
S ............................... ................ i.............. ............... ............. ..........
-- optimal
s if e sstiffness
c r e ncurve
mo l
S-
C
.2...............
2................ ./ ........
.-
i. . . . . . . . ...............
. . . . . ..i. . . . . . . . !...
i...
0 .2 . .......... ............... ................ .... ............
. . . . .....
g0.15...
0
0.05
o........
....
........".......
. ...........
i...............
................................
...............
. . . . . . . . ................
.
0 10 20 30 40 50 60 70 80 90 100
Output Torque (N'm)
0.
S0.
0.1 ..-....
0 10 20 30 40 50 60 70
Output Torque (N'm)
134
Chapter 3: Harmonic-Drive Model 4
0.45 o4. ............ ........... .......... .......... ......... ......... .............. ... ..... - ...... ....
. .............. .... .
i ........... ..........
..........
.......... T........... '."
......... ......... ... ......... ....... ......... ......... .........
.0.25 -
... ............
. ness curve in model
0
0
0.5
1 ..... ...... .....
_ _ _ _
10
_ _
exenena
__ _ _ _
stfns
_
S
_ _
data5
15
_
Given the friction and stiffness properties, the only remaining parameters that need
to be determined are the magnitudes of the kinematic-error components in equation (3.33).
The amplitudes of the three error components were determined by matching the amplitude
of the dynamic position error predicted by the model to experimental values. In theory, if
kinematic error on the testing apparatus could be measured under truly static conditions,
direct estimates of these error components could be made. However, since data on the
apparatus can only be collected for non-zero operating velocities, dynamic vibrations which
induce transmission compliance contaminate the kinematic-error data. As presented in
section 2.2 of this document, by operating the harmonic-drives at slow velocities where
resonance vibration was minimal, measurements were made of the kinematic error in the
three harmonic drives. In table 3.4, these experimental estimates are compared to the error
values which provided optimal dynamic performance in the theoretical model. By
comparing these results, it can be seen that, in most cases, values near the experimental
measurements provided accurate dynamic response simulations, but, in a few instances,
differences between experimental and theoretical values could deviate by up to a factor of
two. From this discrepancy, two possible conclusions can be drawn: (1) slow operating
velocities can significantly influence the magnitude of the position-error measured in the
experimental apparatus, and (2) the model developed in this section does not fully capture
135
Chanter 3: Harmonic-Drive Model 4
Tabl Kieai-ro
3.4: copnnsinmdl4v .
mesre0ale
the influence of kinematic error on dynamic vibrations in the transmission. Since additional
dynamic simulations with the theoretical model illustrated that dynamic vibrations could
cause significant amplification of transmission position error at small operating velocities, it
is likely that the first of these two conclusions is most responsible for the difference
between experimental and theoretical parameters.
WI Position, velocity and torque response in the model shows good agreement to
actual data at operating speeds removed from regions of resonant behavior.
136
Chapter 3: Harmonic-Drive Model 4
Tal
3.5 Sumar of siulte tiersos ltSnapni
n 4o ...
........
t. WrJit 3
E.1.1 E.2.1 E.3.1
::. :. : ::•:
otlo. . ::
..... :..
...
.......
.: : : •. :•
E.1.4 E.2.4 E.3.4
tyve
fuptVlc E.1.5 E.2.5 E.3.5
VeilocitvIe
-O-tput 2i E.1.6 E.2.6 E.3.6
Conversely, by noting the differences between predicted and actual behavior, the
limitations of this model can also be observed:
of output
SAmplitude
torque fluctuation near resonance, especially on joints 1
and 3, is somewhat lower than the experimental torque-fluctuation amplitude.
These limitations on the accuracy of theoretical predictions lend insight into possible areas
for model improvement. First, since the predicted velocities near system resonance are
much higher that the observed velocities, it is clear that this model fails to capture the
enhanced frictional dissipation which appears due to resonant vibration. As discussed in
light of experimental results, this dissipation probably results from a Coulomb-like
frictional torque which increases with the amplitude of torque fluctuations in the
137
Chapter 3: Harmonic-Drive Model 4
transmission. Since, in the three drives tested, the resonance phenomenon dominates
transmission behavior in a wide portion of the observed operating rangei an accurate
harmonic-drive model cannot ignore these resonance losses. As a second consequence of
this resonance-dissipation absence, the velocity jump that occurs in experimental response
when resonance vibrations cease does not appear in the theoretical response.
Consequently, any improved representation which is*able to replicate friction loss at
resonance is likely to show this unusual velocity behavior as well. Furthermore, since a
friction model at resonance is likely to supplement dynamic torque fluctuations with
variations in frictional torque, the low torque-amplitudes observed on the joints 1 and 3 in
this model may also improve. As a final note, the experimental current response on all
three joints displays electrical noise which is not reproduced in the current response of the
model. Since this noise had no effect on the dynamic performance of the actual system, it
was completely ignored in the simulation.
Another conclusion which can be made from this model is that reliable
measurements of transmission stiffness are difficult to obtain. In particular, when stiffness
138
Chapter 3: Harmonic-Drive Model 4
curves collected from static measurements were used in the model, the resulting dynamic
response held few similarities to actual behavior. As explained above, this discrepancy can
be blamed on the inaccuracy of static stiffness measurements due to friction in the
harmonic-drive. From the substantial differences between the theoretical and actual static
stiffness profiles discussed previously, it can be concluded that interaction between
stiffness and friction in the actual system is far more complex than the model recognizes.
Based on this result, an improved harmonic-drive model, in addition to improving
predictions at dynamic resonance, should also better represent the static stiffness behavior
in the transmission. In order to improve model performance, dynamic response data had to
be used to isolate the natural frequency of the system and derive the corresponding
transmission flexibility. Although this new stiffness data provided good dynamic results,
slight changes in its value could cause significant shifts in the operating velocities which
excited resonant behavior. This unfortunate result demonstrates that, unless an extremely
accurate model is developed which can reliably predict transmission stiffness from static
measurements, it is likely that static stiffness data will never be useful for accurate
predictions of harmonic-drive dynamic behavior.
Perhaps the most important conclusion that can be drawn from the results in this
section is that accurate harmonic-drive modeling is a very complex task. The scope of the
model presented in this section represents the maximum useful complexity that can be
included in a model that ignores more complex geometric interaction of the harmonic-drive
components. Since the resulting model performance was completely unreliable for large
portions of the harmonic-drive operating ranges, the forecast is gloomy for designers who
wish to make accurate predictions using a simple transmission representation. In a final
attempt to create a model that can make accurate predictions while still containing a
manageable amount of complexity, a representation will be developed below which
incorporates the three transmission properties used in model 4 into a simplified
representation of harmonic-drive gear-tooth interaction.
139
Chavter 3: Harmonic-Drive Model 5
140
Chanter 3: Harmonic-Drive Model 5
which reduce the complexity required in a harmonic-drive model which mimics this
behavior. First, by realizing that during normal harmonic-drive operation, anywhere from
141
Chapter 3: Harmonic-Drive Model 5
ten to greater than fifty teeth can be in contact between the circular spline and flexspline, the
effect of individual gear teeth can be ignored in a dynamic model. This assumption was
confirmed by the kinematic error results in section 2.2.2.1.4 which discovered that the
error of individual teeth was not detectable. Consequently, in order to capture the
aggregate effects of gear-tooth meshing, a model which describes sliding on a single
continuous gear-tooth surface should suffice. Second, since the gear-teeth are very small
compared to the pitch-diameter of the transmission, the localized motion of the gear-teeth
during meshing shows little trace of the overall rotational motion. Based on this
assumption, a model which describes this gear-tooth action can be simplified by using
translational rather than rotational motion. Using these two assumptions, a dynamic model
that captures the aggregate effect of the gear-tooth meshing mechanism is illustrated in
figure 3.25.
wave-generator,
rotation
FiSr
3.5 of haionSdrv Scenai moe Sihga-oohgoe
142
Chapter 3. !Iarmonic-Drive Model 5
actual transmission. These planes are confined to horizontal and vertical motion by a series
of ideal rollers and rigid surfaces. In particular, note that the output ports of the wave-
generator, flexspline, and circular spline are all restricted solely to horizontal motion so
they can be interfaced easily with an overall system model. Although this cannot be
illustrated in the schematic, by assuming that all of the surfaces which constrain motion are
infinitely long, this model can be used to describe all ranges of transmission operation.
By comparing the actual transmission in figure 3.24 to the model in 3.25, several
parallels can be drawn between actual and modeled behavior. First, horizontal motion in
the model corresponds directly to rotation of the harmonic-drive components and vertical
motion in the model corresponds to radial movement in the harmonic drive. Second, the
low-angle plane and inclined bearing surface on the bottom of the model schematic
correspond to the wave-generator in the actual drive. The mechanical advantage and
kinematic properties of the wave-generator eccentricity are captured by the low-angle of this
wedge. Since the wave-generator motion pushes the flexspline gear teeth into the circular
spline, a third parallel can be drawn between the two central wedges and the flexspline in
the actual system. Similarly, the wedge at the top of the model schematic which mates with
the flexspline wedge corresponds directly to the circular-spline in the actual system. So
that this model effectively replicates the gear-tooth rubbing action in the actual
transmission, the angle of both the flexspline and circular spline mating wedges is identical
to the traditional harmonic-drive gear-tooth angle of thirty degrees. The angle of the wave-
generator wedge at the bottom of the schematic can be calculated from the gear-ratio of the
particular harmonic drive as discussed later.
A better understanding of the likeness between the model in figure 3.25 and the
actual transmission can be gained by considering the case where the flexspline is fixed to
ground. Under this condition, as discussed previously, clockwise rotation of the wave-
generator in an actual transmission produces clockwise rotation of the circular spline. This
relative motion results since the wave-generator pushes the flexspline gear teeth .Ito the
circular spline gear teeth, and the forces at the gear-tooth interface drive the circular spline
in the clockwise direction. Identical behavior can be seen in the schematic in figure 3.25 if
the flexspline port is fixed to ground. Specifically, the clockwise rotation of the wave-
generator can be mimicked by sliding the wave-generator wedge in the model horizontally
to the right. Due to this motion, the two flexspline wedges in the model, which are not
allowed to move horizontally, slide vertically upward and push against the circular-spline
wedge. This vertical motion directly parallels the engagement of the gear teeth that occurs
in the actual transmission. Continuing the analogy, the vertical motion of the flexspline
143
Chapter 3: Harmonic-Drive Model 5
wedges in the model causes the circular-spline wedge to move horizontally to the right.
Again, this horizontal motion mimics the clockwise rotation that results in the circular
spline of the actual transmission. By fixing the circular-spline port in this model to ground,
a similar exercise can be performed to relate model behavior to the meshing action in an
actual harmonic drive under a different operating configuration.
Given this geometric framework which describes the motion and force constraints
present in ha:monic drives, remaining dynamic effects can be incorporated by introducing
friction, compliance, and kinematic-error properties. Specifically, since it is likely that
most frictional loss in harmonic-drives is generated from rubbing at the gear-tooth
interface, all friction in this model is localized at the contact surface between the flexspline
and circular spline wedges. However, unlike the previous model, Coulomb friction
properties are also included at this contact surface to provide a mechanism for
demonstrating increased losses at system resonance. Compliant behavior, on the other
hand, is not as easy to pinpoint as friction. In particular, past research has indicated that, in
addition to deflection in both the flexspline and circular spline gear-teeth, significant
deflection may also occur radially in the wave-generator as well. Since no information was
available to supervise the distribution of the overall transmission compliance into these
individual locations, a single stiffness element was included at the flexspline for simplicity.
In order to maximize the influence of friction at the gear-tooth interface on static stiffness
measurements in the model, this spring was oriented vertically rather than horizontally.
Lastly since kinematic error in harmonic drives is caused mainly by gear-tooth
manufacturing errors on both the circular spline and flexspline, two separate error function
should be included which vary the positions of the flexspline and circular-spline gear-tooth
surfaces. However, since the combined effect of these two error elements can be closely
approximated by a single error function, only one kinematic-error element is included in the
model. For the sake of simplicity, this element is located at the flexspline. Also, since
numerical problems can arise if the velocity at the input of this element becomes zero, the
kinematic error is oriented vertically rather than horizontally . If future models require
improved accuracy, this error function can be distributed into different locations in the
model.
144
Chapter 3: Harmonic-Drive Model 5
better understanding how the load is carried by the harmonic-drive gear-teeth during these
torque oscillations, a new model can be developed to capture this behavior. In particular,
as figure 3.26 illustrates, when the driving torque across the gear-teeth is positive, the load
is carried by the gear-teeth surfaces which push the circular spline in the positive direction.
. . *. ... .. .*
.. ....
........
. ....
•.........
...:g" : ............
145
Chapter 3: Harmonic-Drive Model 5
However, when the driving torque becomes negative due to inertial vibrations, the
transmission velocity is still in the same direction but the opposite face of the gear teeth
must bear the resulting load. In order to capture this behavior in a dynamic model, the
representation shown in figure 3.27 can be used whenever the normal force at the gear
tooth interface becomes negative. Since this new model is essentially a mirror image of the
previous one, the normal force at the gear-tooth interface becomes positive and resulting
coulomb friction can be easily determined. Using this representation along with the
previous one, an attempt was made to mimic both positive and negative gear-tooth loading
in the harmonic-drive model by rapidly switching between models during simulation.
Unfortunately, since the numerical simulation used to solve the equations of motion was
not very tolerant of discontinuous changes in the dynamic qvuations, a solution could not
be reached. Instead, the single representation in figure 3.25 was used exclusively, and
whenever the normal force at the gear-tooth surface became negative, the coulomb friction
146
Chapter 3: Harmonic-Drive Model 5
was set to zero. This inaccurate assumption most likely represents the main shortcoming of
this model and identifies an important area for further improvement.
Given the variable names provided in these two figures, the geometric constraints in
the model will be used to define the kinematic equations governing model motion. First,
due to the inclined face of the wave-generator wedge, the horizontal motion of the wave-
generator and flexspline can be directly related to the vertical motion of the flexspline base
by the equations
0
toothbase - tan awg (Gwg - Of), and (3.49)
Similarly, due to the gear-tooth angle, the circular-spline and flexspline horizontal motion
can be used to determine the vertical motion of the flexspline tip:
147
Chaoter 3: Harmonic-Drive Model S
8
toothaxfmlc bhawac
(X oot
V148
Chapter 3: Harmonic-Drive Model S
Tbh normal
149
Chapter 3: Harmonic-Drive Model 5
Using these equations and the geometry of the model, the equations specifying the motion
at both the wave-generator surface and the gear-tooth surface can be derived:
( 0wg_surface COS aOwg C0 wg + sin awg (0oooth_base - COS xwg (, fs, (3.54)
Otooth_surface = COS
0 tooth Otooth-tip + sin atooth Ocs - (3.55)
sin aooth Ofs , and
0
W toothsurface = COS atooth Qtoothjtip + sin aO•toh Wcs - (3.56)
sin atooth Cofs.
Now that the kinematic relationships are complete, the angle of the wave-generator wedge
required to reproduce the motion in the actual transmission can be found. Specifically,
given the gear-tooth angle, atoth, and the catalog gear ratio, N, the wave-generator angle
can be determined exactly from the equation:
If these kinematic constraints are indeed accurate, they should produce the same
relationship between the positions and velocities of the wave-generator, flexspline, and
circular spline as observed in the actual system. As a confirmation of this fact, while
ignoring kinematic error and compliance in the model, equations (3.49), (3.50), (3.51),
(3.52), and (3.57) can be combined to reproduce the ideal three-port transmission
equations presented in section 3.3. Based on this result, the position and velocity
relationships presented above provide a complete and accurate description of the kinematic
behavior of harmonic drives.
Now that the kinematics of the harmonic-drive model with gear-tooth geometry are
better understood, the constitutive relationships which govern the behavior of the individual
model elements will be derived. First, as seen in previous models, the friction at the gear-
tooth interface is composed of (1) a velocity-independent, constant torque, (2) a velocity-
dependent, damping torque, and (3) a cyclic torque function. However, as a new and
important improvement to this model, a Coulomb friction function is also included at the
gear tooth surface. The resulting aggregate friction function is
150
Chapter 3: Harmonic-DriveModel 5
where
In these equations, g. represents the Coulomb friction coefficient at the gear-tooth interface
and the remaining parameters, as seen before, are defined in table 3.3.
Unlike this friction relationship which contains a new term for the Coulomb losses,
the stiffness equation used in this model resembles closely the cubic relationship seen in
previous models:
3
Tk = k, (0err - Otooth-tip) + k 2 0err - 0tooth-tip) (3.63)
Following the example set by model 4, the kinematic error function used in model 5
will be a harmonic series of three error components. However, slightly different than
model 4, since the kinematic-error function is located on the flexspline, it should vary as a
function of the flexspline gear-tooth meshing position rather than the wave-generator
position. Considering that every tooth on the flexspline will mesh once when the wave-
generator makes one complete revolution relative to the flexspline, the position which
governs the kinematic-error function should be ewg - efs. Adopting the syntax used for
model 4 in equation (3.33), this kinematic error function becomes
where
151
Chapter 3: Harmonic-Drive Model 5
Note that, unlike previous models, since the magnitude of the kinematic-error function in
model 5 does not translate directly to rotation of the flexspline or circular spline, the values
of the position-error amplitudes, Al, A2 , and A3 , must be adjusted as discussed later.
Given the position-error function and its derivative, the position, velocity and torque
constraints can be derived as done for model 4. Specifically, by definition, the position
and velocity constraints across the kinematic-error element are
By requiring that power be conserved, the following equation can also be derived:
Using the kinematic-error relationship in equation (3.70) and the geometric constraint in
equation (3.50), this power relationship can be convened easily into the desired torque
constraint on the kinematic error element:
Terrin
- ' awgI Tkk
I + tan erfn (3.72)
At last, with the torque, velocity, and position relationships for the kinematic error fully
defined, the constitutive relationships required in the model are complete.
Having calculated both the kinematic and constitutive constraints on the model, the
task that remains is to balance the forces across the elements in the transmission. This can
be accomplished by summing the forces separately for the four different geometric shapes
in the model. For the wedge representing the wave-generator, balancing the forces in the
horizontal and vertical directions yields the equations
152
Chapter 3: Harmonic-Drive Model 5
Tfs-normal= Ttoothnormal COS Octooth- Ttooth friction sin actooa, and (3.77)
Using the same procedure, the force-balance equations for the remaining circular-spline
wedge complete the dynamic constraints:
To solve these equations, (3.74), (3.76), (3.78), and (3.80) can be combined with the
constitutive torque-relationship for the kinematic error element in order to determine
TinpuL normal, Twg.normal, Toutput-normal, and Ttooth_normal in terms of known quantities.
Then, by substituting these values along with equation (3.77) into equations (3.73),
(3.75), and (3.79), the expressions for the torques on the wave-generator, flexspline, and
circular spline can be found. When the positions and velocities of the three harmonic-drive
ports are known, these equations, combined with the kinematic and constitutive
relationships presented above, can be used to calculate all torques experienced by
transmission.
153
Chapter 3: Harmonic-Drive Model 5
gear-tooth surface velocity can be easily found from the kinematic relationships presented
above. On joint configuration 1, for example, this scaling factor is
B= 1 , where (3.81)
(N + 1) sin atooth
Given this scaling factor, the values for the tooth-surface friction coefficients required to
generate the same losses observed experimentally can be calculated readily. By assuming
that any power dissipated by a friction source located at the wave-generator should be
identical to the power dissipated when the friction is located at the gear-tooth interface, the
experimental friction results can be directly related to friction parameters at the gear-tooth
interface. In the case of velocity independent friction, this analysis yields the result
bconstan- =bconstantwg
B (3.82)
where bomumt is the constant friction at the gear-tooth interface from equation (3.59),
bcovmta wg is the constant friction measured at the gear-tooth interface, and B is the scaling
factor presented in equation (3.81). Since the amplitude of the cyclic friction in equation
(3.61) is also independent of velocity, it can be calculated in an identical manner.
However, the dynamic friction coefficients in equation (3.60), since they are linear and
cubic functions of velocity, must be scaled by factors of B2 and B4, respectively.
Since no experimental data were collected to measure the Coulomb friction in the
transmission, the remaining friction equation, (3.62), cannot be calculated using the same
method used to derive the other friction terms. Instead, tabulated handbook values for
typical Coulomb friction coefficients were used. Specifically, for lubricated steel-to-steel
contact, listed values ranged from 0.03 to 0.16, and using this range as a guide, friction
coefficients were then adjusted to optimize model performance. The resulting coefficients
that delivered good model performance were 0.15 for joint 1, 0.05 for joint 2 and 0.10 for
joint 3. However, simulated results illustrated that small changes in these friction values
within the quoted range could greatly influence the dynamic response. This unfortunate
result identifies a serious barrier to accurate harmonic-drive modeling since reliable
Coulomb friction values are notoriously difficult to identify.
Having determined useful values for the friction in the harmonic-drive model, the
second property that needs to be evaluated is transmission flexibility. As was the case in
154
Chanter 3: Harmonic-Drive Model 5
the previous model, static stiffness measurement were of little help in determining the ideal
compliance for the model. Instead, the resonances that appeared in the experimental time-
response data were used to identify the natural frequencies of each harmonic-drive system.
Given these frequencies, the inertias seen by the spring, as reflected through the geometry
of the model, were used to calculate the linear component of the ideal compliance required
to produce the same resonant behavior. With the linear coefficient fixed, the cubic stiffness
coefficient was adjusted to reproduce the curvature in the measured static-stiffness profiles
and provide good dynamic-response behavior around resonance. In particular, cubic-
stiffness values were adjusted to increase or decrease the velocity range over which
vibrations were sustained in each system resonance. Using these stiffness parameters,
static stiffness simulations were performed using the harmonic-drive model, and results are
shown in figures 3.30, 3.31, and 3.32.
0.35 1 1 - I .1 .... I
S............... T............... !............... !................ :....... ...
• ............ ...
T ........ - '"- ". .. ....
! .............. .......... ..I
........... ............... ..............................
................. . .. ....................
0
0 10 20 30 40 50 60 70 80 90 100
Output Torque (N*m)
155
Chapter 3: Harmonic-Drive Model 5
...........
...... ...... . . . ............
............. ..........
...............
.......
..........
.
0
0 10 20 30 40 50 60 70
Output Torque (N'm)
UA)
CL
0
0.35 1 1
S....
........... ............
........... ......'"•
.... ....... ................
..... .... ......... .........
........
E"0.32
01 . . . .....
...........
5.......... . . . . . ..... ................. :.."........i
S 0 .2 . ..... *.......... ...... .......... ......... . ..... .......... .... ...... ............. ..................... i....
--. ......... .
0 ....... ........... ,. .......... ..'... .........-
....... .... ..............
.......... . . . ......
....... .... ..... ........-..
... ............
.. ".. ...... .... .... ....
o. .......
. ... ....- . ... ....
o. 1 ....... ist
n tar
stiffness
l dur
te
~
0.1 ~~ ~ ~
............ .......... optiandr stiffness curve
0.0 ...... ........ ..... .......... ......... i•......... .......... g arto th mo e
...... .. |.......... ........... <.......... ..........
. .... stf
lllll i in i lll l m | I
I.. .aI)lll
0 5 10 15
Output Torque (N~m)
156
Chanter 3: Harmonic-Drive Model 5
friction must be overcome. Since the magnitude of the Coulomb friction force at the gear-
tooth surface increases as the load on the surface increases, the applied load required to
overcome this friction force increases linearly as the load increases. During static stiffness
measurements, this variation in Coulomb-friction force will appear as an increase in the
measured stiffness during joint loading and a decrease in stiffness when the joint is
unloaded. This behavior can be observed in all of the theoretical stiffness profiles in
figures 3.30, 3.31, and 3.32. As illustrated most clearly on joints I and 3, this apparent
increase in stiffness during loading produces theoretical stiffness measurements that more
closely resemble the experimental stiffness profiles. This result is positive evidence that the
gear-tooth mechanism captured in the model more closely mimics the physical behavior in
the actual system. However, since in all cases, the overall shapes of the predicted stiffness
profiles show poor agreement to measured values, the actual friction mechanism in
harmonic drives is probably much more complicated than this model recognizes.
The second hallmark which can be used to evaluate the accuracy of the theoretical
stiffness predictions in hysteresis loss. In the model, hysteresis in the calculated stiffness
profiles is governed by the velocity-independent friction component as well as the Coul )mb
friction at the gear-tooth interface. In particular, the constant friction term provides a fixed-
width difference between the loading and unloading stiffness curves while the Coulomb
friction adds a hysteresis-width component that increases linearly with the applied load.
The combined effect of these two friction factors can be seen in the theoretical stiffness
curves in figures 3.30, 3.31, and 3.32. From the poor agreement between the predicted
and experimental hysteresis loss it can be concluded that, although the specified friction
values provided good dynamic response, the interaction between friction and compliance in
the actual system is more complex that the model acknowledges.
Now that values for the friction and compliance effects in the transmission model
have been isolated, the only parameters that remain unknown are the kinematic-error
amplitudes of equation (3.64). As discussed previously, since low-velocity measurements
of these position-error components proved somewhat unreliable due to dynamic vibrations,
the values were selected in the previous model by matching predicted and experimental
dynamic position-error. In theory, since the resulting values used in model 4 should
describe the absolute kinematic error in the transmission, the same numbers should also
apply to the new harmonic-drive model with gear-tooth geometry. When these kinematic-
error amplitudes were reflected appropriately through the geometry of model 5 and
incorporated into the simulation, however, the resulting dynamic-response predictions fell
short of performance expectations. Specifically, the torque-fluctuation amplitudes that
157
Chapter 3: Harmonic-Drive Model 5
resulted from the given kinematic-error values were much smaller than the experimental
values. Because of these diminished torques, the amount of dissipation at resonance due to
Coulomb friction was lower than observed in the actual system. To remedy this problem,
it was discovered that, with slight adjustments to the calculated kinematic-error values,
model predictions improved greatly. The final values which delivered optimal model
performance are listed in Table 3.6 and compared to the experimental predictions from
section 2.2. Based on this evolution of reliable kinematic-error values, a few useful
conclusions can be derived. First, since the amplitudes which provided good results in
model 4 performed less admirably in model 5, it is clear that the mechanisms used to
introduce kinematic-error effects in both these models are incomplete. For example, as
discussed previously, if the model including gear-tooth geometry could also incorporate the
transmission behavior when the normal force at the gear-tooth interface reverses direction,
it is likely that the influence of kinematic error on dynamic behavior could change
dramatically. Additionally, from the difference between low-velocity measurements and
theoretical values shown in table 3.6, the reliability of kinematic-error measurements which
are not taken under truly static conditions is further questioned. In conclusion, the
unfortunate result which can be derived from these observations is that the accurate
kinematic error estimates are elusive. However, through careful parameter adjustment
based on a good understanding of system dynamic operation, values can be chosen that
deliver relatively good model performance.
.... .. ...........
By inserting the harmonic-drive model developed in this section into the overall
system models for the two joint configurations, dynamic step-response simulations were
158
Chapter 3: Harmonic-DriveModel 5
performed to replicate the experimental results. Using the experimental parameters for
friction, stiffness, and kinematic error discussed previously, time-response data was
collected over a range of motor-current steps for the three harmonic-drive testing-station
models. These simulated results are presented in appendix F in the identical format used
for the experimental time-response results of appendix D and the model 4 simulation results
in appendix E. In particular, graphs are shown which illustrate the position, velocity, and
torque response for the observed range of harmonic-drive operating behavior. Table 3.7 is
extracted from a similar index in appendix F and summarizes the figure numbers for all of
the simulated dynamic-response plots. By directly comparing these plots to the
experimental results in appendix D, the capacity of the dynamic model to predict actual
harmonic-drive behavior can be evaluated. As an explanatory note, the electrical noise
present in current-sensor readings was not reproduced in the model simulation since it had
no influence on the dynamic behavior of the system.
i ii
......
.... --
----------.
.F.1.1 F.2.1 F.3.1
`-`.
':::::x11_ý
.... .................
OF.1.2 F.2.2 F.3.2
[F] The surprising jumps in velocity that appear in the experimental response
when the velocity exceeds the resonance range are also observed in the
simulated response.
159
Chapter 3: Harmonic-Drive Model 5
By noting the difference between the experimental and theoretical results, the shortcomings
of this model can be recognized as well:
F21 -Thepredicted dynamic position error at system resonance was greater than
experimental observations for all three joints.
From the improvements in simulated response shown by this model over its predecessors,
it is clear that gear-tooth geometry and Coulomb friction are essential components for
predicting harmonic-drive behavior during resonance. However, from the discrepancies
between the theoretical and actual response, it also clear that a more detailed representation
is required before reliable predictions of harmonic-drive performance can be made. In
particular, as discussed previously, a model that can accommodate the shift in system
dynamics that occurs when the torque at the gear-teeth changes direction should be able to
improve predicted results. Currently, the model in this section avoids this complication by
setting the Coulomb friction at the gear-tooth interface to zero whenever the normal force
on the gear-tooth surface becomes negative. It is likely that a new model, which captures
more accurately the actual friction behavior at the gear-teeth, will improve predicted
performance by (1) increasing loss at system resonance and (2) limiting dynamic position
error and velocity fluctuations while preserving torque amplitudes. If the new model can
indeed deliver these effects, the model predictions should resemble more closely the actual
dynamic response and, consequently, improve the accuracy and reliability of the
representation.
160
Chapter 3: Harmonic-Drive Model 5
caused substantial deviations in the resulting dynamic response. Since experimental results
illustrate that resonant behavior governs system performance for a wide range of operating
velocities, it is not surprising that small variations in model stiffness can have a profound
effect on dynamic response. Second, variations in the coefficient of Coulomb friction, g.,
were also shown to influence dynamic response especially at system resonance.
Surprisingly, decreasing the value of this coefficient often resulted in increased resonance
losses due to heightened system vibration and torque fluctuations. Lastly, the output
friction included in the system model had a surprisingly important relationship on dynamic
behavior. Because increases in this variable resulted in higher normal-forces at the gear-
tooth interface, resonance losses from Coulomb friction were very dependent on its value.
A model in which Coulomb friction acts for both positive and negative normal torques may
not show this same dependence. These observations emphasize the necessity of accurate
characterization of harmonic-drive stiffness and friction properties for reliable model
performance.
161
Chapter3: Harmonic-Drive Model 5
The model presented in this section represents a first attempt to incorporate gear-
tooth-rubbing action into a dynamic model of a harmonic drive. As the simulated results
illustrated, since this model included Coulomb friction losses at the gear-tooth interface,
energy dissipation increased during system resonance and dynamic response integrity
improved greatly in these areas. However, even with this new level of detail in the model,
discrepancies between predicted and actual response indicate that the dynamics of
harmonic-drive operation are still even more complex. However, since this is the only
model which can boast even partially describing harmonic-drive response in all ranges of
operation, it may be useful in any application in which understanding harmonic-drive
operation through resonant regions is required. Since the mathematical complexity of this
model is relatively minimal, it is an ideal candidate for making quick estimates of operating
behavior, but, if more accurate information is required, I hope that this model can provide a
template for more complex representations which better capture the subtleties of harmonic-
drive behavior. In the very least, the model presented in this section can prescribe a
definite lower bound on the level of detail required in harmonic-drive models to deliver
reliable results.
One of the most discouraging conclusions that can be drawn from the results in this
section is that accurate characterization of transmission properties can be very difficult but,
unfortunately, also essential for good model performance. In fact, the only way in which
accurate values could be determined for the stiffness and coulomb friction in this model
was through careful measurement of dynamic response over a wide range of operation. By
comparing these dynamic response measurements to the simulated results, the stiffness and
friction values could be painstakingly adjusted to optimize model performance. Static
stiffness tests and catalog predictions provided little guidance in this task. Based on this
result, it seems unlikely that the model in this section, or any similar representation for that
matter, will ever be able to make accurate predictions of dynamic performance using solely
catalog values or simple experimental observations. Consequently, for applications in
which an accurate harmonic-drive model is required,there are no shorcuts around a detailed
experimental analysis of the actual harmonic-drive system.
Perhaps the most useful feature of the model in this section is its intimate conceptual
relationship to the actual behavior in harmonic drives. Since this model depicts the
interaction of the three different transmission components, harmonic-drive properties such
as friction, compliance, and kinematic error can be placed at locations which have a definite
162
Chanter 3: Harmonic-Drive Model 5
163
164
As presented in this document, my research was targeted at understanding the
operating behavior of harmonic-drive transmissions and developing mathematical
representations to describe the observed performance. From my experimental
investigation, several properties of harmonic drives where characterized and related to the
overall dynamic behavior of the transmission. By developing models to simulate the
experimental observations, the importance of many harmonic-drive properties was
categorized and the complexity of the harmonic-drive operating mechanism was exposed.
This section summarizes the important insights gained from both the experimental and
theoretical investigations and presents recommendations for future harmonic-drive research
and implementation.
165
Chapter 4: Conclusions
Static and dynamic friction in harmonic drives were the third suspects to be
experimentally audited. In particular, starting-torque measurements were made on three
harmonic drives which confirmed past results demonstrating that factors such as wave-
generator angle and gear-tooth preload can greatly affect static friction. Next, experimental
tests performed at non-zero velocities illustrated that the dynamic friction in harmonic
drives has three components: (1) a velocity-independent torque probably due to coulomb-
like friction resulting from gear-tooth preload, (2) a non-linear dynamic torque caused by
the lubrication effects in the transmission, and (3) an increased energy loss during system
resonance resulting from high torque fluctuations which increase coulomb-like friction at
the gear-teeth. The increased losses at resonance, in particular, were sighted for causing
unusual and unpredictable behavior during harmonic-drive operation. Lastly, to add an
extra element of complexity, the velocity-independent friction torque was observed to vary
considerably with the output rotation of the transmission.
Using the results collected for harmonic-drive kinematic error, compliance, and
friction properties, the overall dynamic response of the transmission was scrutinized.
From this investigation, several behavior patterns common to the three harmonic drives
tested were identified and related to known transmission properties. First, all dynamic
response showed significant velocity and torque fluctuations due to the kinematic
inaccuracies of the transmission. Second, due to the relatively low stiffness of harmonic
drives, these fluctuations excited dynamic resonances over a substantial range of the
harmonic-drive operating envelope. Torque and velocity oscillations were observed to
increase dramatically in these resonant regions. Third, since coulomb-like friction
increased energy dissipation at resonance, large increases in motor current produced little or
no increase in operating velocity around resonance. Additionally, due to this energy loss at
166
Chapter 4: Conclusions
resonance, dramatic jumps in velocity were observed when the system finally exceeded
these resonant regimes. Fourth, since the harmonic-drive stiffness is non-linear, the
natural frequency of vibration increased with increased torque levels at resonance. This
rise in resonant frequency increased the velocity range which excited system resonance and
accentuated the velocity jump that occurred when these resonance ranges were exceeded.
Lastly, due to friction variation with output rotation, the location of the velocity jump at
resonance exit was observed to vary unpredictably.
167
Chapter 4: Conclusions
can influence the tightness of gear-tooth meshing, stiffness and friction are vulnerable to
these influences as well. In general, using this argument, any variation in an environmental
or operating conditions which can affect loading at tbh- gear-tooth interface will be
manifested in the stiffness and friction of the transmission. I am optimistic that future work
in this area can not only quantify this interdependence of transmission properties but can
also develop mathematical representations to predict experimental measurements under
static and dynamic conditions.
168
Cha•ier 4: Conclusions
model was incapable of describing the losses which occur at system resonance, model
predictions at velocities which excited resonance vibration were still inaccurate. Based on
these results, it was concluded that this model provides a slightly more accurate
representation than it predecessor for describing the dynamic behavior of systems which
experience little resonance vibration.
From the results of the five transmission models, the complexity of the task of
harmonic-drive modeling has been exposed. Not only are friction, compliance, and
kinematic error essential factors for describing harmonic-drive behavior, but the gear-tooth
meshing dynamics also play an essential role in determining the overall performance of the
transmission. Additionally, experimental results indicate that the actual values of these
transmission properties can vary significantly with environmental and operating conditions,
which makes accurate characterization of these effects difficult if not impractical. To make
matters worse, not only are reliable measurements of transmission properties difficult to
obtain, but, for stiffness and coulomb friction in particular, precise estimates are critical for
good model performance. In fact, the only way in which accurate values could be
determined for the parameters in this model was by measuring carefully the dynamic
response over a wide range of operation and adjusting of model parameters to match this
observed response. This unfortunate result reveals that reliable harmonic drive models
cannot be developed unless experimental data from the actual harmonic-drive system are
available. In the face of these modeling obstacles, it seems unlikely that any model which
does not capture the physics of harmonic-drive operation to a much greater level of detail
169
Chapter 4: Conclusions
will ever be able to deliver accurate and reliable performance over all operating regions.
Given this conclusion, I hope the models developed in this investigation .will act as a
springboard from which to launch more detailed future explorations.
4.3 Recommendations
Given the conclusions of my experimental and theoretical investigation, I can offer
a few recommendations to designers wishing to model systems employing harmonic
drives. The first step in understanding the extent to which harmonic-drive dynamics may
influence system performance is it to determine the prominence of resonance vibration
within the operating range of the transmission. Due to the unreliability of static
measurements and catalog predictions of stiffness values, resonance vibration can be
accurately located only during actual system operation. However, an order-of-magnitude
estimate for the natural frequency of the system can be made by using linearized catalog
stiffness values and system inertia measurements. If a harmonic of the rotational velocity
of the wave-generator relative to the flexspline coincides with the estimated natural
frequency of the system at any point in the operating envelope of the transmission, it is
likely that resonant vibration will occur. Since resonant vibration greatly influences
operating behavior and is difficult to model, the best solution to this dilemma is to revise
the machine specifications to avoid these ranges of operation.
170
Chatier 4: Conclusions
Lastly, for future harmonic-drive modelers, I propose a few areas for investigation
that seem most fruitful. In particular, as described in detail previously, a harmonic-drive
representation that can describe the transition that occurs when the torque at the gear teeth
changes direction should show improved results over the models in this document. After
this improvement has been implemented, a model which predicts the influence of gear-tooth
preload on transmission friction and compliance is also likely to improve dynamic
response. Finally, since the kinematic error in harmonic-drives probably influences the
gear-tooth preload itself, a representation which incorporates this behavior will more
accurately represent the actual physical mechanism inside the harmonic-drive transmission.
171
172
This section contains citations for literature referenced in this thesis as well as
citations for related harmonic-drive research. For convenience, the list of references on
harmonic-drive research has been subdivided into categories describing the nature of the
literature and popular areas of study. Papers which could fit reasonably into more than one
category are listed only once.
r] Harmonic Drive Gearing: Cup Type HDUC and HIUC Component Sets,
HiD Systems, Inc., Hauppauge, NY.
173
Chavier 5: Reference List
11 Hidaka, T., Ishida, T., Zhang, Y., Sasahara, M., and Y. Tanioka,
"Vibration of a Strain-Wave Gearing in an Industrial Robot"
174
Chapter 5 Reference List
1 Ostapski, Wieslaw,
"Analysis of Some Aspects of the Harmonic Drive Dynamic",
Modeling, Simulation & Control A: GeneralPhysics (Matter & Waves),
Electrical& Electronics Engineering,
Vol. 9, No. 1, 1987, pp 21-33.
16 Revin, A. I.,
Dynamic Model of a Harmonic GearPower Transrn.ssion,
Proc. of Bauman MHTS (Moscow Higher Technical School), No. 408,
1984, pp 27-32. (Russian)
175
Chapter 5: Reference List
Russian EngineeringJournal,
Vol. 58, No. 5, 1978, pp 11-15.
21 Hsia, Lih-Min,
HarmonicGearDrive PositionErrorAnalysis,
Technical Report, Jet Propulsion Laboratory, August, 1986.
F5 Popov, P. K.,
Investigation of Angular Positioning Errorsof the Output Shaft of a
Harmonic Drive,
Thesis, MVTU, Moscow, 1972, 15 pages. (Russian)
176
Chapter 5: Reference List
31 Popov, P. K.,
"Torsional Stiffness of a Harmonic Drive",
Russia, M.V.S.S.O., Izvestia V. U. Z. Mashinostroenie,
No. 4, 1972, pp. 43-47. (Russian)
177
Chapter 5: Reference List
[7 Shamsutdinov, F. A.,
"The Transverse Stiffness of the Flexible Wheel of a Harmonic Gear
Transmission",
Russia, M.V.S.S.O., Izvestiya V. U. Z. Mashinostroenie,
No. 2, 1984, pp 36-40. (Russian)
f8 Shuvalov, S. A.,
"Calculation of Harmonic Drives with Allowance for Pliancy of the Links",
Russian EngineeringJournal,
Vol. 54, No. 6, 1974, pp 47-52.
178
Chapter 5: Reference List
4 Chiba, Jiro,
"Adaptive Control Cuts Position Servo Vibration",
PCIM,
June 1990.
179
Chapter 5: Reference List
180
Chapter 5: Reference List
59 Peter, J.,
"Investigation of the Engagement of Harmonic Drives - Part II",
Acta Technica Academiae Scientiarum Hungaricae,
Vol. 94, No. 3-4, 1982, pp 223-233.
181
Chapter 5: Reference List
F Shen Yun-Wen,
"On the Profile of the Harmonic Drive Gearing",
Chilun (Gearing),
Vol. 10, No. 4, 1986, pp 51-56. (Chinese)
65 Nye, T. W.,
"Harmonic Drives: Determining Wear Life Based on Stiffness
Considerations",
Proceedingsof the 1989 InternationalPower Transmissionand Gearing
Conference: New Technologiesfor Power Transmission in the 90's,
Published by ASME, New York, NY, 1989, pp 867-877.
182
ChapterS: Reference List
F9 Istomin, S. N.,
"Choosing the Geometrical Parameters of the Flexible Wheel of a
Hermetically Sealed Harmonic Gear Transmission",
Soviet EngineeringResearch,
Vol. 5, No. 7, 1985, pp 15-17.
183
Chapter 5: Reference List
HarmonicDrive,
Thesis, MAT!, Moscow, 1974, 168 pages. (Russian)
73 Diamond, Albert,
"Sizing Harmonic Drives",
Machine Design,
March 22, 1990
74 Emel'yanov, et al.,
"Rig Tests on the Feed Drive of an NC Machine Tool With Harmonic
Gear",
Soviet EngineeringResearch,
Vol. 2, No. 11, 1982, pp 71-73.
77 Gulyas, M.,
"Harmonic Drive in Special Purpose Machines and Materials Handling
Systems",
Proceedingsof the Second Seminaron Single PurposeMachines,
Published by OMIKK-Technoinform, Budapest, 1983, pp 43-55.
(German)
184
Chapter 5: Reference List
[8] Press, William H., Flannery, Brian P., Teukolsky, Saul A., and William T.
Vetterling,
Numerical Recipies in C: The Art of Scientific Computing,
Published by Cambridge University Press, New York, 1988.
185
186
This appendix describes the calibration procedures and results for all of the
experimental sensors.
A.1 Encoders
0. ..
. 0
~~~~~~~~~~~~......
. . . . . . . --- ---- ........
1000 1, 2, or4
2500 1, 2, or 4
1000 1,2.or 4
500 1, 2, or 4
A.2 Resolvers
The resolver position information was returned as a 16-bit number representing a
fraction of a complete rotation of the resolver shaft. Consequently, to convert this signal
into degrees of resolver rotation, the number of resolver counts was multiplied by (360.0
187
Appendix A: Sensor Calibration
degrees / 216 counts). The resulting value was then converted to degrees of output rotation
by factoring in the sensor gear-ratio of 6.66, on joint 2, or 7.2, on joint 3.
The resolver converter electronics could also return an analog voltage proportional
to the resolver velocity. Since this information was never used, the resolver velocity was
not calibrated.
A.3 Tachometers
The tachometers on joint 1 and 2 were calibrated by simultaneously collecting raw
velocity data, measured in A to D counts, and encoder position data, measured in degrees.
The tachometer data was then integrated to yield raw position data which was compared to
the actual encoder position data to derive a calibration factor. Figures A. 1 and A.2 show
the linear fits of integrated tachometer position data versus encoder data for the two
tachometers. The resulting calibration factors and their associated uncertainty are
summarized in table A.2.
18000 S................................................... T
............... ......... i......................... .•.... .............. . .. . . . ..
..
S1 6000
S~.......................... .................................................... ...................... ..... ... .....................................
S14000
• 12000.....
8000 ..........
...
0
0 200 400 600 800 1000 1200
188
Appvendix A: Sensor Calibration
32000 -
CL 20000....
S16000 ...
128000
......... ................ ........................................... .. .... ...
W 00
.. . . .. .. .... 0. .. . . . . . . .. . . . . . . . ..
.. ...... 1.504
4 Xi 16853 0.4%
189
Appendix A: Sensor Calibration
sensor iiiiiiiiiiiiiiiiii
Jon 3 torqueii
iii! 0.0065 0.o572 1.8 %
190
Appendix A: Sensor Calibration
800-'- - -
700 .........
................... .............. ....... ....... .......
......
5.....................................
........ ..............................
e
100 .......................................................
191 ...
...................
..........
Anvendix A: Sensor Calibration
-00
a-200 ...........................
........... . . . . . . .
0. ...... ....... ......... ........ .............. ............
c'-500
450
00...
.....4........
...... ...
...............
............
..... ...................
......... .......... ........... ..........
o*300 .......
0
- 2 5 0 ...... ......... .........
C6
2500 ...................... ý ............................ ...........
0 ...... ..... ..... ...... ............ .......................................
LL15 ..... ...........
1 00. ..... .
...... ........
.......... ..... ................
....
0
0 400 800 1200 1600 2000 2400
Torque Sensor Raw Data (counts)
192
Appendix A: Sensor Calibration
0
-1
00 ........... ................................ ................ ..... .....
-1 0 ................ ................................ ............. ......... ....
W................................ .......................... ... ...... .. .
...
............. ............. .. ............... ....
....... .......... ....
-2 00 ............................... ................ .............................
CO
1600
12
80
..
........................ ... ...................
................ ... .... .......
00
U.
193
Appendix A: Sensor Calibration
0
..................................... .....-
60......... ............. ......
I- 20 .............
0 ........................ .......... ..........
06
0....... ....................... ............
2....................-1............0
0).............
-120 ........................ ....
.
-1194
The computer code that was used to numerically solve the differential equations
describing the dynamic behavior of the different harmonic-drive models is listed in this
appendix. The header-file of variable definitions used by this program is also presented
after the main program listing. In addition, this program relies on a function, not listed
here, which can integrate ordinary differential equations. The function I used was a
Runga-Kutta solver with adaptive stepsize control taken largely from reference (85),
Numerical Recipes in C. During execution, this program reads the input parameters from a
file, calculates and integrates the specified equations of motion, and outputs time-response
data for the variable names specified in another file. I hope that this code can serve as
template for rapid development of future harmonic-drive simulations.
195
Appendix B: Simulation Code
/* Declare functions which calculate the total equations of motion for the
"* harmonic-drive system when the transmission is ideal, in this case, the
"* system reduced to a first-order system. */
void derivs idealconfigl(double, double *, double );
void derivs idealconfig2(double, double , double *);
/* Declare the individual functions which calculate the input parameters, dynamic
* equations and energy values for the different harmonic-drive models.
void calculate rotaryparams(void);
void rotaryhd model(double, double, double, double, double, double,
double *, double *, double *);
void rotary hd energy(double);
void calculategeartoothparams(void);
void gear tooth hd model(double, double, double, double, double, double,
double *, double *, double );
void gear toothnd enerqy(double);
196
Appendix B: Simulation Code
/* Declare the functions which (1) calculate stiffness data, (2) mooe. the
"* saturation limits of the current amplifiers, and (3) determine tne static
"• and dynamic friction values, respectively. /
void collectstiffnessdata(void);
double iregulate(double, double, double);
double calculatefriction(double, double, double, double, double,
double, double, double, double, double);
'* This main program opens and reads the input data files, initializes appropriate
* variables and prints/plots the results. */
main(int argc, char **argv)
int h, i, j;
int num of trialsiNUMOFAXIS) = (SHLDNUMOFDATARUNS,
ELBNUMOFDATARUNS,
WRISTNUMOFDATARUNS);
int goodsteps(l], badstepsil];
float input_amps(NUMOFAXIS] (SHLD NUM OF DATARUNS] =
(1.6, 2.0, 2.4, 2.8, 3.2, 3.6, 4.0, 4.4, 5.0, 6.0),
(1.4, 1.8, 2.2, 2.6, 3.0, 3.4, 3.8, 4.5, 0.0, 0.0),
(0.28, 0.36, 0.40, 0.44, 0.48, 0.52, 0.56, 0.60, 0.0, 0.0));
char unix command(200);
char vectorname[lO0);
char 'inputdatafile[NUMOFAXIS) - ("hd_shldinput data.dat",
"hdelbinput_data.dat",
"hd wristinput_data.dat");
FILE *infile;
FILE *outfile[NUM OF OUTPUTVAR];
extern int kount;
/* Open the input data file for the given joint, read the simulation
"* parameters from the file, and then close the input file. '1
open_file(&infile, indatapath, input_datafile(joint]);
read_input_data(infile);
closefile(infile);
/* Calculate necessary parameters for the given model from the values read
* from the input file. */
(*(calculateparams(model])) ();
197
Appendix B: Simulation Code
printinput parameters );
/* Set the pointer to the appropriate function containing the equations 'of
" motion. */
if(joint =- WRIST)
if(ideal_flag)
ptr-derivs = derivsideal config2;
else
ptr-derivs = derivsconfig2;
else
if(idealflag)
ptr derivs = derivsideal configl;
else
ptr-derivs = derivsconfigl;
/* Now if the largest line number in the index[] vector is greater than
"* or equal to the first line of the energy variables in the second
"* input data file, then energy is required as an output, so that the
" energy flag shuuld be set so that energy is calculated. /
if(index[numselected_output_var-li >= TOTALENERGY)
calc energy = TRUE;
/* If the energy is to be calculated, increase the order of the system
"* to ORDER+l in order to accomodate the additional total energy variable
" which is to be integrated as a member of the state variable. /
if(calc_energy)
system-order = ORDER + 1;
198
Appendix B: Simulatnan Code
if(debug)
/* :nitialize all output files, export the desired output data to them,
*and close the files. If the -plot option is selected, then save
,the specified data in the plot_data directory with the appropriate
*filename. */
for Ci - 0; 1 < num-selected-output var; i++)
If(plot)
else
close-file(outflletil))
/* NJow plot the results using xgraph, if slow printable plotting is desired,
* or gnuplot, if the fast plot option is selected. "
if(!plot)
if(fast plot)
else
I' Functions which describe the dynamics of the two joint configurations.
199
Appendix B: Simulation Code
I* Function derivs -configio) can be used to solve the equations of motion for :ý'c
"*harmonic drive joint configuration in which: (1) the motor drives the wave-
"*generator and is mounted to the circular spline, (2) the flexspline is attached
"*to ground, and (3) the circular spline is the output port. '
void derlvs-configl(t, x, dxdt)
double t;
double x[], dxdtf];
/* Set the local variable definitions to the value of the state var.~able
* determined on the previous step. *
v(POSIN] x[POS ININDEX];
v]POSOUT] = x(POSOUTINDEX];
v(VELIN] - x]VELININDEX];
v[VELOUT] = x[VEL OUTINDEX];
1* From these positions and velocities, call the appropriate narmonic drive
* model to calculate the resulting torques. */
(*(hd_model~model]))(v[POS_IN], 0.0, v[POSOUT],
v[VEL_ IN], 0.0, v[VELOUT],
6(v]TORQUE-WGI), 6(v[TORQUEFS]), &)v[TORQUE_CS;));
1* Set the value of torque sensor to the torque read on the flexspline. ~
v[TORQUESENSOR] -v[TORQUEFS];
200
Appendix B: Simulation Code
/* Function derivs -config2() can be used to solve the equations of moti~on for
" the
harmonic drive joint configuration in which: (1) the motor drives t-he
"* wave-generator and is mounted to the circular spline, (2) the fIexspllne :s
"* attached to the output link, and (3) the circular spline attached to 9roLuno.
void derivs-config2(t, x, dxdt)
double t;
double K]], dxdtil;
/* Set the local variable definitions to the value of the state variable
* determined on the previous step. ~
v[POSIN] - x(POSININDEX];
v]POSOUT] - x[POSOUTINDEX);
v(VELIN] -x(VELININDEX);
v(VEL OUT] x(VELOUTINDEX];
/* From these positions and velocities, call the appropriate harmonic drive
model to calculate the resulting torques. */
*
(*(hd_model~model]))(v[POSIN], v]POSOUT], 0.0,
v(VELIN], v(VEL OUT], 0.0,
&(v(TORQUEWG], &(v]TORQUEFS]), &(v(TOROUE_CS]));
/* The torque sensor is mounted between the motor and circular spline and
"*
ground, and therefore sees the reaction torque from the motor, the motor
"*
damping torque, and the circular spline torque. */
v[TORQUESENSOR] = v]TORQUE_CS] - v]TORQUEMOTOR] + v]TORQUEB_ IN];
201
Ampendix B: Simulaticn Code
/* Function calculate total energy() calculates the total power in zhe system
"* of equations presented in derivs configl() or derivs_config2(). This power
"* is determined on every time step and stored in the state vector whicn is
" is solved by the Runge-Kutta solver. This can be used to verify that energy
"* is being conserved in both dynamic models. /
void calculatetotal energy(x, dtdt)
double x[], dxdtf];
/* Now calculate the input damping power loss depending on the joint
* configuration. Also calculate the appropriate torque on the output
* port of the harmonic drive. /
if(joint == WRIST)
else
202
Appendix B: Simulation Code
" the system, excluding the harmonic drive. Since this function relies on
" values calculated in the harmonic drive energy function, it should be called
" after the specific harmonic-drive energy function is called. Due to the
"* similarities between the equations in derivsconfigl() and derivs_confi92(),
"* this energy function is valid for both models. /
void calculate configenergy(delta t)
double deltat;
/* Now print out the energy values if the print energy option is selected. /
if (printenergy)
203
Appendix B: Simulation Code
204
Appendix B: Simulation Code
print variables );
/* Set the local variable definitions to the value of the state variabie
* determined on the previous step. /
v[POSIN] = x[POS ININDEX];
v[POS_OUT] - x[POS OUTINDEX];
v[VELIN] = x(VELIN_INDEX];
v[VEL_OUT] - x[VELOUTINDEX];
/* Calculate the motor torque as a function of current. For the ideal
* model, assume that the current amps function ideally as well.
v[CURRENT SENSOR] - irequested;
v[TORQUEMOTOR] = motor_kt * v]CURRENTSENSOR];
/* The torque sensor is mounted between the motor and circular spline and
"* ground, and therefore sees the reaction torque from the motor, the motor
"* damping torque, and the circular spline torque. */
v(TORQUESENSOR] - v(TORQUECS] - v(TORQUEMOTOR] + v[TORQUE_B_IN];
205
Appendix B: Simulation Code
if (debug)
print-variables 0:
I. Store the velocities and positions and the fs, cs, and wq ports
in the global data vector to begin calculations. ~
*
v[POS_WG] - pos~wg;
v[POS_FSI - pos_fs;
v[POS_CS] - p05 Cs;
vjVEL WG] - vel wg;
v[VELFS) - vel fs;
vfVEL CS] - vel-cs;
/* Determine the kinematic error function and its derivative such Chat
(d/dt)v[ERFNI - vIVEL._NC]I vtDERF). *
*
v[ERFNI - ((error_amplitudeO Sin(((1.0 * v[POSWG]) + error phase0) *CONy))
(error-amplitudel *Sin(((2.0 * v[POS NC]) + error phasel) *CONV)) +
(error amplitude2 Sin(((4.0 * v[POS WGJ) + error_ohase2) *CONVMl;
v[DERF) - ((1.0 * CONV * error amplitude0 *
Cos(( (.0 * v[POS NC]G) + error phaseOl CONy)) +
(2.0 * CONV * error amplitudel *
Cos(((2.0 * v(POS_-WG]) + error phasel) *CONV)) +
(4.0 * CONV * error amplitude2
CosM(4.0 * v[POSNC]) + error-phase?) *CONy));
/* Calculate the static, dynamic, and cyclic friction which acts between the
*circular spline and flexspline. */
206
Apnendix B: Simulation Code
/* For the sake of continuity in the simulation, ramp the constant friction
"* in the harmonic drive from zero to the desired value if at the
" beginning of the step-response trial. '/
if(!(stiffness flag 11 stiffnessflag2))
if(currenttime < 0.05)
v[TORUEHDFRICTION] - ((((current time - 0.05)/0.05) + 1.0)
v(TORQUEHDFRICTION]);
/* Store the velocities and positions and the fsr cs. and wq ports
207
A~ondi
B:SimlatonCode
/* Determine the kinematic error function and its derivative such that
* (d/dt)v(ERFN] v[VELWCRELATIVE] v[DERF]./
v(ERFN) - (N*tanl *((error amplitudec *
SinM1(.0 *v(POS_WC_-RELATIVE]) + error_phaseD) *CONV))
(error amplitudel*
Sin(U2.0 *v[POS - C_-RELATIVE]) + error pnasel) *CONy))
(error_amplitude2*
Sin((4.0 * v(POSNC_-RELATIVE]) error_phase2) * CONVM));
v[DERF] (N*tanl ((.0 * CONV ,error amplitudeD
Cos((C1.0 v[POSNC_-RELATIVE)) + error_phase0) CONV))
(2.0 * CONV *error-amplitudel*
Cos((2.0 * v(POS NCRELATIVE)) + error_phasel) *CONV))+
(4.0 * CONV * error amplitude2 *
Cos((4.0 * v[POSNCRELATIVEI) + error-phase2) *CONV)));
/* Find the radial movement at the output of the position error element.
v[POSERR] - v(POSTOOTHBASE) + v[ERFN);
v[VELERR] - v(VELTOOTHBASE] + (vIVELNCRELATIVE] v[DERF]);
208
Appendix B: Simulation Code
/* The total friction at the gear-tooth surface is due to the coulomb friction
"* (mu * N), the constant friction (btoothsurface constant), the velocity-
"* dependent friction, and the cyclic friction. Calculate the component due
"* to everything except the coulomb friction. */
v[TORQUE TOOTHSURFACEFRICTION) -
calculate friction(v(VELTOOTHSURFACE],
b tooth surface_constant, btooth_surfacel, btooth_surface2,
stiction torquetooth, stictionvel.tooth,
v(POSOUT],
cyclic friction_amp_tooth, cyclicfrictionphase_tootn,
(maxfriction vellmodel)[joint)));
/* Set the sign of mu to control the direction of the coulomb friction force
"* along the gear-tooth surface. The sign of mu should be identical to the
"* sign of the surface rubbing-velocity. ./
mu - copysign(mu_save, v[VELTOOTHSURFACEI);
/* If stiffness data is being collected, the friction acts in the opposite
* direction as the joint is being loaded. 1
if(stiffness flag)
mu = -mu save:
/I For the sake of continuity in the simulation, ramp the constant friction
" on the gear-tooth surface from zero to the desired value if at the
"* beginning of the step-response trial. This is a quick fix that should
" be implemented in more general terms for all of the constant friction
* components. */
if(!(stiffnessflag II stiffnessflag2))
if(current time < 0.05)
v(TORQUETOOTH_SURFACEFRICTIONI - ((((currenttime - 0.05)10.05) + 1.0) "
vITORQUE TOOTH SURFACE FRICTION));
209
Appendix B: Simulation Code
/* Find the gear-tooth normal force. Note that this equation might need to
" be calculated iteratively if the velocity is zero since the coulomb
"* friction becomes a different function of the normal force. For now, 1
" will start my simulations at a non-zero initial velocity so that the
" velocity will never be zero. */
vITORQUETOOTHSURFACENORMAL] = ((v[TORQUEOUTPUTNORMAL] -
(cos2 * v[TORQUE TOOTHSURFACEFRICTIONM)) /
(sin2 + (mu , cos2)));
mu - 0.0;
v[TORQUETOOTHSURFACENORMAL] = ((viTORQUEOUTPUT NORMAL] -
(cos2 * v[TORQUE TOOTHSURFACE _FRICTlONIII /
(sin2 + (mu * cos2)));
/* Given this correct tooth-surface normal force, the coulomb friction and
" total friction at the tooth interface can be calculated. */
v]TORQUE TOOTH SURFACECOULOMB] = mu * v(TORQUEY'OOTHSURFACENORMAL];
v[TORQUE TOOTHSURFACETOTALLOSSI = (v(TORQUETOOTH SURFACE COULOMB]
-v[TORQUE TOOTHSURFACEFRICTION]);
/* Now the torques on the wave-generator, flexspline, and circular spline can
* be determined. */
v(TORQUE-WG] = ((v]TORQUECYCLIC]) +
(v[TORQUEWGFRICTION] * cosl) +
(v[TORQUE WG NORMAL] * sinl));
v[TORQUEFS) - ((v[TORQUE TOOTHTIPNORMAL]) -
(v[TORQUE WG FRICTION] * cosl) -
(v(TORQUE WG NORMAL! * sinl));
v[TORQUE CS] - ((v]TORQUE TOOTHSURFACENORMAL] * cos2) -
(v]TORQUETOOTHSURFACETOTALLOSS] * sin2));
/ *********************************.*******.**t***..............
210
Appendix B: Simulation Code
./ . ..............
....................
/* Calculate the cyclic torque energy by integrating the power over time. /
v[CYCLICENERGY] = (v[CYCLICENERGY] +
(deltat * velwgrad * (-v(TORQUECYCLIC])));
/* Calculate the energy lost to friction by integrating the power dissipation
* over time. */
v[HDFRICTIONENERGY] - (v[HDFRICTIONENERGY] +
(delta t * (v[TORQUEHD_FRICTION! I vel hd friction rad)));
/* Calculate total potential energy in the system assuming the cyclic torq.e
"* can be treated like a spring. */
v[TOTALPOTENTIALENERGY) = v[SPRINGENERGY) + v]CYCLICENERGY;;
211
Appendix B: Simulation Code
/* Calculate the cyclic torque energy by integrating the power over time. "/
v]CYCLIC_ENERGY] - (v[CYCLICENERGY] +
(delta t * vel_wg_rad * (-v[TORQUECYCLICI)));
/* Calculate the energy lost to friction by integrating the power dissipation
* over time. */
v[WGFRICTIONENERGY] = (v(WGFRICTION_ENERGY]
(deltat * (v[TORQUE WGFRICTION] * vel_wg_surfacerad)));
v[TOOTHTIPFRICTIONENERGY] = (v[TOOTHTIPFRICTIONENERGY] +
(delta t *
(v]TORQUETOOTHTIPFRICTION]
(veltoothbaserad - veltooth_tip_ rad))));
v[TOOTHSURFACEFRICTIONENERGY] - (v(TOOTHSURFACEFRICTION ENERGYI +
(deltat *
(v[TORQUETOOTHSURFACEFRICTION] *
vel tooth surface rad)));
v[TOOTHSURFACECOULOMBENERGY] = (v[TOOTHSURFACECOULOMBENERGY] +
(delta t *
(v[TORQUETOOTHSURFACECOULOMBI *
vel tooth surface rad)));
v[TOOTH SURFACETOTALLOSSENERGY] = (v[TOOTHSURFACECOULOMBENERGY] +
v(TOOTHSURFACEFRICTIONENERGY]);
/1 Calculate total potential energy in the system assuming the cyclic torque
* can be treated like a spring. -/
212
Appendix B: Simulation Code
/- Calculate the coulomb friction in the harmonic drive based on the width
* of the hysteresis loop in the harmonic-drive stiffness profile. /
b hd constant = (hysteresiswidth / 2.0);
/* Calculate the reduction ratio between the input rotation and the veiocity
* seen by the friction element mounted between the flexspline and circular
"* spline ports. */
if(joint -- WRIST)
reduction-factor - N;
else
reduction-factor = (N + 1.0);
/* Convert the stiction and cyclic friction from input rotation units to
* output rotation. Don't upset the signs of these values. /
stiction torque_hd = reductionfactor * stictiontorque;
stiction vel hd - (stiction vel / .eductionfactor);
cyclic_frictionamp hd - reduction factor * cyclicfriction amplitude;
cyclic_friction phasehd = cyclic frictionphase;
1* Function calculate gear tooth params() is used to take the input parameter
" values and determine: (1) trigonometric and geometric relationships, (2) the
213
Appendix 8: Simulation Code
"* tooth surface friction coefficients, and (3) the ideal spring stiffness
"*coefficients. */
void calculate geartoothparams()
double stiffness(factorA;
double stiffness factor B;
double wg to tooth factor;
double damping factor A;
double dampingfactorB;
double dampingfactor C;
double dampingfactor D;
double inertia in eff, inertiaout_eff;
double templ, temp2, frequency;
/* The wgangle can be found from the gear ratio and the toothangle
* using the relationship: (1/(N+l)) = tanl tan2. */
wg angle - atan(l.0 / (tan2 * (N + 1.0)));
tanl - tan(wg_angle);
sinl - sin(wgangle);
cosl - cos(wgangle);
/* From the hysteresis width on the stiffness curve and the values of the
"* other constant friction components in the harmonic drive, determine
"* the appropriate constant friction component on the gear tooth surface. /
btooth surface_constant - ((1.0 / ((stiffnessfactorA * cos2) + sin2))
((hysteresiswidth / 2.0) -
((stiffness_ factor A - tanl) - b_toothtip_constant) -
/* Scale the stiction torque and cycliý. friction torque from the input
* side (wg) to the gear tooth surface. */
wgto_tooth factor - (1.0 / ((tanl cos2) + (tanl * sin2 * tan2)));
stiction_torque_tooth - wgto_toothfactor * stictiontorque;
stictionvel tooth - (stiction vel / wg to tooth factor);
cyclic_frictionamp tooth - wg-totoothfactor * cyclic friction amplitude;
cyclic_frictionphasetooth - cyclicfrictionphase;
214
Appendix B: Simulation Code
"*
done differently for configuration 1 and configuration 2 since the
"*
dampers outside the harmonic drive depend on different velocities./
damping_factorA = ((cos2 / tan2) + sin2);
damping_factorB -(cosi + (sini tanin);
damping_factor_C - (1.0 / (N + 1.0));
damping factor_D - (1.0 + (1.0 / N));
lf(joint -- WRIST)
b-tooth-surfacel =((b-totall -
b-in -
(b~wgl *(Power((damping factor_B damoing factorD), 2M
(b-out *(Power((-1.0 / N), 2))))/
(Powercidamping factorA *(1.0 / N)), 2));
b-tooth-surface2 =((b total? -
(b~wg2 * (Power((damping factor_B * damping factorD), 4M)))
else
b-tooth-surfacel =((b-totall -
b in -
(b~wgl *(Power((damping factor_-B), 2))M
(b-out *(Power((damping factor_-C), 2M)))
(Power((damping_factor_A * damping factor_C), 2));
b-tooth-surface2 =((b-total2 -
(b~wg2 * (Power((damping factor_-B), 4M)))
(Power((damping~factozA *damping factor_C), 4));
215
Appendix B: Simulation Code
/* Functions which model the friction and amplifiers and collect stiffness aata. "/
/* Function irequlate() takes as arguments the present amp current, the requested
* amp current, and the current motor velocity and returns the actual current ::a:
* the amps can produce. */
double iregulatecipresent, irequested, omega)
double ipresent;
double irequested;
double omega;
double voltageupperlimit;
double voltage_lowerlimit;
double current upper limit;
double current lower-limit;
double currentl, current2;
double ireturn;
/* Make sure that the current requested is within the amp's operating range. */
if (irequested > IMAX)
else
/* Now from this maximum or minimum allowed voltage at present current, calculate
216
Appendix B: Simulation Code
"* the allowable current range thatthe requested current can achieve by assuming
"* that the amp voltage is equal tothe voltage across the motor (- Kb * velocity
" or back EMF) plus the voltage drop due to current flowing through the amp
"* resistance. */
currentl - (voltageupper limit - (motorkb * omega))/amp resistance;
current2 -(voltage_lower_limit - (motor-kb * omega))/amp resistance;
if(currentl < current2)
else
/* Now, if the requested current, irequested, is between the upper and lower
* current bounds, return the value, otherwise return the upper or lower
* limit. */
/* Now add a relative damping factor to the current. Make it behave so that
"* the current can't change from the old to the new value instantaneous.y, buL.
"* instead, a part of the new value, ireturn, is averaged with part of tne olo
"* value, ipresent, to create the actual new current somewhere between the two
"* values. This step is necessary to keep the amps current from fluctuating
" rapidly between each subsequent step and possibly going unstable. The AMP
"* DAMPING factor is a number between 0 and 1.0 that represents the amount
"* of damping in the amps current response. Note that since step size can vary
" significantly, and that this damping factor is calculated on every step, tne
"* amount of damping will vary with the changing step size. */
ireturn - (AMPDAMPING * ipresent) + ((1.0 - AMPDAMPING) * ireturn);
/* Give a run-time warning if the current exceeds the amp's allowed range. I/
if ((ireturn > IMAX) II (ireturn < -IMAX))
return(ireturn);
217
Appendix B: Simulation Code
"* bconstant or else the resulting friction can become negative. Additionally.
"* if the velocity is very small, this friction model is undefined, and a warning
" message is given. Finally, if the velocity is too large, the cubic term of the
"* velocity-dependent damping may decrease and become negative. To prevent this,
"* the velocity dependent damping is calculated for omega max at all velocities
"* greater than omega max. */
double calculatefrictlon(omega,
b_constant, bl, b2,
stictiontorque, stictionvel,
theta,
friction_amplitude, friction_phase,
omegamax)
double omega;
double b-constant, bl, b2;
double stiction torque, stictionvel;
double theta;
double friction amplitude, frictionphase;
double omegamax;
/* If the stiffness data is being collected and the joint is being loaded,
"* return the velocity-independent friction as if the velocity were positive.
"* Note: do not use the stiction torque. *I
if(stiffness flag)
friction = (-bconstant +
(-frictionamplitude
Sin((theta + frictionphase) * CONV)));
return(friction);
friction - (b constant +
(friction amplitude
Sin((theta + frictionphase) * CONV)));
return(friction);
/* Make sure that the cylic friction amplitude is less than the constant friction
* coefficient. */
if (friction-amplitude > b_constant)
218
Appendix 8: Simulation Code
int i;
double torque[STIFFNESSPOINTS];
double displacement[STIFFNESSPOINTS);
double pos_increment;
double actual disp, torque wg, torque_fs, torque_cs;
char unixcommand(200J;
int points - STIFFNESS POINTS;
static double maxdisp[NUM OFAXISJ - (0.1, 0.4, 0.35);
static char *stiffness filename(NUMOFAXIS) - ("shld stiffness data.dat",
"elbstiffnessdata.dat",
"wrist stiffnessoata.dat");
FILE *outfile, *gnufile;
219
Appendix B: Simulation Code
if(joint -- WRIST)
else
if(joint -- WRIST)
else
displacementliI - fabs~actual-disop;
if(debug)
print_variables (;
if (last plot)
220
Appendix B: Simulation Code,
else
exit (0);
mnt arg_index;
for (arg index - 1; arg index < argnum; arg~index++)
else
221
Appendix B: Simulation Code
/* Function usagec) prints out the proper usage of the command-line options. I/
void usage()
/* Function read input datac) parses all of the values in the input data file.
void read-input-datalf-ile pointer)
FILE *file-pointer;
char temp_name(5OI;
double temp_number;
while(fscanf (file pointer, 1%s %lg %*s\n", temp name, &temp_number) !- OF)
if (0 -- strcmp("Iinertia_in", -temp_name))
inertia-in - temp number;
else if (0 -- strcmp("inertia -out", temp name))
inertia out - temp number;
else if (0 __ strcmp(--kl_output_constant"-, temp name))
kl -output_constant - temp number;
else if (0 -- strcmp("kl -output amplitude", temp-name))
kl output amplitude - temp number;
else if (0 -- strcmpt11kl_output phase". temp name))
kl -output-phase - temp number;
else if (0 -- strcmp(*"k2_output -constant", temp name))
k2 -_output constant - temp number;
else if (0 -- strcmp(1'k2_output_amplitude11, temp_name))
k2 -output amplitude - temp number;
else if (0 -- strcmp1"k2_output _phase11, temp name))
Ic?output phase - temp number;
else if (0 -- strcmp(I-hysteresis -width", temp name))
hysteresis_width = temp number;
else if (0 -- strcmp("b-totall", temp name))
b -totall - temp number;
else if (0 -- strcmp("b -total2", temp_name))
b -total2 . temp_number;
else if (0 -- strcmp("lb in", temp name))
b~in - temp number;
else if (0 -- strcmpj11b~wgconstant", temp_name))
b~yq_constant - temp number;
else if (0 -- strcmp("b~wgl", temp name))
b~wql - temp number;
else if (0 -- strcmp("b~wg2", temp name) I
b~wq? - temp number;
else if (0 -- strcmp("b~t 0th tip constant", temp name))
b~tooth_tip_constant - temp number;
else if (0 -= strcmp(I"b -tooth tipi", temp name))
btooth~tipl - temp number;
else if (0 -- strcmp("b tooth tip2", temp name))
b-tooth tip2 - temp number;
222
Avvendsx B: Simuiation Code
/* Initial conditions, */
else if (0 -= strcmp(I"initial pos in"., temp name))
xstart~sav[POSIN_INDEXI - temp number;
else if (0 -- strcmp(I"initial~pos~out', temp name))
xstart-savIPOS_-OUT_-INDEX] - temp-number;
else if (0 =- strcmp("initial_vel_in", temp-name))
xstart -savIVEL-IN_-INDEX) - temp number;
else if (0 =- strcmp("initial -vel_out", temp name))
xstart sav[VEL-OUT_ INDEX) - temp_number;
/* Time conditions */
else if (0 -= strcmp("initial~time**, temp name))
initial -time - temp number;
else if (0 =-strcmp("final -time", temp_name))
final -time -temp number;
else if (0 -- strcmp("'step*", temp-name))
step = temp number;
223
Appendix B: Simulation Code
else
/* Function readinputfilenames() reads the data from the second input file
"* as follows. First it checks to see if the first value on each line is 1.
"* If the value is 1, then it stores the line number in the vector index[]
"* so that it can be used later to remember which output variables to store.
"* Additionally, this function also stores the other informaiton from each line
"* that has a I value to be used later for graph titles and output Filenames.
"* As this function steps through each line of the input file, it counts the
"* number of output variables selected and stores the result in the variable
"* numselected outputvar. */
void readinputfilenames(file_pointer)
FILE *file_pointer;
int linenumber - 0;
int tempflag;
char temp title[50];
,char temp units[50];
if(temp flag)
line-number++;
/* Function print input parameters() prints out the names and vaiues of
* the input parameters. */
void printinputparameters()
printf("\n");
printf(" Model Parameters:\n");
printf(" Input Inertia - %4.101f\n", inertia in);
printf(" Output Inertia = %4.101f\n", inertiaout);
printf(" kl constant - %4.101f\n", kl_outpuc_constant);
printf(" kl amplitude - %4.101f\n, kl_outputamplitude);
printf(" kl phase - %4.10f\n", kl_outputphase);
printf(" k2 constant = %4.10f\n", k2_output constant);
printf(" k2 amplitude - %4.10f\n", k2_output amplitude);
printf(" k2 phase - %4.101f\n", k2_output phase);
printf(" Hysteresis Width - %4.101f\n", hysteresis width);
printf(" b totall - %4.10f\n", b totall);
printf(" b total2 - 14.101f\n", b total2);
224
Appendix B: Simulation Code
printf("\n") ;
printf("\n pos_in: %20.151t"., v[POS IN]);
printt("\n pos_out: %20.151f"l, v(POS OUT]);
printf("\n pos_wg: %20.151f", v(POSWGJ);
printf("\n pos_fs: %20.151fl", v(POSFS]);
printt("\n pos_cs: %20.151f"l, v[POSCSI);
printf("\n pos_n_wg %20.151fll, vIPOSNWGJ);
printt ("\n") ;
printt("\n vel in: %20.151fll, v(VEL_ IN]);
printt("\n vel out: %20.151f"l, v)VELOUTJ);
printf!"\n vel-wg: %20.151fl', viVEL_-WGI);
printf("\n vel fs: %20.151fll, v)VEL FSJ);
printf("\n vel cs: %20.151f", vIVEL CS));
printf C"\n") ;
printf("\n torque-motor: %20.151f", v[TORQUEMOTOR));
printf("\n torque~bin: %20.151lf", v[TORQUE_B_-IN]);
printf("\n torque_b -out: %20.151f", vITORQUE_BOUT));
printf("\n torque wg: %20.151f", v[TORQUE_WGI);
pritt(\ntorue a:%20.151f",- v(TORQUE FSI);
225
Appendix B: Simulation Code
mnt i;
ifC!fast_plot)
fprintf(outfile, "1\" %s\n". filename);
for Ci -1; i <- points; i++)
fprintf(outfile. "%12.101f %12.10lf\n", xvectortil, yvectorlil);
printf C"done.\n$);
1* Function generate gnuplot file)) is called when the fast _pIlt OPtLIOr
"* is selected and creates a command file in the data directory Which
"* is executed in gnuplot by a system command./
void generate_gnuplot-file(filename)
char filename));
int 1;
FILE gqnufile;
226
Appendix B: Simulation Code
227
Appendix B: Simulation Code
1* Define the number of data runs to be made if the -plot option is selected.
tdefine SHLD NUM OF DATA RUNS 10
tdefine ELBNUMOF_DATARUNS 8
#define WRISTNUMOFDATARUNS 8
/* Define the indicies for the vector to hold all variables that can be
* selected for output data plots. */
228
Appendix B: Simulation Code
#define POSWG 0
#define POS FS 1
#define POS CS 2
#define POSTOOTHBASE 3
#define PUS ERR 4
#define POSTOOTHTIP 5
#define POSWGSURFACE 6
#define POSTOOTHSURFACE 7
#define POSWGRELATIVE 8
#define POSK 9
#define POS_N_WG 10
#define POS_N_FS 11
#define POS_N_CS 12
#define VEL WG 13
#define VEL FS 14
#define VELCS 15
#define VEL TOOTH BASE 16
#define VELERR 17
#define VEL TOOTH-TIP 18
#define VELWG SURFACE 19
#define VELTOOTH-SURFACE 20
#define VEL WG RELATIVE 21
#define VEL K 22
#define VEL N WG 23
#define VEL N FS 24
#define VEL N CS 25
#define VEL HE FRICTION 26
#define TORQUE WG 27
#define TORQUEFS 28
#define TORQUECS 29
#define TORQUEINPUTNORMAL 30
#define TORQUECYCLIC 31
#define TORQUE WG NORMAL 32
#define TORQUE_WGFRICTION 33
#define TORQUEERRIN 34
#define TORQUEK 35
#define TORQUETOOTHTIP NORMAL 36
#define TORQUETOOTHTIPFRICTION 37
#define TORQUETOOTHSURFACENORMAL 38
#define TORQUETOOTHSURFACETOTALLOSS 39
#define TORQUETOOTHSURFACECOULOMB 40
#define TORQUETOOTHSURFACEFRICTION 41
#define TORQUECYCLIC FRICTION 42
#define TORQUEOUTPUTNORMAL 43
#define TORQUE_N_WG 44
#define TORQUE_N_FS 45
#define TORQUE N CS 46
#define TORQUE-HDFRICTION 47
/* Sensor variables. */
#define INPUT POSITION SENSOR 55
#define OUTPUT POSITIONSENSOR 56
*define INPUTVELOCITYSENSOR 57
#define OUTPUT VELOCITY SENSOR 58
#define CURRENT SENSOR 59
229
Appendix B: Simulation Code
#define TORQUESENSOR 60
#define POSITIONERROR 61
I/ Define a vector to hold all of the possible output data values for
230
Appendix B: Simulation Code
/* Define the parameters that are calculated from the input parameters
* for the rotary harmonic drive model. */
double klconstant, kl_amplitude, kl_phase;
double k2_constant, k2_amplitude, k2_phase;
double stiction torque_hd, stictionvelhd;
double cyclic_frictionamphd, cyclic_frictionphase_hd;
double bhd constant, bhdl, bhd2;
/* Define the parameters that are calculated from the input parameters
* for the tooth-rubbing harmonic drive model. */
double tanl, cosl, sinl, tan2, cos2, sin2;
double wg_angle;
double stiction torquetooth, stiction veltooth;
double cyclic_frictionamp-tooth, cyclic_frictionphase tooth;
double btoothsurfaceconstant, b tooth surfacel, btooth_surface2;
231
232
For most of the operating range of the harmonic-drive testing equipment, the
amplifiers which supply current to the DC motors can sustain current values up to the their
maximum rated current. However, at high motor velocities, back EMF generated by the
motor can saturate the amplifier voltage which results in a reduction of the maximum
current the amplifier can deliver. This saturation behavior can be captured by developing a
few simple equations which incorporate the velocity and current limits of the amplifiers
with some simplified dynamics of the motor-amplifier circuit.
Identical Aerotech model 4020-LS linear current-amplifiers were used for the three
harmonic-drive testing stations. The voltage and current operating range of these amplifiers
is shown in figure C.. Using this plot, the upper and lower voltage limits of the amplifier
i'K
A
...................
233
Appendix C: Motor Amplifier Model
can be calculated for any desired current within the given current range, as illustrated in
figure C.2. Given this upper and lower voltage-bound, an upper and lower current-bound
can also be found by considering the simplified amplifier circuit shown in figure C.3.
Summing the voltage around the circuit yields the equation
where Kb is the back EMF constant of the motor, ct)o or is the motor velocity, and R is the
resistance across the motor terminals. Resistance values which provided good results are
listed in Table C. 1 for the three harmonic-drive testing stations. By solving equation C.2
for the amplifier current, iamp, and substituting the voltage upper and lower-bounds as
determined in figure C.2, an upper and lower bound for the amplifier current can be
2 341 3.1 i m.
234
Appendix C. Motor Amplifier Model
V.
derived:
Vwpp max - KbO)motor and
inWmx-R ,ad(C.3)
Therefore, for a given motor operating speed and a desired current value, the instantaneous
current limit on the amplifiers can be determined simply by (1) using the amplifier operating
specifications to calculate the voltage limits for the requested current, and (2) substituting
these voltage limits into equations (C.3) and (C.4). This model was incorporated
successfully into a numerical simulation which calculated the amplifier current on every
time-step.
235
236
This appendix contains a complete set of all of the dynamic-response plots collected
from the three harmonic-drive testing stations. For each harmonic-drive, input and output
position were measured and used to derive input and output velocity data. The output
torque on each transmission was monitored by the torque sensors while the input torque
was determined from the motor-current sensor. Additionally, from the input and output
position measurements, the dynamic position-error across the transmission was calculated.
A summary of all of the plots shown in this appendix and their corresponding plot numbers
is given in table D.l. For each of the plots listed in this table, except the dynamic position-
error figure, time-response curves are plotted for several different motor-current step-
commands. Depending on the visual appearance of the data, the different response curves
on each of these plots are either plotted on a single axis or plotted on separate axes. For the
sake of clarity, the output-velocity response data for each harmonic drive is plotted in both
of these ways. In the case of the dynamic position-error plots, a single time-response
curve is shown, rather than several, in order to generalize about behavior for different step
inputs. A complete discussion of all of these plots can be found in section 2.6 of this
document.
237
Appendix D: Experimental Data Plots
238.....
Appendix D: Experimental Data Plots for Joint I
..........
........
.............. 6.0 am ps
80000
...................
440 amps
S4 4 . .4......... ..
............ ... ... ........ ...f............ ... ... .. ...
40000/404 amps
~...............
.... . . . . . . . . . . .
. . . . . . 2.8 am ps
40000 -1..
~~~~~~2.
...............
................
....... ................... amps........
. ...
........ ... ...
239. ....
..... 28a p
Appndx
-Exermeta Data Plots for Joint 1
25000 f 1 rT
6.0 amps
20000 -
CD 15000 44ap
4.0 amps
CD
50000
~ 2.4 amps
Pp
. 2.0 amps
50 0 ........ L iz J. L ii J.
0~. 1a2p3s
1. amps
240
A ppendix D: Experimental Data Plots for JointI
4.4 amps
, L
E 4.0 amps
3.6 amps
~ 3.2 zmps
* * - - -1 .6 amps
0 12 34
time (seconds)
241
D Exprimetal Data Plots for Joint I
Appenix
....
..... ...6.0 am ps
........................ ....
I.. .........
...
......
/..
.................
.................
.... 4 . 4..... .............
1 1.40 amps
200 "/....
:....
/ .2 ..... ...
... amps..
0
0... 1 ... ...... 4...
2..... ... 3........
time..
(seconds
.........
....... .. .... * ....
24...........
2.........
...... ......
..............
...
A opendix D: Exp~erimental Data Plots for Joint I
140
120 am.
5. ap
100. .. ...
......
. .......
14 -! - : .......
0 10 ........... 4. 0 a p
....
Z. 80 '"....:-
0
..
....
.
00
40 . . . ....
2.4 amps
20
1.6 amps
0
0 12 3 4
time (seconds)
Fiur D..5 J5ll I I~ t se te po 0Lt(tV
243
Appvendix D: Experimental Data Plots for Joint 1
150.-
158
158 ± ±
0 --
. 4 ...... 4I.....
Le&
:3 15
....... ... ... ................. ...
E ....
....
...
.... 24...
............. ........ .............. ....
Appendix D: Experimental Data Plots for Joint 1
8.
a I~~4 4 * .............. ........ ........ ...
co
A.i....L. .......
...........................
a.......
1...
......... .... ............. T..... ........
0e1 wua-m rS rrmt~IdA~
V: ......... ........................
...... ... ... ...
I....
... .... 4...
..........
a.............. ....
............ ................. .........
jt
o8
. . . . .. .a. . .
........ ... ..............- .......... ...
- 1 2 3ll 4,
tie(ecns
Y .......... .. ....... .........
.. . . . . . .245. .
Appendfix D: Experimental Data Plots for Joint I
0.1..........................
....
.-
0.1....... . . . . . . ..
0 0
246.. ........
Appendix D: Experimental Data Plots for Joint 2
160000 .. 4.5.amp
140000 .........
3.8 amps
100 -. ................ .................................................
0100000 .-......
.. .......
100000-~4
0~
224 amps
........ .......
.............
... ............
247...
..........
......
.. ........
Appendfix D: Experimental Data Plots for Joint 2
50000
.... ?.......... .... 4.
4..4....... ....... .................
....... ..... 4.5 am ps
40000
1...4.
i ...... ......... . . . . . . . . . . .
0
3.4 amps
300
300....... ....................
...................... /.........
....... ...
00
.. ..................
. .4.
. ......... 4. 4. ..........
........
10000MI-
2.8 amps
'~ T 1.2 amps
0
0 12 3 4
time (seconds)
248
Appendix D: Experimental Data Plots for Joint 2
10,. 11
0........
..... . ...
... ... . ... .. . ...
.. .. . ..
. . . . . . . . .... .
.. ... ... ... ... . . ... ...
.........-.................................. I
...........-.....
6 .................. .. d....
.. ...
..... ........ ..... ........
E : :; : : :249
Appendix D: Experimental Data Plots for Joint 2
800
700..... t .. ....
t....I........ ... ...
..................... ....
....... ... .... ....
1 ...........
0....
0............................................. ......
500
... ......
300-/
....... ....... 30a p
................................
5000
Appendix D: Experimental Data Plots for Joint 2
3500
...... ... ..
a ........
0. ap
30 1...
:...4.........
..... 2........... 4........j
3....
) ... .....
...... ....
3. amps
250 .. 2.1
Appendix D: Experimental Data Plots for Joint 2
350w-----
...- ...... ... -- ... ... .... .... ...
+ 4
i ............
35......... .......
.. J ...
......L....
r ;1
..........-
i...-f
-..............
r....j....
............
......... .............
....
...............
...........
...
...
....
.a........ .........
...
.................
358-J- 1 5 4 -
.............. ........
................
.4.......
............... 4 -4- 4 - - - .-..
.
..................... -
lip.....
V 358
;_-
....
....
............
T
....
....
.
....
T
........
....
..
..
I7I
3582 -rr r - 1
0 12 3 4
time (seconds)
D4. 0ilm
on [tLtV'0~ y ýt' -e, ole
252
Appendix D: Experimental Data Plots for Joint 2
20
±...............................................
. ..
- .Q ;........."
7wU.7. .....
..... .... . .
........... ..........................
o0 .............
. . . . .....
.............
O
20.. .............
... .....
2 ....
.... ......
.4................
0V
~~ ~
. ~ ...
~ ........
~ ~ ........
tm (seconds)............
........ ........ ....
....
.......
.........T ............
.... .... .... ........ .... ... ... ...
..
....
.... .........25 3
. ....... . ....
.... ....
..
Appendix D: Experimental Data Plots for Joint 2
0.15
.;.
0.1~ ~ *..* .**4
--- . . . . . .. . . .
4)a)............
.....
0..... .....................
. ...... . .............. .................. ..
....... ............ 4.....4.. .4. .
0005
05
Appendix D: Experimental Data Plots for Joint 3
160000 -Y-== ~ r! ! T ~ T
0.60 amps
140000 ................
....
.... ...............
....
......... 0.56 am ps
100000
0 0.52 amps
0
.
.... .... .... .... ............
:1..................
0.40 amps
00
60
....
---..
.....
i..... .... ....
........
.... .............
..... .. ps....
0......6...am
5............
Appendix D: Experimental Data Plots for Joint 3
50000 -T..-11
....
;....L... ......
................ 4 % ... . ... . f 0.60 amps ..
40 00...
. .. . ... .......- ....
:...
*40000.
................... .
..... .. . .
....
. ...... 4. ............ ....... ..
0.48 amps
~~~ ~ ......
......
... ~ ... ~ ~ ~ ~ 04
~...........
....................... amps
.~~~...~ .......
~ ~ ~..........
~......
... ~ .......
~ ~......~ ~ 04 amps...
amps
0.28....
0.
0 .... 1...... 2... 4..........
~ ~ ~~~~~tm
.5~ ...
(seconds).... ........
. ... . .6 ..
. .. .... .... .. .. . 8 a p
..........
I... ...... .... ..... .2....56..
Appendix D: Experimental Data Plots for Joint 3
0.8
0L 11,1111"
70.
.. . ......... .
.......
...
0.
-T-a
( U. .- . . .. . . . . ... ... . . .. . . .. ... .. . . . . .. . . .. ..
-CO 0 .1*
0 . ........ .......
0.
CL
0
0... h...
. ....... ... ... . ......
.....-..
.... ...
..... .
.
(U ............ . .
.... I........
..... ..... . .. .
......... ..... . ....
. .
4Co* ..... .
C25
Appenix
D Exprimetal Data Plots for Joint 3
-200
....
... ....
....
............. .... 0.28 am ps
-600
..........
..................
. . ...........
i.......
-800. .... ......
... ..... ....... 0...
48..
amps...
a...
* . ....... ...
-12000.52 amps
4 .........
4... ... .. ....
-1800 ............
:................... ...... ....... .......
4........... 4.4.4.. ..... .....
0.......
1......... 4
..............
.... ....
............
time...........(seconds)
.. ...........
-1600~ ~
........
. .... 0 ....... ...
S...;...O...L....L....L........L...
........
....
....
....
....
....
....
........
258... ............
Appendix D: Experimental Data Plots for Joint 3
.... ... ... ... ... ... ... ....... ........... ... ... ...
....... . ... ... ... ... ... ... ... .. ..
0........ ........
o..........
0.60 amps
.. .... ... .4
0 1...
.... 2..
... ...... ....... ......
...... 4......
.....tim (seconds)..
.... ........
...........
0 0.0.0.0.
2259
Appvendix D: Experimental Data Plots for Joint 3
0 ii1t r r.ri-
9.~ ~.4.....I
4warw wn'i*r*4
&.. . . ..........-..
.. ........ - -... .......... ....
............. ....
. 4. . ................. ...
.......... ......... --- 4... 4... 4 4
3£
0C7 .. .. ...t...±......j... ....... .... ...i................. ............. .... I.................
S4..
. ... 4.. .
I.......
E
.-........
L..... .....
... ..... ... ..... .4
............ $
........ .... ....... ..
-70
708K
0-.
. .............. .... ... ... ............ ... ... ........................
6cc
-7038 ?
E.....................4..4.4
708
-260
Appendix D: Experimental Data Plots for Joint 3
'U.
......
-U---------
.. . ... .. . ................
0. ... ..........;............. ........ ...
E.. ..... 1........
........ ....... ...
..
.. .. ... ... ...
. ............ ..
... . . ... . . .T. . .
0........ ....
T.... ....---
0. .. . ... .... ........ ... .. . ... ... ........ .... .... ........ ...
........ ........ .... ... ... ... ... ........ ... ... ...
. ..... . ... ........ .... .... ........ ... ... ........ ... ... . .
....
E......... .... . .... ...
0... .....
0MP 7 ... ............
...
...
........
........
.a............. ...
....
........
0 ... ............ ... ....
j.... ........ .. ....
.............6............
.... ....
........ ....
.... ...
.... 2....
F~ ~~tm (seconds).....F'7--T-7- ...
2 61.
............
........
............. ............. ....
...
Appendix D: Experimental Data Plots for Joint 3
0.15 .1r r1 r 1
S . . .......... ......
-~~.. 0.......4
. 4
S-0.15.. .....
0
~ 0.15 ....... .. . - . ..
i
-0 .25
0 12 34
time (seconds)
262
This appendix contains a complete collection of all of the dynamic-response plots
generated for each testing station using the harmonic-drive model with friction, compliance,
and kinematic error. This model is derived and these results are discussed in section 3.6 of
this document. The format of these illustrations is identical to the experimental time-
response results presented in appendix D so that direct comparisons can be made easily.
Using the numerical simulation described in section 3. 1, values for the position, velocity,
and torque in the model were calculated, and from the input and output position data, the
dynamic position-error across the transmission was derived. A summary of all of the plots
shown in this appendix and their corresponding plot numbers is given in table E. 1. For
each of the plots listed in this table except the dynamic position-error figure, time-response
curves are plotted for several different motor-current step-commands. Depending on the
visual appearance of the data, the different response curves on each of these plots are either
plotted on a single axis or plotted on separate axes. For the sake of clarity, the output-
velocity response data for each harmonic drive is plotted in both of these ways. In the case
of the dynamic position-error plots, a single time-response curve is shown, rather than
several, from which generalizations about the behavior of dynamic position-error in all
operating ranges are made.
263
Appendix E: Model 4 Simulated Plots
EA .3 E.2.3 E.3.3
264
Appendix E: Model 4 Simulated Plots for Joint I
...................................
...... ...
...
80004 4. amps
. . . ...
0~.... ......
..........
. . . . . . . . . . . .
. . .................. 4.. a p
42.8 amps
200003.4 amps
0 ....
40.... ....... . .265......
....................... ...
Appendix E: Model 4 Simulated Plots for Joint 1
25000 ! 1 1 1
. .... ..
..... ....
.......... .......
.... ...........
.......... 6.0 am ps
K 5.0 amps
20000 -.... .
4.4 amps
0
*15000 -
*...
... . .. 4.0 amps
3.6 amps
CL
. ... .................
.. ........
..... 3.2 am ps
hm~mngma~f~h~ifltl~eakm2.8 amps
500.....-
......
.. ......
.........................
IAM I M A.. .......... ........ 2.4.am p
2.0 amps
5 00 * .....
4~--
J... ........ ........
0....i.... ...
266
Appendix E: Model 4 Simulated Plots for Joint I
6.0 amps
5.0 amps
E
* 4 ~4.0 amps
. .- 3.6 amps
0 3.2 amps
E
S- 2.4 amps
2 ~ 2.0 amps
2- 1.6 amps
o 12 34
time (seconds)
Firur
E..3 Jon smlie tprson 0nu-urn
267
Appendix E: Model 4 Simulated Plots for Joint I
..... .-..
. . . . . . .
. .
.oo..).
~~~~~~~....................
. . .
...... ..............
............................. ... 4
~~~~~. . . ....
.......................... .....
S4 ........
4.... ...
S4.4 amps
S .. ............... . .
30:.
. . .
I. . .
...
. .
....
........
.
~~~~~~.
4amps
.........
...... amps........
36...
2268
A ppendix E: M ode! 4 Simulated Plots for Joint 1
6.0 amp
5.0 amps
0 4.4 amps
140
010
4.0 amps
40
20
1.6 amps
200
12 3 4
time (seconds)
269
Appendix E: Model 4 Simulated Plots for Joint 1
..................
..........................
U!
t..L. . .......... £...
. ....
....... . .
.... .... .... 1 . ... ........
158i1
E........
1-.NO. .4 .
.............
Ea
A(PP .........
... ...... .. .. ....:................. ......... ........ ... ...
15
V -....... .....
. . . 4..
...
................... .... .... 4........ .
E
cc
15 ££ j i Lr- -
a. .......... .......
................. .......................
E. .......................
EU ... ..........
0c
NY
0o 1ls I3 4
15e(ecns
... ........................
.......... ... ... ......... ... .... .... 6 ....
T
80 ~ 1 r f 1
.............. ....... ........ ... ...... ......
E -- "-- I
-- -I' - . 46LM&A I I
0d
v 0.
CL
E0i".-ý~it"A
Oa~~~. I 4 . ..................
cc0.
-8
co~~tm p (seconds)1111
Go
.... ........ * SB
MINTA*S
C~i . . . . . . .271
Appendix E: Model 4 Simulated Plots for Joint I
0.2
...... ......... ............. ..............
.~.. ........ ... .... ................
Z
-0.05
0 ..... .
272
Appendix E: Model 4 Simulated Plots for Joint 2
4.5amps
160000
.. .. .........
140000
820000 /..3..0
CO
400
2L
0000 . 18 amps.
S 052.66amps
273..
...
...............
............
...
Appendix E: Model 4 Simulated Plots for Joint 2
4.5 amps
..............
........
.... ............... .... 3.8 am ps
.......
40000
.. ~ I .. . .
.... ...... . 4..
... ....... ... ............... ...
~) 300003.4 amps
20000
. .... . . .t..- . .
................. ..
.........-..... .
..... ...
4 .... .....
:---
4 -*- *-4 - ' ~.-..-. 3.8 amps
4.5 amps
.................. .. ....
-- 3.4 amps
0
E 2.2 amps
1.8 amps
. .. .. . .4. .. .. $ 4....
1.4 amps
........ ...
.I.4................
4..... .........
012 3 4
time (seconds)
E23 eilr
Jon2 nu-m-n iuie tpesposs
275
Appendix E: Model 4 Simulated Plots for Joint 2
45amp
1000
.oo........ ..............
4....---- .....4....... ... ..4. .........
4... *.. . .
.......
4............ 4.... . . ... .... ........ -...
4 .I ...
I...... /.....
~.
.
.... I .
........ 4. 4....
...... ... ... ... .. . ...... . .... ......... -..
..................
.. .............. ........... -
500 /X 30ap
............
4.... ... . .. .
.... .......... 4......... *.. .
..........
II..........
..... 30 a p
300....................
. . . . .276.
Appendix E: Model 4 Simulated Plots for Joint 2
3500
.4 4.... *.*.............
-- 4.5 amps
30
200~~.. 4 ...........
........ ...........
.0............. .......
2 0 .......
:........
aD....... .......
................ .... .. ...
100amp
4 ,49 4 ..
..... T ..............-
4........
4......
.............. ....... I........
358 L
Zt2I~'y
1 i J . . ......
... .....
L . .....
L..... ..............
a E.
. . .
4*****
......... ...........................
~.
. .........
.......
4
..... I.
... . I ... I...
....
.............
4
....
..............
..........
-...........
358.......T........j.....TT.
T... I....j..
.. .... t
..............
.. .....
.... t L I T.T.T.T.
......................................
0 . ......... ...
aa ..........L.. . .1 . . . .
.... ...... 1---1--
-- - .... ........
o.. a. ...........
......
3278
Appendix E: Model 4 Simulated Plots for Joint 2
20 ...
...............
.......
.......
....... .......
!.......
........
........ .......
.........
.......
....... .............
........ .......
.......-------
. . . . . ........
.
. . . . . . . ........................ .
........
........
aE .......
.......
cc 0
........ . . . MIIT.T.WT
...
.............
t .......
.......
....... .......
.
................. .......
.......
. ... ........ L.......
........
........ ......... ................
....... .......
....... ........
................
. ......
------------- ......
........................
.........
.......
........
................
........
........
.......
.......
........
.......
................ ........
.......
................... ........
....... .......
....... .......
....... .......
........
. ......
..........
CL
......
a
0
. ....... :.......
...... .......
..................
........ ......
........ ....... ..........
...............
..... ....... ..................
....... .......
.....
.....
........ .........................
........
.
........ .......
....... ........
........ .......
........ ....... .......
........
................
.....
................ .......
.. ........ ....... .........
.......
................
....... ................
.......
.......
.......
........
......... 1 .......
....... ................
........
........
........
....... ........
........
........ ........
........
....... ........
........
.................
.......
.......
........
........
.......
.......
.......
....
........
.........
E ...
Mkb.:ý " .""
cd 0
.....
.......
cr) .......
....... .......
.......
........
.......
.......
......................... .......
....... ......
........
................
:f --...... .......
........ .......
....... .......
........................
.......
......................... ......
.......
.................
--------
-- ....... 7r......
.........
..................................
........
..............
_E .......
.......
.................. ........ .........................
........
........
.......
........
........
........................ .......
.......
........................
v.........................
?L 7.......
...................
....
.......
.........
.....
................. ...............
.......
........ ........
.......
........
........
........
.. ..............
CO) .. .......
.. ......
......... . . ...
. ..
..... 4 .......
4 .......
............... 4 .......
4 ....... 4z.......
t.................
*.......
4 .......
................
............... .......
................
.......
4
.......
............. 6.......
........ t....... ........
....... .......
........ .......
.......
CL T .. ....
.
........... . -I---i .......
......... t .......
I ....... ........
........
0
I ......................
....
............
-*
....
.......... ......... I .......
.......
....................... ....... .. .......
....... ........................
................
r........
...................
.. . ........ ................
...... ..
....
........
.....
. ................
E 0
rMUMMM
cm ........
....... .......
........
............... ........
........
.......
.......
........ ..........
.. ........
.......
.......
I ...
................
....... ....... ......
; ..... ................
.......
........
.................
.......
.....
. . ................
. . ....
. .. .......
........
.......
.......
......
Al . .......
. ....... .......
.......
0
C4 . .......... ....... ........ I .......
......
..... . .........
.......
. ........ .......
................
.............
.... ......... .......
........ ........
............... .......
....... .......
................ . ......
... . ........
. .......
..................... ...............
.......
....... ....... ......
.......
.......
......................
........................
. ....... . ..
. ....... ..........
.........
....... ...
................ .......
............. ................
....... .......
....... ........
................
.......
.......
................
.......
....... .......
......................... .......
............... ........
....... .......
................
.......
.......
.......
0
L........
.............. ........ ........
....... ........ ; .......
.........................
.......
................ .......
. .......
........
....... .......
....... ...
................ .......
.......
.......
................. : .......
.........
....... . .......
........
........
........
.............
....... ........
.....
......... ......
.......
.......
-20 L E I
-- ----- ........................
i i i __L I
0 2 3 4
time (sewnds)
279
Appendix E: Model 4 Simulated Plots for Joint 2
0.3 ! T
0 12 3 4
time (seconds)
280
Appendix E: Model 4 Simulated Plots for Joint 3
140000-
0.560 amps
0.48 amps
4000
281........
Appendix E: Model 4 Simulated Plots for Joint 3
50000T T7!] TT 1
0.60 amps
.4.....4.
....... .. . .. . . . . .. ... . . . .. . . . . . ... . .
...... '....
... ... ... ... ... ... ... ....... ... ....... ......
400 . .
..... ......... ................... ........... amps
056..
30000-
t ... ........
0.6 ap
20000.....:
.... .- : : 044; amps:
.
.............. 0 0 a p
...........
............................. 0.28.amps
........
....
....
...
28...
2....
...
...
... ...
.......
Appendix E: Model 4 Simulated Plots for Joint 3
0.7 T 1
0.6
0.60 amps
0.56 amps
-0.44 amps
.a . .
E 0.4 - -------- -b-----0.40 amps
* -P't-P
9Y-1 Y-1--¶-r.t-r.-,--t-0.36 amrps
0
E 0.3- ..... .......
...... ----- .......... ......................
&L 0.28 amps
0.1
0 12 3 4
time (seconds)
E13 0I;Lt
Jon0 1P tC rei 111ia~ tpIS OI
283
Appendix E: Model 4 Simulated Plots for Joint 3
0 . r1 7- T1f1
..................
..... 0.28 am ps-
0.40amp
....
.... ........... 0.44s
. *........ amps
......
o...................... ....... ...... ...... . ...... ..................... 8ap
0.52 amps
-18000
4 . I.....I ....
005 1a2p3s
-1284
Appendix E: -Model 4 Simulated Plots for Joint 3
.....
......... ....... ....... .....
........ ....... ...........
...... ... .......... .. .... ....... ....... .....
........ ....... ... ... 28 a....
ps...
................. ... .. ....... ...... .... ...... . . ........ ....... ... ................. .... ...
-50 .....
.....
. ......... . .
.. . .. . 4.......
..-. ....
..............
....
. ......... . 8 a p
g. .. ...... ....... .......... ......... ....... ....... ....... ........ ......... .. ........ ....... ........ ..... ......... .
o0..2 ..................... ...... ...... ....... ....... ....... ....... ....... ...... ...... ........ ...... ...... ....... ......
....-.......- .... -.....-.. ..-....-..... ....... :....... ........ ....... _.......... ..... ...... ..... .. .
.... •
.. ... ......
... ........•,.,..•.•......
.... •...... .....•. ......... ,....•....
. . ..4...... . ...... . . . .. . . ..... . ...... ... ... , ..... • .... ... ......
. . . . .. . . .. . . .. . ... . . .. . . .. . . .
0 .3 a m ps.•...
.......
.
0~~~~~
0 -250
......
.
""".....
..........
..
...
...
... . .....
.. ... ~
"......•.......3 ...... .3......"....... '.......3 .......' .......' ...... ....
.......
"."' ......
:.......
:...... .3...... '. ......-.
... ' ... 3... ... •...
-45 ..............................
............................... ... 56.a..
....
,.....
....
,.....
' . . . ' ...... ........ ....... ....... ....... ....
...
...
...
...
....... ...
.. ...... ...... ...... ....... ....... .'
...... ...... ....... ....... ........ ...... .i....... ....... !....... !...
"I'
' ' .... ...... • ....... '3 .... ...... • ... .. "....... ...... . . .. ...... .•...... .......• .......÷ ....... ....... . ....... ....... ....... ..
-5 0 ."..
.... ....... ..... ...........
............. &.........
.'"..".. ......
.................
..... ...... .... ..
S•.~~~ ~~.....
• .......
. .. ....... ......
.. T...............
....... . ...... ......
T ..............
...... . .-
S.
. .....
. ..; . - . . . . .... . ...
S...•.... ... ... ----..
.... •.... ,... . ... •...... . ....... .. ....... ÷
.... .... " .. - '. ----'-
..
................. 0.52 am ps
...... ....... •...... ....... ..,..... ...... •....... . ...... .,,...... ...... .,
.......• ...... . .......$...... .•...... . ...... .,...... . ...... . ....... ..... .
0.... .....
•.... ...3......'3 ... .... .. ...... ;........ ...... .3......'...... ;...... '3...... ...... -....... "...... ;......;.... ... 3...
S...
....
•.....•.... • .... •-... -. ... ..... .., ...... •....... •....... •....... •....... .i
...... ....... ."•......•. ...... ."•...... .;,....... ......
.... .. ........... ........... .... ..... .... •... .
4650
..
....
..
. ..
...
.. ...4..
..
... . ... ... .... 6 a p
.(...se.c
............... o...n...d........s...
)....
....
.......
... ...
...
.........
... ....
L....
.... ...
....L...
L........
L....
L....
......
..
L.....
..
Appendix E. Model 4 Simulated Plots for Joint 3
r -b . I................ .... .
70o
.75. . . . . . ...
. ...
.. . ........... ....... ... ....... ....
:... ... .. ...
... ......................
............ ... ~..............
............-............. .............
1
....
..... . i .... ... ........
. ... ... . * ..... ....
.... . . .... ...
L.... ....
a 9
E........... .. ...
. TIa - - ..... .....
070
L........ a.........
...
.........
...
..........
........
..................
-7 J
0........... . L....
............... ...... . L......
............................
708
4 . ... ....... . .... ... ... ........ .... ... ... .... ........
&........ -...
L
9.
....
.......
........... . .
......... .. . 9......
E7
08
-....
.............
> r....r............. . ...r....r.............. ... .
...........
L.......286..
Appendix E: Model 4 Simulated Plots for Joint 3
. . ..... 4 ...........
4 4.4..........
E .................. 4 .
I I.....
... .... 4 .... .... .. ..... ........ ... ... .......
.-...
0.
. I
.4.4........... .4...
.4........ 4 . 4
O 4 ~... 4...... . . . .
t..... ....
................ ......
2 E . 141.....
V ........
..ln... . ......
.........
0......-........
F.......r -
.........;.... ............................;...................:-
.
....
0 ... .......... 1......... ................................3 ..................... 4
............ time (seconds)
............. ........ .............. .... ............
E ......... .... ......... ... ..................... ............
S
Cd2-W
Appendix E: Model 4 Simulated Plots for Joint 3
0~0.2-. T-+-r
....- .4... .
i ........
o0.1
0 . . ....
4.... ...... .. . ... *....... . . . . . .
time (seconds)
288
This appendix contains a complete collection of all of the dynamic-response plots
generated for each testing station using the harmonic-drive model including gear-tooth
geometry. This model is derived and these results are discussed in section 3.7 of this
document. The format of these illustrations is identical to the experimental time-response
results presented in appendix D so that direct comparisons can be made easily. Using the
numerical simulation described in section 3.1, values for the position, velocity, and torque
in the model were calculated, and from the input and output position data, the dynamic
position-error across the transmission was derived. A summary of all of the plots shown
in this appendix and their corresponding plot numbers is given in table F. 1. For each of
the plots listed in this table except the dynamic position-error figure, time-response curves
are plotted for several different motor-current step-commands. Depending on the visual
appearance of the data, the different response curves on each of these plots are either
plotted on a single axis or plotted on separate axes. For the sake of clarity, the output-
velocity response data for each harmonic drive is plotted in both of these ways. In the case
of the dynamic position-error plots, a single time-response curve is shown, rather than
several, from which generalizations about the behavior of dynamic position-error in all
operating ranges are made.
289
Appendix F: Model 5 Simulated Plots
290
Appendix F: Model 5 Simulated Plots for Joint I
6.0 amps
......... ........................
..............................
80000 -7..
60 0 :
......... / ...
....
.... L.....
....
....
...................
000......... ......... ... .......
... amps.
40..
................
....
..
..
...
........ ...
....... ...
......
...
...
.~~~~
. . 6 amps.
20000.................. ...............
.........
......
291. ... ........ 44 a p
App~endix F: Model 5 Simulated Plots for Joint 1
25000 TI T
~~~
... ............ ... ... .......
........ amps
6.0.....
5.0 amps
...........
150004.4 amps
CID
4.0 amps
3.6 amps
OIMMUM
R... mMuIww 32. amps
2.4 amps
. . . .. amp
~ * ~' -~-~' 2.6 amps
0 .~ . . . . .~ .. H .0 a p
0 12 3 4
time (seconds)
292
Appendix F: Model S Simulated Plots for Joint 1
- - 5.04amps
Cu 4 ~6.0 amps
4.4 amps
o ~ .~.-3.2 amps
E
~ 2.8 amps
2.4 amps
2 2.0 amps
o0 2 34
time (seconds)
293
Appendix F: Model 5 Simulated Plots for Joint I
I I
4......... 4. 4 ... ............... ... ... .
4...
40 ...................
......... ....
.... I
.....4.............................
............
............. ...... ...
- ... . . . . . . . . .
4.0 amps
..
...................... ...... 4..........
.............. 28....amps.
100
~ ....-
............... am ps.....
-... 6...
........
........ 29 4 ....
.... .................. ...
A vpendix F: Model S Simulated Plots for Joint 1
160 T1 i1 -- 1 T T r 1
t9~ft$fffff~ffftff 6. amp
6.0 amps
140
120 ...
0)
CD
80
4.4 amps
60 ........ 4. amps....
.............
40
40
295
Appendix F: Model S Simulated Plots for Joint I
150
. .......... ....4....... 4
158mli
01 5
...........-........................
....
...
S4. .
... 4 . 4.... 4. . . 4 . 4......
CL.... 4............
....................... 4..........4
.... ... ...................... ... .........
ccJ
Q .......
1 . ................... i........ ................. ...J...........
..... ..... 1....
158
E ...........
. . . . . . . S..O
C29
Appendix F: Model S Simulated Plots for Joint 1
80
...................................... ....-.......
E0 .4 4 4..
ccJ
0.
-0 0.
0 12 3 4
time (seconds)
298
Appendix F: Model 5 Simulated Plots for Joint 2
4...4.5 amps
. 4... ....
.......
..................
. .....
140000-
.................
........... ..... .... 3.8 am ps
120000........
. ... ... .................
1200
E 0....0........0........-....
0 00 ................. .34..... ... .. .4 aamps....
p
I.....
......... ...
.... amps.........
30..
26 amp
2E2aps
000
0.. 1....... 2..... ............ 4
(seconds)...
....................
,, 0"0*'.* ..I 0-
..... ... ....... S
0 r ,**T
299amp
Appendix F: Model S Simulated Plots for Joint 2
5000 -
4.5 amps
......................
........... ....... ....... ....... ....... ....... ...... ....... ....... ....... 3 amps..
.... ; ............3.8 am p
............ ......
... 4 ...
. ........................... ............
0t
0........ .......
.......... ........................ ....... ..... ....... ....... ....... .......
4'0000 .... 34 amps..........
a... .. .. .. .. ...........
.
.... ...... .... ..... I. . 4. .. .
0.
30000 :.......
....
....- .....
3.4 amps
20000
Appendix F: Model S Simulated Plots for Joint 2
.. "....
.. ... ... . ... ........
. .. 3.8 ams
4.5 amps
-- 3.4 amps
. . .
. . . . ... . . . . ...
i3 ~ -3.0 amps
* ~ 2.6 amps
E . . 2.2 amps
2....................................
1.8 amps
1.4 amps
......
..... ....... ... ....
.....................
..... .....
.......... ...............
I.... 4....
0 12 3 4
time (seconds)
301
Appendix F: Model S Simulated Plots for Joint 2
.....
............ .... 4.5 am ps
4 .... ... .-.. 4.4..... . .......
10
0 ..... .4..... .4......... .....
900 ....
4......... ..............
........
4... 4.
....... ... .
....... ..... I..4....
....
... ... .... . .... ......
g0 -. . . . . ....................... . . . .
..........................................
a / 3.4amps
...........
30 .........
....
.... m..
.....
.... 1. amps.
0. ...... ... ... ........... ......
0... 1.. 2... 4....... ....
~~
300 ~ ~ ~ ~ .......
~ tm ........ .. 26a p
.. ......
..............
(seconds)........
........
................ .... .... ...... .... ... .............. ... S...
.. ...
....
........
L....
..... 2........
.....
..... .30
Appendix F. Model 5 Simulated Plots for Joint 2
05
................. ...
0
. . . . .
. . I I .
.0 ... .............
...
0 ......... .
0...t... ........ ....
350 Ir V
I...... ........
....... .1.u......
4 .. . .. .. . .. .
.........................
................. ....... ....... ............. ........ ....
358 I
* 358
.4 . ... . .. ... . . .... . ...... .. ... .... . .... .. ... .. . .
E ~ .....~
~~......................
~ ~ ..... L L............
0
.................. ................ ..... . .... ... ...
01
358............
13...
....
....
...
...
....
...
. ..
......-.....
~...a....a ......
.-..
...........
.......
-...............
......
...
....
......
...
...l
E
....
..
...
......
.... ......
... .
- J..
............
............. ....
....
...............
__-Appendix F: Model 5 Simulated Plots for Joint 2
20 ......
....
......
................
.......
......
......
......
...... .......
........ ......
.................
............
.... ........ ....... ...
........
........ ......
....
U) : .... .......
........ ........
.......
........
.................
............... ........ ..................................
........
.......
0. . . . . . . . .
E
cc 0-
........................
..... ....... ........
........
.......
.. .... !........
.........
.........
...............
. .......
...............
...... ........
.......
......................... .......
........
................ ...............
........
...... ........
................
................. ........
...................
........
..........................
......................................... ...
. ..
.......
... .......
.........................
........ .......
.................
........ . ........
.......
......
...
....... .......I.................
........ .......
.......
................
...............
...... ..............
........
................
....
...
...
............
... ........ -
........
......
E .....
fla 0
Cq
.7 ...............
...........
CO ....... ................
........ .......
......... ........
........ .......
.................
...
.......
.........................
.......
.......
.................
.......
.................
.................
.
.......
... *........
............................
.......
....... .........
........
.......
.......
.............................
.4 ........
.......
....... ...... ........
........ ....... .................
.......
. .......................
....... . .......
.... ..... .......
....... ........ i .......
........
........ i ....... ........
........ ....... .......
....... ........
E 0
uilaizd&L SWAMI
-
.... ..........
......... .............
.....
.................
.................................. ..
..............
C6
..... ..........................
.......... ..........
......................
...............
.. ........
........
.............
.......
........
........
........
.......
................
.....................
.......
..........
------
...... .............
................
.......
........
................
.................
........
................
........ -
...............
..
2i ........ ....... ................
................
.................
..............
............................. ....
....
......... ...
.... .................. .......
........ ....
........................
......................
Cr E
cc 0
vi .................
. ................
................... ......... I........
...........................
........
........ ... ....... ........
.......
.......
........
.................
.......
.......
....... -
........
......
CL .t::::. ....
...... ........
........
.......
........ .........
........ ....... ......
........
......
...............................
0 ..............
...... A ......
.... ...
........... ........
............ ..............
................
..........................
................
.......
........ .......
.. .......
........
..........................
. .......
. ...
E
cc Ow
CY ................
....... .........
..
........ .................
.......
......
.
.......
........
.................. ....... .......
................
........
..........................
.................................
1.........
...
...
................
..
.........................
E 0 11 WIN
C4 .....
. .......
.................
........
.......
.............
...........
......
............
........
........ ................
....... ...
....
................
......
.......................... .......
........
........
......
.... I......
.. ...........
.................
........
........
....... ?................
................. ..........
...........
........
..
L .......
........
........
..................
....... ; .......
........... ........
....... .......
........
....................
.......
........
........ .......
.................
................ .............
...........
....... ...
......
E
cc 0 .......
.......
.......
...... - - -- ---
...............
........
.................
.......... ........
....................... ........
...............
................. ........
.......
.......
................
................
........
........
.......
.......
........
........
........
.......
.......
.......
...............
....... I...................
.................
........
................ .......
...............
........
........
..............
........
.......
-20
0 2 3 4
time (seconds)
305
Appendix F: Model 5 Simulated Plots for Joint 2
0.2
e 0.05...... ..
0 12 3 4
time (seconds)
306
Appendix F: Model.S Simulated Plots for Joint 31
140000
0.560 amps
.... ....
120000.....
.4. ...... 4......4... ............
0.52 amps
100000 ......
.-.......................
.....
.......................
.... ... .... .... .... ... ... 0.48 am ps
r. ...............
....-.... ...
....................
............. ?......04 ap
8OOOC .... 4..... ..... .... .. 4......
0.84.p
..
....... 4.... ...... ....
0.40 amps
40000...........
0360 amps
20000. ..
. ... ............. 028..
amps.......
... .. .. .....
... ....
.. ..... ...........
030
Appendix F: Model 5 Simulated Plots for Joint 3
50000 f1
0.60 amps
S4.........../
40000 ...
......... .............
. 4 . 4...........
g30000-
tt... 7 . 0.52 amps
0.. . . . . . . ..... .. ...... 4 ........ ...... ... ........ ...............
CD
4. .1. . 4. 4. . . 4
CL 0248 amps
.S 20000308
Appvendix F: Model S Simulated Plots for Joint 3
0.7 T1 1T !rff
0.6 I V
NV 0.60 amps
-0.56 amps
-0.52 amps
0.5- ...............
0..45ap
-0.44 amps
-0.36 amps
3 ........ ..................... ... ... ... ... ...........
0.8.mp
0.1.......................................................
Q...L.....
........ ....... ..... . ....
..... ............
. .... ..........
...
012 3 4
time (seconds)
309
Appendix F: Mode! S Simulated Plots for Joint 3
00
4......... ...............
-2.0 t-.
.......... . .I
0.56 amps
----
1. I ~
L J... ....
A ........ £
0 ....
.............
1.. 2...
.......... ...
...........
... ~tm
....... .......
.(... ec..
o...
nds).............
-40 ..... ... .... ...
..
....
.....
.... ..... ..
3.0..
. .......
....... ...
.... ..
Appendix F: Model 5 Simulated Plots for Joint 3
~0.36 amps
-150.
. ... .....................
...
r .... ...7....-
....
L
L i.......
.......
L.......
.....
L1 .......
.......
.....
.... .. ... ....... ....
...... ....... ........ ..• ..... ........ ....... ÷ ....... ...... .......
....
r . ...7 .............. -r....... :....... ...-.....
.
mp
•., .. ......
•[ ...... ...... ......... ........ ...... .t• ...... .!.=.•.....
....... ... ....... •............. ..... 4....... ......
0
CD .... .. 0.48 amps
5D.....
..
....
.....
....
......
..... ... .. .. .. . ....... ........
(D • ~.......
. .. ...... ..T .... .. ....
...................... ...... ...................
... ... .... ......
i.... .... ...... ...... .....
.... ...... -....... ...... •
...... ...... -....... •
. ..... .. ....
. .....
..
..... ...._.. .....-
......-.....................
. ....... ..... ..........
..... ...... . .. ............ ......
....... ....
..... ..
0 . ... ..... .........
-5 • ......
........ ."...... • ......
......... • ..... •"...... ......• ......• • .............. ....... ....... .......
".......
0--- ...
.......
.,i......
... ,•....
.i......
. ....
.......
, ......
. ...... .......
.....
•.......
•......
.......
.....
i.......
....-
•
...... ....... ......... ...... ".......•....... ..... . .. ....... ....... ....... ............ ....... ......... ...
....... ....... .......n a n a p
S...•..
...... .... -............
i i ......
• . .....
• ......
i . ...... ....$ ...... ...... ...... ....... i ...... ....... ....... •,....... ...
...
....
......
......
31...1....
...
...
Appendix F: Model S Simulated Plots for Joint 3
U)
EL ..... .... .............. .... . .... . .........
CM
* -70
*'WON
4-705 - .. J
di 64
................... r
C, -0
. . __ ._ ........ ........
......
2... 3 ...
........... 4...
time..(seconds)
0 0 0...... 0....... .... .... .....0... .... .. 0...
E .............. 312.
Appendix F: Model S Simulated Plots for Joint 3
.
.4
................ .. .... ... ........ ...
....... ............
.......................*......
i*..F.. T---
E A N .................... ..............
... ... ... ........ ... ... ........
cm......... ...
........ ........
W0........
.... ......... ........ ........ ........ .... ....*.... ......... ... ... ... .......
05 ........ ............ ... ... ............ ... ... ............ ........ ........
E... ......... ............... .................. .................
E
co
313...........
..
0......... ............. ... .......
Appendix F: -Model 5 Simulated Plots for Joint 3
0.4
OD0. 1 2...4
03.4