Required field
Required field
Not a valid email address
Required field
Required field

Documentation

# Validation Case: Harmonic Analysis of a Straight Beam

This harmonic analysis of a straight beam validation case belongs to solid mechanics. It aims to validate the following parameter:

• Harmonic analysis solver in Simscale

The harmonic analysis results from SimScale were compared to the numerical results in SHLL101$$^1$$.

## Geometry

Find below the straight beam geometry used for this validation case: Figure 1: Geometry used in the present validation case

The dimensions of the geometry are tabulated below. Furthermore, the center point C is located at the coordinates (0, 0, 0):

## Analysis Type and Mesh

Tool Type: Code_Aster

Analysis Type: Harmonic analysis

Mesh and Element Types: A total of four meshes are used in this validation case. The standard algorithm was used to create the meshes for cases A and B. The table below contains further information about the cases:

Find below the tetrahedral mesh used for case A: Figure 2: First-order mesh used for case A, created with the Standard meshing algorithm.

## Simulation Setup

Material:

• Material behavior: Linear elastic
• $$(E)$$ Young’s modulus = 165.8 $$GPa$$
• $$(v)$$ Poisson’s ratio = 0.3
• $$(\rho)$$ Density = 13404.106 $$kg/m^3$$
• Rayleigh damping settings (for the cases that model damping):
• Alpha coefficient = 0.001 $$s$$
• Beta damping = 0 $$1/s$$

Boundary Conditions:

• Constraint:
• The free end of the cylinder at $$x$$ = 0 is a Fixed support
• The other free end of the cylinder, at $$x$$ = 10 $$m$$, receives a bending or a tension load, according to table 2.
• Tension load: Force boundary condition, with $$F_x$$ = 3000 $$N$$
• Bending load: Force boundary condition, with $$F_y$$ = 3000 $$N$$

## Reference Solution

The analytical solution for the displacements, velocity, and acceleration on the point (10, 0, 0) are presented in SHLL101$$^1$$.

Note

The point (10, 0, 0), where the reference results were calculated, is the center point of the annulus.

As this point is not within the solid walls of the pipe, a total of 4 points were created on the free end to assess the results: Figure 3: Points at the free end ($$x$$ = 10) where data was collected

The coordinates of the points are as follows:
– Left: (10, 0, 0.08944)
– Right: (10, 0, -0.08944)
– Top: (10, 0.08944, 0)
– Bottom: (10, -0.08944, 0)

SimScale results are the average of the values obtained at these four points. When tension is being applied, the results are evaluated in the x-direction. Otherwise, if a bending load is applied, the results are taken in the y-direction.

## Result Comparison

In Table 3, we compare SimScale results with the reference SHLL101$$^1$$, modeling tension load, and no damping.

Table 4 shows a similar comparison, but now for the cases with a bending load and no damping.

Now, comparing the cases with tension load and the Rayleigh damping model:

Lastly, the results for bending load and damping are compared in Table 6:

SimScale results show great agreement with the analytical solution for all configurations.

Inspecting the displacement magnitude for case B4 in the post-processor: Figure 4: Harmonic analysis results for case B4, showing the displacement contours (magnified by a factor of 100).

Last updated: July 22nd, 2021