DESIGN ANALYSIS OF A RUBBER MOUNT SYSTEM FOR A PUSH- TYPE CENTRIFUGE Carl Nel

Similar documents
THE INFLUENCE OF AMPLITUDE AND PHASE DIFFERENCES IN BI-DIRECTIONAL GROUND MOTION ON THE BEHAVIOUR OF IRREGULAR STRUCTURES.

Start with the equation of motion for a linear multi-degree of freedom system with base ground excitation:

Rigid body simulation

Application to Plane (rigid) frame structure

Modeling of Dynamic Systems

Other Uncertainty Analysis Methods

SCHOOL OF COMPUTING, ENGINEERING AND MATHEMATICS SEMESTER 2 EXAMINATIONS 2011/2012 DYNAMICS ME247 DR. N.D.D. MICHÉ

5.76 Lecture #21 2/28/94 Page 1. Lecture #21: Rotation of Polyatomic Molecules I

One Dimensional Axial Deformations

Modeling and Simulation of a Hexapod Machine Tool for the Dynamic Stability Analysis of Milling Processes. C. Henninger, P.

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

PART 8. Partial Differential Equations PDEs

SCALARS AND VECTORS All physical quantities in engineering mechanics are measured using either scalars or vectors.

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

AMSC 660 Scientific Computing I. Term Project The Assignment

Numerical Methods. ME Mechanical Lab I. Mechanical Engineering ME Lab I

Dynamics 4600:203 Homework 08 Due: March 28, Solution: We identify the displacements of the blocks A and B with the coordinates x and y,

Kinematics in 2-Dimensions. Projectile Motion

Structural Dynamics and Earthquake Engineering

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Study on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component

S Advanced Digital Communication (4 cr) Targets today

ENGN 40 Dynamics and Vibrations Homework # 7 Due: Friday, April 15

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

829. An adaptive method for inertia force identification in cantilever under moving mass

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

Lesson 4: Relative motion, Forces, Newton s laws (sections )

MEEM 3700 Mechanical Vibrations

Mathematical Preparations

PLATE BENDING ELEMENTS

Numerical analysis of buried pipes subjected to lateral soil movements

EN40: Dynamics and Vibrations. Homework 7: Rigid Body Kinematics

Polynomial Regression Models

Numerical Heat and Mass Transfer

Elshaboury SM et al.; Sch. J. Phys. Math. Stat., 2015; Vol-2; Issue-2B (Mar-May); pp

Instituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

Multidisciplinary Design Optimization of Elastomeric Mounting Systems in Automotive Vehicles

where v means the change in velocity, and t is the

CHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)

Statistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

Optimal Design of Monorail Crane Spring Based on MATLAB Quanwei Wang, Wenan Kang, Shan Nie

MODAL ANALYSIS AND TESTING OF ROTATING MACHINES FOR PREDICTIVE MAINTENANCE: EFFECT OF GYROSCOPIC FORCES

So far: simple (planar) geometries

Frequency-weighted MPC of Trailing Edge Flaps. Tests on a V27.

CHAPTER 4. Vector Spaces

Week 9 Chapter 10 Section 1-5

Week 11: Chapter 11. The Vector Product. The Vector Product Defined. The Vector Product and Torque. More About the Vector Product

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Torsten Mayer-Gürr Institute of Geodesy, NAWI Graz Technische Universität Graz

Grid Generation around a Cylinder by Complex Potential Functions

9. Complex Numbers. 1. Numbers revisited. 2. Imaginary number i: General form of complex numbers. 3. Manipulation of complex numbers

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

), it produces a response (output function g (x)

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 141. Lecture 14. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 14, Page 1

Investigation on the Wheel/Rail Contact and Longitudinal Train/Track Interaction Forces

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

Notes on Analytical Dynamics

Research on the Fuzzy Control for Vehicle Semi-active Suspension. Xiaoming Hu 1, a, Wanli Li 1,b

Chapter 3. Estimation of Earthquake Load Effects

Indeterminate pin-jointed frames (trusses)

MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava

coordinates. Then, the position vectors are described by

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES

Kernel Methods and SVMs Extension

Lesson 5: Kinematics and Dynamics of Particles

e i is a random error

Stereo 3D Simulation of Rigid Body Inertia Ellipsoid for The Purpose of Unmanned Helicopter Autopilot Tuning

Identification of Instantaneous Modal Parameters of A Nonlinear Structure Via Amplitude-Dependent ARX Model

CHARACTERISATION OF VIBRATION ISOLATORS USING VIBRATION TEST DATA

Lagrange Multipliers. A Somewhat Silly Example. Monday, 25 September 2013

Second Order Analysis

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES

Conservation of Angular Momentum = "Spin"

LECTURE 21 Mohr s Method for Calculation of General Displacements. 1 The Reciprocal Theorem

PHYS 705: Classical Mechanics. Calculus of Variations II

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Evaluation of the accuracy of the Multiple Support Response Spectrum (MSRS) method

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Pulse Coded Modulation

Measurement and Compensation for Volumetric Positioning Errors of CNC Machine Tools Considering Thermal Effect

A Tutorial on Data Reduction. Linear Discriminant Analysis (LDA) Shireen Elhabian and Aly A. Farag. University of Louisville, CVIP Lab September 2009

C PLANE ELASTICITY PROBLEM FORMULATIONS

A particle in a state of uniform motion remain in that state of motion unless acted upon by external force.

WAVE PROPAGATION, REFLECTION AND TRANSMISSION IN CURVED BEAMS

Part C Dynamics and Statics of Rigid Body. Chapter 5 Rotation of a Rigid Body About a Fixed Axis

International Journal of Mathematics Trends and Technology (IJMTT) Volume 48 Number 2 August 2017

CHAPTER 14 GENERAL PERTURBATION THEORY

11. Dynamics in Rotating Frames of Reference

Canonical transformations

Simple Method for Pushover Curves of Asymmetric Structure with Displacement-Dependent Passive Energy Dissipation Devices

The Minimum Universal Cost Flow in an Infeasible Flow Network

Spin-rotation coupling of the angularly accelerated rigid body

Transcription:

ICSV4 Carns Australa 9-2 Jul, 27 DESIGN ANALYSIS OF A RUBBER MOUNT SYSTEM FOR A PUSH- TYPE CENTRIFUGE Carl Nel School of Mechancal Enneern, North-West Unverst Potchefstroom, South Afrca carlbnel@netactve.co.a Abstract The lack of balance, whch occurs n all rotatn machnes, could cause troublesome vbraton f the attachments to the foundaton of the machne are unsutable. Ths s partcularl true for centrfues, where uneven loadn of the basket, for eample, due to temporar rreulartes n the feed or laers of salt sedment on the basket, ma lead to partcularl bad balancn. The am of ths nvestaton was to evaluate the use of rubber mounts and to predct feasble mount stffness coeffcents to be used. These mounts have to produce ood vbraton solaton to protect the supportn buldn from dnamc forces transmtted, but must also lmt the statc and dnamc mount dsplacements to allow acceptable moton of the centrfue assembl. It was therefore necessar to use a s-deree of freedom mathematcal model, mplemented n computer prorams to compute these statc and dnamc mount dsplacements and mount forces, for 2 possble desns. An objectve functon as a measure of vbraton transmtted, was also used. In addton, the mount sstem natural frequences were also computed, to prevent resonance wth the normal basket anular speed and also the forced frequenc of the aal pusher mechansm.. INTRODUCTION Ths nvestaton reardn vbraton solaton was done for centrfues located at the frst level of supportn renforced concrete buldn. A s-deree of freedom model was used. A drect approach n terms of forces transmtted to the support structure s consdered. 2. MATHEMATICAL MODEL A mount sstem model smlar to that emploed b Nel and Hens s consdered [3], [5], [6]. Hence the centrfue assembl (centrfue and base plate) s dealsed as a rd bod of mass m attached to a rd support structure b means of an arbtrar number n of elastc mounts

wth arbtrar postons and orentatons wth respect to the centrfue assembl lobal coordnate sstem. The orn of the fed lobal co-ordnate sstem s located at, the center of ravt (c..) of the centrfue assembl as shown n fure. It s assumed that the stffness coeffcent k, k, k n the three co-ordnate drectons are ndependent of each other f the rotatonal stffness of the mount s nelected. For elastomerc materals, as are usuall used at rubber mounts, a hsteretc dampn model wth a comple stffness matr s assumed here [2], [4], [5], see also fure, so that j k j k j k + + + ) ( ) ( ) ( η η η () where j =, and η, η and η are the mount loss factors n each of the three co-ordnate drectons. A transformaton s requred to relate the translatonal dsplacements of each mount to translatonal and rotatonal dsplacements of the centrfue assembl c.. Assumn small centrfue dsplacements, G s defned b u = = = = G U (2) where u s the translatonal dsplacement matr at mountn pont and U are the centrfue assembl c.. dsplacements.,, are translatonal and,, rotatonal dsplacements. Wth F vector comprsn the three forces and the three moments due to mount actn on the centrfue assembl alon, and and about, and respectvel, t now follows that F = - [ ] T G K G U (3) Assumn aan small dsplacements, the rd bod equatons of moton [] s M Ü = F (4) wth F the matr of resultant total forces and moments on the centrfue assembl. The sstem mass matr s M = I I I m m m (5)

where m s the total centrfue assembl mass and I, I and I are the elements of centrfue assembl nerta epressed n terms of the lobal co-ordnate sstem. The centrfue assembl c.. acceleraton matr s Ü = && && && && && && (6) where & &, & &, & & are translatonal acceleratons and &, &, & anular acceleratons. Addn toether the effects of all the mounts of the centrfue assembl, t follows that n F = - n T G K G U (7) = = = - K U where K s s a comple dnamc matr whch ncludes hsteretc loss effects at the centrfue c.. From ths and equaton (4) t now follows that M U & + K U = F e (8) where F e represents all forces and moments on the centrfue assembl other than the mount reacton forces. In order to take advantae of the frequenc doman approach, F e s assumed to comprse of v snusodal forces and moments ncluded b F wth correspondn frequences ω h and phase anles a h h,2,..., v, =. It ma then be shown that dh v U = h= 2 [ K ] ω h M [ d ] h F (9) where U and F dh represent dsplacement and force phasors respectvel [2]. The centrfues assembl shakn forces and moments are F F e = e F F M M M e e F F = M M M e e e e e e e e () The model ma be emploed to determne an objectve functon. The dnamc dsplacements at each mount n the, and drectons are computed n the tme doman. It s thus possble to compute for each force component at each mount the averae value of the

sum of the dnamc forces at each tme ncrement dt over the perod T. As crteron of vbraton transmtted, an objectve functon s defned b T n T 2 2 2 ϕ = O ( f ( X ) + f ( X ) + f ( X )) dt = () n whch X s a vector of desn varables chosen here to comprse stffness coeffcents. It s now requred to mnme wth respect to varables X, subject to nequalt constrants w = U ma( U ) U (2) s + d c where U c s a vector of specfed mamum acceptable mount translatonal dsplacement ampltudes. The mamum dnamc translatonal dsplacement elements n U d are n the same drecton as the statc dsplacement elements n U s. These constrants also ve an ndcaton of acceptable centrfue assembl moton and acceptable mount dsplacements. 3. CHARACTERISATION OF PARAMETERS The manufacturer suppled the centrfue and base plate mass and centre of mass propertes. These propertes were used n the computaton of the moments of nerta for the total centrfue and base plate assembl for each lobal as. The total prncpal moments of nerta are 5272.7, 6395. and 7468.5 km 2 correspondn to the centrfue lobal, and aes as shown n fure. The total centrfue assembl mass s 98 k. The lobal coordnates for each of the mounts and those for the basket centre were determned. The centrfue assembl s supported b 4 mounts postoned at the follown lobal co-ordnates: Table 4. Mount lobal co-ordnates [m]. Mount.55.8 -.795 2 -.496.8 -.795 3 -.496 -.8 -.795 4.55 -.8 -.795 The lobal co-ordnates of the basket centre are: Table 5. Basket centre lobal co-ordnates [m]. c c c -.925..42 The basket radus s 4 mm and the reatest possble ecess mass actn on ths radus s 2.3 k (provded b the manufacturer). The normal runnn speed of the basket s 2 rev/mn. The stffness coeffcents of the 2 estn ver soft ppe compensators used were nelected, because these mantudes were rearded as ver small compared to the stffness

coeffcents of tpcal rubber mounts. The dnamc mount stffness coeffcents of standard mounts avalable were chosen for 2 dfferent possble desns A and B as shown n Table 6. Table 6. Mount stffness coeffcents [kn/m]. Mount stffness coeffcents [kn/m] Desn A Desn B (Resonance) (Feasble) k, k 2, k 3, k 4 2 k, k 2, k 3, k 4 2 k, k 2, k 3, k 4 533 456 Snce mount loss factors n the rane of.5 to.5 do not apprecabl affect the dnamc behavour, and rubber mounts normall have mount loss factors n ths rane, mount loss factor values of. were assumed for all the mounts n all 3 orthoonal drectons. 4. COMPUTER IMPLEMENTATION The mathematcal model was mplemented n computer prorams n a MATLAB envronment. The mantudes of the characterstcs of the parameters descrbed above were used as nput data. The centrfue assembl shakn forces and moments were computed frst, and then used n the computaton of dnamc and statc mount dsplacements, dnamc and statc mount forces transmtted to the support structure throuh the mounts and also objectve functon values. 5. CENTRIFUGE ASSEMBLY SHAKING FORCES AND MOMENTS The mantudes of the parameters requred for computaton of the centrfue assembl shakn forces and moments and mathematcal equatons mplemented n computer prorams were also used, to compute these forces and moments. Fure 2 shows ths force and moment tme hstores for the centrfue anular speed at 2 rev/mn. The frequenc observed from these wave forms, corresponds to the centrfue basket anular speed. These tme doman forces and moments were Fourer analsed b usn the MATLAB fft.m functon, and the Fourer coeffcents at correspondn frequences and phase anles were then used n the computaton of the responses n the, and drectons at each mount b usn the frequenc doman computer proram (see for eample fure 3). Table 7 shows the force and moment ampltudes at phase anles. The force mantude of the aal pusher mechansm was consdered relatvel small compared to centrfual forces, and therefore nelected. The mantude of ths ver low forced frequenc force s mportant, and was taken nto account when the mount sstem natural frequences were evaluated (see Table 9), to avod resonance. Table 7. Shakn force [N] and moment [Nm] ampltudes at phase anles [rad] and anular speed of 25.46 rad/s. F e F e M e M e M e Ampltudes [N] and [Nm] 452.8 452.8 6.8 343.8 343.8 Phase anle [rad] -.578.578-3.46

6. RESULTS AND EVALUATION As crteron or measure of vbraton transmtted to the support structure, the objectve functon (equaton ) was used. Ths snle objectve functon value was then computed for each of the 2 dfferent desns A and B and also a rd mounted sstem desn, as shown n Table 8 for a basket speed at 2 rev/mn. Lare objectve functon values ndcate bad vbraton solaton, and lower objectve functon values mproved vbraton solaton [5], [6]. Table 8. Objectve functon [kn] at 2 rev/mn. Desn A Desn B Rd desn (Resonance) (Feasble) Objectve functon ψ [kn] 32.8 47. 57.3 B usn the MATLAB e.m functon, the eenvalues and eenvectors were also computed for the 2 dfferent desns A and B. Table 9 shows these natural frequences and mode shapes. The stffness coeffcents used for desns A and B are shown n Table 6. These natural frequences are mportant n order to avod resonance wth the basket normal anular speed (2 H), but also wth the aal pusher mechansm force frequenc ( to.25 H). Table 9. Natural frequences [H] and correspondn mode shapes. Natural frequences [H] and mode shapes Desn A 5.4 7.65 2.2 2.59 6.78 2. Desn B 2.45 2.87 3.52 5.83 6.88 8.39 7. CONCLUSIONS The use of the ver low stffness coeffcents renders the best vbraton solaton, but has the larest statc mount dsplacements. The use of ver low stffness coeffcents also results n low sstem natural frequences, whch could easl correspond to the pusher mechansm force frequenc ( to.25 H), thus resonance. The use of ver low stffness coeffcents also results n unacceptable centrfue assembl moton, thus lare mount dnamc dsplacements. Desn A wth the moderatel hh stffness coeffcents renders the poorest vbraton solaton (see Table 8 for comparson of objectve functon mantudes), but also has smaller statc mount dsplacements compared to desn B, but results n the larest dnamc mount dsplacements. Ths vbraton solaton s unacceptable. The use of hher stffness coeffcents n eneral tends to create natural frequences whch could easl correspond, or be near to the centrfue basket operaton speed, see also Table 9, whch s the case for desn A. When one of the mount sstem natural frequences corresponds to the forced basket frequenc, then resonance takes place, whch s undesrable and danerous and could lead to falures and mount durablt problems. The objectve functon used has the larest values when resonance takes place (see Tables 8 and 9), whch ndcates that ts use s effectve for the measure of vbraton transmtted.

Desn B renders the best compromse between vbraton solaton and acceptable mount dsplacements. The natural frequences of desn B are also far awa enouh from the basket s normal operatonal speed, and also the forced frequenc of the pusher mechansm (see Table 9). The vbraton solaton obtaned for desn B s also snfcantl better compared to a desn where mount stffness coeffcents were chosen to represent a centrfue whch s rdl attached (see Table 8 for comparson of objectve functon mantudes). REFERENCES [] A.F. D Soua and V.K. Gar, Advanced dnamcs Modeln and analss, Prentce-Hall, 984. [2] D.J. Ewns, Modal Testn: Theor and practce, Research Studes Press, 984. [3] P.S. Hens, C.B. Nel, and J.A. Snman, Optmsaton of enne mountn confuraton. Proceedns of ISMA 9. Tools for Nose and Vbraton Analss, Volume II. Katholeke Unverstet Leuven, Belum, September 994. [4] C.B. Nel, and P.S. Hens, An optmsaton approach to mountn charactersaton, Proceedns of Nose and Vbraton 95, Unverst of Pretora, November 995. [5] C.B. Nel and P.S. Hens, Epermental verfcaton of an optmsaton proram for a front-wheel-drve enne mount sstem, Proceedns of ISMA2, Nose and Vbraton Enneern, Voume. III, Katholeke Unverstet, Leuven, Belum, September 996. [6] C.B. Nel, Optmsaton of enne mount sstems for multple operatonal condtons at front-wheel-drve vehcles, Proceedns of ISMA25, Internatonal Conference on Nose and Vbraton Enneern, Katholeke Unverstet, Leuven, Belum, September 2.