A peridynamic theory for linear elastic shells

A state-based peridynamic formulation for linear elastic shells is presented. The emphasis is on introducing, possibly for the first time, a general surface based peridynamic model to represent the deformation characteristics of structures that have one geometric dimension much smaller than the other two. A new notion of curved bonds is exploited to cater for force transfer between the peridynamic particles describing the shell. Starting with the three dimensional force and deformation states, appropriate surface based force, moment and several deformation states are arrived at. Upon application on the curved bonds, such states yield the necessary force and deformation vectors governing the motion of the shell. By incorporating a shear correction factor, the formulation also accommodates analysis of shells that have higher thickness. In order to attain this, a consistent second order approximation to the complementary energy density is considered and incorporated in peridynamics via constitutive correspondence. Unlike the uncoupled constitution for thin shells, a consequence of a first order approximation, constitutive relations for thick shells are fully coupled in that surface wryness influences the in-plane stress resultants and surface strain the moments. Our proposal on the peridynamic shell theory is numerically assessed against simulations on static deformation of spherical and cylindrical shells, that of flat plates and quasi-static fracture propagation in a cylindrical shell. (C) 2016 Elsevier Ltd. All rights reserved.

Cold gas in cluster cores: global stability analysis and non-linear simulations of thermal instability

We perform global linear stability analysis and idealized numerical simulations in global thermal balance to understand the condensation of cold gas from hot/virial atmospheres (coronae), in particular the intracluster medium (ICM). We pay particular attention to geometry (e.g. spherical versus plane-parallel) and the nature of the gravitational potential. Global linear analysis gives a similar value for the fastest growing thermal instability modes in spherical and Cartesian geometries. Simulations and observations suggest that cooling in haloes critically depends on the ratio of the cooling time to the free-fall time (t(cool)/t(ff)). Extended cold gas condenses out of the ICM only if this ratio is smaller than a threshold value close to 10. Previous works highlighted the difference between the nature of cold gas condensation in spherical and plane-parallel atmospheres; namely, cold gas condensation appeared easier in spherical atmospheres. This apparent difference due to geometry arises because the previous plane-parallel simulations focused on in situ condensation of multiphase gas but spherical simulations studied condensation anywhere in the box. Unlike previous claims, our non-linear simulations show that there are only minor differences in cold gas condensation, either in situ or anywhere, for different geometries. The amount of cold gas depends on the shape of tcool/tff; gas has more time to condense if gravitational acceleration decreases towards the centre. In our idealized plane-parallel simulations with heating balancing cooling in each layer, there can be significant mass/energy/momentum transfer across layers that can trigger condensation and drive tcool/tff far beyond the critical value close to 10.

FAST AND ACCURATE BILATERAL FILTERING USING GAUSS-POLYNOMIAL DECOMPOSITION

The bilateral filter is a versatile non-linear filter that has found diverse applications in image processing, computer vision, computer graphics, and computational photography. A common form of the filter is the Gaussian bilateral filter in which both the spatial and range kernels are Gaussian. A direct implementation of this filter requires O(sigma(2)) operations per pixel, where sigma is the standard deviation of the spatial Gaussian. In this paper, we propose an accurate approximation algorithm that can cut down the computational complexity to O(1) per pixel for any arbitrary sigma (constant-time implementation). This is based on the observation that the range kernel operates via the translations of a fixed Gaussian over the range space, and that these translated Gaussians can be accurately approximated using the so-called Gauss-polynomials. The overall algorithm emerging from this approximation involves a series of spatial Gaussian filtering, which can be efficiently implemented (in parallel) using separability and recursion. We present some preliminary results to demonstrate that the proposed algorithm compares favorably with some of the existing fast algorithms in terms of speed and accuracy.

General Relationship Between Contact Stiffness, Contact Depth, and Mechanical Properties for Indentation in Linear Viscoelastic Solids Using Axisymmetric Indenters of Arbitrary Profiles

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. Although the same expression is well known for indentation in elastic and in elastic-plastic solids, we show that it is also true for indentation in linear viscoelastic solids, provided that the unloading rate is sufficiently fast. Furthermore, the same expression holds true for both fast loading and unloading. These results should provide a sound basis for using the relationship for determining properties of viscoelastic solids using indentation techniques.

Identification of Nonlinear Systems Through Time-Frequency Filtering Technique

In the previous paper, a class of nonlinear system is mapped to a so-called skeleton linear model (SLM) based on the joint time-frequency analysis method. Behavior of the nonlinear system may be indicated quantitatively by the variance of the coefficients of SLM versus its response. Using this model we propose an identification method for nonlinear systems based on nonstationary vibration data in this paper. The key technique in the identification procedure is a time-frequency filtering method by which solution of the SLM is extracted from the response data of the corresponding nonlinear system. Two time-frequency filtering methods are discussed here. One is based on the quadratic time-frequency distribution and its inverse transform, the other is based on the quadratic time-frequency distribution and the wavelet transform. Both numerical examples and an experimental application are given to illustrate the validity of the technique.

Relationship Between Contact Stiffness, Contact Depth, and Mechanical Properties for Indentation in Linear Viscoelastic Solids Using Axisymmetric Indenters

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. Although the same expres

Non-Linear Wave-Induced Transient Response of Soil Around a Trenched Pipeline

Submarine pipelines are always trenched within a seabed for reducing wave loads and thereby enhancing their stability. Based on Biot’s poroelastic theory, a two-dimensional finite element model is developed to investigate non-linear wave-induced responses of soil around a trenched pipeline, which is verified with the flume test results by Sudhan et al. [Sudhan, C.M., Sundar, V., Rao, S.N., 2002. Wave induced forces around buried pipeline. Ocean Engineering, 29, 533–544] and Turcotte et al. [Turcotte, B.R., Liu, P.L.F., Kulhawy, F.H., 1984. Laboratory evaluation of wave tank parameters for wave-sediment interaction. Joseph H. Defree Hydraulic Laboratory Report 84-1, School of Civil and Environmental Engineering, Cornell University]. Non-linear wave-induced transient pore pressure around pipeline at various phases of wave loading is examined firstly. Unlike most previous investigations, in which only a single sediment layer and linear wave loading were concerned, in this study, the influences of the non-linearity of wave loading, the physical properties of backfill materials and the geometry profile of trenches on the excess pore pressures within the soil around pipeline, respectively, were explored, taking into account the in situ conditions of buried pipeline in the shallow ocean zones. Based on the parametric study, it is concluded that the shear modulus and permeability of backfill soils significantly affect the wave-induced excess pore pressures around trenched pipeline, and that the effect of wave non-linearity becomes more pronounced and comparable with that of trench depth, especially at high wave steepness in shallow water.

Nonparametric Bayesian Density Modeling with Gaussian Processes

We present the Gaussian process density sampler (GPDS), an exchangeable generative model for use in nonparametric Bayesian density estimation. Samples drawn from the GPDS are consistent with exact, independent samples from a distribution defined by a density that is a transformation of a function drawn from a Gaussian process prior. Our formulation allows us to infer an unknown density from data using Markov chain Monte Carlo, which gives samples from the posterior distribution over density functions and from the predictive distribution on data space. We describe two such MCMC methods. Both methods also allow inference of the hyperparameters of the Gaussian process.