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

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

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Effect of loading frequency on the settlement of granular layer

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

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

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

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

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

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

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

In this section is given an overview of the common elasticity models.

Finite Element Modelling of truss/cable structures

Dynamic Analysis Based On ANSYS of Turning and Grinding Compound Machine Spindle Box Zanhui Shu and Qiushi Han

Operating conditions of a mine fan under conditions of variable resistance

Pulse Coded Modulation

Indeterminate pin-jointed frames (trusses)

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

Difference Equations

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Uncertainty in measurements of power and energy on power networks

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

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

Design and Analysis of Landing Gear Mechanic Structure for the Mine Rescue Carrier Robot

Module 3: Element Properties Lecture 1: Natural Coordinates

APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS

6.3.7 Example with Runga Kutta 4 th order method

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

Uncertainty and auto-correlation in. Measurement

A Fast Computer Aided Design Method for Filters

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

Linear Feature Engineering 11

De-noising Method Based on Kernel Adaptive Filtering for Telemetry Vibration Signal of the Vehicle Test Kejun ZENG

COEFFICIENT DIAGRAM: A NOVEL TOOL IN POLYNOMIAL CONTROLLER DESIGN

MMA and GCMMA two methods for nonlinear optimization

Numerical Heat and Mass Transfer

DESIGN AND RESEARCH INTO THE NONLINEAR MAIN VIBRATION SPRING IN DOUBLE-MASS HIGH ENERGY VIBRATION MILLING

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

MODEL OF HYDROPNEUMATIC THREE POINT HITCH

The Study of Teaching-learning-based Optimization Algorithm

Over-Temperature protection for IGBT modules

on the improved Partial Least Squares regression

Color Rendering Uncertainty

CHARACTERISATION OF VIBRATION ISOLATORS USING VIBRATION TEST DATA

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Lecture 8 Modal Analysis

Some modelling aspects for the Matlab implementation of MMA

SIO 224. m(r) =(ρ(r),k s (r),µ(r))

Evaluation of Validation Metrics. O. Polach Final Meeting Frankfurt am Main, September 27, 2013

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Digital Signal Processing

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

Vibration analysis on the rolling element bearing-rotor system of an air blower

Parking Demand Forecasting in Airport Ground Transportation System: Case Study in Hongqiao Airport

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Chapter 13: Multiple Regression

Lifetime prediction of EP and NBR rubber seal by thermos-viscoelastic model

The Tangential Force Distribution on Inner Cylinder of Power Law Fluid Flowing in Eccentric Annuli with the Inner Cylinder Reciprocating Axially

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Tensor Smooth Length for SPH Modelling of High Speed Impact

Second Order Analysis

A New Refinement of Jacobi Method for Solution of Linear System Equations AX=b

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

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

Description of the Force Method Procedure. Indeterminate Analysis Force Method 1. Force Method con t. Force Method con t

Introduction to Regression

NON LINEAR ANALYSIS OF STRUCTURES ACCORDING TO NEW EUROPEAN DESIGN CODE

EEE 241: Linear Systems

How Differential Equations Arise. Newton s Second Law of Motion

Valuated Binary Tree: A New Approach in Study of Integers

This column is a continuation of our previous column

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

Integrals and Invariants of Euler-Lagrange Equations

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

EN40: Dynamics and Vibrations. Homework 4: Work, Energy and Linear Momentum Due Friday March 1 st

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

The Analysis of Convection Experiment

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

Dynamic Similarity Design of Geared Rotor System in Five-Shaft Integrally Centrifugal Compressor

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

THE CURRENT BALANCE Physics 258/259

DUE: WEDS FEB 21ST 2018

Turbulence classification of load data by the frequency and severity of wind gusts. Oscar Moñux, DEWI GmbH Kevin Bleibler, DEWI GmbH

Frequency dependence of the permittivity

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

Kernel Methods and SVMs Extension

5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

Appendix B. The Finite Difference Scheme

Solution of Equilibrium Equation in Dynamic Analysis. Mode Superposition. Dominik Hauswirth Method of Finite Elements II Page 1

Handout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs.

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites

MEEM 3700 Mechanical Vibrations

Constitutive Modelling of Superplastic AA-5083

1 Derivation of Point-to-Plane Minimization

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

The Quadratic Trigonometric Bézier Curve with Single Shape Parameter

A Particle Filter Algorithm based on Mixing of Prior probability density and UKF as Generate Importance Function

The internal structure of natural numbers and one method for the definition of large prime numbers

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Transcription:

Study on Non-Lnear Dynamc Characterstc of Vehcle Suspenson Rubber Component Zhan Wenzhang Ln Y Sh GuobaoJln Unversty of TechnologyChangchun, Chna Wang Lgong (MDI, Chna [Abstract] The dynamc characterstc of rubber component s one of mportant factors that nfluences the NVH performances of vehcle In common case, smulaton model bult on ADAMS software consders only the lnear characterstc But n some specal cases, rubber components nstalled n the vehcle suspenson always work n non-lnear condton Amed at the research dffculty of ths knd of non-lnear dynamc characterstc, a numercal ft method had been put forward based on energy analyss theory n ths paper To acheve the nonlnear dynamc characterstc of rubber component n the vehcle dynamcs smulaton model based on the ADAMS, a user wrtten subroutne has been compled, and a custom solver was created, hence the non-lnear analyss of vehcle dynamcs can be acheved Key words: Vehcle Suspenson Rubber Component ADAMS Software Non-lnear Dynamc Characterstc Energy Analyss INTRODUCTION Wth economy growng and traffc envronment mprovng n Chna, the performance of modern vehcle s requred more and more strctly In order to reduce the transference of vbraton and nose, and to mprove the comfortablty of vehcle, the rubber component s used wdely n vehcle structure, especally n some poston where connect and jon the parts of suspenson For example, Aud1 car has more than 1 rubber components nstalled n front suspenson, t s to utlze the far characterstc of rubber components to reduce vbraton and nose Snce many rubber components were used n vehcle, ther dynamc characterstc affect vehcle performance, especally on the dynamc characterstc of suspenson The dynamc characterstc of rubber components s mostly obtaned by testng In common case, the dynamc characterstc of rubber components consdered as lnear But n some specal places, for example n the jont of low control arm and lateral stablty-bar of Aud1, the dynamc characterstc of rubber component s often workng n non-lnear range Based on the energy analyss theory ntroduced n reference [1], ths paper uses dynamc tral data to deal wth the nonlnear dynamc characterstc of rubber components, and puts forward the numercal fttng method based on energy analyss, and realzes non-lnearzaton of vehcle dynamc smulaton LINEAR DYNAMIC CHARACTERISTIC OF RUBBER COMPONENT The dynamc characterstc of rubber components can be classed by lnear and non-lnear characterstc [2], and t s always descrbed by the relatonshp of force and deformaton In the process of vehcle dynamc smulaton, the dynamc model of rubber components serously affects the analyss precson The dynamc characterstc of rubber components n the lnear range, can de obtaned by analyss the force-deformaton curve of rubber component whch got by actng perodcal exctng force In the lnear range, the force-deformaton relatonshp of rubber component can be descrbed wth followng equatons Exctng dsplacement of rubber component can be descrbe as: Xt ( = X cos( ωt 1 Where: X ampltude of exctng dsplacement; ω exctng frequency Correspondng deformaton force of rubber component s: Ft ( = Fcos( ωt+ δ 2 2 Internatonal ADAMS User Conference Page 1 of 5

Where: F ampltude of deformaton force; δ phase dfference between deformaton force and exctng dsplacement In fact, the phase of deformaton force precedes to that of exctng dsplacement, and the phase dfference δ s δ < 9 Expandng the equaton of deformaton force as: F( t = Fk( t + Fc( t 3 Where: F Fk ( t = Fcosδ cos( ωt = cosδxcos( ωt X 4 Fc ( t = F snδ sn( ωt The frst equaton of (4 denotes the sprng force caused by the dynamc stffness of rubber F component, then the dynamc stffness K s cosδ X The second equaton of (4 denotes the dampng force caused by the dampng of rubber component, whch can be expressed as: F t = F C X t F = sn δ CX ( t ω X [ ] c ( snδ ωsn( ω ω X The phase dfference between the dampng force and exctng dsplacement s 9, and equvalent dampng F coeffcent C s snδ ω X Dynamc stffness K and equvalent dampng coeffcent C of rubber component can be obtaned after F X and δ be measured Commonly, X and δ are functons of exctng frequencyω So dynamc stffness K and equvalent dampng coeffcent C are also functons of exctng frequencyω [1] After removng the tems relatng wth tme from the equaton of exctng dsplacement and deformaton force, the curve of dynamc exctng dsplacement-deformaton force s: 2X 2 2 X 2 2 X F cosδ + F = sn δ 6 X X F F Ths equaton defnes an ellpse curve, seeng Fg1 Ths s the feature of the lnear dynamc characterstc of rubber component δ can be calculated wth the functon of exctng dsplacement and deformaton force or the coordnates of ABCD 5 Fg1 A typcal ellpse curve of lnear dynamc characterstc of rubber C D ANALYSIS OF NON-LINEAR DYNAMIC CHARACTERISTIC When the exctng force acted on the rubber component exceed a certan range, the curve of exctng dsplacement-deformaton force wll no longer be an ellpse, and that changes gradually to crescent, seeng Fg2 At ths tme, the dynamc characterstc of rubber component analyss method wth the lnear analyss assumpton appears obvous shortcomngs On the bass of B O Y A X 2F 2 Internatonal ADAMS User Conference Page 2 of 5

dynamc testng of rubber component, a numercal fttng method descrbed the non-lnear dynamc characterstc of rubber component whch base on the energy analyss s put forward n ths paper Obvously, the curve of loadng and unloadng s closed The area enclosed by the curves denotes the energy consumed by the dampng of rubber component that provded by external exctng dsplacement Xt ( = X cos( ωt The energy consumed by dampng of rubber component n a crcle s: W = F( t dx( t 2π = F X sn( ωtcos( ωt+ δ d( ωt 7 = πfx snδ Loadng-unloadng curve of rubber component s also a closed curve n the non-lnear dynamc characterstc testng So we can assume: Loadng curve: Fa( x = ax + b x ( x, x t Unloadng curve: Fr( x = cx + d x ( x, x t Where: s nteger and dependent by the precson of vehcle dynamc analyss Now the energy consumed by dampng of rubber component n the non-lnear dynamc characterstc testng s: xt xt W= F( xdx F( xdx x ( x, x t a x x r 8 For ce St f f ness Damp ng Load ng Fg2 Fttng prncple of non-lnear characterstc of rubber component O If the stffness of rubber component s represented wth the mean value of loadng and unloadng curve (seeng Fg2, the non-lnear dynamc characterstc of rubber component can be denoted as: d( [ Fa( x + Fr( x ]/2 K = dx x ( x, x t 9 1 = ( a + c x The functon of dampng relatve to exctng dsplacement can also be obtaned: dw C = = Fa( x Fr( x dx = a c x + b d ( ( Unl oad ng Def or mat on 1 2 Internatonal ADAMS User Conference Page 3 of 5

Force (KN Fg3 Dynamc testng result Deformaton (mm Fg4 Fttng result of dynamc test Deformaton (mm TEST DATA PROCESS AND EXAMPLE ANALYSIS Accordng to the method educed above, a numercal fttng works for the dynamc testng of suspenson rubber bushng of a domestc car s carred out The fttng result s presented by Fg 3~Fg6 Fg3 s the curve of non-lnear exctng dsplacement-deformaton force of rubber component under external exctng, the frequency of external exctng s 11 Hz Fg4 s the curve of dynamc characterstc after fttng It s clear that the fttng result agrees wth the testng result consstently Force (KN Force (KN Force (KN Fg5 Dynamc stffness characterstc Frequency (Hz Deformaton (mm Fg6 Dynamc dampng characterstc Frequency (Hz Deformaton (mm The 3D-surface of exctng frequencydeformaton and dynamc stffness are showng on Fg5 whle the exctng frequency changed from 1 Hz to 2 Hz, t s obtaned by fttng the testng result of rubber component Fg6 shows the 3D-surface of exctng frequencydeformaton and dampng A parameter table of non-lnear dynamc characterstc of rubber component can be bult wthn the frequency range when the external exctng frequency changed By means of judgng the exctng frequency of model, the non-lnear dynamc characterstc of rubber component can be obtaned by the nterpolaton method based on the table As a result, the effect of lnear and non-lnear dynamc characterstc of rubber component to vehcle can be consdered n the process of vehcle dynamc analyss 2 Internatonal ADAMS User Conference Page 4 of 5

Fg7 Analyss result of lnear system Fg8 Analyss result of non-lnear system Accordng to ths method, an ADAMS user subroutne was compled and an ADAMS custom solver was created, and a free vbraton system that conssts of a part supported by a rubber component s analyzed wth the solver The stffness of rubber component s usng the dynamc stffness whle the exctng frequency s 11Hz Accordng to the lnearzaton assume, The stffness of rubber component s 14KN/mm, the dampng rato s 2, and mass of the part s 245Kg The calculaton result s shown n Fg7 Meanwhle, Fg8 shows the analyss result when the nonlnear dynamc characterstc of rubber component was taken nto account In the Fg7 and Fg8, the top drawng s the dsplacement-force curve, the bottom drawng s the vbraton frequencyampltude curve The calculaton result ndcates that the system wll vbrate wth nherent frequency about 11 Hz, and t s a typcal lnear and one degree of freedom vbraton system Consderng the non-lnear stffness characterstc of rubber component, the response frequency s abundance, whch shows the typcal characterstc of non-lnear and one degree of freedom vbraton system, and manfests the effect of non-lnear stffness characterstc of rubber component to system vbraton CONCLUTION The method of energy analyss s smple and practcal, whch can be appled to study the nonlnear dynamc characterstc of rubber component of vehcle suspenson The analyss result of one degree of freedom vbraton system ndcates that the method of dealng wth the non-lnear dynamc characterstc of rubber component of vehcle suspenson put forward n ths paper can be used reasonably to the dynamc analyss of vehcle The study method makes up the shortage of lnearzaton assume of rubber component of vehcle suspenson, and provdes a feasble method for the smulaton of non-lnear dynamc of vehcle REFERENCE 1 Guo Konghu, Influence of tre dynamc sde-slp characterstc on vehcle shmmy, Automoble Technology, 19954 2 [Japan]Hu Yuan Chun Yan, Mu Chuanwen Translate, Antvbraton Rubber and Its Applcaton, Chna Ralway Publshng Company,1982 3 Ca Yanchun, L Chengde, The Transent Response and Analyss of Automoble Non-lnear System, Chanchun Automoble Research Insttute, 1989511 4 Wen Qang, Yu Zhuopng, Zhang Ljun, Zhang Janqang, Wu Langshou, Testng Research of Rubber Bearng of Suspenson and Its Statc and Dynamc Stffness Characterstc Shangha Automoble, 19988 2 Internatonal ADAMS User Conference Page 5 of 5