Study On Improvement of Damping in Jointed Cantilever Beams Using FEM

ISSN: 2319-8753
International Journal of Innovative Research in Science,

Engineering and Technology
(An ISO 3297: 2007 Certified Organization)
Vol. 3, Issue 6, June 2014
Study on Improvement of Damping in Jointed

Cantilever Beams Using FEM
R.C. Mohanty 1, Rajendra Kumar Mohanty 2, B.K. Nanda 3
1
Associate Professor, Department of Mechanical Engineering, Centurion University of Technology and Management,
Jatni-752050, Orissa, India
2
Research Scholar, Department of Mechanical Engineering, Centurion University of Technology and Management,
Jatni-752050, Orissa, India
3
Professor and Vice Chancellor, Veer Surendra Sai University of Technolgy, Burla-768018, Orissa, India
ABSTRACT: Damping in built-up structures is produced by the energy dissipation due to micro-slip along the
frictional interfaces. The analysis of the problem has been carried out using Finite Element Method (FEM). A finite
element model of the linear elastic system has been formulated using the Euler-Bernoulli beam theory to investigate the
damping phenomena in riveted connections. The discrete element system having two degrees of freedom per node
representing v and v x has been used for the analysis. The solution considers one-dimensional beam elements with
each one consisting of two nodes having two degrees of freedom, i.e. transverse displacement and rotation at each node.
The generalized stiffness and mass matrices for this element has been derived. Extensive experiments have been
conducted for the validation of the analysis. From this study, it is established that the damping capacity increases and
the natural frequency decreases due to the joint effects.
KEYWORDS: Shape function; stiffness matrix; interface pressure; micro-slip; damping capacity.
I. INTRODUCTION
Usually, the fabricated or built-up structures possess low material damping which can be adequately compensated by
incorporating suitable joints to construct the same. The contribution of damping due to joints to such structures depends
largely upon its behavior under dynamic conditions. Gropper [1] and Menq et al. [2] have shown that joints usually
exhibit two types of motions; micro-slip and macro-slip. The micro-slip usually occurs with small relative displacement
between the contacting members resulting in energy losses at the joints causing the damping of the structures. Beards
[3] has performed a series of experiments on the damping of the joints to control the vibration of a structure and
established that the damping in joints is much larger compared to the material damping and there exists an optimum
clamping force at which the energy dissipation is maximum. Shin et al. [4] have shown that the damping effect is
negligible when the joint is very tight. They further established that the natural frequency decreases and the damping
capacity increases with the loosening of the joint. Damisa et al. [5] have studied the effect of non-uniform interface
pressure distribution on the mechanism of slip damping for layered beams. However, their analysis is limited to static
load and linear pressure profile. Further, Olunloyo et al. [6] have also studied the damping mechanism with polynomial
and hyperbolic forms of interfacial pressure distribution.
Many boundary and initial value structural engineering problems are solved with the application of finite
element method (FEM). However, its application in damping analysis is relatively recent. Gaul and Lenz [7] have
worked in details on the finite element models of different slip mechanisms to study the dynamic response of
assembled structures incorporating both the micro- and macro-slip of several bolted/riveted joints. Chen and Deng [8]
have proposed a finite element method for understanding and characterizing the non-linear damping behavior of
structural joints. Hartwigsen et al. [9] have found out the contact area of bolted joint interfaces using FEM analysis and
they further conducted the experiments to verify the same. Sainsbury and Zhang [10] have used finite element analysis
through Galerkin element method (GEM) to study the dynamic analysis of damped sandwich beam structures. . Lee et
Copyright to IJIRSET www.ijirset.com 13030

ISSN: 2319-8753

al. [11] used the finite element model of a jointed beam to obtain the natural frequencies and mode shapes. Oldfield et
al. [12] have investigated the damping of an assembled structure using the dynamic frictional contact analysis under
harmonic loading by the finite element analysis.
Although, a lot of analytical, computational and experimental works have been carried out by several researchers in the
recent past on the damping of bolted structures, but no substantial work has been reported till date on the damping
capacity of riveted structures. Hence, an attempt has been made in the present investigation to analyze the damping
mechanism in assembled structures using the finite element models and to study the effect of various parameters on the
damping of layered and jointed cantilever beams having thickness ratio of 1.5. Further, the finite element model
analysis has been validated with experiments.
The finite element method is widely employed to solve many boundary and initial value engineering problems and its
application in damping analysis has been studied by many researchers [1-4]. Among the work done by many authors,
Nanda and Behera [5] have studied damping in layered and jointed bolted structures using the classical energy
approach. However, a very little amount of work has been reported till date on the damping characteristics of riveted
joints. Further, the bolted and riveted joints exhibit different damping behavior because of varying interface pressure
distribution characteristics and zone of influence under similar dynamic conditions. A finite element approach based on
the energy dissipation mechanism has been presented to investigate the damping in jointed and riveted structures.
II. FINITE ELEMENT FORMULATION

A standard discretized beam element has been modeled using two nodes with each one consisting of two degrees of
freedom (translation and rotation) as shown in Fig. 1. Assuming that each layer has the same transverse displacement
and rotation at any section, the deformation within the eth element v(x,t) is approximated using the standard cubic shape
function S(x) as;
v( x, t ) S( x)d e (t ) (1)
T
where the vector d e (t ) v1 (t ) 1 (t ) v2 (t ) 2 (t ) contains the nodal displacements and the matrix
3x 2 2 x3 2 x2 x3 3x2 2 x3 x 2 x3
S( x ) 1 2 3 x 2 3 defines the cubic shape functions.
l l l l l2 l l l 2
Fig.1 Beam element and its nodal displacements

The contact is developed at the interfaces due to the clamping of the two members by rivets. Minakuchi et al [13] have
found that the interface pressure distribution due to this contact is parabolic having a circular influence zone
circumscribing the rivet with diameter equal to 4.125 times the rivet diameter. The interface pressure under each rivet
in a non-dimensional polynomial is expressed as;
10 8 6 4 2
p R R R R R
C1 C2 C3 C4 C5 C6 (2)
S RB RB RB RB RB

ISSN: 2319-8753

where p, s, R and RB are the interface pressure, surface stress, any radius within the influencing zone and radius of the
connecting rivet respectively. The constants C1 to C6 have been evaluated from the numerical data of Minakuchi et al.
[13] by curve fitting with MATLAB as; 0.14581228 E 05, 0.57951686 E 04, 0.60446153 E 03,
0.30852824 E 02, 0.95814172E 01 and 0.53777732E 00 respectively.
The maximum frictional force at the interfaces is given by;
FrM N (3)
where and N are the coefficient of friction and normal load respectively. The normal load is calculated from Eq. (2).
Fig. 2 Mechanism of dynamic slip

When the beam is excited at the free end, the relative dynamic slip under a connecting rivet as shown in Fig. 2 is given
by;
v( x, t ) dS
ur ( x, t ) 2 h 2 h d e (4)
x dx
where is the dynamic slip ratio.
Considering the influence of q number of equi-spaced rivets, the overall maximum relative dynamic slip of a jointed
cantilever beam is found to be;
urM 2 hX sum (5)

q
dS
where X sum dx d
i 1
e
The energy loss per cycle of vibration is given by;
E f 2 FrM urM (6)
The strain energy stored per cycle of vibration due to free end excitation is given as;
1 T
En D KD (7)
2
where K and D are the global stiffness matrix and displacement vector respectively. The matrix K is obtained
assembling the element stiffness matrices of all beam elements, while the stiffness matrix of a beam element is given
by;

ISSN: 2319-8753

12 6l 12 6l
6l 4l 2 6l 2l 2
EI
ke 3 (8)
l 12 6l 12 6l
2
6l 2l 6l 4l 2
The logarithmic decrement () is evaluated from the ratio of the energy loss to the input strain energy per cycle of
vibration as;
Ef 4 NhX sum
(9)
2 En D T KD
A generalized equation for a jointed cantilever beam consisting of m number of layers is obtained modifying Eq. (9)
as;
4 m 1 NhX sum
(10)
DT KD
Further, the product of . is evaluated from Eq. (9) as;
D T KD
(11)
4 NhX sum
III. EXPERIMENTAL SET-UP AND EXPERIMENTS
Fig. 3 Schematic diagram of the experimental set-up

A schematic diagram of the experimental set-up for conducting experiments is shown in Fig. 3. Initially, the Youngs
modulus of elasticity (E) is found out from a sample solid beam specimen by static bending deflection test and its

ISSN: 2319-8753

average value is found as 196.8 GPa. The specimens are prepared from mild steel flats using power driven rivets as
presented in Table 1. The distance between the consecutive rivets, width and length of specimens are varied as per the
zone of influence. The tests are conducted after mounting the specimens rigidly to the support. A digital storage
oscilloscope coupled with a transducer records the oscillations and frequency of vibration. The logarithmic decrement
has been found out using the expression ln a1 an 1 n , where a1, an+1 and n are the amplitudes of the first cycle,
last cycle and number of cycles respectively. The average value has been taken by repeating each test at least for five
times to minimize the experimental errors.
Table 1 Details of mild steel specimens used in the experiment
Thickness x width (mm Number of layers Diameter of rivet Number Cantilever lengths
x mm) ( mm) of rivets (mm)
8 330
(6+6) x 41.25 2 9 371.25
(4+4+4) x 41.25 3 10 10 412.5
(3+3+3+3) x 41.25 4 11 453.75
(6+6) x 24.75 2 15 371.25

(4+4+4) x 24.75 3 6 16 396
(3+3+3+3) x 24.75 4 17 420.75
18 445.5
The energy dissipation principally depends upon the kinematic coefficient of friction () and dynamic slip ratio () at
the interfaces. These two parameters exhibit complicated behavior under dynamic conditions. Moreover, they are inter-
dependent and inversely related. However, their product (.) being frequency and amplitude dependent is established
to be constant for a particular specimen irrespective of any surface condition. In order to authenticate this, a number of
jointed beams with varying surface roughness at the interfaces have been tested and the logarithmic decrement is found
to be almost the same. Hence, this product has been determined from the experimental results of logarithmic decrement
at a particular rivet diameter using Eq. (11), whose variation with frequency and amplitude is presented in Fig. 4. These
values are subsequently used to find out the damping with respect to other diameters using Eq. (10).
Fig.4 Variation of . with frequency and amplitude of vibration
IV. DISCUSSION OF RESULTS

The finite element approach is an efficient tool for the dynamic analysis of jointed structures. As the direct formation of
damping matrix for interfacial slip damping is very difficult, an alternate energy approach has been used to derive the

ISSN: 2319-8753

damping as presented in Eq. (10). In this regard, the expressions for dynamic slip, input strain energy and energy loss in
a jointed beam are formulated using the stiffness matrix, shape function and displacement vector.
Experiments are conducted with a number of specimens to obtain the logarithmic decrements corresponding to 10 mm
diameter rivet. From these results, the product . corresponding to each value of logarithmic decrements is
computed using Eq. (11). Further, the logarithmic decrements for the specimens connected with different rivet
diameters are found out numerically using Eq. (10). The plots presented in Figs. 5 and 6 give a good comparison
between both the results with a maximum variation of 2.3 %. The higher values of the experimental results are due to
material and support damping which have been neglected in the theoretical analysis. The influence of the following
parameters on the damping is explained below:
Effect of the rivet diameter:- The use of the rivets of larger diameter increases the normal force and energy dissipation
at the interfaces. Moreover, the stiffness increases with the rivet diameter introducing more input strain energy as
evident from Eq. (7). But the energy dissipation due to friction occurs at a higher rate compared to the input strain
energy, thereby increasing the logarithmic decrement as presented in Fig. 5.
Fig. 5 Variation of logarithmic decrement with number of layers at 0.1 mm amplitude with 12 mm overall beam
thickness and 371.25 mm cantilever length
Number of Layers
Fig. 6 Variation of logarithmic decrement with number of layers at 0.1 mm amplitude with 12 mm overall beam
thickness and 10 mm diameter connecting rivet

ISSN: 2319-8753

Fig. 7 Variation of logarithmic decrement with beam length with 12 mm overall beam thickness and 10 mm diameter
connecting rivet
Effect of the number of layers:- The use of more number of layers with the same overall beam thickness is
accompanied by more friction surfaces producing higher energy loss. However, the stiffness as well as strain
energy is reduced with the increased number of layers, thereby increasing the logarithmic decrement as shown in
Figs. 5 and 6.
Effect of the specimen length:- The stiffness is decreased with the increase in the specimen length, thus reducing
the input strain energy into the system as evident from Eq. (7). Moreover, the longer specimens require more
number of rivets, thereby increasing the overall dynamic slip and energy loss at the interfaces. This further
enhances the logarithmic decrement as evident from Fig. 6 and 7.
V. CONCLUSION
Considering one-dimensional beam elements, a finite element procedure based on energy approach has been
enumerated in the present analysis. The results so obtained establish that a substantial increase in the damping capacity
can be achieved by using rivets of larger diameter with more number of layers and cantilever length. The aircraft and
aerospace structures, bridges, frames, machine members, machine tools, automobiles, etc. can be effectively designed
with this concept to maximize their damping capacity.
NOMENCLATURE
D = global displacement vector
2h = thickness of each layer of the cantilever specimen
I = second moment of area of the jointed beam cross-section
K = global stiffness matrix
l = length of individual elements
L = free length of the layered and jointed beam
S = shape function matrix
de = nodal displacement vector of an element
ke = element stiffness matrix
v1, v2 = transverse displacement at the nodes of an element

ISSN: 2319-8753

1 , 2 = slope at the nodes of an element
REFERENCES
[1] Groper, M., Microslip and Macroslip in Bolted Joints, Experimental Mechanics, June, pp. 171174, 1985.
[2] Menq, C. H., Bielak, J., and Griffin, J.H., "The Influence of Micro-Slip on Vibratory Response, Part-I: A New Micro-Slip Model,
Journal of Sound and Vibration, Vol. 107, No. 2, pp. 279-293, 1996.
[3] Beards, C. F., Damping in Structural Joints, The Shock and Vibration Digest, Vol. 24, pp. 3-7, 1992.
[4] Shin, Y. S., Iverson, J. C., and Kim, K. S., Experimental Studies on Damping Characteristics of Bolted Joints for Plates and Shells,
Journal of Pressure Vessel Technology, Vol. 113, No. August, pp. 402-408, 1991.
[5] Damisa, O., Olunloyo, V. O. S., Osheku, C. A., and Oyediran, A. A., Dynamic Analysis of Slip Damping in Clamped Layered Beams
with Non-Uniform Pressure Distribution at the Interface, Journal of Sound and Vibration, Vol. 309, pp. 349-374, 2008.
[6] Olunloyo, V. O. S., Damisa, O., Osheku, C. A., and Oyediran, A.A., Further Results on Static Analysis of Slip Damping with Clamped
Laminated Beams, European Journal of scientific Research, Vol. 17, No. 4, pp. 491-508, 2007.
[7] Gaul, L., and Lenz, J., Nonlinear Dynamics of Structures Assembled by Bolted Joints, Acta Mechanics, Vol. 125, No. 169-181, 1997.
[8] Chen, W., and Deng, X., Structural Damping Caused by Micro-Slip along Frictional Interfaces, International Journal of Mechanical
Science, Vol. 47, pp. 1191-1211, 2005.
[9] Hartwigsen, C. J., Song, Y., Macfarland, D. M., Bergman, L. A., and Vakakis, L. A., Experimental Study of Non-Linear Effects in a
Typical Shear Lap Joint Confifuration, Journal of Sound and Vibration, Vol. 277, pp. 327-351, 2004.
[10] Sainsbury, M. G., and Zhang, Q. J., The Galerkin Element Method Applied to the Vibration of Damped Sandwich Beams, Computers &
Structures, Vol. 7, pp. 239-256, 1999.
[11] Lee, S. Y., Ko, K. H., and Lee, J. M., Analysis of Dynamic Characteristics of Structural Joints using Stiffness Influence Coefficients,
KSME Int. J., Vol. 14, No. 12, pp. 1319-1327, 2000.
[12] Oldfield, M., Ouyang, H., and Mottershead, J. E., Simplified Models of Bolted Joints under Harmonic Loading, Computers and
Structures, Vol. 84, pp. 25-33, 2005.
[13] Minakuchi, Y., Koizumi, T., and Shibuya, T., Contact Pressure Measurement by means of Ultrasonic Waves Using Angle Probes,
Bulletin of JSME, Vol. 28, No. 243, pp. 1859-1863, 1985.

Study On Improvement of Damping in Jointed Cantilever Beams Using FEM

Uploaded by

Copyright:

Available Formats

Study On Improvement of Damping in Jointed Cantilever Beams Using FEM

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Study On Improvement of Damping in Jointed Cantilever Beams Using FEM

Uploaded by

Copyright:

Available Formats

ISSN: 2319-8753

International Journal of Innovative Research in Science,

Vol. 3, Issue 6, June 2014