Required field
Required field
Required field
Required field

Documentation

# Validation Case: Frequency Analysis of a Ring

This validation case belongs to solid mechanics. The aim of this test case is to validate the following parameters:

• Frequency analysis

The simulation results of SimScale were compared to the results from [SDLS109]$$^1$$.

## Geometry

The geometry used for the frequency analysis is as follows:

The ring has a length $$L$$ of 0.05 $$m$$, a thickness $$t$$ of 0.048 $$m$$, and a medium radius $$Rm$$ of 0.369 $$m$$.

## Analysis Type and Mesh

Tool Type: Code_Aster

Analysis Type: Frequency Analysis

Mesh and Element Types:

First and second order meshes were computed using the SimScale standard mesh algorithm:

A mesh independence study was also performed for the second order meshes and IRAM – Sorensen algorithm to ensure the optimal fineness parameter level.

## Simulation Setup

Material:

• Linear Elastic Isotropic:
• $$E =$$ 185 $$GPa$$
• $$\nu =$$ 0.3
• $$\rho =$$ 7800 $$kg.m^{-3}$$

Boundary Conditions:

• Constraints:
• Body is free in space.

Computing Algorithm:

The following available computing algorithms were compared in the different cases:

The main characteristics of each algorithm are summarized below [U4.52.02]$$^2$$:

• IRAM – Sorensen: Uses a sub-space decomposition method to compute the natural frequencies and modes. Suitable for real and complex, symmetrical or non-symmetrical matrices.
• Lanczos: Uses a sub-space decomposition method to compute the natural frequencies and modes. Suitable for real, symmetrical only matrices.
• Bathe – Wilson: Uses a sub-space decomposition method to compute the natural frequencies and modes. Suitable for real, symmetrical only matrices.

A fourth algorithm is also available in SimScale, called QZ. This algorithm suffers from high memory consumption, which limits its application to cases with less than 1000 degrees of freedom. Therefore, it is not suitable for the current validation case.

Note

Complex matrices appear in the frequency analysis of materials with frequency damping. As this model is not available in SimScale, the difference between the algorithms comes down to robustness and speed.

## Frequency Reference Solution

The reference solution is of numerical type, as developed in [SDLS109]$$^1$$. The solution is presented in terms of all the natural frequencies and their corresponding shapes in the frequency range [200, 800] $$Hz$$. This solution was achieved by a convergence analysis using hexahedral elements, and as reported in the reference, a precision of 5% of the computed frequencies is estimated.

The consulted reference solution is:

## Frequency Results Comparison

Below can be found the results of the mesh independence study. For each natural frequency (F1 through F8), the variation of the result (in percent) with respect to the previous solution is plotted against the number of nodes in the mesh. At the final, finer mesh, the solution precision is 0.6% or lower.

Comparison of computed natural frequencies with the reference solution for each case can be seen below:

Results are mesh dependent instead of algorithm dependent because all algorithms produce similar results. The difference between algorithms can be seen when looking at the running times. Cases using IRAM – Sorensen and Lanczos algorithm are much faster than Bathe – Wilson. The recommendation is then to stay with the default algorithm (IRAM – Sorensen), because of its known robustness.

Following are the referenced natural mode shapes as seen on the online post-processor:

References

• SDLS109 – Fréquences propres d’un anneau cylindrique épais – Code_Aster validation case
• [U4.52.02] – Opérateur CALC_MODES – Code_Aster utilization manual

Note

If you still encounter problems validating you simulation, then please post the issue on our forum or contact us.

Last updated: June 30th, 2020