On the Damped Beams with Hysteresis


This paper discusses the hysteretic behavior of beams with external elements made from auxetic materials. The damping force is modeled by using the nonlocal theory. Unlike the local models, the damping force is modeled as a weighted average of the velocity field over the temporal and spatial domains, determined by a kernel function based on distance measures. The hysteresis operator is continuous and it is defined in con-nection with the Euler-Bernoulli equation. The problem is solved by reducing it to a system of differential inclusions.

Share and Cite:

M. Poienariu, M. Ionescu, I. Girip, L. Munteanu and V. Chiroiu, "On the Damped Beams with Hysteresis," World Journal of Mechanics, Vol. 1 No. 1, 2011, pp. 6-14. doi: 10.4236/wjm.2011.11002.

1. Introduction

The term auxetic is coming from the Greek word auxetos, meaning that which may be increased. Instead of getting thinner like an elongated elastic band when stretched, an auxetic material gains volume, expanding laterally [1]. For non-homogeneous linear elastic materials, it is demonstrated that even for the simplest loading case, i.e. quasi-static uniaxial, the Poisson’s ratio is spatial and temporal dependent and not a constant [2]. The Poisson’s ratio can be considered an auxetic descriptor only on such spatial and temporal domains where it gains negative values. Generally, the Poisson’s ratio is a function of the spatially and temporally changing moduli and stresses, and varies accordingly. For instance, we have


Auxetic materials have a negative Poisson’s ratio for large spatial; and temporal domains and are expected to have interesting mechanical properties, such as high energy absorption, fracture toughness, indentation resistance and enhanced shear moduli, which may be useful in some applications [2-6]. Scientists have been aware of the existence of auxetic materials for over a hundred years, though without very special attention, and treating them as an accident or a curiosity. In the case of an isotropic material, the range of Poisson’s ratio is from –1.0 to 0.5, based on thermodynamic considerations of strain energy in the theory of elasticity. The auxetic behavior is found in materials from molecular and microscopic levels up to the macroscopic level. Love [7] presented an example of a cubic single crystal pyrite as having the Poisson’s ratio of –0.14, and he suggested that the effect may be caused by twinned crystals. Subsequently, an auxetic behavior has been observed in other naturally singlecrystal materials [8-10], in real materials with a high degree of anisotropy, such as conventional honeycomb network, re-entrant honeycomb and hexagonal structures (figure 1), reticulated metal foams, the skin covering a cow’s teats, certain rocks and minerals, living bone tissue, etc. In the auxetic materials the shear modulus is significantly greater then the elastic Young modulus It is well known that all major classes of materials (polymers, composites, metals and ceramics) can exist in auxetic form. A specific feature exhibeted by auxetic materials in comparison with other foams is their significant damping capacity at various loading levels, with increase up to 16 times compared to conventional foams [11-13]. The size-dependent theories of the damping capacity of auxetic foams have received increasing attention in recent years due to the need to model and analyse very small-sized mechanical interactions inn the rapid developments of micro-nanotechnologies [14-16]. The nonlocal theory becomes significant when dealing with such problems. The advantages of this theory over the classical local theory are in the high accuracy which can be obtained when employing the energy loss per unit volume versus the number of hysteretic cycles. The nonlocal estimation of the energy loss shows a better approximation of the experimental results than the local energy loss [15,16].

The starting point in the nonlocal damping theory of structures is the damping force which depends at a given point on the past history of a velocity field over a certain domain, through a kernel function. The nonlocal damping force is modeled as a weighted average of the velocity field over the spatial domain, and it is determined by a kernel function based on distance measures [17]. The damping capacity of the auxetic foam, in particular the conventional grey open-cells polyurethane foam [18], is tested by adding it under the form of the external elements to conventional materials [19-21].

Such models represent a generalization of the viscous damping, and examples include the beams and nanobeams with external patches [22,23], the auxetic composites [24-26], the micropaddles with the periodic auxetic cores [27], the multifunctional nanofoils based on the carbon nanotubes and the auxetic foams [14]. A nonlocal damping model including the time and spatial hysteresis effects for beams and plates is analysed in [14] by using the Galerkin method. The hysteresis in connection with PDEs and applicative problems are analyzed in [28-34]. Details on the nonlinear semigroups and differential equations in Banach spaces can be found in [35].

In this paper we apply the nonlocal theory to analyze the hysteretic behavior of the Euler-Bernoulli beams with external elements made from the auxetic material. The hysteresis operator is continuous and it is defined in connection with the motion equation of the beam.

2. Motion Equation with Hysteresis

Let us consider a beam of length L, in which a number of external damping elements of thickness are attached at , , , as shown in figure 2. The length L, the transverse displacement, and the coordinate x are non-dimensionalized by the length of the beam, so that dimensionless quantities are given by .

The bending stiffness is, where E is the Young’s modulus of elasticity and the moment of inertia, is a constant moment of inertia, and the dimensionless moment of inertia. The cross-sectional area A is non-dimensionalized by, and the time t by, where is the first natural frequency. The potential energy is non-dimensionalized by, and the density of the bar is non-dimensionalized by. The virtual work due to the damping force per unit length is non-dimensionalizing by and the damping force by. The dimensionless non-conservative work is expressed as

, where is the external damping force per unit length defined over the spatial subdomains,.


Figure 1. Conventional honeycomb network, re-entrant honeycomb and hexagonal structures, with negative Poisson’s ratio.

Figure 2. The beam with damping elements.

The state of the beam is characterized by two variables confined to a set: the external damping force acting over the spatial subdomains, , and the transverse displacement. Using the extended Hamilton’s principle, by including the non-conservative damping force, the equation of motion for the beam is written as, by omitting the prime




The comma in (2.1) represents the differentiation with respect to the specified variable. The Poisson ratio is an important factor in defining the kernel, which is assumed to be separable in space and time


The expression (2.2) and (2.3) represent the general form of the nonlocal damping force. The Heaviside function denotes the presence of damping. We have (constant) if x is within the patch, and otherwise. The spatial kernel function

 is normalized to satisfy the condition.

A particular case of (2.3) refers to the spatial hysteresis with where is the delta function in time. In this case, the force depends only on the instantaneous value of the velocity or strain rate, but depends on the spatial distribution of the velocities. Also, the velocities at different locations within a certain domain can affect the damping force at a given point.

The external damping force is related to the peculiar deformation mechanisms of the auxetic materials. In fact, the auxetic materials expand in all directions when pulled in only one, leading to an unusual volume change of the solid itself. The auxetic behaviour does not contradict the classical theory of elasticity. A homogeneous, isotropic material has a Poisson’s ratio range between –1.0 and 0.5, while the magnitude of the Poisson’s ratio can be even larger in the case of anisotropic solids. Negative Poisson’s ratio materials can exhibit slow decay of stress according to Saint-Venant’s principle. The choise of the function in (3) depends on the Poisson’s ratio. To understand this, some previous results are further presented [20,21]. The idea was to make a comparison between the non and auxetic grey open-cells polyurethane foams from the point of view of the damping capacities. The Poisson’s ratio values are all strain dependent, tending towards near-zero, negative or positive values depending on the pore dimensions. For example, the conventional foam has a Poisson’s ratio positive of 0.25 at compressive strain of 10%, which decreases sharply with the increase of compressive loading, to become slightly negative from 60 to 80% of tensile strain. The auxetic foam exhibits a negative Poisson’s ratio of –0.185 at compressive strain from 10 to 25%, showing a sharp increase for rising compressive strain, reaching then a zero value at 55% of compressive strain and a positive Poisson’s ratio of 1.33 at 80% [36]. Figure 3 shows the compressive constitutive law of the conventional and auxetic foams for cylinders having the diameter of 30 mm and length of 170 mm.

We determine the best form for so that the variation of the energy loss per unit volume with respect to the number of compressive cyclic loadings is the closest to the experimental data [36,37]. The polyurethane sample is obtained from cylinders having the diameter of 30 mm and length of 170 mm. Then, it is compressed inside the mould obtaining a final nominal diameter of 20 mm and length of 100 mm. Figure 4 defines the compressive loading and displays the variation of the Poisson’s ratio with respect to compressive strain for this

Figure 3. Comparison between compressive stress–strain curves of the considered types of foams.

Figure 4. Poisson’s ratio versus compressive strain for auxetic foam [21].

sample [21]. For any given loading cycle, the dissipated energy is, where and are the minimum and maximum strains, respectively, and. It is interesting to notice the higher energy absorption characteristics of the auxetic foam for the same strain ranges compared to the non-auxetic one.

Accordingly to [3], the calculations of the damping capacities for both non and auxetic foams are carried out on four forms of the spatial kernel: the exponential decay function, the error function, , the hat function, for and 0 otherwise, and the triangular function, for

and 0 otherwise, respectively. In these expressions, and are characteristic parameters of the damping material, and is the influence distance parameter. The numerical results are quite satisfying. The following conclusions can be drawn: 1) for nonauxetic foams the exponential decay function and triangular functions yield to similar values of energy loss per unit volume (15 mJ/cm3) after 20000 cycles; 2) all functions show a good correlation of the energy loss per unit volume with experimental data for auxetic materials; 3) while the hat and triangular functions yield to similar values of energy loss per unit volume (mJ/cm3) after 20000 cycles, the error function and respectively, exponential decay function are closer to an average value of 27 mJ/cm3, and respectively, 29 mJ/cm3. One feature of all results is the high energy dissipated by auxetic foam compared to the conventional foam for the same loading level.

Therefore, we choose for the spatial kernel function the exponential decay with and


The reacting character of the damping of the auxetic foams suggests to consider another particular case of (2.3), i.e. the temporal hysteresis, in which may be the Dirac delta function, which reflects the reacting character of the damping force


In this case we


The ability to dissipate energy is one of the main reasons for using auxetic materials. Stress relaxation refers to the behavior of stress reaching a peak and then decreasing or relaxing over time under a fixed level of strain. For the kernel function depending on the past time histories, we consider


with the time relaxation constant of the auxetic foams and a constant. Both constants depend on the Poisson’s ratio (figure 5).

In this paper we deal with the Euler-Bernouli problem (2.1) and (2.2) that contain a continuous hysteresis operator



The initial conditions attached to (2.8) are given by


The boundary conditions are written, for example, for a simple supported beam


In view of inserting of the delayed relay operator in (2.1) as expressed in (2.8), let us fix any pair with (figure 6). This operator is not closed. For any continuous function is uniformly continuous, hence it may have a finite number of oscillations between the thresholds. For any and, we set and


Figure 5. Poisson’s ratio versus and.

Figure 6. Relay operator.


with [29-31]. We define the rate-independent hysteresis operator

in the sense of Visintin [29-31]. For any increasing continuous function, if.we have.

Visintin has defined the completed relay operator, where is the set of parts of (figure 7). This operator is defines as follows. For any and any, we set if and only if is measurable in



In conformity to (3.4), w is constant in a neighbourhood of t in the first case, w is nonincreasing in a neighbourhood of t in the second case, and w is nondecreasing in a neighbourhood of t in the last case.

The completed relay operator can be regularized by a regularized relay operator

Figure 7. Completed relay operator.

The regularized relay operator may be written an the weak form





The hysteresis motion Equation (2.8) may be rewritten, by taking


We also assume that, and provide a unified formulation of our problem for both the completed relay operator and its regularization with. In this way, the problem can be formulated as Find and such that a.e. in and



This problem may then be interpreted as a weak formulation of the problem (2.8), (2.9). Henceforth we shall write in place of.

The hysteresis relay operator can be either continuous or discontinuous, but in this paper only the continuous case is studied.

3. Results and Conclusions

In what follows, the example refers to a simply supported aluminum beam with constant diameter d = 0.0025, with a single auxetic element, , , and respectively, two auxetic elements , , , , with thickness. The auxetic material has the density, the Poisson’s ratio and the elastic modulus E = 259.93 MPa.

Tables 1 and 2 tabulate the lower estimates for the first five eigenvalues for the beam with a single and respectively, two damping elements. It is observed that the temporal hysteresis gives larger damping ratios then the spatial hysteresis. These results agree to the previous results reported in [21,22].

It is the purpose of this section to present the per-formance of the method of including the relay operator into the Euler-Bernoulli beam with auxetic external elements. To do this, we solve numerically the system of differential inclusions associated to the (2.20) and (2.21). As shown in figures 8 and 9, the hysteresis loops for the ratios for a cyclic variation of u are displayed for both cases of spatial and the temporal hysteresis, respecttively. The calculations are carried out for and in the middle of the beam, for the first four vibration modes in the case of the spatial hysteresis and respectively, for first four vibration modes in the case of the temporal hysteresis.

According to figures 8 and 9, it should be mentioned that the temporal hysteresis gives larger hysteresis loops then the spatial hysteresis. Moreover, the spatial loops look to be in some way to be stretched in the vertical direction.

Table 1. First five eigenvalues for a simply supported beam (spatial hysteresis).

Table 2. First five eigenvalues for a simply supported beam (temporal hysteresis).

Figure 8. The temporal hysteresis for four vibration modes.

Figure 9. The spatial hysteresis for four vibration modes.

A possible drawback of this method is the difficulty in extending the method to deal with discontinuous relay hysteresis operator. The main idea is to explain some features of the experimental results in which the singularities of the hysteresis loops and the presence of the unclosed loops are reported.

For instance, in the figure 10 we see that the displacement drift, the force relaxation and non-closure of hysteretic loops are present [38]. We remember that the displacement drift appears when cycled between two unequal forces, the force relaxation appears when cycled between two unequal displacements. We must specify that such responses exhibit multiple reversals of small amplitude and are experimentally put into evidence [39].

Figure 10. The displacement drift and force relaxation [38].

4. Conclusions

In this paper, the continuous relay hysteresis operator is illustrated, and a weak formulation is studied for the nonlocal Euler-Bernoulli equation of the beam with external elements made from auxetic materials. The damping force is modeled as a weighted average of the velocity field over the temporal and spatial domains, determined by a kernel function based on distance measures. The spatial and temporal hysteresis are defined when employing the presence of the auxetic elements. The problem is solved by reducing it to a system of differential inclusions. From the numerical results presented in this study, we can conclude that the temporal hysteresis is associated with the larger hysteresis loops then in the case of the spatial hysteresis.

5. Acknowledgements

This research was financially supported by the National Authority for Scientific Research (ANCS, UEFISCSU), Romania, through PN2 project ID_1391/2008.

Conflicts of Interest

The authors declare no conflicts of interest.


[1] K. E. Evans, M. A. Nkansah, I. J. Hutchinson and S. C. Rogers, “Molecular Network Design,” Nature, Vol. 353, 1991, pp. 124-125. doi:10.1038/353124a0
[2] R. S. Lakes, “Experimental Microelasticity of Two Porous Solids,” International Journal of Solids and Structures, Vol. 22, No. 1, 1986, pp. 55-63. doi:10.1016/0020- 7683(86)90103-4
[3] R. S. Lakes, “Foam Structures with a Negative Poisson’s Ratio,” Science, Vol. 235, 1987, pp. 1038-1040. doi:10. 1126/science.235.4792.1038
[4] R. S. Lakes, “Experimental Micro Mechanics Methods for Conventional and Negative Poisson’s Ratio Cellular Solids as Cosserat Continua,” Journal of Engineering Materials and Technology, Vol. 113, No. 1, 1991, pp. 148-155. doi:10.1115/1.2903371
[5] D. W. Overaker, L. M. Cuiti?o and N. A. Langrana, “Effects of Morphology and Orientation on the Behavior of Two-Dimensional Hexagonal Foams and Application in a Re-Entrant Foam Anchor Model,” Mechanics of Materials, Vol. 29, No. 1, June 1998, pp. 43-52. doi:10.1016/ S0167-6636(98)00004-0
[6] Y.-C. Wang and R. Lakes, “Analytical Parametric Analysis of the Contact Problem of Human Buttocks and Negative Poisson’s Ratio Foam Cushions,” International Journal of Solids and Structures, Vol. 39, No. 18, September 2003, pp. 4825-4838. doi:10.1016/S0020-7683 (02)00379-7
[7] A. E. H. Love, “A Treatise on the Mathematical Theory of Elasticity,” Dover, 4th Edition, New York, 1926.
[8] D. J. Gunton and G. A. Saunders, “Stability Limits on the Poisson Ratio,” Journal of Material Science, Vol. 7, 1972, pp. 1061-1068. doi:10.1007/BF00550070
[9] Y. Li, “The Anisotropic Behavior of Poisson’s Ratio, Young’s Modulus, and Shear Modulus in Hexagonal Materials,” Physica Status Solidi A, Vol. 38, No. 1, November 1976, pp. 171-175. doi:10.1002/pssa.2210380119
[10] R. H. Baughman, J. M. Shacklette, A. A. Zakhidov and S. Stafstrom, “Negative Poisson’s Ratios as a Common Feature of Cubic Metals,” Nature, Vol. 392, 1998, pp. 362-365. doi:10.1038/32842
[11] F. Scarpa, J. A. Giacomin, A. Bezazi and W. A. Bullough, “Dynamic Behaviour and Damping Capacity of Auxetic Foam Pads,” SPIE Proceeding, 2006.
[12] ?. Donescu, V. Chiroiu amd L. Munteanu, “On the Young’s Modulus of a Auxetic Composite Structure,” Mechanics Research Communications, Vol. 36, No. 3, 2009, pp. 294-301. doi:10.1016/j.mechrescom.2008.10.0 06
[13] ?. Donescu, L. Munteanu, P. P. Delsanto and V. Mo?negu?u, “On the Advanced Auxetic Composites,” Research Trends in Mechanics, Vol. 3, 2009, pp. 78-93.
[14] Y. Lei, M. I. Friswell and S. Adhikari, “A Galerkin Method for Distributed Systems with Nonlocal Damp- ing,” International Journal of Solids and Structures, Vol. 43, No. 11-12, June 2006, pp. 3381-3400. doi:10.1016/j. ijsolstr.2005.06.058
[15] M. I. Friswell, S. Adhikari and Y. Lei, “Non-Local Finite Element Analysis of Damped Beams,” International Journal of Solids and Structures, Vol. 44, No. 22-23, 2007, pp. 7564-7576. doi:10.1016/j.ijsolstr.2007.04.023
[16] M. I Friswell, S. Adhikari and Y. Lei, “Vibration Analysis of Beams with Non-Local Foundations Using the Finite Element Method,” International Journal of Numerical Methods in Engineering, Vol. 71, No. 11, 2007, pp. 1365-1386. doi:10.1002/nme.2003
[17] W. Flugge,“Viscoelasticity,” 2nd Editon, Springer-Verlag, Berlin, 1975.
[18] F. Scarpa, P. Pastorino, A. Garelli, S. Patsias and M. Ruzzene, “Auxetic Compliant Flexible PU Foams: Static and Dynamic Properties,” Physica Status Solidi B, Vol. 242, No. 3, 2005, pp. 681-694. doi:10.1002/pssb.200460 386
[19] L. Munteanu, D. Dumitriu, ?. Donescu and V. Chiroiu, “On the Complexity of the Auxetic Systems,” Proceed- ings of the European Computing Conference, Vol. 2, 2009, pp. 631-636.
[20] V. Chiroiu, St. Donescu, L. Munteanu and V. Mo?negu?u, “The Dynamics of Beams with Auxetic Patches,” Proceedings of the International Conference on Advanced Materials for Application in Acoustics and Vibration AMAAV’09, Cairo, January 2009.
[21] V. Chiroiu, L. Munteanu and St. Donescu, “On the Beams with External Auxetic Patches,” Advances in Mechanical Engineering, 2009, pp. 1-10. doi:10.1155/2009/ 430379
[22] L. Munteanu, V. Chiroiu, D. Dumitriu, D. Baldovin, St. Donescu and C. Chiroiu, “On the Eigenvalues Optimization of Euler-Bernoulli Beams with Nonlocal Damping Pathes,” Révue Roumaine des Sciences Techniques, série de Mécanique Appliquée, Vol. 54, No. 1, 2009, pp. 53-66.
[23] L. Munteanu, V. Chiroiu and P. P. Teodorescu, “Nano- beams with Damping Elements,” In: M. Zaharescu, M. Ciurea, I. Kleps, D. Dascalu, Eds., Series in Micro and Nanoengineering – Nanostructuring and Nanocharacterization, Vol. 16, 2010, pp. 241-256.
[24] L. Munteanu, P. P. Delsanto, D. Dumitriu and V. Mosnegutu, “On the Characterization of Auxetic Materi- als,” In: D. Popa, V. Chiroiu, I. Toma, Eds, Research Trends in Mechanics, Vol. 2, 2008, pp. 205-234.
[25] L. Munteanu, V. Chiroiu, D. Dumitriu and M. Beldiman, “On the Characterization of Auxetic Composites,” Proceedings of the Romanian Academy, Series A: Mathematics, Physics, Technical Sciences, Information Science, Vol. 9, No. 1, 2008, pp. 33-40.
[26] V. Chiroiu, L. Munteanu, D. Dumitriu, M. Beldiman and C. Secara, “On the Arhitecture of a New Cellular Elastic Material,” Proceedings of the Romanian Academy, Series A: Mathematics, Physics, Technical Sciences, Information Science, Vol. 9, No. 2, 2008, pp. 105-115.
[27] L. Munteanu, V. Chiroiu, T. Sireteanu and L. Tenek, “Vibrations of a Micropaddle with Periodic Auxetic Core,” Révue Roumaine des Sciences Techniques, série de Mécanique Appliquée, Vol. 55, No. 2, 2010, pp. 101- 114.
[28] M. Beldiman, L. Munteanu and M. Poienariu, “On the Multifunctional Nanofoils Based on Carbon Nanotubes and Auxetic Foams,” In: L. Munteanu, V. Chiroiu, T. Sireteanu, Eds., Research Trends in Mechanics, Vol. 4, 2010, pp. 9-24.
[29] A. Visintin, “Hysteresis and Semigroups,” In: A. Visintin, Ed., Models of Hysteresis, Longman, Harlow, 1993, pp. 192-206.
[30] A. Visintin, “Differential Models of Hysteresis,” Springer- Verlag, Berlin, 1995.
[31] A. Visintin, “Quasi-Linear Hyperbolic Equations with Hysteresis,” Annales de l’ Institut Henri Poincaré, Non- linear Analysis, Vol. 19, No. 4, 2002, pp. 451-476.
[32] Y. Kōmura, “Nonlinear Semi-Groups in Hilbert Space,” Journal of the Mathematical Society of Japan, Vol. 19, No. 4, 1967, pp. 493-507. doi:10.2969/jmsj/01940493
[33] J. Kopfová, “Nonlinear Semigroup Methods in Problems with Hysteresis,” Discrete and Continuous Dynamical Systems Supplement, pp. 580-589, 2007
[34] M. C. Crandall and T. M. Liggett, “Generation of Semigroups of Nonlinear Transformations on General Banach Spaces,” American Journal of Mathematics, Vol. 93, No. 2, April 1971, pp. 265-298. doi:10.2307/2373376
[35] V. Barbu, “Nonlinear Semigroups and Differential Equations in Banach Spaces,” Noordhoff, Leyden, 1976.
[36] A. Bezazi and F. Scarpa, “Mechanical Behaviour of Conventional and Negative Poisson’s Ratio Thermoplastic Polyurethane Foams under Compressive Cyclic Loading,” International Journal of Fatigue, Vol. 29, No. 5, 2007, pp. 922-930. doi:10.1016/j.ijfatigue.2006.07.015
[37] F. Scarpa, J. A. Giacomin, A. Bezazi and W. A. Bullough, “Dynamic Behaviour and Damping Capacity of Auxetic Foam Pads,” SPIE Proceedings of Smart Structures and Materials, 2006.
[38] V. Preda, M. F. Ionescu, V. Chiroiu and T. Sireteanu, “A Preisach Model for the Analysis of the Hysteretic Phenomena,” Révue Roumaine des Sciences Techniques, série de Mécanique Appliquée, Vol. 55, No. 3, 2010, pp. 75-86.
[39] A. E. Charalampakis and V. K. Koumousis, “A Bouc- Wen Model Compatible with Plasticity Postulates,” Jour- nal of Sound and Vibration, Vol. 322, 2009, pp. 954-968, doi:10.1016/j.jsv.2008.11.017

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.