Microsoft Word - numero_49_art_72ex15_2452 S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 814 Focused on Russian mechanics contributions for Structural Integrity Analysis of spatial vibrations of piezoceramic eccentric cylindrical shells interacting with an annular fluid layer S.A. Bochkarev, S.V. Lekomtsev, A.N. Senin Institute of Continuous Media Mechanics, Ural Branch Russian Academy of Sciences, Perm (Russian Federation) bochkarev@icmm.ru, http://orcid.org/0000-0002-9722-1269 lekomtsev@icmm.ru, http://orcid.org/0000-0002-8331-2979 senin.a@icmm.ru, https://orcid.org/0000-0002-7537-0001 ABSTRACT. The work is devoted to a numerical study of the behavior of vertically oriented eccentric (non-coaxial) electroelastic shells, containing a fluid in the gap between them. The solution of the problem is found in the framework of a three-dimensional formulation using the finite element method. The shells are made of the piezoelectric material polarized in the radial direction. Thin-walled structures are considered in the framework of the classical theory based on the Kirchhoff – Love hypotheses, as well as the equations of linear electroelasticity. The distribution of the electric potential through the thickness is assumed linear. The fluid is considered within the framework of a potential theory. The study of the low natural frequencies and vibration modes is carried out for different variants of the boundary conditions at the edges of the shells, the level of fluid in the gap and the eccentricity of the inner shell. The influence of the electric boundary conditions specified on the parts of the shells surfaces covered with electrodes is estimated. KEYWORDS. Piezoceramics; Potential fluid; Partial filling; Eccentric cylindrical shells; Finite element method; Natural vibrations. Citation: Bochkarev, S.A., Lekomtsev, S.V., Senin A.N., Analysis of spatial vibrations of piezoceramic eccentric cylindrical shells interacting with an annular fluid layer, Frattura ed Integrità Strutturale, 49 (2019) 814-830. Received: 01.04.2019 Accepted: 22.05.2019 Published: 01.07.2019 Copyright: © 2019 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION n the last few decades, the wide-spread use of structural elements in the form of nested shells with parallel axes, interacting with a fluid or gas contained in the gap between them, gave an impetus to the appearance of many publications, in which such systems are the subject of research. An extensive bibliography, covering recent studies into the dynamic behavior of coaxial shells interacting with both stationary and flowing fluids, is given in review [1] and monograph [2]. Smaller reviews extending this list of cited papers and including some unlisted pioneering studies are given in recent articles of the authors [3, 4]. However, in the context of this study, the matter of particular current interest is a few works that take into account the axis misalignment of the inner and outer shells [5–9]. In [5], it was shown that the I http://www.gruppofrattura.it/VA/49/2452.mp4 S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 815 axial eccentricity affects the vibration process only in the case of a narrow annular gap between infinite shells. The analysis of the added liquid masses for two infinite cylinders with eccentricity, containing either an ideal or a viscous fluid in the gap between them, was performed in [6] and [7]. As noted in [7], the added masses and viscous damping coefficients, which grow with increasing eccentricity, are of greater practical interest, because they are widely used in studying hydroelastic vibrations and stability of structures in various engineering applications. The analysis of natural frequencies of cylindrical shells with an eccentricity immersed in a rigid cylindrical container was performed in works [8, 9]. The first of these papers considers the system, in which the non-viscous compressible fluid is contained only in the annular gap between the container and the shell, whereas the second paper examines the situation when the inner shell also contains a fluid. It has been demonstrated that with increasing eccentricity, the frequency of vibrations decreases, which is most significant for low circumferential and axial modes. Particular attention should also be paid to publications dealing with the analysis of coaxial shells partially filled with a fluid [10–15]. In all these papers, the studies of the dynamic behavior of coaxial shells are made with account of the free-surface sloshing effect, but do not include numerical simulation, which would allow one to estimate the effect of this boundary condition. Note that here we cite only those works which analyze the effect of partial filling of shells on the dynamic characteristics of the system. It is shown that an increase in the fluid level causes a decrease in the natural vibration frequencies. A more general review of papers aimed at studying hydroelastic vibrations of partially filled coaxial cylindrical shells is given in monograph [16]. In the above mentioned works, the solution of the problem is searched for in the context of a two-dimensional (flat) or axisymmetric problem formulation. In the case when non-coaxial structures are set in the vertical or horizontal position and partial filled with fluid, their symmetry with respect to the circumferential coordinate is disturbed, which necessitates the use of more complex spatial models. As far as the authors know, such studies, including those that take into account the elasticity of both shells, are not presented in the literature. Electroelastic materials embedded in or attached to engineering structures have been used for rather long time in various fields of technology in order to improve the performance characteristics of a product. Intellectual systems of passive or active control of dynamic behavior built on the basis of their peculiar properties serve to reduce the level of mechanical vibrations or acoustic noise. The bibliography of basic works dealing with the analysis of various approaches used in modeling the structures with piezoelectric elements, and presenting examples of their practical use can be found in monograph [17]. In [18], the authors come to a conclusion that successful practical application of piezoelectric elements coming in contact with a fluid or gaseous medium is possible only with an adequate description of their coupled response. In the case of thin-walled structures this can be achieved with the model proposed in [19]. The model is based on a single- layer representation of the elastic body displacements and a layered description of piezoelectric properties. This approach was used for modeling the plates and single shells, including those interacting with a quiescent fluid, and also the elastic structures with external patches made completely or partially from the electroelastic material [20–23]. For a linear distribution of electrical properties through the body thickness, this method, when used in the framework of finite element modeling, allows us to simplify the solution by eliminating the electrical components at the element level. In the literature, there are no papers describing the application of the above method to coaxial shells. Thus, the purpose of this work is to analyze the natural frequencies and the corresponding modes of vibrations of elastic and piezoelastic eccentric (non-coaxial) cylindrical shells at different values of the annular gap width and different levels of ideal compressible fluid contained in this gap. The solution of the problem is found using a modified version of the finite element algorithm developed in the previous works for evaluating the hydroelastic interaction of single shell structures with arbitrary cross section completely or partially filled with fluid [24–26]. STATEMENT OF THE PROBLEM AND CONSTITUTIVE RELATIONS e consider a structure consisting of two vertically located eccentric (non-coaxial) shells of length L, radiuses (1)R and ( 2)R , and thicknesses ( 1 )h and ( 2 )h (Fig. 1) made of electroelastic material (piezoceramics). Hereinafter, the superscripts “(1)” and “(2)” characterize the inner and outer shells, respectively. The axis of rotation of the inner shell is displaced laterally with respect to rotation axis of the outer shell by the value a  ( 2 ) ( 1)( )a R R , and the space between them is filled with an ideal compressible fluid to the level H  (0 )H L . It is necessary to investigate the natural frequencies and modes of vibrations of the system under different kinematic and electric boundary conditions, and geometric parameters (the size of the annular gap, the height of the fluid, eccentricity). W S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 816 Figure 1: Schematic representation of the computational domain. An ideal compressible fluid is considered in the framework of a potential theory, in which the wave equation in the case of small perturbations takes the following form [27]:     2 2 1 tc   , (1) where  is the velocity potential, c is the speed of sound in the liquid,  is the differential Nabla operator,  t is the time derivative. Using the perturbation velocity potential as an unknown function allows to take into account both the flow or the rotation of the inviscid compressible fluid with a slight modification of the numerical algorithm. In the case of complete filling of the shell with a fluid, the following boundary conditions are prescribed for the velocity potential:    0, : 0x L x . (2a) In the case of partial filling, it is assumed that the free surface of the fluid freeS remains stationary and is free from the dynamic pressure and surface tension. The corresponding boundary condition is given as [28]  : 0x H  . (2b) On the wetted surfaces  ( ) ( )i if sS S S the impermeability conditions hold true      ( ) ( ) i i w tn , (3) where ( )in is the vector of the normal to the shell surface; ( )iw is the normal component of shell displacements; fS and ( )i sS are the surfaces that bound the volumes of the fluid fV and shell ( )i sV . In what follows,  1, 2i . The hydrodynamic pressure p exerted on the elastic structures by the fluid is calculated using the Bernoulli equation: S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 817    ( )i fp t   , (4) where f is the fluid density and the sign of the right-hand parts of Eqns. (3)–(4) depends on the direction of the vectors of the normal to the external surfaces of the shells. Eqn. (1) together with the boundary conditions (2)–(3) are converted by the Bubnov – Galerkin method to a weak form [26]                   ( 1) ( 2 ) 2 (1) ( 2 ) 2 2 ˆ ˆ ˆ1ˆd d d d 0, 1, f f n n n n f V V S S w w F V F V F S F S n m t tc t , (5) where ̂ and ( )ˆ iw are the approximations of the velocity potential and the normal components of the vectors of shell displacements; nF and fm are the basis functions and their number, respectively. In the general case, the behavior of the electroelastic body is described by the equations of state, piezoelectric effect, and the electrostatic relations [29, 30]   Tcε e E , (6)  D eε dE , (7)   0D , (8)  E . (9) Here: , ε are the stress and linear strain tensors; E and D are the vectors of electric field intensity and electrical induction; c, e, d are the matrices of elastic constants, piezoelectric and dielectric coefficients;  is the electrostatic potential. In the case when a thin-walled structure with a radially polarized field is in the plane stress state, Eqns. (6)–(9) can be simplified. In particular, of all components the vectors of the electric field strength and nonzero electrical induction retain only 3E and 3D , and Eqns. (6) and (7) are written as follows [31, 32]:                                                            11 11 31 11 12 22 22 32 21 22 12 12 3 66 0 0 0 0 0 0 0 , 0 0 0 0 0 0 e c c e c c E c c c , (10)                                            11 22 3 31 32 12 33 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0D e e d E    , (11) where           23 3 33 66 66 3 3 33 3 33 33 33 33 33, , , , , 1, 2lm lm l m m m mc c c c c c c e e e c c d d e c l m (12) Here we proceed from the assumption that the shell surface is covered with thin weightless electrodes, which can either be shorted, which corresponds to the boundary condition 3 0E (short circuited) or disconnected 3( 0,D open circuited). From Eqns. (6)–(9) the following integral relation can be obtained [21] S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 818        T T d d 0 s sV V V VE eε E dE . (13) In the shell simulation, it is assumed that their curved surfaces can be represented with a reasonable degree of accuracy as a set of flat segments [33]. The strains in each segment are determined according to the classical theory of thin plates [34]. The corresponding relations in the Cartesian coordinates  , ,x y z , which are tied in with the lateral surfaces of thin- walled structures, are written as                                                     T( ) ( ) ( ) ( ) ( ) ( ) 11 22 12 T( ) ( ) ( ) ( ) T( ) ( ) ( ) ( ) 11 22 12 T2 ( ) 2 ( ) 2 ( ) T( ) ( ) ( ) ( ) 11 22 12 2 2 , , , , , , , , , , , , 2 . i i i i i i i i i i i i i i i i i i i i i z u v u v x y y x w w w x yx y ε ε ε (14) Hereinafter, the over-line means that the quantity is defined in the coordinates  , ,x y z , and ( ) ( ), i iu v and ( )iw is the meridional, circumferential and normal components of the displacement vectors of the middle surfaces of the shells in these coordinates. Taking into account the accepted simplifications, the physical relations interconnecting the vectors of generalized forces and moments ( )iT , generalized deformations ( )i and electric field strength ( )iE , can be written in the matrix form as                          T( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 11 22 12 11 22 12 ( ) ( ) T( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 11 22 12 11 22 12( ) ( ) , , , , , , , , , , , , , , i i i i i i i i i i i i i i i i i i i i i i i i T T T M M MT D ε G E T A B D ε B C (15) where the coefficients entering into the stiffness matrix  ( )iD are calculated as          ( ) ( ) 2 ( ) ( ) ( ) 2 ( ) 2 , , 1, , d i i h i i i i h z z zA B C c , (16) and the structure of the matrix ( )iG will be presented below. A mathematical formulation of the problem of the dynamics of elastic shells is based on the variational principle of virtual displacements, which, with account of relation (4) and the work made by the inertial forces, can be written in the matrix form as                           ( ) ( ) ( ) T T T( ) ( ) ( ) ( ) ( ) ( ) ( ) T T( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) , , , , , , , 0, 0, , 0, 0, 0 , i i i s s i i i i i i i s V V S i i i i i i i i i x y z dV dV dS u v w p ε T u u u P 0 u P (17) where ( )is are the densities of shell materials, ( )iu are the vectors of the displacements and rotation angles of shells, and ( )iP are the vectors of loads acting on the shell surfaces. S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 819 NUMERICAL IMPLEMENTATION ollowing [21], let us decompose the shells through its thickness into N layers and represent the component of the electric field 3E for each layer as    ( ) ( ) ( )3 i i i k kk E V h , (18) where   ( ) ( ) ( ) 1 i i i k k kh z z is the thickness of the k -th layer;    ( ) ( ) ( ) ( )( 2 2)i i i iz h z h is the coordinate measured from the middle surface of the shell;   ( ) ( ) ( ) 1 i i i k k kV   is the difference in the electrostatic potentials between the upper and lower surfaces of the layer, which, as well as the components of the vectors of shell displacements ( )iu and the velocity potential  , turns to be an unknown variable. The numerical solution of the problem is found using the finite element method (FEM) [33]. To describe the velocity potential ˆ, the basic functions nF , and membrane displacements of shells ( ) ( )( , )i iu v we use the Lagrange shape functions with linear approximation, and for bending deflections of shells ( )iw we use the Hermite non-conformal shape functions The discretization of the fluid and shell regions is made using spatial prismatic (8-node brick) and flat quadrangular (4-node plane) finite elements, respectively. Expressing the unknown variables in terms of their nodal values, we arrive at the following matrix relations     ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ), , ,i i i i i i i i ie e e e eu N u F ε B u E B Φ  , (19) where F and ( )iN are the shape functions for the velocity potential of the fluid and the vectors of the nodal displacements of the shell, ( )ieu and e are the vectors of the nodal values, ( )iB are the gradient matrices which are determined by Eqn. (14) and relate the strain vectors to the shell displacements,     T( ) ( ) ( ) ( )1 , , ,i i i ie k NE E EE     T( ) ( ) ( ) ( )1 , , ,i i i ie k NV V VΦ    ( ) ( ) ( ) ( )1diag 1 , ,1 , 1 .i i i ik Nh h hB Using the introduced relations the matrix ( )iG in Eqn. (15) is generated in the following way                                     ( ) ( ) ( ) 11 1 1 ( ) ( ) ( ) 21 2 2 ( ) ( ) ( ) ( ) 11 1 1 ( ) ( ) ( ) 21 2 2 ˆ ˆ ˆ ˆ ˆ ˆ 0 0 0 0 0 0 i i i k N i i i k N i i i i k N i i i k N G G G G G G G G G G G G G , (20) where        2 2( ) ( ) ( ) ( ) ( ) ( ) ( )3 1 31ˆ , , 1, , 1, 2. 2 i i i i i i i lk k l lk k k lk k G h e G z z e k N l           (21) From Eqns. (5), (13) and (17) taking into account (18)–(20), we derive with the aid of the standard FEM operations a coupled system of equations, which can be used to describe the interaction of electroelastic shells with an ideal compressible fluid. The system can be represented in a matrix form as F S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 820                                                            (1) (1)(1)(1) (1) (1) ( 2 )( 2 ) ( 2 )( 2 ) ( 2 ) ( 2 ) (1) ( 2 ) 0 0 0 00 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 sfs sfs f fs fs u uCM Cu uM M C C                                          (1) (1) (1) T(1) (1) (1) ( 2 ) ( 2 ) ( 2 ) T ( 2 )(1) ( 2 ) 0 0 0 0 0 0 0,0 0 0 0 0 0 0 0 0 0 s s s s s f       K K u K K K K u K K K    (22) where ( ) , iu  and ( )i are the generalized vectors of the shell displacements and rotation angles, the perturbation velocity potential of the fluid and the difference in the electrostatic potentials of the shells. The typical stiffness (K), mass (M) and damping (C) matrices for individual finite elements are generated as follows:         ( ) ( ) ( ) ( ) ( ) T T T T ( ) T ( ) 2 T T T( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) T( ) ( ) ( ) ( ) ( 1 d , d , d , d , d , d , d , i f f i i i s s i s i i f f fs w V V S i i i i i i i i i i s s s sf f w V V S i i i i i s V V V S x x y y z z c V V S V                                        F F F F F F K M F F C F N K B D B M N N C N F K B G B K           ( ) T) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 1 33 33 331 d , diag , , , , . i s i i i V i i i i i i i k Nk N V h d h d h d     B H B H     (23) Here ( )iwN are the shape functions for the normal component of the nodal displacement vector of the shells. Since some matrices in (23) are written in the Cartesian coordinates  , ,x y z , the derivation of the solution in the global coordinates, which have been used to specify the problem geometry and to write the constitutive equations for the fluid, involves the transformation of the displacements and rotation angles using the matrix of the directional cosines γ of size 3 3 in the following way:        ( ) ( )0 0 i iu u   . (24) During the numerical implementation, at each node of the finite element of the shell it is necessary to determine six unknowns (three displacements and three angles of rotation). Therefore, the transformation of conventional matrices is performed using the diagonal matrix         diag , , , , , ,L :              T T T( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) T( ) ( ) ( ) ( ) ( ) ( ) , , , , . i i i i i i i i i i i s s s s sf sf i i i i i i s s fs fs  K L K L M L M L C L C K L K C C L (25) Expressing the quantities associated with the electrostatic potential in terms of the shell displacements S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 821      1 T( ) ( ) ( ) ( )i i i is K K u , (26) yields, instead of (22), the system of equations with respect to the displacements of electroelastic bodies and the velocity potential of the fluid                                                                     (1) (1) (1) (1) ( 2 ) ( 2 ) ( 2 ) ( 2 ) (1) ( 2 ) 1 T(1) (1) (1) (1) 1 T( 2 ) ( 2 ) ( 2 ) ( 2 ) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 s sf s sf f fs fs s s s s s s f       M u C u M u C u M C C K K K K K K K K K               (1) ( 2 ) 0, u u  (27) which can be written in a more compact form       T(1) ( 2 )0, , ,Mx Cx Kx x u u  . (28) Representing the motion of shells and fluid in the exponential form     ( ) ( ), ( , , ), ( , , ) exp( )i i x y z x y z tu u  , (29) we obtain instead of (28) the following expression   2 0 Mx Cx Kx . (30) Here,  ( )iu and  ( )i are some functions of the coordinates;   i   is the characteristic coefficient, in which  is the natural frequency of vibrations,  is the quantity responsible for the system damping, and  i 1 is the imaginary unit. The system of Eqns. (30) is reduced to an asymmetric generalized eigenvalue problem                     0 0 0 0   C K M x I I x , (31) where I is the unit matrix. The calculation of the eigenvalues of the system of Eqns. (31) is carried out using the ARPACK procedures, which are based on the implicitly restarted Arnoldi method [35]. RESULTS few numerical examples were considered to investigated the behavior of vertically oriented coaxial cylindrical shells, containing an ideal compressible fluid in the annular gap between them. The system parameters are listed in Tab. 1. The shells were made of PZT-5H piezoceramics, the physical and mechanical characteristics of which are given in Tab. 2. As the boundary conditions for thin-walled bodies, we chose rigid clamping        ( 0)x y zu v w at one end (cantilever, CF) or at both ends (clamped-clamped, CC). The calculations were A S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 822 performed at different values of the dimensionless annular gap k, filling level  and displacement (eccentricity) of the inner shell         ( 2 ) ( 1) ( 1) ( 2 ) ( 1), ,k R R R H L a R R  . (32) Note that the radius of the outer shell ( 2 )( )R is fixed, and the radius of the inner shell (1)( )R varies and depends on the parameter k. Geometric data Fluid properties Parameter Value Parameter Value L, m 1 f , kg/m3 1000 ( 2 )R , m 0.1 c, m/s 1500  (1) ( 2 )h h h , m  45 10 Table 1: Geometrical parameters of the shells and physical properties of the fluid. Parameter Value Parameter Value c11, GPa 126 e15, C/m2 17 c12, GPa 79.5 e24, C/m2 17 c13, GPa 84.1 e31, C/m2 -6.5 c22, GPa 126 e32, C/m2 -6.5 c23, GPa 84.1 e33, C/m2 23.3 c33, GPa 117 d11  1010 , F/m 150.3 c44, GPa 23 d22  1010 , F/m 150.3 c55, GPa 23 d33  1010 , F/m 130 c66, GPa 23.3 s , kg/m3 7500 Table 2: Properties of piezoceramics PZT-5H. Verification of the numerical model The assessment of reliability of the constructed algorithm for the numerical investigation of an electroelastic thin-walled body was made by comparing its performance with the results of work [21]. Here the simply supported at both ends  ( 0)v w shell made of PZT-5H piezoceramics was investigated (Tab. 3). The boundary conditions for the fluid were specified as follows: 0  at  0x and  .x L Geometric data Fluid properties Parameter Value Parameter Value L, m 5 f , kg/m3 1000 R, m 1 c, m/s 1500 h, m 0.02 Table 3: Geometrical parameters of electroelastic shells and physical properties of the fluid [21]. A comparison between the natural frequencies of the empty shell and the shell filled with an ideal compressible fluid is given in Tab. 4. The number of half-waves in the circumferential and meridional directions is denoted by j and m, respectively. According to the data presented in the table, there is good agreement between the results obtained by analytical and numerical solution of the problem. The assessment of the reliability of the results obtained for coaxial shells  ( 0) made of isotropic elastic material and containing fluid in the gap between them was carried out by comparison with the results of analytical solution presented in [36]. In this case, the matrices ( )iK and ( )i sK in Eqn. (27) are excluded from the calculation, and the stiffness matrices S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 823  ( )iD are generated in the known fashion using Young's modulus sE and Poisson’s ratio  . Geometrical and physicomechanical parameters of the problem are summarized in Tab. 5, and a comparison of the obtained results is given in Tab. 6. Without fluid With fluid Wave numbers Short circuited Open circuited Wave numbers Short circuited Open circuited j m [21] Present [21] Present j m [21] Present [21] Present 1 1 88.2310 88.1261 93.7900 93.7266 1 1 44.7600 44.6100 48.2710 48.1062 2 216.260 217.239 228.227 229.386 2 107.390 107.797 115.614 115.968 0 275.956 280.586 275.956 280.586 3 154.131 154.925 166.065 166.705 3 311.647 313.207 330.649 332.755 4 189.327 190.099 204.211 204.721 4 366.968 368.502 392.134 393.952 5 217.190 217.515 234.469 234.471 2 1 37.3160 36.8135 40.1090 39.6043 2 1 19.9160 19.5781 21.4640 21.3157 2 114.210 114.107 122.677 122.616 2 62.5970 62.3712 67.6230 67.4905 3 193.791 194.308 207.916 208.565 3 108.575 108.650 117.248 117.433 4 259.346 260.087 278.434 279.316 4 149.022 149.255 160.912 161.282 5 308.464 308.553 331.551 331.658 5 182.817 182.814 197.400 197.590 Table 4: A comparison of the natural vibration frequencies  (Hz) of a simply supported shell with or without fluid at different combinations of electrical boundary conditions. Geometric data Material properties Fluid properties Parameter Value Parameter Value Parameter Value L, m 0.3  (1) ( 2 )s s sE E E , GPa 69 f , km/m 3 1000 (1)R , m 0.1  (1) ( 2 )   0.3 c, m/s 1483 ( 2 )R , m 0.15   (1) ( 2 ) s s s   , kg/m3 2700  ( 1) ( 2 )h h h , m 0.002 Table 5: Physicomechanical and geometrical parameters of the system, consisting of elastic coaxial shells and a fluid. j m [36] Phase mode Present j m [36] Phase mode Present 1 1 391.1 out-of-phase 390.3 2 1 435.6 out-of-phase 434.9 2 846.7 out-of-phase 846.3 2 907.1 out-of-phase 905.2 3 1397.5 out-of-phase 1396.5 1 996.8 in-phase 995.2 1 1736.6 in-phase 1732.5 3 1401.3 out-of-phase 1399.7 4 1908.5 out-of-phase 1907.3 * 1822.2 mixed phase 1819.1 5 2317.2 out-of-phase 2312.5 * 1892.6 mixed phase 1891.2 * 2623.4 mixed phase 2579.5 * 2265.3 mixed phase 2262.2 3 1 403.0 out-of-phase 403.2 4 1 382.5 out-of-phase 383.3 1 671.3 in-phase 671.4 1 561.9 in-phase 562.9 2 858.30 out-of-phase 857.7 2 791.0 out-of-phase 792.0 2 1344.8 in-phase 1343.8 2 1075.5 in-phase 1076.7 3 1352.4 out-of-phase 1351.5 3 1267.5 out-of-phase 1268.8 4 1810.7 out-of-phase 1811.1 3 1676.9 in-phase 1678.0 3 2010.6 in-phase 2008.6 4 1729.2 out-of-phase 1731.3 Table 6: A comparison of the natural vibration frequencies  (Hz) of clamped-clamped coaxial shells containing an ideal compressible fluid in the annular gap between them. Note that in [36], it was first shown that for two elastic coaxial shells, along with the in-phase (the direction and number of meridional half-waves m coincide) and out-of-phase (opposite directions) vibration modes there also exist mixed modes (the number of half-waves varies). From the data presented in Tab. 6, we may conclude that the results obtained in the framework of the model used in this paper agree fairly well with the data of the analytical solution presented in [36]. S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 824 Natural frequencies of electroelastic shells Fig. 2 shows the graphs of the lowest natural frequencies of vibrations  as a function of the dimensionless displacement of the inner shell  at different filling levels  for the gap  1 10.k The results were obtained for two configurations of electroelastic bodies under different boundary conditions. As is evident from the presented data, in the presence of eccentricity the frequencies decrease regardless of the choice of electrical and kinematic boundary conditions at the edges of the shells. According to Refs. [6, 7, 9], the reason for such a behavior is the growth of the added fluid mass. In this case, a decrease in the level of a fluid and, consequently, a decrease in its the added mass contributes to a stronger dependence of the lowest vibration frequency on the axial misalignment of the shells and eventually leads to their greater relative decrease. The cantilevered structures are less sensitive to the eccentricity of the inner shell, so that over a fairly wide range of the parameter  a decrease in the frequencies is inessential. Similar dependencies for other values of the annular gap k do not demonstrate qualitative differences, which is supported by the data given in Tab. 7. (a) (b) Figure 2: The dependences of the lowest vibration frequencies  (Hz) on the eccentricity  obtained at different filling levels and electrical boundary conditions (shaded symbols — short circuit mode, unshaded symbols — open circuit mode) for clamped-clamped (a) and cantilevered (b) shells,  1 10 .k Due to an additional contribution of the associated stiffness, the electrical boundary conditions corresponding to the open-circuited configuration, lead to an increase in frequencies compared to the case of short-circuited configuration. The relative differences between these two variants ( ) as a function of eccentricity  are presented in Figs. 3–5 for different values of the annular gap k. It can be seen from the above data that the curves are of a qualitatively different character, which is associated with the boundary conditions at the edges of the shells, the size of the annular gap and the level of filling with a liquid. The higher is the fluid level in the annular gap between the shells, the larger is the relative difference in frequencies  . Fig. 6 shows the dependences, which allow us to estimate the variation in the lowest natural frequencies of vibrations with the level of the fluid in the gap at  1 10k and various kinematic boundary conditions. The curves presented in the figure reflect the growth of frequencies with a decrease in the filling level. It can be seen that they are of qualitatively similar character for shells with different boundary conditions. However, in the case of rigid clamping at both ends (Fig. 6a), there is a range of the fluid levels, at which the frequencies practically do not change regardless of the displacement of the inner shell . For cantilevered shells (Fig. 6b), such a behavior is not observed. In this case, the differences in the frequencies obtained at different electrical boundary conditions are insignificant. Figs. 7, 8 show the mode shapes of the shells ( 1 10 ,k  open circuited, CC) obtained for different variants of the inner shell displacement and the fluid level in the annular gap. In the figures, displaying the transverse and longitudinal sections of the system, the dotted lines represent the shells in the undeformed state, and the solid line — the shells in the deformed state. The liquid is shown by light grey. 0.0 0.2 0.4 0.6 0.8  0 10 20 30 40  ●, ○  = 0.25 ▲, Δ  = 0.50 ■, □  = 1.00 0.0 0.2 0.4 0.6 0.8  0 12 24 36 48  ●, ○  = 0.25 ▲, Δ  = 0.50 ■, □  = 1.00 S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 825 k     Clamped-Clamped shells Clamped-Free shells Short circuited Open circuited  , % Short circuited Open circuited  , % 1/2 1.00 0.00 32.7722 36.6289 11.7682 8.58148 9.65366 12.4942 0.50 28.0167 30.9768 10.5657 7.78148 8.71879 12.0455 0.95 15.9592 17.8405 11.7885 5.22212 5.88609 12.7145 0.50 0.00 44.0368 49.0209 11.3179 28.3814 31.4965 10.9758 0.50 38.5264 42.7049 10.8459 26.7499 29.7033 11.0407 0.95 21.5175 24.0634 11.8313 17.2344 19.2967 11.9662 0.25 0.00 89.2623 98.0454 9.83972 39.1212 43.3674 10.8539 0.50 81.4270 90.0862 10.6343 39.0705 43.3060 10.8407 0.95 45.0888 49.9111 10.6952 37.8442 41.8567 10.6027 1/10 1.00 0.00 15.4716 17.2113 11.2443 3.81831 4.25959 11.5569 0.50 13.3385 14.7800 10.8069 3.47482 3.87567 11.5359 0.95 7.13355 7.97309 11.7689 2.24936 2.52362 12.1930 0.50 0.00 21.7425 24.0692 10.7013 15.8505 17.4368 10.0074 0.50 18.7429 20.7607 10.7661 14.1066 15.5642 10.3326 0.95 9.70433 10.8434 11.7380 7.85436 8.76154 11.5500 0.25 0.00 53.0114 57.4618 8.39530 36.5821 40.5511 10.8496 0.50 43.6699 47.9169 9.72512 35.7320 39.4371 10.3692 0.95 20.7597 22.9591 10.5944 20.0218 22.0848 10.3036 1/100 1.00 0.00 4.97846 5.50687 10.6138 1.22342 1.35737 10.9483 0.50 4.27499 4.72525 10.5323 1.11144 1.23392 11.0193 0.95 2.25033 2.51305 11.6749 0.71178 0.79643 11.8917 0.50 0.00 7.08603 7.80343 10.1242 5.42042 5.92045 9.22498 0.50 6.04697 6.68328 10.5228 4.69234 5.15888 9.94261 0.95 3.06832 3.42595 11.6555 2.50227 2.78760 11.4030 Table 7: Natural vibration frequencies  (Hz) of the system of coaxial shells at different values of the annular gap k, eccentricity , and filling level  and varying electrical and kinematic boundary conditions. (a) (b) Figure 3: The relative difference in the frequencies  as a function the eccentricity  at different levels of the fluid in the annular gap  for clamped-clamped (a) and cantilevered shells (b),  1 10 0.k 0.0 0.2 0.4 0.6 0.8  8.0 8.8 9.6 10.4 11.2 , %  = 1.00  = 0.50  = 0.25 0.0 0.2 0.4 0.6 0.8  7.5 8.4 9.3 10.2 11.1 , %  = 1.00  = 0.50  = 0.25 S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 826 (a) (b) Figure 4: The relative difference in the frequencies  as a function the eccentricity  at different levels of the fluid in the annular gap  for clamped-clamped (a) and cantilevered shells (b),  1 10 .k (a) (b) Figure 5: The relative difference in the frequencies  as a function the eccentricity  at different levels of the fluid in the annular gap  for clamped-clamped (a) and cantilevered shells (b),  1 2 .k From the presented data it is clear that the amplitude of the displacements on the wetted areas of the shell surfaces exceeds the displacement amplitude on the areas not contacting with the fluid. A decrease in the level of the fluid also leads to a decrease in the maximum amplitudes of the displacements occurring in the shells. Moreover, in the case of the eccentricity, the displacement maximum occurs in the areas where the distance between the thin-walled bodies is minimal  ( 2 ) (1)( )R R a and gradually decreases towards the region where the spacing between them is maximum  ( 2 ) (1)( )R R a . The behavior described above is the result of a change in the added mass of the fluid in the corresponding regions of the shells. 0.0 0.2 0.4 0.6 0.8  10.5 11.0 11.5 12.0 12.5 , %  = 0.50  = 1.00  = 0.25 0.0 0.2 0.4 0.6 0.8  9.5 10.0 10.5 11.0 11.5 , %  = 0.50  = 1.00  = 0.25 0.0 0.2 0.4 0.6 0.8  9.4 10.0 10.6 11.2 11.8 , %  = 1.00  = 0.25  = 0.50 0.0 0.2 0.4 0.6 0.8  8.0 8.8 9.6 10.4 11.2 , %  = 1.00  = 0.50  = 0.25 S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 827 (a) (b) Figure 6: Variation in the lowest natural frequencies  (Hz) with the level of the fluid in the annular gap between the shells  at different values of the eccentricity  and various electrical boundary conditions (shaded symbols — short-circuit mode, unshaded symbols — open-circuit mode) for clamped-clamped (a) and cantilevered (b) shells: k = 1/10. (a) (b) (c) Figure 7: The lowest vibration mode of the electroelastic shell interacting with a fluid in the transverse ( 2)x L and longitudinal ( 0)y sections at different levels of the fluid in the annular gap  (CC, open-circuit mode,  1 10, 0)k  : (a)  1.0; (b)  0.5; (с)  0.25. 0.0 0.2 0.4 0.6 0.8  0 10 20 30 40  ●, ○  = 0.00 ▲, Δ  = 0.50 ■, □  = 0.95 0.0 0.2 0.4 0.6 0.8  0 28 56 84 112  ●, ○  = 0.00 ▲, Δ  = 0.50 ■, □  = 0.95 S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 828 (a) (b) (c) Figure 8: The lowest vibration mode of the electroelastic shell interacting with a fluid in the transverse ( 2)x L and longitudinal ( 0)y sections at different levels of the fluid in the annular gap  (CC, open-circuit mode,  1 10 , 0.5)k  : (a)  1.0; (b)  0.5; (с)  0.25. CONCLUSION he dynamic behavior of eccentric (non-coaxial) electroelastic cylindrical shells, containing a quiescent fluid in the annular gap between them, was studied in a three-dimensional formulation with the use of the proposed mathematical model, the numerical implementation of which was based on the finite element method. The analysis of the obtained results allowed us to estimate the influence of the size of the annular gap between the shell, the level of the fluid in it, the eccentricity of the inner shell, and the electrical and kinematic boundary conditions on the lowest natural frequencies of vibrations and the corresponding mode shapes. It was demonstrated that both an increase in the misalignment of the shells and the fluid level in the annular gap and a decrease in the size of the gap itself are responsible for a decay of the lowest vibration frequencies. The growth of the lowest frequencies is provided by setting on the electrode covered surfaces of the shells the electrical boundary condition peculiar to the open-circuited configuration. ACKNOWLEDGMENTS he study was supported by the grant of the Russian Scientific Foundation (project No. 18-71-10054). REFERENCES [1] Brown, S.J. (1982). A survey of studies into the hydrodynamic response of fluid-coupled circular cylinders, J. Press. Ves. Tech., 104(1), pp. 2-19. DOI: 10.1115/1.3264181 T T S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 829 [2] Païdoussis, M.P. (2016). Fluid-Structure Interactions: Slender Structures and Axial Flow, vol. 2, 2nd edn., London, Elsevier Academic Press. [3] Bochkarev, S.A., Lekomtsev, S.V. and Senin, A.N. (2018). Analysis of spatial vibrations of coaxial cylindrical shells partially filled with a fluid, Comput. Continuum Mech., 11(4), pp. 448-462. DOI: 10.7242/1999-6691/2018.11.4.35 [4] Bochkarev, S.A., Lekomtsev, S.V., Matveenko, V.P. and Senin, A.N. (2019). Hydroelastic stability of partially filled coaxial cylindrical shells, Acta Mech. (in press) [5] Buivol, B.N. and Guz, A.N. (1966). Oscillations of two cylindrical eccentrically arranged shells in a stream of inviscid liquid, Rep. AS Ukr. SSR, no. 11, pp. 1412-1415. [6] Chung, H. and Chen, S.-S. (1977). Vibration of a group of circular cylinders in a confined fluid, J. Appl. Mech., 44(2), pp. 213–217. DOI: 10.1115/1.3424026 [7] Mateescu, D., Païdoussis, M.P. and Sim, W.-G. (1994). Spectral solutions for unsteady annular flows between eccentric cylinders induced by transverse oscillations, J. Sound Vib., 177(5), pp. 635–649. DOI: 10.1006/jsvi.1994.1458. [8] Jeong, K.-H. (1999). Dynamics of a concentrically or eccentrically submerged circular cylindrical shell in a fluid-filled container, J. Sound Vib., 224(4), pp. 709-732. DOI: 10.1006/jsvi.1999.2209. [9] Jeong, K.-H., Lee G.-M. and Chang, M.-H. (2001). Free vibration analysis of a cylindrical shell eccentrically coupled with a fluid-filled vessel, Comput. Struct., 79(16), pp. 1517-1524. DOI: 10.1016/S0045-7949(01)00031-1. [10] Baghdasaryan, G.Ye. (1968). The vibrations of coaxial cylindrical shells with a clearance partially filled with fluid, Mechanics. Proc. AS Arm. SSR, 21(4), pp. 40-47. [11] Tani, J. and Haiji, H. (1986). The free vibration of free-clamped fluid-coupled coaxial cylindrical shells, Trans. Jpn. Soc. Mech. Eng. Part C, 52(484), pp. 3137-3144. DOI: 10.1299/kikaic.52.3137. [12] Tani, J., Sakai, T. and Chiba, M. (1989). Dynamic stability of fluid-coupled coaxial cylindrical shells under vertical excitation, Trans. Jpn. Soc. Mech. Eng. Part C, 55(512), pp. 870-877. DOI: 10.1299/kikaic.55.870. [13] Okazaki, K., Tani, J. and Sugano, M. (2002). Free vibrations of a laminated composite coaxial circular cylindrical shell partially filled with liquid, Trans. Jpn. Soc. Mech. Eng. Part C, 68(671), pp. 1942–1949. DOI: 10.1299/kikaic.68.1942. [14] Gavrilova, E. (2003). A study of the vibration of fluid coupled coaxial cylindrical shell, In: Benaroya H., Wei T.J. (eds) IUTAM Symposium on Integrated Modeling of Fully Coupled Fluid Structure Interactions Using Analysis, Computations and Experiments. Fluid Mechanics and its Applications, 75, pp. 363-374. DOI: 10.1007/978-94-007-0995-9_26. [15] Baghdasaryan, G.Ye. and Marukhyan, S.A. (2011). Dynamic behavior of a coaxial cylindrical shells, with a gap partially filled with fluid, Mechanics. Proc. NAS Arm., 64(3), pp. 10-21. [16] Ibrahim, R.A. (2005). Liquid Sloshing Dynamics: Theory and Applications, Cambridge, Cambridge University Press. [17] Abramovich, H. (2016). Intelligent Materials and Structures, Berlin, De Gruyter. [18] Kaljevic, I. and Saravanos, D.A. (1997). Steady-state response of acoustic cavities bounded by piezoelectric composite shell structures, J. Sound Vib., 204(3), pp. 459–476. DOI: 10.1006/jsvi.1996.0911 [19] Saravanos, D.A. (1997). Mixed laminate theory and finite element for smart piezoelectric composite shell structures, AIAA J., 35(8), pp. 1327–1333. DOI: 10.2514/2.264. [20] Lammering, R. and Mesecke-Rischmann, S. (2003). Multi-field variational formulations and related finite elements for piezoelectric shells, Smart Mater. Struct., 12(6), pp. 904–913. DOI: 10.1088/0964-1726/12/6/007. [21] Larbi, W., Deü, J.-F. and Ohayon, R. (2007). Vibration of axisymmetric composite piezoelectric shells coupled with internal fluid, Int. J. Num. Meth. Eng., 71(12), pp. 1412–1435. DOI: 10.1002/nme.1987. [22] Deü, J.-F., Larbi, W. and Ohayon, R. (2008). Piezoelectric structural acoustic problems: Symmetric variational formulations and finite element results, Comp. Meth. Appl. Mech. Eng., 197(19–20), pp. 1715–1724. DOI: 10.1016/j.cma.2007.04.014. [23] Larbi, W., Deü, J.-F. and Ohayon, R. (2012). Finite element formulation of smart piezoelectric composite plates coupled with acoustic fluid, Compos. Struct., 94(2), pp. 501–509. DOI: 10.1016/j.compstruct.2011.08.010. [24] Bochkarev, S.A., Lekomtsev, S.V. and Matveenko, V.P. (2013). Numerical modeling of spatial vibrations of cylindrical shells, partially filled with fluid, Comput. technologies, 18(2), pp. 12-24. [25] Bochkarev, S.A., Lekomtsev, S.V. and Matveenko, V.P. (2015). Natural vibrations of loaded noncircular cylindrical shells containing a quiescent fluid, Thin-Walled Struct., 90, pp. 12-22. DOI: 10.1016/j.tws.2015.01.001 [26] Bochkarev, S.A., Lekomtsev S.V. and Matveenko, V.P. (2016). Natural vibrations and stability of elliptical cylindrical shells containing fluid, Int. J. Struct. Stab. Dyn., 16(10), 1550076. DOI: 10.1142/S0219455415500765. [27] Vol’mir, A.S. (1979). Shells in Liquid and Gas Flow. Hydroelastisity problems, Moscow, Nauka. [28] Amabili, M. (1996). Free vibration of partially filled, horizontal cylindrical shells, J. Sound Vib., 191(5), pp. 757-780. S.A. Bochkarev et alii, Frattura ed Integrità Strutturale, 49 (2019) 814-830; DOI: 10.3221/IGF-ESIS.49.15 830 DOI: 10.1006/jsvi.1996.0154. [29] Grinchenko, V.T., Ulitko, A.F. and Shul’ga, N.A. (1989). Mechanics of Coupled Fields in Elements of Structures. vol. 5. Electroelasticity, Kiev, Naukova Dumka. [30] Rogacheva, N.N. (1994). The Theory of Piezoelectric Shells and Plates, London, CRC Press. [31] Sheng, G.G. and Wang, X. (2010). Thermoelastic vibration and buckling analysis of functionally graded piezoelectric cylindrical shells, Appl. Math. Model., 34(9), pp. 2630-2643. DOI: 10.1016/j.apm.2009.11.024. [32] Yao, G. and Li, F.-M. (2014). The stability analysis and active control of a composite laminated open cylindrical shell in subsonic airflow. J. Intel. Mat. Sys. Struct., 25(3), pp. 259-270. DOI: 10.1177/1045389X13491020. [33] Zienkiewicz, O.C. (1971). Finite Element Method in Engineering Science, New York, McGraw-Hill. [34] Reddy, J.N. (2015). An introduction to nonlinear finite element analysis, 2nd edn, Oxford, Oxford University Press. [35] Lehoucq, R.B. and Sorensen, D.C. (1996). Deflation techniques for an implicitly restarted Arnoldi iteration, SIAM J. Matrix Anal. Appl., 17(4), pp. 789–821. DOI: 10.1137/S0895479895281484. [36] Jeong, K.-H. (1998). Natural frequencies and mode shapes of two coaxial cylindrical shells coupled with bounded fluid, J. Sound Vib., 215(1), pp. 105–124. DOI: 10.1006/jsvi.1998.1648. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice