Multilayered polymeric film vibrations under un-symmetrical loading/Daugiasluoksnes polimerines pleveles virpesiai esant nesimetriniam apkrovimui.
Kibirkstis, E. ; Dabkevicius, A. ; Ragulskis, K. 等
1. Introduction
Multilayered materials can be used for producing different designs
of packages applicable to particle-type and fluid-type food and
industrial products. Graphical images are formed on multilayered
materials by using flexographic printing. The combined material that
finds widest application consists from paper which is coated with films
produced from synthetic materials [1].
In this paper the expression of supplementary stiffness from static
tension is taken into account, which is based on the expression for
plate bending presented in [2]. That expression is transformed to the
nodal displacements of the layers. Thus in this paper only the basic
final relationships from the previous papers, published in journal
"Mechanika" are presented and further development of the model
described there is performed. This paper presents the development of
finite element matrices of multilayered polymeric film when
supplementary stiffness from static tension is taken into account.
In flexographic printing machines, when printing is performed on
tensioned materials while the tape for printing is being transported
between different printing sections, the printing material may be loaded
unsymmetrically [3]. This may take place because of the incorrect
mounting of supports of the directing circular elements of the polymeric
film tape, the wear of bearings, rolls, necks and other reasons.
In the printing device often one side of the printing material
experiences higher loading than the other [4]. During the process of
package production, this has a negative influence on the quality of
graphical image transmission, the strength, barrier and other qualities
of produced packages. Therefore, it is important to investigate
vibrations occurring in the process of un-symmetric loading as well as
their eigenmodes, and to predict possible graphical defects occurring in
the process of printing, thus enhancing the quality of off-prints and
avoiding defected products [5].
Taking international experience into account for example on the
basis of references [6-11] it can be concluded that the applied method
of digital speckle photography is simple, fast and accurate.
By investigating the vibrations of polymeric films and their
eigenmodes, it is possible to determine whether the printing material
which is being transported to the printing device is loaded
symmetrically or unsymmetrically. This method of investigation may also
be applied for the identification and diagnostics of physical and
mechanical properties of paper, paperboard or polymeric films [12].
The purpose of investigations presented in this paper is the
analysis of vibrations as well as the eigenmodes of films produced from
polymeric materials and having a number of layers when the loading of
the material is un-symmetric. The shapes of eigenmodes of the film
vibrations enable us to diagnose the character of the film tension, thus
enhancing the quality of the process of printing.
2. Description of the problem
With the purpose of evaluating the mechanical qualities exhibited
by the polymeric material in the process of non-continuous motion in the
machine used for printing, when the polymeric film experiences uneven
tension in its longitudinal direction, the method of experimental
investigation was designed and applied in order to perform the analysis
for the problem when the polymeric film is loaded non symmetrically. In
this case the distribution of loading is uneven (the direction of
loading of the film produced from polymeric materials is assumed
longitudinal) and is achieved by adding a higher load on one of its
ends. The diagram of symmetric loading is presented in Fig. 1, a and of
un-symmetric loading in Fig. 1, b.
[FIGURE 1 OMITTED]
Layers of the analyzed multilayered polymeric film are shown in
Fig. 2.
[FIGURE 2 OMITTED]
The axes of coordinates (Fig. 2) are denoted as x, y, z; the
corresponding displacements are denoted as u, v, w; [E.sub.1] denotes
the modulus of elasticity of the outer layers, [[mu].sub.1]--the
Poisson's ratio of the outer layers, [[rho].sub.1]--the density of
the material of the outer layers, [h.sub.1]--the thickness of the outer
layers; [E.sub.2]--the modulus of elasticity of the inner layer,
[[mu].sub.2]--the Poisson's ratio of the inner layer,
[[rho].sub.2]--the density of the material of the inner layer,
[h.sub.2]--the thickness of the inner layer; subscripts on u, v, w
denote the surface of the layer (1 corresponds to the lower surface and
2 corresponds to the upper surface of the layer).
Multilayered structures consist from:
1) load carrying layers of high strength and stiffness which
protect from thermal, chemical and other external actions;
2) interconnecting layers which are soft and they redistribute the
load between the carrying layers.
Structures with layers having various qualities ensure reliable
operation of systems under unfavourable environmental conditions. By the
proper choice of layers structures with high strength and stiffness and
at the same time with relatively small mass are designed. For stiff
layers the assumptions of plate theory are valid. For soft layers those
assumptions can not be applied.
Thus the models of the layer of multilayered polymeric films can be
of two basic types:
1) the layer of plate type, in which the following assumption is
made: the displacements of the lower surface of the layer in the
direction of the z axis coincide with the displacements of the upper
surface of the layer in the direction of the z axis (Fig. 2, b);
2) the layer of elastic body type, in which the following
assumption is made: the displacements of the lower surface of the layer
in the direction of the z axis may not coincide with the displacements
of the upper surface of the layer in the direction of the z axis (Fig.
2, c).
Thus having the sub-element of the layer of plate type and the
sub-element of the layer of elastic body type the finite element for
multilayered material of desired type may be obtained.
Three layered structures consisting from two load carrying layers
and the soft filler between them find wide application in engineering
and packaging technology. Such structures are near to optimal in
achieving minimum weight while at the same time satisfying the
requirements for strength and stiffness. Thus further calculations are
performed for the model consisting from the external layers of plate
type and an internal layer of elastic body type (Fig. 2, a).
The node of the finite element consisting from three sub-elements
has 10 degrees of freedom, which include: vertical displacement of the
lower layer, displacements of the lower surface of the lower layer in
the directions of the axes x and y, displacements of the upper surface
of the lower layer in the directions of the axes x and y, vertical
displacement of the upper layer, displacements of the lower surface of
the upper layer in the directions of the axes x and y, displacements of
the upper surface of the upper layer in the directions of the axes x and
y.
2.1. Description of the model of the outer layer of a multilayered
polymeric film
The rotations about the axes of coordinates x and y are denoted as
[[THETA].sub.x] and [[THETA].sub.y]. The displacements because of
bending are given by u = z[[THETA].sub.y] and v = -z[[THETA].sub.x].
Thus
[[THETA].sub.y] = [[u.sub.2] - [u.sub.1]]/[h.sub.1],
[[THETA].sub.x] = [[v.sub.1] - [v.sub.2]]/[h.sub.1].
The mass matrix
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)
where the density of the material is denoted as [[rho].sub.1] and
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)
where the shape functions as usual are denoted as [N.sub.i]. The
stiffness matrix
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)
where the modulus of elasticity is denoted as [E.sub.1], the
Poisson's ratio is denoted as [[mu].sub.1] and also
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)
where stresses [[sigma].sub.x], [[sigma].sub.y], [[tau].sub.xy] are
determined from
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)
where the solution of the static problem is denoted as {[delta]}
2.2. Description of the model of the inner layer of a multilayered
polymeric film
The mass matrix
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)
where the density of the material of the layer is denoted as
[[rho].sub.2], the thickness of the layer is denoted as [h.sub.2] and
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)
The derivatives of w with respect to x and y are expressed as
The stiffness matrix
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)
where the vector of nodal displacements is denoted as {[delta]} and
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)
where
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)
where the bulk modulus is denoted as K = [E.sub.2]/3(1 -
2[[mu].sub.2]), the shear modulus is denoted as G = [E.sub.2]/2(1 +
[[mu].sub.2]). And the stresses in the middle of the layer are
determined from the static problem:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)
3. Analysis of the vibrations of the multilayered polymeric film
The length and width of the polymeric film are chosen equal to 0.2
m. The following boundary conditions are assumed: all displacements are
zero on the lower and upper boundaries, only on the upper boundary
special conditions are assumed such as v = 0 on the left side and v =
0.02 on the right side with linear variation between them.
For the outer layer the following data are assumed: Poisson's
ratio [[mu].sub.1] = 0.3, Young's modulus [E.sub.1] = 0.28 x
[10.sup.9] Pa, density [[rho].sub.1] = 800 kg/[m.sup.3], thickness
[h.sub.1] = 10 [micro]m.
[FIGURE 3 OMITTED]
For the inner layer the following data are assumed: Poisson's
ratio [[mu].sub.2] = 0.3, Young's modulus [E.sub.2] = 2.8 x
[10.sup.9] Pa, density [[rho].sup.2] = 800 kg/[m.sup.3], thickness
[h.sub.2] = 100 [micro]m.
There are rather great uncertainties in the determination of the
physical parameters of the layers of the multilayered polymeric film.
They have direct influence to the obtained results and thus the
correspondence of the experimentally and numerically obtained
eigenmodes. The performed investigations [3-4] showed that the
correspondence between the shapes of the experimental and numerical
eigenmodes is satisfactory, but in order to achieve desirable accuracy
between the eigenfrequencies special experiments for the determination
of precise physical parameters of the layers are required.
Transverse displacement of the lower and upper planes is
investigated. In Figs. 3-9 the obtained results are presented.
The numerical results presented in the figure are used in the
process of interpretation of the experimentally obtained eigenmodes as
described further in this paper. They enable to determine which
eigenmode is obtained in the process of experimental investigations.
4. Description of the method and experimental testing devices of
multilayered materials
In order to determine the character of stress distribution in a
multilayered material in the dynamic motion modes when the loading of
the ' polymeric film tape is un-symmetric, a special setup for
experimental investigation was developed (Figs. 4-5).
In the process of experimental tests, one side of the investigated
polymeric material (Fig. 4) was fastened between the fastening devices
in the electro-dynamical exciter of vibrations. At the same time, weight
load P was un-symmetrically fastened to the other side of the polymeric
film (which was also glued to the fastening element), see Figs. 1, b and
4.
[FIGURE 4 OMITTED]
[FIGURE 5 OMITTED]
Experimental setup shown in Fig. 5 was designed and produced. In
this experimental setup one end of the polymeric film 1 (Fig. 5, b) was
fastened (glued) between the aluminum elbow 2 and the keeping strip 4.
This end was fastened to the vibroexciter. To the other end of the
polymeric film which was also glued to the aluminum elbow the weight
(load) P was attached. The position of this load enables to apply
symmetric as well as non-symmetric loadings. In order to achieve minimum
distortion in the distribution of stresses in the film in the whole
tested area of the film several methods of fastening of the polymeric
film to the keeping aluminum elbows were investigated. It was determined
experimentally that the smallest effect to the non-uniformity of
stresses in the film has the special type of fastening when the film is
glued to the aluminum angle by using a two-sided sticking tape. When
fastening the polymeric film by presses, screwing by screws, gluing with
glues or fastening in some other ways the undesirable non-uniformity of
stresses through the whole width of the keeping elements was observed
(the width is up to 200 mm). Because of this non-uniformity of stresses
the results of investigations become inaccurate. Thus for this purpose
foamed acrylic two-sided sticky tape Tesa 52017 (thickness 0.38 mm,
adhesion 25.67 kg/m) was used.
As shown in the structural diagram of the setup for experimental
investigations (Fig. 5, a) the accelerometer KD--32 is fastened to the
operating element of the vibroexciter --to the vibrating membrane. Thus
the mass of the accelerometer has no effect to the loading of the film
and at the same time to the characteristics of vibrations of the film.
Frequencies of the experimentally determined eigenmodes are given
in Fig. 6 and Table 1. Thus in the process of experimental analysis the
resonance of the accelerometer and of the investigated multilayered film
could not take place.
[FIGURE 6 OMITTED]
The measurements were carried out at such frequencies generated by
the forced vibrator, which coincided with eigenfrequencies of the
investigated film. Measurements were carried out when the oscillation
frequency was stabilized.
The camera Edmund Optics EO-1312C USB was used. In the experimental
tests square samples of the film (0.2 x 0.2 m) were produced.
5. Experimental results
In Fig. 6 some of the results are indicated: in Fig. 6, a--the
eigenmode 2, in Fig. 6, b--the eigenmode 4, in Fig. 6, c--the eigenmode
5 and in Fig. 6, d--the eigenmode 7 of the film PET+PAP+LDPE are shown.
There is correspondence of experimental results with the vibrations
of the polymeric film investigated numerically. Figs. 6, (a-d)
correspond to Fig. 3. In the numerical images of the first four
eigenmodes (see Fig. 3, A-D) horizontal lines dominate, which may be
compared with the images obtained in Figs. 6, a and b. Higher eigenmodes
(see Fig. 3, E and F) have two vertically located lines, which may be
compared with Fig. 6, c. And finally the eigenmode represented in Fig.
3, G has three vertically located lines, which may be compared with Fig.
6, d. Thus one is to compare Fig. 3, B, a with Fig. 6, a, Fig. 3, D, a
with Fig. 6, b, Fig. 3, E, a with Fig. 6, c and Fig. 3, G, a with Fig.
6, d.
The performed investigations of polymeric films were in the range
of frequencies and amplitudes of vibrations which are presented as the
data for Fig. 6 and in the Table 1. The frequency ranges of excitation were from 168 Hz up to 272 Hz.
In the Table 1 some of the results of the experimental
investigations of the film with the method of speckle photography when
the loading of the film was non symmetrical are presented. Intervals of
the values of frequencies and of amplitudes of longitudinal vibrations
of a tape of film, in which after excitation of resonance vibrations the
standing waves (eigenmodes) of transverse vibrations of the tape occur,
are given in the second and third columns of the Table 1.
In our previous papers [4, 12] in which vibrations of paper are
analyzed for the investigation of the resonant vibrations a special
experimental setup was designed. The eigenmodes were determined by using
the projection moire techniques for the tape of paper and for the
symmetric load. The several first experimental eigenmodes well
correspond to the numerically obtained ones. But the higher eigenmodes
do not exhibit such correspondence. Thus similar conclusions previously
obtained for the paper on the basis of the current investigation can be
considered to be valid for polymeric films.
6. Conclusions
Experimental and numerical investigations of vibrations of a
polymeric film were performed. Mutual correspondence of the results is
discussed. Non-symmetric loading takes place in a number of engineering
problems. Thus the applied numerical model and the designed experimental
setup provide the basis for investigation of such problems. Because the
results of numerical and experimental investigations have well
corresponded.
Acknowledgements
The authors are grateful for the financial support by the Project
of the European Union Structural Funds "Postdoctoral Fellowship
Implementation in Lithuania". This Project belongs to the framework
of the Measure for Enhancing Mobility of Scholars and Other Researchers
and the Promotion of Student Research (VP1-3.1-SMM-01). This Project
also belongs to the Program of Human Resources Development Action Plan.
References
[1.] Griffiths, A.J.; Williams, K.P.; Owen, N. 2010. Engineering
value recovery from municipal solid waste. Journal of Mechanical
Engineering Science 224(3): 559-570.
http://dx.doi.org/10.1243/09544062JMES1857.
[2.] Zienkiewicz, O.C. 1975. The Finite Element Method in
Engineering Science. Moscow: Mir, 520 p.
[3.] Kibirkstis, E.; Kabelkaite, A.; Dabkevicius, A.; Ragulskis, L.
2008. Investigation of vibrations of packaging materials, Journal of
Vibroengineering 10(2): 225-235.
[4.] Kibirkstis, E.; Kabelkaite, A.; Dabkevicius, A.; Ragulskis, L.
2007. Investigation of vibrations of a sheet of paper in the printing
machine, Journal of Vibroengineering 9(2): 40-44.
[5.] Aubry, E.; Vaclavik, I.; Ragot, P.; Perrenoud, A. 2006. Method
for active compensation for transverse vibrations of a moving paper web,
Experimental Techniques 30(1): 47-50.
http://dx.doi.org/10.1111/j.1747-1567.2006.00010.x.
[6.] Cloud, G. 2006. Optical methods in experimental mechanics,
Experimental Techniques 30(6): 27-30.
http://dx.doi.org/10.1111/j.1747-1567.2006.00103.x.
[7.] Sjodahl, M. 1997. Accuracy in electronic speckle photography,
Applied Optics 36(13): 2875-2885.
http://dx.doi.org/10.1364/AO.36.002875.
[8.] Larsson, L.; Sjodahl, M.; Thuvander, F. 2004. Microscopic 3-D
displacement field measurements using digital speckle photography,
Optics and Lasers in Engineering 41(5): 767-777.
http://dx.doi.org/10.1016/S0143-8166(03)00028-9.
[9.] Synnergren, P.; Sjodahl, M. 1999. A stereoscopic digital
speckle photography system for 3-D displacement field measurements,
Optics and Lasers in Engineering 31: 425-443.
http://dx.doi.org/10.1016/S0143-8166(99)00040-8.
[10.] Sjodahl, M.; Benckert, L.R. 1993. Electronic speckle
photography: analysis of an algorithm giving the displacement with
subpixel accuracy, Applied Optics 32: 2278-2284.
http://dx.doi.org/10.1364/AO.32.002278.
[11.] Jin Seung, K. 1989. Range and accuracy of speckle
displacement measurement in double-exposure speckle photography, The
Journal of the Optical Society of America A 6: 675-681.
http://dx.doi.org/10.1364/JOSAA.6.000675.
[12.] Kibirkstis, E.; Kabelkaite, A.; Dabkevicius, A.; Bivainis,
V.; Ragulskis, L. 2009. Non-destructive diagnostics of uniformity of
loading of the sheet of paper, Ultrasound 64(3): 24-28.
E. Kibirkstis, Kaunas University of Technology, Studentu 56, 51424
Kaunas, Lithuania, E-mail: edmundas.kibirkstis@ktu.lt
A. Dabkevicius Kaunas University of Technology, Studentu 56, 51424
Kaunas, Lithuania, E-mail: arturas.dabkevicius@ktu.lt
K. Ragulskis, Lithuanian Academy of Sciences and Kaunas University
of Technology, Lithuania, E-mail: kazimieras3@hotmail.com,
kazimieras3@yahoo.com
V. Miliunas, Kaunas University of Technology, Studentu 56, 51424
Kaunas, Lithuania, E-mail: valdas.miliunas@ktu.lt
V. Bivainis, Kaunas University of Technology, Studentu 56, 51424
Kaunas, Lithuania, E-mail: vaidas.bivainis@ktu.lt
L. Ragulskis, Vytautas Magnus University, K. Donelaicio g. 58,
44248 Kaunas, Lithuania, E-mail: l.ragulskis@if.vdu.lt
http://dx.doi.org/10.5755/j01.mech.18.5.2706
Received September21, 2011
Accepted October 19, 2012
Table 1
Frequency ranges of vibrations and their amplitude
Ranges for eigenmodes of multilayered polymeric film
Eigenmode Frequency Amplitude
number range, Hz range, m
I 168-176 1.2-5.9 x [10.sup.-6]
II 178-187 1.2-5.4 x [10.sup.-6]
III 189-198 0.9-4.8 x [10.sup.-6]
IV 200-210 0.7-3.3 x [10.sup.-6]
V 212-227 0.6-2.6 x [10.sup.-6]
VI 230-248 0.4-1.5 x [10.sup.-6]
VII 250-272 0.2-1.0 x [10.sup.-6]