Plane wave finite element model for the 2-D phononic crystal under force loadings

Direct numerical simulations of the dynamic response of the phononic crystal to point or distributed forces via conventional numerical methods are still formidable currently due to large number of DOFs involved in the simulations. In this study, based on the plane wave expansion (PWE) method, a plan...

Full description

Saved in:
Bibliographic Details
Published in:European journal of mechanics, A, Solids A, Solids, 2022-01, Vol.91, p.104426, Article 104426
Main Authors: Lu, Jian-Fei, Cheng, Jie, Feng, Qing-Song
Format: Article
Language:English
Subjects:
Computer simulation
Crystals
Direct numerical simulation
Dynamic response
Eigenvalues
Eigenvectors
Finite element method
Fourier transforms
Horizontal orientation
Mathematical models
Matrix methods
Numerical methods
Plane wave expansion
Plane waves
The eigenvalue problem
The phononic crystal
The point force
Transfer matrices
Two dimensional models
Wavelengths
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Direct numerical simulations of the dynamic response of the phononic crystal to point or distributed forces via conventional numerical methods are still formidable currently due to large number of DOFs involved in the simulations. In this study, based on the plane wave expansion (PWE) method, a plane wave finite element (PWFE) model for the 2-D phononic crystal which is infinite in the horizontal direction and finite in the vertical direction and subjected to a point harmonic force is proposed. To this end, the point harmonic force is decomposed into the wavenumber domain (WND) components by the Fourier transform first. The variables of the 2-D phononic crystal under the WND force component are expanded into a series of plane waves along the horizontal direction, and the variables associated with each plane wave are discretized by the FE method in the vertical direction. By using the virtual work principle, FE equations for the plane waves of the phononic crystal are established. With the FE equations for the plane waves, the eigenvalue equation for the displacement of the plane waves and wavenumber is formulated. Superposing the eigenvectors of the above eigenvalue equation yields the displacement for the plane waves due to the WND force component. Synthetization of the contributions from all the plane waves and inversion of the Fourier transform with respect to the wavenumber via the contour integration method generate the response of the phononic crystal to the point force. Comparison of present results with those due to the transfer matrix method validates the proposed model. Also, presented numerical results suggest that the dynamic response of the phononic crystal to a point harmonic force exhibits remarkable resonance phenomenon. •A plane wave finite element (PWFE) model for the 2-D phononic crystal subjected to a point harmonic force is proposed.•The nonlinear eigenvalue equation for the plane wave displacement and wavenumber is derived and linearized.•The model superposition method for constructing the plane wave displacement of the phononic crystal is proposed to avoid direct solution of the FEM equations for the plane wave displacement.•The contour integration method is used to perform the inversion of the Fourier transform with respect to the wavenumber.
ISSN:0997-7538
1873-7285
DOI:10.1016/j.euromechsol.2021.104426