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 connection with the Euler-Bernoulli equation. The problem is solved by reducing it to a system of differential inclusions.


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][3][4][5][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][9][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][12][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][15][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][20][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][25][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][29][30][31][32][33][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.Jx ′ the dimensionless moment of inertia.The cross-sectional area A is non-dimensionalized by 2  L , and the time t by    , Uxt , the equation of motion for the beam is written as, by omitting the prime

Motion Equation with Hysteresis
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 ( ) ,, Cxt ξτ − , which is assumed to be separable in space and time 3) The expression (2.2) and ( 2.
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 ( ) , Uxt 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 ( ) 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 ex-ample, 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 ( ) cx ξ − 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. .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 ( ) with µ the time relaxation constant of the auxetic foams and 0 g a constant.Both constants depend on the Poisson's ratio (Figure 5).
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 : 12 ,R γγ ∈ with 12 γγ < (Figure 6).This operator is not closed.For any continuous function  ( . 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 k ε The regularized relay operator k ε γ may be written an the weak form a.e. in 0, 10, The hysteresis motion Equation (2.8) may be rewritten, by taking We also assume that 0 ε ≥ , and provide a unified formulation of our problem for both the completed relay operator k γ and its regularization k ε γ with 0 ε > .In this way, the problem can be formulated as Find ,0, in 0,, ,,,,dd This problem may then be interpreted as a weak formulation of the problem (2.8), (2.9).Henceforth we shall write ( ) , uw in place of ( ) , uw εε .The hysteresis relay operator can be either continuous or discontinuous, but in this paper only the continuous case is studied.

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 1 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 ( ) , uw for a cyclic variation of u are displayed for both cases of spatial and the temporal hysteresis, respecttively.The calculations are carried out for 0 uuu = % and 0 wWw = % 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.One element -4.97 ± 20.34i -0.35 ± 73.56i -0.050 ± 161.53i -0.023 ± 289.71i -0.0054 ± 455.73i Two elements -5.26 ± 22.84i -0.99 ± 74.57i -0.095 ± 161.40i -0.045 ± 291.54i -0.0108 ± 445.65iA 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].

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.

1 1 ω
, where 1 ω is the first natural frequency.The potential energy is non-dimensionalized by 0 EIL , and the density of the bar ρ is non-dimensionalized by 62 01 EIL ω .The virtual work due to the damping force per unit length is non-dimensionalizing by 0 EIL and the damping force by 3 0 LEI .The dimensionless non-conservative work is expressed as damping force per unit length defined over the spatial subdo-

Figure 2 .L
Figure 2. The beam with damping elements.
3) represent the general form of the nonlocal damping force.The Heaviside function ( ) Hx denotes the presence of damping.We have ( ) 0 HxH = (constant) if x is within the patch, and ( ) 0 Hx= otherwise.The spatial kernel function ( ) case of (2.3) refers to the spatial hysteresis with ( ) ( )

Figure 4 Figure 3 .Figure 4 .
Figure 3.Comparison between compressive stress-strain curves of the considered types of foams.
, respectively.In these expressions, α and β are characteristic parameters of the damping material, and 0 l 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/cm 3 ) 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/cm 3 ) after 20000 cycles, the error function and respectively, exponential decay function are closer to an average value of 27 mJ/cm 3 , and respectively, 29 mJ/cm 3 .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 1 of the damping of the auxetic foams suggests to consider another particular case of (2.3), i.e. the temporal hysteresis, in which ( ) cx ξ − may be the Dirac delta function ( ) x δξ − , which reflects the reacting character of the damping force 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 ( )
and the elastic modulus E = 259.93MPa.