Dynamic Model of a Badminton Stroke
|
|
- Sandra Thomas
- 6 years ago
- Views:
Transcription
1 ISEA 28 CONFERENCE Dynamic Model of a Badminton Stroke M. Kwan* and J. Rasmussen Department of Mechanical Engineering, Aalborg University, 922 Aalborg East, Denmark Phone: / Fax: mak@ime.aau.dk Topics Abstract The purpose of this paper is to contribute to the understanding of dynamic racket behavior in badminton. The paper describes the development of a dynamic model of a badminton stroke based on experimental data. Motion capture and strain gage experiments are performed to clarify the movement and accelerations of a badminton smash stroke. Subsequent data processing reveals that due to the speed and acceleration of the movement, neither of the two methods provides the complete picture of the racket s dynamic behavior, but computational processing of the combination of experimental data together with physical modeling of the racket mechanics allows for an understanding of the basic mechanisms underlying racket movement and deformation. The paper concludes that... Keywords 1 Introduction Although badminton racket technology has advanced considerably, the relationship between design and response are still rather poorly understood. A review of the literature turns up very little on badminton rackets, whereas design of tennis rackets has been studied quite extensively (Brody, 1995; Cross, 1999; Lees, 23). Some tennis studies could be applicable to badminton as well as other racket sports in general, but given the distinct styles of play, the design goals of a tennis racket and a badminton racket are probably quite different. As the world s fastest racket sport, the game of badminton is very much about speed. With smashes reaching up to 332 km/h, the racket is subjected to significant dynamic effects. A dynamic model is therefore necessary to properly analyze racket performance. During the stroke, the racket experiences transverse deflection due to the high rotational and translational accelerations and it is likely that this deflection plays an important role in the transfer of momentum to the shuttle cock. Skilled players perceive significant differences between different high-quality rackets that superficially appear similar in terms of elasticity and mass properties. This indicates that small differences in the racket design have a high influence on the dynamics properties. An understanding of racket dynamics is therefore essential for further improvement of racket design. 2 Motion Capture Experiment *** JR: The section below should be elaborated more: Precisely how are the numbers processed? Where does the inaccuracy appear? What has been done in terms of filtering and such? Motion capture data were taken for ten trials of a smash stroke performed by an advanced player using a Qualisys system of eight high-speed cameras at the maximum frame rate of 24 Hz. Six markers were placed on the racket, two on the handle and four on the head, illustrated in Figure. A 12 mm spherical marker was fixed at the base, and reflective tape was used for the other five markers. Position data of the markers were passed through a low-pass Butterworth filter with a cut-off of 35 Hz and subsequently interpolated using a 6th order B-spline. An optimization-based method designed Paper No. -1- Copyright of ISEA 28
2 ISEA 28 CONFERENCE for over-determinate systems was then applied to determine the kinematics of the system (Andersen et al.). The racket was modeled as a flexible multi-body system with 7 dof, where the handle was considered a rigid segment and the revolute joint between the handle and shaft allowed transverse deflection. The results of the kinematic analysis are shown in Figures and. Figure 1: Illustration of racket anatomy and racket marker placement for motion capture. Due to the high accelerations involved and the limitations of the motion capture system used, the accuracy of marker positions may have been compromised. As seen in Figure, the error between measured and calculated marker positions are below 4 mm during most of the stroke. However, the errors increase up to nearly 2 cm around the time just before impact where accelerations are presumably the highest. The effects of the error can be seen in the calculated angular velocities and accelerations in Figure. Angular kinematics are given with respect to the local reference frame, defined in Figure. The angular accelerations are noisy and somewhat high, with peak accelerations around 15 rad/s 2. 2 p p 2 p 3 Error [cm] 1 p 4 p 5 p Figure 2: Error between calculated and measured marker positions during pre-impact part of stroke. Impact occurs just before t=.3 s. Paper No. -2- Copyright of ISEA 28
3 ISEA 28 CONFERENCE α [rad/s 2 ] ω [rad/s] Figure 3: Angular velocities [rad/s] and accelerations [rad/s 2 ] of the racket about its axial x-axis (blue), y-axis (green), and z-axis (red) during pre-impact part of stroke. Impact occurs just before t=.3 s. The resulting load on the racket was then computed, despite the questionable accuracy of the kinematics derived above. The load is proportional to the acceleration normal to the racket, containing contributions from both linear and angular accelerations. The normal acceleration generated at the tip of the racket is shown in Figure. At the tip of the racket, the normal acceleration is dominated by the rotational component, implying that the angular accelerations contribute much more to the speed of the racket tip than the linear accelerations. The acceleration is especially noisy during backswing, from t= to.2 s, followed by an increase (in the negative direction) as the racket accelerates forward. As the racket continues forward, unloading of the racket begins (t=.25 s) as the normal acceleration decreases back to zero Linear Rotational Total Acceleration [m/s 2 ] Figure 4: Contributions to total normal acceleration of racket tip (red) from translational accelerations (blue) and rotational effects (green) during pre-impact part of stroke. Impact occurs just before t=.3 s. Paper No. -3- Copyright of ISEA 28
4 ISEA 28 CONFERENCE 3 Strain Gage Experiment The motion capture experiment revealed that the high accelerations in the motion combined with the spatial and temporal resolution limitations of the system makes it difficult to obtain reliable values for the accelerations, and thus for the dynamic forces in the system. An alternative is to measure a more direct representation of the forces in the system, namely the strain in racket shaft. Two strain gages were placed on the end of the shaft 8 mm from the top cap in a half-bridge configuration. Ten trials of a smash stroke were performed by the same player and the same racket from the motion capture trials. Through a series of simple bending tests for a clamped racket with a concentrated end load, a linear relationship between the strains and the measured deflection at the tip of the racket was determined. Using this relationship between strain and end deflection for a concentrated load, the end deflection assuming a triangular distributed load was calculated as a function of the strain, shown in Figure Deflection at Tip [m] Strain [µ m/m] Figure 5: Linear relationship between measured strain at shaft base and measured deflection at tip for concentrated load (blue). Calculated tip deflection for an increasing triangular distributed load (red). Applying the linear relationship to the smash stroke data gives the tip deflection curves seen in Figure. Over a stroke, the maximum strain of around 8 µm/m corresponds to a tip deflection of about 5 cm. Paper No. -4- Copyright of ISEA 28
5 ISEA 28 CONFERENCE.6.4 Deflection at Tip [m] Figure 6: Deflection at racket tip over one full stroke. Impact occurs around t=.46 s. During the stroke, the racket is first swung backward and then forward before hitting the shuttle. The backswing generally takes under.2 s, and the forward stroke lasts around.1 s. In Figure, a negative deflection of the racket can be seen (around t=[.14,.34] s) during the backswing (positive normal acceleration), followed by a much greater positive deflection during the forward stroke (negative normal acceleration) until impact around t=.46 s. Figure 7: Illustration of racket deflection during backswing (left) and forward swing (right) before impact. For the maximum benefit from the elastic deformation of the racket, impact should occur when the deflection of the racket returns to zero, where the tip velocity will be highest. However, the elastic deformation will still be advantageous as long as the racket tip is moving forward (towards the shuttle), indicated by a negative deflection slope in this case, so there is a narrow window of about.1 s where slight deviations in the timing of the impact can be tolerated. As perhaps a more reliable source than the motion capture data, the strain gauge data was also used to derive racket kinematic. The racket was assumed to experience rotation only, due to an applied torque at the handle. The magnitude of the torque τ was computed from the strain ɛ, given by τ = EIɛ/r where EI is the flexural rigidity of the shaft and r is the shaft radius. The angular acceleration of the racket due to the torque was determined from α = τ/(i cm + md 2 ) where I c m is the racket moment of inertia, m is the total mass, and d is the distance from the strain gauge to the racket center of mass. The resulting normal acceleration of the racket tip is given in Figure. Paper No. -5- Copyright of ISEA 28
6 ISEA 28 CONFERENCE 4 2 Acceleration [m/s 2 ] Figure 8: Estimated normal acceleration of racket tip during pre-impact part of stroke, calculated from strain measurements. Impact occurs around t=.35 s. Comparing with the normal acceleration derived from motion capture data in Figure, the shape of the curves during the forward swing phase are quite similar, with the peak acceleration near 8 m/s 2. 4 Computational model *** JR: Possibly a shortened version of the Lagrange model leading to the conclusion that the model does not capture the correct dynamics, possibly due to lack of damping and incorrect boundary conditions. The simplified racket model is composed of two uniform beam elements, representing the shaft and the head of the racket. Each element has uniform mass per unit length, ρa, and uniform flexural rigidity, EI. The uniform beam assumption is a good approximation for the shaft, but is a rather gross simplification of the racket head. The handle is assumed to be completely rigid. The origin of the body-fixed coordinate system xyz of the racket, located at the base of the handle, can be expressed in the ground coordinate system XY Z by the vector r, as seen in Figure. The rotation of the racket coordinate system is related to global by the transformation matrix T, such that xyz = T XY Z. Figure 9: Beam model (one element). Paper No. -6- Copyright of ISEA 28
7 ISEA 28 CONFERENCE In addition to the six degrees of freedom from rigid body motion, the racket is allowed to bend in the transverse direction. The finite beam element of length L is represented by two nodes, one at each end. The displacement vector q contains the linear and rotational displacements at each node, q(t) = [v 1 θ 1 v 2 θ 2 ] T, and the displacements along x between the nodes are described by shaping functions, N(x) = [N 1 N 2 N 3 N 4 ]. The deflection of the beam w is then given by w(x, t) = N(x)q(t). Dynamic equations describing the flexible motion of the racket due to both inertial and centrifugal effects are derived using Hamilton s principle (e.g. Yang, 24). Applying Lagrangian mechanics, the kinetic energy of a beam element is given by T = 1 L 2 ρa r r dx (1) where r is the position vector (in local coordinates) of a point P on the racket from the global origin. Since the racket experiences only tranverse deflection, the position vector of point P from the local origin is p = R + x, w,, where R is the length of the rigid handle. Then r = T r + p (2) r = T r + ω p + p (3) where ω contains the angular velocities of the racket about its body-fixed axes. The potential energy consists of elastic energy and energy from centrifugal force. Neglecting axial deformation, the potential energy is given by V = 1 L 2 EI (w ) 2 dx + 1 f c (x) (w ) 2 dx (4) 2 where f c (x) is the centrifugal force (per unit length) that arises from rotation, given by L L L f c (x) = ρa a c (ζ) dζ = ρa x x r i dζ (5) where i is a unit vector along the x-axis of the racket. The centrifugal acceleration a c does not include deflection terms, such that r = T r + ω p (6) where p = R + x,, and r = T r + p represent the undeformed position vectors of point P from the local and global origin, respectively. Applying Lagrange s equations, L q d dt ( ) L = (7) q where the Lagrangian L = T V, the non-linear equations of motion take the form: [m] { q} + ([k] [k 1 ] + [k 2 ] + [k 3 ]) {q} = [c] T + [c 1 ] T (8) where m and k are the mass and stiffness matrices of the beam element. Dynamic stiffening is represented by the stiffness matrices k 1, k 2, and k 3, and the forcing functions due to rotational and translational accelerations are given by the vectors c and c 1, respectively. The equations were solved using a fourth order Runge-Kutta integration. Mass, geometric, and stiffness properties for the shaft and head segments were determined experimentally, shown in Table 1. Using these values, a natural frequency of 14.9 Hz was calculated for the beam model in the clamped condition, which corresponded well with the racket s natural frequency of 15.1 Hz. Paper No. -7- Copyright of ISEA 28
8 ISEA 28 CONFERENCE Table 1: Material properties of beam segments. Mass Length ρa EI Segment [kg] [m] [kg/m] [N-m 2 ] Shaft Head Using kinematic inputs from the motion capture data, the beam model produces the deformation curve seen in Figure. In the clamped case, the predicted deflections are somewhat high, but the frequency of the dynamic response is also too high. One possible explanation is the clamped boundary condition, which may not be very representative of the handheld case (Cross; Brody). To loosen up the boundary conditions, the end of the racket was changed from clamped to pinned with a torsional spring of stiffness k t. A very stiff spring (k t >1e6 N-m/rad) approaches the clamped condition, and lowering the stiffness leads to a reduction in the response frequency. Figure shows the elastic deflection at the tip, over one stroke with k t =5 N-m/rad..1.8 Clamped Pinned spring Deflection at Tip [m] Figure 1: Deflection at beam tip during pre-impact part of stroke for clamped condition (blue) and pinned with spring k t =5 N-m/rad (green), using motion capture data. Impact occurs just before t=.3 s. Figure shows the deflections of the racket tip using kinematic inputs from the strain gauge data, for the clamped condition and the pinned case with a torsional spring of stiffness k t =5 N-m/rad. In both Figures and, the deflections of the racket tip are still rather high, when compared to the values found from strain gauge measurements, shown in Figure. This could be caused by an overestimation of the accelerations from the experimental data and/or a lack of damping in the model. Paper No. -8- Copyright of ISEA 28
9 ISEA 28 CONFERENCE.2.15 Clamped Pinned spring Deflection at Tip [m] Figure 11: Deflection at beam tip during pre-impact part of stroke for clamped condition (blue) and pinned with spring k t =5 N-m/rad (green), using motion capture data. Impact occurs around t=.35 s. 5 Conclusions 1. Our current mocap technology is only partially useful. Is it possible that state-of-the-art technology will be enough? What frame rate and spatial resolution will be sufficient, based on gaps in our mocap results? 2. The strain gage experiment reveals surprisingly large elastic deflections. Is this reasonable? 3. What can we say about the timing of the stroke? How sensitive is elastic advantage to timing? Motion capture could be better with spherical markers drilled into racket to improve spatial resolution and higher frame rate to improve temporal resolution and avoid disappearing markers. (how high is necessary?) The elastic energy of the deflected racket is advantageous as long as impact occurs while the racket is continuing to rotate forward. The ideal timing is to wait long enough for the racket to recover from the backwards deflection until it is vertical again, although ±.5 s is the window of opportunity where the elastic energy would still be beneficial. what is this frequency dependent on, the stroke (kinematics) or the racket (material)? References [A1] Andersen, M.S. [B1] Brody, H. How would a physicist design a tennis racket? In Physics Today, 48(3):26-31, [C1] Cross, R. Impact of a ball with a bat or racket. In American Journal of Physics, 67:692-72, [L1] Lees, A. Science and the major racket sports: a review. 21:77-732, 23. In Journal of Sports Sciences, Paper No. -9- Copyright of ISEA 28
10 ISEA 28 CONFERENCE [YJ1] Yang J.B., Jiang L.L., and Chen D.C.H. Dynamic modelling and control of a rotating Euler-Bernoulli beam. In Journal of Sound and Vibration, 274: , 24. Churchyard. The dynamic behavior of a rotating flexible beam based on finite beam elements has been well documented, e.g. Yang. Montagny et al. developed such a model and used a badminton racket as a case study, which appears to be the only study related to badminton racket dynamics. However, these rotating beam models are generally limited to planar motion and rotation about a fixed point. Here the racket is allowed to translate and rotate in 3D space. Appropriate boundary conditions to describe the handheld racket are also discussed. Here the racket is represented as a two-segment Bernoulli-Euler beam, and the focus is on understanding what dynamic effects take place over the stroke before impact, i.e. deformation of the racket. *** Paper No. -1- Copyright of ISEA 28
Available online at ScienceDirect. The 2014 conference of the International Sports Engineering Association
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 72 ( 2014 ) 496 501 The 2014 conference of the International Sports Engineering Association Dynamic contribution analysis of
More informationModule 4 : Deflection of Structures Lecture 4 : Strain Energy Method
Module 4 : Deflection of Structures Lecture 4 : Strain Energy Method Objectives In this course you will learn the following Deflection by strain energy method. Evaluation of strain energy in member under
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationCO-ROTATIONAL DYNAMIC FORMULATION FOR 2D BEAMS
COMPDYN 011 ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 5-8 May 011 CO-ROTATIONAL
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.
GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system
More informationLecture 15 Strain and stress in beams
Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME
More informationReview of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis
uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationShape Optimization of Revolute Single Link Flexible Robotic Manipulator for Vibration Suppression
15 th National Conference on Machines and Mechanisms NaCoMM011-157 Shape Optimization of Revolute Single Link Flexible Robotic Manipulator for Vibration Suppression Sachindra Mahto Abstract In this work,
More informationFinite element analysis of rotating structures
Finite element analysis of rotating structures Dr. Louis Komzsik Chief Numerical Analyst Siemens PLM Software Why do rotor dynamics with FEM? Very complex structures with millions of degrees of freedom
More informationPart D: Frames and Plates
Part D: Frames and Plates Plane Frames and Thin Plates A Beam with General Boundary Conditions The Stiffness Method Thin Plates Initial Imperfections The Ritz and Finite Element Approaches A Beam with
More informationMechanics of Inflatable Fabric Beams
Copyright c 2008 ICCES ICCES, vol.5, no.2, pp.93-98 Mechanics of Inflatable Fabric Beams C. Wielgosz 1,J.C.Thomas 1,A.LeVan 1 Summary In this paper we present a summary of the behaviour of inflatable fabric
More informationPROJECT 1 DYNAMICS OF MACHINES 41514
PROJECT DYNAMICS OF MACHINES 454 Theoretical and Experimental Modal Analysis and Validation of Mathematical Models in Multibody Dynamics Ilmar Ferreira Santos, Professor Dr.-Ing., Dr.Techn., Livre-Docente
More informationContents. Dynamics and control of mechanical systems. Focus on
Dynamics and control of mechanical systems Date Day 1 (01/08) Day 2 (03/08) Day 3 (05/08) Day 4 (07/08) Day 5 (09/08) Day 6 (11/08) Content Review of the basics of mechanics. Kinematics of rigid bodies
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationStructural Dynamics A Graduate Course in Aerospace Engineering
Structural Dynamics A Graduate Course in Aerospace Engineering By: H. Ahmadian ahmadian@iust.ac.ir The Science and Art of Structural Dynamics What do all the followings have in common? > A sport-utility
More informationFREE VIBRATIONS OF UNIFORM TIMOSHENKO BEAMS ON PASTERNAK FOUNDATION USING COUPLED DISPLACEMENT FIELD METHOD
A R C H I V E O F M E C H A N I C A L E N G I N E E R I N G VOL. LXIV 17 Number 3 DOI: 1.1515/meceng-17- Key words: free vibrations, Coupled Displacement Field method, uniform Timoshenko beam, Pasternak
More informationDYNAMICS OF MACHINERY 41514
DYNAMICS OF MACHINERY 454 PROJECT : Theoretical and Experimental Modal Analysis and Validation of Mathematical Models in Multibody Dynamics Holistic Overview of the Project Steps & Their Conceptual Links
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.
D : SOLID MECHANICS Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. Q.2 Consider the forces of magnitude F acting on the sides of the regular hexagon having
More informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More informationUsing the finite element method of structural analysis, determine displacements at nodes 1 and 2.
Question 1 A pin-jointed plane frame, shown in Figure Q1, is fixed to rigid supports at nodes and 4 to prevent their nodal displacements. The frame is loaded at nodes 1 and by a horizontal and a vertical
More informationModal Analysis: What it is and is not Gerrit Visser
Modal Analysis: What it is and is not Gerrit Visser What is a Modal Analysis? What answers do we get out of it? How is it useful? What does it not tell us? In this article, we ll discuss where a modal
More informationControl of constrained spatial three-link flexible manipulators
Control of constrained spatial three-link flexible manipulators Sinan Kilicaslan, M. Kemal Ozgoren and S. Kemal Ider Gazi University/Mechanical Engineering Department, Ankara, Turkey Middle East Technical
More informationDynamic and buckling analysis of FRP portal frames using a locking-free finite element
Fourth International Conference on FRP Composites in Civil Engineering (CICE8) 22-24July 8, Zurich, Switzerland Dynamic and buckling analysis of FRP portal frames using a locking-free finite element F.
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationMulti Linear Elastic and Plastic Link in SAP2000
26/01/2016 Marco Donà Multi Linear Elastic and Plastic Link in SAP2000 1 General principles Link object connects two joints, i and j, separated by length L, such that specialized structural behaviour may
More informationChapter 11. Displacement Method of Analysis Slope Deflection Method
Chapter 11 Displacement ethod of Analysis Slope Deflection ethod Displacement ethod of Analysis Two main methods of analyzing indeterminate structure Force method The method of consistent deformations
More informationConsider an elastic spring as shown in the Fig.2.4. When the spring is slowly
.3 Strain Energy Consider an elastic spring as shown in the Fig..4. When the spring is slowly pulled, it deflects by a small amount u 1. When the load is removed from the spring, it goes back to the original
More informationName: Fall 2014 CLOSED BOOK
Name: Fall 2014 1. Rod AB with weight W = 40 lb is pinned at A to a vertical axle which rotates with constant angular velocity ω =15 rad/s. The rod position is maintained by a horizontal wire BC. Determine
More information2. Rigid bar ABC supports a weight of W = 50 kn. Bar ABC is pinned at A and supported at B by rod (1). What is the axial force in rod (1)?
IDE 110 S08 Test 1 Name: 1. Determine the internal axial forces in segments (1), (2) and (3). (a) N 1 = kn (b) N 2 = kn (c) N 3 = kn 2. Rigid bar ABC supports a weight of W = 50 kn. Bar ABC is pinned at
More informationIntroduction to Mechanical Vibration
2103433 Introduction to Mechanical Vibration Nopdanai Ajavakom (NAV) 1 Course Topics Introduction to Vibration What is vibration? Basic concepts of vibration Modeling Linearization Single-Degree-of-Freedom
More informationMechanics of Materials II. Chapter III. A review of the fundamental formulation of stress, strain, and deflection
Mechanics of Materials II Chapter III A review of the fundamental formulation of stress, strain, and deflection Outline Introduction Assumtions and limitations Axial loading Torsion of circular shafts
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More informationLaboratory 4 Topic: Buckling
Laboratory 4 Topic: Buckling Objectives: To record the load-deflection response of a clamped-clamped column. To identify, from the recorded response, the collapse load of the column. Introduction: Buckling
More informationDynamics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Dynamics Semester 1, / 18
Dynamics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Dynamics Semester 1, 2016-17 1 / 18 Dynamics Dynamics studies the relations between the 3D space generalized forces
More informationTable of Contents. Preface... 13
Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...
More informationA consistent dynamic finite element formulation for a pipe using Euler parameters
111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,
More informationTransactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
An integral equation formulation of the coupled vibrations of uniform Timoshenko beams Masa. Tanaka & A. N. Bercin Department of Mechanical Systems Engineering, Shinshu University 500 Wakasato, Nagano
More informationVibration of a Euler Bernoulli uniform beam carrying a rigid body at each end
Vibration of a Euler Bernoulli uniform beam carrying a rigid body at each end S. Naguleswaran Department of Mechanical Engineering, University of Canterbury, Christchurch, New Zealand E-Mail: sivapatham.naguleswaran@canterbury.ac.nz
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationMethods of Analysis. Force or Flexibility Method
INTRODUCTION: The structural analysis is a mathematical process by which the response of a structure to specified loads is determined. This response is measured by determining the internal forces or stresses
More informationChapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14
Table of Contents Chapter 1: Research Objectives and Literature Review..1 1.1 Introduction...1 1.2 Literature Review......3 1.2.1 Describing Vibration......3 1.2.2 Vibration Isolation.....6 1.2.2.1 Overview.
More informationEQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION
1 EQUIVALENT SINGLE-DEGREE-OF-FREEDOM SYSTEM AND FREE VIBRATION The course on Mechanical Vibration is an important part of the Mechanical Engineering undergraduate curriculum. It is necessary for the development
More informationChapter 8 Rotational Motion
Chapter 8 Rotational Motion Chapter 8 Rotational Motion In this chapter you will: Learn how to describe and measure rotational motion. Learn how torque changes rotational velocity. Explore factors that
More informationChapter 5 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationCOPYRIGHTED MATERIAL. Index
Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,
More informationResearch Article On Equations of Motion of Elastic Linkages by FEM
Mathematical Problems in Engineering Volume 13, Article ID 6876, 1 pages http://dx.doi.org/1.1155/13/6876 Research Article On Equations of Motion of Elastic Linkages by FEM MichaB HaT Institute of Machine
More informationMembers Subjected to Torsional Loads
Members Subjected to Torsional Loads Torsion of circular shafts Definition of Torsion: Consider a shaft rigidly clamped at one end and twisted at the other end by a torque T = F.d applied in a plane perpendicular
More informationMODAL DERIVATIVES BASED REDUCTION METHOD FOR FINITE DEFLECTIONS IN FLOATING FRAME
Modal derivatives based reduction method for finite deflections in floating frame 11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th
More informationPLAT DAN CANGKANG (TKS 4219)
PLAT DAN CANGKANG (TKS 4219) SESI I: PLATES Dr.Eng. Achfas Zacoeb Dept. of Civil Engineering Brawijaya University INTRODUCTION Plates are straight, plane, two-dimensional structural components of which
More informationLecture 4: PRELIMINARY CONCEPTS OF STRUCTURAL ANALYSIS. Introduction
Introduction In this class we will focus on the structural analysis of framed structures. We will learn about the flexibility method first, and then learn how to use the primary analytical tools associated
More information41514 Dynamics of Machinery
41514 Dynamics of Machinery Theory, Experiment, Phenomenology and Industrial Applications Ilmar Ferreira Santos 1. Recapitulation Mathematical Modeling & Steps 2. Example System of Particle 3. Example
More informationHandout 7: Torque, angular momentum, rotational kinetic energy and rolling motion. Torque and angular momentum
Handout 7: Torque, angular momentum, rotational kinetic energy and rolling motion Torque and angular momentum In Figure, in order to turn a rod about a fixed hinge at one end, a force F is applied at a
More informationDynamic Analysis of Pelton Turbine and Assembly
Dynamic Analysis of Pelton Turbine and Assembly Aman Rajak, Prateek Shrestha, Manoj Rijal, Bishal Pudasaini, Mahesh Chandra Luintel Department of Mechanical Engineering, Central Campus, Pulchowk, Institute
More informationAdvanced Vibrations. Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Exact Solutions (Lecture 10) By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Exact Solutions Relation between Discrete and Distributed
More informationME 563 Mechanical Vibrations Lecture #1. Derivation of equations of motion (Newton-Euler Laws)
ME 563 Mechanical Vibrations Lecture #1 Derivation of equations of motion (Newton-Euler Laws) Derivation of Equation of Motion 1 Define the vibrations of interest - Degrees of freedom (translational, rotational,
More informationarxiv: v2 [physics.class-ph] 15 Dec 2016
Beyond the Kinetic Chain Process in a Stroke using a Triple Pendulum Model. Sun-Hyun Youn Department of Physics, Chonnam National University, Gwangju 61186, Korea arxiv:1612.02110v2 [physics.class-ph]
More informationModeling of Resonators
. 23 Modeling of Resonators 23 1 Chapter 23: MODELING OF RESONATORS 23 2 23.1 A GENERIC RESONATOR A second example where simplified discrete modeling has been found valuable is in the assessment of the
More informationSECTION A. 8 kn/m. C 3 m 3m
SECTION Question 1 150 m 40 kn 5 kn 8 kn/m C 3 m 3m D 50 ll dimensions in mm 15 15 Figure Q1(a) Figure Q1(b) The horizontal beam CD shown in Figure Q1(a) has a uniform cross-section as shown in Figure
More informationChapter 4 Deflection and Stiffness
Chapter 4 Deflection and Stiffness Asst. Prof. Dr. Supakit Rooppakhun Chapter Outline Deflection and Stiffness 4-1 Spring Rates 4-2 Tension, Compression, and Torsion 4-3 Deflection Due to Bending 4-4 Beam
More informationScienceDirect. The Stability of a Precessing and Nutating Viscoelastic Beam with a Tip Mass
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 144 (2016 ) 68 76 12th International Conference on Vibration Problems, ICOVP 2015 The Stability of a Precessing and Nutating
More informationCHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS
61 CHAPTER 4 DESIGN AND ANALYSIS OF CANTILEVER BEAM ELECTROSTATIC ACTUATORS 4.1 INTRODUCTION The analysis of cantilever beams of small dimensions taking into the effect of fringing fields is studied and
More informationMechanics of Materials
Mechanics of Materials 2. Introduction Dr. Rami Zakaria References: 1. Engineering Mechanics: Statics, R.C. Hibbeler, 12 th ed, Pearson 2. Mechanics of Materials: R.C. Hibbeler, 9 th ed, Pearson 3. Mechanics
More informationResponse Spectrum Analysis Shock and Seismic. FEMAP & NX Nastran
Response Spectrum Analysis Shock and Seismic FEMAP & NX Nastran Table of Contents 1. INTRODUCTION... 3 2. THE ACCELEROGRAM... 4 3. CREATING A RESPONSE SPECTRUM... 5 4. NX NASTRAN METHOD... 8 5. RESPONSE
More informationAssumptions: beam is initially straight, is elastically deformed by the loads, such that the slope and deflection of the elastic curve are
*12.4 SLOPE & DISPLACEMENT BY THE MOMENT-AREA METHOD Assumptions: beam is initially straight, is elastically deformed by the loads, such that the slope and deflection of the elastic curve are very small,
More informationModelling and Finite Element Analysis of Double Wishbone Suspension
Modelling and Finite Element Analysis of Double Wishbone Suspension Amol Patil, Varsha Patil, Prashant Uhle P.G. Student, Dept. of Mechanical Engineering, S S B T S College of Engineering, Jalgaon, Maharastra,
More informationUNIT IV FLEXIBILTY AND STIFFNESS METHOD
SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) Subject with Code : SA-II (13A01505) Year & Sem: III-B.Tech & I-Sem Course & Branch: B.Tech
More informationIntroduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.
Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.
More informationMore Examples Of Generalized Coordinates
Slides of ecture 8 Today s Class: Review Of Homework From ecture 7 Hamilton s Principle More Examples Of Generalized Coordinates Calculating Generalized Forces Via Virtual Work /3/98 /home/djsegal/unm/vibcourse/slides/ecture8.frm
More informationE 490 FE Exam Prep. Engineering Mechanics
E 490 FE Exam Prep Engineering Mechanics 2008 E 490 Course Topics Statics Newton s Laws of Motion Resultant Force Systems Moment of Forces and Couples Equilibrium Pulley Systems Trusses Centroid of an
More information0000. Finite element modeling of a wind turbine blade
ArticleID: 16033; Draft date: 2015-07-27 0000. Finite element modeling of a wind turbine blade Mohammad Sheibani 1, Ali Akbar Akbari 2 Department of Mechanical Engineering, Ferdowsi University of Mashhad,
More informationPLANAR KINETIC EQUATIONS OF MOTION (Section 17.2)
PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2) We will limit our study of planar kinetics to rigid bodies that are symmetric with respect to a fixed reference plane. As discussed in Chapter 16, when
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationEngineering Science OUTCOME 2 - TUTORIAL 3 FREE VIBRATIONS
Unit 2: Unit code: QCF Level: 4 Credit value: 5 Engineering Science L/60/404 OUTCOME 2 - TUTORIAL 3 FREE VIBRATIONS UNIT CONTENT OUTCOME 2 Be able to determine the behavioural characteristics of elements
More informationPreliminaries: Beam Deflections Virtual Work
Preliminaries: Beam eflections Virtual Work There are several methods available to calculate deformations (displacements and rotations) in beams. They include: Formulating moment equations and then integrating
More informationPLANAR KINETICS OF A RIGID BODY: WORK AND ENERGY Today s Objectives: Students will be able to: 1. Define the various ways a force and couple do work.
PLANAR KINETICS OF A RIGID BODY: WORK AND ENERGY Today s Objectives: Students will be able to: 1. Define the various ways a force and couple do work. In-Class Activities: 2. Apply the principle of work
More informationPROJECT 2 DYNAMICS OF MACHINES 41514
PROJECT 2 DYNAMICS OF MACHINES 41514 Dynamics of Rotor-Bearing System Lateral Vibrations and Stability Threshold of Rotors Supported On Hydrodynamic Bearing and Ball Bearing. Ilmar Ferreira Santos, Prof.
More informationFinal Exam April 30, 2013
Final Exam Instructions: You have 120 minutes to complete this exam. This is a closed-book, closed-notes exam. You are allowed to use a calculator during the exam. Usage of mobile phones and other electronic
More informationCHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES
CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES 14.1 GENERAL REMARKS In structures where dominant loading is usually static, the most common cause of the collapse is a buckling failure. Buckling may
More informationRotation. I. Kinematics - Angular analogs
Rotation I. Kinematics - Angular analogs II. III. IV. Dynamics - Torque and Rotational Inertia Work and Energy Angular Momentum - Bodies and particles V. Elliptical Orbits The student will be able to:
More informationM.S Comprehensive Examination Analysis
UNIVERSITY OF CALIFORNIA, BERKELEY Spring Semester 2014 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials Name:......................................... M.S Comprehensive
More informationDeflections and Strains in Cracked Shafts due to Rotating Loads: A Numerical and Experimental Analysis
Rotating Machinery, 10(4): 283 291, 2004 Copyright c Taylor & Francis Inc. ISSN: 1023-621X print / 1542-3034 online DOI: 10.1080/10236210490447728 Deflections and Strains in Cracked Shafts due to Rotating
More informationDynamics. describe the relationship between the joint actuator torques and the motion of the structure important role for
Dynamics describe the relationship between the joint actuator torques and the motion of the structure important role for simulation of motion (test control strategies) analysis of manipulator structures
More informationGyroscopic matrixes of the straight beams and the discs
Titre : Matrice gyroscopique des poutres droites et des di[...] Date : 29/05/2013 Page : 1/12 Gyroscopic matrixes of the straight beams and the discs Summarized: This document presents the formulation
More informationChapter 7. ELASTIC INSTABILITY Dr Rendy Thamrin; Zalipah Jamellodin
Chapter 7 ESTIC INSTIITY Dr Rendy Thamrin; Zalipah Jamellodin 7. INTRODUCTION TO ESTIC INSTIITY OF COUN ND FRE In structural analysis problem, the aim is to determine a configuration of loaded system,
More informationME C85/CE C30 Fall, Introduction to Solid Mechanics ME C85/CE C30. Final Exam. Fall, 2013
Introduction to Solid Mechanics ME C85/CE C30 Fall, 2013 1. Leave an empty seat between you and the person (people) next to you. Unfortunately, there have been reports of cheating on the midterms, so we
More informationRigid bodies - general theory
Rigid bodies - general theory Kinetic Energy: based on FW-26 Consider a system on N particles with all their relative separations fixed: it has 3 translational and 3 rotational degrees of freedom. Motion
More informationIf the number of unknown reaction components are equal to the number of equations, the structure is known as statically determinate.
1 of 6 EQUILIBRIUM OF A RIGID BODY AND ANALYSIS OF ETRUCTURAS II 9.1 reactions in supports and joints of a two-dimensional structure and statically indeterminate reactions: Statically indeterminate structures
More informationFree vibration analysis of beams by using a third-order shear deformation theory
Sādhanā Vol. 32, Part 3, June 2007, pp. 167 179. Printed in India Free vibration analysis of beams by using a third-order shear deformation theory MESUT ŞİMŞEK and TURGUT KOCTÜRK Department of Civil Engineering,
More informationGeneralized Single Degree of Freedom Systems
Single Degree of Freedom http://intranet.dica.polimi.it/people/boffi-giacomo Dipartimento di Ingegneria Civile Ambientale e Territoriale Politecnico di Milano March, 8 Outline Until now our were described
More informationArtificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik. Robot Dynamics. Dr.-Ing. John Nassour J.
Artificial Intelligence & Neuro Cognitive Systems Fakultät für Informatik Robot Dynamics Dr.-Ing. John Nassour 25.1.218 J.Nassour 1 Introduction Dynamics concerns the motion of bodies Includes Kinematics
More informationAvailable online at ScienceDirect. The 2014 conference of the International Sports Engineering Association
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 72 ( 2014 ) 593 598 The 2014 conference of the International Sports Engineering Association Effect of moment of inertia and
More informationThree-Dimensional Biomechanical Analysis of Human Movement
Three-Dimensional Biomechanical Analysis of Human Movement Anthropometric Measurements Motion Data Acquisition Force Platform Body Mass & Height Biomechanical Model Moments of Inertia and Locations of
More informationMechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002
student personal identification (ID) number on each sheet. Do not write your name on any sheet. #1. A homogeneous, isotropic, linear elastic bar has rectangular cross sectional area A, modulus of elasticity
More informationMarch 24, Chapter 4. Deflection and Stiffness. Dr. Mohammad Suliman Abuhaiba, PE
Chapter 4 Deflection and Stiffness 1 2 Chapter Outline Spring Rates Tension, Compression, and Torsion Deflection Due to Bending Beam Deflection Methods Beam Deflections by Superposition Strain Energy Castigliano
More informationProceedings The Behavior of Badminton Shuttlecocks from an Engineering Point of View
Proceedings The Behavior of Badminton Shuttlecocks from an Engineering Point of View Christoffer Johansson 1, Kelken Chang 2, Christer Forsgren 3 and Magnus Karlsteen 1, * 1 Department of Physics, Chalmers
More informationSAMCEF For ROTORS. Chapter 1 : Physical Aspects of rotor dynamics. This document is the property of SAMTECH S.A. MEF A, Page 1
SAMCEF For ROTORS Chapter 1 : Physical Aspects of rotor dynamics This document is the property of SAMTECH S.A. MEF 101-01-A, Page 1 Table of Contents rotor dynamics Introduction Rotating parts Gyroscopic
More informationNatural vibration frequency of classic MEMS structures
Natural vibration frequency of classic MEMS structures Zacarias E. Fabrim PGCIMAT, UFRGS, Porto Alegre, RS, Brazil Wang Chong, Manoel Martín Pérez Reimbold DeTec, UNIJUI, Ijuí, RS, Brazil Abstract This
More informationLecture 11: The Stiffness Method. Introduction
Introduction Although the mathematical formulation of the flexibility and stiffness methods are similar, the physical concepts involved are different. We found that in the flexibility method, the unknowns
More information