Etd 137
Etd 137
Etd 137
by Brent L. Weight
A thesis submitted to the faculty of Brigham Young University in partial fulfillment of the requirements for the degree of
Master of Science
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
___________________________ Date
___________________________ Date
As chair of the candidates 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
Accepted for the College _______________________________________ Douglas M. Chabries Dean, College of Engineering and Technology
ABSTRACT
This thesis adds to the knowledge base of constant-force mechanisms (CFMs). It begins by reviewing past work done in the area of CFMs and then develops new nondimensionalized parameters that are used to simplify the calculations required to design a CFM. Comparison techniques are then developed that utilize these non-dimensionalized parameters to compare mechanisms based on stiffnesses, percent constant-force, actual lengths, normal displacements, and feasible design orientations. These comparison techniques 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 successfully developed a constant-force electrical contact (CFEC).
ACKNOWLEDGEMENTS
This research was supported by funding from the Utah Center of Excellence Program and the National Science Foundation through grant No. DMI-9624574. These foundations are gratefully acknowledged for their help. 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 extends beyond the engineering classroom and laboratory. 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 regards to this work. 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 discovering the affects of friction on constant-force mechanisms is appreciated and 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
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
CHAPTER 2
2.1 Constant-Force Behavioral Model..................................................................... 11 2.1.1 Slider-Crank Model ......................................................................... 12 2.1.2 Principle of Virtual Work................................................................. 13 2.1.3 Non-Dimmensionalization............................................................... 15 2.2 Type Synthesis ................................................................................................... 16 2.3 Original Optimization Problem ......................................................................... 17 2.3.1 Displacement Vector ........................................................................ 18 2.3.2 Objective Function........................................................................... 18 2.3.3 Other Variables ................................................................................ 19 2.4 Original Results ................................................................................................. 19 2.5 Model Verification ............................................................................................. 20 2.6 Model Exceptions .............................................................................................. 21 2.6.1 Eccentricity ...................................................................................... 21 2.6.2 Class 1B Mechanisms ...................................................................... 21 2.7 Outstanding Issues ............................................................................................. 22
CHAPTER 3
Classification Refinement............................................................................ 25
3.1 Original Classification ....................................................................................... 25 3.2 Classification Refinement.................................................................................. 27 3.2.1 Flexible Segment Configuration ...................................................... 27 3.2.2 Sub-Classes ...................................................................................... 28 ix
4.1 Stress Feasibility ................................................................................................ 32 4.1.1 Stress Derivation.............................................................................. 32 4.1.2 Stress Parameters Development....................................................... 37 4.1.3 Stress Feasibility Example............................................................... 41 4.2 Force Feasibility ................................................................................................ 43 4.2.1 Force Derivation .............................................................................. 44 4.2.2 Force Parameter Development......................................................... 45 4.2.3 Relating Moments of Inertia ............................................................ 46 4.2.4 Exceptions in the Force Parameter ............................................... 49 4.2.5 Force Feasibility Example ............................................................... 50
CHAPTER 5
5.1 Length Comparison............................................................................................ 54 5.1.1 Length Parameter for class 1A-lpp .................................................. 54 5.1.2 General values .............................................................................. 56 5.2 Stiffness Comparisons........................................................................................ 56 5.2.1 Stiffness Comparison Issues ............................................................ 56 5.2.2 Comparison Derivation.................................................................... 58 5.2.3 Stiffness Comparisons Results......................................................... 61 5.2.4 Stiffness Comparison Conclusions .................................................. 64 5.3 Percent Constant-Force Comparison ................................................................. 65 5.3.1 Percent Constant-Force Inversion.................................................... 65 5.3.2 Percent Constant-Force Comparison ............................................... 66 5.3.3 Percent Constant-Force Results ....................................................... 69 5.4 Manufacturing Orientations ............................................................................... 71 5.4.1 In-plane Orientation......................................................................... 71 5.4.2 Out-of-plane Orientation.................................................................. 71 5.4.3 Orientation Comparison................................................................... 72 5.5 Flexible Segment Design Space Comparison .................................................... 72 5.5.1 Thickness Limits.............................................................................. 72 5.5.2 Width Limits .................................................................................... 74 5.5.3 Flexible Segment Design Area ........................................................ 76 5.5.4 Determining Fabrication Orientations ............................................. 77 5.6 Normal Displacement Comparison.................................................................... 78 5.6.1 Normal Displacement Derivation .................................................... 78 5.6.2 Normal Displacement Results ......................................................... 80 x
6.1 CFM Model........................................................................................................ 83 6.1.1 General Model ................................................................................. 84 6.1.2 Model Inputs .................................................................................... 85 6.1.3 Model Outputs ................................................................................. 86 6.1.4 Model Verification ........................................................................... 87 6.2 Optimization ...................................................................................................... 87 6.2.1 Objective Functions ......................................................................... 87 6.2.2 Variables........................................................................................... 88 6.2.3 Functions.......................................................................................... 88 6.2.4 Optimization Problem - Stiffer Mechanisms ................................... 89 6.2.5 Optimization Problem - Smaller Normal Displacement.................. 89 6.2.6 Optimization Problem - In-plane/Out-of-plane ............................... 89 6.2.7 Combinations of Optimization Problems ........................................ 90
CHAPTER 7
New Mechanisms......................................................................................... 93
7.1 New Mechanisms............................................................................................... 93 7.1.1 Sub-class Expansion ........................................................................ 94 7.2 Class 1A ............................................................................................................. 94 7.2.1 Optimization Details ........................................................................ 95 7.2.2 Optimization Observations .............................................................. 95 7.2.3 Optimization Results........................................................................ 96 7.3 Class 1B ............................................................................................................. 98 7.3.1 Optimization Details ........................................................................ 98 7.3.2 Optimization Observations .............................................................. 99 7.3.3 Optimization Results...................................................................... 100 7.4 Class 2A ........................................................................................................... 102 7.4.1 Optimization Details ...................................................................... 102 7.4.2 Optimization Observations ............................................................ 103 7.4.3 Optimization Results...................................................................... 105 7.5 Class 2B .......................................................................................................... 108 7.5.1 Optimization Details ...................................................................... 108 7.5.2 Optimization Observations ............................................................ 108 7.5.3 Optimization Results...................................................................... 109 7.6 Class 3A.......................................................................................................... 112 7.6.1 Optimization Details ...................................................................... 112 7.6.2 Optimization Observations ............................................................ 113 7.6.3 Optimization Results...................................................................... 114 xi
8.1 General Design Approach................................................................................ 117 8.1.1 Background, Assumption, and Limitations ................................... 118 8.1.2 Principal Equations ........................................................................ 118 8.1.3 Variable Types................................................................................ 120 8.1.4 Variable Value Types...................................................................... 121 8.1.5 Basic Design Steps......................................................................... 122 8.1.6 Other Flexible Segments................................................................ 124 8.1.7 Final Segment Lengths .................................................................. 125 8.2 Feasible Configurations ................................................................................... 127 8.2.1 Method 1 - Stress Feasibility ......................................................... 128 8.2.2 Method 2 - Force Feasibility.......................................................... 131 8.2.3 Combined Stress and Force Feasibilities ....................................... 133 8.3 Other Mechanism Considerations.................................................................... 134 8.3.1 Percent Constant-Force.................................................................. 134 8.3.2 Flexible Segment Configuration .................................................... 134 8.3.3 Manufacturing................................................................................ 134 8.3.4 Normal Displacement .................................................................... 135 8.4 Variable Manipulation...................................................................................... 135 8.4.1 Variable Relationships ................................................................... 136 8.4.2 General Guidelines ........................................................................ 137 8.5 Design Examples ............................................................................................. 138 8.5.1 Example 1 - Basic Design Case..................................................... 138 8.5.2 Example 2 - Mechanism Elimination - Method 1.......................... 141 8.5.3 Example 3 - Feasible Mechanisms Only ....................................... 145
CHAPTER 9
9.1 Introduction...................................................................................................... 150 9.2 Electrical Contacts ........................................................................................... 152 9.2.1 Electrical Contact Industry Practice and Standards ....................... 152 9.3 Constant-Force Mechanisms as Electrical Contacts ........................................ 154 9.4 CFEC Configurations ...................................................................................... 155 9.4.1 Simulated Pin Joints....................................................................... 156 9.4.2 Parameter Definitions .................................................................... 159 9.5 Case Study Details ........................................................................................... 160 9.6 Model Development ........................................................................................ 161 9.7 Model optimization.......................................................................................... 163 9.7.1 Optimization .................................................................................. 164 xii
9.8 Model Validation.............................................................................................. 165 9.8.1 Dimensional Analysis .................................................................... 167 9.8.2 Testing............................................................................................ 167 9.9 Results.............................................................................................................. 168 9.10 Case Study Conclusions and Recommendations ........................................... 172
CHAPTER 10
10.1 Review of Research Purposes........................................................................ 175 10.2 Summary........................................................................................................ 175 10.3 Conclusions.................................................................................................... 176 10.4 Recommendations for Further Research........................................................ 178 10.4.1 New Configurations ..................................................................... 179 10.4.2 Physical Implementation.............................................................. 180 10.4.3 Other Mechanisms ....................................................................... 180 10.4.4 Other Flexible Segment Configurations ...................................... 180 10.4.5 Applications ................................................................................. 181 10.4.6 Constant-Force Electrical Contacts ............................................. 181
APPENDIX A
A.1 PRBM Overview............................................................................................. 187 A.2 Small-Length Flexural Pivots ......................................................................... 189 A.3 Cantilever Beam with Force at End ................................................................ 191 A.4 Cantilever Beam with End Moment Loading ................................................. 193 A.5 Fixed-Guided Beam ........................................................................................ 194
APPENDIX B
B.1 MatLab Code................................................................................................... 197 B.1.1 File Run Order............................................................................... 197 B.1.2 Variable Name Mapping ............................................................... 199 B.1.3 CFMModel Code........................................................................... 199 B.1.4 UserModel Code ........................................................................... 206 B.1.5 OptdesModel Code........................................................................ 208 B.2 Optimization Code .......................................................................................... 210 B.2.1 anasubC Code ............................................................................... 210 B.3 Optdes/Matlab Link ........................................................................................ 215
xiii
Optimization Plots.................................................................................... 217 Mechanism Tables.................................................................................... 221 Glossary of Terms and Variables.............................................................. 229
xiv
LIST OF FIGURES
Figure 1.1 Figure 1.2 Figure 1.3 Figure 2.1 Figure 2.2 Figure 3.1 Figure 3.2 Figure 3.3 Figure 4.1 Figure 4.2 Figure 5.1 Figure 5.2 Figure 5.3 Figure 5.4 Figure 5.5 Figure 5.6 Figure 5.7 Figure 5.8 Figure 7.1 Figure 7.2 Figure 7.3 Figure 7.4 Figure 7.5 Figure 7.6 Figure 7.7 Figure 7.8 Figure 7.9 Figure 8.1 Figure 8.2 (a) A compliant parallel mechanism and (b) its PRBM counter part .......... 4 Two element chain of constant-force generator mechanisms developed by Nathan (1985) .................................................................. 6 Rigid-body CFM developed by Jenuwine and Midha (1994) ..................... 7 Compliant and rigid-body slider crank model and parameters .................. 12 Fifteen original configurations ................................................................... 17 Fifteen original configurations ................................................................... 26 Refined classification system which identifies flexible segment configurations ........................................................................ 28 Summary of classification scheme ............................................................ 30 Values of for each pivot point ............................................................. 33 Graph of vs. d for Class 1A-spp-a .......................................................... 40 values for sub-class b ............................................................................. 61 Normalized values for sub-class b ......................................................... 62 Normalized values for sub-class a ......................................................... 63 Normalized values for sub-class b without stress equalization ............. 64 Parameters for 'abs calculation ................................................................. 68 (a) In-plane and (b) out-of-plane manufacturing orientations ................... 71 Graphical representation of flexible segment design area ......................... 76 Normal displacement definition ................................................................. 78 Explore plot of R versus 'ex and for (a) Class 1A-lpp-a and (b) Class 1A-spp-a ............................................................................... 96 Explore plot for (a) Class 1A-lpp-b and (b) Class 1A-spp-b ..................... 97 Explore plots for (a) Class 1B-psp-a (b) Class 1B-psp-b .......................... 99 Explore plots for (a) Class 1B-plp-a (b) Class 1B-plp-b ......................... 100 2-D explore plots for slp-a of R and K1 with contours of (a) 'ex and (b) ............................................................................... 103 2-D explore plots for ssp-a with R vs. K1 and contours of (a) ' ex and (b) .............................................................................. 104 Optimum plot for lps-a which shows 'ex, and optimal K2 for given R values ......................................................................... 109 2-D explore plots for sps-a with R vs. K2 and contours of (a) 'ex and (b) ............................................................................... 110 Optimum plot of design variables K1 and K2 versus R for sss-b with curves of 'ex and ......................................................... 113 Interdependencies of coupled primary equations .................................... 120 Flow chart of CFM design steps .............................................................. 124 xv
Figure 8.3 Figure 8.4 Figure 8.5 Figure 8.6 Figure 9.1 Figure 9.2 Figure 9.3 Figure 9.4 Figure 9.5 Figure 9.6 Figure 9.7 Figure 9.8 Figure 9.9 Figure 9.10 Figure 9.11 Figure 9.12 Figure 9.13 Figure A.1 Figure A.2 Figure A.3 Figure A.4 Figure B.1 Figure C.1 Figure C.2 Figure C.3 Figure C.4 Figure C.5
Definition of flexible and rigid segment lengths ..................................... Summary of M value requirements for different types of identified mechanisms ....................................................................... Overlap of guaranteed stress and force feasible mechanisms .................. Summary of variable effects on force and stress magnitudes for an increase in variable magnitudes .................................................... (a) Pogo type connector and (b) cantilever type connector ..................... Typical compression slider-crank constant-force configuration ............. Simulation of pin joints with a circular cam ............................................ Limitations and solutions to limitations of simulated pin joint method .. Selected constant-force electrical connector configuration ..................... Selected CFEC configuration in PDA dock. ........................................... Important parameters for general CFEC. ................................................. Key Points for the finite element model .................................................. CFEC final design .................................................................................... CFEC prototype as compared to a dime .................................................. (a) General testing setup and (b) close-up of contact with fixture and probe ............................................................................................ Graph of force versus displacement from test data .................................. Average and predicted force comparison ................................................ (a) A small-length flexural pivot and (b) its PRBM ................................ (a) Cantilever beam with force at end and (b) its PRBM ........................ (a) Cantilever beam with force at end and (b) its PRBM ........................ (a) Fixed-guided beam and (b) its PRBM ................................................ Summary of Matlab file paths and run order ........................................... 2-D Explore plots for slp-b with R vs. K1 and contours of (a) 'ex and (b) ............................................................................... 2-D Explore plots for ssp-b with R vs. K1 and contours of (a) 'ex and (b) ............................................................................... Optimum plot for lps-b which shows 'ex, and optimal K2 for given R value ........................................................................... 2-D explore plots for sps-b with R vs. K2 and contours of (a) 'ex and (b) ............................................................................... Optimum plot of design variables K1 and K2 versus R for sss-b with curves of 'ex and .........................................................
126 129 133 136 152 155 156 158 158 159 160 163 166 167 169 169 171 189 191 193 194 198 217 218 218 219 219
xvi
LIST OF TABLES
Table 2.1 Table 4.1 Table 4.2 Table 4.3 Table 4.4 Table 5.1 Table 5.2 Table 5.3 Table 5.4 Table 5.5 Table 6.1 Table 7.1 Table 7.2 Table 7.3 Table 7.4 Table 7.5 Table 7.6 Table 7.7 Table 7.8 Table 7.9 Table 7.10 Table 7.11 Table 7.12 Table 7.13 Table 7.14 Table 8.1 Table 8.2 Table 8.3 Table 8.4 Table 8.5 Table 8.6 Table 9.1 Table 9.2 Table 9.3 Table 9.4 Table 9.5 Table 10.1 Original results ........................................................................................... 20 Primary pivot and parameter C for each configuration ............................. 38 Power function values for all classifications ............................................. 40 for each configuration and sub-class ...................................................... 46 Values for Kx-1 and x for each sub-class and configuration ..................... 48 Length parameter formulas and values ...................................................... 57 Percent constant-force values for the original mechanisms. ...................... 70 Maximum thickness and minimum width constants .................................. 73 Summary of manufacturing orientation possibilities ................................. 77 Summary of dNmax values .......................................................................... 81 Summary of different possible optimization problems .............................. 90 Optimization results ................................................................................... 96 Parameter summary for new mechanisms in Class 1A. ............................. 98 Optimization results for Class 1B configurations .................................... 101 Parameter summary for new mechanisms of Class 1B ............................ 102 Optimization results for Class 2A ............................................................ 105 Summary of Class 2A parameters ........................................................... 107 Summary of Class 2A thickness and width ratios ................................... 107 Optimization results for Class 2B ............................................................ 110 Summary of Class 2B parameters ............................................................ 112 Summary of Class 2B thickness and width ratios .................................... 112 Optimization results for Class 3A ............................................................ 114 Summary of Class 3A parameters ........................................................... 115 Summary of Class 3A thickness and width ratios ................................... 115 Largest stiffness parameters and percent increase ................................... 115 Needed values to calculate flexible and rigid segment lengths ............... 127 Summary of stress and force feasibility methods .................................... 130 Example 1 design requirements ............................................................... 139 Final design values ................................................................................... 141 Summary of example 2 design requirements ........................................... 142 Design problem summary for example 3 ................................................. 146 Summary of case-study constraints ......................................................... 161 Parameter summary of final design ......................................................... 166 Dimensional analysis of CFEC prototypes .............................................. 168 Parameter summary of prototype 3 .......................................................... 171 Summary of testing and prediction comparisons ..................................... 172 Mechanisms with highest stiffness .......................................................... 177 xvii
Table B.1 Table D.1 Table D.2 Table D.3 Table D.4 Table D.5 Table D.6
Mapping of parameter names to names in Matlab and OptdesX code .... Combined mechanism table ..................................................................... Combined mechanism table ..................................................................... Combined mechanism table sorted by M ................................................. Combined mechanism table sorted by M ................................................. Combined mechanism table sorted by ................................................. Combined mechanism table sorted by .................................................
xviii
CHAPTER 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 independent of the pseudo-rigid-body model. In parallel with this theoretical work, an attempt to incorporate constant-force behavior into electrical contacts is made. Before new work is done on compliant constant-force mechanisms, general background information is given on constant-force mechanisms and compliant mechanisms. Additionally, a literature review in these areas is presented. Next, an entire chapter is devoted to a review of the work on compliant constant-force mechanisms done prior to 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 is presented, after which conclusions and recommendations are given.
1.2 Background
1.2.1 Constant-Force Mechanisms
Constant-force mechanisms (CFMs) can be rigid-body mechanisms with linear and/or torsional springs or they can be compliant mechanisms. In general, they use the principles of mechanical advantage and stored strain energy to produce a near-constant output force over a large range of input displacement. This is accomplished by determining specific geometric ratios that allow for equal increases in stored strain energy and
mechanical advantage. In this way, the output remains constant throughout the displacement.
(a)
Figure 1.1
(b)
(a) A compliant parallel mechanism and (b) its PRBM counter part
The conversion between rigid-body and compliant mechanisms is done by representing 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 information on the pseudo-rigid-body model see Appendix A.
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 constant-torque spring motors which are capable of producing 50 revolutions on one winding (Wahl, 1963). 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 produce a drive system that maintained a constant normal force on the workpiece while following a prescribed path. Successful drive systems have been developed and demonstrated.
Supported Mass
Figure 1.2
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 mechanism, as seen in Figure 1.3, uses rigid-links and linear springs to achieve a constant-force and has been successfully implemented in concrete testing equipment. Compression slider-crank compliant CFMs have been proposed (Murphy et al., 1994, Howell et al., 1994, Midha et al., 1995). Millar et al. (1996) developed non-dimensionalized 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 configurations. Most recently, Evans and Howell (1999) implemented the compliant CFM into a
r4
r3
springs
r8
r5
rigid link F
r2
r7 r6
r1
Figure 1.3
robot end-effector that successfully demonstrated constant-force behavior while cutting glass. Parkinson et al. (1997) used a parametric optimization approach to develop compliant mechanisms. In one study, they examined a fully compliant constant-force mechanism that was developed through optimization techniques. Herder and Tuijthof (2000), have developed 4 and 6 degrees-of-freedom spatial 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 balanced 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 force generated at the output. 7
Chapter 2 reviews in detail the work done on the compliant slider-crank CFMs. The chapter provides valuable information and a foundation for the proposed research. The notation developed, equations derived, and the optimization problem used are presented and discussed.
sections. Methods have been developed that take this into consideration allowing for accurate calculation of the degrees of freedom (Howell and Midha, 1995).
10
CHAPTER 2
Work was done prior to this thesis on the development of a behavioral model for compliant constant-force mechanisms (CFMs). This behavioral model is used extensively in this work and is used as the basis for further development. All other equations developed in this thesis are derived from this behavioral model or the PRBM and thus have the same accuracy as these two models. 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 (2001) unless otherwise noted.
k2 r3
slider deflection
r2
k3
2
x
F
k1 r1
(b)
Figure 2.1
r1 + r2 r3 2 = acos -------------------------2r 1 r 2 12
(2.1)
r 2 sin 2 3 = asin --------------------r3 r 1 = r 2 cos 2 + r 3 cos 3 These equations allow for all unknown variables to be determined.
(2.2)
(2.3)
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.
(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 energy, ( V ), with respect to the generalized coordinate, ( q ). This results in
dV W = ------ q dq Summing the virtual works in Equation (2.4) to Equation (2.6) results in dVk -------- q k dq k
(2.6)
Having established equations for virtual work, the principle of virtual work can be applied. The principle of virtual work can be stated as (Paul, 1979): 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 written in terms of the generalized coordinate, equation (2.7) reduces to
where A and B are vectors that change the linear and angular displacements into terms of the generalized coordinate. If q k is assumed to be zero (hence the fictitious displacement), 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
Fi A +
Mj B
dVk -------- ( qk ) = 0 dq k
W =
F i z i +
Mj j
(2.7)
(2.8)
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 equations were summed, and the principle of virtual work was applied. The force, F , was solved for resulting in
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 = ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------r 2 r 3 ( sin 2 cos 3 sin 3 cos 2 )
(2.9)
This equation tells how the force, F , is related to the link lengths, spring constants, and angles of the mechanism as defined in Figure 2.1b.
2.1.3 Non-Dimmensionalization
Inspection of the model shows that it relies on many independent variables. It is 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 parameters were chosen. They were r3 R = ---r2 k2 K 1 = ---k1
2
(2.10)
(2.11)
K3 = ----k1
(2.12)
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 equation (2.9). This results in k1 F = ---- r2 where ( R cos 3 [ 2 + K 1 ( 2 + 2 3 ) ] + cos 2 [ K 1 ( 2 + 2 3 ) + K 2 ( 2 3 ) ] ) = -----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------R sin ( 2 3 )
(2.13)
(2.14)
Close examination shows that equation (2.14) is dimensionless. Therefore, F 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 length is known as the geometric parameter. Thus, the creation of non-dimmensionalized parameters reduces the number of independent variables in the model, making the model easier to use.
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)
Figure 2.2
on the number of flexible segments and their location in each configuration. These classifications and configurations can be seen in Figure 2.2.
(2.15)
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 total x would be 4 inches. Originally, two deflection parameters, d, were chosen, 16 (16%) and 40 (40%) deflection. 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 corresponds to the vector r 1 .
(2.17)
18
This parameter shows how constant the static force is for the mechanism throughout the slider displacement. As defined, a perfectly constant-force mechanism would have a level of constant-force equal to 1. Due to the definition in (2.17), is always greater than or equal to 1. Therefore, the optimization objective can be stated Minimize
(2.18)
19
unique solutions to the problems, but represent what was felt to be the best combinations that provided the minimum value.
.
Table 2.1
Original results
K1 1.0000 1.0000 0.1906 0.2237 1.0000 1.0000 K2 0.1208 0.1208 12.6704 9.3816
Class 1A 1B 2A 2B 3A
d 16 40 16 40 16 40 16 40 16 40
R 0.8274 0.8853 1.0000 1.0000 0.3945 0.4323 0.7591 0.8441 2.6633 2.0821
1.0030 1.0241 1.0564 1.1576 1.0015 1.0058 1.0029 1.0235 1.0002 1.0049
0.4537 0.4773 2.0563 2.1513 0.9575 1.0466 0.5230 0.5438 3.4016 3.6286
20
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 from all equations in the behavioral model.
(2.19)
21
It should be noted that Equation (2.19) is used exclusively with the Class 1B mechanisms. With these modifications, the force equation becomes k2 Fin = ---- r2
(2.20)
This configuration was then run through the optimization routine and the results can be found in Table 2.1.
22
Third, the prior work used only one method for comparing mechanisms, the parameter . No attempt was made to determine other methods of comparison or other areas in which improved performance could be achieved. 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 configuration 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 comparison and differentiation between the mechanisms, and eventually lead to an ability to design CFMs for a wide variety of situations. 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 this refined classification, groups of mechanisms or individual mechanism can be referred to quickly and easily.
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)
Figure 3.1
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 3A. Additionally, specific constant-force parameters ( R , K 1 , K 2 , , and ) are associated 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
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
last two pivots. However, notice that the specific set of parameters to be used with the mechanisms has not yet been identified.
3.2.2 Sub-Classes
Sub-classes will be used to distinguish between different sets of constant-force 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 mechanisms are defined. 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 percentage 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 the 40% deflection mechanisms. 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 classification will be Class 2B-lps-c.
29
3A-sss-b
Level Common Elements Applicable
Figure 3.3
Class PRBM
Group of Mechanisms
Sub-Class Parameters
Specific Mechanism
30
CHAPTER 4
To overcome the challenges associated with CFM design and analysis, it is desirable 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 the comparison of different configurations revealing strengths and weaknesses of each. 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 configurations 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 work to make comparisons between different configurations, helping to develop a better understanding of their strengths and weaknesses. Additionally, values for the parameters developed in this chapter will be summarized and tabulated for new configurations developed as part of this work.
(4.1)
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
EI K = K ----l
(4.3)
where is the PRBM characteristic radius factor, K is the stiffness coefficient, E is the materials modulus of elasticity, I is the moment of inertia, and l is the length of the flexible segment. Constant-force mechanism configurations use two different types of flexible segments, small-length flexural pivots, and fixed-pinned beams. There are some assumptions associated with each of these flexible beams. For fixed-pinned beams, the assumptions are: 1. The length of the flexible segment, l , is
ri l = --
(4.4)
33
where is the PRBM characteristic radius factor and r i is an effective link length associated with the flexible segment. For fixed-pinned beams, r i can be either r 2 or r 3 . 2. is typically assumed to be 0.85 .
3. K is approximated as 2.65 For small-length flexural pivots, the following assumptions are made: 1. The flexural pivot length ( l ) is much smaller than the corresponding PRBM link length. Mathematically, 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, or an average link length, r ave , defined below. 2. Commonly, the value for is 0.10. This value will be used unless stated otherwise. 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 ) of r 2 and r 3 . Thus, r i used to find l in Equation (4.5) for k 2 is
r2 + r3 r i = r ave = --------------2
(4.6)
34
Equations (4.4) and (4.5) can be generalized to l = r 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 configuration of the mechanism. With this in mind, a new parameter, , can be developed where r tot = ------ri and r tot = r 2 + r 3
(4.7)
(4.8)
(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 encountered in a CFM are: r2 r3 = R+1 1 = -- + 1 R
(4.10)
(4.11)
r ave = 2
(4.12)
Rearranging Equation (4.8) and substituting Equations (4.2), (4.3), (4.7), and (4.8) into equation (4.1) yields K Ec c = -------------------------r tot 35
(4.13)
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 geometry of the flexible segment. The stress can be related to the safety factor, SF, and the yield strength, S y , as c SF = S y Substituting equation (4.13) into equation (4.14) and rearranging, results in K c Sy 1 ---------------------- = ---- -----r tot E SF
(4.14)
(4.15)
Equation (4.15) can then be separated into a non-dimensionalized stress factor, , a geometric parameter, A , and a material parameter, , where K = ------------------ c A = ------r tot Sy = ---E and equation (4.15) becomes A = -----SF
(4.16)
(4.17)
(4.18)
(4.19)
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.
(4.20)
(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 distance c , and the PRBM length of the mechanism. In determining the stress feasibility of a mechanism, it is necessary to look only at 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 dependent. 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. As an example, take a Class 2A-ssp-a mechanism. The second flexible segment has the highest stress and therefore, the c value for this flexible segment is c p . Addition-
Table 4.1
Configuration Class 1A-spp Class 1A-lpp Class 1B-psp Class 1B-pl p Class 2A-ssp Class 2A-slp Class 2B-sps Class 2B-lps Class 3A-sss
Primary Pivot 1 1 2 2 2 2 3 3 1
C sub-class a sub-class b 4.647 5.090 1.474 1.721 1.405 1.739 7.339 9.081 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, c o < 1.739c p or c o < 0.3478
(4.23)
(4.24)
Equation (4.24) indicates that if c o becomes larger than this value, then the flexible segment 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 classification 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 to be quickly calculated using = Md n where M is the multiplier of the power function and n is the exponent. 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 calculated for each classification at any displacement. It should be remembered that the val-
(4.25)
39
vs d
10 8 6 4 2 0 0 10
Percent Deflection, d
Stress Parameter,
y = 2.351x 2 R =1
0.5004
20
Figure 4.2
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-
Table 4.2
Configuration Class 1A-spp Class 1A-lpp Class 1B-psp Class 1B-pl p Class 2A-ssp Class 2A-slp Class 2B-sps Class 2B-lps Class 3A-sss
40
tain parameters. The next section demonstrates the usefulness and practicality of this method.
(4.26)
41
Therefore, the percent deflection or d for this problem is 12. From Table 4.2, M and n are found to be
M = 3.743 n = 0.511 The parameter can be calculated using equation (4.25) and is = 3.743 ( 12 0.511 ) = 13.31
(4.27) (4.28)
(4.29)
Using the value of c = h 2 for a rectangular cross section, the maximum thickness, h max , can be calculated by combining and rearranging Equations (4.17) and (4.19)
2r tot h max < --------------SF For 1010 Steel, Equation (4.30) becomes 2 ( 0.00087 )5 h max < -----------------------------13.31 ( 1.5 ) or h max < 0.00043 inches and for Beryllium Copper 2 ( 0.0092 )5 hmax < --------------------------13.31 ( 1.5 ) or h max < 0.0046 inches
(4.30)
(4.31)
(4.32)
(4.33)
(4.34)
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 maximum 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 thickness 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 value for h max in Equation (4.34), Equation (4.22) becomes h 1 < 9.08 ( 0.0046 ) = 0.0418
(4.35)
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 values for the exponential relationship for would have to be generated.
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.
(4.36)
where k 1 is the PRBM spring constant for the first pivot, r 2 is a PRBM link length, and is one of the non-dimensionalized CFM parameters. The equation for the PRBM spring constant is EI1 k 1 = K ------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 equation (4.7) into equation (4.37), results in EI1 k 1 = K ------r 2
(4.37)
(4.38)
44
(4.39)
Using the parameter presented in Equation (4.8), Equation (4.39) can be further generalized to 2 K EI 1 F = ---------------------------2 r tot where = R+1 for r 2 . Therefore, substituting Equation (4.41) into Equation (4.40) results in K EI 1 ( R + 1 ) 2 F = ------------------------------------------2 r tot
(4.40)
(4.41)
(4.42)
(4.43)
45
(4.44)
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 be found in Table 4.3. Equation (4.44) is used to determine if a specific configuration, 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 to determine an unknown parameter.
Table 4.3
Configuration Class 1A-spp Class 1A-lpp Class 1B-psp Class 1B-pl p Class 2A-ssp Class 2A-slp Class 2B-sps Class 2B-lps Class 3A-sss
sub-class a 0.4537 15.1508 0.4537 2.9008 2.0563 82.2520 2.0563 15.7482 0.9575 18.6332 0.9575 18.6332 1.2259 37.9347 1.2259 7.2631 3.4016 456.4868
sub-class b 0.4773 16.9649 0.4773 3.2482 2.1500 86.0000 2.1500 16.4658 1.0466 21.4708 1.0466 21.4708 1.2154 41.3322 1.2154 7.9136 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, the general spring constant equation K EI k = --------------l
(4.45)
will be used. Rearranging this equation and substituting values according to the methods above, Equation (4.45) becomes K EI k = -----------------r tot The constant-force parameters K 1 and K 2 , where k2 K 1 = ---k1 k3 K 2 = ---k1
(4.46)
(4.47)
(4.48)
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 constant of any flexible segment can be related to I1 through Equation (4.47) and (4.48). Generalizing, these equations become Ki 1 1 1 K E 1 I 1 i i K E i I i ------------------------------------------ = ------------------------ 1 r tot i r tot
1 i
(4.49)
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
1 i
(4.50)
(4.51)
This parameter depends upon the sub-class and the configuration. Table 4.4 gives the values 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 = ----------------1 Ki 1
(4.52)
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 and stress equations developed to determine feasibilities.
Table 4.4
Configuration Class 1A-spp Class 1A-lpp Class 1B-psp Class 1B-pl p Class 2A-ssp Class 2A-slp Class 2B-sps Class 2B-lps Class 3A-sss
Primary Pivot 1 1 2 2 2 2 3 3 1
K1
K2
K2
48
(4.53)
(4.54)
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 r tot
(4.56)
where I 2 is the first moment of inertia of the flexible segment associated with pivot two ( k 2 ) and is found from Table 4.3.
49
(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 parameter in Table 4.4. The equation becomes 23.17 I1 = ------------------------------- I 3 = 6.86I 3 3.368 ( 1.003 )
(4.58)
50
Substituting the formula for a rectangular cross section and Equation (4.58) into Equation (4.56) and rearranging for b results in
2 12Fr tot b = ----------------------6.86Eh 3
(4.59)
If the original thickness restriction of h > 0.001 inches is used, then the minimum width ( b min ) can be found from
(4.60)
Plugging values into Equation (4.60) results in 12 ( 0.1 )5 2 bmin > -----------------------------------------------------------------6.86 ( 7.263 ) ( 18.5e6 )0.0013 or b min > 32.54 inches
(4.61)
(4.62)
This exceeds the acceptable width and therefore is not an acceptable design. However, 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 > -----------------------------------------------------------------------6.86 ( 3.0987 ) ( 18.5e6 )0.0046 3 or b min > 0.334 inches
(4.63)
(4.64)
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 requirements. It should also be noted that there are several combinations of h and b that will satisfy 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 manner, the complete geometry of all of the flexible segments can be found while still satisfying the constraints of the problem. This example demonstrates the usefulness and ease with which force feasibility checks can be performed.
52
CHAPTER 5
MECHANISM COMPARISONS
To effectively use the CFMs, it is necessary to understand the advantages and disadvantages of each configuration. One way to understand this is to compare the mechanisms 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 understand each configuration including comparisons of:
Actual lengths Maximum stiffness Percent constant-force Manufacturing orientations Normal displacement
53
(5.1)
where is a length parameter. A value for for each mechanism can be found by examining them individually. The calculation of the value for the parameter for the class 1Alpp mechanism is presented as an example. Following which, the general equations and values for of all configurations will be presented.
(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 = --- The second beam is a rigid link with two pin joints and its length is expressed as l 3 = r3 54
(5.3)
(5.4)
(5.5)
To find the value for the parameter L , Equation (5.5) can be substituted into Equation (5.1) and solved for . This results in r2 ---- + r 3 = --------------r tot where r tot = r 2 + r 3
(5.6)
(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 into Equation (5.6) and rearranging results in R + 1 = ------------------( R + 1 )
(5.8)
Equation (5.8) expresses in terms of R and . Substituting the values for R and that correspond to the class 1Aa-lpp mechanism results in 0.8274 ( 0.85 ) + 1 = ---------------------------------------- = 1.027 ( 0.8274 + 1 )0.85
(5.9)
This indicates that the actual length of the mechanism is 1.027 times longer than the PRBM length of the mechanism.
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 Equation (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.
necessary to derive a method to hold conditions such as size and stress constant while comparing the output force. The force equation, Equation (4.44), is EI 1 F = ----------2 r tot
(5.10)
2 where is configuration dependent, I1 and r tot are geometry dependent, and E is mate-
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 ltot
(5.11)
Table 5.1
Value
sub-class b 1.027 1.094 1.0 1.088 1.035 1.088 1.050 1.119 1.050
Configuration
Class 1A-spp Class 1A-lpp Class 1B-psp Class 1B-plp Class 2A-ssp Class 2A-slp Class 2B-sps Class 2B-lps Class 3A-sss
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 underlying issues that have not been addressed. By making the assumptions listed above, the amount of stress found in each configuration at the total displacement has not been taken into account. By holding I1 constant, 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 configuration 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 predetermined stress level is reached.
(5.12)
58
Substituting Equation (5.12) into Equation (5.11), and setting i to the primary pivot, the force equation becomes L 2 E p Ip F = -------------------------2 1 K p 1 ltot
(5.13)
where p , 1 , and K p 1 are configuration dependent, and Ip is geometry dependent. 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 -----SF
(5.14)
To make the comparison, the right hand side of the side of Equation (5.14) is set equal to 1 and the inequality is removed. This results in 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 hLA = ------ = --------l tot 2l tot
(5.16)
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 maximum stresses in two configurations, i and j, can be expressed as i h i L i = j h j Lj 59
(5.17)
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 substituting it into Equation (5.13) gives
3 j i E p i2 b i ---------j h j i i F i = -------------------------------------------------2 12 1 K ( p 1 ) l tot
i 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 length ( l tot ), and width (b). Equation (5.18) can be used to make comparisons between two configurations. Ultimately, comparisons between all of the configurations, not just two, are to be made. This is accomplished by setting all of the configuration b parameters, as well as E, l tot , and ba to 1, and removing the 12. This results in a stiffness intensity parameter where 1 3 p - = --------------------- 1 K p 1
Equation (5.19) no longer results in a force value, but is a dimensionless parameter 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
(5.18)
(5.19)
value
Percent Deflection (d )
Figure 5.1
25
Normalized Values
20 15 10 5 0 0 20 40 60
Percent Deflection (d )
Mech 1A-spp-b Mech 1A-lpp-b Mech 1B-psp-b Mech 1B-plp-b Mech 2A-ssp-b Mech 2A-slp-b Mech 2B-sps-b Mech 2B-lps-b Mech 3A-sss-b
Figure 5.2
The data in Figure 5.1 can be normalized using the maximum value at maximum deflection, c . That is, a normalized parameter value is expressed as = ----c
(5.20)
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 better understanding of the force behavior. For example, the graph shows that the force output rises from a value of 1 at 16% deflection to 22 at 2% deflection, an increase of 2200%. This information helps show the trade off between deflection and stiffness of the mechanisms, 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
25
Values
20 15 10 5 0 0 5 10 15 20
Percent Deflection (d )
Mech 2A-ssp-a Mech 2A-slp-a Mech 2B-sps-a Mech 2B-lps-a Mech 3A-sss-a
Normalized
Figure 5.3
example, mechanism 1B-plp-b has only 40% of the force or stiffness of mechanism 1Alpp-b when stress levels and deflections are equivalent. The same procedures were followed for the mechanisms in sub-class a. The normalized 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. To illustrate the importance of equating the stresses, normalized values when 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
Non-Equivalent Stress
Normalized Stress (Percent of Maximum Force) 1.2 1 0.8 0.6 0.4 0.2 0 0 20 40 60 Percent Deflection (d )
Figure 5.4
Mech 1A-spp-b Mech 1A-lpp-b Mech 1B-psp-b Mech 1B-plp-b Mech 2A-ssp-b Mech 2A-slp-b Mech 2B-sps-b Mech 2B-lps-b Mech 3A-sss-b
maximum deflection. However, this information is misleading because the stresses in the mechanism are not considered.
64
nisms that are stiffer than the mechanisms currently defined. The parameter also makes it quick to determine if a mechanism is stiffer than another. 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 indicating what deflection range must be used to obtain a desired force.
(5.21)
To help this parameter be more intuitive, it can be inverted and multiplied by 100, resulting in min ( F ) ' = 100 -----------------max ( F )
(5.22)
Multiplying this parameter by a hundred gives the percent constant-force as a percentage 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 and maximum output force of a CFM.
(5.23)
range at the beginning in which friction and other factors distort the output force. In these cases, to ensure comparability, ' ex must be determined through another method. 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 ---------------------------- + F 0 d max d min
(5.24)
67
Fmax
Force (F)
Fmin F0
dmin
Figure 5.5
dmax
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
F0 = Fmax
68
(5.25)
(5.26)
' ex
This equation gives a level of constant-force for the entire deflection of the mechanism. If the minimum force ( F min ) is already on the y-axis, then Equation (5.27) reduces to Equation (5.22).
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 intuitive way as 99% constant rather than the value 1.0030. 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
Fmin 1 ----------F0 Fmax = 100 ----------- = 1 -------------------F max d min 1 ---------d max
(5.27)
(5.28)
(5.29)
a greater variation in force than the other mechanisms. This knowledge will allow designers 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. The parameter ' ex is useful in establishing a consistent manner to measure percent 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 constantforce '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 usefulness of this parameter to determine the percent constant-force from experimental data will be exhibited as testing results are explored.
Table 5.2
1.0030 1.0564 1.0015 1.0721 1.0002 1.0241 1.1576 1.0058 1.1914 1.0049
' ex
99.70 94.66 99.85 93.27 99.98 97.65 86.39 99.42 83.93 99.51
70
(a)
Figure 5.6
(b)
be fabricated from a uniform piece of material with either simple milling or stamping type operations.
(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 presented in this work. These minimum values, C , can be determined as
no )
(5.31)
By selecting the minimum values, the inequality in Equation (5.30) will be satisfied 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 segments is the minimum thickness of the mechanism. This limitation is defined by the limitations of the manufacturing process and materials.
Table 5.3
Configuration Class 1A-spp Class 1A-lpp Class 1B-psp Class 1B-pl p Class 2A-ssp Class 2A-slp Class 2B-sps Class 2B-lps Class 3A-sss
Primary D equal D min D equal D min C C Pivot 1 1 2 2 7.522 6.242 0.062 2 5.090 0.057 4.647 1.721 2.543 0.499 1.474 0.618 2 1.980 1.739 1.314 0.250 1.405 1.157 0.417 3 9.081 6.861 0.009 7.339 5.103 0.013 3 1 1.384 8.258 1.832 33.745 0.691 0.060 1.043 4.359 1.541 19.533 1.359 0.236
73
(5.32)
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 Equation (4.52), this relationship can be defined as Ki 1p I i = ----------------- I p K p 1 i Furthermore, using Equation (5.32), Equation (5.33) can be expanded to Ki 1p 3 b i h i3 = ----------------- b p h p K p 1 i
(5.33)
(5.34)
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. The minimum width is also controlled by limitations on the thicknesses and, in 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 are all equal, Dequal. 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 bi Ki 1 p ---- = ------------------------ = D min bp Kp 1 i C3
(5.35)
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 primary flexible segment. 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 the thicknesses, Equation (5.34) becomes bi Ki 1 p ---- = ----------------- = D equal bp Kp 1 i
(5.36)
This ratio indicates what ratio the widths of the flexible segments must be when the thicknesses are equal for all flexible segments.
75
1 bx bp
Dmin
Feasible Design Space Limited by Manufacturing Process Limited by Mechanism Direction to move for Improved Design Area
hx hp
Figure 5.7
76
Table 5.4
1
NA NA
77
k2 r3
slider deflection
normal deflection
r2
k3
2
x
F
k1 r1
Figure 5.8
(5.37)
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 can be described in terms of 2 as y = r 2 sin 2 where, by using the law of cosines,
2 2 2 r1 + r2 r3 ------------------------- 2 = acos 2r 1 r 2
(5.38)
Using Equations (4.8), (4.10), and (4.11), r 2 and r 3 can be related to r tot by
79
dr 1 = r tot 1 -------100
(5.39)
(5.40)
(5.41)
(5.42)
Substituting Equations (5.40), (5.41), and (5.42) into Equation (5.39), and rearranging, results in
where r tot drops out of the equation. Substituting Equations (5.40) and (5.43) into Equation (5.38) gives r tot 1y = ------------ sin acos -------R+1 200
Equation (5.45) gives the normal displacement as a percentage of the total displacement 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 parameter help compare different mechanisms, but it can be used during optimization to look for better mechanisms.
80
2 = acos
(5.43)
(5.44)
(5.45)
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 for any d can be found using Equation (5.45).
Table 5.5
Mechanism 1a 1b 2a 2b 3a
sub-class a d Nmax R 0.8853 26.96 1.0000 27.13 0.4323 23.23 0.8441 26.77 2.0821 22.82
sub-class b d Nmax R 0.8853 39.79 1.0000 40.00 0.4323 30.03 0.8441 39.60 2.0821 32.44
81
82
CHAPTER 6
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 of CFMs connected with this research.
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.
(6.1)
84
where ( R cos 3 [ 2 + K 1 ( 2 + 2 3 ) ] + cos 2 [ K 1 ( 2 + 2 3 ) + K 2 ( 2 3 ) ] ) = -----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------R sin ( 2 3 ) or k2 F = ---- r2 where ( R cos 3 [ K 1 ( 2 + 2 3 ) ] + cos 2 [ K 1 ( 2 + 2 3 ) + K 2 ( 2 3 ) ] ) = -----------------------------------------------------------------------------------------------------------------------------------------------------------------------------R sin ( 2 3 ) for Class 1B mechanisms. As before, a displacement vector is generated from zero displacement to the maximum deflection, dr tot .
(6.2)
(6.3)
(6.4)
Flexible pivot types: pin, small-length, or long flexible Percent deflection: d Associated link for Class 1-B mechanisms
The inputs that have no direct affect on the configurations, but are required to analyze actual designs are:
Material Properties: E, Sy Spring Constant for the first flexible pivot: k1 Total Length: ltot Width of Flexible segments: b
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.
Force Parameter:
6.2 Optimization
The optimization is performed to determine the best combination of constant-force parameters that provides the most desirable CFM performance. The performance is measured based on the comparison methods described in Chapter 5.
87
on the desired output. Different optimization problems are defined below, at which point the objective functions are clearly stated. 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 values for the desired mechanism. The only variables that are used as design variables are the constant-force parameters. These parameters allow the optimization to modify the link length ratio, R, and the spring constant ratios, K, in search of mechanisms with more desirable performance.
6.2.3 Functions
The analytical functions for the optimization problem include all of the model outputs described above. The selection of analytical functions as design functions depends upon the objective of the optimization problem. For each optimization problem described below, the design functions are clearly outlined.
88
Maximize
(6.5)
(6.6)
Minimize
(6.7)
(6.8)
89
criteria defined in Section 5.5.4 and tabulated in Table 5.4. Formally, the optimization problem is written as: '
Maximize
(6.9)
subject to
C1 0.5 D equeal 2
and/or
D min 1
(6.10)
(6.11)
Table 6.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 identify the particular problems used in the search for better mechanisms within each configuration, as well as the results achieved.
91
92
CHAPTER 7
NEW MECHANISMS
This chapter presents the results of the optimization, including explore plots and optimum plots, increases in stiffness, and summarization of parameters.
7.2 Class 1A
Class 1A was the first class of mechanisms examined. This class consists of both the lpp and spp mechanisms. 94
95
.0777 .0629
1.96
99.3
1.59
93.9
1
2-psi .0335 .0482 2-psi .848 1.22 88.6 88.6
93.9
99.3
1-Xi
83.2
1-Xi 83.2
1 2 1 2
1 2 1 2
.0188 .00403
.477
77.8
.107
72.4
.500
1.00
1.50 R
2.00
2 2.50
72.4
77.8
3.00
.500
1.00
1.50 R
2.00
2 2.50
3.00
(a)
Figure 7.1
(b)
Explore plot of R versus 'ex and for (a) Class 1A-lpp-a and (b) Class 1A-spp-a
Table 7.1
Optimization results
Configuration lpp
lpp
spp
spp
Sub-Class a a99 a95 a90 b b99 b95 b90 a a99 a95 a90 b b99 b95 b90
1.000 1.046 1.239 1.490 1.000 N/A 1.052 1.129 0.039 0.041 0.049 0.059 0.039 N/A 0.041 0.044
R 0.8274 0.8018 0.7106 0.6185 0.8853 N/A 0.8595 0.8226 0.8274 0.8018 0.7104 0.6187 0.8853 N/A 0.8595 0.8226
% 4.6% 23.9% 49.0% 5.2% 12.9% 4.9% 24.9% 50.8% 5.1% 13.0%
%'ex -0.7% -4.7% -9.7% -4.7% -9.7% -0.7% -4.7% -9.7% -2.7% -7.8%
96
.0996 .0810
95.3
2.51
78.5
2.04
78.5
95.3
2-Psi .0624
61.8
2-Xi
1-Xi' 61.8
.0438
2 2
1 1
45.1
.0251
28.4
.640
28.4
45.1
11.6
.172
11.6
.00651
.500
.500
(a)
Figure 7.2
(b)
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 about 4% of the total stiffness of the lpp-a. The Class 1A-lpp-b and Class 1A-spp-b mechanisms had similar trends as the subclass a mechanisms. However, at 90% constant, these mechanisms only showed a 13% increase in stiffness. Additionally, no feasible mechanism was found for both configurations 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 prior work. 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 prevent 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
This class of CFMs consists of the psp and plp configurations.
Table 7.2
Configuration Sub-Class Primary Pivot lpp a a99 1 a95 a90 lpp b b99 1 b95 b90 spp a a99 1 a95 a90 spp b b99 1 b95 b90
98
.0151
.0373 .0328
94.6
12
.0146
92.9
82.5
86.2
1 2
91.2
78.8
1-Xi
89.4
1 2
.0132
.0193
87.7
.0128
.0148
85.9
.500
1.00
1.50 R
2.00
2.50
3.00
(a)
Figure 7.3
67.8
(b)
99
.710
94.6
86.2 82.5
.574
90.3
1.47
1.80
1 2-Psi .303 .439 86.0 78.8 2-Xi 81.6 1 2 71.5 .167 77.3 .449 75.2 1-Psi .788 1.13
1-Xi
2 2 .0312 67.8 73.0 1.00 2.00 R 3.00 2 4.00 .110 5.00 1.00 1 1 2.00 R 3.00 2 1 4.00 5.00
(a)
Figure 7.4
(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 constantforce.
100
For the Class 1B-plp-a mechanisms, feasible mechanisms for percent constantforce 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 constant-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 mechanisms is found in Table 7.4
Table 7.3
Configuration Sub-class psp a a99 a95 a90 psp b b99 b95 b90 pl p a a99 a95 a90 a85 a80 pl p b b99 b95 b90 b85 b80
'ex 94.6 N/A N/A N/A 86.3 N/A N/A N/A 94.6 N/A N/A 90 85 80 86.3 N/A N/A N/A 85 80
0.015
R 1.0000
% -
%'ex -
0.016 N/A N/A N/A 0.378 N/A N/A 0.644 0.698 0.710 0.409
1.0000 N/A N/A N/A 1.0000 N/A N/A 0.4388 0.3170 0.2530 1.0000
0.500 0.596
0.7919 0.6043
22.1% 45.5%
-1.5% -7.3%
101
7.4 Class 2A
The Class 2A mechanisms consist of the slp and ssp configurations. These mechanisms are combinations of the Class 1A-spp and Class 1B-psp and pl3p configurations.
Table 7.4
pl p
102
20.0 .00991
.233
.341
.420
.420 .341
.264 .233
16.0
20.0
16.0
93.3
94.5
94.5
94.0
93.3
92.2
91.1
90.1
88.0
12.0
86.4
K1
8.00
K1
4.00 80.1 .500 1.00 1.50 R 2.00 2.50 3.00
4.00
84.3
8.00
12.0
.500
1.00
1.50 R
2.00
2.50
(a)
Figure 7.5
(b)
2-D explore plots for slp-a of R and K1 with contours of (a) 'ex and (b)
103
.0148
.0145
.0140
.0135
.0128
.0122
12.0
66.6 K1 8.00
K1 8.00
12.0
.0118
69.7 75.4
85.4
85.4
80.1
75.4
69.7
68.4
67.6
67.1
.0118
65.6 4.00
4.00
.00973
(a)
Figure 7.6
(b)
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 mechanism becomes infeasible. The 2-D explore plots for slp-b showed similar relationships among the parameters. These plots can be found in Appendix C under Figure C.1. The ssp-a and ssp-b configurations were examined in the same way as the slp configurations. 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. However, the contours for are different. In this mechanism, the area in which is the highest corresponds to the area with the highest ' ex . Additionally, in this area, as K1 decreases, increases, indicating that K1 should be minimized.
104
The ssp-b explore plots, shown in Figure C.2 in Appendix C, display the same trends as those found in the ssp-a configuration.
Table 7.5
Configuration Sub-Class slp asI a99IO a95Io a90I a90IO slp bIO b99 b95IO b90Io ssp asI a99I a95I a95IO a90IO ssp bsI b99 b94I b90IO b90o
'ex 99.8 99 95 93 97.7 99.8 N/A 96 90 99.8 99 96.2 95 90 99.4 N/A 94 90 91.2
0.145 0.219 0.338 0.534 0.233 0.181 N/A 0.294 0.331 0.017 0.042 0.048 0.022 0.028 0.019 N/A 0.030 0.022 0.052
R 0.3950 0.5057 0.8237 1.5278 0.5437 0.4323 N/A 0.6248 0.7267 0.3950 0.7864 0.9182 0.9563 1.4688 0.4323 N/A 0.3000 0.8714 0.9325
K1 0.1906 0.2640 1.6370 15.0000 0.3521 0.2237 N/A 0.3924 0.8283 0.1906 0.1000 0.1000 0.5113 0.4050 0.2237 N/A 0.1000 0.5343 0.1000
% %'ex 51.4% -0.8% 133.6% -4.8% 268.7% -6.8% 61.0% -2.1% 62.9% 83.4% 141.4% 176.9% 24.3% 59.5% 54.5% 14.7% 170.8% -3.8% -9.8% -0.8% -3.6% -4.8% -9.8% -5.4% -9.5% -8.2%
105
resents a fairly large increase in stiffness, with only a marginal decrease in percent constant-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 increase of 83% in stiffness. 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 orientation 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 subclass has a ' ex value of 96.2 and a 176.9% increase in stiffness over the original subclass. However, forcing the optimization to look for sub-classes suitable for in-plane and out-of-plane orientations resulted in improvements of 24% for 95 percent constant-force and 60% increase for 90 percent constant-force. 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 constantforce. However, a sub-class, not suitable for in-plane orientation and only possibly suitable for out-of-plane fabrication orientation, was found for 90 percent constant-force with a 170% increase in stiffness.
106
In all, the results of the optimization are promising. Many new sub-classes are presented in Table 7.5. The important design parameters for these new sub-classes are summarized in Table 7.6 while the thickness and width ratios determined from the optimization are summarized in Table 7.7.
Table 7.6
Configuration slp Sub-Class Primary Pivot aIo 2 a99IO 2 a95Io 1 a90I 1 a90IO 2 bIO 2 b99 b95IO 1 b90Io 1 aIo 2 a99I 2 a95I 2 a95IO 2 a90IO 2 bIo 2 b99 b94I 2 b90IO 2 b90o 2 'ex 99.8 99.0 95.0 93.3 97.7 99.4 96.2 90.0 99.8 99.0 96.2 95.0 90.0 99.4 R 0.3950 0.5057 0.8237 1.5278 0.5437 0.4323 0.6248 0.7267 0.3950 0.7864 0.9182 0.9563 1.4688 0.4323
slp
K1 0.1906 0.2640 1.6370 15.0000 0.3521 0.2237 0.3924 0.8283 0.1906 0.1000 0.1000 0.5113 0.4050 0.2237
ssp
ssp
0.030 1.038 6.3534 0.5372 23.08 0.022 1.027 5.5051 0.5168 39.74 0.052 1.026 5.4994 0.5165 39.93
Table 7.7
Configuration slp
slp
ssp
ssp
Sub-Class Primary Pivot aIo 2 a99IO 2 a95Io 1 a90I 1 a90IO 2 bIO 2 b99 b95IO 1 b90Io 1 aIo 2 a99I 2 a95I 2 a95IO 2 a90IO 2 bIo 2 b99 b94I 2 b90IO 2 b90o 2
107
7.5 Class 2B
The lps and sps configuration make up the 2B class of CFMs. The original configurations in this class of CFMs exhibit the lowest values for percent constant-force of all the original mechanisms.
108
14.1
.508
99.8
11.3
8.52
3-K2
95.9
97.9
5.72
93.9
2 1 3
2.91
.121
92.0
3
.0246 .100 90.0
.500
1.00
1.50
1 R
2.00
2.50
3.00
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 generated for the Class 2A configurations were generated for the sps-a configuration. The highest ' 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.
109
.00837
.00551
.00387
.00284
.00216 4.50
93.7
93.7
87.6
84.2
80.8
8.00
76.2
6.00
K2
70.6 4.00
K2
2.00
.500
1.50
2.50 R
3.50
4.50
2.00
4.00
6.00
8.00
.500
1.50
2.50 R
3.50
(a)
Figure 7.8
(b)
2-D explore plots for sps-a with R vs. K2 and contours of (a) 'ex and (b)
Table 7.8
Configuration Sub-Class lps aIo a99I a95I a90 a99IO a95IO a90IO lps bsI b99IO b95IO b90IO b95I b90I sps aIO a99I a90IO sps bsI b98Io b86IO
0.026 0.241 0.253 0.286 0.173 0.211 0.278 0.035 0.222 0.235 0.252 1.053 1.137 0.028 0.063 0.057 0.037 0.031 0.062
R 0.7591 2.5750 2.4050 2.0960 2.1623 2.1623 2.1623 0.8441 1.9335 1.9339 1.9290 0.8561 0.8177 0.7591 0.7328 1.0000 0.8441 0.6528 1.0000
%'ex K2 % 1.0029 6.4256 826.5% 6.2% 4.3690 873.1% 1.9% 2.0313 1000.0% -3.4% 5.1740 565.4% 6.2% 3.6510 711.5% 1.9% 2.4156 968.5% -3.4% 1.0230 4.9463 539.8% 18.3% 4.5230 576.4% 13.5% 3.9830 626.2% 7.5% 0.0100 2934.2% 13.5% 0.0100 3177.2% 7.5% 1.0029 0.1845 123.8% 6.7% 1.0022 104.6% 1.7% 1.0235 0.2549 -17.2% 17.6% 1.0000 68.6% 3.2%
110
.00169
.0418
.0418 .0270
.0180
.0119
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 unconstrained by manufacturing orientations. The lps-b configurations had increases in stiffness from 539% to 3177%. This represents a significant improvement. In fact, the b95I and b90I sub-classes have higher stiffness 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, five new sub-classes were added to this configuration. The sps-a configuration showed moderate increases in stiffness of about 100%. The level of constant-force was increased from 93 to 99 while stiffness for the same subclass 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 and in-plane orientation, 95 percent constant-force, and an increase of 104% in stiffness. 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 original value of 84 to a high of 98 percent constant-force was achieved with sub-class b98Io . This is a substantial increase in the level of constant-force. 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.
Table 7.9
'ex 93.1 99.0 95.0 90.0 99.0 95.0 90.8 83.7 99.0 95.0 90.0 95.0 90.0 93.1 99.3 94.6 83.7 98.4 86.3
Configuration lps
lps
sps
sps
Sub-Class aIo a99I a95I a90 a99IO a95IO a90IO bIo b99IO b95IO b90IO b95I b90I aIO a99I a90IO bIo b98Io b86IO
Primary Pivot 3 3 3 1 3 3 3 3 3 3 3 1 1 3 1 3 3 3 3
R 0.7591 2.5750 2.4052 2.0960 2.1623 2.1623 2.1623 0.8441 1.9336 1.9339 1.9292 0.8561 0.8178 0.7591 0.7328 1.0000 0.8441 0.6528 1.0000
Table 7.10
Configuration lps
lps
sps
sps
Sub-Class aIo a99I a95I a90 a99IO a95IO a90IO bIo b99IO b95IO b90IO b95I b90I aIO a99I a90IO bIo b98Io b86IO
112
.0865
99.9
2 1
95.9
97.9
1-Xi 94.0
.0194
92.0
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.
113
Table 7.11
Configuration sss
sss
K1 K2 % %'ex 1.0000 12.6700 0.3497 8.3820 291.7% -1.0% 1.7047 2.9708 5.4% -5.0% 0.4140 15.3420 317.8% -4.4% 0.3969 0.6459 52.5% -9.0% 1.0000 9.3816 0.6930 7.1497 13.7% -0.5% 0.3329 5.4511 157.8% -3.2% 1.8187 3.3349 -19.8% -11.3%
114
Table 7.12
Configuration Sub-Class Primary Pivot 'ex a 1 100.0 sss a99 1 99.4 a95IO 2 95.0 a95I 3 95.6 a90IO 2 91.1 b 1 99.5 sss b99 1 99.0 b90I 3 96.3 b88IO 2 88.4 R 2.6633 1.4903 1.3464 1.6573 1.5518 2.0821 1.8709 1.2420 1.1990
K1 K2 1.0000 12.6700 0.3497 8.3820 1.7047 2.9708 0.4140 15.3423 0.3970 0.6459 1.0000 9.3816 0.6930 7.1497 0.3330 5.4510 1.8187 3.3349
Table 7.13
Configuration Sub-Class Primary Pivot a sss 1 a99 1 a95IO 2 a95I 3 a90IO 2 b sss 1 b99 1 b90I 3 b88IO 2
Table 7.14
Configuration Sub-Class 'ex 99.7 lpp a 99 a99 95 a95 90 a90 97.6 lpp b 95 b95 90 b90 lps b95M 95 b90M 90
% K2 R 0.8274 0.8018 4.6% 0.7106 23.9% 0.6185 49.0% 0.8853 0.8595 5.2% 0.8226 12.9% 0.8561 0.0100 5.3% 0.8177 0.0100 13.7%
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 subclasses with 16% deflection belong to Class 1A-lpp. Two of the 40% deflection mechanisms 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 additional 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
CFMs have the potential of being used in a wide variety of applications. To be used in these applications, the CFMs must first be designed to meet all of the applications requirements. An understanding of a suitable approach to design is desired to simplify the process used by designers. This chapter looks at the design aspect of CFMs including general design procedures, infeasible mechanism elimination methods, secondary issues, and trends between variables. The chapter ends by presenting different design example problems and solutions.
117
(8.1)
Substitution of Equation (4.52) (moment of inertia ratios), Equation (5.1) (length parameter equation), and the equation for the moment of inertia for a rectangular cross-section into Equation (8.1) results in
3 2 p b p h p E F = -------------------------------2 12ltot 1 K p 1
(8.2)
This equation will be referred to as the force design equation. The only exception to Equation (8.2) is for Class 1B mechanisms. The force design equation for this class of mechanisms is 2 EI 2 F = ----------------2 l tot
(8.3)
The second principal design equation is a modified form of the stress feasibility equation introduced in Equation (4.19) in Chapter 4 as A = -----SF hSubstitution of = Md n and A = --------- into Equation (8.4) results in 2l tot h------ = Md n --------2l tot SF This equation will be referred to as the stress design equation.
(8.4)
(8.5)
119
Displacement Equation
disp
ltot
Sy
SF
ltot
Figure 8.1
The final principal equation is a combination of the percent deflection equation found in Equation (2.16) and Equation (5.1), the length parameter equation. The third equation is dl tot disp = ----------100
(8.6)
This equation is an intermediate equation which helps to provide needed values for variables 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 to determine a design. A graphical representation of the interdependencies of these equations are given in Figure 8.1.
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 isolated 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 number of possible classifications, and an iterative process may be required. Methods for eliminating infeasible mechanisms will be discussed later.
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 additional values are needed to solve the principal equations.
Stress and Force Feasibility (Section 8.2) Percent constant-force (Section 8.3.1) Flexible segment configuration (Section 8.3.2) Manufacturing orientations such as in-plane and out-of-plane (Section 8.3.3) Normal Displacement (Section 8.3.4)
These points can be used to reduce the number of mechanisms by eliminating any mechanism that is not feasible, greatly simplifying the selection process.
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 purpose. 3. Add constraint values, starting with the most important ones, until the equations 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 direction in which the constraint value should be adjusted to achieve the desired performance. At this point, the process will become iterative until a suitable design is chosen. 4. Use unknown variables as free variables if needed. If, after adding all constraint variables, there are still not enough variable values to solve the equations, use reasonable assumptions for unknown values until there are enough to solve the system. 5. Remove variables if system is over-constrained. If, at any point, the system is over-constrained, remove the least important values first. 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 will be given in Section 8.5.
123
Yes
Yes
No
Figure 8.2
(8.7)
This equation relates the moment of inertia of the primary pivot with the moment of inertia of any other flexible segment.
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 thicknesses) or in-plane (equal widths) orientation is used. However, if either one of these orientations, 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 stress levels and satisfy Equations (8.8) and (8.9).
125
f2 f1 f3 r2 r3
Figure 8.3
1. Calculate r tot using the length parameter 2. Calculate r 2 , r 3 , and r ave where (from Equations (4.8) to (4.12)) r tot r i = ------ and for r2 r3 = R+1 1 = -- + 1 R
(8.11)
(8.12)
(8.13)
r ave = 2 3. Use Table 8.1 to calculate the flexible segment lengths 4. Use the relations in Table 8.1 to calculate the lengths of the rigid sections 5. Replace flexible segments of length zero with pin-joints
(8.14)
126
Table 8.1
To get: Do spp lpp psp plp ssp slp sps lps sss
f1 Multiply r2 by 0.1
0 0 0.1
1--------0.85
0 0 0.1 0.1 0.1
1--------0.85
1--------0.85
0.1
127
Since all 3 principal equations contain classification variables, some approximations must be made. For this reason, the results of these methods are only as accurate as the approximations and must be verified.
128
Mechanism M value
Stress Infeasible
Mmax
Stress Feasible
Mmin
Figure 8.4
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 may only be necessary to identify stress feasible mechanisms. To use this method, the following must be satisfied:
In the deflection equation, 2 out of the 3 isolated and coupler variables must have known or constraint values. 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 = 1.12 for Mmin and = 1 for Mmax. 2. Calculate the desired M value. Calculate the M value that corresponds with the type of mechanism to be identified. Table 8.2 defines which constraint values to use to calculate the different M values. 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 identified. Table D.3 in Appendix D contains all of the mechanisms and their parameters sorted by M.
Table 8.2
Stress Feasiblity - Method 1 M min M max minimum maximum maximum minimum minimum maximum maximum minimum maximum minimum minimum minimum Calculated Calculated
maximum for desired type minimum for desired type of of configurtaion or 0.5256 configuration or .4699 maximum for desired type minimum for desired type of of configuration or 1.12 configuration or 1.00 Mechanisms that can undergo deflection without violating design constraints Full Constrained Design Space Partial Constrained Design Space
If Force is feasible
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 approximations made during the calculations. After choosing a mechanism, the actual n, M, and values should, together with the stress design equation, be used to verify that the constraints have not been violated. 5. Finish the design. Once mechanisms have been eliminated, a mechanism must be chosen and the design process continued.
In the deflection equation, 2 out of the 3 isolated and coupler variables must have known or constraint values. 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 variable values that correspond to the lpp configuration. 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 the lpp configuration. 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. 4. Eliminate mechanisms. Any sub-class with the same deflection percentage 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 eliminated from consideration. Table D.5 in Appendix D contains all of the mechanisms and their parameters sorted by . 5. Determine thickness. The actual thickness, which will be different than the thickness in step 1, must be determine. This can be accomplished by using the classification variable values that correspond to the chosen mechanism in the force design equation. 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 con132
Force Feasible
Method 2
Figure 8.5
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 required to calculate min .
133
8.3.3 Manufacturing
The types of manufacturing processes that can be used to make a given mechanism are summarized in the sub-class name. If a certain process is to be used, a configuration suitable for this process must be identified and used. 134
(8.15)
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 subclass a, 40 for sub-class b). For deflection percentages below the maximum value, Equation (5.45) can be used. For design purposes, the normal displacement information can be used to determine 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 during the basic design. The location of their use in the basic design steps depends on the type of variable values that result from Equation (8.15).
Sy
SF
h 3
ltot 2
disp -
M -
L -
Effect on Stress
Key n n Decrease, Increase in Magnitude by power n (default is 1) Decrease, Increase in Magnitude Limit
Figure 8.6
Summary of variable effects on force and stress magnitudes for an increase in variable magnitudes
136
Adjust thickness to change the force while minimizing the change in stress Adjust length to change stress while minimizing the change in force
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 design equation. Therefore, generally, it is best to maximize the length of the mechanism for the reasons stated directly above 4. Minimize safety factor. Use the smallest possible value rather for the safety factor. This allows for the largest design space. 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 the stress, and not the stress or force itself. 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 can be achieved without switching flexible segment configurations. However, a decrease in ' may be necessary. 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.
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 constantforce for the lpp configuration. Additionally, solving the deflection equation for the maximum 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. 2. Fill in known values This example has four known values for the three principal equations. Examination of Figure 8.1 indicates that more variable values are needed. 3. Fill in constraint values 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-
Table 8.3
139
tions. Using the maximum length results in a value for d of 29.3. This value can also be added to the design equations. 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 the geometry and mechanism constraints.
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 iterations or optimization methods.
Table 8.4
141
The first basic design step is to choose a configuration. In this particular problem, a specific mechanism can not be chosen. Therefore, it is necessary to eliminate as many mechanisms as possible using one of the two methods described in Section 8.2. 2. Determine which elimination method can and should be used The deflection equation is missing one of the coupler variables and the classification 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 variable b. This indicates that only method 1 can be used for this problem. 3. Calculate the percent displacement In this case, we are looking for a fully compliant mechanism. All of these mechanisms have a value of 1.05. Using this value and Equation (8.6) results in ) 0.0075 ( 1.05 ) disp ( d = 100 ------------------ = 100 ------------------------------- = 5.25 l tot 0.15 4. Calculate Mmax to identify feasible mechanisms
Table 8.5
(8.16)
Units mN Pa Pa mm mm mm mm -
Low 1.5 90 -
High 2.0 -
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. 2l totmax ( 2 )5.02 10 3 ( 0.15 ) Mmax = ------------------------------------------ = ----------------------------------------------------------------- = 0.3792 1.1 ( 1.05 ) ( 0.0015 ) ( 5.25 0.5 ) SFmin h min d n
a p p ro x
(8.17)
5. Use Table D.3 to eliminate configurations Examination of Table D.3 shows that the smallest M value for any sss configuration 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 conclusion 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 compliant mechanism without changing any of the other design constraints? 1. Choose a mechanism The sss configuration with he lowest M value was identified above. This mechanism also has the highest value of all of the sss-a mechanisms. Therefore, this mechanism, the sss-a95I, will be used. 2. Calculate the needed length A non-linear method must be used to calculate the needed length due to the interdependencies 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 more than double the initial length value. Question: What thickness would be required to be able to use a fully compliant mechanism and all of the original design constraints? 1. Choose a mechanism The sss-a95I mechanism will be used for the same reasons listed above. 2. Calculate the maximum thickness for the primary pivot Using the appropriate values for M (1.7796) and n (.4888) and Equation (8.5), the maximum thickness can be calculated as 2l tot 5.02 10 3 ( 0.15 ) h max = --------------------- = 2 ------------------------------------------------------------------- = 0.326 mm ( 1.05 )1.7796 ( 5.25 0.4888 )1.1 Md n SF
(8.18)
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 would be required to achieve the desired force? 1. Calculate the required width to achieve the force using the force design equation, Equation (8.2) The equation is
2 12Fl tot 1 K p 1 12 ( 0.003 )0.15 2 ( 26.57 )15.3423 b p = ------------------------------------ = -------------------------------------------------------------------------------------------------------- = 2.15 mm 3 2 p h p E 1.05 2 ( 16.03 ) ( 0.000326 3 )2.279 ( 1.1 10 11 )
(8.19)
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 rectangular cross section can be used to calculate the thicknesses of the other beams. These thickness are 0.111 mm and 0.091 mm for the first and second flexible segments respectively. Question: Is the above design feasible? The only thing that has not been verified are the stresses in the other flexible segments. 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 suitable for in-plane orientation, this inequality must be satisfied.
145
The mechanism was not identified in the problem. However, the value for d, using a value of 1, is found to be 25. This indicates that a 40% mechanism (b sub-classes) must be used. To identify all of the possible mechanisms, both the stress and force feasibilities must be used. 2. Use force feasibility method to elimination mechanisms 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 generate the needed force without violating the design constraints (assuming the same width and lengths are used).
Table 8.6
Low 0.01 -
High 0.035 -
146
The application of this method reduces the number of sub-class b mechanisms down to 5 mechanisms, 3 lpp configurations and 2 lps configurations. 3. Use method 1 to identify guaranteed 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 Appendix D, 11 guaranteed mechanisms can be identified. 4. Determine overlap of two methods Method 1 resulted in 11 guaranteed mechanisms being identified, 5 of which were identified also using method 2. Therefore, any one of these 5 mechanisms (lps-b90I, lppb90, lps-b95I, lpp-b95, and lpp-b) can be used to satisfy the design requirements of this example.
147
148
CHAPTER 9
The objective of this chapter is to demonstrate the potential benefits and viability of applying constant-force mechanism technology to electrical contact design. The successful 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 CFECs into industry applications. The chapter begins by discussing current electrical contact industry practices and standards. It then presents a discussion of constant-force mechanism technology. Different configurations are explored to discover one suitable for application as a CFEC in a presented case study. (The CFEC case study presented in this chapter is focused on limitations 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 meets the requirements of the case study.
149
9.1 Introduction
The reliability of high-cycle electrical connectors is of great concern to designers, and methods to improve this reliability are always being evaluated. According to Deshpande 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 minimal contact resistance under all types of use conditions and accommodate expected geometric variations in manufacture and assembly. The factor that contributes most significantly to the reliability of electrical contacts is the contact-surface mating conditions. Two physical parameters that greatly affect mating conditions are contact surface finish, and contact normal force at mating. When contact surface finish remains corrosion free, either by being corrosion resistant or by being self-cleaning, greater reliability is achieved. When contact normal force is maintained 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 geometric variations. Thus a desirable contact system would maintain an optimal contact force regardless of variations in assembly or use. In addition to achieving high levels of reliability, the electrical contact industry is 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 designs that can be produced at low cost (Brush Wellman, 1999). The recent introduction and advancements in design of simple constant-force 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 disassociation of contact normal force and contact deflection. Traditional mechanics describe force (F) and deflection (d) of springs as that represents the stiffness of the connector. A CFEC uses constant-force technology to separate the contact normal force and the contact deflection, resulting in a relationship where force is relatively independent of 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 dismissed or undiscovered can be explored. The disassociation of contact normal force and contact deflection could lead to three advantages. First, current electrical contacts require tight manufacturing and assembly tolerances to ensure that the contact deflection is within acceptable limits. The decoupling 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 manufacturing and assembly. The second advantage may be in applications were the user interacts with the electrical contacts, such as docking stations. In this case, the decoupling of the force and 151 = kd where k is the spring constant
Undeflected Contact
Housing Springs
Housing
Figure 9.1
(a) (a) Pogo type connector and (b) cantilever type connector
(b)
deflection helps prevent any variation introduced by the user, such as different docking methods, from affecting the performance of the electrical contacts. 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 system environment. Examples of this include connectors in aircraft, vehicles, and machinery.
of these are divided into subgroups for presentation here, but it is not intended to be an exhaustive list of all design issues. 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 commonly used to achieve higher yield strengths. Unfortunately, they are more difficult to use and more expensive than phosphor bronze. Manufacturability is an important aspect of electrical contact design. Electrical 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:
Minimum Material Thickness - There is a minimum material thickness that is suitable for stamping. Minimum Bend Radius - There is a minimum bend radius allowed during stamping operations. A general rule used is 4 times the thickness of the material. 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 suitable for use. Design for assembly is vital to achieve a low cost and reliable part. Some of the practices associated with assembly are:
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). Single Assembly - Assembly of all of the contacts into the housing at one time is the standard method. This greatly simplifies the process. 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, examination of traditional contacts shows that they are not suitable for use as constant-force mechanisms, requiring that a new configuration be developed for the CFEC.
Manufacturability - The stamping of the necessary geometry would be difficult. Material - The deflections and size constraints would cause extremely high stresses compared to the strengths of common electrical contact materials. Assembly - The assembly of pin-joints makes the use of traditional slider-crank configurations in electrical contacts unlikely.
154
Displacement
Figure 9.2
Electrical Continuity - Pin joints would introduce gaps and areas of high resistance in the electrical path making the contact inefficient and unreliable.
Evaluation of the latest configurations shows that pin-joints and small-length flexural 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.
155
Displacement
Figure 9.3
156
If the flexible beam is loaded, but doesnt slide around the cam, then the beam could buckle. Slider crank change points should be avoided. At these points, it may be difficult 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 Figure 9.4a and Figure 9.4b. 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. To combine the strengths of constant-force mechanisms and bent beam electrical contacts, many different possible configurations were evaluated. Using the industry practice and standards criteria and a screening process, the configuration determined most viable 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 inside of a PDA dock is shown in Figure 9.6. This configuration is easy to manufacture and assemble, and has electrical continuity. Additionally, the beam and cam combination provide the necessary increases in
157
Displacement
Displacement
Eccentric
Cam
Cam
Displacement
Displacement
Displacement
Cam
Cam
Figure 9.4
Deflection
Cam
Figure 9.5
158
kM Doc
n otio
le kab Do c
t Un i
Docking Station
Figure 9.6
mechanical advantage and the strain energy storage device necessary for constant-force behavior.
159
Deflection R1
Cross Section
h package
l2 3
l3 b
l1
l4 R2
R3 l5
l6 b package
R4
Figure 9.7
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 column 3 lists the actual constraint values for the case-study.
Table 9.1
Summary of case-study constraints General Constraint Symbol hpc bpc hc bc Rc Ec Syc SFc Fave, c Fmin,c Fmax, c Constraint Value for Case-Study
hpackage bpackage h b Bend Radius E Sy SF Fave Fmin Fmax Percent Constant (')
6 mm 12 mm
0.2 mm 1.0 mm
0.7 mm
110e9 Pa 552e6 Pa
1.0
'c
161
tion phase, the model is called many times to calculate function values and derivatives. Therefore a simple model is preferred, but it must also be accurate. Originally, a combination of a numerical CFM model based on the PRBM and a numerical bent beam model based on Eulers Method were combined to model the CFEC configuration. However, the model was not accurate enough due to violations of assumptions 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 actually 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 between the FEA and optimization programs. 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 175 beam elements are used to model the CFEC. 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 given a large width and height to ensure that it is rigid.
162
100 C 4 20 D B 21 F E 7 A 1 J 22 23 I 11 9 H 10 5
6 2 3
8 G
12
Figure 9.8
It is also necessary to constrain segment J at key point 12 in the x and y directions, 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 the top of the mechanism at key point 100. 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 by the optimization software (Optdes-X).
mechanism from the layout presented in Figure 9.5 that satisfies all of the design constraints. 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 values 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 functions calculated from variables and other model results.
9.7.1 Optimization
Using the design constraints of the case study, the optimization problem can be formally stated in as: Maximize ' Subject to the constraints: ' > ' c F ave > F min
c
(9.1)
(9.2)
(9.3)
(9.4)
c
(9.5)
(9.6)
(9.7)
164
With the following constraints on the design variables: R 1 > Rc R 2 > Rc R 3 > Rc R4 > Rc
(9.8)
(9.9)
(9.10)
(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.
165
4.7672 2.9398
78 .92 R0
x2
[9 5 ]
[10]
.7 R0
R0 .7
1.0
0.4184
0.4956
1.35
0.5087
2.72
5.1221
5.3912
Figure 9.9
Table 9.2
Parameter summary of final design Final Design 5.9 mm 5.4 mm 0.2 mm 1.0 mm 0.7 mm 110e9 Pa 552e6 Pa 1.29 478 mN (48.8 gf) 423 mN (43.2 gf) 577 mN (58.9 gf) 73.2 Constraint Value for Case-Study
Model Function Symbol hpackage bpackage h b Bend Radius E Sy SF Fave Fmin Fmax Percent Constant (')
2.04
.5 R0
4.5521
.7 R0
6 mm 12 mm
0.2 mm 1.0 mm
0.7 mm
110e9 Pa 552e6 Pa
1.0
166
Figure 9.10
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 affected the performance of the prototype. A force transducer was attached to a computer-controlled actuator. The computer 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 maintained a near constant-force throughout the deflection and behaved as predicted.
Table 9.3
Dim
Design
Weight Value 5.401 2.925 1.490 0.183 0.199 0.700 0.738 90.000 2.287 5.334 0.394 0.475 0.537 96.100 10.300 0.952
A 5.2809 3 B 2.9116 1 C 1.45 1 D 0.2 2 d2 0.2 2 E 0.8 1 F 0.8 1 G ang 90 2 G 2.24 1 H 5.5188 3 I 0.4184 3 J 0.6 1 K 0.8 1 L 95 2 M 10 2 N 1.0278 1 Total Weighted Variation Rank
3 Variation 0.068 0.005 0.028 0.170 0.010 0.125 0.078 0.000 0.021 0.100 0.175 0.208 0.329 0.023 0.060 0.074 1.473 1
168
(a)
Figure 9.11
(b)
(a) General testing setup and (b) close-up of contact with fixture and probe
1000
800
400
200
0 0.000 -200
0.100
0.200
0.300
0.400
0.500
0.600
0.700
0.800
0.900
Displacement (mm )
Figure 9.12
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 mechanisms. This is easily addressed by applying a deflection preload of 0.05 mm. 169
The second phenomenon observed is a difference in force between the compression 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 the expansion stroke. 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 within the mechanism and testing system. 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 friction can be removed by averaging the compression and expansion strokes . 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 predictions are listed in Table 9.4 while Figure 9.13 shows the predicted and adjusted-measured forces for prototype 3.
170
600 500 Force (mN) 400 300 200 100 0 0 0.2 0.4 0.6 0.8 Predicted Polypropylene Teflon
Deflection (mm)
Figure 9.13
Table 9.4
Parameter summary of prototype 3 Prototype 3 5.6 mm 5.5 mm 0.2 mm 1.0 mm 0.7 mm 110e9 Pa 552e6 Pa 1.29 448 mN (45.7 gf) 418 mN (42.6 gf) 500 mN (51.0 gf) 79.4 Constraint Value for Case-Study
Function Symbol hpackage bpackage h b Minimum Bend Radius E Sy SF Fave Fmin Fmax Percent Constant (')
6 mm 12 mm
0.2 mm 1.0 mm
0.7 mm
110e9 Pa 552e6 Pa
1.0
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 percent-constant force. Table 9.5 contains a summary of the comparison between the testing and predicted values.
Table 9.5
172
The successful demonstration of the uncoupling of the force and deflection in an electrical contact promises to create new possibilities in electrical contact designs, possibilities in lowering required manufacturing tolerances, reduction of system sensitivity to variations introduced by the user, and increased system robustness in applications where movement and/or vibrations exist. 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 understood as well as a better understanding of the phenomenon observed.
173
174
CHAPTER 10
Seek to understand compliant constant-force mechanism behavior from both a stress and force viewpoint. Develop different ways of comparing constant-force mechanisms. Define new mechanisms with improved performance. Outline a design method that can be used independent of the pseudorigid-body model. Incorporate constant-force behavior into electrical contacts.
10.2 Summary
The purposes listed above were attained through the main contributions of this research including:
Non-dimensionalized parameters were developed. These parameters can be used to assess stress and force feasibilities.
175
The classification system for CFMs was refined and expanded to include naming methods for groups of mechanisms, families of mechanisms, and individual mechanisms. Methods for comparing constant-force mechanisms were established and defined. A detailed optimization model was developed that incorporated the new non-dimensionalized parameters. 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. A design methodology, including detailed examples, was outlined and set forth to assist designers in the design of CFMs. Included were the primary design equations, methods to identify stress and force feasible mechanisms, and description of variable interactions and trends. A constant-force electrical contact was developed and tested. The CFEC performed as intended and the details and results were presented.
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 comparison 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 constant-force increased from 84 to 99, an 18% increase. 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 mechanisms 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 pivots. 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-
Table 10.1
' Configuration Sub-Class 99.7 lpp a 99 a99 95 a95 90 a90 97.6 lpp b 95 b95 90 b90 lps b95M 95 b90M 90
% K2 R 0.8274 0.8018 4.6% 0.7106 23.9% 0.6185 49.0% 0.8853 0.8595 5.2% 0.8226 12.9% 0.8561 0.0100 5.3% 0.8177 0.0100 13.7%
177
ology allows the designer to make educated decisions and choices based on the information 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 that will satisfy all of the design constraints. 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 configurations were developed, and one configuration was chosen. From this configuration, a design was generated, prototyped, and tested. Testing results showed that the CFEC displayed constant-force behavior and characteristics similar to the models predictions. Additionally, the successful demonstration of the uncoupling of the force and deflection in an electrical contact promises to create new possibilities in electrical contact designs, possibilities in lowering required manufacturing tolerances, reduction of system sensitivity to variations introduced by the user, and increased system robustness in applications where movement and/or vibrations exist. This could lead to benefits in terms of performance, robustness, and package size.
178
Cam- Replacement of link with circular cam Offset - Mechanisms that utilize the eccentricity of the mechanism Orthoplaner - Configurations that start in an orthoplaner (in plane) position Linear Spring - Configurations that include a linear spring Perpendicular Force - The movement of the force from a horizontal position to a vertical position Inverted Mechanisms - These are mechanisms that work under tension Series/Parallel configurations - Different arrangements of CFM in series and parallel to achieve desired displacement/force Pre-loads - Different pre-loads on the torsional and linear springs 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 presented by these types of mechanisms.
179
Sharp rise as the mechanism is initially displaced Difference in force between compression and expansion strokes Sensitivity to tolerances
The effects of physical implementation of CFMs should be studied more fully. These phenomenon and others should be studied, allowing CFMs to be used in future applications.
ment without affecting stress, should exhibit an even greater stiffness than the lpp configuration. It is recommended that further work should be done to understand the behavior of the lpl configuration, its stiffness, and its susceptibility to stress.
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 (2001) on the dynamics of CFMs to develop CFMs for dynamic applications.
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. Boyle, C. L. , 2001, A Closed-Form Dynamic Model of the Compliant ConstantForce Mechanism using the Pseudo-Rigid-Body Model, M.S. Thesis, Brigham Young University, Provo, Utah. Brush Wellman Inc., 1999, Connector Engineering Design Guide: Material Selection in the Design of Spring Contacts and Interconnections, Brush Wellman Inc. Chang, L. H., and Fu, L. C., 1997, Nonlinear adaptive control of a flexible manipulator for automated deburring, Proceedings - IEEE International Conference on Robotics asnd Automation, Vol. 4, pp. 2844-2849. Deshpande, A. and Subbarayan, G., 2000, LGA Connectors: An Automated Design Technique for Shrinking Design Space, Journal of Electronic Packaging, Vol. 122, pp. 247-254. Evans, M. S., and Howell, L. L., 1999, Constant-Force End-Effector Mechanism, Proceedings of the IASTED International Conference, Robotics & Applications, Santa Barbara, CA, USA, pp. 250-256. Harper, C. A. 1996, Electronic Packaging & Interconnection Handbook -Second Edition, McGraw-Hill, New York. Herder, J.L., Tuijthof, G.J.M., Two Spatial Gravity Equilibrators, Proceedings of the 2000 Design Engineering Technical Conferences and Computer in Engineering Conference, DETC2000/MECH-14120. Herder, J. L., van den Berg, F. P. A., Statically Balanced Compliant Mechanisms (SBCMs), An Example and Prospects, Proceedings of the 2000 Design Engineering Technical Conferences and Computer in Engineering Conference, 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 DETC94, ASME Design Engineering Technical Conferences, DETC98/MEMD-71. Howell, L. L. and Midha, A., 1994, The Development of Force-Deflection Relationships for Compliant Mechanisms, Proceedings of the 1994 ASME Mechanisms Conference, DE-Vol. 71, pp. 501-508. 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 MultipleOutput 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, DEVol. 70, pp. 61-66. Nathan, R. H., 1985, A Constant Force Generation Mechanism, Journal of Mechanisms, Transmissions, and Automation in Design, Trans. ASME, Vol. 107, December, pp. 508-512. 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 University, Provo, Utah. Williman, J., 1995, Small Torque, Engineering (London), Gillard Welch Associates, 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 different 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).
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 beginning ideas to refined designs. 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 of compliant mechanism in displacement, force, velocity, and acceleration. Thus, designers 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 presented. The strain energy stored in each flexible member is represented by a torsional 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 each of these spring constants will be given.
188
nP
nP
P
l 2
nP l
Figure A.1
nP l L
Torsional Spring
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 in each flexible element will be given.
(A.1) (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 ( EI ) K = -----------l l T = K The x and y coordinates of the end of the beam can be found through
(A.3)
(A.4)
(A.5)
l l a = -- + L + -- cos 2 2 and
(A.6)
l b = L + -- sin 2 The stress equations for this flexible segment are ( Pa + nPb )c nP top = ---------------------------------- -----I A ( Pa + nPb )c nP bot = ------------------------------ -----I A
(A.7)
(A.8)
(A.9)
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
nP
nP
nP L
Figure A.2
nP L
(1-)L
Torsional Spring
(a) Cantilever beam with force at end and (b) its PRBM
191
(A.13)
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 where c 1.24
(A.14)
(A.15)
The value of c, the parametric angle coefficient, varies between 1.256 and 1.179. Values for different values of n are tabulated in Howell (2000). The torsional spring constant for a cantilever beam with a force on the end can be found from EI K = K ----l K 2.65
(A.16)
(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 ) 192
(A.18)
M
o
M
L
Torsional Spring
L
(1-)L
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.20)
193
Centerline
Torsional Spring
L
(1-)L 2 (1-)L 2
Figure A.4
The parametric angle coefficient is c = 1.5164 and the stiffness coefficient is K = 2.0643
(A.21)
(A.22)
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 The spring constant is found to be EI K = 2K ----l
(A.24)
(A.25)
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 fixedguided is four times that of the cantilever beam. Since the end of the beam is constrained, a reactionary moment, Mo is created. The formula for this moment is Pa M o = -----2 or Pl M o = ---- [ 1 ( 1 cos ) ] 2
(A.26)
(A.27)
195
The maximum stress occurs at the ends of the beam where the moment is largest. It has a value of Pac max = --------2I
(A.28)
196
APPENDIX B
197
Independent of Optimization
Linked to Optimization
Run Order
UserModel
OptdesLink
OptdesModel
CFMModel
Figure B.1
The UserModel file is used by the user to run the model independent of the optimization 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 values from CFMModel. The file OptdesModel is similar in purpose to UserModel. However, this file 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 create an output file that is read by the optimization code. 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 finished and the output file has been generated.
198
Table B.1
Variable Name Matprop(1) Matprop(2) Matprop(3) SF stress(x) lengths(1) lengths(2) lengths(3) lengths(4) geoparam(x,1) geoparam(x,2) geoparam(x,3) pivotparam(x,1) pivotparam(x,2) pivotparam(x,3) pivotparam(x,4) pivotparam(x,5) pivotparam(x,5) cfparam(1) cfparam(2) cfparam(3) cfparam(4) cfparam(5) cfparam(6) cfparam(7) configparam(1) configparam(2) configparam(3) configparam(4) configparam(5) configparam(6) ri r2 r3 r4 Alternate Name
fp dN * note: in AnasubC, array values are all shifted down by 1due to C code requirments.
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
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)); 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
%Calcu-
%Forces Orthoplaner position by changing r4 from user defined to ideal for orthoplaner
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 ri=sqrt((r2+r3)^2-r4^2)*(1-di); %Calculate initial displacement (rest) if (ri<rimin) %& (d<0) %Force ri to be no longer than maximum allowable displacement with given geometery ri=rimin; di=1-ri/(r2+r3); if d>0 %Don't allow any compression because ri=rimin d=0; end elseif (d>0) & (ri*(1-d)<rimin) %checks to see if compression mechanism is going too far d=1-rimin/ri; %Resets percd to maximum allowable end %Solve for theta i's if (di==0) & (r4==0)%if initially flat theta2i=0; theta3i=2*pi; theta4i=pi; elseif (di==1) %if orthoplaner theta2i=pi/2; theta3i=3*pi/2; elseif (r4==0) theta2i=acos((ri*ri+r2*r2-r3*r3)/(2*ri*r2)); %Solves for initial angles (not pre-load) theta3i=2*pi+asin(-r2*sin(theta2i)/r3); %Solves for initial angle 3 (not pre-load) else del = 1e-6; tol = 1e-3; theta2i = Newton(type,'theta2',thetaT,[ri,r2,r3,r4],limit,del,tol); theta3i = asin((r4-r2*sin(theta2i))/r3)+2*pi; end
theta2o=theta2i+pivotparam(1,3); %Initial angle + Pre-load-->Pre-load is Counter Clock-wise resulting in negative torque (radians) theta3o=theta3i+pivotparam(2,3); %Initial angle + Pre-load-->Pre-load is Counter Clock-wise resulting in negative torque (radians) theta4o=2*pi-theta3i+pivotparam(3,3); %d is the percent deflection of mechanism. if (d>0) %If d is >0, then is percent of rest length. deltax=0.001:0.001*ri*d:(ri)*d; %Set up percentage vector [r1]=ri-deltax; else % if d <0, then is percent of total allowable(total length - rest length) deltax=0.001:0.001*(sqrt((r2+r3)^2-ri^2)):(sqrt((r2+r3)^2-r4^2)-ri)*d*-1; [r1]=ri+deltax; %Set up r1 vector end [trash,sz]=size(deltax); %Find value of lengths
202
%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
%Equation to solve for non-dimensionalized parameter phi % Note: Equation for 1B mechanism is not the same as for others. if class1B==0
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
%Equations to calculate alpha alpha(1,i)=gamma(pivotparam(1,1))*zeta(1)*Ktheta(pivotparam(1,1))*(theta2(i)-theta2o)/roe(1); alpha(2,i)=gamma(pivotparam(2,1))*zeta(2)*Ktheta(pivotparam(2,1))*(theta3o+theta2(i)-theta3(i)-theta2o)/ roe(2); alpha(3,i)=gamma(pivotparam(3,1))*zeta(3)*Ktheta(pivotparam(3,1))*(theta3o-theta3(i)-theta4o)/roe(3); %Equations to calculate stress stress(1,i)=alpha(1,i)*A(1)*matprop(1); stress(2,i)=alpha(2,i)*A(2)*matprop(1); stress(3,i)=alpha(3,i)*A(3)*matprop(1); %Equations to calculate which pivot has the highest stress if (stress(1,i)>=stress(2,i)) & (stress(1,i)>=stress(3,i)) if (maxstresspivot==1) maxstresspivot=1; else maxstresspivot=1; changepivot=changepivot+1; end elseif (stress(2,i)>stress(1,i)) & (stress(2,i)>stress(3,i)) if (maxstresspivot==2) maxstresspivot=2; else maxstresspivot=2; changepivot=changepivot+1; end elseif (stress(3,i)>stress(2,i)) & (stress(3,i)>stress(2,i)) if (maxstresspivot==3) maxstresspivot=3; else maxstresspivot=3; changepivot=changepivot+1; end 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
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Return important parameters back to the interace function. pivot=maxstresspivot; kapivot=ka(pivot) ka1=ka(1) Kpivot=K(pivot) 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; %Average array values
205
%Multiply in the average Beta value for i=1:sz fparam(i)=fparam(i)*average_Beta; end 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
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)
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); %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)=cfc_Data1(4,2); %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)
[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) save cfc_Results.txt lengths geoparam pivotparam cfparam configparam averageforce di d SF pc ka C DesignParam -ASCII %Result File Setup %ri r2 r3 r4 %b1 h1 l1 %b2 h2 l2 %b3 h3 l3 %pinjoint1 0 theta2pre k1 theta2i theta2final %pinjoint2 Reference theta3pre k2 theta3i theta3final %pinjoint3 0 theta4pre k3 0 0 %0 0 0 ks 0 0 %R K1 K2 K3 phi phis xi %alphamax pivot changepivot average_Beta fp dN %averageforce %di %d %SF %m,n %ka1 ka2 ka3 %C1 C2 C3 %hparam1 bparam1 D1 Area1 %hparam2 bparam2 D2 Area2 %hparam3 bparam3 D3 Area3
209
210
/*File declarations for linking to other models */ FILE *out, * in, *flag; /* char str[80];*/
/* get AV values from OptdesX (Variable names 16 chars max) */ /* be sure to use the ADDRESS of the variables in the function calls */ /*Checks to see if Hold existes */ flag=fopen("Hold.txt","r"); if(flag!=NULL) {fclose(flag); system("rm Hold.txt"); } fclose(flag); 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(&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"); /*Create Data file for input into Matlab */ fprintf(out,"%g %g %g %g %g -1 -1\n",E,Sy,ltot,pivotparam[0][3],geoparam[0][0]); fprintf(out,"%g %g %g %g %g %g %g\n",cfparam[0],cfparam[1],cfparam[2],cfparam[3],di,d,ds); fprintf(out,"%g %g %g %g %g -1 -1\n",pivotparam[0][0],pivotparam[1][0],pivotparam[2][0],r4type,r4); fprintf(out,"-1 %g -1 -1 -1 -1 -1\n",pivotparam[1][1]); fprintf(out,"%g %g %g -1 -1 -1 -1\n",pivotparam[0][2],pivotparam[1][2],pivotparam[2][2]); fprintf(out,"%g -1 -1 -1 -1 -1 -1\n",fit);
fclose(out); /*Close matlab data file */ hold=1; /*system("matlab ,main.m &"); */ /*Force anasubC to hold till Matlab is finished*/ 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"); /*Read in data from Result file*/ in = fopen("cfc_Results.txt","r"); /*variable=getValue(in,#));*/ 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];
/* send functions to OptdesX (Function names 16 chars max) */ afdscaC(averageforce, "Average Force" ); afdscaC(cfparam[0], "R" ); afdscaC(cfparam[1], "K1" ); afdscaC(cfparam[2], "K2" ); afdscaC(cfparam[3], "K3"); afdscaC(cfparam[4], "Phi"); afdscaC(cfparam[5], "Phi S"); afdscaC(cfparam[6], "Xi" ); afdscaC(dN,"dN"); afdscaC(alphamax,"Max. Alpha"); afdscaC(pivot,"max S P"); afdscaC(changepivot,"pivot changes"); afdscaC(averagebeta,"average Beta"); afdscaC(fp,"f"); afdscaC(k1,"k1"); afdscaC(k2,"k2"); afdscaC(k3,"k3"); afdscaC(ks,"ks"); 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"); afdscaC(di,"% Defl. Initial"); afdscaC(d,"d"); afdscaC(ds,"% Spring Pre-Load"); afdscaC(SF,"Safety Factor"); afdscaC(m,"Power m"); afdscaC(n,"Power n"); afdscaC(L,"L"); afdscaC(b1,"b1"); afdscaC(h1,"h1"); afdscaC(l1,"l1"); afdscaC(b2,"b2"); afdscaC(h2,"h2");
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(hp1,"h param 1"); afdscaC(bp1,"b param 1"); afdscaC(D1,"D1"); afdscaC(Area1,"Design Area 1"); afdscaC(hp2,"h param 2"); afdscaC(bp2,"b param 2"); afdscaC(D2,"D2"); afdscaC(Area2,"Design Area 2"); afdscaC(hp3,"h param 3"); afdscaC(bp3,"b param 3"); afdscaC(D3,"D3"); afdscaC(Area3,"Design Area 3"); 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 { FILE *out, * in, *flag; out = fopen("cfc_postdata1.txt","w"); /*variable=getValue(in,#));*/ fprintf(out,"%lf %lf %lf %lf 0 0 0\n",geoparam[0][0],geoparam[0][1],geoparam[0][2],geoparam[0][3]);
214
%lf
%lf
fprintf(out,"%lf %lf %lf 0 0 0 0\n",pivotparam[0][0],pivotparam[0][1],pivotparam[0][2]); fprintf(out,"%lf %lf %lf 0 0 0 0\n",pivotparam[1][0],pivotparam[1][1],pivotparam[1][2]); fprintf(out,"%lf %lf %lf 0 0 0 0\n",pivotparam[2][0],pivotparam[2][1],pivotparam[2][2]); fprintf(out,"%lf %lf %lf 0 0 0 0\n",pivotparam[3][0],pivotparam[3][1],pivotparam[3][2]); fprintf(out,"%lf %lf %lf 0 0 0 0\n",pivotparam[4][0],pivotparam[4][1],pivotparam[4][2]); 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"); }
215
% Insert first function here x=1; optdesmodel(x); save Hold.txt temp -ASCII end
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 been written to the data file. 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 for the file Hold to be created. 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
K1 8.00
K1 8.00
.049
65.8 4.00
.034 4.00
.018
62.7 1.00 2.00 R 3.00 4.00 5.00
1.00
2.00 R
3.00
4.00
5.00
(a)
Figure C.1
(b)
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
.022 .024
.027
.031
12.0
12.0
65.6 4.00 4.00 62.4 1.00 2.00 R 3.00 4.00 5.00 1.00 2.00 R 3.00 4.00
(a)
Figure C.2
(b)
2-D Explore plots for ssp-b with R vs. K1 and contours of (a) 'ex and (b)
20.0
.987 .792
16.0
95.7
97.4
12.0
2-Psi .597
93.9
3-K2
8.01
.403
1-Xi
92.1
3 2
4.01
.208
90.4
3 2
.0172
.0133
88.6
.500
Figure C.3
Optimum plot for lps-b which shows 'ex, and optimal K2 for given R value
218
71.3
79.3 84.2
93.6
93.6
88.5
84.2
79.3
75.4
.00876
.00770
.00682
.00612
.00557 .00512
.0406 .0249
.0149
8.00
71.3
8.00
67.7 6.00
6.00
K2
4.00
62.8
K2
2.00
2.00
58.6
4.00
.0107
(a)
(b)
Figure C.4
2-D explore plots for sps-b with R vs. K2 and contours of (a) 'ex and (b)
.0808 .0660
96.9
99.2
1-Xi
92.4
2 90.1
.0218
1 2
.00708
87.8
.500
1.50
2.50 R
3.50
4.50
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) deflection mechanisms.
221
Table D.1
Configuration Sub-Class Primary Pivot a 1 lpp a99 1 a95 1 a90 1 b 1 lpp b95 1 b90 1 spp a 1 a99 1 a95 1 a90 1 spp b 1 b95 1 b90 1 psp a 2 psp b 2 pl p a 2 a90 2 a85 2 a80 2 pl p b 2 b85 2 b80 2 slp aIo 2 a99IO 2 a95Io 1 a90I 1 a90IO 2 slp bIO 2 b95IO 1 b90Io 1 ssp aIo 2 a99I 2 a95I 2 a95IO 2 a90IO 2 ssp bIo 2 b94I 2 b90IO 2 b90o 2 lps aIo 3 a99I 3 a95I 3 a90 1 a99IO 3 a95IO 3 a90IO 3 lps bIo 3 b99IO 3 b95IO 3 b90IO 3 b95I 1 b90I 1 sps aIO 3 a99I 1 a90IO 3 sps bIo 3 b98Io 3 b86IO 3 sss a 1 a99 1 a95IO 2 a95I 3 a90IO 2 sss b 1 b99 1 b90I 3 b88IO 2 Configuration Sub-Class Primary Pivot
' ex 99.7 99.0 95.0 90.0 97.6 95.0 90.0 99.7 99.0 95.0 90.0 97.6 95.0 90.0 94.6 86.3 94.6 90.0 85.0 80.0 86.3 85.0 80.0 99.8 99.0 95.0 93.3 97.7 99.4 96.2 90.0 99.8 99.0 96.2 95.0 90.0 99.4 94.1 90.0 91.2 93.1 99.0 95.0 90.0 99.0 95.0 90.8 83.7 99.0 95.0 90.0 95.0 90.0 93.1 99.3 94.6 83.7 98.4 86.3 100.0 99.4 95.0 95.6 91.1 99.5 99.0 96.3 88.4 ' ex
R 0.8274 0.8018 0.7106 0.6185 0.8853 0.8595 0.8226 0.8274 0.8018 0.7104 0.6187 0.8853 0.8595 0.8226 1 1 1 0.4387 0.3170 0.2529 1 0.7919 0.6043 0.3950 0.5057 0.8237 1.5278 0.5437 0.4323 0.6248 0.7267 0.3950 0.7864 0.9182 0.9563 1.4688 0.4323 0.3000 0.8714 0.9323 0.7591 2.5750 2.4052 2.0960 2.1623 2.1623 2.1623 0.8441 1.9336 1.9339 1.9292 0.8561 0.8178 0.7591 0.7328 1.0000 0.8441 0.6528 1.0000 2.6633 1.4903 1.3464 1.6573 1.5518 2.0821 1.8709 1.2420 1.1990 R
K1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 1 0.1906 0.2640 1.6370 15.0000 0.3521 0.2237 0.3924 0.8283 0.1906 0.1000 0.1000 0.5113 0.4050 0.2237 0.1000 0.5343 0.1000 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.0000 0.3497 1.7047 0.4140 0.3970 1.0000 0.6930 0.3330 1.8187 K1
K2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.0029 6.4256 4.3688 2.0313 5.1740 3.6510 2.4156 1.0230 4.9463 4.5230 3.9830 0.0100 0.0100 1.0029 0.1845 1.0022 1.0235 0.2549 1.0000 12.6700 8.3820 2.9708 15.3423 0.6459 9.3816 7.1497 5.4510 3.3349 K2
1 3.50 3.45 3.28 3.10 3.61 3.56 3.49 18.27 18.02 17.10 16.19 18.85 18.60 18.23 13.95 15.06 18.24 25.28 15.44 14.32 16.25 17.27 13.95 17.86 19.18 19.56 24.69 14.32 13.00 18.71 19.32 3.37 6.84 6.52 5.93 6.05 6.05 6.05 3.53 5.62 5.62 5.61 3.55 3.48 17.59 17.33 20.00 18.44 16.53 20.00 36.63 24.90 23.46 26.57 25.52 30.82 28.71 22.42 21.99 1
2 20.00 20.00 3.83 2.75 2.52 2.40 3.83 3.43 3.07 6.76 5.70 4.24 3.17 5.44 6.34 4.98 4.55 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 2
3 23.17 13.88 14.16 14.77 14.62 14.62 14.62 21.85 15.17 15.17 15.18 21.68 22.23 23.17 23.65 20.00 21.85 25.32 20.00 13.75 16.71 17.43 16.03 16.44 14.80 15.34 18.05 18.34 3
0.4537 0.4439 0.4067 0.3653 0.4773 0.4630 0.4421 0.4537 0.4439 0.4066 0.3654 0.4773 0.4630 0.4421 2.0560 2.1500 2.0561 3.4511 4.4878 5.4987 2.1501 2.4473 2.9399 0.9573 1.1290 4.1829 26.3159 1.3688 1.0467 1.4438 2.5248 0.9573 0.6718 0.7015 1.5750 1.3429 1.0467 0.6423 1.7057 0.7255 1.2248 1.4238 1.2692 1.0408 1.4476 1.2475 1.0851 1.2126 1.5344 1.4816 1.4161 0.4683 0.4472 1.2248 0.5711 1.0292 1.2129 0.6380 1.0750 3.4016 3.4388 4.6011 4.6831 1.4903 3.6285 3.0645 3.1967 5.4789
2.901 2.759 2.279 1.832 3.248 3.065 2.812 15.152 14.412 11.895 9.576 16.964 16.010 14.686 82.242 86.000 15.747 13.677 14.904 16.527 16.466 15.046 14.487 18.628 25.595 139.117 1681.489 32.622 21.473 38.116 75.281 18.628 21.440 25.810 60.277 81.851 21.473 10.854 59.736 27.087 7.257 34.840 28.177 19.101 27.717 23.885 20.776 7.895 25.282 24.419 23.263 3.089 2.829 37.901 17.147 41.166 41.247 17.429 43.001 456.482 213.264 253.313 330.679 97.043 344.684 252.583 160.686 264.926
1.000 1.046 1.239 1.492 1.002 1.052 1.129 0.039 0.041 0.049 0.059 0.039 0.041 0.044 0.015 0.016 0.378 0.644 0.698 0.710 0.409 0.500 0.596 0.145 0.219 0.338 0.534 0.233 0.181 0.294 0.331 0.017 0.042 0.048 0.022 0.028 0.019 0.030 0.022 0.052 0.026 0.241 0.253 0.287 0.173 0.211 0.278 0.035 0.222 0.235 0.252 1.053 1.137 0.028 0.063 0.057 0.037 0.031 0.062 0.020 0.079 0.021 0.084 0.031 0.029 0.033 0.076 0.024
1.097 1.098 1.103 1.109 1.094 1.095 1.097 1.027 1.028 1.029 1.031 1.027 1.027 1.027 1.000 1.000 1.088 1.123 1.134 1.141 1.088 1.098 1.110 1.086 1.092 1.107 1.126 1.095 1.088 1.099 1.103 1.036 1.028 1.026 1.026 1.020 1.035 1.038 1.027 1.026 1.122 1.085 1.087 1.091 1.090 1.090 1.090 1.119 1.093 1.093 1.093 1.118 1.120 1.050 1.050 1.050 1.050 1.050 1.050 1.050 1.050 1.050 1.050 1.050 1.050 1.050 1.050 1.050
M 0.4501 0.4373 0.3923 0.3479 0.4788 0.4683 0.4535 2.3511 2.2841 2.0482 1.8177 2.5006 2.4460 2.3687 5.6291 5.4974 1.0777 0.8380 0.8214 0.8272 1.0525 0.9468 0.8593 2.0988 1.6933 2.3414 4.3532 1.5975 1.8373 1.9768 2.1734 6.2080 5.6676 5.6339 5.6304 5.7279 5.7929 6.3534 5.5051 5.4994 3.7184 1.2616 1.3261 1.1869 1.4370 1.4370 1.4370 3.1623 1.7882 1.7881 1.7902 0.4670 0.4516 3.7183 2.1056 2.8145 3.1623 3.9715 2.7487 8.2162 4.2382 5.6882 1.7796 5.7582 5.2859 4.7814 2.3326 5.5108 M
n 0.5004 0.4994 0.4952 0.4898 0.5033 0.5000 0.4949 0.5004 0.4994 0.4951 0.4898 0.5033 0.5000 0.4949 0.5062 0.5164 0.5062 0.5115 0.5176 0.5239 0.5164 0.5178 0.5233 0.5132 0.5097 0.5002 0.5175 0.5090 0.5341 0.4578 0.4794 0.5132 0.5066 0.5062 0.5062 0.5073 0.5341 0.5372 0.5168 0.5165 0.5137 0.4656 0.4699 0.5256 0.4761 0.4761 0.4761 0.5331 0.4250 0.4249 0.4254 0.4995 0.4942 0.5137 0.4963 0.5062 0.5331 0.5572 0.5164 0.5324 0.5168 0.5068 0.4888 0.5076 0.5883 0.5768 0.4923 0.5172 n
222
Table D.2
Configuration Sub-Class a lpp a99 a95 a90 b lpp b95 b90 spp a a99 a95 a90 spp b b95 b90 psp a psp b pl p a a90 a85 a80 pl p b b85 b80 slp aIo a99IO a95Io a90I a90IO slp bIO b95IO b90Io ssp aIo a99I a95I a95IO a90IO ssp bIo b94I b90IO b90o lps aIo a99I a95I a90 a99IO a95IO a90IO lps bIo b99IO b95IO b90IO b95I b90I sps aIO a99I a90IO sps bIo b98Io b86IO sss a a99 a95IO a95I a90IO sss b b99 b90I b88IO Configuration Sub-Class
d Nmax 26.96 26.90 26.57 26.04 39.79 39.68 39.47 26.96 26.90 26.57 26.04 39.79 39.68 39.47 27.13 40.00 27.13 24.03 21.34 19.11 40.00 39.24 36.46 23.23 24.97 26.95 26.28 25.39 30.18 36.92 38.58 23.23 26.85 27.09 27.12 26.43 30.18 23.08 39.74 39.93 26.77 23.10 23.63 24.60 24.39 24.39 24.39 39.60 33.91 33.91 33.95 39.66 39.44 26.77 26.67 27.13 39.60 37.47 40.00 22.82 26.37 26.71 25.93 26.21 32.44 34.51 39.35 39.54 d Nmax
C1
C2
C3
1.71 1.13 1.00 1.00 1.00 1.47 1.00 1.00 5.07 2.54 2.18 2.09 1.33 4.63 6.67 2.30 2.15 9.07 0.68 0.79 1.00 1.00 1.00 1.00 7.33 0.99 0.99 1.00 1.00 1.00 1.74 1.00 1.00 1.40 2.35 1.00 1.000 1.000 1.462 0.340 1.260 1.000 1.000 0.589 1.610 C1
1.00 1.00 1.92 4.82 1.00 1.00 1.09 1.47 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.332 0.745 1.000 0.289 1.000 1.041 0.936 0.381 1.000 C2
1.00 1.00 1.00 0.84 1.00 1.00 1.00 1.00 1.00 1.00 1.00 0.13 0.12 1.00 0.52 1.00 1.00 1.00 1.00 7.099 2.222 2.693 1.000 3.104 4.338 3.502 1.000 2.398 C3
D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class a a99 a95 a90 b b95 b90 a a99 a95 a90 b b95 b90 a b a a90 a85 a80 b b85 b80 2.54 1.00 0.51 1.00 aIo a99IO 1.43 1.00 1.00 1.00 a95Io 1.00 7.04 1.00 1.00 a90I 1.00 119.69 1.00 1.07 1.00 1.00 1.00 1.00 a90IO 1.98 1.00 0.63 1.00 bIO b95IO 1.00 1.28 1.00 1.00 b90Io 1.00 3.14 1.00 1.00 7.52 1.00 0.06 1.00 aIo 11.20 1.00 0.68 1.00 a99I a95I 10.43 1.00 1.01 1.00 a95IO 2.00 1.00 0.22 1.00 a90IO 2.00 1.00 0.84 1.00 6.24 1.00 0.06 1.00 bIo b94I 15.38 1.00 0.05 1.00 b90IO 2.00 1.00 0.17 1.00 b90o 10.35 1.00 1.05 1.00 6.86 1.00 0.01 1.00 aIo 0.32 1.00 1.00 1.00 a99I 0.50 1.00 1.00 1.00 a95I 1.00 0.82 1.00 1.37 a90 0.47 1.00 0.47 1.00 a99IO 0.66 1.00 0.66 1.00 a95IO 1.00 1.00 1.00 1.00 a90IO 6.05 1.00 0.02 1.00 bIo 0.55 1.00 0.56 1.00 b99IO 0.60 1.00 0.61 1.00 b95IO 0.68 1.00 0.68 1.00 b90IO 1.00 0.00 1.00 0.73 b95I 1.00 0.00 1.00 0.98 b90I 1.31 1.00 0.25 1.00 aIO 1.00 0.14 1.00 0.98 a99I 1.00 1.00 1.00 1.00 a90IO 1.16 1.00 0.42 1.00 bIo 6.01 1.00 0.46 1.00 b98Io 1.00 1.00 1.00 1.00 b86IO 1.000 1.832 33.744 1.000 0.775 0.094 a 1.000 0.435 12.492 1.000 1.052 1.139 a99 0.500 1.000 2.000 0.160 1.000 0.102 a95IO 0.039 0.022 1.000 0.997 0.895 1.000 a95I 1.974 1.000 1.979 0.986 1.000 0.066 a90IO 1.000 1.541 19.533 1.000 1.365 0.239 b 1.000 0.995 13.377 1.000 1.215 0.311 b99 0.148 0.055 1.000 0.724 0.996 1.000 b90I 0.500 1.000 2.000 0.120 1.000 0.145 b88IO D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class
Configuration lpp
lpp
spp
spp
psp psp pl p
pl p
slp
slp
ssp
ssp
lps
lps
sps
sps
sss
sss
Configuration
223
Table D.3
Configuration lpp lpp lpp lpp pl p pl p pl p pl p lps lps lps lps lps lps slp slp sss spp spp slp sps spp slp spp sps sps lps sss slp psp ssp ssp ssp sss ssp sss ssp sss
Sub-Class Primary Pivot a90 1 a95 1 a99 1 a 1 a85 2 a80 2 a90 2 a 2 a90 1 a99I 3 a95I 3 a90IO 3 a95IO 3 a99IO 3 a90IO 2 a99IO 2 a95I 3 a90 1 a95 1 aIo 2 a99I 1 a99 1 a95Io 1 a 1 a90IO 3 aIO 3 aIo 3 a99 1 a90I 1 a 2 a95IO 2 a95I 2 a99I 2 a95IO 2 a90IO 2 a90IO 2 aIo 2 a 1 Configuration Sub-Class Primary Pivot lps b90I 1 lpp b90 1 lps b95I 1 lpp b95 1 b 1 lpp pl p b80 2 pl p b85 2 pl p b 2 lps b95IO 3 lps b99IO 3 lps b90IO 3 slp bIO 2 slp b95IO 1 slp b90Io 1 sss b90I 3 spp b90 1 spp b95 1 spp b 1 sps b86IO 3 lps bIo 3 sps bIo 3 sps b98Io 3 sss b99 1 sss b 1 psp b 2 ssp b90o 2 ssp b90IO 2 sss b88IO 2 ssp bIo 2 ssp b94I 2 Configuration Sub-Class Primary Pivot
'ex 90.0 95.0 99.0 99.7 85.0 80.0 90.0 94.6 90.0 99.0 95.0 90.8 95.0 99.0 97.7 99.0 95.6 90.0 95.0 99.8 99.3 99.0 95.0 99.7 94.6 93.1 93.1 99.4 93.3 94.6 95.0 96.2 99.0 95.0 90.0 91.1 99.8 100.0 'ex 90.0 90.0 95.0 95.0 97.6 80.0 85.0 86.3 95.0 99.0 90.0 99.4 96.2 90.0 96.3 90.0 95.0 97.6 86.3 83.7 83.7 98.4 99.0 99.5 86.3 91.2 90.0 88.4 99.4 94.1 'ex
224
Table D.4
Configuration lpp lpp lpp lpp pl p pl p pl p pl p lps lps lps lps lps lps slp slp sss spp spp slp sps spp slp spp sps sps lps sss slp psp ssp ssp ssp sss ssp sss ssp sss
26.04 26.57 26.90 26.96 21.34 19.11 24.03 27.13 24.60 23.10 23.63 24.39 24.39 24.39 25.39 24.97 25.93 26.04 26.57 23.23 26.67 26.90 26.95 26.96 27.13 26.77 26.77 26.37 26.28 27.13 27.12 27.09 26.85 26.71 26.43 26.21 23.23 22.82 Configuration Sub-Class d Nmax lps b90I 39.44 lpp b90 39.47 lps b95I 39.66 lpp b95 39.68 b 39.79 lpp pl p b80 36.46 pl p b85 39.24 pl p b 40.00 lps b95IO 33.91 lps b99IO 33.91 lps b90IO 33.95 slp bIO 30.18 slp b95IO 36.92 slp b90Io 38.58 sss b90I 39.35 spp b90 39.47 spp b95 39.68 spp b 39.79 sps b86IO 40.00 lps bIo 39.60 sps bIo 39.60 sps b98Io 37.47 sss b99 34.51 sss b 32.44 psp b 40.00 ssp b90o 39.93 ssp b90IO 39.74 sss b88IO 39.54 ssp bIo 30.18 ssp b94I 23.08 Configuration Sub-Class d Nmax
Sub-Class a90 a95 a99 a a85 a80 a90 a a90 a99I a95I a90IO a95IO a99IO a90IO a99IO a95I a90 a95 aIo a99I a99 a95Io a a90IO aIO aIo a99 a90I a a95IO a95I a99I a95IO a90IO a90IO aIo a
d Nmax
C1
C2
C3
1.00 0.84 0.68 1.00 0.79 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.00 1.13 1.00 0.340 0.289 1.000
1.00 1.92
1.00 1.74 9.07 1.000 0.745 1.00 4.82 2.09 2.18 2.54 1.462 1.33 1.260 5.07 1.000 C1 1.00 1.00
1.00 1.00 1.00 1.000 2.693 1.00 1.000 3.104 1.00 1.332 7.099 C2 C3 0.12 0.13
D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class a90 a95 a99 a a85 a80 a90 a 1.00 0.82 1.00 1.37 a90 0.32 1.00 1.00 1.00 a99I 0.50 1.00 1.00 1.00 a95I 1.00 1.00 1.00 1.00 a90IO 0.66 1.00 0.66 1.00 a95IO 0.47 1.00 0.47 1.00 a99IO a90IO 1.00 1.00 1.00 1.00 a99IO 1.43 1.00 1.00 1.00 0.039 0.022 1.000 0.997 0.895 1.000 a95I a90 a95 aIo 2.54 1.00 0.51 1.00 1.00 0.14 1.00 0.98 a99I a99 a95Io 1.00 7.04 1.00 1.00 a 1.00 1.00 1.00 1.00 a90IO 1.31 1.00 0.25 1.00 aIO 6.86 1.00 0.01 1.00 aIo 1.000 0.435 12.492 1.000 1.052 1.139 a99 1.00 119.69 1.00 1.07 a90I a 2.00 1.00 0.22 1.00 a95IO a95I 10.43 1.00 1.01 1.00 a99I 11.20 1.00 0.68 1.00 0.500 1.000 2.000 0.160 1.000 0.102 a95IO a90IO 2.00 1.00 0.84 1.00 1.974 1.000 1.979 0.986 1.000 0.066 a90IO aIo 7.52 1.00 0.06 1.00 1.000 1.832 33.744 1.000 0.775 0.094 a D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class 1.00 1.00 0.00 0.00 1.00 1.00 0.98 b90I b90 0.73 b95I b95 b b80 b85 b 1.00 b95IO 1.00 b99IO 1.00 b90IO bIO b95IO b90Io 1.000 b90I b90 b95 b 1.00 b86IO 1.00 bIo 1.00 bIo 1.00 b98Io 0.311 b99 0.239 b b b90o b90IO 0.145 b88IO bIo b94I D 3min Sub-Class
Configuration lpp lpp lpp lpp pl p pl p pl p pl p lps lps lps lps lps lps slp slp sss spp spp slp sps spp slp spp sps sps lps sss slp psp ssp ssp ssp sss ssp sss ssp sss Configuration lps lpp lps lpp lpp pl p pl p pl p lps lps lps slp slp slp sss spp spp spp sps lps sps sps sss sss psp ssp ssp sss ssp ssp Configuration
1.00 0.60 1.00 0.55 1.00 0.68 1.98 1.00 1.00 1.000 0.148
1.00 0.61 1.00 0.56 1.00 0.68 0.63 1.00 1.00 1.000 0.724
1.00 1.00 1.00 7.33 1.00 6.05 1.40 1.00 1.16 2.35 1.00 6.01 1.000 0.936 3.502 1.000 1.000 1.041 4.338 1.000 2.15 2.30 1.610 4.63 6.67 C1
0.995 1.541
1.00 10.35 1.00 1.05 1.00 1.00 2.00 1.00 0.17 1.00 1.000 2.398 0.500 1.000 2.000 0.120 1.000 1.00 6.24 1.00 0.06 1.00 15.38 1.00 0.05 1.00 1.00 C2 C 3 D 1equal D 2equal D 3equal D 1min D 2min
225
Table D.5
K1 R Sub-Class Primary Pivot 'ex a90 1 90.0 0.6185 0 a95 1 95.0 0.7106 0 a99 1 99.0 0.8018 0 99.7 0.8274 0 a 1 a80 2 80.0 0.2529 1 a85 2 85.0 0.3170 1 90.0 0.4387 1 a90 2 a90I 1 93.3 1.5278 15.0000 a 2 94.6 1 1 a95Io 1 95.0 0.8237 1.6370 a90 1 90.0 2.0960 0 a90IO 3 90.8 2.1623 0 a95I 3 95.0 2.4052 0 a99I 3 99.0 2.5750 0 a90IO 2 97.7 0.5437 0.3521 a99IO 2 99.0 0.5057 0.2640 a95IO 3 95.0 2.1623 0 a99IO 3 99.0 2.1623 0 aIo 2 99.8 0.3950 0.1906 a95I 3 95.6 1.6573 0.4140 a99 1 99.4 1.4903 0.3497 a99I 1 99.3 0.7328 0 a90 1 90.0 0.6187 0 a90IO 3 94.6 1.0000 0 a95 1 95.0 0.7104 0 a95I 2 96.2 0.9182 0.1000 a99I 2 99.0 0.7864 0.1000 a99 1 99.0 0.8018 0 a 1 99.7 0.8274 0 a90IO 2 91.1 1.5518 0.3970 aIO 3 93.1 0.7591 0 a90IO 2 90.0 1.4688 0.4050 aIo 3 93.1 0.7591 0 a95IO 2 95.0 0.9563 0.5113 a95IO 2 95.0 1.3464 1.7047 a 1 100.0 2.6633 1.0000 aIo 2 99.8 0.3950 0.1906 94.6 1 1 a 2 K1 R Configuration Sub-Class Primary Pivot 'ex lps b90I 1 90.0 0.8178 0 lpp b90 1 90.0 0.8226 0 lps b95I 1 95.0 0.8561 0 lpp b95 1 95.0 0.8595 0 b 1 97.6 0.8853 0 lpp pl p b80 2 80.0 0.6043 1 pl p b85 2 85.0 0.7919 1 pl p b 2 86.3 1 1 slp b90Io 1 90.0 0.7267 0.8283 slp b95IO 1 96.2 0.6248 0.3924 lps b90IO 3 90.0 1.9292 0 lps b95IO 3 95.0 1.9339 0 lps b99IO 3 99.0 1.9336 0 slp bIO 2 99.4 0.4323 0.2237 sss b90I 3 96.3 1.2420 0.3330 sps b86IO 3 86.3 1.0000 0 ssp b90o 2 91.2 0.9323 0.1000 spp b90 1 90.0 0.8226 0 spp b95 1 95.0 0.8595 0 spp b 1 97.6 0.8853 0 sps bIo 3 83.7 0.8441 0 lps bIo 3 83.7 0.8441 0 sss b99 1 99.0 1.8709 0.6930 sps b98Io 3 98.4 0.6528 0 ssp b94I 2 94.1 0.3000 0.1000 sss b 1 99.5 2.0821 1.0000 sss b88IO 2 88.4 1.1990 1.8187 ssp b90IO 2 90.0 0.8714 0.5343 ssp bIo 2 99.4 0.4323 0.2237 psp 86.3 1 1 b 2 K1 R Configuration Sub-Class Primary Pivot 'ex Configuration lpp lpp lpp lpp pl p pl p pl p slp pl p slp lps lps lps lps slp slp lps lps slp sss sss sps spp sps spp ssp ssp spp spp sss sps ssp lps ssp sss sss ssp psp
K2 0 0 0 0 0 0 0 0 0 0 2.0313 2.4156 4.3688 6.4256 0 0 3.6510 5.1740 0 15.3423 8.3820 0.1845 0 1.0022 0 0 0 0 0 0.6459 1.0029 0 1.0029 0 2.9708 12.6700 0 0 K2 0.0100 0 0.0100 0 0 0 0 0 0 0 3.9830 4.5230 4.9463 0 5.4510 1.0000 0 0 0 0 1.0235 1.0230 7.1497 0.2549 0 9.3816 3.3349 0 0 0 K2
1 3.10 3.28 3.45 3.50 25.28 18.24 5.93 6.05 6.52 6.84 15.44 15.06 6.05 6.05 13.95 26.57 24.90 17.33 16.19 20.00 17.10 19.18 17.86 18.02 18.27 25.52 17.59 24.69 3.37 19.56 23.46 36.63 13.95 1 3.48 3.49 3.55 3.56 3.61 17.27 16.25 5.61 5.62 5.62 14.32 22.42 20.00 19.32 18.23 18.60 18.85 18.44 3.53 28.71 16.53 13.00 30.82 21.99 18.71 14.32 1
2 2.40 2.52 2.75 3.17 3.83 4.24 5.44 5.70 6.76 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 2 3.07 3.43 3.83 4.55 4.98 6.34 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 20.00 2
3 14.77 14.62 14.16 13.88 14.62 14.62 16.03 16.71 23.65 20.00 16.44 23.17 23.17 17.43 13.75 3 22.23 21.68 15.18 15.17 15.17 18.05 20.00 21.85 21.85 15.34 25.32 14.80 18.34 3
0.3653 0.4067 0.4439 0.4537 5.4987 4.4878 3.4511 26.3159 2.0561 4.1829 1.0408 1.0851 1.2692 1.4238 1.3688 1.1290 1.2475 1.4476 0.9573 4.6831 3.4388 0.5711 0.3654 1.0292 0.4066 0.7015 0.6718 0.4439 0.4537 1.4903 1.2248 1.3429 1.2248 1.5750 4.6011 3.4016 0.9573 2.0560 0.4472 0.4421 0.4683 0.4630 0.4773 2.9399 2.4473 2.1501 2.5248 1.4438 1.4161 1.4816 1.5344 1.0467 3.1967 1.0750 0.7255 0.4421 0.4630 0.4773 1.2129 1.2126 3.0645 0.6380 0.6423 3.6285 5.4789 1.7057 1.0467 2.1500
1.832 2.279 2.759 2.901 16.527 14.904 13.677 1681.489 15.747 139.117 19.101 20.776 28.177 34.840 32.622 25.595 23.885 27.717 18.628 330.679 213.264 17.147 9.576 41.166 11.895 25.810 21.440 14.412 15.152 97.043 37.901 81.851 7.257 60.277 253.313 456.482 18.628 82.242 2.829 2.812 3.089 3.065 3.248 14.487 15.046 16.466 75.281 38.116 23.263 24.419 25.282 21.473 160.686 43.001 27.087 14.686 16.010 16.964 41.247 7.895 252.583 17.429 10.854 344.684 264.926 59.736 21.473 86.000
1.492 1.239 1.046 1.000 0.710 0.698 0.644 0.534 0.378 0.338 0.287 0.278 0.253 0.241 0.233 0.219 0.211 0.173 0.145 0.084 0.079 0.063 0.059 0.057 0.049 0.048 0.042 0.041 0.039 0.031 0.028 0.028 0.026 0.022 0.021 0.020 0.017 0.015 1.137 1.129 1.053 1.052 1.002 0.596 0.500 0.409 0.331 0.294 0.252 0.235 0.222 0.181 0.076 0.062 0.052 0.044 0.041 0.039 0.037 0.035 0.033 0.031 0.030 0.029 0.024 0.022 0.019 0.016
1.109 1.103 1.098 1.097 1.141 1.134 1.123 1.126 1.088 1.107 1.091 1.090 1.087 1.085 1.095 1.092 1.090 1.090 1.086 1.050 1.050 1.050 1.031 1.050 1.029 1.026 1.028 1.028 1.027 1.050 1.050 1.020 1.122 1.026 1.050 1.050 1.036 1.000 1.120 1.097 1.118 1.095 1.094 1.110 1.098 1.088 1.103 1.099 1.093 1.093 1.093 1.088 1.050 1.050 1.026 1.027 1.027 1.027 1.050 1.119 1.050 1.050 1.038 1.050 1.050 1.027 1.035 1.000
M 0.3479 0.3923 0.4373 0.4501 0.8272 0.8214 0.8380 4.3532 1.0777 2.3414 1.1869 1.4370 1.3261 1.2616 1.5975 1.6933 1.4370 1.4370 2.0988 1.7796 4.2382 2.1056 1.8177 2.8145 2.0482 5.6339 5.6676 2.2841 2.3511 5.7582 3.7183 5.7279 3.7184 5.6304 5.6882 8.2162 6.2080 5.6291 M 0.4516 0.4535 0.4670 0.4683 0.4788 0.8593 0.9468 1.0525 2.1734 1.9768 1.7902 1.7881 1.7882 1.8373 2.3326 2.7487 5.4994 2.3687 2.4460 2.5006 3.1623 3.1623 4.7814 3.9715 6.3534 5.2859 5.5108 5.5051 5.7929 5.4974 M
n 0.4898 0.4952 0.4994 0.5004 0.5239 0.5176 0.5115 0.5175 0.5062 0.5002 0.5256 0.4761 0.4699 0.4656 0.5090 0.5097 0.4761 0.4761 0.5132 0.4888 0.5168 0.4963 0.4898 0.5062 0.4951 0.5062 0.5066 0.4994 0.5004 0.5076 0.5137 0.5073 0.5137 0.5062 0.5068 0.5324 0.5132 0.5062 n 0.4942 0.4949 0.4995 0.5000 0.5033 0.5233 0.5178 0.5164 0.4794 0.4578 0.4254 0.4249 0.4250 0.5341 0.4923 0.5164 0.5165 0.4949 0.5000 0.5033 0.5331 0.5331 0.5768 0.5572 0.5372 0.5883 0.5172 0.5168 0.5341 0.5164 n
226
Table D.6
Configuration lpp lpp lpp lpp pl p pl p pl p slp pl p slp lps lps lps lps slp slp lps lps slp sss sss sps spp sps spp ssp ssp spp spp sss sps ssp lps ssp sss sss ssp psp
Sub-Class a90 a95 a99 a a80 a85 a90 a90I a a95Io a90 a90IO a95I a99I a90IO a99IO a95IO a99IO aIo a95I a99 a99I a90 a90IO a95 a95I a99I a99 a a90IO aIO a90IO aIo a95IO a95IO a aIo a Configuration Sub-Class lps b90I lpp b90 lps b95I lpp b95 b lpp pl p b80 pl p b85 pl p b slp b90Io slp b95IO lps b90IO lps b95IO lps b99IO slp bIO sss b90I sps b86IO ssp b90o spp b90 spp b95 spp b sps bIo lps bIo sss b99 sps b98Io ssp b94I sss b sss b88IO ssp b90IO ssp bIo psp b Configuration Sub-Class
d Nmax 26.04 26.57 26.90 26.96 19.11 21.34 24.03 26.28 27.13 26.95 24.60 24.39 23.63 23.10 25.39 24.97 24.39 24.39 23.23 25.93 26.37 26.67 26.04 27.13 26.57 27.09 26.85 26.90 26.96 26.21 26.77 26.43 26.77 27.12 26.71 22.82 23.23 27.13 d Nmax 39.44 39.47 39.66 39.68 39.79 36.46 39.24 40.00 38.58 36.92 33.95 33.91 33.91 30.18 39.35 40.00 39.93 39.47 39.68 39.79 39.60 39.60 34.51 37.47 23.08 32.44 39.54 39.74 30.18 40.00 d Nmax
C1
C2
C3
1.00 1.00 1.00 1.00 0.79 0.68 1.00 1.13 1.00 1.00 1.71 0.340 1.000 1.00 1.00 2.18 2.54
0.84 1.00 1.00 1.00 1.00 1.00 1.000 2.222 0.52 1.00 -
1.260 1.74 1.33 9.07 2.09 1.462 1.000 5.07 C1 1.00 1.00
D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class a90 a95 a99 a a80 a85 a90 a90I 1.00 119.69 1.00 1.07 a 1.00 7.04 1.00 1.00 a95Io 1.00 0.82 1.00 1.37 a90 1.00 1.00 1.00 1.00 a90IO 0.50 1.00 1.00 1.00 a95I 0.32 1.00 1.00 1.00 a99I a90IO 1.00 1.00 1.00 1.00 a99IO 1.43 1.00 1.00 1.00 0.66 1.00 0.66 1.00 a95IO 0.47 1.00 0.47 1.00 a99IO aIo 2.54 1.00 0.51 1.00 0.039 0.022 1.000 0.997 0.895 1.000 a95I 1.000 0.435 12.492 1.000 1.052 1.139 a99 1.00 0.14 1.00 0.98 a99I a90 1.00 1.00 1.00 1.00 a90IO a95 a95I 10.43 1.00 1.01 1.00 11.20 1.00 0.68 1.00 a99I a99 a 1.974 1.000 1.979 0.986 1.000 0.066 a90IO 1.31 1.00 0.25 1.00 aIO a90IO 2.00 1.00 0.84 1.00 6.86 1.00 0.01 1.00 aIo a95IO 2.00 1.00 0.22 1.00 0.500 1.000 2.000 0.160 1.000 0.102 a95IO 1.000 1.832 33.744 1.000 0.775 0.094 a 7.52 1.00 0.06 1.00 aIo a D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class 1.00 0.00 1.00 0.98 b90I b90 1.00 0.00 1.00 0.73 b95I b95 b b80 b85 b 1.00 3.14 1.00 1.00 b90Io b95IO 1.00 1.28 1.00 1.00 0.68 1.00 0.68 1.00 b90IO 0.60 1.00 0.61 1.00 b95IO 0.55 1.00 0.56 1.00 b99IO bIO 1.98 1.00 0.63 1.00 0.148 0.055 1.000 0.724 0.996 1.000 b90I 1.00 1.00 1.00 1.00 b86IO b90o 10.35 1.00 1.05 1.00 b90 b95 b 1.16 1.00 0.42 1.00 bIo 6.05 1.00 0.02 1.00 bIo 1.000 0.995 13.377 1.000 1.215 0.311 b99 6.01 1.00 0.46 1.00 b98Io 15.38 1.00 0.05 1.00 b94I 1.000 1.541 19.533 1.000 1.365 0.239 b 0.500 1.000 2.000 0.120 1.000 0.145 b88IO 2.00 1.00 0.17 1.00 b90IO 6.24 1.00 0.06 1.00 bIo b D 1equal D 2equal D 3equal D 1min D 2min D 3min Sub-Class
Configuration lpp lpp lpp lpp pl p pl p pl p slp pl p slp lps lps lps lps slp slp lps lps slp sss sss sps spp sps spp ssp ssp spp spp sss sps ssp lps ssp sss sss ssp psp Configuration lps lpp lps lpp lpp pl p pl p pl p slp slp lps lps lps slp sss sps ssp spp spp spp sps lps sss sps ssp sss sss ssp ssp psp Configuration
227
228
APPENDIX E
The stress factor, a non-dimensionalized parameter where K = ------------------- . This parameter is used in the stress feasibil ity equation.
229
Length parameter used to calculate the actual length from the PRBM length where ltot = r tot . This parameter is mechanism dependent.
Ratio of the small-length flexural pivot length over the associated PRBM length
(Stiffness Intensity Factor) This parameter summarizes the stiffness at a given primary pivot stress level and mechanism size. This allows for comparisons in maximum stiffnesses between mechanisms.
(Normalized Stiffness Intensity Factor) This parameter divides the stiffness intensity factor by the stiffness intensity factor of the lpp mechanims for the given sub-class. This normalizes the values with 1 being the highest.
230
Defines the ratio between the total PRBM length r tot and a link length
The original percent constant-force where the maximum force is divided by the minimum force.
'
Percent constant-force using the minimum and maximum known forces (no extrapolation).
' ex
The extrpolated percent constant-force calculated by dividing the extraoplated force at zero displacement by the maximum force at full displacement.
CFEC
CFM
Constant-force mechanism
Classification Varaibles
Those variables that dependent on the mechanism and are looked up in tables. 231
Configuration
Coupled Variables
Those variables that are shared between tow or more of the principal design equations.
Dequal
Ratio of a flexible segment width and the primary pivot width when thicknesses are equal.
Dmin
Ratio of a flexible segment width and the primary pivot width when thicknesses are maximized.
dNmax
The maximum normal displacement percentage which allows for quick calculations of the maximum normal deflection.
Displacement Equaiton
Main equation used to determine the percent displacement when doing design of CFMs.
Main equation used to calculate the force of a CFM. Knowledge of the PRBM is not required to use this equation.
232
Force Feasibility
EI 1 F = ----------- Equation used to quickly calculate the force 2 r tot knowing only the cross section of the first flexible segment, the material properties, the configuraiton, and the total PRBM length.
Those mechanisms that are guaranteed to produce the needed force without violating the stress requirments. However, no guarantee is given that the flexible segment thicknesses will satisfy the design constraints.
Those mechanisms that are guaranteed to procude the desired displacement in any part of the design space. However, no guarantee is given that the need force can be achieved.
In-plane Orientation
In this orientation, all motion takes place in the plane of fabrication. This requires that all link and flexible segment widths be the same.
Isolated Varaibles
Multiplier in the power curve fit of . Power in the power curve fit of .
Out-of-plane Orientation
All motion is normal to the plane of fabrication. This requires that all link and flexible segment thicknesses be the same.
Principal Equations
The three primary equations required to design CFMs. These equations are combinations of the stress feasibility equation, force feasibility equation, and other equations developed.
Primary Pivot
Main equation used to determine the stress within a CFM. Knowledge of the PRBM is not required to use this equation. A ------ Equation used to determine the stress feasibility SF of a given design.
Stress Feasibility
Those mechanisms that are capable of producing the desired displacement in part of the design space, but not capable in
234
other parts of the design space. However, no guarantee is given that the need force can be achieved.
Sub-class
Distinguishes between different sets of constant-force parameters within a given class and configuration thus specifying a specific mechanism.
Those variables values that have no information given by the design requirements.
235
236