The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

The time finite element as a robust general scheme for solving nonlinear dynamic equations including chaotic systems

Schemes that can be proven to be unconditionally stable in the linear context can yield unstable solutions when used to solve nonlinear dynamical problems. Hence, the formulation of numerical strategies for nonlinear dynamical problems can be particularly challenging. In this work, we show that time finite element methods because of their inherent energy momentum conserving property (in the case of linear and nonlinear elastodynamics), provide a robust time-stepping method for nonlinear dynamic equations (including chaotic systems). We also show that most of the existing schemes that are known to be robust for parabolic or hyperbolic problems can be derived within the time finite element framework; thus, the time finite element provides a unification of time-stepping schemes used in diverse disciplines. We demonstrate the robust performance of the time finite element method on several challenging examples from the literature where the solution behavior is known to be chaotic. (C) 2015 Elsevier Inc. All rights reserved.

Seismic characterization and dynamic site response of a municipal solid waste landfill in Bangalore, India

Seismic design of landfills requires an understanding of the dynamic properties of municipal solid waste (MSW) and the dynamic site response of landfill waste during seismic events. The dynamic response of the Mavallipura landfill situated in Bangalore, India, is investigated using field measurements, laboratory studies and recorded ground motions from the intraplate region. The dynamic shear modulus values for the MSW were established on the basis of field measurements of shear wave velocities. Cyclic triaxial testing was performed on reconstituted MSW samples and the shear modulus reduction and damping characteristics of MSW were studied. Ten ground motions were selected based on regional seismicity and site response parameters have been obtained considering one-dimensional non-linear analysis in the DEEPSOIL program. The surface spectral response varied from 0.6 to 2g and persisted only for a period of 1s for most of the ground motions. The maximum peak ground acceleration (PGA) obtained was 0.5g and the minimum and maximum amplifications are 1.35 and 4.05. Amplification of the base acceleration was observed at the top surface of the landfill underlined by a composite soil layer and bedrock for all ground motions. Dynamic seismic properties with amplification and site response parameters for MSW landfill in Bangalore, India, are presented in this paper. This study shows that MSW has less shear stiffness and more amplification due to loose filling and damping, which need to be accounted for seismic design of MSW landfills in India.

Global Response Sensitivity Analysis of Randomly Excited Dynamic Structures

The response of structural dynamical systems excited by multiple random excitations is considered. Two new procedures for evaluating global response sensitivity measures with respect to the excitation components are proposed. The first procedure is valid for stationary response of linear systems under stationary random excitations and is based on the notion of Hellinger's metric of distance between two power spectral density functions. The second procedure is more generally valid and is based on the l2 norm based distance measure between two probability density functions. Specific cases which admit exact solutions are presented, and solution procedures based on Monte Carlo simulations for more general class of problems are outlined. Illustrations include studies on a parametrically excited linear system and a nonlinear random vibration problem involving moving oscillator-beam system that considers excitations attributable to random support motions and guide-way unevenness. (C) 2015 American Society of Civil Engineers.

SMO: An Integrated Approach to Intra-array and Inter-array Storage Optimization

The polyhedral model provides an expressive intermediate representation that is convenient for the analysis and subsequent transformation of affine loop nests. Several heuristics exist for achieving complex program transformations in this model. However, there is also considerable scope to utilize this model to tackle the problem of automatic memory footprint optimization. In this paper, we present a new automatic storage optimization technique which can be used to achieve both intra-array as well as inter-array storage reuse with a pre-determined schedule for the computation. Our approach works by finding statement-wise storage partitioning hyper planes that partition a unified global array space so that values with overlapping live ranges are not mapped to the same partition. Our heuristic is driven by a fourfold objective function which not only minimizes the dimensionality and storage requirements of arrays required for each high-level statement, but also maximizes inter statement storage reuse. The storage mappings obtained using our heuristic can be asymptotically better than those obtained by any existing technique. We implement our technique and demonstrate its practical impact by evaluating its effectiveness on several benchmarks chosen from the domains of image processing, stencil computations, and high-performance computing.

Identification of nonlinear dynamic systems: Time-frequency filtering and skeleton curves

The nonlinear behavior varying with the instantaneous response was analyzed through the joint time-frequency analysis method for a class of S. D. O. F nonlinear system. A masking operator an definite regions is defined and two theorems are presented. Based on these, the nonlinear system is modeled with a special time-varying linear one, called the generalized skeleton linear system (GSLS). The frequency skeleton curve and the damping skeleton curve are defined to describe the main feature of the non-linearity as well. Moreover, an identification method is proposed through the skeleton curves and the time-frequency filtering technique.

Chaotic dynamic analysis of viscoelastic plates

The dynamic buckling of viscoelastic plates with large deflection is investigated in this paper by using chaotic and fractal theory. The material behavior is given in terms of the Boltzmann superposition principle. in order to obtain accurate computation results, the nonlinear integro-differential dynamic equation is changed into an autonomic four-dimensional dynamical system. The numerical time integrations of equations are performed by using the fourth-order Runge-Kutta method. And the Lyapunov exponent spectrum, the fractal dimension of strange attractors and the time evolution of deflection are obtained. The influence of geometry nonlinearity and viscoelastic parameter on the dynamic buckling of viscoelastic plates is discussed.

Dynamic Characteristics of Spherically Converging Detonation Waves

The spherically converging detonation wave was numerically investigated by solving the one-dimensional multi-component Euler equations in spherical coordinates with a dispersion-controlled dissipative scheme. Finite rate and detailed chemical reaction models were used and numerical solutions were obtained for both a spherical by converging detonation in a stoichiometric hydrogen-oxygen mixture and a spherically focusing shock in air. The results showed that the post-shock pressure approximately arises to the same amplitude in vicinity of the focal point for the two cases, but the post-shock temperature level mainly depends on chemical reactions and molecular dissociations of a gas mixture. While the chemical reaction heat plays an important role in the early stage of detonation wave propagation, gas dissociations dramatically affect the post-shock flow states near the focal point. The maximum pressure and temperature, non-dimensionalized by their initial value, are approximately scaled to the propagation radius over the initial detonation diameter. The post-shock pressure is proportional to the initial pressure of the detonable mixture, and the post-shock temperature is also increased with the initial pressure, but in a much lower rate than that of the post-shock pressure.

Comparison of the Dynamic Response of a Saturated Soil Between Two Models

By comparing the dynamic responses of saturated soil to Biot's and Yamamoto's models, the properties of the two models have be pointed out. First of all, an analysis has been made for energy loss of each model from the basic equations. Then the damping of elastic waves in coarse sand and fine sand with loading frequency and soil's parameters have been calculated and the representation of viscous friction and Coulomb friction in the two models has been concluded. Finally, the variations of loading wave damping and stress phase angles with water depth and soil's parameters have been obtained as loading waves range in ocean waves.

Saturated duration for dynamic structural plastic response and fundamental periods of elastic vibration

The relationship is determined between saturated duration of rectangular pressure pulses applied to rigid, perfectly plastic structures and their fundamental periods of elastic vibration. It is shown that the ratio between the saturated duration and the fundamental period of elastic vibration of a structure is dependent upon two factors: the first one is the slenderness or thinness ratio of the structure; and the second one is the square root of ratio between the Young's elastic modulus and the yield stress of the structural material. Dimensional analysis shows that the aforementioned ratio is one of the basic similarity parameters for elastic-plastic modeling under dynamic loading.

Dynamic behavior of a cylindrical crack in a functionally graded interlayer under torsional loading

In this paper the problem of a cylindrical crack located in a functionally graded material (FGM) interlayer between two coaxial elastic dissimilar homogeneous cylinders and subjected to a torsional impact loading is considered. The shear modulus and the mass density of the FGM interlayer are assumed to vary continuously between those of the two coaxial cylinders. This mixed boundary value problem is first reduced to a singular integral equation with a Cauchy type kernel in the Laplace domain by applying Laplace and Fourier integral transforms. The singular integral equation is then solved numerically and the dynamic stress intensity factor (DSIF) is also obtained by a numerical Laplace inversion technique. The DSIF is found to rise rapidly to a peak and then reduce and tend to the static value almost without oscillation. The influences of the crack location, the FGM interlayer thickness and the relative magnitudes of the adjoining material properties are examined. It is found among others that, by increasing the FGM gradient, the DSIF can be greatly reduced.

Generalization of Response Number for Dynamic Plastic Response of Shells Subjected to impulsive Loading

A dimensionless number, termed as response number in Zhao [Archive of Applied Mechanics 68 (1998) 524], has been suggested for the dynamic plastic response of beams and plates made up of rigidly perfect plastic materials subjected to dynamic loading. Many theoretical and experimental results can be reformulated into new concise forms with the response number. The concept of a new dimensionless number, response number, termed as Rn(n), is generalized in Zhao [Forschung im Ingenieurwesen 65 (1999) 107] to study the elastic, plastic, dynamic elastic as well as dynamic plastic buckling problems of columns, plates as well as shells. The response number Rn(n) is generalized to the dynamic behaviour of shells of various shapes in the present paper.

Dynamic properties and asymmetry of tension and compression of metal matrix composites under impact loading

The mechanical behaviors of 2124, Al-5Cu, Al-Li and 6061 alloys reinforced by silicon carbide particulates, together with 15%SiCw/6061 alloy, were studied under the quasi-static and impact loading conditions, using the split Hopkinson tension/compression bars and Instron universal testing machine. The effect of strain rate on the ultra tensile strength (UTS), the hardening modulus and the failure strain was investigated. At the same time, the SEM observations of dynamic fracture surfaces of various MMC materials showed some distinguished microstructures and patterns. Some new characteristics of asymmetry of mechanical behaviors of MMCs under tension and compression loading were also presented and explained in details, and they could be considered as marks to indicate, to some degree, the mechanism of controlling damage and failure of MMCs under impact loading. The development of new constitutive laws about MMCs under impact loading should benefit from these experimental results and theoretical analysis.