Logo Logo

Open Access Journal

  • About Journal
  • Aims and scope
  • Editorial Board
  • For Authors
  • Special Issues
  • History
  • Contact
  • Statistics
  • Deklaracja dostępności


Issues:
Search In print
JAMCM
Vol. 24, 2025
Issue 1 Issue 2
Vol. 23, 2024
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 22, 2023
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 21, 2022
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 20, 2021
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 19, 2020
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 18, 2019
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 17, 2018
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 16, 2017
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 15, 2016
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 14, 2015
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 13, 2014
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 12, 2013
Issue 1 Issue 2 Issue 3 Issue 4
SRIMCS
Vol. 11, 2012
Issue 1 Issue 2 Issue 3 Issue 4
Vol. 10, 2011
Issue 1 Issue 2
Vol. 9, 2010
Issue 1 Issue 2
Vol. 8, 2009
Issue 1 Issue 2
Vol. 7, 2008
Issue 1 Issue 2
Vol. 6, 2007
Issue 1
Vol. 5, 2006
Issue 1
Vol. 4, 2005
Issue 1
Vol. 3, 2004
Issue 1
Vol. 2, 2003
Issue 1
Vol. 1, 2002
Issue 1
Article - HTML version

Vibration analysis of composite circular and annular membranes



Urszula Siedlecka

,

Izabela Zamorska

,

Stanisław Kukla


Journal of Applied Mathematics and Computational Mechanics
Year 2016, Volume 15, Issue 1, pages 149-159
DOI: 10.17512/jamcm.2016.1.15

PDF
Download
Full Text


VIBRATION ANALYSIS OF COMPOSITE CIRCULAR AND ANNULAR MEMBRANES

Urszula Siedlecka, Izabela Zamorska, Stanisław Kukla

Institute of Mathematics, Czestochowa University of Technology
Częstochowa, Poland
urszula.siedlecka@im.pcz.pl, izabela.zamorska@im.pcz.pl, stanislaw.kukla@im.pcz.pl

Abstract. In this paper a solution to the free vibration problem of composite circular and annular membranes is presented. The vibrations of membranes whose material densities and/or thicknesses varied step-wise with the radial co-ordinate are considered. This approach is applied to approximate the solution to the vibration problem of a membrane with continu- ously varying density and/or thickness with the radial co-ordinate. The obtained analytical solutions are used in numerical investigations into the effect of parameters characterizing the composite membranes on their eigenfrequencies.

Keywords: composite membrane, free vibration, Green’s function

1. Introduction

The vibration analysis of membranes is of particular interest in the design of various acoustical devices. The solution to the vibration problem of a non-homogenous membrane in a closed form can be derived only for the cases of some functions describing the change in the material density and thickness of this membrane. Free vibration problems of circular and annular membranes when the density varies with the radial co-ordinate are the subject of papers [1-5]. The solution to the vibration problem of a membrane comprising two concentric annular membranes has been derived in an exact form by Laura et al. in paper [1]. Gottlieb [2] gives the explicit values of the radial spectrum of an annular membrane with a stepped density which contains inverse fourth power logarithmic terms in the density function. The exact solutions to both the axisymmetric and antisymmetric modes of non-homogenous circular and annular membranes with polynomial variation of the density are given by Jabareen and Eisenberger in paper [3]. The eigenfrequencies for variable density membranes are obtained by a power series expansion and for multiple-connected regions by the dynamic stiffness method.

Approximate methods have been applied to solve the vibration problems of non-homogeneous membranes in numerous papers (for example in references [4-9] various approximate methods are used). An application of the boundary point collocation method to determine the eigenfrequencies of membranes with varying mass density is presented by Cap in paper [4]. Gutierrez et al. [5] present numerical results for two lower free vibration frequencies of circular and annular membranes whose densities varied linearly, quadratically or cubically with the radial coordinate. The frequencies were calculated by using the differential quadrature method, the finite element technique, an optimized and/or improved Rayleigh quotient method and a lower bound based on the Stodola-Vianello method. The multi-symplectic methods for free vibration of the membrane are proposed by Wei-Peng et al. in paper [6]. Numerical results presented in the paper verify the efficiency of the methods. The free vibration frequencies of an annular membrane for axisymmetric modes by the discrete singular convolution method, based on the regularized Shannon’s delta kernel, were determined by Civalek and Gürses in reference [7]. The numerical technique for problems of free vibrations of non-homogenous membranes is presented by Reutskiy in paper [8]. The method is based on the mathematical modeling of the response of a system to external excitation over a range of frequencies. Buchanan [9] has studied vibrational properties and has demonstrated the accuracy of finite element formulation for circular and annular membranes with density variation assumed as a linear function of the radius. The free vibration of an annular membrane consisting of three concentric annular membranes was considered in the paper [10]. The solution of the problem has been derived by using the properties of Green’s functions.

In this paper, a solution to the free vibration problem of composite circular and annular membranes whose densities and/or thicknesses varied step-wise with the radial co-ordinate is derived. This approach has been applied for numerical computation of eigenfrequencies of membranes with a continuous change in density and thickness in the radial direction.

2. Formulation and solution to the free vibration problem of a composite membrane

Consider the vibration problem of an annular composite membrane consisting of m annular homogenous membranes as shown in Figure 1. The thickness and the material density of the composite membrane change step-wise at circles with radii . The j-th annular membrane includes an elastic support distributed along the circle , where .

Free vibrations of the homogenous annular membranes with an elastic support are governed by the following differential equations:

, (1)

where is the displacement of the j-th annular membrane, is the tension per unit length, is the mass per unit area, is the thickness of the j-th annular membrane, is the Laplacian operator, and are polar coordinates and is time. The radial variable for the j-th annular membrane is in the interval: , , where , .

Fig. 1. The sketch of the composite annular membrane under study

Functions satisfy the continuity conditions

, (2)
, (3)

and the boundary conditions

, (4)

Considering the free vibration of the membrane we assume functions in the form

, , (5)

where is the natural frequency of the composite membrane. Taking equation (6) into account in differential equation (1), in continuity conditions (3)-(4) and in boundary conditions (5), we obtain a differential equation, continuity conditions and boundary conditions for functions :

for , (6)
, (7)
, (8)

(9)

where .

The general solution of the differential equation (7) is given by

for , (10)

where and are arbitrary constants. Substituting the functions (10) into conditions (7)-(9), we obtain a set of equations which can be written in a matrix form

(11)

where , and . For (an annular membrane) the non-zero elements of the matrix A are (index n is omitted)

, , , ,

, , ,

, ,

for and

, .

For a nontrivial solution of the equation (11) the determinant of the matrix A is set equal to zero yielding the frequency equation of the composite membrane

(12)

The equation (12) is then solved numerically with respect to by using an approxi- mate method. Note that for a fixed we obtain a sequence , , of the roots of the equation (12).

The frequency equation for a circular membrane () is obtained similarly. In this case, the first row in the matrix A for the annular membrane should be changed by assuming: , . The remainder elements of the matrix A for the circular membrane are the same as for annular membrane.

For the computed values , , the corresponding eigenfunctions (the mode shapes of vibration) are given by equation (10) where the coefficients , are determined by solving equation (11), in which should be assumed. The eigenfunctions satisfied the following orthogonality conditions

(13)

where and

(14)

3. Forced vibration of a composite membrane

The differential equation to vibration of an annular membrane forced by an outer force is

, (15)

The functions satisfy the continuity conditions (2)-(3) and the boundary conditions (4). Moreover, the initial conditions are

(16)

(17)

We seek a solution of the problem in the form of a series

(18)

where were derived in Section 2 (equation (10)) and functions will be now determined.

Substituting the function in the form (18) into equation (15) and using the orthogonality condition (13), we obtain an equation with unknown functions in the form

(19)

where and , for . This equation is complemented by initial conditions which follow from (16) and (17). Using (18) and (13) one obtains the initial conditions in the form

(20)
(21)

The solution of the equation (19) with conditions (20), (21) is as follows:

(22)

Finally, the solution of the forced vibration problem of the composite membrane can be written in the form

(23)

4. Numerical examples

The numerical computations presented here concern the frequency analysis of composite annular and circular membranes for various values of parameters which characterize their non-uniformity. The calculations of the non-dimensional free vibration frequencies were performed using frequency equation (12). The roots of this equation were determined by the application of the false position method [11].

In the first example, the eigenfrequencies of circular and annular membrane were computed with the density function given by the formula: . The results of calculations: (;) for various values of and various numbers of annular membranes are shown in Table 1 - for the annular membrane and in Table 2 - for the circular membrane. The eigenfrequencies for these membranes were determined earlier by using the power series method in paper [3] by Jabareen and Eisenberger. The results for , obtained by using the presented method, are compared with the results given in reference [3]. The calculated free vibration frequencies are consistent for both annular (Table 1) and circular (Table 2) membranes for and those presented in reference [3]. The differences do not exceed 5∙10–4 for and 3∙10–2 for . Moreover, it can be shown that the differences decrease as the number of annular membranes is increased.

The frequencies of the circular composite membrane of radius consisting of three parts (a circular inner membrane and two annular membranes) as functions of the ratio , are presented in Figure 2 for The material density of the inner annular membrane is much greater than the density of the other parts of the membrane, i.e. , and . The calculations were performed for: and Figure 2 shows that both the material density (or membrane thickness) and the location of the inner annular membrane cause significant changes in the eigenfrequencies of the composite membrane.

Table 1

Values of for an annular membrane,

a = 1.0

a = 2.0

m

Ref [3]

m

Ref [3]

10

15

20

10

15

20

W01

3.0156

3.0156

3.0156

3.0156

2.5691

2.5691

2.5691

2.5691

W02

6.1639

6.1640

6.1640

6.1640

5.2653

5.2654

5.2654

5.2654

W03

9.2537

9.2937

9.2937

9.2937

7.9442

7.9442

7.9443

7.9443

W04

12.4161

12.4163

12.4163

12.4163

10.6161

10.6163

10.6164

10.6164

W05

15.5350

15.5354

15.5354

15.5354

13.2845

13.2851

13.2852

13.2852

W11

3.3349

3.3349

3.3349

3.3349

2.8360

2.8360

2.8360

2.8360

W12

6.3801

6.3801

6.3802

6.3802

5.4503

5.4504

5.4504

5.4504

W13

9.4532

9.4532

9.4532

9.4532

8.0825

8.0826

8.0826

8.0826

W14

12.5414

12.5415

12.5416

12.5416

10.7255

10.7258

10.7258

10.7259

W15

15.6377

15.6381

15.6381

15.6381

13.3746

13.3753

13.3753

13.3754

W21

4.0734

4.0734

4.0734

4.0734

3.4492

3.4492

3.4492

3.4492

W22

6.9718

6.9718

6.9718

6.9718

5.9526

5.9526

5.9526

5.9526

W23

9.9155

9.9156

9.9156

9.9156

8.4816

8.4817

8.4817

8.4817

W24

12.9120

12.9122

12.9122

12.9123

11.0488

11.0491

11.0492

11.0492

W25

15.9440

15.9444

15.9445

15.9445

13.6434

13.6441

13.6442

13.6443

W31

4.9359

4.9359

4.9359

4.9359

4.1603

4.1603

4.1603

4.1603

W32

7.7955

7.7956

7.7956

7.7956

6.6424

6.6424

6.6425

6.6425

W33

10.6269

10.6271

10.6271

10.6271

9.0889

9.0891

9.0892

9.0892

W34

13.5092

13.5095

13.5095

13.5096

11.5663

11.5668

11.5669

11.5669

W35

16.4474

16.4480

16.4481

16.4481

14.0838

14.0848

14.0850

14.0850

Table 2

Values of for a circular membrane,

a = 1.0

a = 2.0

m

Ref [3]

m

Ref [3]

10

15

20

10

15

20

W01

2.0108

2.0108

2.0108

2.0108

1.7598

1.7598

1.7598

1.7598

W02

4.5548

4.5549

4.5549

4.5549

3.9802

3.9802

3.9802

3.9802

W03

7.1194

7.1195

7.1196

7.1196

6.2147

6.2149

6.2150

6.2150

W04

9.6897

9.6900

9.6901

9.6901

8.4537

8.4542

8.4543

8.4543

W05

12.2622

12.2629

12.2631

12.2631

10.6943

10.6956

10.6957

10.6958

W11

3.0678

3.0678

3.0678

3.0678

2.6273

2.6273

2.6273

2.6273

W12

5.6990

5.6991

5.6991

5.6991

4.9267

4.9268

4.9268

4.9268

W13

8.3004

8.3006

8.3006

8.3006

7.1990

7.1993

7.1993

7.1994

W14

10.8914

10.8919

10.8920

10.8920

9.4608

9.4616

9.4617

9.4618

W15

13.4773

13.4785

13.4786

13.4787

11.7170

11.7188

11.7191

11.7192

W21

4.0223

4.0223

4.0224

4.0224

3.4110

3.4110

3.4110

3.4110

W22

6.7453

6.7454

6.7454

6.7454

5.7872

5.7874

5.7874

5.7874

W23

9.3940

9.3943

9.3944

9.3944

8.1011

8.1015

8.1016

8.1016

W24

12.0150

12.0157

12.0158

12.0159

10.3905

10.3916

10.3918

10.3919

W25

14.6217

14.6234

14.6236

14.6236

12.6665

12.6691

12.6694

12.6696

W31

4.9283

4.9283

4.9284

4.9284

4.1549

4.1549

4.1550

4.1550

W32

7.7381

7.7383

7.7383

7.7383

6.6029

6.6031

6.6032

6.6032

W33

10.4357

10.4362

10.4362

10.4363

8.9581

8.9587

8.9588

8.9588

W34

13.0894

13.0904

13.0906

13.0906

11.2758

11.2774

11.2776

11.2776

W35

15.7197

15.7220

15.7223

15.7224

13.5726

13.5763

13.5767

13.5768

Fig. 2. Four free vibration frequencies of the circular composite membrane consisting of three segments for n = 0, as functions of the ratio for various and

5. Conclusions

In this paper, a solution to the free vibration problem of annular and circular composite membranes has been derived. Numerical analysis has shown the effect of parameters characterizing the composite annular or circular membranes on their eigenfrequencies. The frequency analysis was performed for a circular membrane with an inner highlighted annular membrane whose material density (or thickness) is much greater than the material density (or thickness) of the remaining part of the membrane. The numerical results show that the location of the highlighted annular membrane has a significant effect on eigenfrequencies of the composite membrane. It follows that the application of this annular membrane can be used to introduce a specific change in the dynamical characteristic of the composite membrane. Numerical examples show that the presented method can be used to determine approximate natural frequencies of composite membranes with continuously varying densities or thicknesses.

References

[1] Laura P.A.A., Rossit C.A., La Malfa S., Transverse vibrations of composite, circular annular membranes: exact solution, Journal of Sound and Vibration 1998, 216(1), 190-193.

[2] Gottlieb H.P.W., Exact solutions for vibrations of some annular membranes with inhomogeneous radial densities, Journal of Sound and Vibration 2000, 223(1), 165-170.

[3] Jabareen M., Eisenberger M., Free vibrations of non-homogeneous circular and annular membranes, Journal of Sound and Vibration 2001, 240(3), 409-429.

[4] Cap F.F., Eigenfrequencies of membranes of arbitrary boundary and with varying surface mass density, Applied Mathematics and Computation 2001, 124, 319-329.

[5] Gutierrez A.A., Laura P.A.A., Bambill D.V., Jederlinic V.A., Hodges D.H., Axisymmetric vibrations of solid circular and annular membranes with continuously varying density, Journal of Sound and Vibration 1998, 212(4), 611-622.

[6] Wei-Peng H., Zi-Chen D., Wen-Cheng L., Multi-symplectic methods for membrane free vibration equation, Applied Mathematics and Mechanics 2007, 28(9), 1181-1189.

[7] Civalek Ö., Gürses M., Discrete singular convolution for free vibration analysis annular membranes, Mathematical and Computational Applications 2009, 14(2), 131-138.

[8] Reutskiy S.Y., The methods of external and internal excitation for problems of free vibrations of non-homogeneous membranes, Engineering Analysis with Boundary Elements 2007, 31, 906-918.

[9] Buchanan G.R., Vibration of circular membranes with linearly varying density along a diameter, Journal of Sound and Vibration 2005, 280, 407-414.

[10] Zamorska I., Kukla S., Siedlecka U., Frequency analysis of composite annular membranes, Scientific Research of the Institute of Mathematics and Computer Science 2012, 11(1), 129-135.

[11] Dahlquist G., Björck A., Numerical Methods in Scientific Computing, Society for Industrial and Applied Mathematics, Philadelphia 2008.


Journal of Applied Mathematics and Computational Mechanics
p-ISSN: 2299-9965, e-ISSN: 2353-0588
Editorial address: Department of Mathematics, Czestochowa University of Technology, Armii Krajowej 21, 42-200 Częstochowa, Poland
E-mail: jamcm@pcz.pl