Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 83
PROCEEDINGS OF THE EIGHTH INTERNATIONAL CONFERENCE ON COMPUTATIONAL STRUCTURES TECHNOLOGY Edited by: B.H.V. Topping, G. Montero and R. Montenegro
Paper 75
On the Dynamics of a Duffing Oscillator with an Exponential NonViscous Damping Model D.J. Wagg and S. Adhikari
Faculty of Engineering, University of Bristol, United Kingdom D.J. Wagg, S. Adhikari, "On the Dynamics of a Duffing Oscillator with an Exponential NonViscous Damping Model", in B.H.V. Topping, G. Montero, R. Montenegro, (Editors), "Proceedings of the Eighth International Conference on Computational Structures Technology", CivilComp Press, Stirlingshire, UK, Paper 75, 2006. doi:10.4203/ccp.83.75
Keywords: damping, Duffing, nonviscous, bifurcation.
Summary
In this paper a Duffing oscillator with nonviscous damping function is considered.
The nonviscous damping function is an exponential damping model with a decaying memory property to the damping term of the oscillator.
Viscous damping is the most common model for modelling of
vibration damping. This model, first introduced by Lord Rayleigh
[1], assumes that the instantaneous generalized
velocities are the only relevant variables that determine damping.
Viscous damping models are used widely for their simplicity and
mathematical convenience even though the behaviour of real
structural materials may not be adequately represented by viscous
models. For this reason it is well recognized that in general a
physically realistic model of damping will not be viscous. Damping
models in which the dissipative forces depend on any quantity
other than the instantaneous generalized velocities are
nonviscous damping models.
Mathematically, any causal model which
makes the energy dissipation functional nonnegative is a possible
candidate for a nonviscous damping model.
Clearly a wide range of
choice is possible, either based on the physics of the problem, or
by a priori selecting a model and fitting its parameters from
experiments.
Among various damping models, the 'exponential
damping model' is particularly promising and has been used by many
authors. With this model the damping force is expressed as
Here is the viscous damping constant, μ is the relaxation parameter and is the displacement as a function of time. In the context of viscoelastic materials, the physical basis for exponential models has been well established, see for example [2]. A selected literature review including the justifications for considering exponential damping model may be found in reference [3]. The Duffing oscillator with viscous damping has been widely studied, and comprehensive reviews of the associated literature can be found in [4]. For these type of systems the damping is nearly always assumed to be viscous. One exception is a study of study of a Duffing vibration isolator which has combined Coulomb and viscous damping [5]. In this study the authors demonstrate that the addition of a Coulomb damping term can lead to new dynamical behaviour in the system. In Section 2 of the paper we present the mathematical model of the system, which is then partly nondimensionalised and reduced to an equivalent system of three first order differential equations. The reduced system has six system parameters, of which two relate to the external forcing, two capture the stiffness behaviour, one represents the viscous damping and one the nonviscous damping. In Section 3 we show numerical simulations of the system in both phase space, and parameter space. The reduced set of three ODE's can be simulated using a fixed step, 4th order RungeKutta algorithm. The numerical parameter variation is used to create bifurcation diagrams, in this case using the external forcing frequency as the bifurcation parameter  this allows easy comparison with work on the viscous Duffing system. Structural changes in the bifurcation diagrams are then observed for increasing levels of nonviscous damping parameter. The two cases considered relate to small and large cubic stiffness nonlinearity. In the case of small nonlinearity, hardening spring behaviour is found, and the result of adding nonviscous effects is to increase the amplitude of the resonance peak. In the case of strong nonlinearity, more complex dynamics are observed close to the first natural frequency. Here the effect of adding nonviscous damping is to cause subtle changes in the dynamics. These changes include subtle movement or shifts in the location of the different regions of dynamical behaviour along the frequency axis. Changes in the periodicity of solutions can also be observed, although it should be noted that our numerical approach is not guaranteed to capture all coexisting solutions. For both cases, significant changes only occur for large values of nonviscous damping. For general nonlinear system modelling of mechanical or structural applications, this type of nonviscous damping may give improved accuracy of the simulations when compared to experimental data. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
