COUPLED BELT-PULLEY MECHANICS IN SERPENTINE
BELT DRIVES
DISSERTATION
Presented in Partial Fulfillment of the Requirements of the Degree Doctor of Philosophy in the Graduate School of The Ohio State University
By
Lingyuan Kong, B.S., M.S. *****
The Ohio State University 2003
Dissertation Committee:
Prof. Robert G. Parker, Advisor Approved by
Prof. Stephen E. Bechtel
Prof. Chia-Hsiang Menq ______Advisor Prof. Rajendra Singh Department of Mechanical Engineering
ABSTRACT
Belt vibration and slip are primary concerns in the design of serpentine belt drives.
Belt-pulley coupling is essential for the analysis. This work investigates issues to advance the understanding of belt-pulley mechanics.
Closed-form eigensolution approximations for an axially moving beam with small bending stiffness are given. This model is the first order approximation for the transverse vibration of each span in a serpentine belt drive. Perturbation techniques for algebraic equations and the phase closure principle are used. The eigensolutions are interpreted in terms of propagating waves.
For a complete serpentine belt drive, a hybrid continuous-discrete model is built.
Incorporation of belt bending stiffness introduces linear belt-pulley coupling. This model can explain the transverse span vibrations caused by crankshaft pulley fluctuations at low engine idle speeds where other coupling mechanisms do not. For the steady state analysis, a novel transformation of the governing equations to a standard ODE form for general-purpose BVP solvers leads to numerically exact steady solutions. A closed-form singular perturbation solution is developed for the small bending stiffness case. A coupling indicator based on the steady state is defined to quantify the undesirable belt- pulley coupling. A spatial discretization is developed to find the free vibration ii eigensolutions. In contrast to prior formulations, this discretization is numerically robust and free of missing/false natural frequency concerns. New dynamic properties induced by bending stiffness are characterized. Dynamic response calculations using the discretized model follow naturally. The effects of major design variables are investigated. This provides knowledge to help optimize structural design, especially to reduce large belt transverse vibration.
Finally, to better predict the belt-pulley contact interactions applicable to serpentine belt drives an improved model is established for the steady state mechanics. Bending stiffness is considered while other factors in the literature such as belt-pulley friction and belt inertia are retained. An iterative solution based on general-purpose BVP solvers is presented to determine the belt deflections and the distributions of speed, tension, and friction along the belt as well as the belt-pulley contact points and adhesion/slip zones on the pulleys. Key design criteria like maximum transmissible moment and power efficiency are examined.
iii
Dedicated to my wife, Hong Chi
iv
ACKNOWLEDGMENTS
I wish to express my sincere thanks to Dr. Robert G. Parker for his assistance and guidance throughout the research project upon which this dissertation is based.
Especially, I am grateful for his help in continuously correcting and modifying my writing, from which I benefited greatly. I also thank Dr. Stephen E. Bechtel, Dr. Chia-
Hsiang Menq, and Dr. Rajendra Singh, who served as members of my dissertation committee. I would like to acknowledge the generous financial support given to this project by Mark IV/Dayco Corporation. Furthermore, I thank my colleagues for their suggestions regarding my dissertation. Finally I thank my wife who constantly supported me throughout this endeavor.
v
VITA
March, 1972…..……………….Born Guizhou, China
1993……………………………B.S., Shanghai Jiaotong University
1999……………………………M.S., Tsinghua University
1999 – present …………………Graduate Research Assistant, The Ohio State University
PUBLICATIONS
1. L. Kong and R. G. Parker, 2003, “Equilibrium and Belt-Pulley Vibration Coupling in Serpentine Belt Drives,” ASME Journal of Applied Mechanics, 70(5), pp.739-750.
FIELDS OF STUDY
Major Field: Mechanical Engineering
vi
TABLE OF CONTENTS
Abstract ...... ii Dedication ...... iv Acknowledgments...... v Vita...... vi List of Figures ...... ix List of Tables...... xiii Chapters
CHAPTER 1 INTRODUCTION ...... 1 1.1 Motivation and Objectives ...... 1 1.2 Literature Review...... 6 1.2.1 Vibration of axially moving materials ...... 7 1.2.2 Serpentine belt drives...... 8 1.2.3 Belt-pulley steady state contact mechanics...... 12 1.3 Scope of Investigation...... 13 CHAPTER 2 APPROXIMATE EIGENSOLUTIONS OF AXIALLY MOVING BEAMS WITH SMALL BENDING STIFFNESS...... 17 2.1 Introduction ...... 18 2.2 Model Equations ...... 22 2.3 Application of the Phase Closure Principle...... 26 2.4 Other Boundary Conditions ...... 32 CHAPTER 3 MODELING AND STEADY STATE ANALYSIS OF SERPENTINE BELT DRIVES WITH BENDING STIFFNESS...... 36 3.1 Introduction ...... 37 3.2 System Model...... 40
vii
3.3 Numerical Solution ...... 49 3.4 Numerical Results and Discussion...... 52 3.5 Approximate Closed-Form Solution ...... 64 CHAPTER 4 DYNAMIC ANALYSIS OF SERPENTINE BELT DRIVES WITH BENDING STIFFNESS...... 79 4.1 Introduction ...... 80 4.2 Linearization of Equations of Motion ...... 83 4.3 Extended Operator Formulation...... 88 4.4 Galerkin Discretization ...... 94 4.5 Results and Discussion...... 96 CHAPTER 5 STEADY STATE MECHANICS OF BELT-PULLEY SYSTEMS ...... 113 5.1 Introduction ...... 114 5.2 Nonlinear Equations of a Moving Curved Beam...... 116 5.3 BVP Solver Based Method for Problem with Unknown Boundaries ...... 119 5.4 Steady State Analysis of Belt Pulley Drives: an Iteration Method ...... 122 5.4.1 Regular Moment Transmission Problem ...... 125 5.4.2 Maximum Transmissible Moment Problem...... 131 5.5 Results and Discussion...... 134 5.5.1 Example of regular moment transmission problem ...... 136 5.5.2 Example of maximum transmissible moment problem ...... 140 CHAPTER 6 SUMMARY AND FUTURE WORK...... 145 6.1 Summary ...... 145 6.2 Future Work ...... 150 REFERENCE……………… …………………………………………………………..158
viii
LIST OF FIGURES
Figure Page
Figure 1.1 Installed serpentine belt drive...... 3 Figure 1.2 Sketch of a seven-pulley serpentine belt drive...... 4 Figure 2.1. Waves in a finite moving beam with small bending stiffness ...... 26 =− =2 −222 − Figure 2.2. Phase angle of D RrrrrA ( 3 231)/( ) from equation (2.21),
=−εω 2 =−−+εω2 22 ω D1 2/1vv, Dvv2 {[/(1)][/(1)]}, =− −2 +εω2222 − − ω + Dvvv3 (1 ) {[ /(1 )] [ /(1 v )] } ...... 28 Figure 2.3. Comparison of fundamental natural frequency for a simply supported beam.
−−−, perturbation; , exact...... 30 Figure 3.1. A prototypical three-pulley serpentine belt drive system...... 38 Figure 3.2. Detail of tensioner region and pulley 2, defining alignment angles...... 43 Figure 3.3. Steady state deflections of spans 1 and 3 for varying belt bending stiffness. = = γ = = = = s 0, ks 4 , 400 , P1 P2 P3 1...... 54 Figure 3.4. Steady state deflections of spans 1 and 3 for varying span tensions. ε = 0.05 , = = γ = s 0, ks 4 , 400 ...... 55 ε = = Figure 3.5. Steady state deflections of spans 1 and 3 for varying speed. 0.05 , ks 4 , γ = = = = 400 , P1 P2 P3 1...... 58 Figure 3.6. Steady state deflections of spans 1 and 3 for varying tensioner spring ε = = γ = = = = stiffness. 0.05 , s 0 , 400 , P1 P2 P3 1...... 59 Figure 3.7. Steady state deflections of spans 1 and 3 for varying longitudinal belt ε ======stiffness. 0.05 , s 0 , ks 4 , P1 P2 P3 1...... 60
ix
Figure 3.8. System coupling indicator Γ for varying belt bending stiffness. s = 0 , = γ = = = = k s 4, 400 , P1 P2 P3 1...... 61 Figure 3.9. System coupling indicator Γ for varying span tensions. ε = 0.015 , s = 0 , = γ = k s 4, 400 ...... 62 Γ ε = = Figure 3.10. System coupling indicator for varying speed. 0.015 , k s 4, γ = = = = η = 400 , P1 P2 P3 1, 0.78 corresponds to tensioner orientation β = o β = o η = 1 135.79 , 2 178.74 in Table 3.2; while 0 corresponds to tensioner β = o β = o orientation: 1 68.53 , 2 111.47 ...... 72 = γ = Figure 3.11. Variation of tension with different tensioner orientation. k s 4, 400 , = = = β = o β = o η = β = o P1 P2 P3 1, a) 1 135.79 , 2 178.74 , 0.78 , b) 1 68.53 , β = o η = 2 111.47 , 0 ...... 73 Figure 3.12. Steady state deflections (dimensionless) of the first and third spans. ε = = = γ = = = = 0.01, s 0.6 , ks 4 , 400 , P1 P2 P3 1...... 74 Figure 3.13. Steady state deflections (dimensionless) of the first and third spans.ε = 0.01, = = γ = = = = s 0.9 , ks 4 , 400 , P1 P2 0.9395 , P3 1.5395...... 75 Figure 4.1. A prototypical three-pulley serpentine belt drive...... 83 Figure 4.2. Rotationally dominant mode (ε = 0 ) for increasing belt bending stiffness. The dimensionless natural frequency for ε = 0 is ω = 4.1205 . a) ε = 0.01, b) ε = 0.04 , ε = ε = = = γ = = = = β = o c) 0.07 , d) 0.1. s 0, ks 4 , 400 , P1 P2 P3 1, 1 135.79 , β = o 2 178.74 ...... 98 Figure 4.3. Span 2 transversely dominant mode (ε = 0 ) for increasing belt bending stiffness. The dimensionless natural frequency for ε = 0 is ω = 3.0951. a) ε = 0.01, ε = ε = ε = = = γ = = = = b) 0.04 , c) 0.07 , d) 0.1. s 0, ks 4 , 400 , P1 P2 P3 1, β = o β = o 1 135.79 , 2 178.74 ...... 99 Figure 4.4. Span 3 transversely dominant mode (ε = 0 ) for increasing belt bending stiffness. The dimensionless natural frequency for ε = 0 is ω =1.9968. a) ε = 0.01,
x
ε = ε = ε = = = γ = = = = b) 0.04 , c) 0.07 , d) 0.1. s 0, ks 4 , 400 , P1 P2 P3 1, β = o β = o 1 135.79 , 2 178.74 ...... 100 = = Figure 4.5. Natural frequency spectrum for varying belt bending stiffness. s 0, ks 4 , γ = = = = β = o β = o 400 , P1 P2 P3 1, 1 135.79 , 2 178.74 ...... 102 Figure 4.6. Natural frequency spectrum for varying belt bending stiffness. −−−, fix
= = steady state; , fix bending stiffness value in (4.39)-(4.45). s 0 , ks 4 ,
γ = = = = β = o β = o 400 , P1 P2 P3 1, 1 135.79 , 2 178.74 ...... 105
Figure 4.7. Natural frequency spectrum for varying belt speed. , ε = 0.1; −−−,
ε = = γ = = = = β = o β = o 0.01. ks 4 , 400 , P1 P2 P3 1, 1 135.79 , 2 178.74 ...... 106
η = β = o Figure 4.8. Natural frequency spectrum for varying belt speed. , 0 ( 1 68.53 ,
β = o −−− η = β = o β = o ε = = 2 111.47 ); , 0.78 ( 1 135.79 , 2 178.74 ). 0.04 , ks 4 , γ = = = = 400 , P1 P2 P3 1...... 107
Figure 4.9. Natural frequency spectrum for varying tensioner effectiveness η . ,
ε = −−− ε = = = γ = = = = 0.1; , 0.01. s 0 , ks 4 , 400 , P1 P2 P3 1...... 108 Figure 4.10. Fifth vibration mode for varying tensioner effectiveness η . a) η = 0
β = o β = o η = β = o β = o ( 1 68.53 , 2 111.47 ), b) 0.5 ( 1 95.53 , 2 138.47 ), c) η = β = o β = o ε = = = γ = 0.78 ( 1 135.79 , 2 178.74 ). 0.1, s 0, ks 4 , 400 , = = = P1 P2 P3 1...... 112 Figure 5.1 Free body diagram of a moving curved beam...... 117 Figure 5.2. Single span boundary value problem with unknown boundaries...... 120 Figure 5.3. Two-pulley belt drive with inclusion of belt bending stiffness...... 127 Figure 5.4. Steady solutions for the system properties specified in Table 5.1. a) EI = 0.0015, b) EI = 0.015, c) EI = 0.05 N ⋅m2 ...... 135
xi
Figure 5.5. Variations of tension and speed in the tight and slack spans for the belt-pulley drive in Table 5.1...... 138 Figure 5.6. Deflections of the free spans for two different belt-pulley models. (a) and (b) correspond to the current model (symbols denote span endpoints); (c) and (d) correspond to the fixed boundary model in chapter 3 [11]. The system is specified in Table 5.1...... 139 Figure 5.7. Steady solutions for the system properties specified in Table 5.3. Full slip occurs on the driver pulley. a) EI = 0.0015, b) EI = 0.015, c) EI = 0.05 N ⋅m2 . .. 142 = Figure 5.8. Comparison of maximum transmitted moment M maxMTR 2 _ max / ini 2 , power η efficiency , Tt _ midspan/T ini , and Ts _ midspan/T ini between the string and beam models for the belt drive in Table 5.3...... 144 Figure 6.1. Model considering the radial stiffness of the belt...... 150 Figure 6.2. Steady motion of an extensible belt considering bending stiffness...... 154 Figure 6.3. a) Self-sustained system, b) friction force characteristics, c) limit cycle. .... 155 Figure 6.4. Two-dimensional friction between belt and grooved pulley...... 157
xii
LIST OF TABLES
Table Page
Table 2.1. Comparison of approximate roots of (2.4) from (2.15) with numerically exact roots...... 24 Table 2.2. Reflection coefficients for different end supports...... 33 Table 3.1. List of symbols defining serpentine belt drives...... 42 Table 3.2. Physical properties of the prototypical system...... 53 Table 3.3 Comparison of approximate analytical and numerical solutions. Case 1 parameters are those used in Figure 3.12, and Case 2 parameters are those used in Figure 3.13...... 70 Table 4.1. Physical properties of the example system, from which nominal dimensionless parameters are calculated...... 97 Table 5.1. Physical properties of the belt drive with two identical pulleys...... 136 Table 5.2. Numerical results for the belt drive specified in Table 5.1...... 137 Table 5.3. Physical properties of the belt drive with two different pulleys...... 137
xiii
CHAPTER 1
INTRODUCTION
1.1 Motivation and Objectives
Serpentine belt drives with a long, flat, multi-ribbed belt are widely used in the automobile industry for passenger vehicles and heavy duty trucks (Figure 1.1). With a single belt, the engine power is delivered from the crankshaft to all of the individual accessories (air conditioner, alternator, power steering, water pump, etc.). To maintain proper belt tension as accessory loads and engine speed vary, a spring-loaded tensioner is introduced (Figure 1.2). There are many advantages of the serpentine belt drives over the
V-belt drives including compactness, longer life, automatic tension loss compensation, ease of assembly, and so on [1-3].
Crankshaft torque pulsations, resulting from cylinder combustion in the engine, and dynamic accessory torques excite rotational vibration of the pulleys. These pulley vibrations can then be further transmitted to the belt spans. This belt-pulley coupling has strong influence on the entire system’s dynamic behavior. Under certain circumstances, 1 the coupling is so strong that some belt spans experience large transverse vibrations although no direct transverse forces are applied to these belt spans. Belt-pulley coupling is undesirable because it allows crankshaft excitation, which drives pulley rotation directly, to indirectly drive transverse belt vibrations that cause noise, fatigue failure, and belt slip. A more complete understanding of belt-pulley vibration coupling will allow engineers to design for reduced span vibrations.
To accurately predict the entire system’s dynamic behavior, both the discrete motions of the pulleys and the continuous motions of the belt need to be considered in the model. This results in a coupled continuous-discrete system. Former models concentrate mainly on the pulley rotational motions while treating the belt spans only as springs linking the pulleys [2,4-7]. These discrete models can not predict the system behavior without losing fidelity, especially when coupling between the belt spans and pulleys is strong. Furthermore, they provide no prediction of the troublesome belt span response.
Thus, a comprehensive hybrid model is preferable, in which both the pulley rotational motions and the transverse span motions are incorporated.
Because the transport speed of the belt is high for high engine speeds, this system has the gyroscopic characteristics of an axially moving medium [8-10]. Many industrial applications fall into this class, like band saws, tape drives, paper-handling machinery, textile processing, high-speed printers, and copiers. Due to the mathematical complexities, hybrid continuous-discrete gyroscopic systems have been rarely studied in the literature. Novel mathematical formulations are needed to present the governing equations, and new analytical and numerical methods are needed to solve the problems
2 arising from these kinds of systems (such as finding the steady state mechanics at a specified speed, solving the corresponding eigenproblems, sensitivity analysis, dynamic response, and so on). Thus, developing mathematical methods for hybrid continuous- discrete systems is an objective of this research.
Figure 1.1 Installed serpentine belt drive. 3
J5 Power Steering J4 θ θ 5 4 J J 7 Idler 3 w (x ,t) θ θ i-1 i-1 7 3 Alternator J θ Tensioner t t β 1 ζ w (x ,t) ζ1 2 i i K Water Pump t β θ 2 6
J Air 6 Conditioner Crankshaft θ 2 J θ 1 1 J2
γ
Figure 1.2 Sketch of a seven-pulley serpentine belt drive.
4
For practicing accessory drive designers in the automotive/heavy vehicle industries, the natural frequencies and dynamic response are of utmost interest. Predicting the system’s dynamic behavior in the design stage is important because changes are difficult to accomplish once prototypes are built. This requires a comprehensive hybrid continuous-discrete, gyroscopic model of serpentine belt drive vibration. Thus, a practically important engineering need creates the opportunity for analytical scientific research.
Reducing belt-pulley coupling is one of the main concerns in the design stage.
Undesirable transverse span vibrations reduce belt life, produce noise, and accelerate belt fatigue. Large belt vibrations can even cause the belt to jump a pulley groove, dramatically increasing belt stress that ultimately snaps the belt and fails the system. The mechanisms causing large transverse vibrations in some belt spans of the serpentine system are not clear. Better understanding of those coupling mechanisms can deepen our knowledge of this system and provide useful information to designers of serpentine systems. Correspondingly, development of a mathematical model that can capture the belt-pulley coupling is one of the main objectives of this research. This dynamic model should explain the large belt transverse vibration, which occurs at engine firing frequency, as observed by our research sponsor. Further analytical investigation of design parameters on this belt-pulley coupling mechanism follow naturally.
Another objective of this research is to advance the investigation of steady state contact mechanics for general belt-pulley drives, where serpentine belt drives are one application. This analysis is somewhat different in spirit from the above investigation of
5 belt-pulley dynamical coupling phenomenon. This analysis focuses closely on the detailed steady state mechanics of the belt, the pulley, and their interactions, such as the distribution of friction and belt tension on the belt-pulley contact zones. For this analytical purpose, fewer modeling assumptions are made, producing a more general and complicated model. For example, belt speed and tension are no longer assumed to be uniform throughout the belt, the belt has slip and stick arcs on the pulleys, and the belt- pulley contact points are not known a priori. Only steady state analysis on multi fixed pulley systems is pursued in this research. This steady state belt-pulley mechanics analysis is useful in practice to better predict the stress distribution of the belt-pulley drives, which is one of the dominant factors used to calculate belt fatigue and predict belt durability. It can also better predict some important design performance criteria such as maximum transmissible moment and power efficiency. The analysis along this path can be extended to more complicated situations like free vibration analysis or nonlinear response prediction of entire serpentine belt drives, which is planned for future work.
1.2 Literature Review
In this section, a brief review of the related work is given. The review is divided into three parts. Part 1 focuses on the discussion of the general class of axially moving materials. The belts in the free spans of serpentine belt drives belong to this class. Part 2 discusses papers that address the entire serpentine belt drive. Part 3 reviews work concentrating on the steady state belt-pulley mechanics.
6
1.2.1 Vibration of axially moving materials
The moving belt of the serpentine belt system belongs to the class of axially moving materials. Extensive research has been done for these materials. The topics studied include steady state, free vibration, dynamic response, stability, and so on [8-15]. Most research in this field concerns the transverse vibration of axially moving materials and treats the continua as moving strings or beams with the boundary conditions of the materials as simple supports. This modeling excludes pulley rotations at the ends of the belt span. This results in purely continuous systems and no interaction between the span and the pulleys is considered. See papers [10-12] for comprehensive reviews for these axially moving materials.
For axially moving (gyroscopic) systems, the natural frequencies are velocity dependent and the eigenfunctions are complex [13], which means that the different components of a system do not have the same phase in the vibration modes (that is, they do not pass through the steady state position at the same time even in single mode response). The property results from the transport speed of the material, which makes the system gyroscopic. The axial velocity introduces two convective acceleration terms that are not present in the equivalent stationary system: one is the centripetal acceleration term and the other is the Coriolis acceleration term. The canonical form of a gyroscopic
+ + = system is MU tt GU t KU F where M and K are self-adjoint operators and G is skew-self-adjoint [16-18].
7
One prominent characteristic of this study is the incorporation of the belt bending stiffness into the modeling of the entire serpentine belt drive. The belts in serpentine belt drives have small bending stiffness relative to their tensions. The belts can be modeled as either a traveling string or a traveling, tensioned, Euler-Bernoulli beam with small dimensionless bending stiffness. One contribution of this research is to develop a uniformly valid perturbation method for closed-form approximate eigensolutions of moving, tensioned beams with vanishing bending stiffness. Thus, the works in the literature investigating the transition behavior from the axially moving string case to the beam case is highlighted here. Pellicano and Zirilli [19] study the nonlinear vibrations of large amplitudes for axially moving beams with small bending stiffness. Boundary layers due to the consideration of small bending stiffness appear at the end boundaries, which need the special treatment of singular perturbation techniques (like the method of matched asymptotic expansions). Oz et al. [20] and Ozkaya and Pakdemirli [21] examine the transition from axially moving string to beam for an axially accelerating material. O’Malley [22] obtains the perturbation solution for the eigenvalue problem of a stationary beam for clamped boundaries.
1.2.2 Serpentine belt drives
For serpentine belt drives, considerable research has been done with models that only consider pulley rotational motions. The pulleys are linked by the belt spans as springs, and no belt transverse vibrations are considered. This results in a completely discrete model. This model is relatively simple, non-gyroscopic, and has been studied extensively. Many factors have been incorporated to describe the system dynamic
8 behavior such as bearing damping, Coulomb friction of the tensioner arm, and so on [2,4-
6]. The modeling assumptions exclude the span transverse vibrations, which may be large and interact strongly with the pulleys, so this is only an approximation for the whole serpentine belt system.
The above two models represent two distinct types of vibration found in serpentine belt systems: 1) transverse vibration of individual belt spans, 2) rotational vibration of pulleys with the belts acting as linking springs. In contrast to the above models, there are much fewer papers in the literature that address coupled vibration of the belt spans and the pulley rotations. Models that consider the interactions between the continuous (belt span) and discrete (pulley) components are more complicated and less developed. The study of these hybrid models appeared later than the aforementioned two kinds of models. Initially, only the simpler band saw system was investigated. The band saw consists of two belt spans and two pulleys. Mote and Wu [23] first noticed the coupling phenomena between the belt span vibration and the pulley rotations from experiments.
Wang and Mote [24] established an analytical model to obtain a linear belt-pulley coupling mechanism due to the finite bending stiffness of the metal belt. (Although a lot of research has been done related to band saw systems, most of it treats the pulleys as simple supports for the spans, which results in a purely continuous, single span system decoupled from the pulleys [10,11,13]). Hwang and Perkins [25-27] focused mainly on the response of axially moving beam-like elements at translation speeds that exceed the classical “critical speed stability limit.” Further, they applied their theory to high-speed band/wheel systems and studied the response of the entire hybrid system. Leamy and
9
Wasfy [28] developed a dynamic finite element model to study the response of a simple two pulley belt-drive system by modeling the belt as truss elements while treating the pulleys as rotating circular constraints.
Because the serpentine belt system has more spans and pulleys plus the spring- loaded tensioner, modeling the entire system is more complicated than the two pulley system. Ulsoy, et al. [29] considered the possibility of parametric instability and presented a mechanism which may cause large transverse vibration in some spans due to tension fluctuations under certain operating conditions. Note this model incorporates the pulley-tensioner coupling. In [14,15], this behavior is considered further for the case of parametric instability in a single span. The boundaries of the belt spans are specified, thus excluding the interaction between the belt span vibrations and pulley rotations. In addition to giving a comprehensive review of available models for analyzing the free and forced vibrations of the power transmission belt, Abrate [12] first mentioned the challenging research direction of modeling the entire serpentine belt drives with the consideration of geometric nonlinearity due to the introduction of the tensioner device.
Beikmann, et al. [30,31] treated the belt as a moving string and studied a prototypical three-pulley model that has all essential components of a real serpentine belt drive. He captured a linear coupling mechanism between the tensioner rotation and the transverse vibrations of the two spans adjacent to the tensioner. This coupling results from tensioner rotation moving the boundary points of the two adjacent spans. Zhang and Zu [32] and
Zhang, et al. [33] further refined this linear model by incorporating belt damping and gave a complex modal analysis of the hybrid model for the serpentine belt drive system.
10
Parker [34] developed a spatial discretization of this model extended to n pulleys.
Beikmann, et al. [35] demonstrate a nonlinear coupling mechanism between all spans and pulleys due to finite belt stretching. He studied the two-to-one internal resonance using a numerical method. For the same nonlinear model, Zhang, et al. [36] investigated a one- to-one internal resonance combined with primary external resonance using the multiple scales method.
In practice, transverse vibrations of the belt generally appear on those spans not adjacent to the tensioner and not satisfying the conditions of parametric instability. None of the aforementioned models can predict this behavior. This suggests that belt span transverse vibrations can be excited through as yet unidentified belt-pulley coupling mechanisms. Finding the mechanisms accounting for the belt-pulley coupling, investigating their roles under different operating conditions, and predicting the system’s dynamic behaviors are the main concerns of this research.
After the models are established, solving them analytically or numerically also presents challenges due to the complication of the hybrid continuous-discrete modeling.
Up to now, mathematical methods and tools dealing with hybrid systems are rare in the literature, and only a few papers in the literature relate to such systems. The methods presented in these papers are quite useful. Parker and Mote [37,38] developed an extended operator form to study the vibration of disk-spindle systems where the conventional operator form is expanded to act on an extended variable that includes multiple unknown functions and discrete unknown variables. With proper definition of the inner product, the system is cast as a canonical gyroscopic one and powerful
11 analytical methods can be applied naturally (e.g., Galerkin discretization, perturbation, etc.) For serpentine drives, it is preferable to study the system in an extended operator formulation since the system involves both discrete (pulleys) and continuous (belt spans) components. Beikmann [31] presented the eigenvalue problem in the extended operator form and proved that these operators satisfy the symmetric/skew-symmetric properties for a conservative gyroscopic system. In his free vibration analysis, however, this extended operator form is not used (instead, a modification of Holzer’s method was used). In this research, the extended operator form and its associated properties are used to the full extent.
1.2.3 Belt-pulley steady state contact mechanics
Serpentine belt drives belong to the broad class of belt-pulley drives, which have been widely used to transmit power for hundreds of years. The power is transmitted from the driver pulley to the driven pulleys through friction between the belt and the pulleys.
The belt-pulley contact mechanics are important in industrial applications as they impact belt tension, belt life, power transmission efficiency, maximum transmissible moment, and noise. Considerable research has been done in this field. The earliest work can be traced back to Leonard Euler’s study of a belt wrapped around a fixed pulley or capstan
[39]. Fawcett [40] gives a comprehensive review of belt-pulley contact mechanics up to
1981. Two different theories have been used to describe the belt behavior. One is known as creep theory, which assumes that the belt is elastically extensible, friction is developed due to the relative slip motion between the belt and pulley, and a Coulomb law describes the belt-pulley friction. Another model is the shear theory, which addresses shearing
12 deformation of the belt and assumes that the belt is inextensible. The shear theory is developed recently in [41,42]. Alciatore and Traver [43] give a comparison between these two different theories. In this research the creep theory is adopted with the refinement of incorporating belt bending stiffness.
Creep theory has been widely used to calculate and design belt-pulley drives in industry. Johnson [44] gives a review of the classic creep theory. Recently, by considering inertial effects, Bechtel et al. [45] update the classic creep theory to include belt inertia and present a complete solution for a two-pulley belt drive. Independently,
Rubin [46] investigates the effects of the same inertia terms and presents a method to find solutions for general multi-pulley systems. Although the derivations in [45,46] seem different, the analysis and main conclusions are essentially the same. The main contribution of these two papers is that they include belt inertia terms and determine the relative errors of prior creep theories that neglect these terms. Belt bending stiffness is ignored in both studies.
1.3 Scope of Investigation
The main scope of this project is to advance the knowledge and understanding of serpentine belt drive mechanics. A prominent characteristic of this research is the incorporation of bending stiffness into the modeling of the axially moving belt. The background information and literature review is introduced in chapter 1. Chapter 2 analytically explores the transition behavior of modal properties from an axially moving string to a tensioned, axially moving, Euler-Bernoulli beam with small bending stiffness.
An approximate perturbation method is developed. Closed-form, uniformly valid 13 approximate eigensolutions are derived for combinations of different boundary conditions. Chapters 3 and 4 investigate the steady state and dynamic behavior of the entire serpentine belt drive when the belt bending stiffness is considered. This part of the analysis is based on the assumption of fixed belt-pulley contact boundaries. Belt-pulley mechanics on the contact arc are not considered in detail. Chapter 3 first establishes the model. A computational method based on boundary value problem solvers is developed to obtain the numerically exact solution of the nonlinear steady state equations. An approximate analytical solution of closed-form is also obtained for the case of small bending stiffness. Based on these solutions, the effects of design variables on the steady state deflections and span-pulley coupling are investigated. Chapter 4 then studies the linear vibration about the obtained nontrivial steady state. A robust spatial discretization is presented to solve the corresponding eigenproblem. New dynamic characteristics of the system induced by belt bending stiffness are discussed. Belt-pulley dynamic coupling is investigated thoroughly through evolution of the vibration modes. Chapter 5 focuses on the investigation of steady state contact mechanics for belt drives of fixed pulley systems.
More complicated factors such as belt inertia, elastic extension, and Coulomb friction are retained, and belt bending stiffness is also included. A numerical iteration method based on general purpose BVP solvers is developed. This overcomes the obstacles induced by inclusion of the bending stiffness (e.g. non-uniform distribution of the tension and speed in the free belt spans, unknown belt-pulley contact points a priori, etc).
14
The benefits of this research include:
1. A uniformly valid approximate perturbation technique is developed for the eigenvalue
problem of axially moving beams with small bending stiffness. The closed-form
approximate eigensolutions for different boundary conditions can be used as a basis
to analyze nonlinear and parametric excited moving media problems. The knowledge
applies to most of the belts in belt-pulley drives, that is those with small bending
stiffness relative to tension.
2. The developed model for the entire serpentine belt drives is suitable for the
explanation and understanding of the belt-pulley coupling phenomena observed in
experiment. Belt-pulley coupling induced by the consideration of bending stiffness is
thoroughly investigated and can be used in troubleshooting large belt vibration caused
by this coupling mechanism.
3. Knowledge of the natural frequencies and vibration modes is essential for almost all
dynamic analyses. This study characterizes the special properties that are applicable
for serpentine belt drives with bending stiffness, which is a fundamental advance for
serpentine belt vibration research.
4. The research of steady state contact mechanics between the belt and pulley reveals
how the belt bending stiffness alters the distribution of belt tension, speed, and
friction in the belt-pulley drives. Further some important design criteria like
maximum transimmible moment and power efficiency are also analytically obtained.
15
5. Effective methods, both theoretical and numerical, are developed to address the
steady state, eigenvalue, and dynamic response problems arising from the above
investigation. Some of these methods can be generalized to other coupled continuous-
discrete systems.
6. A comprehensive simulation tool for the dynamics of serpentine drives based on this
analysis has been developed for the analysis of practical serpentine belt drives.
Applied in practice by our research sponsor, it provides practical design guidance to
designers of serpentine drive systems.
16
CHAPTER 2
APPROXIMATE EIGENSOLUTIONS OF AXIALLY
MOVING BEAMS WITH SMALL BENDING STIFFNESS
For the simplest case, without considering the dynamic interaction between the pulley and belt, each free span of a serpentine belt drive can be modeled as an axially moving string or a tensioned, axially moving beam with small bending stiffness. In this chapter, a perturbation method based on the phase closure principle is developed to find the closed-form approximate eigensolutions of an axially moving beam with vanishing bending stiffness. This method is suited for different combinations of boundary conditions and uniformly valid approximate eigenfunctions are obtained. Clear physical interpretation based on wave propagation and reflection is available from this perturbation method.
17
2.1 Introduction
Before studying the belt interaction with the pulleys, it is beneficial to first understand the belt transition behavior from the known axially moving string case to the beam case. In the simplest sense, the transverse vibration of each span of serpentine belt drives is typically modeled as either a traveling string or a traveling, tensioned, Euler-
Bernoulli beam. This chapter focuses on the transition of modal properties between these two cases for an axially moving material.
Closed-form solutions for the natural frequencies and vibration modes are available for the string model [8,47]. For the axially moving beam model, due to the beam’s dispersive property, only numerical solutions are available [13,48]. Because most axially moving media have small bending stiffness relative to their tension, they can be modeled as an axially moving beam with small dimensionless bending stiffness. The transition of modal properties from the axially moving string case to the beam case is desirable from both practical and theoretical viewpoints. Finding closed-form approximate solutions of the eigenvalue problem for such transitional systems is the main objective of this chapter.
Before highlighting the works in this field, we first review work related to the phase closure principle, which is one of the main tools used in this chapter.
The vibration of elastic structures can be described in terms of waves propagating and attenuating in structures. The phase closure principle [49] states that if the phase change for propagating (or evanescent) waves is an integer multiple of 2π as they return to their start point after traveling forward and back along a finite structure, then the frequency at which the waves travel is a natural frequency and the corresponding
18 vibration mode is the superposition of the component waves. In the field of acoustics, solids, and fluids, wave propagation and attenuation in waveguides and wave reflection/transmission at a boundary point has been studied extensively. The phase closure principle links the knowledge of wave motion in these fields to computation of natural frequencies and modes of finite structures. Mead [50] applies the method to find the eigensolutions of stationary beams. Exact frequency equations are established that differ from the conventional ones but have identical roots. These frequency equations have clear physical interpretation and deepen understanding of the beam vibration modes.
Mace [51] develops a numerical matrix method based on wave propagation, reflection, and transmission at a point support (or a change of cross-section or material property) to calculate the natural frequencies and modes for beams. Tan and his coworkers [52,53] extend this method to some complex beam structures, like those consisting of several different uniform segments.
For axially moving continua it is well known that the vibration modes can be viewed as the superposition of pairs of opposite–going propagating waves. The phase speeds in the opposite directions are different due to the convective effect of the medium’s axial speed. Lengoc and McCallion [54] study the relation between wave propagation and natural frequency, but their work is limited to non-dispersive system like taut strings. Lee and Mote [55] investigate the energy transfer due to the interaction between the translating continua and its boundary supports. The phase closure principle is used to obtain the natural frequencies of an axially moving string. Chakraborty and
Mallik [56] study the free vibration of a traveling beam simply supported at both ends.
19
The frequency equation is derived based on the phase closure principle. This work applies only to beams with finite bending stiffness and zero tension; transition behavior from a taut string to a tensioned beam is not investigated.
More commonly, researchers investigate axially moving continua mathematically without consideration of the physical wave propagation. They study transition behaviors for moving beams with vanishing bending stiffness by using perturbation techniques directly on the differential equations. Because the main concern of this chapter is the modal properties, only those aspects of related works are reviewed here. Pellicano and
Zirilli [19] study axially moving beams with simple supports at both ends. While not specifically addressed, their natural frequencies can be extracted from the results. These natural frequencies depend only on the displacement boundary condition of each end, suggesting that the remaining two beam boundary conditions do not affect the natural frequencies. Oz et al. [20] and Ozkaya and Pakdemirli [21] examine the transition from axially moving string to beam for an axially accelerating material. By letting the accelerating terms vanish, the free vibration solutions for constant belt speed follow. In
[20], multiple scales perturbation is applied to find the approximate natural frequencies.
Problems are apparent because no boundary conditions are considered in the derivation, indicating that different boundary conditions yield the same natural frequencies. The problem considered in [21] is similar to that in [20] and similar techniques are used. The improvement is that the spatial boundary layer terms arising from small bending stiffness are considered. Two sets of boundary conditions are considered. The solutions in [21] incorrectly imply that the approximate natural frequencies for these two sets of boundary
20 conditions are the same. Further, for clamped boundaries the zero speed solution fails to give the approximate solution given by O’Malley [22]. In contrast, the present analysis gives different natural frequencies for these two kinds of boundary conditions, and the zero speed results converge to the exact solution (simply supported) and that given by
O’Malley (fixed-fixed), although the adopted methods differ. O’Malley’s work [22] treats stationary beams with two clamped ends. When we extended this method to axially moving beams, the procedure became cumbersome and no explicit solutions were obtained.
In this chapter, a different perturbation method is developed to find closed-form, approximate eigensolutions of axially moving beams with small bending stiffness. Wave propagation considerations lead to an algebraic equation with a small dimensionless bending stiffness parameter. Taking advantage of the simplicity of the propagation and attenuation properties of the waves, which are determined by the roots of an algebraic equation, the phase closure principle is used to find the natural frequencies. The complex vibration modes are obtained naturally from the superposition of all component waves in the beam. Approximate eigensolutions for different boundary conditions are presented.
The perturbation solutions are confirmed by comparison with numerically exact ones. For the special cases mentioned above where the exact or approximate solutions are available, the derived approximate solutions agree with them.
Instead of considering spatial and temporal variations for the governing partial differential equation (like in [19,21]), this approach focuses on perturbation of algebraic equations. No boundary layers or secular terms need to be considered explicitly in the
21 derivation. Although the method is simple, no completeness of the solutions is sacrificed.
For example, the evanescent wave components (if not zero) automatically generate boundary layer terms for those beams where small bending stiffness creates edge effects at the boundaries. Unlike prior perturbations where assumed mode spatial expansions are only suited for certain boundary conditions, this method handles different boundary conditions with a consistent treatment.
2.2 Model Equations
The dynamic equation for an axially moving beam is
+−−+=<<2 mwmcwPmcwEIwtt2() xt xx xxxx 0,0 xL (2.1) where m is the belt mass per unit length, w(x,t) is the transverse displacement, c is the belt transport speed, P is the tension, and EI is the bending stiffness. The following non- dimensional variables are introduced,
x wPEIm xwttˆˆ===,,ˆ ,ε 2 = , vc ˆ = (2.2) L LmLPL22P
Substitution of (2.2) into (2.1) leads to the dimensionless equation (after dropping the hat)
+−−+22ε =<< wtt2vw xt(1 v ) w xx w xxxx 0, 0 x 1 (2.3)
Assuming we= i()rx−ω t , where r is the wavenumber and ω is the wave propagation frequency, equation (2.3) yields
εωω24rvrvr+−(1)2 2 2 + − 2 =0 (2.4)
Note that the small parameter ε <<1 multiplies the highest power of r. The roots of such a polynomial equation have two possible forms [57].
22
In the first form, the roots of (2.4) are expressed using the straightforward expansion
=+εε +2 + rx01 x x 2... (2.5)
Substitution of (2.5) into (2.4) leads to the ε 0 order result
ω x = (2.6) 0 v ±1 and the ε 1 order result
−+=2 ω = 2(1vx )0 2 v 0 or x1 0 (2.7)
= ε 2 The first equation in (2.7) is discarded because it contradicts (2.6), so x1 0 . The order equation gives
ωω3 −=1 4 , when x0 −x4 2 (1)v ++v 1 ==0 x2 (2.8) 22(1)vvxω +−2 ωω3 − 0 1 = 4 , when x0 2(1−−vv ) 1
Equations (2.5)-(2.8) provide two of the four roots of (2.4).
The remaining two roots are expressed as the singular expansion
y rx=++..., λ >0 (2.9) ε λ 0
Substitution of (2.9) into (2.4) gives
yyy εωω2422( ++xvxvx...) +− (1 )( ++ ...) + 2 ( ++ ...) − 2 =0 (2.10) εεελλλ000
The dominant terms in (2.10) are y4 / ε 4λ−2 and (1)/− vy2 22ε λ . Balancing these leads to
4λλλ−=22⇒ =1 (2.11)
23
Equation (2.10) then becomes
εε−−24+− 22 + 1 3 + − 2 + ω += [yvyxyvxyvy(1 ) ] [400 2(1 ) 2 ] ... 0 (2.12) with the solutions
y =±iv1 − 2 (2.13)
vω 2xv [(1−−−22 ) 2(1 v )] + 2 vω = 0 ⇒ x= (2.14) 001− v2
Case Results r1 r2 r3 r4
Exact 0.7222 -6.4335 2.8556+60.2421 i 2.8556-60.2421 i ε = 0.01 v = 0.80 ω =1.30 Approx. 0.7222 -6.4313 2.8889+59.9999 i 2.8889-59.9999 i
Exact 0.6482 -1.1501 0.2510+19.2243 i 0.2510-19.2243 i ε = 0.05 v = 0.28 ω = 0.83 Approx. 0.6482 -1.1501 0.2522 +19.2000 i 0.2522 -19.2000 i
Exact 1.2949 -1.7342 0.2197+10.0070 i 0.2197-10.0070 i ε = 0.10 v = 0.15 ω =1.50 Approx. 1.2947 -1.7324 0.2302 +9.8869 i 0.2302 -9.8869 i
Table 2.1. Comparison of approximate roots of (2.4) from (2.15) with numerically exact roots. 24
In summary, the wave dispersion equation (2.4) has the four roots
ωω1133 ωω rOrO=−εε2 +( 323)( =−+ εε +) 1212(1)++vv44 12(1) −− vv (2.15) vvωω11−−22 vv riOriO=+ +()εε =− + () 3411−−vv22εε
Table 2.1 compares the approximate and numerically exact roots for three cases.
Consistent with (2.15), r1 and r2 are best approximated by perturbation. Physically, the real parts of the roots represent the phase change between two points unit distance apart, and the imaginary parts represent the variation of the wave amplitudes for two points unit
distance apart. Specifically, r1 represents the wave propagating in the positive direction,
r2 the wave propagating in the negative direction, r3 the evanescent wave attenuating in
the positive direction, and r4 the evanescent wave attenuating in the negative direction.
The beam motion is the superposition of the four components
− ω =+++ir12 x ir xir3 x ir 4 x it w(,)xt [ ce12 ce ce 3 ce 4 ] e (2.16)
where cc14~ are complex coefficients.
ε << For small bending stiffness 1, the imaginary parts of the evanescent waves r3
and r4 become very large. Consequently the r3 component can exist only close to the
= = boundary x 0 , and the r4 component exists only close to the boundary x 1 (Figure
2.1). They can be viewed as part of the reflected waves as the propagating waves ( r1 and
= r2 ) travel forward and back along the beam between the boundary points A( x 0 ) and
B( x =1) (Figure 2.1). The phase closure principle can now be applied to the propagating waves to find the eigensolutions for different boundary conditions. 25
x
r3 r r1 v 1
A x=0 Beam with small ε = EI/PL2 x=1 B
r2 r2
r4
Figure 2.1. Waves in a finite moving beam with small bending stiffness
2.3 Application of the Phase Closure Principle
To apply the phase closure principle to the propagating waves ( r1 and r2 ), one needs to find four different phase changes: first, the propagating wave leaves from
= boundary A and arrives at boundary B with a phase change Re (rr11) (because the span is normalized to unit length); second, it reflects at boundary B with a phase change
φ ()RB ; third, it travels from boundary B to boundary A with another phase change -
26
=− Re ()rr22 (minus sign due to leftward propagation); finally, it reflects at boundary A
φ with the fourth phase change ()RA and returns to the start point boundary A.
Mathematically, the phase closure principle requires
+−+=φφπ =±± rRrR12( BA)()2,0,1,2, nn ... (2.17)
Consider the case of a simply supported beam with boundary conditions
== == wt(0,)(0,)0,(1,)(1,) wxx t wt w xx t 0 (2.18)
= At x 0 (point A in Figure 2.1), there is no r4 evanescent component as noted above, and
(2.16) becomes
− ω =++ir12 x ir x ir3 x it w(,)xtA [ ce12 ce ce 3 ] e (2.19)
Substitution into the boundary conditions at x = 0 yields
11c 1 c 1 rr22− 1 =− c ⇒ 1 = 23c (2.20) 22 22 22− 22− 2 rr13 c 3 r2 c3 rr31rr12
For the two propagating waves, the relative phase due to the reflection at the left boundary is given by the phase of
c rr22− {−− (1v2 ) +εω22222 [vv /(1 − )] − εω [ /(1 − vivv )] } + 2 εω / 1 − 2 R ==−1 32 =− (2.21) A 22− 2 22222 2 crr231{−− (1v ) +εω [vv /(1 − )] − εω [ /(1 + vivv )] } + 2 εω / 1 −
One can prove mathematically that the phase angle of RA is
φ =+πε3 (ROA )() (2.22) as shown graphically in Figure 2.2, where π is from the leading minus sign preceding
2 −−222 (r32rrr)/( 31) in (2.21).
27
2 Im D2= O( ε )
D1= O( ε) Denominator of D = -RA Ο(ε3)
Numerator of D = -RA = D3 O(1) Re
=− =2 −222 − Figure 2.2. Phase angle of D RrrrrA ( 3 231)/( ) from equation (2.21), =−εω 2 =−−+εω2 22 ω D1 2/1vv, Dvv2 {[/(1)][/(1)]}, =− −2 +εω2222 − − ω + Dvvv3 (1 ) {[ /(1 )] [ /(1 v )] } .
= At x 1(point B), there is no r3 evanescent component in (2.16). In seeking the relative phase between propagating waves at B, it is notationally convenient to introduce
ξ =− = irk x 1 and express the coefficients in (2.16) using bkkce . This gives
ξξξ − ω =++ir12 ir ir3 it w(,)xtB [ be12 be be 4 ] e (2.23)
Substitution of (2.23) into the x =1 boundary conditions yields
11b 1 b 1 rr22− 2 =− b ⇒ 2 = 14b (2.24) 22 21 22− 22− 1 rrb24 4 r 1 b4 rr42rr21
For the two propagating waves, the relative phase due to the reflection at the right boundary B is given by the phase of
28
b rr2 − 2 R ==2 14 (2.25) B 2 − 2 b1 rr42
Similar to the handling of RA , the phase angle of RB is
φ =+πε3 (ROB )() (2.26)
Substitution of (2.15), (2.22), and (2.26) into (2.17) leads to
21ω 1 1 −++==±±εω23[ ] ... 2nn π , 0, 1, 2,... (2.27) 12(1)(1)−+−vvv244
This is an algebraic equation (for ω ) with the small parameter ε multiplying the highest power. Application of the previously discussed algebraic perturbation technique leads to three different roots. Only the root from the straightforward expansion form is retained.
The two roots from the singular expansion form are discarded because they yield complex roots, and physically the natural frequency ω must be real for subcritical
ωω=+ εωεω +2 + speeds. Substitution of 01 2... into (2.27) leads to
ωπ=−+222242 επ +++ = n nv[1 nvv ( 6 1) / 2 ...], n 1,2,3,... (2.28)
Figure 2.3 compares the fundamental ( n =1) natural frequencies obtained from
(2.28) with numerically exact solutions for different belt speed v and different bending stiffness. The approximation results are best for small bending stiffness and low axial belt speeds. This is because for such cases the four roots in (2.15) have the best perturbation approximation. For large bending stiffness or high speed, more terms need to be incorporated in the perturbation approximation.
29
4
3.5 v=0 3 v=0.3 2.5
2 v=0.6 1.5 perturbation
1 v=0.9 exact 0.5 Fundamental natural frequency 0 0 0.02 0.04 0.06 0.08 0.1 0.12 Bending stiffness, ε
Figure 2.3. Comparison of fundamental natural frequency for a simply supported beam. −−−, perturbation; , exact.
When v = 0 , (2.28) becomes
1 ωπ=+nnn[1 επ22 ( ) ] += ..., 1,2,3,... (2.29) n 2
The exact eigensolution for the special case v = 0 is
30
ωπεπ=+22 4 = π = exact (nnwxnxn)(),()sin(),1,2,3,... (2.30)
Expansion of the eigenvalue in (2.30) for small ε yields (2.29).
Computation of the eigenfunctions requires additional consideration of the evanescent waves at the boundaries A and B. From (2.20) and (2.24),
c rr22−+−−[/(1)][/(1)]ωω v 2 v 2 R ==3 12 = (2.31) A 22− 2 222 2 2 crr231{[vωεωωε /(1−−−−vv )] (1 ) / [ /(1 ++ vivv )] } 2 /( 1 −)
brr22−−+−−[ωω /(1 v )] 2 [ /(1 v )] 2 R ==−412 = (2.32) B 22− 2 222 2 2 brr142{[vωεωωε /(1−−−−−−vv )] (1 ) / [ /(1 vivv )] } 2 /( 1 −)
The eigenfunction can be written as
ξ =+irx1 irx24 ++ir3 x ir w()xcececebe12 3 4 (2.33)
Normalization of (2.33) by dividing it by c2 , application of (2.21), (2.31), (2.32), and use
ξ =− = ir1 of x 1 and b11ce give
− =+++ir12 x ir x ir3 x ir 14 ir( x 1) wx() ReAAAB e Re Re Re (2.34)
≠ When v 0 , there are boundary layer terms from the evanescent r3 and r4 terms. But
= =− == when v 0 , RA 1 and RRAB0 , leading to the eigenfunctions
w()xee=−irx1 + ir2 x (2.35)
Substitution of (2.15), (2.29), and v = 0 into (2.35) gives the eigenfunction approximation
wx( )sin(),1,2,3== nπ x n ... (2.36) in agreement with the exact solution in (2.30). The eigenfunctions have no boundary layer terms for v = 0 while they do for v ≠ 0 .
31
2.4 Other Boundary Conditions
The above perturbation method for the eigenvalue problem can be applied to other