MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi Publishing Corporation 651356 10.1155/2014/651356 651356 Research Article Nonlinear Free Vibration Analysis of Axisymmetric Polar Orthotropic Circular Membranes under the Fixed Boundary Condition Zheng Zhoulian 1, 2 http://orcid.org/0000-0001-6414-4765 Guo Jianjun 3 Song Weiju 2 He Xiaoting 2 Lu Faming 2 Xie Chuanxi 2 Sun Junyi 2 Nguyen Thoi Trung 1 Department of Civil Engineering Chongqing University Chongqing 400045 China cqu.edu.cn 2 Chongqing Vocational College of Architectural Engineering Chongqing 400039 China 3 Chongqing Water Resources and Electric Engineering College Chongqing 402160 China 2014 25 2 2014 2014 26 08 2013 20 12 2013 15 01 2014 25 2 2014 2014 Copyright © 2014 Zhoulian Zheng et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper presents the nonlinear free vibration analysis of axisymmetric polar orthotropic circular membrane, based on the large deflection theory of membrane and the principle of virtual displacement. We have derived the governing equations of nonlinear free vibration of circular membrane and solved them by the Galerkin method and the Bessel function to obtain the generally exact formula of nonlinear vibration frequency of circular membrane with outer edges fixed. The formula could be degraded into the solution from small deflection vibration; thus, its correctness has been verified. Finally, the paper gives the computational examples and comparative analysis with the other solution. The frequency is enlarged with the increase of the initial displacement, and the larger the initial displacement is, the larger the effect on the frequency is, and vice versa. When the initial displacement approaches zero, the result is consistent with that obtained on the basis of the small deflection theory. Results obtained from this paper provide the accurate theory for the measurement of the pretension of polar orthotropic composite materials by frequency method and some theoretical basis for the research of the dynamic response of membrane structure.

1. Introduction

The membrane structure is a kind of large span flexible space structure, which has a wide range of applications as discussed elsewhere . However, due to its small weight and stiffness and low natural frequency as well as small damping, it is very sensitive to the action of wind, rain, and other external loads and easy to produce vibration as discussed elsewhere . The deformation of membrane structure, caused by the vibration, leads to the destruction of membrane structure in practical engineering. Since this kind of case is common, studying the vibration of membrane structure and its dynamic response under the external excitation becomes the key scientific issues in the field of membrane structure. In actual engineering, the dynamic characteristics of the membrane structure mainly show up as the vibration problem in large deflection, when the membrane is under impact loading.

From the last century, many scholars at home and abroad have done a lot of research work about it. The problem of large deflection of elastic circular membrane was studied by Chien et al. . A nonlinear relaxation method was employed to solve the nonlinear partial differential equations governing the large deflection response of various axisymmetric circular membranes, by Kao and Perrone . The problem of large deflections of rotationally orthotropic circular membrane was investigated by Cheng and Yang , in which the relations between the load and membrane forces as well as the deflection and the central deflection were obtained. The large deformation of circular membrane under the concentrated force was analyzed by Chen and Zheng . Kang  dealt with a composite rectangular membrane with an oblique interface. Young et al.  made a numerical calculation and experimental study on the dynamic characteristics of thin-film membranes. The nonlinear vibration analysis of a prestretched hyperelastic annular membrane under finite deformations was completed by Soares and Gonalves . Qian  researched the vibration problem of the rectangular membrane with two different directions of force, circular membrane, and elliptic membrane and obtained the approximate solution of the vibration frequency of arbitrary shape membrane. Wu  obtained the analytical solution of natural vibration modes of the intermediate support membrane by the Laplace transform. Lin and Chen  used analytical method to obtain the theoretical solution of free vibration of annular isotropic membrane. Ou yang  studied a Robin problem of the vibration of a perturbed boundary membrane, using parametric deformation method and the condition for the solvability of the problem. The nondegenerate vibrational frequency and vibration mode of this problem are obtained. Lin and Chen  used the finite element method to analyze the free vibration of flat membrane. Qiao et al.  used the Newmark method of dynamic analysis and the Newton-Raphson iterative method of nonlinear analysis to program composition for the nonlinear dynamic response analysis of membrane structure. However, the theoretical research on the nonlinear vibration of orthotropic circular membrane is relatively rare.

In this paper, according to the large deflection theory of membrane, we study the vibration of axisymmetric polar orthotropic circular membrane and use the principle of virtual displacement to create the governing vibration equations of circular membrane with outer edges fixed. The Bessel function and the Galerkin method are used to get a generally accurate expression of the nonlinear free vibration frequency of circular membrane, and the results obtained in this paper are compared with other existing results. The work presented in this paper provides a theoretical basis for the study of the dynamic response of polar orthotropic membrane materials.

2. Basic Equations

The studied circular membrane is axisymmetric polar orthotropic and it is assumed that all outer edges of the membrane are fixed. r denotes the radius of the membrane, h denotes the thickness of the membrane, ρ denotes aerial density of the membrane material, and Nr and Nθ denote initial tension along radial and circumferential direction, respectively, as shown in Figure 1.

Circular membrane with all edges fixed and force diagram of the element.

According to the large deflection theory of membrane, the equations are established, which are based on the hypothesis as follows.

The membrane is soft and springy, which cannot resist moment, and the tension is in the tangent plane at any time.

There is only transverse vibration without considering the damping in vibration process.

The influence of membranes thickness variation in vibration process can be neglected.

The paper only considers geometric nonlinearity for membranes with the hypothesis that the membrane is homogeneous, continuous, and linearly elastic in vibration process.

The corresponding boundary conditions can be expressed as follows: (1)w|r=a=0,ur|r=a=0, where w denotes the deflection of membranes and ur denotes the radial displacement.

The geometrical equation of axisymmetric deformation of the circular membrane in polar coordinates is given by (2)εr=dudr+12(dwdr)2,εθ=ur.

The physical equations are (3)Nr=hσr=Erh1-μrμθ(εr+μθεθ),Nθ=hσθ=Eθh1-μrμθ(εθ+μrεr), where Er and Eθ represent Young’s modulus along radial and circumferential direction, respectively; μrμθ represents Poisson’s ratio, respectively; ε represents the strain; σ represents the stress.

The equation of compatibility is (4)εr=ddr(rεθ)+12(dw(r,t)dr)2 and εr and εθ can be transformed into the below form, with (3); that is, (5)εr=NrErh-NθμθEθh,εθ=NθEθh-NrμrErh.

Substituting (5) into (4) and letting λ2=Eθ/Er=μθ/μr yield (6)rd2dr2(rNr)+d(rNr)dr-λ2Nr=-12Eθh(dwdr)2.

Introducing the stress function σr=(1/r)(dϕ/dr), σθ=d2ϕ/dr2 and letting Nr=hσr=(h/r)(dϕ/dr), Nθ=hσθ=h(d2ϕ/dr2), (6) can be simplified as follows: (7)rd3ϕdr3+d2ϕdr2-λ21rdϕdr=-12Eθ(dw(r,t)dr)2.

According to the principle of virtual displacement, considering the vibration inertia power of membranes as distributed load, the equation of the large deflection vibration of membranes can be given as follows: (8)0a[ψ+h(σr+σ0r)wr]δwrrdr=0, where σ0r denotes pretension of circular membrane; load function is ψ=(1/r)0r(-ρhw¨)rdr.

Taking the small deflection theory of membrane into account, every piece element of membranes has no elongation; therefore, according to the elastic Hooke’s law, the membrane is given a constant tension. The equation of the small deflection vibration of membranes should be (9)0a(ψ+hσ0rwr)δwrrdr=0.

3. Governing Equation

The equation of transverse free vibration of membranes is (10)2w=ρN2wt2, where w denotes the deflection of membranes. On the polar coordinates, 2=2/r2+(1/r)(/r)+(1/r2)(2/θ2); for the axisymmetric circular membrane, 2=d2/dr2+(1/r)(d/dr).

In order to solve the vibration differential equation, through the method of separation of variables, the functions which satisfy the boundary conditions (1) are taken as follows: (11)w(r,t)=W(r)T(t).

Substituting (11) into (6) yields (12)2W·T(t)=ρN·T¨(t)·W. Letting μ2=ρT¨(t)/NT(t), then (13)2W=ρNμ2W. Letting (ρ/N)μ2=γ2, then (14)(2-γ2)W=0,(15)  (2r2+1rr-γ2)W=0. Let the shape function of displacement be (16)W=R(r)cosnθ, where n=0,1,2,3,. When n=0, the shape function is symmetrical. In this paper, the deformation of circular membranes is axisymmetric, so n=0.

Substituting W=R(r)cosnθ into (15), an ordinary differential equation is obtained: (17)d2Rdr2+1rdRdr+(γ2-n2r2)R=0.

A dimensionless parameter (x=γr) is introduced; we have (18)x2d2Rdr2+xdRdr+(x2-n2)R=0.

A general solution for this differential equation is (19)R=C1Jn(x)+C2Nn(x), where Jn(x) and Nn(x) are the first and second kind Bessel functions of the Order-N real variable, respectively: (20)Jn(x)=k=0(-1)k1k!Γ(n+k+1)(x2)2k+n(|argz|<π),Nn(x)=Jn(x)cosnπ-J-n(x)sinnπ(|z|<,|argz|<π).

Equation (11) is only the general solution of the governing differential equation. Only after it satisfies concrete boundary conditions does it meet the practical engineering problems. Due to the clamped edges of the circular membrane, we have (21)w|r=a=0,|w|<,wr|r=a=0, where a denotes the radius of the circular membrane.

N n ( x ) and Jn(x) are infinite. Because |w|<, when x=γr=0, so C2=0.

And (19) can be simplified as (22)R=C1Jn(γir),W(r)=C1Jn(γir)cosnθ, And when n=0, the vibration mode is symmetrical. In this paper, the deformation of the circular membrane is axisymmetric, so n=0; namely, (23)W(r)=C1J0(γir).

Substituting the boundary conditions into (19) and (23) yields (24)C1J0(γia)=0.

Because C1 is not identically equal to zero, there is γi, satisfying (25)J0(γia)=0.

Namely, γia is the zero solution of the J0(x).

Using the MATLAB software to solve the six zero point of the first kind zero-order Bessel function, we have rootBessel = 2.4048, 5.5201, 8.6537, 11.7915, 14.9309, 18.0711.

Obtaining the Bessel function expansion by definition, we have (26)J0(x)=1-x222+x424·(2!)2-x626·(3!)2+.

The front four of J0(x) will be used in the following, namely: (27)J0(x)=1-x24+x464-x62304.

Substituting (27) into (11), the deflection expression is(28)w(r,t)=J0(γir)T(t).

Substituting (27) into (28) yields (29)w(r,t)=(b+cr2a2+dr4a4+er6a6)T(t),w(r,t)hhhhhhhhhhhb,d,>0,c,e<0.

Substituting (29) into (6) and solving the Euler equation yield (30)dϕdr=-Eθ2XT2(t)+C1rλ+C2r-λ, where (31)X=4c2(9-λ2)a4r3+16d2+24ce(49-λ2)a8r7+36e2(121-λ2)a12r11X=+16cd(25-λ2)a6r5+48de(81-λ2)a10r9. when r=0, the stress (σr) is finite, so (32)C2=0.

At the same time, the boundary condition is: ur|r=a=0, yields: (33)rεθ|r=a=0,(rEθd2ϕdr2-μrErdϕdr)|r=a=0.

Substituting (30) into (33) yields (34)C1=Eθ2(λ-μθ)aλ-1YT2(t), where (35)Y=a-2{4c2(3-μθ)(9-λ2)+16(d2+24ce)(7-μθ)(49-λ2)Yhhhh0hh+36e2(11-μθ)(121-λ2)+16cd(5-μθ)(25-λ2)hYhhhh0h+48de(9-μθ)(81-λ2)4c2(3-μθ)(9-λ2)4c2(3-μθ)(9-λ2)}.

Substituting (32) and (34) into (30) yields (36)dϕdr=[-Eθ2X+YEθ2(λ-μθ)aλ-1rλ]T2(t).

Substituting (29) into (8), according to the Galerkin method, yields (37)0aQdW(r)drrdr=0. Here, Q=ψ+h((1/r)(dϕ/dr)+σ0r)(w/r)(38)0aψdW(r)drrdr=0a1r0r(-ρhw¨)rdrdW(r)drrdr0aψdW(r)drrdr=-Sa2ρhT¨(t), where (39)0ah(1rdϕdr+σ0r)wrdW(r)drrdr=(-h2EθF1+YEθh2(λ-μθ)F2)T3(t)+H1hσ0rT(t), where (40)F1=a-2{(81-λ2)((13c2+d2+95e2-1(23c2+85d2+187e2+2cd+4de+125ce)0000000×(9-λ2)-10000000+(16d2+24ce)0000000×(25c2+87d2+2e2+43cd+3de+127ce)0000000×(49-λ2)-10000000+36e2(27c2+89d2+1811e2hhhhhhh000hhhhh+cd+125de+43ce)0000000×(121-λ2)-10000000+16cd(12c2+43d2+94e2hhhhhh00hhhhhhh+85cd+247de+2ce)0000000×(25-λ2)-10000000+48de(13c2+d2+95e2000000000000000+87cd+83de+32ce)0000000×(81-λ2)-1(((23)c2(81-λ2)((13c2+d2+95e2-1+(85)d2+(187)e2},F2=4c23+λ+16d2+24ce7+λ+36e211+λ+16cd5+λ+48de9+λ,H1=c2+2d2+3e2+83cd+245de+3ce.

Substituting all the above into (8) yields (41)d2T(t)dt2+(-H1σ0rSa2ρ)T(t)d2T(t)dt2+[EθF12Sa2ρ-EθY2Sa2ρ(λ-μθ)]T3(t)=0. Letting M=-H1σ0r/Sa2ρ, N=EθF1/2Sa2ρ-EθY/2Sa2ρ(λ-μθ) yields (42)d2T(t)dt2+M·T+N·T3(t)=0.

4. Solution of Free Vibration Frequency

Integrating (42) yields (43)(dT(t)dt)2+M·T2(t)+N2·T4(t)=K.

In (43), the value of K is determined by the initial conditions. Assuming that the initial displacement is T|t=0=T0, T0 is the amplitude of the membrane, so the initial velocity is (44)dT(t)dt|t=0=0.

Substituting T|t=0=T0 and (44) into (43) yields (45)K=MT02+N2T04.

Substituting K=MT02+(N/2)T04 into (43) yields (46)dT(t)dt=T02MT02+N2dT(t)dt·(1-T2T02)·(1+NT022M+NT02T2T02). Letting ξ=T02M/T02+N/2, k2=NT02/(2M+NT02).

Integrating (46) by the method of separation of variables, we can obtain the period of the vibration of the membrane: (47)Z=4ξ0T01(1-(T2/T02))(1+k2·(T2/T02))dT. Letting T/T0=sinθ, then (47) can be simplified as follows: (48)Z=4T0ξ0π/211+k2·sin2θdθ,(0k·sinθ1), where (1+k2·sin2θ)-1/2 could be spread as a power series with respect to k·sinθ: (49)(1+k2·sin2θ)-1/2=1-12k2sin2θ+1·32·4k4sin4θ(1+k2·sin2θ)-1/2=++(-1)n(2n-1)!!(2n)!!(ksinθ)2n.

Substituting (49) into (48) and then solving (48) through integrating item by item yield (50)Z=4T0ξ0π/2[(2n-1)!!(2n)!!1-12k2sin2θ+1·32·4k4sin4θZhhhhhhhhhhh+(-1)n(2n-1)!!(2n)!!(ksinθ)2n]dθZ=2πT0ξ[((2p-1)(2p))21-(12)2k2+(1·32·4)2k2hhZhhhhh++(-1)p((2p-1)(2p))2k2p]Z=2πT0ξn=0(-1)p((2p-1)!!(2p)!!)2k2p, where p=n=0,1,2,3,.

Therefore, the vibration frequency of the circular membrane is (51)ω=2πZω=M+(N/2)T02p=0(-1)p((2p-1)!!/(2p)!!)2(NT02/(2M+NT02))p, where M=-H1σ0r/Sa2ρ, N=EθF1/2Sa2ρ-EθY/2Sa2ρ(λ-μθ), and T0 is the amplitude of the circular membrane. Next, we will solve the free vibration frequency of circular membrane in small deflection.

Substituting (29) into (9) yields (52)d2T(t)dt2+(-H1σ0rSa2ρ)T(t)=0.

Equation (52) can be seen as (43) when N=0; substituting N=0 into (51) yields (53)ω=1aσ0rρ·-H1S.

Directly letting T00 in (51), the situation may degenerate into a small deflection. Equation (51) can be turned into the free vibration frequency of circular membrane in small deflection: (54)ω=1aσ0rρ·-H1S.

So, the frequency formula which is deduced directly out according to the vibration of circular membrane in small deflection is equal to the frequency formula which is based on the vibration of circular membrane in large deflection. This proves that the result of this paper is correct.

5. Computational Examples and Discussion

Take the membrane material commonly applied in engineering as an example. Young’s moduli are Er=1.4×106 kN/m2 and Eθ=0.9×106 kN/m2, respectively; Poisson’s ratios are μr=0.3 and μθ=0.193, respectively; the aerial density of the membrane materials is ρ=1.7 kg/m2; the circular membrane’s thickness is h=1.0 mm. The first order vibration frequency of the membrane is calculated by (51).

We can draw the conclusion from the result of Table 1: The initial displacement (the amplitude) has influenced the vibration frequency of circular membrane on the basis of the large deflection theory. The frequency is enlarged with the increase of the initial displacement, and the larger the initial displacement is, the larger the effect on the frequency is, and vice versa. When the initial displacement approaches zero, the result is consistent with that obtained in this paper, according to the small deflection theory.

The frequency under different initial displacements (r=0.5 m and σ0r=5.0×103 kN/m2).

ω T 0 = 0.10 T 0 = 0.08 T 0 = 0.04 T 0 = 0.01
Hz 117.075 97.090 60.741 43.215

ω T 0 = 0.005 T 0 = 0.001 T 0 = 0.0001 T 0 0

Hz 42.125 41.769 41.754 41.754

In Figure 2, we have drawn the relation curve of frequency and initial displacement. “Result Curve 1” and “Result Contrast 1” are based on large deflection; “Result Curve 2” and “Result Contrast 2” are based on small deflection; “Result Contrast 1” and “Result Contrast 2” are the result for the other method which is used in . It could be concluded that when the material properties and the initial pretension are given, there is nonlinear relationship between the vibration frequency and the initial deflection for the membrane in the large deflection; when the initial deflection (T0) is increasing, not only the vibration frequency is increasing, but also the slope of the curve will be slowly increasing, meaning dω/dT0 is gradually increasing. This also shows that the larger T0 is, the greater impact on ω it has. When the initial deflection approaches zero, the result is consistent with that obtained from the small deflection theory. Figure 2 reflects substantially the geometrical nonlinearity of the vibration of membranes; namely, formula (51) can reflect the geometric nonlinear characteristics of the vibration of membranes.

Frequency(ω)-initial displacement (T0) curve.

We can draw the conclusion from the comparative analysis of the relation curves of the frequency and the initial deflection drawn from this paper and . The results of the two kinds of formula are consistent; they both reflect the nonlinear characteristics of the membrane vibration; therefore, the expression of frequency obtained in this paper is correct. In [12, 17], the Bessel equation was deduced, so using the Bessel function as the vibration function of the membrane is reasonable. The more the number of the expansions of the Bessel function is, the more accurate the frequency is. In practice, the number of the expansions is determined by necessity.

Compared to [16, 18, 19], in this paper, more precise displacement potential function is applied, and the nonlinear vibration and orthotropy of membrane are taken into consideration too. In [18, 19], the thin circular plate is the object of study. In  the nonlinear vibration of thin circular plate has not been considered, and the initial pretension is imposed on the membrane in our paper, so it is different from . In , the vibration displacement potential function of circular membrane is based on the vibration of the thin circular plate; its precision is not enough; the vibration of the plate and the vibration of the thin membrane are different. In this paper, the Bessel function is chosen as the vibration displacement potential function of circular membrane; it is more precise.

This formula can also be used to solve other membrane structures with different geometry and boundary conditions. For example, assume a square membrane structure with 1.0 m in side length; then Young’s moduli are Ex=1.52×106 kN/m2 and Ey=1.29×106 kN/m2, respectively; Poisson’s ratios are μx=0.4 and μy=0.39, respectively; the aerial density of the membrane materials is ρ=1.05 kg/m2; the membrane’s thickness is h=0.82 mm; the initial displacement is T0. We used formula (51) to calculate the vibration frequency ω1, roughly. In Table 2, the experimental results ω2 are given too.

The frequency under various initial displacements and various pretensions (σox=σoy=σο).

 T ο (m) 0.0022583 0.0023742 0.001593 0.0016005 σ o (kN/m2) 1200 2400 3700 4800 ω 1 (Hz) 13.2160 18.5630 22.9094 26.0790 ω 2 (Hz) 15 19 21 25 T 0 (m) 0.0015625 0.0014709 0.0013644 0.001361 σ o (kN/m2) 6100 7300 8500 9800 ω 1 (Hz) 29.3847 32.1331 34.6638 37.2162 ω 2 (Hz) 28 31 33 36

In Figure 3, the calculation data is compared with the experimental data. The calculation results are consistent with the results obtained from the experiment.

Frequency(ω)-pretension (σ) curve.

6. Conclusions

We apply the large deflection theory of membrane and the principle of virtual displacement to derive the governing vibration equations of circular membrane with outer edges fixed, use the Bessel function to establish the accurate vibration mode function of the circular membrane, and then solve them by the Galerkin method. At last, we obtain the generally accurate formula of nonlinear vibration frequency of circular membrane with outer edges fixed. And the frequency of the vibration of circular membrane in small deflection which is degenerated by the vibration of circular membrane in large deflection is equal to the frequency which is deduced directly out, according to the vibration of circular membrane in small deflection. According to [18, 19], the formula of the vibration frequency for orthotropic membrane, which is based on the large deflection theory, is correct. The theoretical result is of great significance to study the dynamic response of membranes with forced vibration in large deflection and is also available to verify the correctness of the other numerical methods.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research presented in this paper was supported by the National Natural Science Foundation of China (Project no. 51178485).

Shen S. Z. A study on some theoretical aspects of long-span spatial structures Journal of Suzhou Institute of Urban Construction and Environmental Protection 2000 13 3 1 8 Chang H. S. Development and strategy of reliability for long-span spatial The World of Building Materials 2012 33 2 111 113 Chien W. Z. Lin H. S. Hu H. C. Yeh K. Y. The Problem of Large Deflection of Elastic Circular Membranes 1953 Beijing China Science Press Kao R. Perrone N. Large deflections of axisymmetric circular membranes International Journal of Solids and Structures 1971 7 12 1601 1612 2-s2.0-0015203744 Cheng C. J. Yang X. The problem of large deflections of rotationally orthotropic circular membranes Communication on Applied Mathematics and Computation 1987 1 1 38 47 Chen S.-L. Zheng Z.-L. Large deformation of circular membrane under the concentrated force Applied Mathematics and Mechanics 2003 24 1 28 31 2-s2.0-0038545560 Kang S. W. Free vibration analysis of composite rectangular membranes with a bent interface Journal of Sound and Vibration 2004 272 1-2 39 53 2-s2.0-1942485328 Young L. G. Ramanathan S. Hu J. Pai P. F. Numerical and experimental dynamic characteristics of thin-film membranes International Journal of Solids and Structures 2005 42 9-10 3001 3025 2-s2.0-12244265524 10.1016/j.ijsolstr.2004.09.031 Soares R. M. Gonalves P. B. Nonlinear vibrations and instabilities of a stretched hyperelastic annular membrane International Journal of Solids and Structures 2012 49 3-4 514 526 2-s2.0-84655161568 10.1016/j.ijsolstr.2011.10.019 Qian G. Z. The solutions of the free vibration problem of planar membrane with two different direction of force Applied Mathematics and Mechanics 1982 3 6 817 824 10.1007/BF01895336 Wu X. Natural vibration analysis of films with middle supports Journal of Vibration and Shock 2006 25 1 140 141 Lin W. J. Chen S. H. Free vibration analysis of annular membranes Journal of Vibration and Shock 2008 47 2 103 108 Ou yang C. The solvability study for a vibrating membrane robin problem with perturbed boundary Mathematics in Practice and Theory 2009 39 23 184 185 Lin W. J. Chen S. H. Free vibration analysis of plane membranes by finite element method Journal of Dynamics and Control 2010 8 3 202 206 Qiao L. Tan F. Yang Q.-S. Dynamic analysis of membrane structures Journal of Vibration and Shock 2011 30 6 109 113 2-s2.0-79960123573 Zheng Z. Song W. Liu C. Xu Y. Long J. Free vibration analysis of orthotropic circular membranes in large deflection Advanced Materials Research 2011 255–260 1730 1734 2-s2.0-79959801615 10.4028/www.scientific.net/AMR.255-260.1730 Lin W.-J. Chen S.-H. Free vibration analysis of circular membranes Journal of Vibration and Shock 2009 28 5 84 86 Li J. J. Large deflection of orthotropic thin circular plates Journal of Shanghai University (Natural Science Edition) 1997 3 1 9 18 Wang J. Nonlinear free vibration of the circular plate with large deflection Journal of South China University of Technology (Natural Science Edition) 2001 29 8 4 6 2-s2.0-0035416019