scholarly journals Lamb's problem on random mass density fields with fractal and Hurst effects

Author(s):  
V. V. Nishawala ◽  
M. Ostoja-Starzewski ◽  
M. J. Leamy ◽  
E. Porcu

This paper reports on a generalization of Lamb's problem to a linear elastic, infinite half-space with random fields (RFs) of mass density, subject to a normal line load. Both, uncorrelated and correlated (with fractal and Hurst characteristics) RFs without any weak noise restrictions, are proposed. Cellular automata (CA) is used to simulate the wave propagation. CA is a local computational method which, for rectangular discretization of spatial domain, is equivalent to applying the finite difference method to the governing equations of classical elasticity. We first evaluate the response of CA to an uncorrelated mass density field, more commonly known as white-noise, of varying coarseness as compared to CA's node density. We then evaluate the response of CA to multiscale mass density RFs of Cauchy and Dagum type; these fields are unique in that they are able to model and decouple the field's fractal dimension and Hurst parameter. We focus on stochastic imperfection sensitivity; we determine to what extent the fractal or the Hurst parameter is a significant factor in altering the solution to the planar stochastic Lamb's problem by evaluating the coefficient of variation of the response when compared with the coefficient of variation of the RF.

Author(s):  
Xian Zhang ◽  
Martin Ostoja-Starzewski

This paper reports the application of cellular automata to study the dynamic responses of Lamb-type problems for a tangential point load and a concentrated moment applied on the free surface of a half-plane. The medium is homogeneous, isotropic and linear elastic while having a random mass density field with fractal and Hurst characteristics. Both Cauchy and Dagum random field models are used to capture these effects. First, the cellular automata approach is tested on progressively finer meshes to verify the code against the continuum elastodynamic solution in a homogeneous continuum. Then, the sensitivity of wave propagation on random fields is assessed for a wide range of fractal and Hurst parameters. Overall, the mean response amplitude is lowered by the mass density field’s randomness, while the Hurst parameter (especially, for β  < 0.2) is found to have a stronger influence than the fractal dimension on the response. The resulting Rayleigh wave is modified more than the pressure wave for the same random field parameters. Additionally, comparisons with previously studied Lamb-type problems under normal in-plane and anti-plane loadings are given. This article is part of the theme issue ‘Advanced materials modelling via fractional calculus: challenges and perspectives’.


2019 ◽  
Vol 8 (1) ◽  
pp. 231-246 ◽  
Author(s):  
Xian Zhang ◽  
◽  
Vinesh Nishawala ◽  
Martin Ostoja-Starzewski ◽  
◽  
...  

2020 ◽  
Vol 87 (4) ◽  
Author(s):  
K. S. Parel

Abstract A fundamental plane frictional receding contact of a thin layer pressed onto a substrate is solved with the insertion of distributions of dislocations. The new formulation overcomes the restriction of standard contact analysis with the distributed dislocation technique, which requires the distributions to be square root bounded to zero or square root singular at the ends of the area of insertion. This new formulation opens up a range of plane frictional receding contacts to theoretical study. Solutions are obtained for the basic plane receding contact problem of a homogeneous linear elastic semi-infinite layer pressed by a line load onto a half-plane of the same material with a Coulomb friction interface. The frictional dependence of receding contact behavior is identified and discussed.


Geophysics ◽  
1967 ◽  
Vol 32 (1) ◽  
pp. 88-98 ◽  
Author(s):  
Prabakar S. Naidu

Following a study of the statistical properties of a random potential field, which is often encountered as a background noise in observed geophysical potential fields, and particularly a study of conditions under which the random process is homogeneous and gaussian, an approximate analytical method has been developed to evaluate the correlation functions and the power spectrum of the random process when the mass density is locally uniform but varies randomly over an entire semi‐infinite space. The method is particularly convenient for it gives a unified approach to evaluating the statistical properties of gravity as well as magnetic fields. The relationship between the autocorrelation function of the random potential field and the random mass density is quite involved, but by assuming that the mass density in successive cells is completely uncorrelated the relationship can be simplified. It is a matter of considerable interest to estimate, from the autocorrelation function, the size of cell characteristic of the random medium which reflects the geological characteristics of the rocks.


2004 ◽  
Vol 14 (2) ◽  
pp. 109-117 ◽  
Author(s):  
K Murawski ◽  
L Nocera ◽  
E N Pelinovsky

2001 ◽  
Vol 09 (04) ◽  
pp. 1329-1345
Author(s):  
FRANZ ZIEGLER

Monitoring of structures which may be subjected to overloads can be based on observing signals emitted in the course of developing defects. Ductile structures react on overloads by the formation of (small scale) plastic zones. Within the linear elastic background concept such an event is considered by the formation of sources of eigenstress. Since a direct description of the source characteristics is rather cumbersome we choose the time convolution of the imposed plastic strains with the complementary Green's stress dyadic to describe the signal of the acoustice mission. Thus, the dynamic generalization of Maysel's formula of thermo-elasticity, to include all kinds of eigenstrains, enters the field of computational acoustics. In that context, the novel contribution of this paper to acoustic emission and monitoring of (layered) structures is the formulation of the full 3-D problems and the introduction of the generalized rays in the background considering an instantaneous oblique force point source at the transducer's location. That means, all the information on the wave guide is contained in the Green's stress dyadic. The expansion into plane waves of cylindrical or spherical waves propagating in a layered elastic half-space or plate proves to be quite efficient for short observation times. Even the divergence effects of dipping interfaces of wedge-type layers are perfectly included by proper coordinate rotations and the exact "seismograms" are observed at a point receiver (where the localized plastic source is assumed to develop, commonly buried and often localized at an interface) from any source located at the hypo center (the site of the transducer in receiving mode, commonly placed at the "outer" surface). This nontrivial technique relies on the invariance of the phase function (arrival time) and of the infinitesimal amplitude of the plane waves in the ray expansion. The concept of the elastic background is illustrated by elastic-viscoplastic waves propagating in thin rods and subsequently extended to the 3-D problem of spherical waves with point symmetry. In that context and in an incremental formulation, the notion of plastic sources is introduced, which emit elastic waves in the background. Finally, the full 3-D problem in a layered half space or layered plate is solved in terms of generalized rays to be received at a transducer in receiving mode. Taking into account the progress in symbolic manipulation with integrated numeric capabilities (e.g., of Mathematica), such a formulation seems timely and may prove to be competitive to the entirely computational Finite Element Method of analysis of signals received from plastic sources. Time signatures of Green's displacement components at the surface of a half-space are illustrated when produced by a vertical, horizontal and inclined line load with a triangular time source function, respectively.


2016 ◽  
Vol 33 (1) ◽  
pp. 28-63 ◽  
Author(s):  
Richard Regueiro ◽  
Zheng Duan ◽  
Beichuan Yan

Purpose – The purpose of this paper is to develop a concurrent multiscale computational method for granular materials in the quasi-static loading regime. Design/methodology/approach – Overlapped-coupling between a micropolar linear elastic one-dimensional (1D) mixed finite element (FE) model and a 1D chain of Hertzian nonlinear elastic, glued, discrete element (DE) spheres is presented. The 1D micropolar FEs and 1D chain of DEs are coupled using a bridging-scale decomposition for static analysis. Findings – It was found that an open-window DE domain may be coupled to a micropolar continuum FE domain via an overlapping region within the bridging-scale decomposition formulation for statics. Allowing the micropolar continuum FE energy in the overlapping region to contribute to the DE energy has a smoothing effect on the DE response, especially for the rotational degrees of freedom (dofs). Research limitations/implications – The paper focusses on 1D examples, with elastic, glued, DE spheres, and a linear elastic micropolar continuum implemented in 1D. Practical implications – A concurrent computational multiscale method for granular materials with open-window DE resolution of the large shearing region such as at the interface with a penetrometer skin, will allow more efficient computations by reducing the more costly DE domain calculations, but not at the expense of generating artificial boundary effects between the DE and FE domains. Originality/value – Open-window DE overlapped-coupling to FE continuum domain, accounting for rotational dofs in both DE and FE methods. Contribution of energy from micropolar FE in overlap region to underlying DE particle energy.


2013 ◽  
Vol 845 ◽  
pp. 324-329 ◽  
Author(s):  
Theng Pin Ng ◽  
Seyed Saeid Rahimian Koloor ◽  
Joy Rizki Pangestu Djuansjah ◽  
Mohammed Rafiq Abdul Kadir

Linear elastic response of the bovine cortical bone has been examined under compression load. Experimental and computational methods were used to observe and predict the response of cortical bone. In computational method, two mechanical behaviors of isotropic and orthotropic were considered to simulate the cortical bone deformation. In experimental process, the specimens were designed to show maximum stiffness and strength by specifying osteon direction along loading axis during tests. The tests were controlled by displacement rate of 0.5 mm/minute and the overall stiffness responses of the structures were recorded to extract mechanical properties and also for validation aims. Finite Element Method (FEM) was used to model the linear response of the structure by using ABAQUS6.9EF. The FE results using orthotropic definition shows a good correlation with experimental data. A discussion was given based on overall stiffness and effective stress variation for both mechanical behaviors. In order to design the optimal implant structure, the presented study was proposed for prediction of bone structure deformation that attached to the orthopedic implants.


Sign in / Sign up

Export Citation Format

Share Document