Free vibration analysis of elastic plate structures by boundary element method

  • Published on
    15-Jul-2016

  • View
    215

  • Download
    3

Embed Size (px)

Transcript

<ul><li><p>Free vibration analysis of elastic plate structures by boundary element method </p><p>MASATAKA TANAKA </p><p>Department of Mechanical Engineering, Shinshu University, 500 Wakasato, Nagano 380, Japan </p><p>KOUJI YAMAGIWA </p><p>Graduate School of Shinshu University </p><p>KENICHI MIYAZAKI </p><p>Seiko-Epson Co., Ltd </p><p>and TAKAHIRO UEDA </p><p>Tateishi Electric Co., Ltd </p><p>This paper is concerned with an integral equation approach to the free vibration problems of elastic plate structures. In the usual boundary integral equation method, the eigen frequency must be determined by means of the so called direct search of the zero-determinant value of the system matrix. To circumvent these difficulties this paper presents a new integral equation approach and its solution procedure, in which an approximate fundamental solution to the static problem is used for the formulation. The resulting set of integral equations are discretized by means of the boundary-domain element method and reduced to the system of algebraic eigenvalue equations. The potential usefulness of the proposed method is demonstrated through some sample computations. </p><p>Key Words: free vibration analysis, plate structure, bending vibration, in-plane vibration, boundary integral equation, boundary-domain element method. </p><p>1. INTRODUCTION </p><p>The plate structures are widely used in various fields because they have an advantage of larger load-carrying capacity in spite of light weight 1-3. It is one of the most important subjects in engineering to develop an effective solution procedure for the free vibration problems of plate structures. Much knowledge has been already accumulated for both the in-plane and out-of-plane vibration problems of single plates with arbitrary shapes through the application of finite element, finite difference or other analytical methods of solution TM. Recently, the boundary element methods have been successfully applied to these problems 5-1 o. However, there have still been few investigations for the free vibration problem of plate structures composed of elastic plates. In this study, we first discuss the in-plane and out-of-plane vibration problems of a single plate, and then propose a new integral equation formulation and its solution procedure by means of the boundary-domain element method. </p><p>As has been well known 5-7, the exact fundamental solutions of the problem under consideration are expressed in terms of the Hankel and Bessel functions. This would often lead to some difficulties in manipulation of the functions and numerical implementation of the </p><p>Accepted December 1987. Discussion closes May 1989. </p><p>formulation. Much care should be paid to circumvent these difficulties 5 6. In the present method of solution, approximate fundamental solutions of the corresponding static problems are used for the formulation. The resulting boundary integral equations include unknown variables in the inner domain of the plate component, and </p><p>Corner Point Ak+l y" </p><p>Boundary /-' I Ak-1 Y) Fig. 1. Notation and coordinate system </p><p>182 Engineering Analysis, 1988, Vol. 5, No. 4 1988 Computational Mechanics Publications </p></li><li><p>Free vibration analysis of elastic plate structures: M. Tanaka et al. </p><p>hence they are solved together with the additional integral equations by means of the boundary-dorrlain element method proposed previously by the senior author ~1. The powerful usefulness of the proposed method is demonstrated through some sample computations. It should be noted that similar ideas to that of the present paper are also successfully applied to the free vibration problem of single plates s. </p><p>2. FREE VIBRATION OF SINGLE PLATE </p><p>First, we consider the out-of-plane bending vibration of a single plate. We use the Cartesian coordinate system O - xlx2x 3 and the notation as shown in Fig. 1. If the plate vibrates harmonically with the angular frequency o9, the governing equation of the problem can be expressed by </p><p>V4 w =1 0)2 phw (1) D </p><p>Where V 4 denotes the biharmonic differential operator, p the mass density of the plate material, h the plate thickness, w the lateral deflection of the plate, D the flexural rigidity of the plate. Suppose that the pseudo pressure q defined as </p><p>q=0)2phw (2) </p><p>is applied to the plate, equation (1) can be formulated into the integral equation by using the fundamental solution of the static bending problem: </p><p>, 1 z w*(x ,x )=~r lnr (3) </p><p>Now, we consider the weighted residual statement of equation (1) expressed as </p><p>f (V4w- q/D)w* df~ = 0 (4) </p><p>in which the approximate fundamental solution given by equation (3) is used as the weight function w*. Applying repeatedly the Gaussian divergence theorem, we can finally obtain the boundary integral equations as follows: </p><p>{" cw+ Jr (V*w-M*T. + T*M,-w*V.)dF </p><p>+ Z [[M*,w- w*M.3 </p><p>= 0) 2 jf~ phw*w c df~ (5a) </p><p>T*M.- if* V.) dF </p><p>+ </p><p>=0)~ I- ph~*w~df~ (5b) d, </p><p>which c is a constant depending only on the geometrical property of the boundary and equal to 1/2 if the boundary is smooth. The asterisked functions in equations (5) are the set of fundamental solutions corresponding to </p><p>equation (3), and those denoted by the tilde their normal derivatives. The deflection ~ multiplied by ~'* and also M*t denotes the relative displacement with respect to the observation point on the boundary, so that the corresponding boundary integrals can be evaluated in the sense of Cauchy's principal value integral as shown in Ref. 1 5. The variables T., M,, V. and M.t stand for the slope, bending moment, shearing force and twisting moment on the boundary, respectively. The quantity denoted by the symbol [1" -l] means the discontinuity jump at a corner point of the plate, and I]-M.t]] implies the concentrated force acting there. The notation ~ in front of the discontinuities indicates the summation over the whole corner points except for the one coinciding with the observation point. As has been shown in Refs 8 and 1 2, the discontinuity []-M,,]] can be approximated by the slope. 7". before and after the corner point under consideration by using the following relation: </p><p>~T. M., = - D(1 - v) ~- (6) </p><p>where v denotes Poisson's ratio and O( )lOt the derivative along tangential direction. Through such a procedure we can eliminate the corner variable JIM.t]] in the numerical implementation of the formulation. Other details with respect to the variables and the implementation of equations (5) can be found in the literature 12-14. </p><p>Since the approximate fundamental solution is used for the formulation, the r.h.s, of the boundary integral equations (5) includes the unknown variable w~ in the inner domain of the plate. Therefore, an additional equation is needed to solve the problem without any iterative solution procedure. For this purpose, we use the integral equation expressed by </p><p>w~+ fr (V*w-M*T,+ T*M,-w*V.)dF </p><p>+Z JIM*w- </p><p>= 09 2 [ phw*wc df~ Jn </p><p>(7) </p><p>The boundary integral equation (5) and the additional integral equation (7) are discretized by means of the boundary-domain element method 11 as in the previous paper for elastic buckling analysis of plate structures x4. The boundary-domain element method proposed in Ref. 11 can use any higher-order finite elements as the domain elements to discretize the domain integrals included in the integral equations (5) and (7). However, it should be noted that the basic idea of the boundary-domain element method for the free vibration problem of single plates was originally proposed by Bezine 8 and the constant domain element was successfully applied to numerical computation of some sample problems. If we employ the same solution scheme as Bezine's and take into account the homogeneous boundary conditions, we can finally obtain the discretized set of equations in the following matrix form: </p><p>[ GE] { XB} = (D2[JE] {Wc} (8) </p><p>{we) + [Go]{x.] = 0)2[JJ{w } (9) </p><p>Engineerin9 Analysis, 1988, Vol. 5, No. 4 183 </p></li><li><p>Free vibration analysis of elastic plate structures: M. Tanaka et al. </p><p>where {X~} denotes the column vector of nodal unknowns on the boundary and {w~} the column vector of nodal unknowns regarding the deflection w~ of the inner domain of the plate. Elimination of {XB} from the above two sets of equations can lead to the algebraic set of eigenvalue equations expressed as </p><p>jointed edge of plate components, we can finally obtain the global set of equations, which are solved under the appropriate boundary conditions prescribed on the periphery of the structure. </p><p>We first express the discretized sets of equations with respect to the integral equations (5) and (11) as follows: </p><p>([JE] -- [G~][G~] - x[SE]){We} = 5 {We} (10) </p><p>from which the eigen frequency and the eigen mode can be computed by using the subroutine library available at most computing centres. </p><p>Next, we shall show briefly a similar solution procedure for the in-plane free vibration problem of elastic plates. By assuming that the plate is in the plane stress state, and applying the approximate fundamental solution of the static problem, that is, Kelvin's solution to the integral equation formulation, we can obtain the boundary integral equation as follows X: </p><p>ck~u,,, + p*,,,u,, dF- u'~,,,p,, dF = po 2 u~mu m df~ dr dr </p><p>(11) </p><p>Also in this case, the r.h.s, of the boundary integral equation includes the unknown variable u,, in the inner domain of the plate. Hence, we use an additional integral equation expressed </p><p>Uk+ I </p><p>(12) </p><p>Discretization of equations (11) and (12) through the boundary-domain element method and taking account of the homogeneous boundary conditions can lead to </p><p>[HE]{XE] = coZ[RE]{Uc} (13) </p><p>{u~} + [H~] {XE} = ~02[RE] {Uc} (14) </p><p>Eliminating the nodal unknows {XE}, we can obtain </p><p>~ 1 ([R~] - [Hc][HE]- [Re]){ud =~- {u~} (15) </p><p>[S~] {~,, } - JR.] {} = o) 2 [AB] { ~;i.} </p><p>From equations (7) and (12) we can obtain </p><p>(16) </p><p>u ~2FA ]~u~ {~}+[S,]{r.}-[R,]{~}= t. ,j~w~, (17) </p><p>In equations (16) and (17) we have introduced the following notation: </p><p>{u}~=kut u2 w l </p><p>{"c}T=["~c .~A {f}r=tpth p2 h V.J </p><p>(18) </p><p>where the symbol [ J denotes a row vector. In order to transform the variables in the local </p><p>coordinate system to those in the global coordinate system, we use the following relationships: </p><p>{M.} =[ ]{M.# (19) </p><p>where the superimposed bar denotes the variables in the global coordinate system, and the transformation matrix [A] can be given by </p><p>[11 ml nl [A]= 12 m2 n2 (20) </p><p>13 m3 n3 0 0 0 </p><p>In equation (20) 1~, m~ and n~ (i = 1, 2, 3) are the directional cosines of the axes x'~, x~ and x; with respect to the global coordinate axis x~, respectively. </p><p>If equations (16) and (17) are expressed in the global coordinate system by using equation (19), it follows that </p><p>{S"] {t} -- [R"] {'M.} -- 2 [A.]{. J (21) </p><p>{~;} + [St]{t} - [/~] {~io} -co- 2[At]{w~},o (22) </p><p>where </p><p>This algebraic set of eigenvalue equations can be easily solved as in the out-of-plane free vibration problem. </p><p>[g.] =[S . ] [A] </p><p>[S,] =[Sl][A ] </p><p>= [R.] [A] </p><p>[R,] =[g , ] [A] (23) </p><p>3. FREE VIBRATION PROBLEM OF PLATE STRUCTURES </p><p>The out-of-plane and in-plane vibrations are combined in the free vibration problem of assembled plate structures. To analyse this problem we introduce a local coordinate system to each plate component of the structure. Then, the discretized sets of equations derived in the previous chapter can hold for each plate component. These equations expressed in the local coordinate system are assembled into the global coordinate system which is used in common for the whole structure. By taking account of the equilibrium and compatibility conditions on the </p><p>Equations (21) and (22) of every plate component are assembled into the global set of the whole structure. Then, we consider the compatibility and equilibrium conditions on the jointed edge of plate components. Suppose that L plate components are jointed together, these conditions can be expressed as follows: </p><p>(1) Compatibility of displacements {u} and slope T,; </p><p>(24) </p><p>fllT~=fl2T~ . . . . . flLT~ (25) </p><p>184 Engineering Analysis, 1988, Vol. 5, No. 4 </p></li><li><p>Free vibration analysis of elastic plate structures: M. Tanaka et al. </p><p>(2) Equilibrium of forces {f} and moment M,; </p><p>{fl} + {f2} +. . . + {fL} = 0 (26) </p><p>fllM~ +f12M2 +. . . +ilL M L = 0 (27) </p><p>It is assumed in equations (26) and (27) that no external forces and moments are applied at the jointed edge. In addition, the coefficient fli is defined such that #i= 1 when the plate component has the same tangential direction as that of the reference plate component and otherwise #i= -1 . </p><p>If the homogeneous boundary conditions expressed as </p><p>w = 0, T, = 0 (clamped) </p><p>w = 0, M, = 0 (simply supported) </p><p>M, = 0, V, = 0 (free) (28) </p><p>are taken into account, the assembled set of equations for the whole structure can be obtained as follows: </p><p>[H.]{X} =o~Z[G.]{u} (29) </p><p>Fig. 2. Element division of square plate </p><p>{u} (30) </p><p>where {u}r=[uc w d and {X} denotes the nodal unknowns on the boundary and the jointed edges of the structure. Elimination of {X} from equations (29) and (30) can finally lead to the following set of algebraic eigenvalue equations: </p><p>([Gc] - [Hc] [HB] -~ [G.]){u} = ~ {u} (31) </p><p>from which we can easily compute the eigen frequency o~ and the corresponding eigen mode. </p><p>4. NUMERICAL COMPUTATIONS AND DISCUSSION </p><p>A new computer program has been developed on the basis of theoretical considerations mentioned before. We will show here the numerical results obtained for some sample problems. It is noted that only constant boundary elements and also constant domain elements have been used for the present computations. </p><p>The first example is the out-of-plane bending vibration of a square plate. The boundary conditions are assumed such that the whole boundary is simply supported or clamped. The element division is illustrated in Fig. 2 where the whole boundary is equally divided and also the inner domain into triangular domain elements of equal size. To check the numerical accuracy computation is performed under the following conditions: </p><p>(1) The number of boundary elements is changed as 12, 16, 20, 32, and 48, while the number of domain elements is kept constant as 72 triangular elements. </p><p>(2) The number of triangular domain elements is changed as 8, 32 and 72, while the number of boundary elements is kept constant as 48 on the whole boundary. </p><p>The results obtained are summarized in Figs 3 and 4 for the clamped boundary conditions. The error in these </p><p>i_ o i,. t_ lid </p><p>Fig. 3. er ror </p><p>40.0 </p><p>50.0 </p><p>20.0 </p><p>I0.0 </p><p>I I I </p><p>\\\\\\\N\\ </p><p>lNNNNNX\/ </p><p> Mode I _ </p><p> Mode 2 </p><p>I - I I - </p><p>0 20 40 60 80 </p><p>Number of Internot Cetts Influence of boundary element mesh on numerical </p><p>figures are calculated such that </p><p>Numerical solution- Analytical solution Error(~o) = </p><p>Analytical solution </p><p>xl00 </p><p>In Tables 1 and 2 comparison is made between the numerical results obtained and the analytical ones 1'2. It can be seen that better accuracy can be achieved for lower eigen modes. </p><p>Next, we compute the in-plane free vibration problem of a rectangular plate as shown in Fig. 5. The dement division used is illustrated in Fig. 6. In Table 3 the numerical results...</p></li></ul>