Finite Element Simulation of Axisymmetric Elastic and Electroelastic Wave Propagation Using Local-Domain Wave Packet Enrichment

2021 ◽  
Vol 144 (2) ◽  
Author(s):  
Amit Kumar ◽  
Santosh Kapuria
Keyword(s):  
Finite Element ◽  
Wave Propagation ◽  
Wave Packet ◽  
Equations Of Motion ◽  
Time Integration ◽  
Coupled System ◽  
Integration Scheme ◽  
Local Domain ◽  
Complex Wave ◽  
Wave Modes

Abstract A local-domain wave packet enriched multiphysics finite element (FE) formulation is employed for accurately solving axisymmetric wave propagation problems in elastic and piezoelastic media, involving complex wave modes and sharp jumps at the wavefronts, which pose challenges to the conventional FE solutions. The conventional Lagrangian interpolations for the displacement and electric potential fields are enriched with the element-domain sinusoidal functions that satisfy the partition of unity condition. The extended Hamilton’s principle is employed to derive the coupled system of equations of motion which is solved using the simple Newmark-β direct time integration scheme without resorting to any remeshing near the wavefronts or post-processing. The performance of the enrichment is assessed for the axisymmetric problems of impact waves in elastic and piezoelectric cylinders and elastic half-space, bulk and Rayleigh waves in the semi-infinite elastic domain and ultrasonic Lamb wave actuation and propagation in plate-piezoelectric transducer system. The element shows significant improvement in the computational efficiency and accuracy over the conventional FE for all problems, including those involving multiple complex wave modes and sharp discontinuities in the fields at the wavefronts.

scholarly journals Numerical Integration of Multibody Dynamic Systems Involving Nonholonomic Equality Constraints

2021 ◽  
Author(s):  
Sotirios Natsiavas ◽  
Panagiotis Passas ◽  
Elias Paraskevopoulos
Keyword(s):  
Dynamic Systems ◽  
Equations Of Motion ◽  
Time Integration ◽  
Nonholonomic Constraints ◽  
Coupled System ◽  
Weak Form ◽  
Integration Scheme ◽  
Motion Constraints ◽  
Multibody Dynamic ◽  
Time Integration Scheme

Abstract This work considers a class of multibody dynamic systems involving bilateral nonholonomic constraints. An appropriate set of equations of motion is employed first. This set is derived by application of Newton’s second law and appears as a coupled system of strongly nonlinear second order ordinary differential equations in both the generalized coordinates and the Lagrange multipliers associated to the motion constraints. Next, these equations are manipulated properly and converted to a weak form. Furthermore, the position, velocity and momentum type quantities are subsequently treated as independent. This yields a three-field set of equations of motion, which is then used as a basis for performing a suitable temporal discretization, leading to a complete time integration scheme. In order to test and validate its accuracy and numerical efficiency, this scheme is applied next to challenging mechanical examples, exhibiting rich dynamics. In all cases, the emphasis is put on highlighting the advantages of the new method by direct comparison with existing analytical solutions as well as with results of current state of the art numerical methods. Finally, a comparison is also performed with results available for a benchmark problem.

A Smoothing Scheme for Seismic Wave Propagation Simulation with a Mixed FEM of Diffusionized Wave Equation

2021 ◽  
pp. 2150061
Author(s):  
Ryuta Imai ◽  
Naoki Kasui ◽  
Masayuki Yamada ◽  
Koji Hada ◽  
Hiroyuki Fujiwara
Keyword(s):  
Wave Propagation ◽  
Wave Equation ◽  
Seismic Wave ◽  
Time Integration ◽  
Coupled System ◽  
Seismic Wave Propagation ◽  
Implicit Time ◽  
Integration Scheme ◽  
Weak Formulation ◽  
Mixed Fem

In this paper, we propose a smoothing scheme for seismic wave propagation simulation. The proposed scheme is based on a diffusionized wave equation with the fourth-order spatial derivative term. So, the solution requires higher regularity in the usual weak formulation. Reducing the diffusionized wave equation to a coupled system of diffusion equations yields a mixed FEM to ease the regularity. We mathematically explain how our scheme works for smoothing. We construct a semi-implicit time integration scheme and apply it to the wave equation. This numerical experiment reveals that our scheme is effective for filtering short wavelength components in seismic wave propagation simulation.

A Newly Developed Algorithm for Treating the Coupled Dynamic Response of Robotic Fixtureless Assembly

1996 ◽  
Author(s):  
Genady Shagal ◽  
Shaker A. Meguid
Keyword(s):  
Finite Element ◽  
Dynamic Response ◽  
Equations Of Motion ◽  
Time Integration ◽  
Integration Scheme ◽  
The Finite Element Method ◽  
Coupled Equations ◽  
Time Integration Scheme ◽  
Implicit Approach ◽  
Cooperating Robots

Abstract The coupled dynamic response of two cooperating robots handling two flexible payloads for the purpose of fixtureless assembly and manufacturing is treated using a new algorithm. In this algorithm, the equations describing the dynamics of the system are obtained using Lagrange’s method for the rigid robot links and the finite element method for the flexible payloads. A new time integration scheme is developed to treat the coupled equations of motion of the rigid links for a given displacement of the flexible payloads. The finite element equations of the flexible payloads are then treated using an implicit approach. The new algorithm was verified using simplified examples and was later used to examine the dynamic response of two cooperating robot arms manipulating flexible payloads which are typical of the automotive industry.

scholarly journals Evaluation of Component Mode Synthesis Methods for the Detection of Modal Interaction Through Rotor Stator Contacts

2009 ◽  
Author(s):  
Alain Batailly ◽  
Mathias Legrand ◽  
Patrice Cartraud ◽  
Christophe Pierre ◽  
Jean-Pierre Lombard
Keyword(s):  
Finite Element ◽  
Equations Of Motion ◽  
Time Integration ◽  
Component Mode Synthesis ◽  
Integration Scheme ◽  
Synthesis Methods ◽  
Reduced Order Models ◽  
Reduced Order ◽  
Free Interface ◽  
Fixed Interface

The study of interactions through direct contact between bladetips and outer casings in modern turbomachines may be very time-consuming when the classical finite element method is used. The construction of reduced-order models using component mode synthesis (CMS) methods generally allows for dramatic increase in computational efficiency and may be used in order to improve the knowledge over these interaction phenomena. Among the available approaches, both a fixed-interface method and a free-interface method are considered here in an original manner to reduce the size of a realistic two-dimensional model. The equations of motion are solved using an explicit time integration scheme with the Lagrange multiplier method where friction is accounted for. This method offers energy momentum conserving which is a critical point to ensure the convergence of the algorithm. Moreover, it is shown that even in a non-linear framework the reduced-order models converge to the finite element solution as the number of modes included in the models increases. Considering the fixed-interface method of Craig-Bampton (CB) and the free-interface method of Craig-Chang-Martinez (CCM), it is shown that a method with fast displacement convergence may be less efficient in terms of motion convergence.

scholarly journals Transient Wave Propagation Dynamics with Edge-Based Smoothed Finite Element Method and Bathe Time Integration Technique

2020 ◽  
Vol 2020 ◽  
pp. 1-16
Author(s):  
Yingbin Chai ◽  
Yongou Zhang
Keyword(s):  
Finite Element ◽  
Wave Propagation ◽  
Time Integration ◽  
Numerical Dispersion ◽  
Integration Scheme ◽  
Transient Wave ◽  
Time Integration Scheme ◽  
Smoothed Finite Element Method ◽  
Transient Wave Propagation ◽  
Edge Based

In this work, the edge-based smoothed finite element method (ES-FEM) is incorporated with the Bathe time integration scheme to solve the transient wave propagation problems. The edge-based gradient smoothing technique (GST) can properly soften the “overly soft” system matrices from the standard finite element approach; then, the spatial numerical dispersion error of the calculated solutions for wave problems can be significantly reduced. To effectively solve the transient wave propagation problems, the Bathe time integration scheme is employed to perform the involved time integration. Due to the appropriate “numerical dissipation effects” from the Bathe time integration method, the spurious oscillations in the relatively large wave numbers (high frequencies) can be effectively suppressed; then, the temporal numerical dispersion error in the calculated solutions can also be notably controlled. A number of supporting numerical examples are considered to examine the capabilities of the present approach. The numerical results show that ES-FEM works very well with the Bathe time integration technique, and much more numerical solutions can be reached for solving transient wave propagation problems compared to the standard FEM.

An Efficient Numerical Method for Dynamic Interaction Analysis of High-Speed Shinkansen Vehicles, Rail, and Bridge

1994 ◽  
Author(s):  
Makoto Tanabe ◽  
Hajime Wakui ◽  
Nobuyuki Matsumoto
Keyword(s):  
Finite Element ◽  
Numerical Method ◽  
High Speed ◽  
Interaction Analysis ◽  
Equations Of Motion ◽  
Time Integration ◽  
Integration Scheme ◽  
Element Formulation ◽  
Exact Time ◽  
Time Integration Scheme

Abstract This paper describes a finite element formulation to solve for the combined dynamic behavior of Shinkansen (bullet train) vehicles, irregular rails, and bridges. A mechanical model for interactions between a wheel and an irregular rail is discussed. The bridge is modeled by use of various finite elements. An efficient numerical method, based on modal analysis and exact time integration, is described for solving the nonlinear equations of motion of the Shinkansen vehicle and bridge. The convergence of the exact time integration scheme is discussed and compared with a previous numerical time integration scheme. A finite element computer program has been developed to analyze the dynamic response of Shinkansen vehicles operating at high speed over irregular rails and a bridge. Numerical examples are presented to demonstrate the effectiveness and validity of the present approach.

Contribution of Lamb wave modes in the formation of higher order modes cluster (HOMC) guided waves

2021 ◽  
pp. 095440622110424
Author(s):  
Z Abbasi ◽  
F Honarvar
Keyword(s):  
Finite Element ◽  
Wave Packet ◽  
Lamb Wave ◽  
Guided Waves ◽  
Element Model ◽  
Higher Order ◽  
Guided Wave ◽  
Cross Sectional ◽  
Higher Order Modes ◽  
Wave Modes

In recent years, Higher Order Modes Cluster (HOMC) guided waves have been considered for ultrasonic testing of plates and pipes. HOMC guided waves consist of higher order Lamb wave modes that travel together as a single nondispersive wave packet. The objective of this paper is to investigate the effect of frequency-thickness value on the contribution of Lamb wave modes in an HOMC guided wave. This is an important issue that has not been thoroughly investigated before. The contribution of each Lamb wave mode in an HOMC guided wave is studied by using a two-dimensional finite element model. The level of contribution of various Lamb wave modes to the wave cluster is verified by using a 2D FFT analysis. The results show that by increasing the frequency-thickness value, the order of contributing modes in the HOMC wave packet increases. The number of modes that comprise a cluster also increases up to a specific frequency-thickness value and then it starts to decrease. Plotting of the cross-sectional displacement patterns along the HOMC guided wave paths confirms the shifting of dominant modes from lower to higher order modes with increase of frequency-thickness value. Experimental measurements conducted on a mild steel plate are used to verify the finite element simulations. The experimental results are found to be in good agreement with simulations and confirm the changes observed in the level of contribution of Lamb wave modes in a wave cluster by changing the frequency-thickness value.

Transonic Flow Computations in Cascades Using Finite Element Method

1981 ◽  
Vol 103 (4) ◽  
pp. 657-664 ◽  
Author(s):  
H. U. Akay ◽  
A. Ecer
Keyword(s):  
Finite Element Method ◽  
Finite Element ◽  
Transonic Flow ◽  
Method Development ◽  
Time Integration ◽  
Integration Scheme ◽  
Complex Flow ◽  
Time Integration Scheme ◽  
Naca 0012 ◽  
Element Method

Analysis of transonic flow through a cascade of airfoils is investigated using the finite element method. Development of a computational grid suitable for complex flow structures and different types of boundary conditions is presented. An efficient pseudo-time integration scheme is developed for the solution of equations. Modeling of the shock and the convergence characteristics of the developed scheme are discussed. Numerical results include a 45 deg staggered cascade of NACA 0012 airfoils with inlet flow Mach number of 0.8 and angles of attack 1, 0, and −1 deg.

Time Integration of Impact Phenomena in Solid Mechanics: A One-Dimensional Model Problem

1996 ◽  
Author(s):  
V. Chawla ◽  
T. A. Laursen
Keyword(s):  
Finite Element ◽  
Solid Mechanics ◽  
Exact Analytical Solution ◽  
Time Integration ◽  
Dynamic Contact ◽  
Momentum Conservation ◽  
Integration Scheme ◽  
Mass System ◽  
Element Discretization ◽  
Time Integration Scheme

Abstract 1D impact between two identical bars is modeled as a simple spring-mass system as would be generated by a finite element discretization. Some commonly used time integrators are applied to the system to demonstrate defects in the numerical solution as compared to the exact analytical solution. Using energy conservation as the criterion for stability, a new time integration scheme is proposed that imposes a persistency condition for dynamic contact. Finite element simulation with Lagrange multipliers for enforcing the contact constraints shows exact energy and momentum conservation.

Plain and Full Floating Bearing Simulations With Rigid Shaft Dynamics

2007 ◽  
Author(s):  
Ian McLuckie ◽  
Scott Barrett
Keyword(s):  
System Development ◽  
Equations Of Motion ◽  
Time Integration ◽  
Forced Response ◽  
Design Parameters ◽  
Integration Scheme ◽  
Damping Coefficients ◽  
Time Integration Scheme ◽  
Combined Solution ◽  
Stiffness And Damping

This paper shows a promising predictive bearing model that can be used to reduce turbocharger bearing system development times. Turbocharger development is normally done by varying design parameters such as bearing geometry in a very time consuming experimentation process. Full Floating Bearings (FFB) are used in most automotive turbochargers and, due to emissions regulations, there has been a push towards downsizing engines and applying turbo charging to generate optimized engine solutions for both gasoline and diesel applications. In this paper the turbocharger rotor is regarded as being rigid, and the equations of motion are solved using the Bulirsch Stoer time integration scheme. These equations are solved simultaneously with the bearing model which is used also to determine nonlinear stiffness and damping coefficients. The bearings are solved using a Rigid Hydro Dynamic (RHD) Finite Difference Successive Over Relaxation (SOR) scheme of Reynolds equation that includes both rotational and squeeze velocity terms. However the solver can also consider bearing and rotor elasticity in a Multi-Body Dynamic (MBD) and Elasto-Hydro Dynamic (EHD) combined solution. Two bearing types have been studied, a plain grooved (PGB) and a full floating bearing (FFB) for comparative purposes. The mathematical models used are generic and suitable for whole engine bearing studies. The results in this paper show they are suitable for determining the onset of turbocharger bearing instability, and also the means by which bearing instability may be suppressed. The current study has investigated forced response with the combined effects of gravity and unbalance. It is worth noting that the effects of both housing excitation and aerodynamic excitation from the compressor and turbine can be easily accommodated, and will be the subject of a future paper. Other topics introduced here that will be explored further in the future include the effect of bearing and rotor flexibility in the MBD and EHD solution and the use of automatically generated stiffness and damping coefficients for any bearing geometry.

Sign in / Sign up

Export Citation Format

Share Document