## Journal of the Brazilian Society of Mechanical Sciences and Engineering

*versão On-line* ISSN 1806-3691

### J. Braz. Soc. Mech. Sci. & Eng. v.25 n.1 Rio de Janeiro jan./mar. 2003

#### http://dx.doi.org/10.1590/S1678-58782003000100013

**Chaotic motion of an elasto-plastic beam**

**M. A. Savi ^{I, *}; P. M. C. L. Pacheco**

^{I}Department of Mechanical and Materials Engineering, Instituto Militar de Engenharia, 22.290.270 Rio de Janeiro, RJ. Brazil. E-mail: savi@ime.eb.br

^{II}Department of Mechanical Engineering, CEFET/RJ, 20.271.110 Rio de Janeiro, RJ. Brazil. E-mail: calas@cefet-rj.br

**ABSTRACT**

The non-linear dynamics of a pin-ended elasto-plastic beam is discussed in this contribution. Ideal plasticity is in focus. Free and forced responses are considered. An iterative numerical procedure, based on the operator split technique, is developed. This work shows that jump phenomenon and sensitivity to initial conditions imply that the system response may become unpredictable even when a periodic response is expected. Chaotic motion is also associated with the unpredictability of the beam response.

**Keywords:** Chaos, non-linear dynamics, elasto-plasticity

**Introduction**

This contribution investigates the non-linear dynamics of a pin-ended elasto-plastic beam. The beam is modeled by the Symonds' model where a pin-ended beam with length 2*L*, and uniform rectangular cross section of area *A = bh'*, is represented by two rigid links, each of length *L*, joined by an elasto-plastic element. The two rigid bars are assumed to have mass per unit length r , the same as for the uniform beam (Shanley, 1947). The beam model is depicted in Fig.1.

An iterative numerical procedure, based on the operator split technique, is developed. Numerical simulations discuss the free and forced responses of the beam. Some dynamical characteristics of the beam imply that the system response may become unpredictable. Chaotic motion is one of these and it may occur for some parameters. There is also jump phenomenon which implies that small changes in force amplitudes cause great variations in the steady state responses.

**Model for an Elasto-Plastic Beam**

Geometric considerations allows one to define the relation between the cell position, *a*, and the angle of rotation, j, and also the semi-extension of the cell centerline, *e*, and the semi-extension of each flange, *e _{1}* and

*e*,

_{2}The constitutive equation for ideal plasticity is given by,

where sign(s) = s / |s|. s is the one-dimensional stress, e and e * ^{p}* is the total and plastic one-dimensional strain, respectively. g represents the rate at which plastic deformations take place.

*E*is the Young modulus. The yield function,

*h*(s), the Kuhn-Tucker conditions and the consistency condition are given by (Savi & Pacheco, 1997):

Here, s* _{y}* is the yield stress.

The force and moment resultants in the cell are taken by considering the same relations of those of a sandwich beam, consisting of two bars each in simple tension or compression (Symonds & Yu, 1985). Hence,

where *N _{1}* and

*N*are the forces on each flange. Assuming that the total strain on each flange is obtained by dividing the semi-extension by the semi-length of the beam, and the area of each flange is a half of the beam cross section area, hence,

_{2}where is the plastic semi-extension on each flange.

Governing equations of the model are obtained by establishing the equilibrium of moment on the half beam. Neglecting the inertia of the elasto-plastic element, and assuming a linear viscous external dissipation,

where µ = 3/rL^{3} and *c* is the linear viscous dissipation parameter. *N* and *M* are given by equations (1,5). Now, consider the following definitions,

Here and Denoting the non-dimensional time derivative by ( )´= d( )/dt, the following system can be written,

Actually, the space of state variables includes more variables than *y _{1}* and

*y*(Poddar

_{2}*et al.*, 1988), however, the analysis is developed on a subspace of dimension 2 (Savi & Pacheco, 1997). The numerical solution procedure here proposed uses the operator split technique (Ortiz

*et al.*, 1983). First, the equations (8) are integrated using any classical scheme, like fourth order Runge-Kutta, assuming that the variables

*n*and

*m*are known parameters.

*n*and

*m*are evaluated by considering an elastic predictor step (trial state), where plastic variable, e

*, remains constant from the previous time instant. The next step of solution procedure consists on a plastic corrector step where the feasibility of trial state is evaluated using the return mapping algorithm (Simo & Taylor, 1985). An iterative process takes place until the convergence is achieved (Savi & Pacheco, 1997).*

^{p}

**Free Vibrations**

In this Section, one discusses the free response of the elasto-plastic beam. This is done by letting ƒ(t vanish in the equations of motion (8). In all simulations, one has taken *L* = 0.10m, *b* = 0.02m, *h'* = 0.04m, *h* = 0.68*h'*, r = 0.216 kg/m, *E* = 120 GPa, s* _{y}* = 0.3 GPa (Poddar

*et al.*, 1988). The procedure converges with time steps less than Dt = 2p / 1000W.

To illustrate the free response, one considers a system with no external dissipation (*c _{0}* = 0). Figure 2 presents results from simulations in the form of phase portrait. Different initial conditions cause different plastification of the cell and, as a consequence, alter the position of the equilibrium points in phase space.

The beam response presents sensitivity to initial conditions. Poddar *et al.* (1988) show the fractal basin boundaries of the system under free vibration. This conclusion was used to explain the discrepancy among some finite element results shown in Symonds & Yu (1985). Figure 3 shows the steady state envelope of maximum and minimum angular displacement, *y _{1}*, for different initial conditions. A null initial angular velocity is considered. Figure 3a considers no external dissipation (

*c*= 0), and reproduces the result obtained by Symonds & Yu (1985). It should be pointed out that a global change on the response occurs when initial conditions are on the range from 0.086 to 0.092. By considering an external dissipation, asymptotic behavior is expected. Hence, the response tends to be a closed curve in phase space, and the envelope becomes a line. When

_{0}*c*= 0.2, global changes are characterized by jumps (Fig.3b). For high values of this parameter, for example

_{0}*c*= 1.5, global changes do not occur anymore (Fig.3c).

_{0}

**Forced Vibrations**

Forced vibration analysis is now in focus. First, one considers a square wave excitation with amplitude d. Actually, the considered square wave uses 1% of the period to load or unload. One starts the analysis by considering bifurcation diagrams (Fig.4). The dissipation parameter is *c*_{0} = 1.5 and different frequency parameters are considered. The first 30 cycles are neglected.

In order to analyze the jump phenomenon for W= 1, one considers two values of forcing amplitude near amplitude d= 0.22, say d= 0.22 and d= 0.23. Figure 5a shows the steady state response for these two situations. For W= 2, high values of forcing amplitudes presents different behavior on bifurcation diagram. After jump, which occurs near d= 0.37, one value of the forcing amplitude is related with many points. This behavior is associated with transient response since the steady state is periodic. Figure 5b shows the steady state response for two forcing amplitudes very close: d= 0.37 and d= 0.38. There is a global change on dynamical response.

In spite of the beam does not exhibit chaotic motion with the considered physical parameters, jump phenomenon and sensitivity to initial conditions introduces difficulties to predict the beam behavior. Poddar *et al.* (1988) show that chaotic motion occurs in Symonds' beam when it is subjected to periodic impulses. Hence, it is interesting to consider different parameters to evaluate the possibility of chaotic motion. With this aim, one considers a dissipation parameter *c _{0}* = 0.2. Figure 6 shows the bifurcation diagram with frequency parameter W= 1. Now, it is possible to identify regions with cloud of points associated with chaos.

The strange attractor of the motion for d= 0.08 is depicted in Fig.7a, while Fig.7b presents the periodic phase plane for d= 0.04. The Fast Fourier Transform (FFT) analysis permits to clearly identify the difference between the two responses (Fig.8). As it is well known, the FFT of a chaotic signal presents continuous spectra over a limited range. The energy is spread over a wider bandwidth. On the other hand, the FFT of a periodic signal presents discrete spectra, where a finite number of frequencies contribute for the response (Moon, 1992; Mullin, 1993).

Chaotic motion may also occur when the beam is subjected to a harmonic excitation. In order to analyze this situation, one considers a harmonic sinusoidal excitation, with a dissipation parameter *c _{0}* = 0.2. Figure 9 shows the bifurcation diagram with frequency parameter W= 1. As in the previous case, it is possible to identify regions with cloud of points associated with chaos. Figure 10 shows the strange attractor of the motion for d= 0.13, while Fig.10b presents the periodic phase plane for d= 0.15. The Fast Fourier Transform (FFT) analysis permits to identify a clearly difference between the two responses (Fig.11).

**Conclusions**

A dynamic analysis of a pin-ended elasto-plastic beam, described by the Symonds model, is considered. Ideal plasticity is focused. The numerical method proposed in this work permits the use of a combination of classical algorithms to evaluate the response of elasto-plastic dynamical systems. Numerical simulations of free and forced vibrations are evaluated. Forced vibrations analysis shows that jump phenomenon occurs for some parameters. In this situation, a very close change in forcing parameter causes a qualitative change in the beam response. Chaotic motion is also possible in the beam response. These behaviors may cause practical problems to predict the response of the beam, even when periodic response is expected.

**Acknowledgements**

The authors would like to acknowledge the support of the Brazilian Research Council (CNPq) and of the Research Foundation of Rio de Janeiro (FAPERJ).

**References**

Moon, F.C., 1992, "Chaotic and Fractal Dynamics", John Wiley. [ Links ]

Mullin, T., 1993, "The Nature of Chaos", Clarendon Press. [ Links ]

Ortiz, M., Pinsky, P.M. & Taylor, R.L., 1983, "Operator Split Methods for the Numerical Solution of the Elastoplastic Dynamic Problem", Computer Methods in Applied Mechanical Engineering, Vol. 39, pp. 137-157. [ Links ]

Poddar, B., Moon, F.C. & Mukherjee, S., 1988, "Chaotic Motion of an Elastic-Plastic Beam", ASME Journal of Applied Mechanics, Vol. 55, pp. 185-189. [ Links ]

Pratap, R., Mukherjee, S. & Moon, F.C., 1994, "Dynamic Behavior of a Bilinear Hysteretic Elasto-Plastic Oscillator, Part I: Free Oscillations", Journal of Sound and Vibration, Vol. 172(3), pp. 321-337. [ Links ]

Pratap, R., Mukherjee, S. & Moon, F.C., 1994, "Dynamic Behavior of a Bilinear Hysteretic Elasto-Plastic Oscillator, Part II: Oscillations Under Periodic Impulse Forcing", Journal of Sound and Vibration, Vol. 172(3), pp. 339-358. [ Links ]

Savi, M.A. & Pacheco, P.M.C.L., 1997, "Non-Linear Dynamics of an Elasto-Plastic Oscillator with Kinematic and Isotropic Hardening", Journal of Sound and Vibration, Vol. 207(2), pp.207-226. [ Links ]

Shanley, F.R., 1947, "Inelastic Column Theory", Journal of the Aeronautical Sciences, Vol. 14, pp. 261-267. [ Links ]

Simo, J.C. & Taylor, R.L., 1985, "Consistent Tangent Operators for Rate-Independent Elastoplasticity", Computer Methods in Applied Mechanical Engineering, Vol. 48, pp. 101-118. [ Links ]

Symonds, P.S. & Yu, T.X., 1985, "Counterintuitive Behavior in a Problem of Elastic-Plastic Beam Dynamics", ASME Journal of Applied Mechanics, Vol. 52, pp. 517-522. [ Links ]

Presented at COBEM 99 – 15th Brazilian Congress of Mechanical Engineering, 22-26 November 1999, São Paulo. SP. Brazil

Paper accepted: July, 2002. Technical Editor: José Roberto de França Arruda

* Current Address: Department of Mechanical Engineering, COPPE / UFRJ, 21945-970 Rio de Janeiro, RJ. Brazil. E-mail: savi@ufrj.br