Natural Vibrations of a Spherical Dome

Let us consider a spherical dome of radius R, clamped along the contour (see figure).

Natural Vibrations of a Spherical Dome

Thickness of wall of the dome h is considerably smaller than its radius R.
Only one quarter of spherical surface will be considered. The bottom edge is fully restrained, the symmetry boundary conditions are applied to the side faces.

Natural Vibrations of a Spherical Dome, the finite element model with restraints

The finite element model with restraints

Let us use the following data: radius R = 300 mm, thickness h = 3 mm (R / h = 100).
The material properties are: the Young's modulus E = 2.1E+011 , Poisson's ratio ν=0.28, the density ρ = 7800 kg / m3.
Analytical solution of this problem is given by:

fi= ki . ω0 / 2π

,

where E – Young’s modulus, ki – coefficient whose value for the first five natural frequencies is: 0.5457, 0.7377, 0.8563, 0.8598, 0.9034.
Thus, f1 = 1564.7 Hz , f2 = 2115.3 Hz , f3 = 2455.4 Hz, f4 = 2465.4 Hz f5 = 2590.4 Hz.

After carrying out calculation with the help of AutoFEM, the following results are obtained*:

Table 1. Parameters of finite element mesh

Finite Element Type

Number of nodes

Number of finite elements

quadratic tetrahedron

21784

10762

Table 2. Result "Frequency"*

 

Numerical solution
Frequency fi*, Hz

Analytical solution
Frequency fi, Hz

Error δ = 100%*| fi* - fi| / | fi |

1

1569.065

1564.7

0.28

2

2107.329

2115.3

0.38

3

2458.378

2455.4

0.12

4

2481.949

2465.4

0.67

5

2575.043

2590.4

0.59

Natural Vibrations of a Spherical Dome, modes of vibration

Conclusion:

The relative error of the numerical solution compared with the analytical solution not exceed 1%.

*The results of numerical tests depend on the finite element mesh and may differ slightly from those given in the table.

 

 

Read more about AutoFEM Frequency Analysis

autofem.com

Return to contents