Lab Home | Phone | Search
Center for Nonlinear Studies  Center for Nonlinear Studies
 Colloquia Archive 
 Postdoc Seminars Archive 
 Quantum Lunch 
 Quantum Lunch Archive 
 CMS Colloquia 
 Q-Mat Seminars 
 Q-Mat Seminars Archive 
 P/T Colloquia 
 Kac Lectures 
 Kac Fellows 
 Dist. Quant. Lecture 
 Ulam Scholar 
 CNLS Fellowship Application 
 Student Program 
 Past Visitors 
 History of CNLS 
 Maps, Directions 
 CNLS Office 
Monday, June 02, 2014
3:00 PM - 4:00 PM
CNLS Conference Room (TA-3, Bldg 1690)


Advances in Computational Dynamics

Jiri Plesek
Institute of Thermomechanics, Prague

Spatial discretization of a continuum introduces dispersion error to the numerical solution of wave propagation. In the introduction, review is made of fundamental approaches used to derive the truncation error of the finite element method in an elastodynamic analysis, valid for elements with linear shape functions. In the second part, recent results accomplished by the author and his coworkers are summarized, namely the extension of dispersion theory to quadratic finite elements, following the lines of reasoning introduced by Belytschko and Mullen (1978) for one-dimensional elements and those of Abboud and Pinsky (1992), concerning the scalar Helmholtz equation.

Next, a detailed study follows of the mass matrix lumping schemes suitable for wave propagation simulations. To this end, a variable parameter is defined whose role is to distribute total mass between the element's corner and midside nodes. Based on that, dispersion analysis is carried out for varying value of this parameter. For example, it is shown in terms of dispersion curves that the Hinton-Rock-Zienkiewicz midside to corner mass ratio is far from optimum and the most accurate representation is obtained when 92% of total mass is coalesced into four midside nodes leaving only 8% share to the corner nodes.

Second, the issue of explicit time integration is addressed. Rigorous stability theory is invoked to assess the stability properties of the central difference method (as a typical representative), employing both the derived dispersion curves to impose the upper bound on the critical time step as well as the complementary estimate offered by the Fried theorem, which establishes the lower bound on the system eigenvalues.

This section closes with two numerical examples. In the first example, an infinite 2D space loaded by a point-wise source is considered to test the spatial finite element discretization. The loading frequency gradually changes from zero to high values to mimic the dispersion response to a broad loading spectrum. With the second example, the analytical solution to the longitudinal impact of two cylindrical bars as in the split-Hopkinson pressure bar test is employed to gauge dispersion for a contact-impact problem defined by the eight-node serendipity elements. Both examples results show superior agreement with the developed theory.

The main difficulty in the contact-impact analysis is a non-smoothness of contacting surfaces. It arises from inequality constraints as well as from geometric discontinuities induced by spatial discretization. The contact analysis based on traditional finite elements utilizes element facets to describe a contact surface. Unfortunately, the facet interfaces are only C0 continuous so that normal surface vectors may experience jump across element boundaries, which may lead to artificial oscillations of contact forces. A remedy to this geometric discontinuity may be provided by isogeometric formulation. In this approach, the known geometry is accurately described by the Non-Uniform Rational B-Splines (NURBS) basis functions. In this work, the mortar-KTS contact algorithm is utilized tohether with the central difference time integration scheme. The present algorithm is studied by means of a numerical example, which involves impact of two elastic tubes. The results clearly demonstrate the superiority of the NURBS discretization over the conventional Lagrange polynomial ansatz, provided the consistent matrix had been used, whereas considerable deterioration was observed for lumped mass matrices.

Host: Ondrej Certik