<figure><img src="image_1.png" /></figure>

HAL
open science

# Free wave propagation in two-dimensional periodic beam grillage

Changwei Zhou, Jean-Pierre Lainé, Mohamed Ichchou, Abdelmalek Zine

► To cite this version:

Changwei Zhou, Jean-Pierre Lainé, Mohamed Ichchou, Abdelmalek Zine. Free wave propagation in two-dimensional periodic beam grillage. 21st International Congress on Sound and Vibration, Jul 2014, Beijing, China. hal-01054912

HAL Id: hal-01054912

https://hal.science/hal-01054912

Submitted on 3 Nov 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# FREE WAVE PROPAGATION IN TWO-DIMENSIONAL PERIODIC BEAM GRILLAGE

Changwei Zhou, Mohamed Ichchou, Jean-Pierre Lainé

Laboratoire de Tribologie et Dynamique des Systèmes, Ecole Centrale de Lyon. 36, Avenue Guy de Collongue, 69134 Ecully, France

e-mail: chang-wei.zhou@ec-lyon.fr

Abdel-Malek Zine

Institut Camille Jordan, Ecole Centrale de Lyon. 36, Avenue Guy de Collongue, 69134 Ecully, France

Two-dimensional periodic structures are widely applied in the engineering domain, such as the orthogonally rib-stiffened panels used in the fuselage of an aircraft or the cellular structures like honeycombs used in sandwich construction. The dynamic behaviour of these structures has drawn a lot of attention among researchers for the last decades. In this paper, A beam grillage is studied as a particular example of a two-dimensional periodic structure. A single unit cell is modelled by the finite element method and periodic structure theory is employed to study the whole beam grillage. The method is known as the Wave Finite Element method (WFEM), which is usually used to study the one-dimensional waveguide or two-dimensional homogeneous waveguide. This work has extended the WFEM to the two-dimensional periodic beam grillage to calculate the wave dispersion characteristics. It also extended the model reduction based on Craig-Bampton method in WFEM to two-dimensional periodic structures. The simulation has been validated by the theoretical result in the literature. It is also compared with the dispersion relation by using the Fourier transformation of the displacement field under harmonic excitation. The WFEM shows a good agreement with little computational cost. The fixed frequency dispersion relation deduced from WFEM is used to predict the direction of wave propagation and to investigate the existence of direction where waves cannot propagate in the beam grillage.

## 1. Introduction

Periodic structures which can be obtained by repeating a single unit cell occur largely in engineering domain. For example an orthogonally stiffened plate or shell used in the fuselage of an aircraft or a ship hull. In this case the structure is called two-dimensional periodic because the whole structure can be formed by repeating the unit cell in two directions. This paper use the WFEM to analyse the two-dimensional periodic structure, in particular a two-dimensional periodic beam grillage. The method is usually used to model the one-dimensional waveguide or two-dimensional homogeneous waveguide [1, 2]. In this work, WFEM has been extended to study the two-dimensional periodic beam grillage. Similar to the one-dimensional waveguide, the internal nodes in an unit cell of two-dimensional periodic structures can be numerous, so it is interesting to apply the model reduction as
it is used in one-dimensional periodic case [3]. The result of WFEM on two-dimensional structures (WFE2DM) can be presented with the fixed frequency contour of dispersion relation, which can be used to determine the direction of wave propagation in the considered structures and in the investigation of the existence of directions within the structure, where waves do not propagate. It is called “wave beaming effect” where local vibrations are of primary concern [4, 5, 6]. The wave beaming effect of beam grillage is studied in this paper by both the dispersion relation and the displacement of an finite beam grillage under excitation.

## 2. Model reduction on the unit cell based on Craig-Bampton method

The structure to be modelled is periodic in both the x and y directions, and properties can also vary through its thickness in the z direction as in the laminate structures. The unit cell which forms the structures is discretized by the Finite Element (FE) method. The mass and stiffness matrices M and K of the unit cell are found using commercial FE package. All the boundary nodes are denoted $q_{bd}$, $q_I$ represents the internal DOFs. So the nodal DOFs of the unit cell become $((q_{bd})^T(q_I)^T)^T$. For free wave propagation, no external forces act on the structures apart from those on its boundaries from adjacent unit cells. So the equation of motion of the unit cell, assuming time-harmonic behaviour and neglecting damping, becomes:

$$[\mathbf{K} - \omega^2 \mathbf{M}] \begin{pmatrix} \mathbf{q}_{bd} \\ \mathbf{q}_I \end{pmatrix} = \begin{pmatrix} \mathbf{f} \\ \mathbf{0} \end{pmatrix} \qquad (1)$$

Matrix $[\mathbf{K} - \omega^2 \mathbf{M}]$ is the dynamic stiffness matrix of the unit cell, with the presence of the damping, it becomes complex and can be expressed by $[\mathbf{K} + i\mathbf{K}' - \omega^2 \mathbf{M}]$. To model the dynamic behaviour in high frequency, the unit cell is often discretized with numerous number of internal DOFs. A reduced model of the unit cell is necessary so as to gain the computation time. In this approach, Craig-Bampton hybrid coordinates $((q_{bd})^T(P_c)^T)$ are related to the physical coordinates $((q_{bd})^T(q_I)^T)$ as follows:

$$\begin{pmatrix} \mathbf{q}_{bd} \\ \mathbf{q}_I \end{pmatrix} = \begin{bmatrix} \mathbf{I}_n & \mathbf{0} \\ \mathbf{\Psi}_{bd} & \mathbf{\Psi}_c \end{bmatrix} \begin{pmatrix} \mathbf{q}_{bd} \\ \mathbf{P}_c \end{pmatrix} \qquad (2)$$

With $[\mathbf{\Psi}_{bd}]$ represent the static boundary modes, $[\mathbf{\Psi}_c]$ represent the fixed base mode shapes. The approach has been used in one-dimensional WFEM, which shows that if the fixed base mode shapes under natural frequency lower than three time the maximum investigated frequency are retained in the modal basis, the WFEM on reduced model is able to give precise result [3]. The reduced model with less internal DOFs is easier to be applied dynamic condensation. The Eq. (1) can be represented by reduced dynamic stiffness matrix with only the boundary DOFs:

$$D\mathbf{q}_{bd} = \mathbf{f} \qquad (3)$$

Once the reduced dynamic stiffness matrix of an unit cell is obtained, the WFE2DM can be applied to study the wave propagation characteristics of the periodic structures.

## 3. Wave Finite Element method applied on two-dimensional periodic beam grillage

A beam grillage is a simple form of continuum two-dimensional periodic structure. This section explains how WFE2DM can be applied to study the beam grillage.
## 3.1 WFE2DM on periodic beam grillage

The beam grillage consists of beams in two directions which intersect on the joints. An unit cell can be used to represent the beam grillage. The one used in this paper is as shown in Fig. 1, which contains two beams and 1/4 mass at the joints and 1/2 mass at the other end of the beams. So the whole beam grillage with (2*1/4+1/2)=1 mass at the joints can be obtained by repeating this unit cell. The unit cell consists of [q₁, q₂, q₃, q₄], so after the dynamic condensation of the internal DOFs, the

<figure><img src="image_4.png" /><figcaption>Figure 1. Unit cell used in WFE2DM</figcaption></figure>

equilibrium equation becomes:

$$ \begin{bmatrix} D_{11} & D_{12} & D_{13} \\ D_{21} & D_{22} & D_{23} \\ D_{31} & D_{32} & D_{33} \end{bmatrix} \begin{pmatrix} q_1 \\ q_2 \\ q_3 \end{pmatrix} = \begin{pmatrix} f_1 \\ f_2 \\ f_3 \end{pmatrix} \qquad (4) $$

Using the periodic structure theory, the displacement in the unit cell can be related as follows:

$$ \begin{pmatrix} q_1 \\ q_2 \\ q_3 \end{pmatrix} = \begin{pmatrix} I_s \\ \lambda_x I_s \\ \lambda_y I_s \end{pmatrix} q_1 = \Lambda_R q_1 \qquad (5) $$

The relation between λ and wavenumber k is as follows:

$$ \lambda_x = e^{-ik_x L_x}, \lambda_y = e^{-ik_y L_y} \qquad (6) $$

The sum of the nodal forces of all elements connected to the nodes 1 is zero, so

$$ ( I_s \quad \lambda_x^{-1} I_s \quad \lambda_y^{-1} I_s ) f_1 = \Lambda_L f_1 = 0 \qquad (7) $$

So the equilibrium equation becomes:

$$ \Lambda_L * D * \Lambda_R = 0 \qquad (8) $$

Suppose the frequency ω and one between (λₓ, λᵧ) is known, for example λᵧ. This may represent the situation where a known wave is incident on a straight boundary so that λᵧ is given and all possibly solutions for λₓ are sought which correspond to evanescent or propagation waves. λₓ can be calculated by solving the following eigenvalue problem

$$ \frac{1}{\lambda_x} (A_2 \lambda_x^2 + A_1 \lambda_x + A_0) = 0 \qquad (9) $$
with

$$A_0 = D_{23}\lambda_y + D_{21} \qquad (10a)$$

$$A_1 = D_{11} + D_{22} + D_{33} + D_{13}\lambda_y + D_{31}\lambda_y^{-1} \qquad (10b)$$

$$A_2 = D_{12} + D_{32}\lambda_y^{-1} \qquad (10c)$$

Then the quadratic eigenvalue problem on $\lambda_x$ can be solved, which contain 2s solutions, s is the number of the DOFs of $q_1$. The dispersion relation is presented as an iso-frequency contour in coordinate system $(k_x, k_y)$, then the Poynting vector (i.e., the energy flow vector) is normal to the contour[7].

## 3.2 Validation of the WFE2DM by the result in the reference

Langley and Bardell have done the experimental test on a beam grillage constructed from strips of aluminium with bolted joints[4]. The same structure is used here in order to compare the results by WFE2DM with his results. The periodic grillage structure is constructed from an orthogonal array of aluminium strips as shown in 2. The strips were of rectangular cross-section with 19.06 mm × 1.59mm. At the crossing points, the joints are composed of a set of washers, nut and bolt, which is of $6 \times 10^{-3}$ kg. The joints are taken account as the mass point.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** A schematic of the experimental structure by [4]. The excitation is at position (A). All dimensions are in mm.</figcaption></figure>

Unit cell used in WFE2DM is defined in the precedent subsection. In the unit cell, the mass point is modelled by Mass21 and the beam is modelled by Beam188 in ANSYS. The latter is a three-dimensional two-node Timoshenko element, with 6 DOFs per node. Since the cross-section of the element is rectangular, the movements of the 6 DOFs are decoupled. The present analysis is restricted to the out-of-plane bending of the beam components so the DOFs related to in-plane movement are fixed to zero, subsequently only rotx, roty and uz are conserved in each node. To assure the convergence of the problem, the two beams are both meshed with 100 elements. So the $q_1$, $q_2$ as well as $q_3$ is of 3 DOFs, the internal DOFs is of 594 DOFs. In the reduced model, only 8 internal modes which represent the natural frequency under 4500 Hz are taken account in. The size
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Fixed frequency contours of dispersion relation, – Langley, \* WFE2DM on completed model, o WFE2DM on reduced model</figcaption></figure>

of the reduced model is about 1% of the full model. The computation time decreases a lot to apply the dynamic condensation on the reduced model.

The dispersion relations at 100Hz, 600Hz, 800Hz, 1000Hz, 1200Hz, 1400Hz are given in the Fig. 3. The dispersion relation obtained by the proposed WFE2DM conforms with the theoretic result in the reference [4]. The method applied on the reduced model can predict the dispersion relation with less computation time.

## 3.3 The wave beaming effects in beam grillage

The dispersion relation can also be calculated by the two-dimensional Discrete spatial Fourier Transformation (DFT2D). Here the Finite Element model based on ANSYS which consists of 100 unit cells in x direction and 100 unit cells in y direction. An harmonic excitation is applied in the center of the grillage, then the out of plane displacement is calculated and then DFT2D is applied to obtain the dispersion relation at the frequency of the excitation. At 100Hz, the normal directions of the fixed frequency contour contain all the directions. The conclusion is confirmed by the displacement field. However at higher frequency, such as 600Hz, 800Hz and 1000Hz, as shown in the Fig.4, Fig.5 and Fig.6 respectively, the directivity of the propagation can be observed. It exists some “dead zones” for directions. The normal directions of the fixed frequency contours are restrained to certain angle. So is the displacement field under harmonic excitation at these frequencies.

## 4. Conclusion

In this paper we propose the WFE2DM to study a two-dimensional periodic structures such as beam grillage. The model reduction based on Craig-Bampton is applied to decrease computation time and proved to give as accurate result as the full model. The fixed frequency contours of dispersion relation are obtained by the WFE2DM and are validated by the theoretical result existing in the reference. The wave beam effects in the beam grillage are studied, by using the fixed frequency contour of dispersion relation and by numerical simulation of the finite beam grillage. There is a good correlation between the results obtained by the two methods, it exists some dead directions of the propagation while other directions are privileged. Some experimental studies on the two-dimensional periodic structure such as beam grillage will be considered in the future work.
<figure><img src="image_7.png" /></figure>

**Figure 4.** At 600Hz (a) Phase constant surface by the DFT2D, the line with 'o' indicates the result by WFE2DM on reduced model. (b) The displacement field of the beam grillage with 100 periods in each direction

<figure><img src="image_8.png" /></figure>

**Figure 5.** At 800Hz (a) Phase constant surface by the DFT2D, the line with 'o' indicates the result by WFE2DM on reduced model. (b) The displacement field of the beam grillage with 100 periods in each direction
<figure><img src="image_9.png" /><figcaption>(a) Fixed frequency contour of dispersion relation</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>(b) Displacement field of the beam grillage</figcaption></figure>

**Figure 6.** At 1000Hz (a) Phase constant surface by the DFT2D, the line with 'o' indicates the result by WFE2DM on reduced model. (b) The displacement field of the beam grillage with 100 periods in each direction

REFERENCES

1 Mencik, J. -M. and Ichchou, M. N. Multi-mode propagation and diffusion in structures through finite elements. *European Journal of Mechanics A/Solids*, **24**(5), 877-898, (2005).

2 Mace, B. R. and Manconi, E. Modelling wave propagation in two-dimensional structures using finite element analysis, *Journal of Sound and Vibration*, **318**(4-5),884-902.(2008)

3 Zhou, C. W., Laine, J. -P., Ichchou, M. N. and Zine, A. M. Wave finite element method based on reduced model for periodic structures, submitted to *International Journal of Solids and Structures*

4 Langley, R. S. and Bardell, N. S. The response of two-dimensional periodic structures to harmonic point loading: a theoretical and experimental study of a beam grillage, *Journal of Sound and Vibration*, **207**(4), 521-535, (1997).

5 Ruzzene, M., Scarpa, F. and Soranna, F. Wave beaming effects in two-dimensional cellular structures, *Smart Materials and Structures*, **12**(3),363-372,(2003).

6 Kohrs, T. and Björn, A.T. Petersson, Wave beaming and wave propagation in light weight plates with truss-like cores, *Journal of Sound and Vibration*, **321**(1-2),137-165,(2009).

7 Langley, R. S. On the modal density and energy flow characteristics of periodic structures, *Journal of Sound and Vibration*, **172**(2),491-511,(1994).