scholarly journals The Wave Finite Element Method Applied to a One-Dimensional Linear Elastodynamic Problem With Unilateral Constraints

Author(s):  
Carlos Yoong ◽  
Anders Thorin ◽  
Mathias Legrand

The Wave Finite Element Method (WFEM) is implemented to accurately capture travelling waves propagating at a finite speed within a bouncing rod system and induced by unilateral contact collisions with a rigid foundation; friction is not accounted for. As opposed to the traditional Finite Element Method (FEM) within a time-stepping framework, potential discontinuous deformation, stress and velocity wave fronts are accurately described, which is critical for the problem of interest. A one-dimensional benchmark with an analytical solution is investigated. The WFEM is compared to two time-stepping solution methods formulated on a FEM semi-discretization in space: (1) an explicit technique involving Lagrange multipliers and (2) a non-smooth approach implemented in the Siconos package. Attention is paid to the Gibb’s phenomenon generated during and after contact occurrences together with the time evolution of the total energy of the system. It is numerically found that the WFEM outperforms the FEM and Siconos solution methods because it does not induce any spurious oscillations or dispersion and diffusion of the shock wave. Furthermore, energy is not dissipated over time. More importantly, the WFEM does not require any impact law to close the system of governing equations.

2021 ◽  
pp. 107754632098131
Author(s):  
Jamil Renno ◽  
Sadok Sassi ◽  
Wael I Alnahhal

The prediction of the response of waveguides to time-harmonic base excitations has many applications in mechanical, aerospace and civil engineering. The response to base excitations can be obtained analytically for simple waveguides only. For general waveguides, the response to time-harmonic base excitations can be obtained using the finite element method. In this study, we present a wave and finite element approach to calculate the response of waveguides to time-harmonic base excitations. The wave and finite element method is used to model free wave propagation in the waveguide, and these characteristics are then used to find the amplitude of excited waves in the waveguide. Reflection matrices at the boundaries of the waveguide are then used to find the amplitude of the travelling waves in the waveguide and subsequently the response of the waveguide. This includes the displacement and stress frequency response transfer functions. Numerical examples are presented to demonstrate the approach and to discuss the numerical efficiency of the proposed method.


1994 ◽  
Vol 61 (4) ◽  
pp. 919-922 ◽  
Author(s):  
Taein Yeo ◽  
J. R. Barber

When heat is conducted across an interface between two dissimilar materials, theimoelastic distortion affects the contact pressure distribution. The existence of a pressure-sensitive thermal contact resistance at the interface can cause such systems to be unstable in the steady-state. Stability analysis for thermoelastic contact has been conducted by linear perturbation methods for one-dimensional and simple two-dimensional geometries, but analytical solutions become very complicated for finite geometries. A method is therefore proposed in which the finite element method is used to reduce the stability problem to an eigenvalue problem. The linearity of the underlying perturbation problem enables us to conclude that solutions can be obtained in separated-variable form with exponential variation in time. This factor can therefore be removed from the governing equations and the finite element method is used to obtain a time-independent set of homogeneous equations in which the exponential growth rate appears as a linear parameter. We therefore obtain a linear eigenvalue problem and stability of the system requires that all the resulting eigenvalues should have negative real part. The method is discussed in application to the simple one-dimensional system of two contacting rods. The results show good agreement with previous analytical investigations and give additional information about the migration of eigenvalues in the complex plane as the steady-state heat flux is varied.


Sign in / Sign up

Export Citation Format

Share Document