Development and Design of Constant-Force Mechanisms
Development and Design of Constant-Force Mechanisms
BYU ScholarsArchive
2002-11-08
This Thesis is brought to you for free and open access by BYU ScholarsArchive. It has been accepted for inclusion
in Theses and Dissertations by an authorized administrator of BYU ScholarsArchive. For more information, please
contact [email protected], [email protected].
DEVELOPMENT AND DESIGN OF
CONSTANT-FORCE MECHANISMS
by
Brent L. Weight
Master of Science
December 2001
Copyright © 2001 Brent L. Weight
All Rights Reserved
BRIGHAM YOUNG UNIVERSITY
of a thesis submitted by
Brent L. Weight
This thesis has been read by each member of the following graduate committee and by
majority vote has been found to be satisfactory.
____________________________ _______________________________________
Date Larry L. Howell, Chair
___________________________ _______________________________________
Date Spencer P. Magleby
___________________________ _______________________________________
Date Mark S. Evans
BRIGHAM YOUNG UNIVERSITY
As chair of the candidate’s graduate committee, I have read the thesis of Brent L. Weight
in its final form and have found that (1) its format, citations, and bibliographical style are
consistent and acceptable and fulfill university and department style requirements; (2) its
illustrative materials including figures, tables, and charts are in place; and (3) the final
manuscript is satisfactory to the graduate committee and is ready for submission to the
university library.
___________________________ _______________________________________
Date Larry L. Howell
Chair, Graduate Committee
CONSTANT-FORCE MECHANISMS
Brent L. Weight
Master of Science
begins by reviewing past work done in the area of CFMs and then develops new non-
dimensionalized parameters that are used to simplify the calculations required to design a
CFM. Comparison techniques are then developed that utilize these non-dimensionalized
lengths, normal displacements, and feasible design orientations. These comparison tech-
niques are then combined with optimization to define new mechanisms with improved
performance and range of capabilities. This thesis also outlines a design process, methods
to identify mechanisms that are suitable for a given design problem, and relationships and
trends between variables. The thesis concludes by discussing the adaptation of CFMs for
use in electrical contacts and presenting the results of a design case study which success-
This research was supported by funding from the Utah Center of Excellence Pro-
gram and the National Science Foundation through grant No. DMI-9624574. These foun-
I would like to acknowledge the continuous help and support of Dr. Larry Howell.
He has dedicated countless hours to training, guiding, and teaching me. His impact on me
I would also like to acknowledge the help of the Mechanical Engineering Faculty
and staff. They have been an immense help to me and have always been there. Especially
the help of Dr. Spencer Magleby and Dr. Mark Evans for their support and direction with
The help of Christopher Mattson, Dan Carroll, and Ryan Weight in development
and testing of the constant-force electrical contact and Cameron Boyle for his efforts in
acknowledged.
Most of all, I would like to thank my wonderful wife. She is truly a pillar of
strength in my life. I am thankful for her patience with long, odd hours, for her genuine
interest in my research, her playful teasing, and her efforts to motivate me through to the
end.
TABLE OF CONTENTS
CHAPTER 1 Introduction and Background ........................................................................ 1
1.1 Thesis Problem .................................................................................................... 1
1.2 Background.......................................................................................................... 2
1.2.1 Constant-Force Mechanisms.............................................................. 2
1.2.2 Compliant Mechanisms ..................................................................... 3
1.2.3 The Pseudo-Rigid-Body Model (PRBM) .......................................... 3
1.3 Literature Review ................................................................................................ 4
1.3.1 Constant-Force Systems .................................................................... 4
1.3.2 Constant-Force Mechanisms.............................................................. 5
1.3.3 Compliant Mechanisms ..................................................................... 8
ix
3.2.3 Classification Summary................................................................... 29
3.2.4 Mechanism Inversions ..................................................................... 30
x
5.6.3 Normal Displacement Behavior....................................................... 81
xi
7.7 Optimization Results Summary ....................................................................... 115
xiii
APPENDIX C Optimization Plots.................................................................................... 217
xiv
LIST OF FIGURES
Figure 1.1 (a) A compliant parallel mechanism and (b) its PRBM counter part .......... 4
Figure 1.2 Two element chain of constant-force generator mechanisms
developed by Nathan (1985) .................................................................. 6
Figure 1.3 Rigid-body CFM developed by Jenuwine and Midha (1994) ..................... 7
Figure 2.1 Compliant and rigid-body slider crank model and parameters .................. 12
Figure 2.2 Fifteen original configurations ................................................................... 17
Figure 3.1 Fifteen original configurations ................................................................... 26
Figure 3.2 Refined classification system which identifies flexible
segment configurations ........................................................................ 28
Figure 3.3 Summary of classification scheme ............................................................ 30
Figure 4.1 Values of ∆θ for each pivot point ............................................................. 33
Figure 4.2 Graph of α vs. d for Class 1A-spp-a .......................................................... 40
Figure 5.1 Ψ values for sub-class b ............................................................................. 61
Figure 5.2 Normalized Ψ values for sub-class b ......................................................... 62
Figure 5.3 Normalized Ψ values for sub-class a ......................................................... 63
Figure 5.4 Normalized Ψ values for sub-class b without stress equalization ............. 64
Figure 5.5 Parameters for Ξ'abs calculation ................................................................. 68
Figure 5.6 (a) In-plane and (b) out-of-plane manufacturing orientations ................... 71
Figure 5.7 Graphical representation of flexible segment design area ......................... 76
Figure 5.8 Normal displacement definition ................................................................. 78
Figure 7.1 Explore plot of R versus Ξ'ex and Ψ for (a) Class 1A-lpp-a and
(b) Class 1A-spp-a ............................................................................... 96
Figure 7.2 Explore plot for (a) Class 1A-lpp-b and (b) Class 1A-spp-b ..................... 97
Figure 7.3 Explore plots for (a) Class 1B-psp-a (b) Class 1B-psp-b .......................... 99
Figure 7.4 Explore plots for (a) Class 1B-plp-a (b) Class 1B-plp-b ......................... 100
Figure 7.5 2-D explore plots for slp-a of R and K1 with contours of
(a) Ξ'ex and (b) Ψ ............................................................................... 103
Figure 7.6 2-D explore plots for ssp-a with R vs. K1 and contours of
(a) Ξ' ex and (b) Ψ .............................................................................. 104
Figure 7.7 Optimum plot for lps-a which shows Ξ'ex, Ψ and optimal
K2 for given R values ......................................................................... 109
Figure 7.8 2-D explore plots for sps-a with R vs. K2 and contours of
(a) Ξ'ex and (b) Ψ ............................................................................... 110
Figure 7.9 Optimum plot of design variables K1 and K2 versus R for
sss-b with curves of Ξ'ex and Ψ ......................................................... 113
Figure 8.1 Interdependencies of coupled primary equations .................................... 120
Figure 8.2 Flow chart of CFM design steps .............................................................. 124
xv
Figure 8.3 Definition of flexible and rigid segment lengths ..................................... 126
Figure 8.4 Summary of M value requirements for different types of
identified mechanisms ....................................................................... 129
Figure 8.5 Overlap of guaranteed stress and force feasible mechanisms .................. 133
Figure 8.6 Summary of variable effects on force and stress magnitudes for
an increase in variable magnitudes .................................................... 136
Figure 9.1 (a) Pogo type connector and (b) cantilever type connector ..................... 152
Figure 9.2 Typical compression slider-crank constant-force configuration ............. 155
Figure 9.3 Simulation of pin joints with a circular cam ............................................ 156
Figure 9.4 Limitations and solutions to limitations of simulated pin joint method .. 158
Figure 9.5 Selected constant-force electrical connector configuration ..................... 158
Figure 9.6 Selected CFEC configuration in PDA dock. ........................................... 159
Figure 9.7 Important parameters for general CFEC. ................................................. 160
Figure 9.8 Key Points for the finite element model .................................................. 163
Figure 9.9 CFEC final design .................................................................................... 166
Figure 9.10 CFEC prototype as compared to a dime .................................................. 167
Figure 9.11 (a) General testing setup and (b) close-up of contact with fixture
and probe ............................................................................................ 169
Figure 9.12 Graph of force versus displacement from test data .................................. 169
Figure 9.13 Average and predicted force comparison ................................................ 171
Figure A.1 (a) A small-length flexural pivot and (b) its PRBM ................................ 189
Figure A.2 (a) Cantilever beam with force at end and (b) its PRBM ........................ 191
Figure A.3 (a) Cantilever beam with force at end and (b) its PRBM ........................ 193
Figure A.4 (a) Fixed-guided beam and (b) its PRBM ................................................ 194
Figure B.1 Summary of Matlab file paths and run order ........................................... 198
Figure C.1 2-D Explore plots for slp-b with R vs. K1 and contours of
(a) Ξ'ex and (b) Ψ ............................................................................... 217
Figure C.2 2-D Explore plots for ssp-b with R vs. K1 and contours of
(a) Ξ'ex and (b) Ψ ............................................................................... 218
Figure C.3 Optimum plot for lps-b which shows Ξ'ex, Ψ and optimal
K2 for given R value ........................................................................... 218
Figure C.4 2-D explore plots for sps-b with R vs. K2 and contours of
(a) Ξ'ex and (b) Ψ ............................................................................... 219
Figure C.5 Optimum plot of design variables K1 and K2 versus R for
sss-b with curves of Ξ'ex and Ψ ......................................................... 219
xvi
LIST OF TABLES
Table 2.1 Original results ........................................................................................... 20
Table 4.1 Primary pivot and parameter C for each configuration ............................. 38
Table 4.2 Power function values for all classifications ............................................. 40
Table 4.3 β for each configuration and sub-class ...................................................... 46
Table 4.4 Values for Kx-1 and κx for each sub-class and configuration ..................... 48
Table 5.1 Length parameter formulas and values ...................................................... 57
Table 5.2 Percent constant-force values for the original mechanisms. ...................... 70
Table 5.3 Maximum thickness and minimum width constants .................................. 73
Table 5.4 Summary of manufacturing orientation possibilities ................................. 77
Table 5.5 Summary of dNmax values .......................................................................... 81
Table 6.1 Summary of different possible optimization problems .............................. 90
Table 7.1 Optimization results ................................................................................... 96
Table 7.2 Parameter summary for new mechanisms in Class 1A. ............................. 98
Table 7.3 Optimization results for Class 1B configurations .................................... 101
Table 7.4 Parameter summary for new mechanisms of Class 1B ............................ 102
Table 7.5 Optimization results for Class 2A ............................................................ 105
Table 7.6 Summary of Class 2A parameters ........................................................... 107
Table 7.7 Summary of Class 2A thickness and width ratios ................................... 107
Table 7.8 Optimization results for Class 2B ............................................................ 110
Table 7.9 Summary of Class 2B parameters ............................................................ 112
Table 7.10 Summary of Class 2B thickness and width ratios .................................... 112
Table 7.11 Optimization results for Class 3A ............................................................ 114
Table 7.12 Summary of Class 3A parameters ........................................................... 115
Table 7.13 Summary of Class 3A thickness and width ratios ................................... 115
Table 7.14 Largest stiffness parameters and percent increase ................................... 115
Table 8.1 Needed values to calculate flexible and rigid segment lengths ............... 127
Table 8.2 Summary of stress and force feasibility methods .................................... 130
Table 8.3 Example 1 design requirements ............................................................... 139
Table 8.4 Final design values ................................................................................... 141
Table 8.5 Summary of example 2 design requirements ........................................... 142
Table 8.6 Design problem summary for example 3 ................................................. 146
Table 9.1 Summary of case-study constraints ......................................................... 161
Table 9.2 Parameter summary of final design ......................................................... 166
Table 9.3 Dimensional analysis of CFEC prototypes .............................................. 168
Table 9.4 Parameter summary of prototype 3 .......................................................... 171
Table 9.5 Summary of testing and prediction comparisons ..................................... 172
Table 10.1 Mechanisms with highest stiffness .......................................................... 177
xvii
Table B.1 Mapping of parameter names to names in Matlab and OptdesX code .... 199
Table D.1 Combined mechanism table ..................................................................... 222
Table D.2 Combined mechanism table ..................................................................... 223
Table D.3 Combined mechanism table sorted by M ................................................. 224
Table D.4 Combined mechanism table sorted by M ................................................. 225
Table D.5 Combined mechanism table sorted by Ψ ................................................. 226
Table D.6 Combined mechanism table sorted by Ψ ................................................. 227
xviii
CHAPTER 1 INTRODUCTION AND
BACKGROUND
Many efforts have been made to design systems that produce constant force. Some
of these systems use complex control loops and algorithms that result in costly systems or
operate only in tension. There still exists the need for inexpensive compression constant-
force mechanisms.
Through the effort of others and in correlation with other research projects, com-
pliant constant-force mechanisms have been developed. However, despite the many recent
tions. This includes their design, their limits, and how to improve their performance.
Additionally, there exists new configurations that have yet to be explored. Once a better
force mechanisms. Using the results of prior research efforts as a starting point, this work
1
first seeks to understand compliant constant-force mechanism behavior from both a stress
and force viewpoint. This work then examines different ways of comparing constant-force
mechanisms and then use these comparison methods to define mechanisms with improved
performance. Finally, this work attempts to outline a design method that can be used inde-
pendent of the pseudo-rigid-body model. In parallel with this theoretical work, an attempt
this research. This includes the behavioral model developed for these mechanisms and
other important background information. After the review, the new research of this thesis
1.2 Background
and/or torsional springs or they can be compliant mechanisms. In general, they use the
output force over a large range of input displacement. This is accomplished by determin-
ing specific geometric ratios that allow for equal increases in stored strain energy and
2
mechanical advantage. In this way, the output remains constant throughout the displace-
ment.
deflection of one or more of its members. Compliant mechanisms offer several advantages
allow energy storage directly within a flexible member, eliminating the need for additional
energy storage devices (i.e. springs) found in rigid-body mechanisms. The member deflec-
tion also allows for the replacement of pin joints with small-length flexural pivots, or liv-
ing hinges, thereby reducing part count and assembly time. In fact, one of the most
fewer pieces of material providing savings in both production time and manufacturing
cost. This increase in performance with lower maintenance makes them better suited for
harsh environments. Due to these advantages, compliant mechanisms are replacing many
rigid-body mechanisms.
The pseudo-rigid-body model is used to help model and design compliant mecha-
nisms. With this design technique, compliant mechanisms can be converted into function-
analysis can be performed. Once designed and analyzed, the resulting rigid mechanism
3
(a) (b)
Figure 1.1 (a) A compliant parallel mechanism and (b) its PRBM counter part
senting compliant members as rigid links with torsional springs at the pin joints to account
for the moments at the pin joints due to compliant member deflection. An example of this
can be seen in Figure 1.1. The compliant parallel mechanism in Figure 1.1a has the same
motion and force characteristics as the rigid mechanism in Figure 1.1b. The PRBM works
well for both small and large deflections as well as with a variety of compliant member
types (i.e. small length flexural pivots, living hinges, and fixed guided beams). For more
For many years, there has been a search for reliable mechanical methods that pro-
duce a constant force. The first real success in this endeavor was the development of con-
stant-force tension springs. These springs, also known as “Neg’ator” springs, consist of a
4
coil of flat spring material which has been given a heavy forming operation. When
unstressed, the material tends to form a tight coil. These springs exhibit little change in
load with deflection (Wahl, 1963). Constant-force tension springs have been around for
many years and can be found in many common applications such as inertia reel seat belts,
tape measures, and pull starts (Williman, 1995). They have even been used in creating
Much work has also been given to develop drive units that produce a constant
force. These include electrical, hydraulic, and pneumatic systems (Nathan, 1985). Many
of these drive systems use complex algorithms and feed-back loops to achieve the desired
goal. Bossert et al. (1996) developed a complex algorithm for following unknown surfaces
with a robotic arm. As part of the method, the robot kept a constant normal force on the
surface. In other work (Chang and Fu, 1997), a complex deburring model was used to pro-
duce a drive system that maintained a constant normal force on the workpiece while fol-
lowing a prescribed path. Successful drive systems have been developed and
demonstrated.
Recently, much effort has been made to design mechanisms that produce a con-
stant output force. Nathan (1985) proposed a rigid-link constant-force generator. His work
resulted in the creation of a hinged lever that produces a constant unidirectional force for
any position. This work was extended resulting in a chain of parallel mechanisms that
5
Supported Mass
Figure 1.2 Two element chain of constant-force generator mechanisms developed by Nathan
(1985)
would support a mass when moved to any position. A diagram of this mechanism can be
seen in Figure 1.2. This mechanism can be seen in applications such as desk lamp stands
(Nathan, 1985). Jenuwine and Midha (1994) have proposed a rigid-link CFM. This mech-
anism, as seen in Figure 1.3, uses rigid-links and linear springs to achieve a constant-force
1994, Howell et al., 1994, Midha et al., 1995). Millar et al. (1996) developed non-dimen-
sionalized parameters to facilitate their design and tested several mechanisms. Murphy et
al. (1994) used type synthesis on the compression CFM to develop 28 configurations
while Howell et al. (1994) performed dimensional synthesis of several of these configura-
tions. Most recently, Evans and Howell (1999) implemented the compliant CFM into a
6
r4
r3 r5
r8
springs rigid link
F
r7
r2 r6
r1 spring end
attached
to ground
Figure 1.3 Rigid-body CFM developed by Jenuwine and Midha (1994)
robot end-effector that successfully demonstrated constant-force behavior while cutting
glass.
pliant mechanisms. In one study, they examined a fully compliant constant-force mecha-
gravity equilibrators. These mechanisms are similar to the work done by Nathan (1985),
who developed the constant-force generator commonly found in desk lamps, but have a
larger range of motion. Additionally, Herder and Berg (2000) developed a statically bal-
anced compliant mechanism. This system consisted of a compliant gripper on the output
end and a balancing mechanism on the input end. In this fashion, the force required by the
user to deflect the gripper is offset by the balance mechanism, and the user only feels the
The chapter provides valuable information and a foundation for the proposed research.
The notation developed, equations derived, and the optimization problem used are pre-
Compliant mechanisms get their motion and energy from the deflection of their
members. The PRBM uses links and springs to model motion and compliance (Howell,
2001).
The PRBM allows for easy design and synthesis of compliant mechanisms. Com-
pliant mechanism synthesis can be divided into rigid-body replacement synthesis and syn-
thesis for compliance. Rigid-body replacement synthesis deals only with the motion and
path of the mechanisms, while synthesis for compliance takes into account both the
used. Conventional Newtonian methods require free-body diagrams of each link and result
in forces for the entire mechanism. A second method, the principle of virtual work, also
works well with the PRBM. This method looks at the whole mechanism and accounts eas-
The PRBM also allows for the determination of the degrees of freedom of a com-
pliant mechanism. While traditional methods predict many compliant mechanisms are
structures, consideration must be made for the movement made possible by the compliant
8
sections. Methods have been developed that take this into consideration allowing for accu-
9
10
CHAPTER 2 BEHAVIORAL MODEL
DEVELOPMENT AND
PREVIOUS RESULTS
Work was done prior to this thesis on the development of a behavioral model for
in this work and is used as the basis for further development. All other equations devel-
oped in this thesis are derived from this behavioral model or the PRBM and thus have the
This chapter summarizes the derivation of the behavioral model, the extend of its
development, and the validity of the model. The following is summarized from Howell
slider-crank model using the PRBM and the principal of virtual work. Several non-dimen-
sional parameters can be developed and the model simplified. This next section discusses
11
flexible
segments flexible
(small-length segment
flexural pivots) (small-length
flexural pivot) slider
(rigid segment)
rigid segments F
θ'2
(a)
θ3
k2
r3
r2 slider
deflection
∆x
k3
θ2 F
k1
r1
(b)
Figure 2.1 Compliant and rigid-body slider crank model and parameters
The original behavioral model is based upon a simple compliant slider crank
mechanism. The PRBM and standard kinematic equations are used to solve for the posi-
tion of the slider crank given a deflection. The variables used in the equations and the
mechanism orientation is shown in Figure 2.1b. The known variables for the problem are
2 2 2
r1 + r2 – r3
θ 2 = acos -------------------------- (2.1)
2r 1 r 2
12
– r 2 sin θ 2
θ 3 = asin --------------------- (2.2)
r3
Figure 2.1a shows one configuration of the compliant version of the slider crank.
The appropriate lengths of the flexible segments can be determined using the PRBM as
discussed in Appendix A.
The principle of virtual work and the PRBM can be used to determine the static
force for a given deflection. It can be assumed that all force references throughout this
To determine the static force for a given deflection, equations must be developed
relating displacement, compliant member deflection, and static input force. Using the
principle of virtual work and the PRBM, a fictitious or virtual displacement ( δz ) can be
δW = F ⋅ δz (2.4)
δW = M ⋅ δθ (2.5)
13
where δW is the virtual work due to the moment, M , and virtual angular displacement,
δθ . A good equation for conservative forces is found by taking the derivative of potential
dV
δW = – ------- δq (2.6)
dq
∑ F ⋅ δz + ∑ M ⋅ δθ – ∑ --------
dVk
δW = i i
- δq
j j k (2.7)
dq k
i j
Having established equations for virtual work, the principle of virtual work can be
The net virtual work of all active forces is zero if and only if an
ideal mechanical system is in equilibrium.
This principle allows equation (2.7) to be set equal to zero. If all virtual displacements are
∑ F ⋅ A + ∑ M ⋅ B – ∑ --------
dVk
i j- ( δq ) = 0 k (2.8)
dq k
i j k
where A and B are vectors that change the linear and angular displacements into terms of
ment), then the remaining equation can be solved for the unknown force or moment.
The method of virtual work was applied to the slider crank. The variable θ2 was
chosen to be the generalized coordinate. Equations for the virtual work associated with
14
each torsional spring were developed. A fourth equation was developed relating an
unknown static input force applied to the slider in the horizontal direction. These equa-
tions were summed, and the principle of virtual work was applied. The force, F , was
r 3 cos θ 3 [ k 1 θ2 + k 2 ( 2π + θ 2 – θ 3 ) ] + r 2 cos θ 2 [ k 2 ( 2π + θ 2 – θ 3 ) + k 3 ( 2π – θ 3 ) ]
F = ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- (2.9)
r 2 r 3 ( sin θ 2 cos θ 3 – sin θ3 cos θ 2 )
This equation tells how the force, F , is related to the link lengths, spring constants,
2.1.3 Non-Dimmensionalization
beneficial to generalize the model to simplify its use. One method to do this is to try and
remove all the independent variables replacing them with dimensionless parameters. In
the complimentary work done by Millar et al. (1996), three non-dimensionalized parame-
r
R = ----3 (2.10)
r2
k2
K 1 = ----- (2.11)
k1
K3
2 = ------ (2.12)
k1
15
These parameters where substituted into equation (2.9). Furthermore, the trigonometric
identity sin ( α – β ) = sin α cos β – cos α sin β was used to simplify the denominator of
k
F = ----1- Φ (2.13)
r2
where
( R cos θ 3 [ θ 2 + K 1 ( 2π + θ 2 – θ 3 ) ] + cos θ 2 [ K 1 ( 2π + θ2 – θ 3 ) + K 2 ( 2π – θ 3 ) ] )
Φ = ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
- (2.14)
R sin ( θ 2 – θ 3 )
depends only on the non-dimensional parameter Φ and the spring constant k 1 and link
length r 2 . The spring constant is considered to be the stiffness parameter, while the link
parameters reduces the number of independent variables in the model, making the model
easier to use.
Murphy et al. (1994) performed type synthesis on the slider-crank model. This
work resulted in the development of 28 configurations for the CFM. The 28 configurations
consist of different arrangements of pin joints and flexible segments. These 28 configura-
tions were reduced to 15 viable configurations and are divided into 5 classifications based
16
Class 1A
(a) (b)
(c) (d)
Class 1B
(e)
(f) (g)
Class 2A
(h) (i)
Class 2B
(j) (k)
(l) (m)
Class 3A
(n) (o)
The objective of the original project was to find combinations of the non-dimme-
17
2.3.1 Displacement Vector
r 1 = ( r 2 + r 3 ) – ∆x (2.15)
where
d
0 ≤ ∆x < --------- ( r 2 + r 3 ) (2.16)
100
was created where d is the deflection parameter and the vector ∆x contains 50 points (an
arbitrary number). For example, if the mechanism is 10 inches long and d is 40, then the
The vector r 1 was then used to calculate the angles θ2 and θ3 for the 50 positions
using equations (2.1) and (2.2). Subsequently, these values were then used to calculate the
force from equation (2.13). The result of this process is a force vector, F which corre-
An objective function was needed for the optimization routine. This was accom-
max ( F )
Ξ = ------------------- (2.17)
min ( F )
18
This parameter shows how constant the static force is for the mechanism throughout the
Minimize Ξ (2.18)
The design variables for this problem vary depending upon the configuration. The
design variables are taken from the non-dimensionalized parameters R , K 1 , and K 2 . The
parameter R is always included as a design variable while K 1 and K 2 are added when
there is a k 2 and k 3 respectively. The analysis variables are the displacement vectors, k 1 ,
Each of the five classes of configuration found in Figure 2.2 were run through the
optimization code to find ideal values for the parameters R, K1, and K2. The values for K1
and K2 were set to zero for configurations in which pin joints were present at the respec-
tive joints. The results of the original work are summarized in Table 2.1
This table shows that a set of viable non-dimensionalized parameters were found
for all 5 CFM classifications with deflections (d) of both 16 and 40. These configurations
have a percent constant-force close to one. However, these solutions are not the only set of
19
unique solutions to the problems, but represent what was felt to be the best combinations
The validity of the behavioral model was verified through prototyping and testing
of mechanisms from each of the classes of CFMs. Millar et al. (1996) describe some of
these tests and their results. Mechanisms from different classifications were prototyped
from various materials and tested using a compression testing machine. The mechanisms
displayed a constant output force for a large range of deflection. However, it was noted
that the initial deflection of the mechanism resulted in a large force spike. This was attrib-
uted to internal friction (for partially compliant mechanisms) and polymer re-alignment.
However, after the initial deflection, the results were very reliable.
Class d R K1 K2 Ξ Φ
1A 16 0.8274 - - 1.0030 0.4537
40 0.8853 - - 1.0241 0.4773
1B 16 1.0000 1.0000 - 1.0564 2.0563
40 1.0000 1.0000 - 1.1576 2.1513
2A 16 0.3945 0.1906 - 1.0015 0.9575
40 0.4323 0.2237 - 1.0058 1.0466
2B 16 0.7591 - 0.1208 1.0029 0.5230
40 0.8441 - 0.1208 1.0235 0.5438
3A 16 2.6633 1.0000 12.6704 1.0002 3.4016
40 2.0821 1.0000 9.3816 1.0049 3.6286
20
2.6 Model Exceptions
Close examination of previous work reveals some exceptions and alterations to the
behavioral model. These exceptions are important to the general understanding of past
work and will be re-examined in future work. For this purpose, they are briefly explained
at this point.
2.6.1 Eccentricity
Close examination of the slider-crank in Figure 2.1 on page 12 reveals that when
the model is fully extended, the mechanism is horizontal. This shows that eccentricity, or
an offset from the slider in the vertical direction, is omitted. In fact, eccentricity is omitted
The PRBM of the class 1B mechanism consists of a torsional spring in the center
of the mechanisms with pin joints on either side. This requires that k1 equal 0. However,
infinity. To solve this problem, K 1 is set equal to 1, which allows k 1 to equal k 2 . Since k 3
( R cos θ 3 [ K 1 ( 2π + θ 2 – θ 3 ) ] + cos θ2 [ K 1 ( 2π + θ2 – θ 3 ) + K 2 ( 2π – θ 3 ) ] )
Φ′ = ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- (2.19)
R sin ( θ 2 – θ 3 )
21
It should be noted that Equation (2.19) is used exclusively with the Class 1B mechanisms.
k2
Fin = ----- Φ′ (2.20)
r2
This configuration was then run through the optimization routine and the results can be
and a classification system for CFMs. Although these equations and classification system
work, there are several challenges to using them in the exploration and design of CFMs.
First, the equations developed rely heavily upon the PRBM making the design of
CFMs difficult for engineers with little experience with the PRBM. Additionally, although
the equations are non-dimensionalized, they require that a full model be developed for
each configuration to determine the stresses in the flexible segments and the mechanism
output force.
Second, the classification system developed identifies only large groups of mecha-
nisms and provides no way of identifying smaller groups of mechanisms or specific mech-
anisms.
22
Third, the prior work used only one method for comparing mechanisms, the
Finally, a methodical approach to design was not developed. The basic design
steps and issues are not discussed and the advantages and disadvantages of each configu-
ration of mechanisms are not known, making the selection of an appropriate configuration
for a given application difficult. The following chapters address these issues.
23
24
CHAPTER 3 CLASSIFICATION
REFINEMENT
To improve the performance of CFMs and better understand how to design them, it
becomes necessary to look more closely at specific mechanisms. This will allow for com-
parison and differentiation between the mechanisms, and eventually lead to an ability to
The original classification system, presented in Section 2.2, does not differentiate
between specific, single mechanisms. For this reason, the rest of this chapter is devoted to
the presentation of a refined classification system based upon the original system. With
The original classification system developed by Millar et al., 1996, is first pre-
sented in Section 2.2. The system, seen in Figure 3.1, is based first upon the number of
25
Class 1A
(a) (b)
(c) (d)
Class 1B
(e)
(f) (g)
Class 2A
(h) (i)
Class 2B
(j) (k)
(l) (m)
Class 3A
(n) (o)
flexible members in each mechanism, and then upon the arrangements of those flexible
segments. It can be seen that mechanisms that have one flexible segment located at the
first pivot fall into Class 1A while mechanisms with three flexible segments fall into Class
with each class of mechanism. These parameters are valid for every mechanism in that
class.
Therefore, the original classification refers not only to the physical arrangement of
flexible segments and pin-joints, but also to the specific parameter set.
26
3.2 Classification Refinement
The original classification system is simple and easy to use. Therefore, the new
system simply adds a method to distinguish between mechanisms within a class that have
the same pseudo-rigid-body model (PRBM) but have different constant-force parameters
and/or flexible segment types. This refinement allows for various levels of classification.
Each level corresponds to various sets of parameters and groups of mechanisms. At the
most refined level, a specific mechanism can be identified along with a specific set of
between possible flexible segment configurations within each class. This configuration is
denoted by a string of letters representing the order and type of pivots used. The letter “s”
will be used for small-length flexural pivots, the letter “l” will be used for long flexible
segments, and the letter “p” will be used for pin joints. Figure 3.2 shows the refined classi-
fication system using the original classes. For each mechanism, a flexible segment config-
uration has been added and can be seen under the corresponding mechanism.
When writing the classification, the configuration is added after the class. For
example, a Class 2A-ssp is a mechanism that has two small-length flexural pivots located
at the first and second pivot points, and a pin joint at the third pivot point. A Class 1A-lpp
is a mechanism that has a single fixed-pinned beam at the first pivot, and pin joints at the
27
Class 1A
(spp) (pps')
(lpp) (ppl')
Class 1B
(psp)
(plp) (plp’)
Class 2A
(ssp) (pss')
(slp) (psl')
Class 2B
(sps) (sps')
(lps) (spl’)
Class 3A
(sss) (sss')
Figure 3.2 Refined classification system which identifies flexible segment configurations
last two pivots. However, notice that the specific set of parameters to be used with the
3.2.2 Sub-Classes
parameters within a given class and configuration. Each sub-class within a given class will
28
have the same PRBM, but will lead to different CFMs. These new sub-classes will be
denoted using a string of lower case letters and numbers. The first letter will denote the
maximum percent deflection for which the sub-class was designed. The other letters and
numbers will be explained later in this work as other parameters are developed and new
The original work performed actually resulted in different sub-classes even though
they were not thought of in this manner. These sub-classes were distinguished by the per-
centage of deflection for which they worked. They are commonly referred to as 16% and
40% deflection mechanisms. In the new classification system, these sub-classes will be
distinguished as sub-class “a” for the 16% deflection mechanisms and sub-class “b” for
When specifying a specific sub-class, the identifying string of letters and numbers
is added after the configuration. For example, if a new Class 2B-lps mechanism is defined
with a maximum deflection of 25% and a unique set of parameters, then the new classifi-
class refers to a group of mechanisms that share a common PRBM, each configuration
refers to a mechanism with specific flexible segment types, and each sub-class identifies a
29
A configuration contains all of the following:
3A-sss-b
Level Class Configuration Sub-Class
Common Flexible
PRBM Parameters
Elements Segments
Group of Family of Specific
Applicable Mechanisms Mechanisms Mechanism
Figure 3.3 Summary of classification scheme
In some cases, mechanisms are simply inversions of other mechanisms. The slider
is fixed, while the fixed end is allowed to slide. These mechanisms are easily accounted
for. In Figure 3.2, these inversions are shown next to their counterparts and have the same
flexible segment arrangement followed by a prime. It should be noted that all parameter
values that have been given are valid for both mechanisms. It is only necessary to apply
the nomenclature in the same way in either instance. For example, in a 1A-spp, the link
length r 2 is always associated with the link with the flexible segment and k 1 is always
associated with the flexible segment regardless of which end of the mechanism is
grounded.
30
CHAPTER 4 STRESS AND FORCE
FEASIBILITIES
To overcome the challenges associated with CFM design and analysis, it is desir-
able to develop a method that would allow for the quick and simple determination of the
stress and force feasibility for a particular application. This would greatly reduce the
amount of work required to determine which, if any, of the CFM configurations is viable
for a given application. Additionally, a method is needed that will be simple to use and
require only a limited understanding of the PRBM. Finally, the new method should aid in
This chapter adds to the work presented in Chapter 2. It begins by deriving several
new parameters which can be used to help analyze the stress and force feasibilities. These
derivations and design techniques are based upon the pseudo-rigid-body model and the
behavioral model. The steps to the derivation are outlined, the end parameters are defined
and further developed, and parameter values are summarized for the sub-classes and con-
figurations of the original results as presented in Chapter 2. The derivations and results are
then followed up with examples that show how the derivation works and its usefulness.
31
The parameters and methods established in this chapter will be used later in this
understanding of their strengths and weaknesses. Additionally, values for the parameters
developed in this chapter will be summarized and tabulated for new configurations devel-
To analyze the stresses in the CFMs, it is necessary to first look at the stress in one
of the links of the mechanisms. The critical stress, σ c , in a flexible beam under bending
Mc
σ c = -------- (4.1)
I
where M is the bending moment in the flexible segment, c is the distance from the neutral
plane to the top/bottom plane, and I is the moment of inertia of the cross section of the
flexible segment. Using the pseudo-rigid-body model, the bending moment M is found to
be
M = K∆θ (4.2)
where K is the PRBM spring constant and ∆θ is the actual angle of deflection of the
beam. The values for ∆θ for each pivot in the slider-crank are defined in Figure 4.1.
32
∆θ=θ3ο+θ2-θ3-θ2o
r2 r3
∆x
∆θ=θ2-θ2o ∆θ=θ3ο-θ3 F
r1
Figure 4.1 Values of ∆θ for each pivot point
EI
K = γK θ ------ (4.3)
l
where γ is the PRBM characteristic radius factor, K θ is the stiffness coefficient, E is the
material’s modulus of elasticity, I is the moment of inertia, and l is the length of the flex-
ible segment.
ments, small-length flexural pivots, and fixed-pinned beams. There are some assumptions
r
l = ---i (4.4)
γ
33
where γ is the PRBM characteristic radius factor and r i is an effective link length associ-
ated with the flexible segment. For fixed-pinned beams, r i can be either r 2 or r 3 .
3. K θ is approximated as 2.65
1. The flexural pivot length ( l ) is much smaller than the corresponding PRBM
l = µr i (4.5)
where µ is the ratio of l over r i and r i will be either one of the two PRBM link lengths,
2. Commonly, the value for µ is 0.10. This value will be used unless stated other-
wise.
3. The values for K θ and γ in Equation (4.3) are 1. This is consistent with the
PRBM.
4. The link length r i used for k 2 (middle pivot) is taken to be the average ( r ave )
r2 + r3
r i = r ave = ---------------
- (4.6)
2
34
Equations (4.4) and (4.5) can be generalized to
l = ρr (4.7)
1
where ρ is either --- for fixed-pinned beams or µ for small-length flexural pivots.
γ
Furthermore, the link length or average link length, r i , depends upon the configu-
ration of the mechanism. With this in mind, a new parameter, ζ , can be developed where
r tot
ζ = ------- (4.8)
ri
and
r tot = r 2 + r 3 (4.9)
This new parameter defines the ratio between the total PRBM length r tot of the
CFM and the link length of interest. The values for ζ for the different link lengths encoun-
r2 … ζ = R+1 (4.10)
1
r3 … ζ = --- + 1 (4.11)
R
r ave … ζ = 2 (4.12)
Rearranging Equation (4.8) and substituting Equations (4.2), (4.3), (4.7), and (4.8)
γζK θ E∆θc
σ c = --------------------------- (4.13)
ρr tot
35
Equation (4.13) gives the stress in a flexible beam according to the PRBM. The
values for γ , K θ , and ρ depend upon the assumptions for each type of flexible beam used,
ζ depends on the configuration, E depends on the material selected, ∆θ is based upon the
deflection and the sub-class as defined in Figure 4.1, and c and r tot depend on the geom-
The stress can be related to the safety factor, SF, and the yield strength, S y , as
σ c ⋅ SF = S y (4.14)
γζK θ ∆θc S 1
- = ----y- ⋅ ------
---------------------- - (4.15)
ρr tot E SF
γζK θ ∆θ
α = -------------------
- (4.16)
ρ
c
A = ------- (4.17)
r tot
Sy
Ω = ----- (4.18)
E
Ω
αA = ------- (4.19)
SF
36
The parameter α is determined by the configuration and sub-class, A is based
upon the geometry, Ω is determined by the material, and SF is a design parameter. The
following section will develop and refine these parameters so that they are easy to use.
sure of the stress in the specified flexible segment at a given deflection. It is dependent
upon the type of flexible segment, the amount of deflection, and the constant-force param-
eter R. The actual size and material of the flexible segment have no affect on this parame-
ter.
Using Equation (4.16) and the assumptions for each flexible segment as stated in
α = 10ζ∆θ (4.20)
α = 1.91ζ∆θ (4.21)
Equations (4.20) and (4.21) show that small-length flexible pivots have an α
approximately 5 times larger than fixed-pinned beams if ζ∆θ is held constant. According
to the relationship defined in equation (4.19), small-length flexural pivots are much higher
in stress than long flexible segments. This is consistent with what would be expected.
The change in angular deflection, ∆θ , in equation (4.16) does not depend upon the
flexible segment type. It depends on the mechanism displacement and the parameter R .
37
The second parameter, A , is a geometric parameter. This parameter depends upon
the flexible segment that has the highest stress, or the primary pivot. Table 4.1 shows the
primary pivot for each configuration. This table holds true provided the inequality
c o ≤ Cc p (4.22)
remains true where c p is the value for c for the primary pivot, c o is the value(s) for the
other flexible segment in the mechanism, and C is a parameter that is mechanism depen-
dent. In the case of the fully compliant mechanism (sss configuration), two C values are
given. The first one is for the second pivot, the second value is for the third pivot.
has the highest stress and therefore, the c value for this flexible segment is c p . Addition-
Primary C
Configuration
Pivot sub-class a sub-class b
Class 1A-spp 1 - -
Class 1A-lpp 1 - -
Class 1B-psp 2 - -
Class 1B-pl p 2 - -
Class 2A-ssp 2 5.090 4.647
Class 2A-slp 2 1.721 1.474
Class 2B-sps 3 1.739 1.405
Class 2B-lps 3 9.081 7.339
Class 3A-sss 1 1.384 8.258 1.043 4.359
38
ally, the c value for the first flexible pivot is c o . If c p is 0.2, then c o must satisfy Equation
(4.22). Therefore,
or
Equation (4.24) indicates that if c o becomes larger than this value, then the flexible seg-
ment with the highest stress changes and a new value for α must be calculated. Values for
C for each configuration and sub-class can also be found in Table 4.1.
Once the primary pivot has been identified, α can be calculated for each classifi-
cation for a percent deflection, d. Figure 4.2 shows a graph of α vs. d for the 1A-spp-a
mechanism. Additionally, the curve has been fitted with a power function. This allows α
α = Md n (4.25)
This procedure was repeated for all of the configurations in both sub-classes and
similar functions were determined. The values for the parameters M and n in the α power
function for each configuration and sub-class are listed in Table 4.2. This table also
restates the information from Table 4.1. From this table, a value for α can be quickly cal-
culated for each classification at any displacement. It should be remembered that the val-
39
α vs d
10
8
Stress Parameter,
CFM 1A-spp-a
0
0 10 20
Percent Deflection, d
ues for α are good for all deflection percentages, d , up to the maximum percent
deflection of the sub-class. The α value, along with other known parameters, can then be
used with Equations (4.17) to (4.19) to validate the stress feasibility of a design given cer-
40
tain parameters. The next section demonstrates the usefulness and practicality of this
method.
Suppose that a CFM is needed in an application where the overall PRBM length
can not be larger than 5 inches and a deflection of 0.6 inches is needed. The mechanism is
to be made with a single rectangular cross section ( c = h ⁄ 2 ) under the constraints that
h ≥ 0.001 inches . The mechanism must use the Class 2B-lps-a configuration, and must be
Assume SF = 1.5 . Using the information in Table 4.2, determine if the mechanism is
This problem calls for the use of the Class 2B-lps-a mechanism. This mechanism
consists of a long flexible beam at the first pivot, a pin joint at the second, and a small-
length flexural pivot at the third. Table 4.2 indicates that the third pivot has the highest
stress. To solve this problem, it is necessary to determine the value of the parameter α .
First, the percent deflection or d of the problem is determined. This is done by dividing the
0.6
d = ------- ( 100 ) = 12 (4.26)
5
41
Therefore, the percent deflection or d for this problem is 12. From Table 4.2, M
M = 3.743 (4.27)
n = 0.511 (4.28)
Using the value of c = h ⁄ 2 for a rectangular cross section, the maximum thick-
ness, h max , can be calculated by combining and rearranging Equations (4.17) and (4.19)
2Ωr tot
h max < ---------------- (4.30)
αSF
For 1010 Steel, Equation (4.30) becomes
2 ( 0.00087 )5
h max < ------------------------------- (4.31)
13.31 ( 1.5 )
or
2 ( 0.0092 )5
hmax < ---------------------------- (4.33)
13.31 ( 1.5 )
or
42
Equation (4.32) shows that the thickness for the flexible segment at the third pivot
must be less than 0.00043 inches when using 1010 Steel to keep the stress below the max-
imum level. This thickness is below the minimum thickness value defined in the problem,
and therefore, 1010 Steel can not be used for this situation. However, Equation (4.34)
shows that using Beryllium Copper up to a thickness of 0.0046 inches will satisfy the
stress requirements for the problem. Therefore, the flexible segment can have any thick-
ness between 0.001 inches and 0.0046 inches, and still satisfy the requirements.
Also, using Equation (4.22), it is possible to determine the maximum thickness for
the flexible segment on the first pivot in the configuration. From Table 4.2, the value for
C for this configuration is 9.08 . Using the value for C , the equation c = h ⁄ 2 , and the
This indicates that the width of the first flexible pivot must be less than 0.0418 when the
thickness 0.0046 inches is used for the flexible segment of the third pivot. Any value
above 0.0418 will cause the stress parameter equation to become invalid as the flexible
segment with the maximum stress shifts from one segment to the other. At this point, new
The force feasibility equations are similar in purpose to the stress feasibility equa-
tions developed above. These equations contain a set of unique parameters that help in
43
determining if a given mechanism can meet the force demands of a given situation while
still satisfying all the constraints. In many cases, the information used and determined in
the stress feasibility calculations can be applied to the force feasibility calculations.
The derivation for the force feasibility will be presented, results for the original
configurations will be shown, and the example started above will be continued.
The static force equation for the CFMs is given in equation (2.13) as
k1
F = ----- Φ (4.36)
r2
where k 1 is the PRBM spring constant for the first pivot, r 2 is a PRBM link length, and
EI1
k 1 = γK θ -------- (4.37)
l
where the parameters are the same as defined in the above sections.
Following the assumptions explained in Section 4.1.1 for l , and substituting equa-
EI1
k 1 = γK θ -------- (4.38)
ρr 2
44
In turn, this is combined with equation (4.36) to give
γK θ EI 1 Φ
F = ----------------------
- (4.39)
ρr 22
Using the parameter ζ presented in Equation (4.8), Equation (4.39) can be further
generalized to
γζ 2 K θ EI 1 Φ
F = ----------------------------
- (4.40)
ρr tot2
where
ζ = R+1 (4.41)
γK θ EI 1 Φ ( R + 1 ) 2
F = -------------------------------------------
- (4.42)
ρr tot 2
sub-class, E depends upon the material, and r tot and I 1 depend upon the geometry. The
variables that are dependent upon configuration and sub-class can be combined to form a
γK θ ( R + 1 ) 2 Φ
β = ----------------------------------
- (4.43)
ρ
45
and
βEI 1
F = -----------
2
- (4.44)
r tot
The first moment of inertia, I1 , is associated with the flexible segment of the first
pivot ( k 1 ). The parameter β is easily calculated for each specific configuration (with an
exception for Class 1B, discussed in Section 4.2.4). The results of these calculations can
material, length, and cross sectional geometry are suitable to achieve a desired force. This
equation can be readily used without a complex model to run a simple feasibility check or
In order to use the force equation with the stress equation, it becomes important to
relate the moment of inertia of the first flexible segment, I1 , with the moment of inertia of
the flexible segment with the maximum stress. This moment of inertia will be denoted as
sub-class a sub-class b
Configuration
Φ β Φ β
Class 1A-spp 0.4537 15.1508 0.4773 16.9649
Class 1A-lpp 0.4537 2.9008 0.4773 3.2482
Class 1B-psp 2.0563 82.2520 2.1500 86.0000
Class 1B-pl p 2.0563 15.7482 2.1500 16.4658
Class 2A-ssp 0.9575 18.6332 1.0466 21.4708
Class 2A-slp 0.9575 18.6332 1.0466 21.4708
Class 2B-sps 1.2259 37.9347 1.2154 41.3322
Class 2B-lps 1.2259 7.2631 1.2154 7.9136
Class 3A-sss 3.4016 456.4868 3.6286 344.6931
46
I x where x is the primary pivot as defined in Table 4.1. To relate the moments of inertia,
γK θ EI
k = --------------- (4.45)
l
will be used. Rearranging this equation and substituting values according to the methods
γζK θ EI
k = ------------------- (4.46)
ρr tot
k2
K 1 = ----- (4.47)
k1
k3
K 2 = ----- (4.48)
k1
can now be used to relate the spring constants for the first flexible pivot and any other
flexible pivot. If i is used to represent any of the flexible segments, then the spring con-
stant of any flexible segment can be related to I1 through Equation (4.47) and (4.48). Gen-
Ki – 1 γ 1 ζ1 Kθ E 1 I 1 γi ζi Kθ E i I i
------------------------------------------- = -------------------------
1 i
(4.49)
ρ 1 r tot ρ i r tot
47
If it is assumed that the flexible segments are made from the same material and
that r tot is the same for each, then the equation reduces to
Ki – 1 γ 1 ζ 1 Kθ I 1 γi ζ i K θ I i
------------------------------------ = -------------------
1 i
- (4.50)
ρ1 ρi
γi ζ i Kθ
κ i = ---------------- i
(4.51)
ρi
This parameter depends upon the sub-class and the configuration. Table 4.4 gives the val-
ues for κ i for each flexible segment in each configuration. Substituting Equation (4.51)
into Equation (4.50) and rearranging for I1 , the equation of interest becomes
κiIi
I 1 = -----------------
- (4.52)
κ1 Ki – 1
This equation relates the two moments of inertia together allowing quick calculations of
any moment of inertia within the mechanism and can be used in connection with the force
Table 4.4 Values for Kx-1 and κx for each sub-class and configuration
48
4.2.4 Exceptions in the Force Parameter β
In the above derivation, Equations (4.41) and (4.42) are true for all configurations
except the Class 1B mechanisms. From Chapter 2, the force equation for Class 1B mecha-
nisms is defined as
k
F = ----2- Φ (4.53)
r2
where
EI2
k 2 = γK θ -------- (4.54)
l
In this case, k 2 depends upon the length of the flexible segment which is no longer
guaranteed to be related to r 2 . Therefore, the value for ζ is different depending upon the
specific mechanism used and the equation for the parameter β becomes
γK θ ( R + 1 )ζΦ
β = ----------------------------------- (4.55)
ρ
where ζ has the values as defined by Equations (4.10) to (4.12) and Equation (4.53)
becomes
βEI 2
F = -----------
2
- (4.56)
r tot
where I 2 is the first moment of inertia of the flexible segment associated with pivot two
49
4.2.5 Force Feasibility Example
This example is a continuation of the example found in Section 4.1.3. The mecha-
nism is required to have a force of 0.1 pounds and the width of the flexible segments can
be no more than 1.0 inch. If Beryllium Copper is chosen as the material (E=18.5 Mpsi), is
the mechanism feasible? The stress feasibility calculations were performed in a previous
example.
We know that the stress is feasible when Beryllium Copper is used. Therefore, the
question is whether or not the mechanism can generate the required force with the limita-
at Table 4.3. The value for β on a Class 2B-lps-a mechanism with k 1 coming from a
β = 7.263 (4.57)
Now, the moments of inertia between the third flexible segment and the first must
be related. This is done by using Equation (4.52) and looking up the values for each
23.17
I1 = -------------------------------- I 3 = 6.86I 3 (4.58)
3.368 ( 1.003 )
50
Substituting the formula for a rectangular cross section and Equation (4.58) into
12Fr tot 2
b = -----------------------
- (4.59)
6.86βEh 3
If the original thickness restriction of h > 0.001 inches is used, then the minimum
12Fr tot 2
b min > -----------------------
- (4.60)
6.86βEh 3
12 ( 0.1 )5 2
bmin > ------------------------------------------------------------------3- (4.61)
6.86 ( 7.263 ) ( 18.5e6 )0.001
or
This exceeds the acceptable width and therefore is not an acceptable design. How-
ever, if the maximum value for h max (0.0046) given in equation (4.34) is used, then b min
becomes
12 ( 0.1 )5 2
b min > ------------------------------------------------------------------------3- (4.63)
6.86 ( 3.0987 ) ( 18.5e6 )0.0046
or
51
The stress feasibility equation indicates that this combination of b and h allows
for the mechanism to meet the force requirement without violating the width require-
ments. It should also be noted that there are several combinations of h and b that will sat-
isfy both the hmax and b min inequalities. Assuming that the maximum b (1.0 in) is used,
then the h to be used with it for the third flexible pivot becomes 0.0032. This combination
will also give an acceptable value as will any values between these two sets.
The final part of the design process is to find the width and thickness of the first
flexible segment. Depending upon the constraints, it is an easy matter to use Equation
(4.58) and pick a value for one of the dimensions while solving for the other. In this man-
ner, the complete geometry of all of the flexible segments can be found while still satisfy-
ing the constraints of the problem. This example demonstrates the usefulness and ease
52
CHAPTER 5 MECHANISM
COMPARISONS
To effectively use the CFMs, it is necessary to understand the advantages and dis-
advantages of each configuration. One way to understand this is to compare the mecha-
nisms with each other to determine which ones are best suited under given circumstances.
This chapter presents several different comparisons that can be made to help better under-
• Actual lengths
• Maximum stiffness
• Percent constant-force
• Manufacturing orientations
• Normal displacement
53
5.1 Length Comparison
The equations developed in Chapter 4 are in terms of the total PRBM length r tot .
However, it is often desirable to express the equations for CFMs in terms of the actual
length of the mechanism, l tot . This can be done by using the following equation:
where λ is a length parameter. A value for λ for each mechanism can be found by exam-
ining them individually. The calculation of the value for the parameter λ for the class 1A-
lpp mechanism is presented as an example. Following which, the general equations and
The total length of the mechanism is found by adding the actual length of each link
as expressed in
ltot = l2 + l 3 (5.2)
where l 2 and l 3 are the actual link lengths. In this configuration, the first link is a long
flexible beam. The length of the beam is found by dividing the PRBM link length by γ or
r2
l 2 = ---- (5.3)
γ
The second beam is a rigid link with two pin joints and its length is expressed as
l 3 = r3 (5.4)
54
Combination of Equations (5.2) to (5.4) results in
r2
l tot = ---- + r 3 (5.5)
γ
To find the value for the parameter L , Equation (5.5) can be substituted into Equation
r2
---- + r 3
γ
λ = ---------------- (5.6)
r tot
where
r tot = r 2 + r 3 (5.7)
Once again, Equation (4.8) can be used to express the PRBM link lengths in terms
of the overall PRBM length. Substituting Equation (4.8) and the corresponding ζ values
Rγ + 1
λ = -------------------- (5.8)
( R + 1 )γ
Equation (5.8) expresses λ in terms of R and γ . Substituting the values for R and
0.8274 ( 0.85 ) + 1
λ = ----------------------------------------- = 1.027 (5.9)
( 0.8274 + 1 )0.85
This indicates that the actual length of the mechanism is 1.027 times longer than the
55
The length parameter is valuable in that it allows any equations that are expressed
in terms of r tot to be easily expressed in terms of l tot . This is done by substituting Equa-
tion (5.1) into the equation and using the appropriate λ value. This allows the actual
mechanism length to be used throughout the design process. This parameter also allows
for comparisons of the actual lengths of mechanisms that have the same PRBM length.
The technique described above can be used for all of the different configurations
and corresponding λ values can be found. Table 5.1 presents the general equation and
value of λ for each mechanism. Two mechanisms can be compared by comparing λ val-
ues. The mechanism with the higher λ value will have a larger actual length assuming that
stiffness of each configuration under identical sets of circumstances. By making this com-
parison, the mechanisms that provide the most force are identified, allowing them to be
The stiffness referred to in this work is a measure of the maximum force that can
be generated for a given mechanism size and stress level. To make this comparison, it is
56
necessary to derive a method to hold conditions such as size and stress constant while
βEI 1
F = -----------
2
- (5.10)
r tot
rial dependent. Equation (5.10) can be written in terms of l tot by using Equation (5.1),
resulting in
βEL 2 I1
F = -----------------
2
- (5.11)
ltot
λ Value
Configuration λ Formula
sub-class a sub-class b
R + 1.05
Class 1A-spp -------------------- 1.027 1.027
R+1
Rγ + 1
Class 1A-lpp -------------------- 1.097 1.094
γ(R + 1 )
Rγ + 1
Class 1B-plp -------------------- 1.088 1.088
γ(R + 1 )
R + 1.05
Class 2A-ssp -------------------- 1.036 1.035
R+1
R + 1.05γ
Class 2A-slp ----------------------- 1.086 1.088
γ(R + 1)
1.05R + 1.05
Class 2B-sps ------------------------------- 1.050 1.050
(R + 1)
1.05Rγ + 1
Class 2B-lps -------------------------- 1.122 1.119
γ(R + 1)
1.05R + 1.05
Class 3A-sss ------------------------------- 1.050 1.050
(R + 1)
57
A method could be pursued that would fix I 1 , l tot , and E to form a common set of
parameters. This would allow for comparison of the configuration based on F , which
would be directly proportional to β . Although this method would be valid, there are some
By making the assumptions listed above, the amount of stress found in each con-
figuration at the total displacement has not been taken into account. By holding I1 con-
stant, the maximum stress in each mechanism is not equivalent with the maximum stress
in the other mechanisms. Therefore, it would be ideal to compare the forces of each con-
figuration at the same stress value. This can be done by holding E , l tot , and the width b
constant for each mechanism while adjusting the height of the flexible segments until a
To equate the stress levels it is necessary to change the focus of the force equation
from the first flexible segment to the primary pivot, after which, the stresses can be
equated.
κiIi
I 1 = ------------------ (5.12)
κ1 Ki – 1
58
Substituting Equation (5.12) into Equation (5.11), and setting i to the primary pivot, the
βL 2 Eκ p Ip
F = --------------------------
- (5.13)
κ 1 K p – 1 ltot 2
The stress parameter developed in Chapter 4 can be used to equate the stress in
each mechanism to some maximum allowable stress value. The basic stress equation
found in Chapter 4 is
Ω
αA ≤ ------- (5.14)
SF
To make the comparison, the right hand side of the side of Equation (5.14) is set
αA = 1 (5.15)
Equation (5.15) assumes that the same material and safety factor are used for each
mechanism and that the stress level in each mechanism is the same. Furthermore,
cL hL
A = ------- = ---------- (5.16)
l tot 2l tot
To equate the stresses, either h or l tot can be adjusted to give the maximum stress.
However, if ltot is assumed to be constant as described above, then the comparison of the
α i h i L i = α j h j Lj (5.17)
59
Equation (5.17) can be used to maximize the stress in each configuration. Using
the general equation for a rectangular cross section, rearranging Equation (5.17), and sub-
α j λj 3
β i Eκ p λ i2 b i --------- -h
i α i λ i j
F i = --------------------------------------------------
2
(5.18)
12κ 1 K ( p – 1 ) l tot
i i
Equation (5.18) gives the force in configuration i at a stress level equal to the
stress in configuration j given identical material properties ( Ω ), safety factor (SF), overall
Ultimately, comparisons between all of the configurations, not just two, are to be made.
and ba to 1, and removing the 12. This results in a stiffness intensity parameter ψ where
1 3
βκ p ---
α
ψ = ----------------------- (5.19)
λκ 1 K p – 1
that describes the stiffness of the mechanism related to stress. This parameter can be used
to compare all of the configurations assuming identical lengths, materials, widths, and
stresses.
60
ψ Values given Equivalent Stress
2.500
Mech 1A-spp-b
2.000 Mech 1A-lpp-b
Mech 1B-psp-b
ψ value
1.500 Mech 1B-plp-b
Mech 2A-ssp-b
1.000
Mech 2A-slp-b
Mech 2B-sps-b
0.500
Mech 2B-lps-b
0.000 Mech 3A-sss-b
0 20 40 60
Percent Deflection (d )
Examination of Equation (5.19) shows that for a given configuration, ψ varies only
of displacement have been calculated for all of the configurations and sub-classes as intro-
duced in Chapter 3. Figure 5.1 shows a graph of the ψ values for the mechanisms in sub-
class b versus the percent displacement. This graph shows that the Mech 1A-lpp-b has the
highest value for the parameter ψ. The graph also shows how the available force decreases
increase, stresses will increase, limiting the height of the flexible segment. This in turn
decreases the moment of inertia, the corresponding spring constant, and ultimately, the
output force.
61
Normalized Ψ Values given Equivalent Stress
25
Mech 1A-spp-b
Normalized Ψ Values
20 Mech 1A-lpp-b
Mech 1B-psp-b
15
Mech 1B-plp-b
Mech 2A-ssp-b
10
Mech 2A-slp-b
Mech 2B-sps-b
5
Mech 2B-lps-b
Mech 3A-sss-b
0
0 20 40 60
Percent Deflection (d )
Figure 5.2 Normalized Ψ values for sub-class b
The data in Figure 5.1 can be normalized using the maximum value at maximum
ψ
Ψ = ------ (5.20)
ψc
Figure 5.2 shows a graph of the same data, but in normalized form. This graph
shows the values as percentages of the maximum value at full displacement, giving a bet-
ter understanding of the force behavior. For example, the graph shows that the force out-
This information helps show the trade off between deflection and stiffness of the mecha-
nisms, helping to indicate in which deflection ranges a mechanism must operate to obtain
a desired force. Additionally, Figure 5.2 shows the relationship between mechanisms. For
62
Normalized Ψ Values given Equivalent Stress
Mech 1A-spp-a
25 Mech 1A-lpp-a
Mech 1B-psp-a
20
Values
Mech 1B-plp-a
15 Mech 2A-ssp-a
Normalized
Mech 2A-slp-a
10
Mech 2B-sps-a
5 Mech 2B-lps-a
Mech 3A-sss-a
0
0 5 10 15 20
Percent Deflection (d )
example, mechanism 1B-plp-b has only 40% of the force or stiffness of mechanism 1A-
The same procedures were followed for the mechanisms in sub-class a. The nor-
malized results can be found in Figure 5.3. Comparison of Figure 5.2 and Figure 5.3
shows that the order of the mechanisms is the same for both sub-classes.
stresses are not equated are graphed in Figure 5.4. This graph shows that when the stresses
are not equated, mechanism 3A-sss-b is the stiffest mechanism, while mech 1A-lpp-b is
not even in the top three. It also indicates that the maximum stiffness is not a function of
63
Normalized Stress (Percent Non-Equivalent Stress
1.2
of Maximum Force)
1 Mech 1A-spp-b
Mech 1A-lpp-b
0.8
Mech 1B-psp-b
0.6 Mech 1B-plp-b
Mech 2A-ssp-b
0.4 Mech 2A-slp-b
Mech 2B-sps-b
0.2
Mech 2B-lps-b
0 Mech 3A-sss-b
0 20 40 60
Percent Deflection (d )
Figure 5.4 Normalized Ψ values for sub-class b without stress equalization
maximum deflection. However, this information is misleading because the stresses in the
Based on the results from above, the stiffness of each mechanism can be deter-
mined relative to all of the other mechanisms. The results indicate that the three stiffest
mechanisms are, in descending order, 1A-lpp, 1B-plp, and 2A-slp. When large output
forces are needed, these mechanisms should be the first ones considered. The rest of the
The parameter Ψ also adds a valuable tool in looking for new classes of mecha-
nisms. With the parameter, optimization routines can be designed to look only for mecha-
64
nisms that are stiffer than the mechanisms currently defined. The parameter Ψ also makes
The parameter Ψ is also beneficial because, when values are standardized, the
effect of percent deflection on the force can be determined quickly. This aids in design by
A second parameter that can be used for comparing configurations is the percent
max ( F )
Ξ = ------------------- (5.21)
min ( F )
To help this parameter be more intuitive, it can be inverted and multiplied by 100,
resulting in
min ( F )
Ξ' = 100 ------------------- (5.22)
max ( F )
Multiplying this parameter by a hundred gives the percent constant-force as a per-
centage with 100% being perfectly constant. Redefining this parameter make it more intu-
65
itive and easier to understand. It measures the amount of variation between the minimum
To use this parameter, it is important to ensure that the values are used in a consis-
tent manner. The maximum force, as defined in Equation (5.21), is taken to be the maxi-
mum force throughout the percent displacement specified for the sub-class. Due to the
nature of CFMs, this force is usually located at the maximum deflection. The minimum
force is defined similarly to the maximum force, and can generally be found at the small-
est deflection.
across the full range of displacement of the mechanism. This value will be termed Ξ' ex
and is
F0
Ξ' ex = ----------- (5.23)
Fmax
Often, Ξ' ex can not be calculated directly due to limitations in the models. These
limitations may include inflection points at zero displacement or limitations in step size.
For example, in a finite-element model (FEA), it may only be feasible to calculate the
force at five points along the deflection. In this case, d min may be some distance from the
zero displacement. Additionally, results from experimental testing often have a small
66
range at the beginning in which friction and other factors distort the output force. In these
A value for Ξ' ex can be determined by curve fitting a line through the maximum
force at the total displacement and the minimum force at the smallest known deflection.
The y-intercept of this line can be used as the force at zero deflection ( F 0 ).
Using the basic slope-intercept equation of a line, the maximum force is found to
be
F max – F min
F max = d max ----------------------------
- + F0 (5.24)
d max – d min
67
Fmax
Force (F)
Fmin
F0
Sample Curve
Curve Fit
dmin dmax
Percent deflection (d)
Figure 5.5 Parameters for Ξ'abs calculation
where F max is the force calculated at displacement d max (percent or actual displacement)
and F min is the force at the displacement d min . Figure 5.5 shows the parameters used to
calculated Ξ' ex . Rearranging for F0 and factoring d max out of the denominator gives
F max – F min
F 0 = Fmax – ----------------------------- (5.25)
d min
1–d
-----------
max
Fmin
1 – -----------
F max
F0 = Fmax 1 – --------------------
- (5.26)
d min
1 – -----------
d max
68
Finally, dividing F0 by F max gives
Fmin
1 – -----------
F0 Fmax
--------------------
Ξ' ex = 100 ⋅ ----------- = 1– - (5.27)
F max d min
1 – -----------
d max
This equation gives a level of constant-force for the entire deflection of the mecha-
nism. If the minimum force ( F min ) is already on the y-axis, then Equation (5.27) reduces
to Equation (5.22).
robust. It is assumed
F min ≅ F0 (5.28)
The original values for Ξ from Chapter 2 are shown in Table 5.2. These values
were calculated using the method illustrated in Equation (5.21). The values for the new
method, the inverse of the original, can be seen next to the old values. Examination of
these values show that the percent constant-force can now be referred to in a more intui-
From Table 5.2, it can be seen that the class 1B mechanisms have a smaller percent
constant-force than the other mechanisms. This indicates that these mechanisms will have
69
a greater variation in force than the other mechanisms. This knowledge will allow design-
ers to pick a mechanism suitable for the application. Additionally, as new sub-classes are
found, a method has been established to help identify which mechanisms exhibit a higher
percentage of constant-force.
cent constant-force. For the new mechanisms presented in this work, Ξ' is calculated and
presented as Ξ' ex since it can be assumed that inequality in Equation (5.29) is satisfied. In
the cases in which this inequality does not hold true, the extrapolated percent constant-
force Ξ'ex will be presented. This will help ensure the comparability of this parameter,
allowing for useful and accurate conclusions to be drawn. Additionally, the further useful-
ness of this parameter to determine the percent constant-force from experimental data will
70
(a) (b)
CFMs can be fabricated in many ways. The special case orientations are discussed
in this section.
The first special case CFM orientation is an in-plane orientation. In this orienta-
tion, all motion takes place in the plane of manufacturing. To simplify the manufacturing
of CFMs in an in-plane orientation, the widths of the flexible segments must be equal to
one another and to the thickness of the work material. This is shown graphically in Figure
5.6a.
5.6b. In this case, to simplify manufacturing, the thicknesses of the flexible segments must
be equal to one another and to the thickness of the material. This allows the mechanism to
71
be fabricated from a uniform piece of material with either simple milling or stamping type
operations.
A mechanism that is suitable for both types of fabrication orientations may be preferable
over a mechanism that is not suitable for either one of these fabrication orientations, all
When comparing different CFMs, the design space surrounding the flexible seg-
ments in the mechanisms should be compared. The moment of inertia value is set by the
spring constants requirements and can not be changed. However, there are many different
combinations of flexible segment thicknesses and widths that can give the correct moment
of inertia. It is advantageous to have combinations of thicknesses and widths that meet the
The thicknesses of the primary pivot is limited by the stress constraint. However,
as presented in Chapter 4, if the primary pivot is used during design, the other flexible seg-
ments in the mechanism must be limited to values that will maintain a stress level equal to
72
The relationship above is defined mathematically in Equation (4.22) as
c i ≤ C'c p (5.30)
where c p is half the thickness of the primary pivot, c i is half the thickness of the other
pivots in the mechanism, and C' is a parameter for each specific mechanism.
The parameter C' varies with displacement. For simplicity, and to ensure that the
inequality in Equation (5.30) is always satisfied, only the minimum values for C' are pre-
αp Mp d n Mp
p
C = minimum of -----
- = minimum of -------------n- = ------- d ( n p – no ) (5.31)
αo Mo d Mo
o
By selecting the minimum values, the inequality in Equation (5.30) will be satis-
fied at any displacement, ensuring that the primary pivot will always have the maximum
stress. The values for C for each mechanism are repeated in Table 5.3.
The other limitation on the design space for the thicknesses of the flexible seg-
ments is the minimum thickness of the mechanism. This limitation is defined by the limi-
73
5.5.2 Width Limits
The limitations on the width of the flexible segments are also important in deter-
mining the design space of a flexible segment. According to the first moment of inertia,
the width of the beam is inversely related to the thickness of the flexible segment where
b o h o3
I o = ----------- (5.32)
12
for rectangular cross-sections. As for the thicknesses, the moment of inertia of any given
flexible pivot is related to the moment of inertia of the primary pivot. Generalizing Equa-
Ki – 1κp
I i = ------------------ I p (5.33)
K p – 1 κi
Ki – 1κp
b i h i3 = -----------------
-b h3 (5.34)
K p – 1 κi p p
Equation (5.34) shows that the combinations of thickness and width are defined by
the geometry of the primary pivot. The maximum thickness for any pivot is limited first by
the minimum thickness, and second by the design constraints of the problem. Because
general manufacturing limitations for thicknesses vary for different processes, the upper
limit for widths is difficult to define. However, it is known that the limit can be relatively
large.
some cases, perhaps even design constraints. There are two important width ratios limits
74
that limit the width of a flexible segment and the width of the primary pivot. The first is
the minimum ratio, Dmin, which occurs when the thicknesses of the flexible segments are
set at their highest value. The second is the width ratio that occurs when the thicknesses
The minimum width ratio can be calculated by combining the equality portion of
Equation (5.30) and Equation (5.34) and solving for b i ⁄ b p . This results in
b Ki – 1 κ p
-----i = ------------------------
- = D min (5.35)
bp Kp – 1 κi C3
Values for the parameter Dmin are tabulated in Table 5.3. These values represent
the lowest possible ratio between the width of a flexible segment and the width of the pri-
The second ratio, Dequal, is important in cases when the mechanism needs to be
manufactured using an out-of-plane orientation. This ratio is similar to the one above and
is derived from Equation (5.34). Setting the thicknesses equal and solving for the ratio of
b Ki – 1 κp
-----i = -----------------
- = D equal (5.36)
bp Kp – 1 κi
This ratio indicates what ratio the widths of the flexible segments must be when
75
C
1
Dmin
Feasible Design Space
bx
bp Limited by Manufacturing Process
Limited by Mechanism
Direction to move for Improved Design Area
hx 1
hp
The design area for the flexible segments that are not the primary pivot are par-
tially defined by the parameters C and Dmin. The rest of the area is defined by manufactur-
ing limitations. A graphical representation of this design area can be seen in Figure 5.7.
The design area is bounded on the upper left corner by manufacturing limitations.
The limitations on the lower and right side are due to stress limitations. The mechanism
with the largest area, assuming the same limitations in manufacturing, will have the most
To improve the design area, the constants C and Dmin must move according to the
arrows in Figure 5.7. This provides for more suitable combinations of widths and thick-
nesses.
76
5.5.4 Determining Fabrication Orientations
The width and thickness ratios can also be used to determine suitable fabrication
orientations. For in-plane orientation, the parameter Dmin must be less than or equal to 1.
This ensures that the flexible segments can have the same width as the material, a neces-
sary trait for in-plane orientation. The magnitude of C in this case is not significant.
equal to 1 while Dequal is greater than or equal to 0.5 and less than or equal to 2. This
allows the thicknesses to be the same width as the material while the ratio between the dif-
ferent widths is not larger than double. The magnitude of Dmin is not significant in this ori-
entation.
thickness (C≥1), but the width ratio is smaller than 0.5 or larger than 2. In these cases, if
Dequal is greater than 0.10 and less than 10, the mechanism is considered possibly suitable
The different types of fabrication orientations mentioned above, along with the
required values for the thickness and width ratios, are summarized in Table 5.4.
In-plane NA NA ≤1
Out-of-plane ≥1 ≤ 2 (0.5) NA
Possibly out-of-plane ≥1 > 2 (0.5), ≤ 10 (0.1) NA
77
3
k2
r3
normal r2 slider
deflection ∆y deflection
∆x
k3
F
2
k1
r1
Figure 5.8 Normal displacement definition
This indicates that a method is needed to compare the normal displacement of different
configurations.
The normal deflection, ∆y , can be defined as the deflection normal to the input
deflection resulting from the outward deflection of one slider-crank mechanism. Figure
5.8 shows a graphical definition of this deflection. The normal displacement, dN , is the
ratio between the normal deflection ( ∆y ) at displacement d and the total undeflected
∆y
dN = 100 ⋅ ------- (5.37)
r tot
78
It is important to calculate d N at a displacement corresponding to the sub-class
maximum deflection percentage (16 for sub-class a, 40 for sub-class b). This value is
known as d Nmax . To determine d Nmax it is first necessary to develop an equation that will
determine the normal displacement d N in terms of r tot and d. First, the normal deflection
∆y = r 2 sin θ 2 (5.38)
r 12 + r 22 – r 32
θ 2 = acos --------------------------
- (5.39)
2r 1 r 2
Using Equations (4.8), (4.10), and (4.11), r 2 and r 3 can be related to r tot by
r tot
r 2 = ------------
- (5.40)
R+1
r tot
r 3 = ------------- (5.41)
1
1 + ---
R
d
r 1 = r tot 1 – --------- (5.42)
100
79
Substituting Equations (5.40), (5.41), and (5.42) into Equation (5.39), and rear-
ranging, results in
where r tot drops out of the equation. Substituting Equations (5.40) and (5.43) into Equa-
Equation (5.45) gives the normal displacement as a percentage of the total dis-
placement and in terms of R and d. This equation will allow a comparison of the normal
displacement of each of the mechanisms. As mentioned before, not only will this parame-
ter help compare different mechanisms, but it can be used during optimization to look for
better mechanisms.
The R values for the original mechanisms were used to calculate the normal dis-
placement at the sub-class maximum displacement percentage. The values calculated for
d Nmax are summarized in Table 5.5. These values represent the percent of the total length
that the mechanisms will displace at the maximum sub-class displacement. These values
80
can be used to help determine the space required for a particular design using a single
CFM. If the total displacement for a sub-class is not utilized, the normal displacement d N
By observing a general slider crank, it can be seen that the normal deflection ∆y
can never be larger than the smallest link, either r 2 or r 3 . Since the smallest link can
never be larger than 50% of the mechanism, the normal displacement can never be greater
than 50%.
sub-class a sub-class b
Mechanism
R d Nmax R d Nmax
1a 0.8853 26.96 0.8853 39.79
1b 1.0000 27.13 1.0000 40.00
2a 0.4323 23.23 0.4323 30.03
2b 0.8441 26.77 0.8441 39.60
3a 2.0821 22.82 2.0821 32.44
81
82
CHAPTER 6 MODEL AND
OPTIMIZATION
The development and improvement of CFMs is partially based upon the ability to
model the mechanisms and optimize for the correct parameters. It is through the model
and optimization that new sub-classes of constant-force parameters can be found. This
chapter describes the key features and important issues of the modeling and optimization
There are several directions from which to approach the construction of the CFM
model to be linked with the optimization software. Since the optimization will alter key
parameters rather than the mechanism geometries, which is opposite of a design approach,
the model to be used at this point in this research will be constructed differently from a
83
Although the model has the ability to produce valid designs, it was not developed
for primary design of CFMs. Modeling methods and tips for use in design are presented in
Chapter 8.
The CFM model is based on the behavioral model presented in Chapter 2, but has
been adapted to include more general modeling abilities and the parameters developed in
previous chapters.
The modeling methods presented in Chapter 2 are used as a starting point for the
modeling of CFMs in this research. However, some generalization of the model was made
software allows for easy use, quick alterations, data file manipulation, and linking with the
optimization software. Detailed information and a copy of the code used in Matlab can be
found APPENDIX B.
k
F = ----1- Φ (6.1)
r2
84
where
( R cos θ 3 [ θ2 + K 1 ( 2π + θ 2 – θ 3 ) ] + cos θ 2 [ K 1 ( 2π + θ 2 – θ 3 ) + K 2 ( 2π – θ3 ) ] )
Φ = ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
- (6.2)
R sin ( θ 2 – θ 3 )
or
k
F = ----2- Φ′ (6.3)
r2
where
( R cos θ 3 [ K 1 ( 2π + θ 2 – θ 3 ) ] + cos θ2 [ K 1 ( 2π + θ 2 – θ 3 ) + K 2 ( 2π – θ 3 ) ] )
Φ′ = ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- (6.4)
R sin ( θ 2 – θ 3 )
The inputs for the CFM were selected specifically to aid in the search for new sub-
classes. Some of the input parameters allow the optimization routine to alter the parame-
ters that most directly affect the sub-classes. The other input parameters are parameters
that have no direct impact on sub-classes, and in fact are often set at arbitrary values such
as 1. They are included merely to allow the model to be used for calculations.
85
• Flexible pivot types: pin, small-length, or long flexible
• Percent deflection: d
The inputs that have no direct affect on the configurations, but are required to ana-
• Material Properties: E, Sy
These inputs allow the optimization software to have the best access to modifying
the mechanism and also allow the model to be used to generate real designs.
The CFM model determines all the important parameters needed to define a new
configuration. For a couple of these parameters, values are calculated at each point of the
deflection, but only the average value of the parameter is reported - as is done with the
• Comparison Parameters: λ, D
• Stiffness Parameter: Ψ
86
• Primary Pivot
• Θ
• Force Parameter: β
spread sheet, which in turn was used to verify the output of the model. All 6 classes of
mechanisms were used at both sub-classes a and b. In all 12 cases, the model agreed with
the spread sheet and correctly calculated the parameter values presented in Chapters 4 and
5.
6.2 Optimization
parameters that provides the most desirable CFM performance. The performance is mea-
The prior work was concerned with finding mechanisms that have the highest per-
cent constant-force. For this research, several different objectives are possible depending
87
on the desired output. Different optimization problems are defined below, at which point
In several cases, dual objective functions could be used. However, the use of dual
objective functions complicates the optimization process. Where possible, dual objective
functions are reduced to a single objective function with the other objective function
becoming a constraint.
6.2.2 Variables
The analytical variables for the optimization are all of the model inputs described
above. In the search for new configurations, all comparisons are made so that material
properties and mechanism size are not important. Therefore, these variables are set equal
to 1. The variables that control the mechanism type and sub-class are set to the respective
The only variables that are used as design variables are the constant-force parame-
ters. These parameters allow the optimization to modify the link length ratio, R, and the
6.2.3 Functions
The analytical functions for the optimization problem include all of the model out-
puts described above. The selection of analytical functions as design functions depends
upon the objective of the optimization problem. For each optimization problem described
88
6.2.4 Optimization Problem - Stiffer Mechanisms
To develop mechanisms that are stiffer, thus allowing for higher forces with the
same stress limits, it is important to maximize the stiffness parameter Ψ . In this case, to
eliminate the dual objective functions, the parameter Ξ' can be established as a design
Maximize Ψ (6.5)
subject to
The optimization program can look for mechanisms that have a smaller normal
displacement by minimizing the normal displacement parameter. Once again, the percent
Minimize dN (6.7)
subject to
To determine the best configurations that are suitable for in-plane and out-of-plane
manufacturing orientations, it is necessary that the width and thickness ratios satisfy the
89
criteria defined in Section 5.5.4 and tabulated in Table 5.4. Formally, the optimization
subject to
0.5 ≤ D equeal ≤ 2
(6.11)
The optimization problems described above are summarized in Table 6.1. Not all
of the possible optimization problems that can be examined for improvements in CFMs
are represented in the table. In fact, the problems presented above can be combined in
many different ways. One such way is shown in Table 6.1 under the title “Best Method”.
This problem defines an optimization problem that searches for the stiffest mechanism
that is suitable for in-plane and out-of-plane orientation, as well as a certain percent con-
stant-force value.
Normal
Objective Stiffness Stamping Milling Best Overall
Displacement
Objective Function Maximize Ψ Minimize d N Maximize Ξ' Maximize Ξ' Maximize Ψ
Design Functions Ξ' ≥ Ξ'c Ξ' ≥ Ξ'c C ≥ 1 D min ≤ 2 Ξ' ≥ Ξ'c
D equal ≤ 2 C≥ 1
D equal ≤ 2
D min ≤ 1
90
The optimization problems presented in the chapter are presented as examples and
guidelines. The optimization problems used may not be suitable for all mechanism types,
thereby requiring some alterations in the problems. The next chapter, Chapter 7, will iden-
tify the particular problems used in the search for better mechanisms within each configu-
91
92
CHAPTER 7 NEW MECHANISMS
This chapter presents the results of the optimization, including explore plots and
One of the main objectives of this research is to find new and improved CFMs,
principally CFMs stiffer than the original mechanisms. To accomplish this, each configu-
ration was examined through optimization for new mechanisms which were evaluated by
Two main objectives existed for the optimization: find the stiffest mechanisms for
three different levels of constant-force (90, 95, and 99 percent constant-force), and find
the stiffest mechanisms for the same four levels of constant-force that are suitable for fab-
In many cases, suitable mechanisms could not be found for several of these objec-
93
7.1.1 Sub-class Expansion
Two new additions to the sub-class definition are needed. First, the limit for the
level of constant-force used in the optimization is added after the percent displacement let-
ter. This number identifies the mechanism and the constraint used. The number does not
give the actual percent constant-force value of the mechanism, only the constraint value
used.
The second addition identifies the mechanism as suitable for fabrication in an in-
plane orientation, “I”, an out-of-plane orientation, “O”, and/or possibly suitable for out-of-
plane orientation, “o”. These letters, when added after the percent deflection letter and
percent constant-force limit, identify the mechanism as suitable for these orientations,
the classification for the mechanism would be Class 2A-ssp-a90Io. This name summarizes
important information and provides a unique naming method for the mechanisms.
7.2 Class 1A
Class 1A was the first class of mechanisms examined. This class consists of both
94
7.2.1 Optimization Details
The mechanisms that make up this class only have one flexible segment, thus
making it the primary pivot. Additionally, with only one pivot, the mechanism is suitable
The constant-force parameters K1 and K2 are both zero, leaving only the parameter
was performed to find the stiffest mechanisms for each of the three new sub-classes with
An explore plot with R as the variable was generated for both configurations and
deflection ranges. The constant-force parameter R was allowed to vary from 0.5 to 3.0.
Once the explore plot was generated, graphs of the relationship of Ξ' ex and Ψ with respect
to R were generated. These plots are shown in Figure 7.1a. In both configurations, Ξ' ex
the optimization and shows that the highest percent constant-force occurs around an R
value of 0.8. This is consistent with the findings of the prior work. The graphs in Figure
7.1 also shows that in order to increase Ψ, R must be decreased, sacrificing the level of
constant-force.
95
.0777
1.96
99.3
99.3
.0629
1.59
93.9
93.9
1 1
.0482
1.22
88.6
88.6
2-psi
2-psi
1-Xi
1-Xi
.0335
.848
83.2
83.2
1 1
2 2
.0188
.477
77.8
77.8
1 1
2 2
.00403
.107
72.4
72.4
2 2
.500 1.00 1.50 2.00 2.50 3.00 .500 1.00 1.50 2.00 2.50 3.00
R R
(a) (b)
Figure 7.1 Explore plot of R versus Ξ'ex and Ψ for (a) Class 1A-lpp-a and (b) Class 1A-spp-a
The optimization was able to improve the stiffness of both configurations within
Class 1A for both deflection ranges. The results of the optimization are given in Table 7.1.
The Class 1A-lpp-a and Class 1A-spp-a mechanisms both showed a 50% increase in stiff-
96
.0996
95.3
2.51
95.3
.0810
78.5
2.04
78.5
1
2
.0624
61.8
1.58
61.8
1-Xi'
2-Psi
1-psi
2-Xi
2 1
.0438
45.1
1.11
45.1
2 1
.0251
28.4
.640
28.4
1
2
.00651
1 2
11.6
.172
11.6
1 2
.500 1.00 1.50 2.00 2.50 3.00 .500 1.00 1.50 2.00 2.50 3.00
R R
(a) (b)
Figure 7.2 Explore plot for (a) Class 1A-lpp-b and (b) Class 1A-spp-b
ness when Ξ' ex was allowed to drop to 90. However, the spp-a configuration still only has
The Class 1A-lpp-b and Class 1A-spp-b mechanisms had similar trends as the sub-
class a mechanisms. However, at 90% constant, these mechanisms only showed a 13%
increase in stiffness. Additionally, no feasible mechanism was found for both configura-
tions at 99 percent constant-force. This is supported by the explore graph generated for
these mechanisms. As seen in Figure 7.2, the highest Ξ' ex value is around 95, the same as
The optimization has defined stiffer mechanisms for three new sub-classes with
16% deflection and two new sub-classes with 40% deflection. If the design does not pre-
vent the level of constant-force from being sacrificed, these new sub-classes can be used
to obtain a higher stiffness in the mechanism. Table 7.2 summarizes the important design
parameters for the Class 1A CFMs which can be used, in correlation with the equations
97
developed earlier, to quickly design any new mechanism. The width and thickness ratios
are not applicable for these configurations and are emitted from the table.
7.3 Class 1B
tions. These configurations also only contain one flexible pivot, the second pivot, indicat-
ing that this pivot is the primary pivot and all configurations are suitable for both in-plane
For the optimization, K1 must be set equal to 1, while K2 is set equal to 0. Again,
this leaves only R to be used as a design variable. The optimization is not complicated and
new mechanisms for the three new sub-classes are easily defined.
98
.0151
.0373
94.6
86.2
12
.0146
.0328
92.9
82.5
1
1
2
.0142
.0283
91.2
78.8
2
2-Psi
2-Psi
1-Xi
1-Xi
.0137
.0238
89.4
75.2
1
2
.0132
.0193
87.7
71.5
2
1
1
.0128
.0148
2
85.9
67.8
.500 1.00 1.50 2.00 2.50 3.00 1.00 2.00 3.00 4.00 5.00
R R
(a) (b)
Figure 7.3 Explore plots for (a) Class 1B-psp-a (b) Class 1B-psp-b
Four explore plots were generated, one each for every combination of configura-
tion and displacement percentage. Figure 7.3 shows the explore plots for the psp configu-
ration. The peaks in the parameter Ξ' ex are clearly visible and correspond to the values
found in the original works. Additionally, according to Figure 7.3a, Ψ peaks at the same
point as Ξ' ex for the Class 1B-psp-a mechanism, thus indicating that no improvements are
possible. Examination of Ψ in Figure 7.3b shows that increases in Ψ occur only once Ξ' ex
has decreased below 75. This indicates that stiffer mechanisms for the three new sub-
The explore plots for the plp configurations in Figure 7.4 appear similar to the
explore plots in Figure 7.3. In fact, within the same percent deflection and between the
two configurations, the curves for Ξ' ex are identical. The differences are in the Ψ curves.
99
.710
94.6
86.2
1.80
1
.574
90.3
82.5
1.47
2
1
.439
86.0
78.8
1.13
2-Psi
1-Psi
1-Xi
2-Xi
.303
81.6
75.2
.788
1
2
.167
77.3
71.5
.449
2
2 1
2
.0312
73.0
67.8
.110
2 1 1
1.00 2.00 3.00 4.00 5.00 1.00 2.00 3.00 4.00 5.00
R R
(a) (b)
Figure 7.4 Explore plots for (a) Class 1B-plp-a (b) Class 1B-plp-b
This is consistent with expectations as Ξ' ex is based on the PRBM which does not change
between configurations within a certain class, while Ψ is based on the actual configuration
of flexible segments.
Once again the peaks in Ξ' ex correspond with the original work, indicating no
improvement in Ξ' ex . However, the peaks in the Ψ curves do not correspond to the peaks
in the Ξ' ex curve as they did with the psp configurations. This offset between the Ξ' ex and
Ψ curves provides an opportunity for trade-offs between stiffness and percent constant-
force.
The results of the optimization are summarized in Table 7.3. As expected, no feasi-
ble mechanisms for the psp with either percent deflection were found that offered a higher
100
For the Class 1B-plp-a mechanisms, feasible mechanisms for percent constant-
force values of 95 and 99 could not be found. Therefore, stiffer mechanisms were searched
for at percent constant-force levels of 80, 85, and 90. The Class 1B-plp-a90 mechanism
showed an increase of 70% while the Class 1B-plp-a80 mechanism increased in stiffness
by 88%.
Improved mechanisms for Class 1B-plp-b were found for 80 and 85 percent con-
stant-force values. These mechanisms showed an increase of 22% and 45% in stiffness
from the original mechanism. The summary of the parameters for the Class 1B mecha-
101
7.4 Class 2A
The Class 2A mechanisms consist of the slp and ssp configurations. These mecha-
nisms are combinations of the Class 1A-spp and Class 1B-psp and pl3p configurations.
mining which mechanisms to include in this work. This difficulty comes from the intro-
duction of a second design variable because of the use of two flexible segments.
Additionally, the ratios between the thickness and width of these two flexible segments are
The parameters R and K1 are the two design variables for the optimization, while
the thickness and width ratios for the first two pivots are added as design functions. Multi-
102
.00991
.0458
.0962
.140
.233
.341
.420
.420
.341
.264
.233
.153
.140
80.1
86.4
91.1
93.3
94.5
94.5
94.0
93.3
92.2
91.1
90.1
88.0
86.4
20.0
20.0
16.0
16.0
.0962
12.0
12.0
K1
K1
.0578
8.00
8.00
.0458
84.3
4.00
4.00
80.1 .00991
.500 1.00 1.50 2.00 2.50 3.00 .500 1.00 1.50 2.00 2.50 3.00
R R
(a) (b)
Figure 7.5 2-D explore plots for slp-a of R and K1 with contours of (a) Ξ'ex and (b) Ψ
The slp-a configuration was considered first. Two 2-D explore plots were gener-
ated for this mechanism and are shown in Figure 7.5 with R and K1 on the horizontal and
vertical axes. These plots show contours of Ξ' ex in Figure 7.5a and Ψ in Figure 7.5b.
Examination of these plots show that Ξ' ex is fairly independent of K1 but highly depen-
dent on R. Additionally, K1 has very little effect on Ψ for values of R below 1, but a large
These observations are valuable and give insight into the mechanism. The value
for K1, since it has little effect on Ξ' ex , should be made as large as possible. However, by
increasing K1, the ratio of the spring constants is increased, affecting the width and thick-
103
.0122
.0135
.0145
.0148
.0148
.0145
.0140
.0135
.0128
.0122
.0118
69.7
75.4
85.4
85.4
80.1
75.4
69.7
68.4
67.6
67.1
12.0
12.0
66.6
8.00
8.00
K1
K1
.0118
65.6
.00973
4.00
4.00
62.4
1.00 2.00 3.00 4.00 5.00 1.00 2.00 3.00 4.00 5.00
R R
(a) (b)
Figure 7.6 2-D explore plots for ssp-a with R vs. K1 and contours of (a) Ξ' ex and (b) Ψ
ness ratios. Ultimately, K1 can only be increased so much before fabricating the mecha-
The 2-D explore plots for slp-b showed similar relationships among the parame-
The ssp-a and ssp-b configurations were examined in the same way as the slp con-
figurations. The 2-D explore plots for the ssp-a configurations can be found in Figure 7.6.
These plots show that the contours of Ξ'ex are similar to those found in Figure 7.5a. How-
ever, the contours for Ψ are different. In this mechanism, the area in which Ψ is the high-
est corresponds to the area with the highest Ξ' ex . Additionally, in this area, as K1
104
The ssp-b explore plots, shown in Figure C.2 in Appendix C, display the same
stiffer mechanisms and the definition of mechanisms suitable for in-plane and out-of-
The original slp-a sub-class was suitable for in-plane orientations, and possibly
suitable for out-of-plane orientations. A new sub-class was found that had the same level
of constant-force, was suitable for both in-plane and out-of-plane orientations, and had a
50% increase in stiffness. Another sub-class was found for this configuration that allowed
the level of constant-force to decrease to 90, but increased the stiffness by 268%. This rep-
105
resents a fairly large increase in stiffness, with only a marginal decrease in percent con-
stant-force. The other new sub-classes for slp-a are shown in Table 7.5.
In the slp-b configuration, the original sub-class was the best mechanism that
could be found for the given value of Ξ' ex . However, a decrease of 10 in Ξ' ex provided an
The original ssp-a configuration was improved by 141% without changing Ξ' ex .
However, the original sub-class was possibly suitable for out-of-plane orientation, while
the new sub-class was not suitable for out-of-plane orientation. When the out-of-plane ori-
entation constraints were removed from the optimization, the stiffest mechanism that
could be found with Ξ' ex values between 90 and 100 was the a95M sub-class. This sub-
class has a Ξ' ex value of 96.2 and a 176.9% increase in stiffness over the original sub-
class. However, forcing the optimization to look for sub-classes suitable for in-plane and
The final configuration in this class of CFMs is the ssp-b. As with the slp-b, the
original sub-class was the stiffest mechanism that could be found for 99 percent constant-
force. However, a sub-class, not suitable for in-plane orientation and only possibly suit-
able for out-of-plane fabrication orientation, was found for 90 percent constant-force with
106
In all, the results of the optimization are promising. Many new sub-classes are pre-
sented in Table 7.5. The important design parameters for these new sub-classes are sum-
marized in Table 7.6 while the thickness and width ratios determined from the
Configuration Sub-Class Primary Pivot C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min
slp aIo 2 1.71 1.00 - 2.54 1.00 - 0.51 1.00 -
a99IO 2 1.13 1.00 - 1.43 1.00 - 1.00 1.00 -
a95Io 1 1.00 1.92 - 1.00 7.04 - 1.00 1.00 -
a90I 1 1.00 4.82 - 1.00 119.69 - 1.00 1.07 -
a90IO 2 1.00 1.00 - 1.00 1.00 - 1.00 1.00 -
slp bIO 2 1.47 1.00 - 1.98 1.00 - 0.63 1.00 -
b99
b95IO 1 1.00 1.09 - 1.00 1.28 - 1.00 1.00 -
b90Io 1 1.00 1.47 - 1.00 3.14 - 1.00 1.00 -
ssp aIo 2 5.07 1.00 - 7.52 1.00 - 0.06 1.00 -
a99I 2 2.54 1.00 - 11.20 1.00 - 0.68 1.00 -
a95I 2 2.18 1.00 - 10.43 1.00 - 1.01 1.00 -
a95IO 2 2.09 1.00 - 2.00 1.00 - 0.22 1.00 -
a90IO 2 1.33 1.00 - 2.00 1.00 - 0.84 1.00 -
ssp bIo 2 4.63 1.00 - 6.24 1.00 - 0.06 1.00 -
b99
b94I 2 6.67 1.00 - 15.38 1.00 - 0.05 1.00 -
b90IO 2 2.30 1.00 - 2.00 1.00 - 0.17 1.00 -
b90o 2 2.15 1.00 - 10.35 1.00 - 1.05 1.00 -
107
7.5 Class 2B
The lps and sps configuration make up the 2B class of CFMs. The original config-
urations in this class of CFMs exhibit the lowest values for percent constant-force of all
The optimization of the Class 2B configurations was done similar to the Class 2A
configurations with the only difference being a broadening of focus from increasing stiff-
The optimization of the lps found many different sub-classes that looked promis-
ing, making the task of selecting sub-classes to be presented in this work challenging. An
optimum plot was generated to help view the design area, and is shown in Figure 7.7. The
plot shows the values for Ξ' ex and Ψ , as well as the optimal value for K2, for any given R
value. Two peaks in Ξ' ex are present in the graph. The peak centered at an R value of 2.50
corresponds with higher values of Ξ' ex than the peak centered at 0.75. Therefore, mecha-
nisms within this area were considered. The optimum plot for the lps-b mechanism can be
108
14.1
.508
99.8
2
11.3
.411
97.9
8.52
.314
95.9
2-Psi
3-K2
1-Xi
2 2
1
5.72
.218
93.9
3
2.91
.121
92.0
1
3
.0246
.100
90.0
3 1
.500 1.00 1.50 2.00 2.50 3.00
R
Figure 7.7 Optimum plot for lps-a which shows Ξ'ex, Ψ and optimal K2 for given R values
The change in the type of flexible segments from the lps to the sps greatly affects
the optimization problem and the associated results. Explore plots, similar to those gener-
ated for the Class 2A configurations were generated for the sps-a configuration. The high-
est Ξ' ex and Ψ values are located around an R value of 1 and a K2 value of 1, as seen in
Figure 7.8. It is in this area that the search was focused and new sub-classes were found.
Similar plots and results were found for the sps-b configuration and can be seen in Figure
C.4 in Appendix C.
ments are summarized in Figure 7.8. For the lps-a configuration, improvements ranged
109
.00837
.00551
.00387
.00284
.00216
.00169
.0119
.0180
.0270
.0418
.0418
.0270
.0180
.0119
10.0
84.2
90.9
95.2
98.5
98.5
97.0
95.2
93.7
93.7
95.2
97.0
98.5
98.5
97.0
95.2
93.7
90.9
87.6
84.2
80.8
10.0
8.00
8.00
76.2
6.00
6.00
K2
70.6
K2
4.00
4.00
2.00
2.00
.500 1.50 2.50 3.50 4.50
R .500 1.50 2.50 3.50 4.50
R
(a) (b)
Figure 7.8 2-D explore plots for sps-a with R vs. K2 and contours of (a) Ξ'ex and (b) Ψ
110
and manufacturing process constraints. In all, six new sub-classes were added for 16%
deflection, three suitable for both in-plane and out-of-plane orientations, and three uncon-
The lps-b configurations had increases in stiffness from 539% to 3177%. This rep-
resents a significant improvement. In fact, the b95I and b90I sub-classes have higher stiff-
ness parameters than any of the original configurations and sub-classes. Additionally, the
level of constant-force was increased from the original value of 84 to a high of 99. In all,
The level of constant-force was increased from 93 to 99 while stiffness for the same sub-
class was increased 123%. The a90IO sub-class was constrained to greater than or equal to
90 percent constant-force. However, the result was a sub-class suitable for out-of-plane
Very little increase in stiffness was found for the sps-b configuration. For the
b86IO sub-class, stiffness was increased by only 68%. However, a large increase from the
All of the important design parameters for the CFM Class 2B are summarized in
Figure 7.9, with the thickness and width ratios summarized in Figure 7.10.
111
7.6 Class 3A
The last class of CFMs is Class 3A, consisting of only the sss configuration.
The optimization of the sss configuration requires that R, K1, and K2 be design
variables. Additionally, all of the width and thickness ratios must be added as design func-
Configuration Sub-Class Primary Pivot C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min
lps aIo 3 9.07 - 1.00 6.86 - 1.00 0.01 - 1.00
a99I 3 0.68 - 1.00 0.32 - 1.00 1.00 - 1.00
a95I 3 0.79 - 1.00 0.50 - 1.00 1.00 - 1.00
a90 1 1.00 - 0.84 1.00 - 0.82 1.00 - 1.37
a99IO 3 1.00 - 1.00 0.47 - 1.00 0.47 - 1.00
a95IO 3 1.00 - 1.00 0.66 - 1.00 0.66 - 1.00
a90IO 3 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00
lps bIo 3 7.33 - 1.00 6.05 - 1.00 0.02 - 1.00
b99IO 3 0.99 - 1.00 0.55 - 1.00 0.56 - 1.00
b95IO 3 0.99 - 1.00 0.60 - 1.00 0.61 - 1.00
b90IO 3 1.00 - 1.00 0.68 - 1.00 0.68 - 1.00
b95I 1 1.00 - 0.13 1.00 - 0.00 1.00 - 0.73
b90I 1 1.00 - 0.12 1.00 - 0.00 1.00 - 0.98
sps aIO 3 1.74 - 1.00 1.31 - 1.00 0.25 - 1.00
a99I 1 1.00 - 0.52 1.00 - 0.14 1.00 - 0.98
a90IO 3 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00
sps bIo 3 1.40 - 1.00 1.16 - 1.00 0.42 - 1.00
b98Io 3 2.35 - 1.00 6.01 - 1.00 0.46 - 1.00
b86IO 3 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00
112
.0865
99.9
2
.0697
97.9
.0529
95.9
1
2-Psi
1-Xi
.0362
94.0
2
.0194
92.0
2
.00257
90.0
1
.500 1.50 2.50 3.50 4.50
R
Figure 7.9 Optimum plot of design variables K1 and K2 versus R for sss-b with curves of Ξ'ex and
Ψ
tions. The increase in the number of design variables and design functions complicates the
optimization.
The addition of a third design variable limits the usefulness of 2-D explore plots.
Therefore, an optimum plot, shown in Figure 7.9, was generated with R as the independent
variable and K1 and K2 as design variables. This plot clearly shows the local optima, the
peaks of both Ξ' ex and Ψ , and their correlation. Both curves have peaks around an R
value of 1.5. It is in this area that the optimization was concentrated. The explore plot for
sss-b is similar to Figure 7.9 and can be found in Figure C.5 in Appendix C.
113
7.6.3 Optimization Results
Four new sub-classes are added for the sss-a configuration. The a99 sub-class has
the same level of constant-force, but shows an increase of 291% in stiffness. The a95IO
and the a90IO sub-classes are both suitable for in-plane and out-of-plane orientations. The
a95I sub-class is not suitable for out-of-plane orientation, but shows a 317% increase in
stiffness over the original sub-class. The results, plus the results for sss-b, are summarized
in Table 7.11.
For the sss-b configuration, three sub-classes are added. The b99 and b90I sub-
classes showed a 14% and 158% increase in stiffness respectively with the b90I sub-class
finding the stiffest mechanism to be 96 percent constant-force. The final sub-class, b88IO,
is the highest percent constant-force mechanism found suitable for in-plane and out-of-
Once again, the important parameters have been summarized and can be found in
Table 7.12, with width and thickness ratios found in Table 7.13.
114
7.7 Optimization Results Summary
Many improvements in stiffness were made within each class of CFMs. In some
classes, increases in stiffness of as much as 3000% were made. However, the maximum
stiffness of CFMs did not see much increase. Originally, the Class 1A-lps-a mechanism
Configuration Sub-Class Primary Pivot C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min
sss a 1 1.000 1.332 7.099 1.000 1.832 33.744 1.000 0.775 0.094
a99 1 1.000 0.745 2.222 1.000 0.435 12.492 1.000 1.052 1.139
a95IO 2 1.462 1.000 2.693 0.500 1.000 2.000 0.160 1.000 0.102
a95I 3 0.340 0.289 1.000 0.039 0.022 1.000 0.997 0.895 1.000
a90IO 2 1.260 1.000 3.104 1.974 1.000 1.979 0.986 1.000 0.066
sss b 1 1.000 1.041 4.338 1.000 1.541 19.533 1.000 1.365 0.239
b99 1 1.000 0.936 3.502 1.000 0.995 13.377 1.000 1.215 0.311
b90I 3 0.589 0.381 1.000 0.148 0.055 1.000 0.724 0.996 1.000
b88IO 2 1.610 1.000 2.398 0.500 1.000 2.000 0.120 1.000 0.145
115
Review of the optimization results indicates that three 16% deflection mechanisms
and four 40% deflection mechanisms have a stiffness parameter greater than 1. These
mechanisms are defined in Table 7.14, along with the percent increase in the stiffness
parameter. The largest increase of 49% occurred in Class 1A-lpp-a90. All three of the sub-
classes with 16% deflection belong to Class 1A-lpp. Two of the 40% deflection mecha-
nisms belong to Class 1A-lpp, while the other two belong to Class 2B-lps.
These increases in stiffness will provide designers with greater flexibility and addi-
tional options for design. The design space for the existing configurations has been
explored and the best mechanisms for a variety of circumstances have been defined.
116
CHAPTER 8 DESIGN APPROACH
AND METHODS
used in these applications, the CFMs must first be designed to meet all of the application’s
process used by designers. This chapter looks at the design aspect of CFMs including gen-
eral design procedures, infeasible mechanism elimination methods, secondary issues, and
trends between variables. The chapter ends by presenting different design example prob-
Every CFM design problem can be approached with the same general method.
This allows designers to quickly become familiar with the process. The length of the pro-
cess varies depending upon the constraints and requirements of the problem. In some
cases, an immediate and direct solution can be found while in others, iteration must be
117
8.1.1 Background, Assumption, and Limitations
The design method outlined in this chapter utilizes the stress and force parameters,
along with other parameters, developed in earlier chapters. Since these parameters attempt
to remove PRBM details such as individual link lengths, knowledge of the PRBM is not
This design method builds upon the assumptions used throughout this work.
Although the general steps outline in this chapter can always be used, the tabulated param-
Additionally, since the parameters for each mechanism use estimated values for
the PRBM and other parameters, the accuracy of the designs that come from this design
method cannot be guaranteed. For the most accurate design, a complete CFM model must
be constructed and more accurate approximations based on the given design problem must
be utilized. Wittwer (2001) addressed some of these accuracy issues and methods to deter-
There are three principal equations that are necessary in designing a CFM. These
equations are combinations and modifications of equations presented earlier in this work.
The first equation is a modified force feasibility equation. In Chapter 4, the force feasibil-
βEI 1
F = -----------
2
- (8.1)
r tot
118
Substitution of Equation (4.52) (moment of inertia ratios), Equation (5.1) (length parame-
ter equation), and the equation for the moment of inertia for a rectangular cross-section
λ 2 κ p b p h p3 βE
F = --------------------------------
2 κ K
- (8.2)
12ltot 1 p–1
The only exception to Equation (8.2) is for Class 1B mechanisms. The force
βλ 2 EI 2
F = -----------------
2
- (8.3)
l tot
The second principal design equation is a modified form of the stress feasibility
Ω
αA = ------- (8.4)
SF
λh
Substitution of α = Md n and A = ---------- into Equation (8.4) results in
2l tot
Ω- h
------ = λMd n ---------- (8.5)
SF 2l tot
119
Isolated Variables Coupler Variables Classification Variables
found in Equation (2.16) and Equation (5.1), the length parameter equation. The third
equation is
dl tot
disp = ------------ (8.6)
100λ
This equation is an intermediate equation which helps to provide needed values for vari-
ables in both the force and stress design equations. It is referred to as the displacement
equation.
These three equations are coupled through different variables and must be solved
The boxed variables in Figure 8.1 are isolated variables, or those variables that
occur only in one of the three equations. There are a total of five isolated variables in the
120
three equations, two in each of the force design equation and stress design equation, and
one in the deflection equation. For the system to be solvable, at least one of the two iso-
lated variables must be known in both the force design and stress design equations.
The circled variables in Figure 8.1 are shared between two or more of the principal
equations and are termed coupler variables. The number of coupler variables that must be
known to solve the system depends on which coupler variables are known and how many
isolated variables are known. No easy method, other than examining the missing variables
in Figure 8.1, has been determined for deciding how many coupler variables are needed.
The final type of variable in the principal equations are the classification variables.
These variables, outlined by a hexagon in Figure 8.1, can only be determined by choosing
a specific classification of mechanism and substituting the values for the variables into the
equations. If the classification is not known, then efforts must be made to reduce the num-
ber of possible classifications, and an iterative process may be required. Methods for elim-
The three types of variables in the principal equations can have three types of val-
ues: known, constraint, or unknown. A known variable value is one that the designer must
A constraint variable value is one that has no exact value given by the design
121
can be used in the principal equations to help determine design feasibilities and an initial
design.
An unknown variable value is one that has no set value and is not constrained by
upper and lower bounds. These values can be solved for or chosen as free variables if
before the equations can be solved. The following steps outline the order in which types of
1. Identify a suitable mechanism. If the design does not specify a specific mecha-
nism, then one must be chosen. It is ideal to choose the classification that is
best suited for the design problem. The following points should be considered
These points can be used to reduce the number of mechanisms by eliminating any
122
2. Fill in all the known values. Start by identifying all of the known values and
their location in the principal design equations. Figure 8.1 is useful for this pur-
pose.
3. Add constraint values, starting with the most important ones, until the equa-
tions can be solved. From Figure 8.1, determine how many more variables are
required to solve the equations. Use the constraint values of the design problem
to fill in the needed number of variables. The constraint value chosen depends
on the needs of the design objectives. Section 8.4 can be used to identify the
desired performance. At this point, the process will become iterative until a
4. Use unknown variables as free variables if needed. If, after adding all con-
straint variables, there are still not enough variable values to solve the equa-
tions, use reasonable assumptions for unknown values until there are enough to
These steps, along with the steps in the following sections, are summarized in the
flow chart in Figure 8.2. This flow chart, along with the diagram in Figure 8.1, will help
solve the system of equations for basic design cases. Examples of how to use these tools
123
Suitable Yes Add constraint values Use unknown
Identify a suitable Fill in all known
mechanism (start with the most variables as free
mechanism values
known important ones) variables if needed
No
No
Remove variables
Solve principal
if system is
design equations
over-constrained
Determine other
flexible segment
geometries
The principal equations listed above identify only the geometry of the primary
pivot. The following equations can be used after the basic design of the mechanism is
complete to determine, and if needed, verify the feasibility of the other flexible segments’
geometries. Of course, for the Class 1A and 1B mechanisms, this is not necessary.
The important equation for the design of the other flexible segments is
Ki – 1 κ p I
Ii = ----------------------p- (8.7)
Kp – 1 κi
This equation relates the moment of inertia of the primary pivot with the moment of iner-
124
The limits on the other flexible segments can be determined quickly by using the
width and thickness ratios. The equation for determining the thickness limitation is
c i ≤ Cc p (8.8)
The minimum thickness of any of the other flexible segments is found by using
b i ≥ D min b p (8.9)
If an out-of-plane orientation is used, then the width of the other flexible segments
is found to be
b i = D equal b p (8.10)
Equation (8.7) can be used to quickly calculate the geometry of the other flexible
segments in the mechanism. The ratio values and Equations (8.8) and (8.9) can be used to
verify the geometries of the other flexible segments if neither an out-of-plane (equal thick-
nesses) or in-plane (equal widths) orientation is used. However, if either one of these ori-
entations, including the perhaps out-of-plane orientation are used, then the geometries of
the flexible segments do not need to be verified. They are guaranteed to have acceptable
The final step in the design is to determine the length of the flexible and rigid seg-
ments of each link. The following steps, along with the definitions in Figure 8.3, can be
125
f2
f1
f3
slider
(rigid segment)
r2
F
r3
r tot
r i = ------- (8.11)
ζ
and for
r2 … ζ = R+1 (8.12)
1
r3 … ζ = --- + 1 (8.13)
R
r ave … ζ = 2 (8.14)
4. Use the relations in Table 8.1 to calculate the lengths of the rigid sections
126
8.2 Feasible Configurations
mechanism to be chosen. The following two methods can be used to identify mechanisms
that are suitable for a given design problem based on either stress or force feasibility.
These methods utilize the classification variables and design constraints to eliminate
mechanisms that are infeasible, leaving only those mechanisms that can provide the
Table 8.1 Needed values to calculate flexible and rigid segment lengths
1-
lpp --------- 0 0 r2 0
0.85
psp 0 0.1 0 f2/2 f2/2
1-
plp 0 --------- (times r2) 0 r2 0.15*f2
0.85
ssp 0.1 0.1 0 f1/2+f2/2 f2/2
1-
slp 0.1 --------- (times r3) 0 f1/2+0.15*f2 r3
0.85
sps 0.1 0 0.1 f1/2 f3/2
1-
lps --------- 0 0.1 r2 f3/2
0.85
sss 0.1 0.1 0.1 f1/2+f2/2 f2/2+f3/2
127
Since all 3 principal equations contain classification variables, some approxima-
tions must be made. For this reason, the results of these methods are only as accurate as
This method uses the α curve fit parameter M and the stress design equation to
identify stress feasible mechanisms. The mechanisms identified are capable of producing
the desired deflection without violating any of the stress and geometry constraints.
Although these mechanisms are stress feasibly, they are not guaranteed to be able to pro-
There are two types of stress feasible mechanisms. The first type, guaranteed stress
feasible mechanisms, are those mechanisms with an M value less than Mmin, the minimum
M value calculated for the worst of the design constraints. These mechanisms are guaran-
The second set of mechanisms are simply referred to as the stress feasible mecha-
nisms and can be identified using the maximum value of M, Mmax, calculated using the
conservative constraints. These mechanisms with M values less than Mmax, but greater
than Mmin, are stress feasible for only part of the constrained design space. At certain
points in the constrained design space, these mechanisms will not be able to achieve the
128
Stress Infeasible
Mechanism M value
Mmax
Stress Feasible
Mmin
Guaranteed
Stress Feasible
Figure 8.4 Summary of M value requirements for different types of identified mechanisms
Any mechanisms that have an M value greater than Mmax can not obtain the
desired deflection without violating the design constraints. A graphical summary of the
required values for M for each type of mechanism is shown in Figure 8.6.
In most cases, it is most beneficial to identify and use the guaranteed stress feasible
mechanisms. However, at times, guaranteed stress feasible mechanisms may not exist or it
• All of the isolated and coupler variables in the stress design equation
must have known or constraint values.
129
1. Solve the displacement equation. If a specific type of configuration is desired,
use the highest λ value for Mmin and lowest value for Mmax. Otherwise, use
2. Calculate the desired M value. Calculate the M value that corresponds with the
3. Eliminate mechanisms. Any sub-classes with a value for M larger than the M
value calculated can be eliminated from the group of mechanisms being identi-
fied. Table D.3 in Appendix D contains all of the mechanisms and their param-
eters sorted by M.
Mechanisms that can undergo deflection without Mechanisms that can produce the desired force and
Identifies
violating design constraints deflection without violating design constraints
130
4. Choose mechanism and verify. It should be noted that those mechanisms that
have an M value close to Mmin and Mmax may or may not be suitable due to the
actual n, M, and λ values should, together with the stress design equation, be
5. Finish the design. Once mechanisms have been eliminated, a mechanism must
This method uses the corresponding percent deflection lpp configuration and the
parameter Ψ to identify mechanisms that are force feasible, those mechanisms capable of
delivering the required force without violating the stress, width, and length design con-
straints. Even though these mechanisms can produce the desired force, due to the nature of
the method, it can not be guaranteed that the thicknesses of the flexible segments will fall
• All of the isolated and coupler variables in the stress and force design
equations must have known or constraint values, except for h.
131
1. Calculate maximum h. The stress design equation can be used to calculate the
maximum h value for the lpp-a or lpp-b mechanisms. Use the classification
2. Calculate maximum force. Calculate the maximum force possible for the lpp-a
or lpp-b sub-classes using the maximum h value and the other variables in the
force design equation. Use the classification variable values that correspond to
3. Divide forces. Divide the needed force by the maximum force found in step 2.
This gives Ψ min , the minimum Ψ value that can be used to achieve this force.
used in step 1 and a Ψ value greater than Ψmin is capable of producing the
required force without violating design constraints given the same length,
stress, and flexible segment widths used in step 1. All other mechanisms can be
5. Determine thickness. The actual thickness, which will be different than the
Once a mechanism is chosen, there is no need to re-solve the stress design equation
as long as all parameters, except the thickness, remain the same. However, if any con-
132
Guaranteed
Guaranteed Stress and Force
Stress Force Feasible
Feasible Feasible Method 2
Method 1
Methods 1&2
straint values are modified, it becomes necessary to re-check all of the principal equations.
Table 8.2 contains a summary of force feasible method, including the variable values
Either of the two methods can be used to identify feasible mechanisms for a given
problem depending upon the objective and needs of the design problem. In some cases, it
is advantageous to use both methods on the same problem. This, in a sense, overlaps the
two types of feasible regions to generate a region in which both types of feasibilities exist,
as graphically demonstrated in Figure 8.5. The mechanisms within this region are those
mechanisms that have been identified as guaranteed stress and force feasible. All of these
mechanisms will satisfy all of the constraints in all circumstances (assuming approxima-
tions are close). This makes the selection process very simple, allowing the designer to
133
8.3 Other Mechanism Considerations
The methods discussed in Section 8.2 identify those mechanisms that are stress
and force feasible. This section discusses other aspects of CFMs that must be considered
The percent constant-force must be considered when defining the design problem
and choosing a mechanism. When selecting a mechanism, it is generally best to choose the
lowest percent constant-force value possible. These mechanisms tend to have a higher
stiffness than the mechanisms with a higher percent constant-force allowing them to
achieve higher forces for the same relative stress. If a minimum percent constant-force
The requirements for the flexible segment configuration can also be used to elimi-
nate infeasible mechanisms. If the design problem consists of constraints that limit the use
of certain types of flexible pivots or pin-joints, then those configurations that are not suit-
8.3.3 Manufacturing
The types of manufacturing processes that can be used to make a given mechanism
134
8.3.4 Normal Displacement
or verify a design if the design problem contains normal displacement constraints. The
normal displacement can be calculated quickly using the equations developed in Chapter
5. Solving Equation (5.37) for ∆y and using d Nmax and the length parameter equation
results in
d Nmax l tot
∆y max = ---------------------
- (8.15)
100λ
Equation (8.15) can be used to calculate the maximum normal deflection when the
mechanism is deflected the maximum deflection percentage for the sub-class (16 for sub-
class a, 40 for sub-class b). For deflection percentages below the maximum value, Equa-
For design purposes, the normal displacement information can be used to deter-
mine the upper limits on the length and/or deflection percentage when a maximum normal
displacement is specified. The normal displacement equation can be used at any time dur-
ing the basic design. The location of their use in the basic design steps depends on the type
During the iterative process found at times in CFM design, it is valuable to have an
made.
135
Variable Increased E Sy SF b h ltot disp d M Ψ β λ
3 2 - 2
Effect on Force -
L L L
Effect on Stress - - -
L L
Key
n n Decrease, Increase in Magnitude
by power n (default is 1)
Figure 8.6 Summary of variable effects on force and stress magnitudes for an increase in variable
magnitudes
The variables in Equations (8.2), (8.5), and (8.6) affect the force and/or stress mag-
nitudes of a CFM. Some of these variables have a direct effect on the magnitude of the
force or stress, while others affect the limits of the force or stress magnitude. The effect of
an increase in magnitude of these variables on force and stress are summarized in Figure
8.6.
This figure can be used to determine which directions the variables should move to
affect the stress, force, or other variables. Figure 8.6 can also be used to decide whether
the lower or upper limit of a constraint value should be used in the principal equations.
136
8.4.2 General Guidelines
Several general guidelines can be established from the trends illustrated in Figure
8.6. These guidelines provide insights into the response of CFMs to changes in the vari-
ables.
1. Adjust width to change force. For deflection loads, the width increases the
2. Length verses thickness. Length and thickness both have a 1:1 ratio with the
stress design equation. However, any changes in thickness are cubed and any
changes in length are squared in the force design equation. Therefore, the fol-
3. Maximize the length of the mechanism. Most design problems tend to have
have more difficulty satisfying the stress design equation rather than the force
4. Minimize safety factor. Use the smallest possible value rather for the safety fac-
5. Increase the yield strength. Increasing the yield strength is not always an
option, but it is very useful. Changes in yield strength affect only the limit on
137
6. Change the type of mechanism to affect force. By moving to a mechanism with
a higher Ψ , a larger force can be achieved for the same given stress level.
However, the flexible segments thicknesses will adjust and the usefulness of
this method depends upon the constraints of the system. Generally, increases in
In some design cases, the force may be too high, and manufacturing limits on the
widths or thicknesses may be preventing the force from being lowered. In this case, most
of the rules listed above can be reversed to help decrease the magnitude of the force.
A steel Class 1A-lpp mechanism is needed that produces 25 lbs of force and dis-
places 4 inches. The mechanism should also be as high of level of constant-force possible.
The length must be between 12 and 15 inches and the width must be between 3 and 5
inches. The material comes in sizes ranging from 0.04 to 0.07 inches thick with a range of
50 to 200 Kpsi yield strength, with the material becoming more expensive as yield
strength increases. Even though cost needs to be minimized, a safety factor of at least 1.2
is required. Design a mechanism that satisfies the criteria summarized in Table 8.3.
1. Choose a mechanism
138
The configuration of the mechanism has been specified, but the sub-class has not.
However, Table D.1 indicates that sub-classes a and b have the highest percent constant-
force for the lpp configuration. Additionally, solving the deflection equation for the maxi-
mum and minimum mechanism lengths gives the high and low values of d as 36.6 and
29.3. This leads to the conclusion that the lpp-b mechanism must be used.
This example has four known values for the three principal equations. Examina-
tion of Figure 8.1 indicates that more variable values are needed.
The constraint values are chosen based on the objective of the problem. In this
example, the cost (yield strength) needs to be minimized. According to Figure 8.6, the
required yield strength can be decreased by decreasing the stress. Figure 8.6 also indicates
that increasing the length and decreasing the thickness will decrease the stress. Therefore,
the maximum length, and the minimum thickness should be used in the principal equa-
139
tions. Using the maximum length results in a value for d of 29.3. This value can also be
Examination of the force design equation shows that it can now be solved for b
without adding additional constraint values. To achieve the desired force, b must be 8.99
inches, which exceeds the upper constraint. This now indicates that b must be a constraint
value and that one of the previous constraint values must be removed from the system.
According to the Figure 8.6 and the general guidelines, h should be used rather than l tot to
change the force while minimizing the affect on stress (yield strength).
By removing the constraint value for h and adding the maximum b value, the force
design equation can now be solved for the need h to achieve the desired force. This results
in a value for h of 0.049 inches. The final step is to resolve the stress design equation for
the needed value of Sy for the given geometries. This results in a final value for Sy of 167
Kpsi. This is the smallest yield strength that can be used to achieve the desired force using
140
In this example, a suitable design, summarized in Table 8.4, was defined. By using
Figures 8.1 and 8.6, as well as the general design steps and guidelines presented in this
chapter, a suitable design was quickly discovered without the need for large numbers of
Bronze. A force of 3 mN is required from a mechanism that has a length of 150 mm and a
displacement of 7.5 mm. The thicknesses of the flexible segments must be between 1.5
and 2.0 mm thick. It is also necessary that the mechanism is at least 90 percent constant
and has a safety factor of 1.1 or greater. A summary of all of the requirements can be
Question: Is it possible to use a fully compliant configuration for the above design
problem?
141
The first basic design step is to choose a configuration. In this particular problem,
mechanisms as possible using one of the two methods described in Section 8.2.
The deflection equation is missing one of the coupler variables and the classifica-
tion variable, which satisfies the first requirement for both methods. Examination of the
stress design equation shows that all of the isolated and coupler variables have either
known or constraint values, while the force design equation is missing values for the vari-
able b. This indicates that only method 1 can be used for this problem.
In this case, we are looking for a fully compliant mechanism. All of these mecha-
nisms have a λ value of 1.05. Using this value and Equation (8.6) results in
142
Since the only objective is to determine if a fully compliant CFM can be used, only
the stress feasible mechanisms need to be identified using method 1 and Mmax. This value
can be found by using a form of Equation (8.5) and the appropriate values given above.
Examination of Table D.3 shows that the smallest M value for any sss configura-
tion that is suitable for an in-plane orientation is the sss-a95I mechanism with an M value
of 1.7796, well above the maximum M value calculated above. This leads to the conclu-
sion that the fully compliant mechanism can not be used for these design constraints.
Question: What length of mechanism is needed to be able to use the fully compli-
1. Choose a mechanism
The sss configuration with he lowest M value was identified above. This mecha-
nism also has the highest Ψ value of all of the sss-a mechanisms. Therefore, this mecha-
A non-linear method must be used to calculate the needed length due to the inter-
dependencies between the stress design equation and displacement equation. Performing
the calculation with the appropriate values for M (1.7796) and n (.4888) for the sss-a95I
143
mechanism, as well as the other known values results in a length of 418.4 mm. This is
1. Choose a mechanism
The sss-a95I mechanism will be used for the same reasons listed above.
Using the appropriate values for M (1.7796) and n (.4888) and Equation (8.5), the
This leads to the conclusion that the sss-a95I mechanism can be used for the given
length and deflection if the primary pivots thickness is no larger than 0.326 mm.
Question: Assuming that the thickness above is suitable, what width of material
1. Calculate the required width to achieve the force using the force design equa-
The equation is
2 κ K
12Fl tot 1 p–1 12 ( 0.003 )0.15 2 ( 26.57 )15.3423
bp = ------------------------------------
- = --------------------------------------------------------------------------------------------------------
- = 2.15 mm (8.19)
λ 2 κ p h p3 βE 1.05 2 ( 16.03 ) ( 0.000326 3 )2.279 ( 1.1 × 10 11 )
144
Equation (8.19) indicates that the width of the flexible segments must be 2.15 mm
wide.
Question: What are the required thicknesses for the other flexible segments?
Equation (8.7), together with the equation for the moment of inertia for a rectangu-
lar cross section can be used to calculate the thicknesses of the other beams. These thick-
ness are 0.111 mm and 0.091 mm for the first and second flexible segments respectively.
The only thing that has not been verified are the stresses in the other flexible seg-
ments. However, this is not necessary. The highest stress in the mechanism is found in the
primary pivot assuming that the other flexible segments satisfy the inequality in Equation
(8.8), and because all of the flexible segments have the same width and this mechanism is
A CFM is needed that produces 100 lbs of force over 6 inches of deflection and
with 10% or less variation in the force magnitude. The mechanism must be exactly 2 feet
long and 5 inches wide. The mechanism is to be made out of spring steel with a yield
strength of 200 Kpsi and a thickness between 0.01 and 0.035 inches. Identify all of the
configurations that satisfy all of the design requirements as summarized in Table 8.6.
1. Choose a mechanism.
145
The mechanism was not identified in the problem. However, the value for d, using
must be used. To identify all of the possible mechanisms, both the stress and force feasi-
Using the design constraints and the stress design equation, a maximum value for
the thickness of the lpp-b sub-class can be calculated. This turns out to be 0.116 inches.
This value can now be used in the force design equation to determine the maximum force
the lpp-b configuration can generate given the design constraints. The maximum force is
130.3 lbs.
Dividing the needed force of 100 lbs by the maximum force of 130.3 results in a
Ψ min of 0.767. Any mechanism with a value of Ψ greater than Ψ min will be able to gen-
erate the needed force without violating the design constraints (assuming the same width
146
The application of this method reduces the number of sub-class b mechanisms
Using Equation (8.5), and the maximum thickness constraint together with the
other constraints results in a M min value of 1.83. Using this value and Table D.3 in
identified also using method 2. Therefore, any one of these 5 mechanisms (lps-b90I, lpp-
b90, lps-b95I, lpp-b95, and lpp-b) can be used to satisfy the design requirements of this
example.
147
148
CHAPTER 9 CONSTANT-FORCE
ELECTRICAL CONTACT
The objective of this chapter is to demonstrate the potential benefits and viability
cessful development of a CFEC that meets all of the requirements of an electrical contact
will lay a ground work for further exploration and introduction of CFEC’s into industry
applications.
The chapter begins by discussing current electrical contact industry practices and
configurations are explored to discover one suitable for application as a CFEC in a pre-
sented case study. (The CFEC case study presented in this chapter is focused on limita-
tions and common industry practices associated with Personal Digital Assistants (PDA)
docks, but the principles are applicable to a wide range of connector applications.) The
chapter finishes by describing the modeling, optimization, and verification of a CFEC that
149
9.1 Introduction
and methods to improve this reliability are always being evaluated. According to Desh-
pande and Subbarayan (2000), the reliability of high-cycle electrical connectors is related
to electrical signal propagation, and mechanical performance and stability. To achieve this
reliability in practice, the connector contacts must transmit the electrical signal with mini-
mal contact resistance under all types of use conditions and accommodate expected geo-
The factor that contributes most significantly to the reliability of electrical contacts
is the contact-surface mating conditions. Two physical parameters that greatly affect mat-
ing conditions are contact surface finish, and contact normal force at mating. When con-
tact surface finish remains corrosion free, either by being corrosion resistant or by being
above a certain level, greater reliability is also achieved (Harper, 1997). Contact normal
forces must be small enough to minimize plating damage over the life of the contact, yet
large enough to overcome co-planarity differences from adjacent contacts and other geo-
metric variations. Thus a desirable contact system would maintain an optimal contact
being driven to produce innovative products that have faster speeds, smaller packages, and
150
higher contact density. To remain competitive, performance gains must be achieved with
mechanisms, have created the potential for small-scale, low-cost, constant-force electrical
contacts (CFECs). CFECs differ from traditional contacts and springs by the separation or
describe force (F) and deflection (d) of springs as = kd where k is the spring constant
A CFEC uses constant-force technology to separate the contact normal force and
deflection. By removing the traditional constraints imposed by forces and deflections that
are dependent on each other, the addition of new types of electrical connectors previously
The disassociation of contact normal force and contact deflection could lead to
three advantages. First, current electrical contacts require tight manufacturing and assem-
bly tolerances to ensure that the contact deflection is within acceptable limits. The decou-
pling of the force and deflection may allow the tolerances to be loosened, while still
having acceptable performance. This can help to reduce the cost and difficulty of manu-
The second advantage may be in applications were the user interacts with the elec-
trical contacts, such as docking stations. In this case, the decoupling of the force and
151
Undeflected Contact Deflected Contact
Undeflected Contact
Housing Housing
Springs
Solder Point Deflected Contact
(a) (b)
Figure 9.1 (a) Pogo type connector and (b) cantilever type connector
deflection helps prevent any variation introduced by the user, such as different docking
The third advantage that may be seen from decoupling is in applications where the
output force needs to remain relatively constant to ensure performance, but the deflection
does not remain constant due to movement and/or vibrations of the contacts due to the sys-
tem environment. Examples of this include connectors in aircraft, vehicles, and machin-
ery.
Traditionally, electrical contacts for use in PDA docks have consisted of linear
spring assemblies or an arrangement of cantilever beams (See Figure 9.1). These configu-
rations usually require large deflections to obtain the desired force and must be long to
The electrical contact industry has several practices and standards that constitute
essential performance characteristics for electrical contacts. The most basic and important
152
of these are divided into subgroups for presentation here, but it is not intended to be an
The first standard is that electrical continuity must be continuous. The electrical
path created by an electrical contact can not be interrupted or contain high resistive areas.
Additionally, this path should consist of few parts (preferably one piece) which are easily
assembled.
Electrical contacts can be fabricated from any conductive material, but current
industry practice is to use alloys that contain copper. Phosphor bronze is a common alloy
that is easy to use and readily available. Beryllium copper and titanium copper are com-
monly used to achieve higher yield strengths. Unfortunately, they are more difficult to use
contacts are being produced in ever increasing volumes at lower costs. Current industry
practice is to use progressive stamping techniques to shape the metallic beams. Generally,
the contacts are stamped at the desired pitch distance and are left attached to the flashing.
This allows for easier material handling and assembly, but limits the shape and design of
the contact. Some of the limitations imposed on designs due to this manufacturing process
are:
153
Many industry practices are associated with assembly of the electrical contacts. In
this case, assembly deals with packaging the electrical contacts into a housing that is suit-
able for use. Design for assembly is vital to achieve a low cost and reliable part. Some of
• Mount type - Through-hole mounts are generally easier to use and most
commonly used. Surface mounts usually are more difficult to do, but
take up less room on the board (one side only).
• Housing - The plastic housing that holds the contacts generally is one
or two pieces. This part holds the contacts into place.
The design of electrical contacts is well defined and understood. However, exami-
nation of traditional contacts shows that they are not suitable for use as constant-force
Traditional CFM configurations such as the one shown in Figure 9.2 are not suit-
able for use as electrical contacts for several different reasons which include:
154
Displacement
Evaluation of the latest configurations shows that pin-joints and small-length flex-
ural pivot can not be used in electrical contacts, indicating that different configurations
that combine the benefits of both electrical contacts and CFMs must be developed for use
as a CFEC.
tional electrical contacts and current constant-force mechanisms are not acceptable for use
as electrical contacts, they provide a starting point in the search for new configurations
155
Displacement
The slider-crank constraints can be greatly simplified by using a method that simu-
lates the function and motion of a fixed-pinned flexible segment and a rigid link joined by
a pin joint. In this method a circular cam is used to represent the rigid-link. If the simu-
lated joint remains in compression, then the flexible link will follow the cam profile - the
exact path of the replaced rigid body link - as shown in Figure 9.3.
However, there are limitations to the simulated pin joint method. It must be in
compression to ensure that the tip of the beam remains in contact with the cam (See Figure
9.4c). It is also important to ensure that the simulated pin joint has minimal friction so that
there is smooth motion around the cam. This can be partially accomplished by rounding
over the tip of the beam and providing smooth surface finishes on the cam as shown in
Figure 9.4d.
156
If the flexible beam is loaded, but doesn’t slide around the cam, then the beam
could buckle. Slider crank change points should be avoided. At these points, it may be dif-
ficult to get the beam to begin to slide around the cam. This can be done by changing the
initial angle of the beam or changing the eccentricity of the slider-crank as illustrated in
There is also a limit to how far around the circular cam the flexible link can travel.
If the mechanism is deflected beyond the point at which the flexible beam is tangent to the
cam, the tip of the flexible beam will no longer be in contact with the cam (See Fig. 9.4e).
Figure 9.4 shows a graphical summary of the limitations, along with the methods
to overcome these limitations, associated with the simulated pin joint method. Despite
these exceptions, the simulation of pin joints with the use of a circular cam is an important
tool.
contacts, many different possible configurations were evaluated. Using the industry prac-
tice and standards criteria and a screening process, the configuration determined most via-
ble for use in a CFEC is one in which a slider-crank mechanism is attached directly to the
end of a bent cantilever beam as illustrated in Figure 9.5. A concept drawing of the CFEC
This configuration is easy to manufacture and assemble, and has electrical continu-
ity. Additionally, the beam and cam combination provide the necessary increases in
157
Displacement Displacement
Cam Cam
Figure 9.4 Limitations and solutions to limitations of simulated pin joint method
Deflection
Cam
158
n
otio t
kM Un i
Doc kab
le
Do c
Docking Station
mechanical advantage and the strain energy storage device necessary for constant-force
behavior.
Parameters establish the shape and size of the mechanism and are used as inputs in
the model and optimization. Among these parameters are link lengths, angles, and cross
sectional geometries. Some of the parameters are also used to assess performance relative
to design requirements associated with the industry practice and standards. A graphical
summary of each of the important parameters for the selected configuration is shown in
Figure 9.7.
159
Deflection
δ
R1
h Cross Section
l3
l2
h package b
θ3
l1 θ2 R3
l4
l5
R2
θ1
R4
l6
b package
CFEC is demonstrated by using a specific case study of electrical contacts for Personal
Digital Assistants (PDAs) docking stations. Specific requirements for the case study were
gathered based on existing dock designs and by working closely with the engineers of a
for this contact because it is relatively common, cheap, and easy to work with during pro-
duction when compared to other alternatives. The force range for the design is 294 - 588
mN (30 - 60 gf), and the maximum stress in the contact should not exceed the yield
strength. The CFEC is required to have a cross sectional height, h, of 0.2 mm and a width,
b, of 1 mm. The case study also requires that the contact fit inside of a 12 mm wide by 6
mm tall rectangle. The final design constraint is that the output force of the mechanism
160
should be at least 60% constant. This means that there is only a 40% variation between the
minimum force and the maximum force within the possible displacement range. Table 9.1
summarizes the design constraints for the case study. The symbols for the model functions
for each constraint are listed in the first column. The second column contains the general
constraint symbol which represents the fixed design constraints for any problems, while
An accurate model in which the governing parameters can be modified and accu-
rate resulting forces and displacements can be calculated is needed. During the optimiza-
General
Constraint Value
Model Function Symbol Constraint
for Case-Study
Symbol
hpackage hpc ≤ 6 mm
bpackage bpc ≤ 12 mm
h hc 0.2 mm
b bc 1.0 mm
E Ec 110e9 Pa
Sy Syc 552e6 Pa
SF SFc ≥ 1.0
161
tion phase, the model is called many times to calculate function values and derivatives.
numerical bent beam model based on Euler’s Method were combined to model the CFEC
configuration. However, the model was not accurate enough due to violations of assump-
tions in the CFM model. The CFM portion of the CFEC does not act as a pure slider. The
motion is not straight line and the moment assumed not to pass through the slider, is actu-
ally passed to the bent beam. These differences proved too much for the joint numerical
model.
To overcome the model problem, a finite element analysis (FEA) program capable
of nonlinear analysis (ANSYS) was used to model the deflections, contact forces, and
stresses in the CFEC. A parametric model was used so that values could be passed
The FEA model was generated using the input parameters to calculate the location
of the key points shown in Figure 9.8. Once all the key points have been defined, a total of
The cam is replaced with the rigid link (segment A) that it is simulating. This
requires that segment A be pinned to ground at key point 1, and that key points 2 and 3 be
constrained to have the same x and y displacement, thus forming a pin-joint. Segment A is
162
100
5
4
20
D
B
6 21
F 8
2 3 G
E
7
22 9
A
H
1 23 10
I
J 11
12
Figure 9.8 Key Points for the finite element model
as well as rotation about the z -axis. This represents the way the bent beam attaches to
ground as a cantilever, simulating its attachment when soldered to a printed circuit board
(PCB). Finally, 5 vertical displacement load steps in the downward direction are applied to
Once the model has run for the 5 different load steps, the contact force for each
load step and the highest stresses over the total deflection are written to a data file for use
ables, design functions, and constraints that facilitates the development of a constant-force
163
mechanism from the layout presented in Figure 9.5 that satisfies all of the design con-
straints.
As with any CFM, the principle objective function is the parameter Ξ' . The
lengths, angles, and radii described in Figure 9.7 are established as design variables in the
optimization problem with reasonably assumed bounds. The beam height (h), width (b),
modulus of elasticity (E), and safety factor (SF) are set up as analytical variables. The val-
ues for the analytical variables are established by the requirements of the case study as
described in Table 9.1. The remaining constraints of the case study show up in the design
9.7.1 Optimization
Using the design constraints of the case study, the optimization problem can be
SF ≥ SFc (9.5)
164
With the following constraints on the design variables:
R 1 > Rc (9.8)
R 2 > Rc (9.9)
R 3 > Rc (9.10)
R4 > Rc (9.11)
Where the variables with a subscript c denote the constraint values found in Table 9.1.
The optimization and FEA model were linked together in a similar manner as the
CFM model. At first, a feasible starting point was difficult to find so the constraints were
loosened by 10%-15%. Once a starting point was found, the optimization was allowed to
run and the constraints were tightened. This was repeated until a suitable design was
found.
The final design chosen satisfied the design constraints and requirements of the
case study. A detailed drawing of the final design chosen for the case study is shown in
Figure 9.9. The model values for the design and constraint parameters are listed in Table
9.2.
To confirm the behavior of the CFEC and the accuracy of the model, 9 prototypes
of the final design were produced for testing. The photo in Figure 9.10 illustrates the com-
165
4.7672
2.9398
x2
78
.92
R0
[9
5°
]
[10°]
.7
R0
4.5521
.7
R0
.5
2.04
R0
R0
.7
1.0
0.4184
0.4956
1.35
0.5087
2.72
5.1221
5.3912
Constraint Value
Model Function Symbol Final Design
for Case-Study
hpackage 5.9 mm ≤ 6 mm
bpackage 5.4 mm ≤ 12 mm
h 0.2 mm 0.2 mm
b 1.0 mm 1.0 mm
Sy 552e6 Pa 552e6 Pa
SF 1.29 ≥ 1.0
Fave 478 mN (48.8 gf) ≈ 441 mN (45 gf)
Fmin 423 mN (43.2 gf) ≥ 294 mN (30 gf)
166
Figure 9.10 CFEC prototype as compared to a dime
dimensions were to the specified dimensions. An optical comparator was used to take 16
which of the prototypes were closest to the final design. The three prototypes closest to the
final design were chosen for testing. The results of the dimensional analysis for the three
9.8.2 Testing
A rigid test fixture was designed to allow for easy and accurate placement of the
prototype. The cam was fabricated as a separate piece to help ensure tight tolerances and
allow for different materials for the cam to be used, including polypropylene and teflon.
167
The purpose of the different materials was to investigate how different material types
controlled the actuator and collected position and force data. During testing, the contacts
were deflected to 0.75 mm and back. Figure 9.11a shows a photo of the general test setup
and Figure 9.11b shows a close-up photo of the contact in the test fixture.
9.9 Results
The prototypes were each tested using two different cams. Figure 9.12 shows a
graph of testing results for prototype 3 with a polypropylene cam. The mechanism main-
168
(a) (b)
Figure 9.11 (a) General testing setup and (b) close-up of contact with fixture and probe
1000
800
600
Force (mN)
Compression
400
Expansion
200
0
0.000 0.100 0.200 0.300 0.400 0.500 0.600 0.700 0.800 0.900
-200
Displacement (mm )
There are two interesting phenomenon that were observed in every test. First, since
there is no pre-load on the mechanism, the initial force must be zero. However, as the
mechanism goes through the initial displacement (about 0.05 mm), there is a sharp rise
from zero force to the intended constant-force. The phenomenon was observed in every
test and in fact was observed by Millar et al. (1996) during initial testing of constant-force
169
The second phenomenon observed is a difference in force between the compres-
sion and expansion strokes of the testing. During the compression stroke, the mechanism
experienced a higher force than predicted. As the mechanism reverses direction, there is a
sharp decrease in the force to a point below the predicted force which persists throughout
This phenomenon was consistent throughout testing and is found in all types of
electrical contacts and mechanisms. In fact, this same phenomenon was also observed by
Boyle (2001) while studying the dynamics of constant-force mechanisms. In all cases, this
behavior is consistent with the effects of friction, which acts in the direction to oppose
motion. Boyle (2001) was successful at modeling this phenomenon as friction found
The accuracy of the model can be verified by comparing the testing results with
the predicted results. However, the model used for the case study does not account for the
friction, requiring that the effects of the friction be removed from the test data. Assuming
that the difference in force between the compression and expansion strokes and the force
of the mechanism without friction is the magnitude of the friction force, the effects of fric-
Additionally, to make the comparison between test results and model predictions,
new predictions were made based on the actual shape and size of prototype 3. These pre-
dictions are listed in Table 9.4 while Figure 9.13 shows the predicted and adjusted-mea-
170
600
500
200
100
0
0 0.2 0.4 0.6 0.8
Deflection (mm)
Figure 9.13 Average and predicted force comparison
Constraint Value
Function Symbol Prototype 3
for Case-Study
hpackage 5.6 mm ≤ 6 mm
bpackage 5.5 mm ≤ 12 mm
h 0.2 mm 0.2 mm
b 1.0 mm 1.0 mm
Sy 552e6 Pa 552e6 Pa
SF 1.29 ≥ 1.0
Fave 448 mN (45.7 gf) ≈ 441 mN (45 gf)
Fmin 418 mN (42.6 gf) ≥ 294 mN (30 gf)
171
The percent constant-force for prototype 3 can be calculated by using Equation
(9.1). However, since the model does not include the region of quickly rising forces in the
initial deflections, the method introduced in Section 5.3 must be used to calculate the
extrapolated percent constant-force, Ξ' ex . This value is found by using Equation (5.27)
with the lowest average force on the flat part of the curve as Fmin and the highest force as
F max . This results in a value for Ξ' ex that is within 12% of the predicted extrapolated per-
cent-constant force. Table 9.5 contains a summary of the comparison between the testing
for use in electrical connectors. Viable configurations were developed, and one configura-
tion was chosen for use in a case study. A design for the chosen configuration was gener-
ated, prototyped, and tested. The testing results indicated that the model predicted the
could provide a number of benefits in terms of performance, robustness, and package size.
172
The successful demonstration of the uncoupling of the force and deflection in an electrical
introduced by the user, and increased system robustness in applications where movement
Further work to explore the contribution of CFECs to these areas must be done.
The size and force limitations and the affects of tolerances on CFECs should be better
173
174
CHAPTER 10 SUMMARY,
CONCLUSIONS, AND
RECOMMENDATIONS
10.2 Summary
The purposes listed above were attained through the main contributions of this
research including:
175
• The classification system for CFMs was refined and expanded to
include naming methods for groups of mechanisms, families of mecha-
nisms, and individual mechanisms.
• The model was used, together with the comparison methods, to seek
out and identify new and improved mechanisms. This effort resulted in
the discovery of mechanisms with major improvements in stiffness and
feasible manufacturing orientations.
10.3 Conclusions
Before this research was performed, the CFM knowledge base was limited. This
work has resulted in new understanding and insights into CFMs. Important new compari-
son methods were developed that can be used to identify new mechanisms and aid in the
design process.
While it was not known whether or not improvements over the original CFMs
could be made, this work showed that improvements were possible. In most cases, the
level of constant-force could not be improved. However, for the Class 2B mechanisms, the
level of constant-force was improved from 93 percent constant to 99 percent constant for
176
the sub-class a mechanisms. For the sub-class b mechanisms in this class, the percent con-
For every class of mechanism, an increase in the stiffness was made without
increasing the relative stress in the mechanism. Each class of mechanisms had some kind
of increase (ranging from 5% to 3000%) with the largest improvements in Class 2B.
Four 16% deflection mechanisms and five 40% deflection mechanisms were
defined with a stiffness parameter greater than or equal to 1 (see Table 10.1). These mech-
anisms were in the lpp and lps configurations and are the stiffest known CFMs.
It was also observed that the configurations with at least one long flexible segment
had higher stiffness parameters than their counterparts with only small-length flexural piv-
ots. The use of long flexible beams is important in gaining stiffness without increasing
stress.
This work also showed that a reasonable design approach is possible. This design
approach was defined and works for many different design scenarios. The design method-
177
ology allows the designer to make educated decisions and choices based on the informa-
tion summarized in the tables. Additionally, two methods were outlined that identify
mechanisms that are either feasible from a stress or force standpoint. It was also shown
that when these methods are used together, it is possible to identify all of the mechanisms
It can be concluded from this work that the methods and techniques used to design
CFMs can be used to develop models for mechanisms other than the traditional slider
crank. The CFEC developed in this work is such a mechanism. A number of viable config-
urations were developed, and one configuration was chosen. From this configuration, a
design was generated, prototyped, and tested. Testing results showed that the CFEC dis-
sensitivity to variations introduced by the user, and increased system robustness in appli-
cations where movement and/or vibrations exist. This could lead to benefits in terms of
This section recommends and outlines some of the areas in which future work
should be undertaken.
178
10.4.1 New Configurations
Much work was done to understand the limitations of the original configurations.
However, much work is needed to understand the following types or variations of config-
urations.
• lpl- Configuration with two long flexible segments attached with a pin
The model developed for this work had the ability to add pre-loads, linear springs,
and offsets to the traditional mechanisms, as well as consider orthoplaner type CFMs. It is
recommended that further work be done to fully understand the effects and possibilities
179
10.4.2 Physical Implementation
the development of the CFEC, several different phenomenon where observed that affected
• Sensitivity to tolerances
These phenomenon and others should be studied, allowing CFMs to be used in future
applications.
work. With the development of the CFEC, a different type of mechanism was developed
to examine other types of mechanisms for the properties required for constant-force
behavior. Some of these mechanisms include: centrifugal clutches and four-bar mecha-
nisms.
The lpp configuration exhibited the highest stiffness of the configurations exam-
ined. Other configurations were constrained by stress due to small-length flexural pivots.
It can then be reasoned that the lpl configuration, which adds a second long flexible seg-
180
ment without affecting stress, should exhibit an even greater stiffness than the lpp config-
uration. It is recommended that further work should be done to understand the behavior of
10.4.5 Applications
The improved mechanisms and the design methodology developed can be used to
develop CFMs for applications. It is recommended that work be done to implement CFMs
in viable commercial applications. This can be in conjunction with the work of Boyle
robustness must be done. The size and force limitations and the affects of tolerances on
CFECs should be better understood, as well as a better understanding of the effects of the
181
182
REFERENCES
Bossert, D., Ly, U. L., and Vagners, J., 1996, “Experimental Evaluation of a
Hybrid Position and Force Surface Following Algorithm for Unknown
Surfaces,” Proceedings - IEEE International Conference on Robotics asnd
Automation, Vol. 3, pp. 2252-2257.
Brush Wellman Inc., 1999, Connector Engineering Design Guide: Material Selec-
tion in the Design of Spring Contacts and Interconnections, Brush Well-
man Inc.
Chang, L. H., and Fu, L. C., 1997, “Nonlinear adaptive control of a flexible manip-
ulator for automated deburring,” Proceedings - IEEE International Confer-
ence on Robotics asnd Automation, Vol. 4, pp. 2844-2849.
Herder, J. L., van den Berg, F. P. A., “Statically Balanced Compliant Mechanisms
(SBCM’s), An Example and Prospects,” Proceedings of the 2000 Design
Engineering Technical Conferences and Computer in Engineering Confer-
ence, DETC2000/MECH-14144.
183
Howell, L. L., 2001, Compliant Mechanisms, John Wiley and Sons, New York.
Howell, L. L., Midha, A., and Murphy, M. D., 1994, “Dimensional Synthesis of
Compliant Constant-Force Slider Mechanisms,” Proceedings of DETC’94,
ASME Design Engineering Technical Conferences, DETC98/MEMD-71.
Howell, L. L., and Midha, A., 1995, “Determination of the Degrees of Freedom of
Compliant Mechanisms Using the Pseudo-Rigid-Body Model Concept,”
Proceedings of the Ninth World Congress on the Theory of Machines and
Mechanisms, Milano, Italy, Vol. 2, p. 1537-1541.
Howell, L. L., and Midha, A., 1996, “A Loop-Closure Theory for the Analysis and
Synthesis of Compliant Mechanisms,” ASME Journal of Mechanical
Design, Vol. 118, No. 1, pp. 121-125.
Jenuwine, J. G., and Midha, A., 1994, “Synthesis of Single-Input and Multiple-
Output Port Mechanisms with Springs for Specified Energy Absorption,”
Journal of Mechanical Design, Trans. ASME, Vol. 116, No. 3, September,
pp. 937-943.
Midha, A., Murphy, M. D., and Howell, L. L., 1995, “Compliant Constant-Force
Mechanism and Devices Formed Therein”, U.S. Patent 5649454, Issued
July 22, 1997.
Millar, A. J., Howell, L. L, and Leonard, J. N., 1996, “Design and Evaluation of
Compliant Constant-Force Mechanisms,” Proceedings of the 1996 ASME
Design Engineering Technical Conferences and Computers in Engineering
Conference, 96-DETC/MECH-1209.
Murphy, M. D., Midha, A., and Howell, L. L., 1994, “Methodology for the Design
of Compliant Mechanisms Employing Type Synthesis Techniques with
Example,” Proceedings of the 1994 ASME Mechanisms Conference, DE-
Vol. 70, pp. 61-66.
Parkinson, M. B., Howell, L. L., and Cox, J. J., 1997, “A Parametric Approach to
the Optimization-Based Design of Compliant Mechanisms,” Proceedings
of the 23rd Design Automation Conference, DETC97/DAC-3763.
184
Paul, B. 1979, Kinematics and Dynamics of Planar Machinery, Prentice-Hall,
Englewood Cliffs, New Jersey.
Wahl, A., 1963, Mechanical Springs, 2nd. Ed. McGrawHill, New York.
Wittwer, J. W., 2001, “Predicting the Effects of Dimensional and Material Property
Variations in Compliant Mechanisms,” M.S. Thesis, Brigham Young Uni-
versity, Provo, Utah.
Williman, J., 1995, “Small Torque,” Engineering (London), Gillard Welch Associ-
ates, pp. 27-28.
185
186
APPENDIX A PSEUDO-RIGID-BODY
MODEL
The Pseudo-Rigid-Body Model (PRBM) plays an important part in the design and
analysis of compliant mechanisms. This appendix offers a closer look at the PRBM. It
gives a brief overview and then discusses the nomenclature and equations for several dif-
ferent types of flexible segments. However, no attempt is made to show evidence of the
validity of the PRBM or its limitations. For further information and details, the reader is
referred to the work from which this appendix was summarized (Howell, 2000).
ally, the large non-linear deflections have caused significant difficulties in the design of
compliant mechanisms. Techniques such as finite element analysis (FEA) and elliptic inte-
grals provide accurate information, but make design very drawn out and complicated.
Fortunately, the development of the PRBM has greatly increased the speed and
ease in which compliant mechanisms can be designed. The PRBM allows for the approxi-
187
mation of the force-deflection characteristics of flexible segments. Thus, the PRBM is
intended to be an intermediate design tool, allowing for the rapid design and analysis of
first generation compliant mechanisms. Afterwards, techniques such as FEA and other
numerical methods can be used to refine the designs. The PRBM becomes a tool to take
The power of the PRBM comes from its ability to model compliant members using
rigid members that have the same force-deflection characteristics as the original member.
Continuous work in developing the PRBM has shown it to accurately model the behavior
ers can draw from the vast number of traditional mechanism design and analysis tools.
For each type of flexible segment, several parameters are defined. The first of
these is the characteristic pivot. The characteristic pivot is the center of the arc created by
the path of the end of the beam. This pivot lies on the flexible beam and is represented as a
pin joint in the PRBM. Equations for the position of this characteristic pivot is given for
each type of flexible segment and is easily determined. Additionally, the variable Θ is the
pseudo-rigid-body angle and equations relating it to the end angle of the beam are pre-
sented.
spring with a spring constant of K. This spring constant is determined by geometric and
material properties and is also dependent upon the type of flexible member. Formulas for
188
P P
nP nP
θo Θ
P P
l
2
nP nP
l L Torsional l L
Spring
Figure A.1 (a) A small-length flexural pivot and (b) its PRBM
These are the main parameters that will be described below. Each section gives a
diagram for the flexible segment and its corresponding PRBM diagram, the characteristic
pivot, and the torsional spring constant. Additionally, the formula for the maximum stress
L»l (A.1)
( EI )L » ( EI )l (A.2)
189
The characteristic pivot is located at the center of the small-length flexural pivot as
shown in Figure A.1b. For small-length flexural pivots, the basic equations are
Θ = θo (A.3)
( EI )
K = ------------l (A.4)
l
T = KΘ (A.5)
The x and y coordinates of the end of the beam can be found through
l l
a = --- + L + --- cos Θ (A.6)
2 2
and
l
b = L + --- sin Θ (A.7)
2
– ( Pa + nPb )c nP
σ top = ---------------------------------- – ------- (A.8)
I A
( Pa + nPb )c nP
σ bot = ------------------------------- – ------- (A.9)
I A
The equations presented here are sufficient in most cases. At times, the size of the
small-length flexural pivot is small enough that the spring constant can be ignored if other
larger torques are present. This special case small-length flexural pivots are called living
hinges.
190
P P
nP nP
P P
nP nP
L Torsional γL
Spring (1-γ)L
Figure A.2 (a) Cantilever beam with force at end and (b) its PRBM
A second type of flexible segment is a cantilever beam with a force at the end. Fig-
ure A.2a show the cantilever beam with its PRBM and corresponding parameters (Figure
A.2b).
The characteristic pivot is located a distance γL from the free end where γ is the
characteristic radius factor. The value of the characteristic radius factor is a function of
the direction of the applied force and can be expressed in terms of n as follows:
191
For values of n between -0.5 and 1.0,
There is some slight deviation between the pseudo rigid body angle and the actual
angle of the beam. This variation is almost linear and is compensated through
θo = cθ Θ (A.14)
where
c θ ≈ 1.24 (A.15)
The value of cθ, the parametric angle coefficient, varies between 1.256 and 1.179. Values
The torsional spring constant for a cantilever beam with a force on the end can be
found from
EI
K = γK θ ------ (A.16)
l
K θ ≈ 2.65 (A.17)
The value Kθ is called the stiffness coefficient. The approximation given in Equation
(A.17) is accurate in most cases. However, more accurate values are given in Howell
(2000). With the spring constant, force and torque calculations can be made according to
Equation (A.5).
The x and y coordinates of the end of the beam can be found through
a = L – γL ( 1 – cos Θ ) (A.18)
192
M
M
θo Θ
L γL
Torsional (1-γ)L
Spring
Figure A.3 (a) Cantilever beam with force at end and (b) its PRBM
and
b = γL sin Θ (A.19)
The stress in the beam can be calculated from Equation (A.8) and Equation (A.9),
which are the same equations used for the small-length flexible segment.
A cantilever beam is often loaded with an end moment. Figure A.3 show this load-
ing configuration along with its PRBM. The equations for this configuration are identical
to the previous configuration. However, there are some differences in the values of the
γ = 0.7346 (A.20)
193
Centerline
P P
L γL
Torsional (1-γ)L (1-γ)L
Spring 2 2
c θ = 1.5164 (A.21)
K Θ = 2.0643 (A.22)
A common type of flexible segment is the fixed-guided beam. This beam consists
of a beam fixed on one end, while the other end is kept perpendicular to ground during dis-
placement. This is commonly see in mechanisms such as parallel and folded beam mecha-
Close observation of the fixed-guided beam shows that the curvature is zero at the
middle due to symmetry. It is also known that the curvature at the end of a cantilever beam
with an end load is also zero. Therefore, the fixed-guided beam can be modeled as two
194
cantilever beams each half the length of the original beam. Thus, the PRBM is easily
found.
The parameters for this beam are similar to previous parameters. The characteristic
radius factor is
γ = 0.8517 (A.23)
Since the beam has a constant end angle, the parametric angle coefficient is trivial
and
cθ = 0 (A.24)
EI
K = 2γK θ ------ (A.25)
l
From observation, it can be seen that the spring constant for the fixed-guided beam
is twice that for cantilever beams. This indicates that the overall stiffness of the fixed-
Since the end of the beam is constrained, a reactionary moment, Mo is created. The
Pa
M o = ------- (A.26)
2
or
Pl
M o = ----- [ 1 – γ ( 1 – cos Θ ) ] (A.27)
2
195
The maximum stress occurs at the ends of the beam where the moment is largest. It
has a value of
Pac
σ max = ---------- (A.28)
2I
196
APPENDIX B MODEL AND
OPTIMIZATION CODE
The CFM model is written in Matlab. Matlab has the ability to manipulate files,
store matrices, and use functions. However, a significant advantage is that no compiling is
required as in C. Changes can be made and the program can be execute immediately. This
allows for quick and easy changes. Unfortunately, this software requires that Matlab is
The 4 separate Matlab files make up the CFM model. These files each serve a dif-
ferent purpose allowing the model to be either independently or linked to the optimization
code. The different operating paths of the model and the run order of the files are summa-
rized in Figure B.1. The primary file, CFMModel, is the model core. This file contains all
197
Independent of
Optimization Linked to Optimization
UserModel
Run Order
OptdesLink
OptdesModel
CFMModel
The UserModel file is used by the user to run the model independent of the optimi-
zation software. The user can change the input values within this file. This file then calls
and passes all important parameters to CFMModel. UserModel then receives all the output
receives all the needed information from data files created by the optimization code. It
then calls CFMModel and passes and receives all important values. Its final duty is to cre-
The final file is OptdesLink. This file handles the handshaking that goes on
between Matlab and OptdesX. It triggers OptdesModel, telling it when the optimization
code is finished creating data files. It then tells the optimization code when Matlab is fin-
198
B.1.2 Variable Name Mapping
Due to the number of variables required, there relationship to one another, and
requirements for input and output of data, the variable names in the Matlab and OptdesX
code are not necessarily intuitive. In some cases, to simplify the writing of the code, two
names were used. Table B.1 outlines the variable names, any second names that may exist,
and the parameter or item that the variable represents. This table can be used to help
%[cfparam,deltax,phi]=cfmmodel(r2,r3,k1,k2,k3,percentage)
%percentage in
%r2,r3
Table B.1 Mapping of parameter names to names in Matlab and OptdesX code
Variable Name Alternate Name Parameter Variable Name Alternate Name Parameter
199
function[lengths,cfparam,configparam,pivotparam,geoparam,deltax,F,di,d,SF,pc,alpha,Beta,fparam,phi,ka,C,Des
ignParam]=cfmmodel(ltot,r4type,r4,di,d,ds,geoparam,pivotparam,cfparam,matprop,fit)
gamma=[0,1,0.85];
roedata=[1,0.1,1/0.85]; %One in first place is to avoid dividing by zero.
Ktheta=[0,1,2.65];
fgoal=.4518;
%10.2483;
Ldata=[1,1.05,gamma(3)];
%Spot at which f parameter is compared.
%if d==16
% famount=15.9/100;
%elseif d==40
% famount=0.05;
%end
R=cfparam(1);
geoparam(2,1)=geoparam(1,1);
geoparam(3,1)=geoparam(1,1);
K1=cfparam(2);
K2=cfparam(3);
K3=cfparam(4);
K(1)=1;
K(2)=K1;
K(3)=K2;
K(4)=K3;
k1=pivotparam(1,4);
k2=K1*k1;
k3=K2*k1;
pivotparam(1,4)=k1;
pivotparam(2,4)=k2;
pivotparam(3,4)=k3;
Lvalues(1)=R;
Lvalues(2)=1;
Lvalues(3)=1;
if pivotparam(1,1)==2
Lvalues(2)=Lvalues(2)*Ldata(2);
elseif pivotparam(1,1)==3
Lvalues(1)=Lvalues(1)*Ldata(3);
Lvalues(3)=Lvalues(3)+1;
end
if pivotparam(2,1)==3
if pivotparam(1,1)==2
Lvalues(2)=Lvalues(2)*Ldata(3);
Lvalues(3)=Lvalues(3)+1;
elseif (pivotparam(2,2)==3)
Lvalues(2)=Lvalues(2)*Ldata(3);
Lvalues(3)=Lvalues(3)+1;
else
Lvalues(1)=Lvalues(1)*Ldata(3);
Lvalues(3)=Lvalues(3)+1;
end
end
if pivotparam(3,1)==2
Lvalues(1)=Lvalues(1)*Ldata(2);
elseif pivotparam(3,1)==3
Lvalues(2)=Lvalues(2)*Ldata(3);
Lvalues(3)=Lvalues(3)+1;
end
if Lvalues(3)==1
200
Lvalues(3)=R+1;
else
Lvalues(3)=Ldata(3)*(R+1);
end
L=(Lvalues(1)+Lvalues(2))/Lvalues(3);
rtot=ltot/L;
r3=rtot/((1+1/R)); %Calculate r3 value
r2=rtot/(1+R); %Calculate r2 value
d=d/100;
class1B=0; %Boolean for class 1B config.
for i=1:3
roe(i)=roedata(pivotparam(i,1));
zeta(1)=R+1;
zeta(3)=1+1/R;
if (pivotparam(2,1)==2)%Small length in center
zeta(2)=2;
elseif (pivotparam(2,1)==3) %Long flexible in center
if pivotparam(1,1)==1
if (pivotparam(2,2)==2) %Associated with r2
zeta(2)=zeta(1);
else
zeta(2)=zeta(3); %Associated with r3
end
else
zeta(2)=zeta(3);
end
else
zeta(2)=1; %Sets to arbitary 1 (not zero so can divide w/o error
end
geoparam(i,3)=rtot*roe(i)/zeta(i); %Calculate flexible segment lengths
if pivotparam(i,1)==1
I(i)=0;
geoparam(i,2)=0;
else
I(i)=pivotparam(i,4)*geoparam(i,3)/(gamma(pivotparam(i,1))*Ktheta(pivotparam(i,1))*matprop(1)); %Calcu-
late spring constants
geoparam(i,2)=(12*I(i)/geoparam(i,1))^(1/3);
end
A(i)=geoparam(i,2)/(2*rtot); %Calculate A value
ka(i)=gamma(pivotparam(i,1))*zeta(i)*Ktheta(pivotparam(i,1))/roe(i);
end
if (pivotparam(1,1)==1)
class1B=1;
geoparam(1,3)=geoparam(1,2);
pivotparam(1,4)=pivotparam(2,4);
A(1)=A(2);
ka(1)=ka(2);
end
type=0;
if (r4type==1) %Forces Orthoplaner position by changing r4 from user defined to ideal for orthoplaner
r4=(1-R)*r2;
end
ks=K3*k1/(r2*r2);
201
pivotparam(4,4)=ks;
rimin=0;
if (r3+r4>r2) %& ~(r4==0) %Limits ri so that mechanisms doesn't go back on itself
theta3max=asin((r4-r2)/r3); %Limits to where r2 is vertical. Finds theta3max
rimin=r3*cos(theta3max); %Sets start point and limits for non-linear solver
thetaT = pi/4;
limit = pi/2;
elseif (r3+r4<r2) %& ~(r4==0) %Limits to where r3 and r4 are vertical. Finds theta2max
theta2max=asin((r4+r3)/r2);
rimin=r2*cos(theta2max);
thetaT = theta2max/2;
limit = theta2max;
elseif (r3+r4==r2) %OrthoPlaner position
thetaT=pi/4;
limit=2*pi;
end
202
for i=1:sz
darray(i)=deltax(i)/ri*100;
end
%bo=0;
%for i=1:sz
% if (deltax(i)>ri*famount) & (bo==0)
% bo=1;
% fspot=i
% end
%end
fspot=sz; %calculate f at d
deltax(fspot)=ri*d;
if (d>0)
r1(fspot)=ri-deltax(fspot);
else
r1(fspot)=ri+deltax(fspot);
end
changepivot=-1;
maxstresspivot=0;
dd=d*100;
dN=100*(1/(R+1))*sin(acos((-1/200)*((dd^2*R-200*dd*R-200*dd+20000+dd^2)/(dd-100))));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Basic slider crank model
for i=1:sz
if r4==0
theta2(i)=acos((r1(i)^2+r2^2-r3^2)/(2*r1(i)*r2));
tempangle=asin((-r2*sin(theta2(i)))/r3);
else %solves equations using Newton Raphson if r4!=0
del = 1e-6;
tol = 1e-3;
limit=2*pi;
if i==1
theta2temp = Newton(type,'theta2',theta2i,[r1(i),r2,r3,r4],limit,del,tol);
else
theta2temp = Newton(type,'theta2',theta2(i-1),[r1(i),r2,r3,r4],limit,del,tol);
end
theta2(i) = theta2temp;
tempangle = asin((r4-r2*sin(theta2(i)))/r3)+2*pi;
end
if tempangle<0
theta3(i)=2*pi+tempangle;
else
theta3(i)=tempangle;
end
if i==sz
theta2final=theta2(i);
theta3final=theta3(i);
end
203
phi(i)=(R*cos(theta3(i))*((theta2(i)-theta2o)+K1*(theta3o+theta2(i)-theta3(i)-theta2o))...
+cos(theta2(i))*(K1*(theta3o+theta2(i)-theta3(i)-theta2o)+K2*(theta3o-theta3(i))))/(R*sin(theta2(i)-theta3(i)));
else
phi(i)=(R*cos(theta3(i))*(K1*(theta3o+theta2(i)-theta3(i)-theta2o))...
+cos(theta2(i))*(K1*(theta3o+theta2(i)-theta3(i)-theta2o)))/(R*sin(theta2(i)-theta3(i)));
end
phis(i)=K3*((cos(theta2(i))-(1-ds)*cos(theta2i))+R*(cos(theta3(i))-(1-ds)*cos(theta3i)));
F(i)=k1*(phi(i)-phis(i))/r2; %Force equation
%Calculate Beta
if class1B==1
Beta(i)=gamma(pivotparam(2,1))*Ktheta(pivotparam(2,1))*(R+1)*zeta(2)*phi(i)/roe(2);
else
Beta(i)=gamma(pivotparam(1,1))*Ktheta(pivotparam(1,1))*(R^2+2*R+1)*phi(i)/roe(1);
end
%Calculate f parameter by normalizing highest stress parameter (Offset for K already applied)
fparam(i)=ka(maxstresspivot)*(1/alpha(maxstresspivot,i))^3/(ka(1)*L*K(maxstresspivot));
204
%Equations for calculating force on sliding Cam
Fx(i)=((theta2(i)-theta2o)*k1)/((cos(theta3o-theta3(i))*sin(theta2(i)-theta2o)*r2)/sin(theta3o-theta3(i))...
+cos(theta2(i)-theta2o)*r2);
Fy(i)=(-Fx(i)/tan(theta3o-theta3(i)));
Fcam(i)=sqrt(Fx(i)^2+Fy(i)^2);
end
for i=1:3
if alpha(i,sz)==0
C(i)=-1;
else
C(i)=alpha(maxstresspivot,sz)/alpha(i,sz);
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:3
if(pivotparam(i,1)==1)
DesignParam(i,1)=-1
DesignParam(i,2)=-1
DesignParam(i,3)=-1
DesignParam(i,4)=-1
else
DesignParam(i,1)=(K(i)*kapivot/(ka(i)*Kpivot))^(1/3) %maximum h given same stress and b values
DesignParam(i,2)=(K(i)*kapivot/(ka(i)*Kpivot)) %maximum b given same stress and h values
DesignParam(i,3)=(K(i)*kapivot/(ka(i)*Kpivot)*C(i)^3) %D ratio of maximum b to primary stress pivot b
given max h on x
DesignParam(i,4)=(11-DesignParam(i,3))*(C(i)) %A
end
end
[alphamax,temp]=max(alpha(maxstresspivot,:));
alphamax;
Amax=A(pivot);
maxstress=max(stress(maxstresspivot,:));
SF=matprop(3)/alphamax/Amax;
lengths(1)=ri;
lengths(2)=r2;
lengths(3)=r3;
lengths(4)=r4;
lengths(5)=L;
pivotparam(1,5)=theta2i;
pivotparam(2,5)=theta3i;
pivotparam(1,6)=theta2final;
pivotparam(2,6)=theta3final;
205
average_phi=average(phi);
average_phis=average(phis);
average_Beta=average(Beta);
average_alpha=average(alpha(maxstresspivot,:));
Fmin=min(F);
Fmax=max(F);
if d>0
xi=((Fmin/Fmax))*100;
else
xi=abs(Fmin/Fmax);
end
fp=fparam(fspot);
fp=fparam(fspot)/fgoal;
pc=[-1,-1];
if fit==1
X=[1,1];
Options = optimset('TolFun',.0001);
LL=[];
UL=[];
pc=lsqcurvefit('powerb',X,darray,alpha(maxstresspivot,:),LL,UL,Options); %Curve fit alpha
end
ForceAverage=average(F);
cfparam=[R,K1,K2,K3,average_phi,average_phis,xi];
configparam=[alphamax,pivot,changepivot,average_Beta,fp,dN];
dd=d*100;
d=dd;
%This is the user interface for the CFORIGINAL model. The model can also be accessed through
%the optdes interface
%function usermodel
clear all;
close all;
%User Inputs
matprop(1)=1; %E%E
matprop(2)=1;%Sy %Sy
matprop(3)=matprop(2)/matprop(1);%Sy/E
K3=0;
R=0.395;
K1=.1906;
K2=0;
cfparam(2)=K1;
cfparam(3)=K2;
206
ltot=1;
k1=.60106;
b=1;
geoparam(1,1)=b;
geoparam(2,1)=b;
geoparam(3,1)=b;
pivotparam(1,4)=k1;
r4=0;
r4type=0; %0==> User input 1 ==> Orthoplaner
di=0; %0..1
d=16; %percent displacement (0..100) +--> percent of ri to displace - -->percent of (r2+r3-ri) to displace
ds=0; %Percent of rest length to preload spring
%Pin joint types 1=pin 2=small 3=long
pivotparam(1,1)=2;
pivotparam(2,1)=3;
pivotparam(3,1)=1;
%Pin Joint reference links for long length segments
%Should only matter for pivot 2 Can be either link 2 or link 3
pivotparam(2,2)=3;
%Initial angles can be inputed through the pivot parameters. The defaults for
%the normal slider crank are as follows:
%Initial wind up on springs --> pos. # is a counter clockwise rotation results in negative torque
%With no pre-load, the springs will be at rest in rest posistion (theta#i)
pivotparam(1,3)=0;
pivotparam(2,3)=0;
pivotparam(3,3)=0;
fit=0;
cfparam(1)=R;
cfparam(4)=K3;
[lengths,cfparam,configparam,pivotparam,geoparam,deltax,force,di,d,SF,pc,alpha,Beta,fp,phi,ka,C,DesignPara
m]=cfmmodel(ltot,r4type,r4,di,d,ds,geoparam,pivotparam,cfparam,matprop,fit);
%Plot Force output
figure(1)
plot(deltax,force);
axis([0,max(deltax),0,max(force)*1.1]);
%Plot mechanism in initial and fully deflected position
figure(2)
hold on;
r2=lengths(2);
axis equal;
theta2i=pivotparam(1,5);
theta3i=pivotparam(2,5);
r2ix=cos(theta2i)*r2;
r2iy=sin(theta2i)*r2;
x=[0,r2ix];
y=[0,r2iy];
r3=lengths(3);
r3ix=r2ix+cos(pivotparam(2,5))*r3;
r3iy=r2iy+sin(pivotparam(2,5))*r3;
plot(x,y,'g');
x=[r2ix,r3ix];
y=[r2iy,r3iy];
plot(x,y,'b');
x=[r3ix,r3ix];
y=[r3iy,0];
plot(x,y,'r');
if r2iy > r3ix
range = r2iy+1;
else
range = r3ix+1;
end
theta2final=pivotparam(1,6);
theta3final=pivotparam(1,6);
207
r2ix=cos(theta2final)*r2;
r2iy=sin(theta2final)*r2;
x=[0,r2ix];
y=[0,r2iy];
r3ix=r2ix+cos(pivotparam(2,6))*r3;
r3iy=r2iy+sin(pivotparam(2,6))*r3;
plot(x,y,'g-.');
x=[r2ix,r3ix];
y=[r2iy,r3iy];
plot(x,y,'b-.');
x=[r3ix,r3ix];
y=[r3iy,0];
plot(x,y,'r-.');
hold off;
%if r2iy > r3ix & r2iy>range
% range = r2iy+1;
%elseif r3ix>range
% range = r3ix+1;
%end
%axis([0 range 0 range]);%
%plotmech(2,r2,r3,theta2i,theta3i,theta2final,theta3final);
%plot(deltax,cfparam(5))
%postprocess2(r2,r3,percentage,pivotparam,matprop,cfparam,deltax,phi,force)
function optdesmodel(x)
load cfc_Data1.txt;
matprop(1)=cfc_Data1(1,1);
matprop(2)=cfc_Data1(1,2);
matprop(3)=matprop(2)/matprop(1);
ltot=cfc_Data1(1,3);
pivotparam(1,4)=cfc_Data1(1,4);
geoparam(1,1)=cfc_Data1(1,5);
cfparam(1)=cfc_Data1(2,1);
cfparam(2)=cfc_Data1(2,2);
cfparam(3)=cfc_Data1(2,3);
cfparam(4)=cfc_Data1(2,4);
di=cfc_Data1(2,5);
d=cfc_Data1(2,6);
ds=cfc_Data1(2,7); %Percent of rest length to preload spring
208
pivotparam(1,1)=cfc_Data1(3,1);
pivotparam(2,1)=cfc_Data1(3,2);
pivotparam(3,1)=cfc_Data1(3,3);
r4type=cfc_Data1(3,4);
r4=cfc_Data1(3,5);
%Initial angles can be inputed through the pivot parameters. The defaults for
%the normal slider crank are as follows:
%Initial wind up on springs --> pos. # is a counter clockwise rotation results in negative torque
%With no pre-load, the springs will be at rest in rest posistion (theta#i)
pivotparam(1,3)=cfc_Data1(5,1);
pivotparam(2,3)=cfc_Data1(5,2);
pivotparam(3,3)=cfc_Data1(5,3);
fit=cfc_Data1(6,1);
[lengths,cfparam,configparam,pivotparam,geoparam,deltax,force,di,d,SF,pc,alpha,Beta,fp,phi,ka,C,DesignPara
m]=cfmmodel(ltot,r4type,r4,di,d,ds,geoparam,pivotparam,cfparam,matprop,fit);
averageforce=average(force)
209
B.2 Optimization Code
zation routines. It is necessary to construct a C file called anasubC. This file establishes
the analysis variables, analysis functions, and the model or links to the model. The vari-
able names in the anasubC code are mapped to the given parameter names in Table B.1
#include "supportC.h"
#include "math.h"
#include "stdlib.h"
#include "string.h"
#include <stdio.h>
double ri,r2,r3,r4,r4type,d,L,di,ds,E,Sy,SF,rtot;
double pivotparam[6][6];
double hp1,bp1,D1,Area1,hp2,bp2,D2,Area2,hp3,bp3,D3,Area3;
double AA1,BB1,CC1,DD1,AA2,BB2,CC2,DD2,AA3,BB3,CC3,DD3;
double geoparam[3][3];
double cfparam[7];
double lengths[4];
double configparam[5];
double C1,C2,C3;
double R,K1,K2,xi,averageforce,m,n,ka1,ka2,ka3,dN;
double k1,k2,k3,ks,averagephi,averagephis;
double alphamax, pivot, changepivot,averagebeta,fp,ltot;
double fit,b1,h1,b2,h2,b3,h3,l1,l2,l3;
int hold;
/*=============================================================================
Function anapreC
Preprocessing Function
-----------------------------------------------------------------------------*/
#ifdef __STDC__
void anapreC( char *modelName )
#else
void anapreC( modelName )
char *modelName;
#endif
{
/* set model name (16 chars max) */
strcpy( modelName, "Original CF" );
}
/*=============================================================================
Function anafunC
Analysis Function
-----------------------------------------------------------------------------*/
#ifdef __STDC__
void anafunC( void )
210
#else
void anafunC( )
#endif
{
avdscaC(&E,"E");
avdscaC(&Sy,"Sy");
avdscaC(<ot,"ltot");
avdscaC(&pivotparam[0][3],"k1");
avdscaC(&geoparam[0][0],"b");
avdscaC(&cfparam[0],"R");
avdscaC(&cfparam[1],"K1");
avdscaC(&cfparam[2],"K2");
avdscaC(&cfparam[3],"K3");
avdscaC(&di,"% Defl. Intial");
avdscaC(&d,"d");
avdscaC(&ds,"% Spring Pre-load");
avdscaC(&r4,"r4");
avdscaC(&r4type,"r4 type");
avdscaC(&pivotparam[0][0],"Pin 1 Type");
avdscaC(&pivotparam[1][0],"Pin 2 Type");
avdscaC(&pivotparam[2][0],"Pin 3 Type");
avdscaC(&pivotparam[1][1],"Pin 2 Reference");
avdscaC(&pivotparam[0][2],"Theta2 Preload");
avdscaC(&pivotparam[1][2],"Theta3 Preload");
avdscaC(&pivotparam[2][2],"Theta4 Preload");
avdscaC(&fit,"Fit 1Yes 0No");
211
/*Checks to see if Hold existes */
flag=fopen("Hold.txt","r");
if(flag!=NULL)
{fclose(flag);
system("rm Hold.txt");
}
fclose(flag);
out = fopen("cfc_Data1.txt","w");
flag=fopen("Stop.txt","w");
fprintf(flag,"%g \n",1.00);
fclose(flag);
flag=NULL;
while (flag==NULL)
{
flag= fopen("Hold.txt","r");
}
fclose(flag);
system("rm Hold.txt");
fscanf(in,"%lf%lf%lf%lf%lf",&ri,&r2,&r3,&r4,&L);
fscanf(in,"%lf%lf%lf",&b1,&h1,&l1);
fscanf(in,"%lf%lf%lf",&b2,&h2,&l2);
fscanf(in,"%lf%lf%lf",&b3,&h3,&l3);
fscanf(in,"%lf%lf%lf%lf%lf%lf",&pivotparam[0][0],&pivotparam[0][1],&pivotparam[0][2],&pivotparam[0][3],&p
ivotparam[0][4],&pivotparam[0][5]);
fscanf(in,"%lf%lf%lf%lf%lf%lf",&pivotparam[1][0],&pivotparam[1][1],&pivotparam[1][2],&pivotparam[1][3],&p
ivotparam[1][4],&pivotparam[1][5]);
fscanf(in,"%lf%lf%lf%lf%lf%lf",&pivotparam[2][0],&pivotparam[2][1],&pivotparam[2][2],&pivotparam[2][3],&p
ivotparam[2][4],&pivotparam[2][5]);
fscanf(in,"%lf%lf%lf%lf%lf%lf",&pivotparam[3][0],&pivotparam[3][1],&pivotparam[3][2],&pivotparam[3][3],&p
ivotparam[3][4],&pivotparam[3][5]);
fscanf(in,"%lf%lf%lf%lf%lf%lf%lf",&cfparam[0],&cfparam[1],&cfparam[2],&cfparam[3],&cfparam[4],&cfpara
m[5],&cfparam[6]);
212
fscanf(in,"%lf%lf%lf%lf%lf%lf",&alphamax,&pivot,&changepivot,&averagebeta,&fp,&dN);
fscanf(in,"%lf",&averageforce);
fscanf(in,"%lf",&di);
fscanf(in,"%lf",&d);
fscanf(in,"%lf",&SF);
fscanf(in,"%lf%lf",&m,&n);
fscanf(in,"%lf%lf%lf",&ka1,&ka2,&ka3);
fscanf(in,"%lf%lf%lf",&C1,&C2,&C3);
fscanf(in,"%lf%lf%lf%lf",&hp1,&bp1,&D1,&Area1);
fscanf(in,"%lf%lf%lf%lf",&hp2,&bp2,&D2,&Area2);
fscanf(in,"%lf%lf%lf%lf",&hp3,&bp3,&D3,&Area3);
fclose(in);
k1=pivotparam[0][3];
k2=pivotparam[1][3];
k3=pivotparam[2][3];
ks=pivotparam[3][3];
afdscaC(pivotparam[0][4],"Theta2 initial");
afdscaC(pivotparam[1][4],"Theta3 initial");
afdscaC(pivotparam[0][5],"Theta2 final");
afdscaC(pivotparam[1][5],"Theta3 final");
afdscaC(ri,"r initial");
afdscaC(r2,"r2");
afdscaC(r3,"r3");
afdscaC(r4,"r4");
213
afdscaC(l2,"l2");
afdscaC(b3,"b3");
afdscaC(h3,"h3");
afdscaC(l3,"l3");
afdscaC(ka1,"ka1");
afdscaC(ka2,"ka2");
afdscaC(ka3,"ka3");
afdscaC(C1,"C1");
afdscaC(C2,"C2");
afdscaC(C3,"C3");
afdscaC(AA1,"AA 1");
afdscaC(AA2,"AA 2");
afdscaC(AA3,"AA 3");
afdscaC(BB1,"BB 1");
afdscaC(BB2,"BB 2");
afdscaC(BB3,"BB 3");
afdscaC(CC1,"CC 1");
afdscaC(CC2,"CC 2");
afdscaC(CC3,"CC 3");
afdscaC(DD1,"DD 1");
afdscaC(DD2,"DD 2");
afdscaC(DD3,"DD 3");
/*=============================================================================
Function anaposC
Postprocessing Function
-----------------------------------------------------------------------------*/
#ifdef __STDC__
void anaposC( void )
#else
void anaposC( )
#endif
{
out = fopen("cfc_postdata1.txt","w");
/*variable=getValue(in,#));*/
214
fprintf(out,"%lf %lf %lf %lf 0 0 0\n",geoparam[1][0],geoparam[1][1],geoparam[1][2],geoparam[1][3]);
fprintf(out,"%lf 0 0 0 0 0 0\n",averageforce);
fprintf(out,"%lf 0 0 0 0 0 0\n",di);
fprintf(out,"%lf 0 0 0 0 0 0\n",d);
fclose(out);
system("matlab <postprocess.m");
}
The Matlab model code and the OptdesX optimization package are linked together
through "hand shaking" that occurs through data files. One piece of software signals the
other piece of software that it is finished with its part of the problem by creating or editing
a data file. Meanwhile, the software that is not currently doing anything in the problem
waits for the other software to signal that it is finished. Ideally, the two programs would
share a common database and have internal triggers. However, with the software being
while 1
stop=0;
load Stop.txt
while Stop == 0.0
load Stop.txt
end
temp=1.0;
215
temp2=0.000;
save Stop.txt temp2 -ASCII;
Stop is a file that tells MatLab to wait. MatLab continues to read the file until the
value in the file is no longer zero. This value is changed by anasubC upon once all data has
The MatLab function OptdesModel is then called. Once all the data has been saved
from MatLab to a data file, MatLab creates the file Hold. Meanwhile, anasubC is waiting
To run the linked models, it is necessary to start Matlab and run the OptdesLink
function, as well as start OptdesX. The model runs fairly quickly. It is much faster than
launching the MatLab application each time through a system command. Additionally,
this set up allows for easy debugging and model changes. It is not necessary to terminate
the optimization code. The matlab code can be stopped, the model altered, and OptdesLink
restarted.
216
APPENDIX C OPTIMIZATION PLOTS
69.8
75.6
85.4
85.4
80.2
75.6
69.8
68.5
67.7
67.2
.155
.205
.336
.336
.205
.155
.126
.096
.079
12.0
12.0
.065
66.7
.049
8.00
8.00
K1
K1
65.8 .034
4.00
4.00
.018
62.7
1.00 2.00 3.00 4.00 5.00 1.00 2.00 3.00 4.00 5.00
R R
(a) (b)
Figure C.1 2-D Explore plots for slp-b with R vs. K1 and contours of (a) Ξ'ex and (b) Ψ
217
69.7
75.4
85.4
85.4
80.1
75.4
69.7
68.4
67.6
67.1
.022
.015
.016
.016
.016
.017
.019
.022
.024
.027
.031
.036
12.0
12.0
66.6
8.00
8.00
K1
K1
65.6 .036
4.00
4.00
.031
.027
.022
.019
62.4 .015
1.00 2.00 3.00 4.00 5.00 1.00 2.00 3.00 4.00 5.00
R R
(a) (b)
Figure C.2 2-D Explore plots for ssp-b with R vs. K1 and contours of (a) Ξ'ex and (b) Ψ
20.0
.987
97.4
16.0
.792
95.7
1 1
12.0
.597
93.9
2-Psi
1
3-K2
1-Xi
8.01
.403
92.1
2
4.01
.208
90.4
3 3
2
.0172
.0133
2
88.6
Figure C.3 Optimum plot for lps-b which shows Ξ'ex, Ψ and optimal K2 for given R value
218
54.0
67.7
79.3
93.6
84.2
75.4
71.3
71.3
79.3
84.2
93.6
93.6
88.5
84.2
79.3
75.4
.00876
.00876
.00770
.00682
.00612
.00557
.0149
.0249
.0406
.0406
.0249
.0149
.0107
10.0
10.0
.00512
71.3
8.00
8.00
67.7
6.00
6.00
K2
K2
62.8
4.00
4.00
58.6
2.00
2.00
54.0
.500 1.50 2.50 3.50 4.50 .500 1.50 2.50 3.50 4.50
R R
(a) (b)
Figure C.4 2-D explore plots for sps-b with R vs. K2 and contours of (a) Ξ'ex and (b) Ψ
.0808
99.2
.0660
96.9
1
.0513
94.7
2
2-Psi
1-Xi
.0365
92.4
2
.0218
90.1
1 1
2
.00708
87.8
Figure C.5 Optimum plot of design variables K1 and K2 versus R for sss-b with curves of Ξ'ex and Ψ
219
220
APPENDIX D MECHANISM TABLES
The tables on the following pages summarize all of the new mechanisms along
with their parameters. The lpp-a and lpp-b mechanisms are bolded in each graph so that
they can be quickly found during design. Tables D.1 and D.2 list the mechanisms in the
order in which they were presented. Tables D.3 and D.4 list all of the mechanisms sorted
by M while Tables D.5 and D.6 tabulate all of the mechanisms sorted by Ψ. In the sorted
tables, the mechanisms are separated into 16% (sub-class a) and 40% (sub-class b) deflec-
tion mechanisms.
221
Table D.1 Combined mechanism table
222
Table D.2 Combined mechanism table
Configuration Sub-Class d Nmax C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class Configuration
lpp a 26.96 a lpp
a99 26.90 a99
a95 26.57 a95
a90 26.04 a90
lpp b 39.79 b lpp
b95 39.68 b95
b90 39.47 b90
spp a 26.96 a spp
a99 26.90 a99
a95 26.57 a95
a90 26.04 a90
spp b 39.79 b spp
b95 39.68 b95
b90 39.47 b90
psp a 27.13 a psp
psp b 40.00 b psp
pl p a 27.13 a pl p
a90 24.03 a90
a85 21.34 a85
a80 19.11 a80
pl p b 40.00 b pl p
b85 39.24 b85
b80 36.46 b80
slp aIo 23.23 1.71 1.00 - 2.54 1.00 - 0.51 1.00 - aIo slp
a99IO 24.97 1.13 1.00 - 1.43 1.00 - 1.00 1.00 - a99IO
a95Io 26.95 1.00 1.92 - 1.00 7.04 - 1.00 1.00 - a95Io
a90I 26.28 1.00 4.82 - 1.00 119.69 - 1.00 1.07 - a90I
a90IO 25.39 1.00 1.00 - 1.00 1.00 - 1.00 1.00 - a90IO
slp bIO 30.18 1.47 1.00 - 1.98 1.00 - 0.63 1.00 - bIO slp
b95IO 36.92 1.00 1.09 - 1.00 1.28 - 1.00 1.00 - b95IO
b90Io 38.58 1.00 1.47 - 1.00 3.14 - 1.00 1.00 - b90Io
ssp aIo 23.23 5.07 1.00 - 7.52 1.00 - 0.06 1.00 - aIo ssp
a99I 26.85 2.54 1.00 - 11.20 1.00 - 0.68 1.00 - a99I
a95I 27.09 2.18 1.00 - 10.43 1.00 - 1.01 1.00 - a95I
a95IO 27.12 2.09 1.00 - 2.00 1.00 - 0.22 1.00 - a95IO
a90IO 26.43 1.33 1.00 - 2.00 1.00 - 0.84 1.00 - a90IO
ssp bIo 30.18 4.63 1.00 - 6.24 1.00 - 0.06 1.00 - bIo ssp
b94I 23.08 6.67 1.00 - 15.38 1.00 - 0.05 1.00 - b94I
b90IO 39.74 2.30 1.00 - 2.00 1.00 - 0.17 1.00 - b90IO
b90o 39.93 2.15 1.00 - 10.35 1.00 - 1.05 1.00 - b90o
lps aIo 26.77 9.07 - 1.00 6.86 - 1.00 0.01 - 1.00 aIo lps
a99I 23.10 0.68 - 1.00 0.32 - 1.00 1.00 - 1.00 a99I
a95I 23.63 0.79 - 1.00 0.50 - 1.00 1.00 - 1.00 a95I
a90 24.60 1.00 - 0.84 1.00 - 0.82 1.00 - 1.37 a90
a99IO 24.39 1.00 - 1.00 0.47 - 1.00 0.47 - 1.00 a99IO
a95IO 24.39 1.00 - 1.00 0.66 - 1.00 0.66 - 1.00 a95IO
a90IO 24.39 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00 a90IO
lps bIo 39.60 7.33 - 1.00 6.05 - 1.00 0.02 - 1.00 bIo lps
b99IO 33.91 0.99 - 1.00 0.55 - 1.00 0.56 - 1.00 b99IO
b95IO 33.91 0.99 - 1.00 0.60 - 1.00 0.61 - 1.00 b95IO
b90IO 33.95 1.00 - 1.00 0.68 - 1.00 0.68 - 1.00 b90IO
b95I 39.66 1.00 - 0.13 1.00 - 0.00 1.00 - 0.73 b95I
b90I 39.44 1.00 - 0.12 1.00 - 0.00 1.00 - 0.98 b90I
sps aIO 26.77 1.74 - 1.00 1.31 - 1.00 0.25 - 1.00 aIO sps
a99I 26.67 1.00 - 0.52 1.00 - 0.14 1.00 - 0.98 a99I
a90IO 27.13 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00 a90IO
sps bIo 39.60 1.40 - 1.00 1.16 - 1.00 0.42 - 1.00 bIo sps
b98Io 37.47 2.35 - 1.00 6.01 - 1.00 0.46 - 1.00 b98Io
b86IO 40.00 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00 b86IO
sss a 22.82 1.000 1.332 7.099 1.000 1.832 33.744 1.000 0.775 0.094 a sss
a99 26.37 1.000 0.745 2.222 1.000 0.435 12.492 1.000 1.052 1.139 a99
a95IO 26.71 1.462 1.000 2.693 0.500 1.000 2.000 0.160 1.000 0.102 a95IO
a95I 25.93 0.340 0.289 1.000 0.039 0.022 1.000 0.997 0.895 1.000 a95I
a90IO 26.21 1.260 1.000 3.104 1.974 1.000 1.979 0.986 1.000 0.066 a90IO
sss b 32.44 1.000 1.041 4.338 1.000 1.541 19.533 1.000 1.365 0.239 b sss
b99 34.51 1.000 0.936 3.502 1.000 0.995 13.377 1.000 1.215 0.311 b99
b90I 39.35 0.589 0.381 1.000 0.148 0.055 1.000 0.724 0.996 1.000 b90I
b88IO 39.54 1.610 1.000 2.398 0.500 1.000 2.000 0.120 1.000 0.145 b88IO
Configuration Sub-Class d Nmax C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class Configuration
223
Table D.3 Combined mechanism table sorted by M
224
Table D.4 Combined mechanism table sorted by M
Configuration Sub-Class d Nmax C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class Configuration
lpp a90 26.04 a90 lpp
lpp a95 26.57 a95 lpp
lpp a99 26.90 a99 lpp
lpp a 26.96 a lpp
pl p a85 21.34 a85 pl p
pl p a80 19.11 a80 pl p
pl p a90 24.03 a90 pl p
pl p a 27.13 a pl p
lps a90 24.60 1.00 - 0.84 1.00 - 0.82 1.00 - 1.37 a90 lps
lps a99I 23.10 0.68 - 1.00 0.32 - 1.00 1.00 - 1.00 a99I lps
lps a95I 23.63 0.79 - 1.00 0.50 - 1.00 1.00 - 1.00 a95I lps
lps a90IO 24.39 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00 a90IO lps
lps a95IO 24.39 1.00 - 1.00 0.66 - 1.00 0.66 - 1.00 a95IO lps
lps a99IO 24.39 1.00 - 1.00 0.47 - 1.00 0.47 - 1.00 a99IO lps
slp a90IO 25.39 1.00 1.00 - 1.00 1.00 - 1.00 1.00 - a90IO slp
slp a99IO 24.97 1.13 1.00 - 1.43 1.00 - 1.00 1.00 - a99IO slp
sss a95I 25.93 0.340 0.289 1.000 0.039 0.022 1.000 0.997 0.895 1.000 a95I sss
spp a90 26.04 a90 spp
spp a95 26.57 a95 spp
slp aIo 23.23 1.71 1.00 - 2.54 1.00 - 0.51 1.00 - aIo slp
sps a99I 26.67 1.00 - 0.52 1.00 - 0.14 1.00 - 0.98 a99I sps
spp a99 26.90 a99 spp
slp a95Io 26.95 1.00 1.92 - 1.00 7.04 - 1.00 1.00 - a95Io slp
spp a 26.96 a spp
sps a90IO 27.13 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00 a90IO sps
sps aIO 26.77 1.74 - 1.00 1.31 - 1.00 0.25 - 1.00 aIO sps
lps aIo 26.77 9.07 - 1.00 6.86 - 1.00 0.01 - 1.00 aIo lps
sss a99 26.37 1.000 0.745 2.222 1.000 0.435 12.492 1.000 1.052 1.139 a99 sss
slp a90I 26.28 1.00 4.82 - 1.00 119.69 - 1.00 1.07 - a90I slp
psp a 27.13 a psp
ssp a95IO 27.12 2.09 1.00 - 2.00 1.00 - 0.22 1.00 - a95IO ssp
ssp a95I 27.09 2.18 1.00 - 10.43 1.00 - 1.01 1.00 - a95I ssp
ssp a99I 26.85 2.54 1.00 - 11.20 1.00 - 0.68 1.00 - a99I ssp
sss a95IO 26.71 1.462 1.000 2.693 0.500 1.000 2.000 0.160 1.000 0.102 a95IO sss
ssp a90IO 26.43 1.33 1.00 - 2.00 1.00 - 0.84 1.00 - a90IO ssp
sss a90IO 26.21 1.260 1.000 3.104 1.974 1.000 1.979 0.986 1.000 0.066 a90IO sss
ssp aIo 23.23 5.07 1.00 - 7.52 1.00 - 0.06 1.00 - aIo ssp
sss a 22.82 1.000 1.332 7.099 1.000 1.832 33.744 1.000 0.775 0.094 a sss
Configuration Sub-Class d Nmax C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class Configuration
lps b90I 39.44 1.00 - 0.12 1.00 - 0.00 1.00 - 0.98 b90I lps
lpp b90 39.47 b90 lpp
lps b95I 39.66 1.00 - 0.13 1.00 - 0.00 1.00 - 0.73 b95I lps
lpp b95 39.68 b95 lpp
lpp b 39.79 b lpp
pl p b80 36.46 b80 pl p
pl p b85 39.24 b85 pl p
pl p b 40.00 b pl p
lps b95IO 33.91 0.99 - 1.00 0.60 - 1.00 0.61 - 1.00 b95IO lps
lps b99IO 33.91 0.99 - 1.00 0.55 - 1.00 0.56 - 1.00 b99IO lps
lps b90IO 33.95 1.00 - 1.00 0.68 - 1.00 0.68 - 1.00 b90IO lps
slp bIO 30.18 1.47 1.00 - 1.98 1.00 - 0.63 1.00 - bIO slp
slp b95IO 36.92 1.00 1.09 - 1.00 1.28 - 1.00 1.00 - b95IO slp
slp b90Io 38.58 1.00 1.47 - 1.00 3.14 - 1.00 1.00 - b90Io slp
sss b90I 39.35 0.589 0.381 1.000 0.148 0.055 1.000 0.724 0.996 1.000 b90I sss
spp b90 39.47 b90 spp
spp b95 39.68 b95 spp
spp b 39.79 b spp
sps b86IO 40.00 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00 b86IO sps
lps bIo 39.60 7.33 - 1.00 6.05 - 1.00 0.02 - 1.00 bIo lps
sps bIo 39.60 1.40 - 1.00 1.16 - 1.00 0.42 - 1.00 bIo sps
sps b98Io 37.47 2.35 - 1.00 6.01 - 1.00 0.46 - 1.00 b98Io sps
sss b99 34.51 1.000 0.936 3.502 1.000 0.995 13.377 1.000 1.215 0.311 b99 sss
sss b 32.44 1.000 1.041 4.338 1.000 1.541 19.533 1.000 1.365 0.239 b sss
psp b 40.00 b psp
ssp b90o 39.93 2.15 1.00 - 10.35 1.00 - 1.05 1.00 - b90o ssp
ssp b90IO 39.74 2.30 1.00 - 2.00 1.00 - 0.17 1.00 - b90IO ssp
sss b88IO 39.54 1.610 1.000 2.398 0.500 1.000 2.000 0.120 1.000 0.145 b88IO sss
ssp bIo 30.18 4.63 1.00 - 6.24 1.00 - 0.06 1.00 - bIo ssp
ssp b94I 23.08 6.67 1.00 - 15.38 1.00 - 0.05 1.00 - b94I ssp
Configuration Sub-Class d Nmax C1 C2 C 3 D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class Configuration
225
Table D.5 Combined mechanism table sorted by Ψ
226
Table D.6 Combined mechanism table sorted by Ψ
Configuration Sub-Class d Nmax C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class Configuration
lpp a90 26.04 a90 lpp
lpp a95 26.57 a95 lpp
lpp a99 26.90 a99 lpp
lpp a 26.96 a lpp
pl p a80 19.11 a80 pl p
pl p a85 21.34 a85 pl p
pl p a90 24.03 a90 pl p
slp a90I 26.28 1.00 4.82 - 1.00 119.69 - 1.00 1.07 - a90I slp
pl p a 27.13 a pl p
slp a95Io 26.95 1.00 1.92 - 1.00 7.04 - 1.00 1.00 - a95Io slp
lps a90 24.60 1.00 - 0.84 1.00 - 0.82 1.00 - 1.37 a90 lps
lps a90IO 24.39 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00 a90IO lps
lps a95I 23.63 0.79 - 1.00 0.50 - 1.00 1.00 - 1.00 a95I lps
lps a99I 23.10 0.68 - 1.00 0.32 - 1.00 1.00 - 1.00 a99I lps
slp a90IO 25.39 1.00 1.00 - 1.00 1.00 - 1.00 1.00 - a90IO slp
slp a99IO 24.97 1.13 1.00 - 1.43 1.00 - 1.00 1.00 - a99IO slp
lps a95IO 24.39 1.00 - 1.00 0.66 - 1.00 0.66 - 1.00 a95IO lps
lps a99IO 24.39 1.00 - 1.00 0.47 - 1.00 0.47 - 1.00 a99IO lps
slp aIo 23.23 1.71 1.00 - 2.54 1.00 - 0.51 1.00 - aIo slp
sss a95I 25.93 0.340 0.289 1.000 0.039 0.022 1.000 0.997 0.895 1.000 a95I sss
sss a99 26.37 1.000 0.745 2.222 1.000 0.435 12.492 1.000 1.052 1.139 a99 sss
sps a99I 26.67 1.00 - 0.52 1.00 - 0.14 1.00 - 0.98 a99I sps
spp a90 26.04 a90 spp
sps a90IO 27.13 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00 a90IO sps
spp a95 26.57 a95 spp
ssp a95I 27.09 2.18 1.00 - 10.43 1.00 - 1.01 1.00 - a95I ssp
ssp a99I 26.85 2.54 1.00 - 11.20 1.00 - 0.68 1.00 - a99I ssp
spp a99 26.90 a99 spp
spp a 26.96 a spp
sss a90IO 26.21 1.260 1.000 3.104 1.974 1.000 1.979 0.986 1.000 0.066 a90IO sss
sps aIO 26.77 1.74 - 1.00 1.31 - 1.00 0.25 - 1.00 aIO sps
ssp a90IO 26.43 1.33 1.00 - 2.00 1.00 - 0.84 1.00 - a90IO ssp
lps aIo 26.77 9.07 - 1.00 6.86 - 1.00 0.01 - 1.00 aIo lps
ssp a95IO 27.12 2.09 1.00 - 2.00 1.00 - 0.22 1.00 - a95IO ssp
sss a95IO 26.71 1.462 1.000 2.693 0.500 1.000 2.000 0.160 1.000 0.102 a95IO sss
sss a 22.82 1.000 1.332 7.099 1.000 1.832 33.744 1.000 0.775 0.094 a sss
ssp aIo 23.23 5.07 1.00 - 7.52 1.00 - 0.06 1.00 - aIo ssp
psp a 27.13 a psp
Configuration Sub-Class d Nmax C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class Configuration
lps b90I 39.44 1.00 - 0.12 1.00 - 0.00 1.00 - 0.98 b90I lps
lpp b90 39.47 b90 lpp
lps b95I 39.66 1.00 - 0.13 1.00 - 0.00 1.00 - 0.73 b95I lps
lpp b95 39.68 b95 lpp
lpp b 39.79 b lpp
pl p b80 36.46 b80 pl p
pl p b85 39.24 b85 pl p
pl p b 40.00 b pl p
slp b90Io 38.58 1.00 1.47 - 1.00 3.14 - 1.00 1.00 - b90Io slp
slp b95IO 36.92 1.00 1.09 - 1.00 1.28 - 1.00 1.00 - b95IO slp
lps b90IO 33.95 1.00 - 1.00 0.68 - 1.00 0.68 - 1.00 b90IO lps
lps b95IO 33.91 0.99 - 1.00 0.60 - 1.00 0.61 - 1.00 b95IO lps
lps b99IO 33.91 0.99 - 1.00 0.55 - 1.00 0.56 - 1.00 b99IO lps
slp bIO 30.18 1.47 1.00 - 1.98 1.00 - 0.63 1.00 - bIO slp
sss b90I 39.35 0.589 0.381 1.000 0.148 0.055 1.000 0.724 0.996 1.000 b90I sss
sps b86IO 40.00 1.00 - 1.00 1.00 - 1.00 1.00 - 1.00 b86IO sps
ssp b90o 39.93 2.15 1.00 - 10.35 1.00 - 1.05 1.00 - b90o ssp
spp b90 39.47 b90 spp
spp b95 39.68 b95 spp
spp b 39.79 b spp
sps bIo 39.60 1.40 - 1.00 1.16 - 1.00 0.42 - 1.00 bIo sps
lps bIo 39.60 7.33 - 1.00 6.05 - 1.00 0.02 - 1.00 bIo lps
sss b99 34.51 1.000 0.936 3.502 1.000 0.995 13.377 1.000 1.215 0.311 b99 sss
sps b98Io 37.47 2.35 - 1.00 6.01 - 1.00 0.46 - 1.00 b98Io sps
ssp b94I 23.08 6.67 1.00 - 15.38 1.00 - 0.05 1.00 - b94I ssp
sss b 32.44 1.000 1.041 4.338 1.000 1.541 19.533 1.000 1.365 0.239 b sss
sss b88IO 39.54 1.610 1.000 2.398 0.500 1.000 2.000 0.120 1.000 0.145 b88IO sss
ssp b90IO 39.74 2.30 1.00 - 2.00 1.00 - 0.17 1.00 - b90IO ssp
ssp bIo 30.18 4.63 1.00 - 6.24 1.00 - 0.06 1.00 - bIo ssp
psp b 40.00 b psp
Configuration Sub-Class d Nmax C1 C2 C3 D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class Configuration
227
228
APPENDIX E GLOSSARY OF TERMS
AND VARIABLES
γζK θ ∆θ
α = -------------------
- . This parameter is used in the stress feasibil-
ρ
ity equation.
γK θ ( R + 1 ) 2 Φ
tion where β = -----------------------------------
ρ
γi ζi Kθ
different flexible segments. κ i = ----------------
i
ρi
229
λ Length parameter used to calculate the actual length from
mechanism dependent.
Sy
Ω = -----
E
1
ρ Is either --- for fixed-pinned beams or µ for small-length
γ
flexural pivots
230
ζ Defines the ratio between the total PRBM length r tot and a
link length
c
where A = -------
r tot
thickness.
Classification Varaibles Those variables that dependent on the mechanism and are
looked up in tables.
231
Configuration Distiguishes between different possible flexible segment
configurations.
Constraint Variable Value Those variable values that are constrained by the design
requirments.
Coupled Variables Those variables that are shared between tow or more of the
deflection.
Force Design Equation Main equation used to calculate the force of a CFM. Knowl-
232
βEI 1
Force Feasibility F = -----------
2
- Equation used to quickly calculate the force
r tot
PRBM length.
Force Feasible Mechanism Those mechanisms that are guaranteed to produce the
Guaranteed Stress Feasbile Those mechanisms that are guaranteed to procude the
achieved.
In-plane Orientation In this orientation, all motion takes place in the plane of fab-
Isolated Varaibles Those variables that occur only in one of the three equa-
tions.
Known Variable Value A variable value that is set by the design requirements.
233
M Multiplier in the power curve fit of α.
same.
developed.
Primary Pivot The flexible pivot in the CFM with the highest stress.
Stress Design Equation Main equation used to determine the stress within a CFM.
tion.
Ω
Stress Feasibility αA ≤ ------- Equation used to determine the stress feasibility
SF
of a given design.
Stress Feasible Mechanism Those mechanisms that are capable of producing the desired
234
other parts of the design space. However, no guarantee is
Unknown Variable Value Those variables values that have no information given by
235
236