Elsevier

Applied Acoustics

Volume 67, Issue 2, February 2006, Pages 93-105
Applied Acoustics

Simulation of the nonlinear vibration of a string using the Cellular Automata based on the reflection rule

https://doi.org/10.1016/j.apacoust.2005.05.008Get rights and content

Abstract

In this study, the nonlinear dynamic responses of a string are simulated using the Cellular Automata method based on the reflection rule. In the case of nonlinear systems, the velocity of wave propagation is not constant and depends on the amplitude. A new treatment of the dynamic time step is proposed for the Cellular Automata method considering the effect of the propagation velocity. As numerical examples, first, the dynamic responses of a string with linear characteristic are simulated using the Cellular Automata method. A typical resonance curve can be obtained. Second, the dynamic responses of a string with nonlinear characteristic are simulated using the proposed method. Some characteristic types of vibration can be obtained. It is concluded that the linear and nonlinear dynamic responses of a string may be obtained by simulation using the Cellular Automata method.

Introduction

Recently, some machines or structures become lighter in spite of its severe operating condition. As their vibrations may be large, a nonlinear vibration analysis becomes important for precise analysis and design. A nonlinear analysis of the dynamics of a multi-degree-of-freedom system or continuous system can be carried out, for instance, using the method of multiple scales or the Harmonic Balance method [1]. Many studies have been found analyzing the nonlinear phenomena by these methods. These methods are very effective, though its procedures for computation are complex and fixing the order for solution is not easy. One of the other methods for nonlinear dynamic analysis is the numerical integration of a nonlinear differential equation, however, this requires much computation time in the case of a multi-degree-of-freedom system.

A discrete computation method, such as the Cellular Automata (CA) method or the Lattice Gas Automata method, have been recently introduced to analyze some problems in the engineering field [2]. In the CA method, the analytical domain is divided into finite state variables called ‘cells’. The state of each cell is updated according to local rules at every discrete time step. That is, the state of a cell at a given time step depends only on its previous state and that of the neighbor cells. The states of all cells are updated synchronously. Because of such computational characteristics, analysis can be performed for only a desired portion of the total space. A local rule is designed in such a way that the results will satisfy the requirements. These methods, for instance, have been applied to sound wave propagation [3], [4] and the fluid dynamics problem [5].

The dynamic analysis of a string is one of the problems regarding wave propagation. Local rules for wave propagation have been previously proposed [3], [4]. There are some studies regarding wave propagation [6], [7], [8], [9] based on the rule proposed in [3]. The rule is applied for the amplitude of vibration and the decay rate of amplitude, however, it is feasible to only the linear problem. Applying the CA method, the wave will propagate over a distance between cells in one time step, thus the propagation velocity must coincide real velocity estimated from the material properties of the media. For the nonlinear dynamic analysis of a string, the tension may depend on its amplitude of vibration, thus the propagation velocity is not constant. The local rules for a linear problem must be modified for a nonlinear one. One of the authors has been proposed a new CA method for nonlinear dynamic analysis [10] based on the rule in [3]. In that proposed method, the time step for simulation in the CA method is adjusted to real time by considering the effect of the propagation velocity. Both the linear and nonlinear dynamic responses of a string can be demonstrated using the proposed method. On the other hand, the rule proposed in [4] is based on the collision between cells, and the elastic collision between some cells is simulated. Due to the same reasons as mentioned above, this rule can be applied to only the linear problem.

In this study, the nonlinear dynamic responses of a string are simulated using the CA method based on the rule proposed in [4]. The rule in [4] is used for the transverse vibration and a new treatment of the dynamic time step will be proposed for nonlinear dynamic analysis. In our proposed method, the time step, usually constant, is adjusted to real time by considering the effect of the propagation velocity. The basic idea is as same as the one in the previous study [10], but, the concrete procedures are different. As a numerical example, first, the dynamic response of a string having linear characteristics is calculated using the CA method. Second, the dynamic response of a string having nonlinear characteristics is calculated. The nonlinear term related to the tension of the string is calculated using the dynamic response at the past time of one time step. The results are examined from the viewpoint of their physical meaning.

Section snippets

Setting state parameters

A string can be thought of as a chain of masses linked by springs. It is composed of two kinds of cells, the white ones and the black ones, which alternate along a line. Fig. 1 shows a string composed of some cells. Each cell has state parameter such as the amplitude of vibration u, which is defined as real numbers.

General local rules

The rule of update is the following. The time evolution has two phases. First, the black cells are held fixed and the white ones move according to the rule below. Then, the white

Simulation model

The nonlinear response of a uniform string with simple supports, which is the similar one described in [12], is treated. The string has a length l = 0.8 m, width b = 20 × 10−3 m, thickness h = 0.5 × 10−3 m, density ρ = 7.8 × 103 kg/m3, Young’s modulus E = 206 GPa, damping coefficient per unit length C = 0.14 N s/m2 and initial tension N0 = 130 N. The external force per unit length q is assumed to consist of two components, one independent on, and the other harmonically varying with time, i.e., it is given in the

Conclusions

In this study, the nonlinear dynamic responses of a string are simulated using the CA method based on the reflection rule. In the case of nonlinear systems, the velocity of wave propagation is not constant and depends on the amplitude. A new treatment of the dynamic time step is proposed for the CA method considering the effect of the propagation velocity.

As numerical examples, first, the dynamic responses of a string with linear characteristic are simulated using the CA method. A typical

References (12)

There are more references available in the full text version of this article.

Cited by (6)

  • Cellular automata in nonlinear vibration problems of two-parameter elastic foundation

    2011, Archives of Civil and Mechanical Engineering
    Citation Excerpt :

    Considering that the work of load q(x, t) and Winkler foundation reaction -k(x)w(x, t) for the displacements shown in Figure 2 is a linear displacement function and the energy dissipation function is a homogenous quadratic velocity function, the current load in relation (4) is f(i, j + 1) = q (i, j + 1) – k(i)w(i, j + 1) and the damping per unit string length parameter in relation is (7) c = ⅔C, which is due to the zero velocities of the balls directly neighbouring the ball moving in one of the stages of string evolution [22] According to the proposal by Kawamura et al. [28], time-variable velocity v(t) is the basis for the use of evolution-variable time increment Δt (changing the values of parameters δ and χ) so matched that the disturbance wave front covers exactly distance a = l/n (constant throughout the entire analysis) during this time. A numerical analysis of the behaviour of an exemplary two-parameter nonlinear elastic foundation with length l = 80.0 m and a width of 1.0 m is carried out here for several test problems whose schemes are shown in Figure 3.

  • Cellular automata in nonlinear string vibration

    2010, Archives of Civil and Mechanical Engineering
    Citation Excerpt :

    The mathematical Kirchhoff model [33] of nonlinear transverse vibration, neglecting the displacements along the string's axis and averaging tension (N) over its length (l), can be in an elementary way derived from Newton's second law and written as Applying central difference approximation for ∂2w/∂t2, ∂2w/∂x2 and ∂w/∂t one can write Equation (5), after elementary transformations, in the form presented in the papers by Kawamura et al. [31–32] For the above relations it can be shown that (1 − 2δ) = (γ − 1)/(γ + 1), which fully agrees with the relations presented by Kawamura et al. [31–32].

1

Master course student of Toyohashi University of Technology.

View full text