The effect of anisotropic elasticity on the yielding characteristics of overconsolidated natural clay

Quantitative application of elastoplastic theory to the yielding behaviour of natural soils has always been uncertain. Part of the reason is that the theory was developed for reconstituted materials with isotropic structure, in contrast to natural soils that are usually anisotropic. The approach considered in this study assumes that pre-yielding behaviour is governed by the theory of linear anisotropic elasticity and that yield loci in the mean effective stress ( p') – deviator stress (q) plane are aligned approximately along the coefficient of earth pressure (K0) line. The assumption of a rotated yield locus associated with anisotropic elastic behaviour within the state boundary surface indicates that the elastic wall within the state boundary surface is inclined. The form of the state boundary surface has been determined mathematically in terms of anisotropic elastic and Cam-Clay soil parameters. Stress path tests were conducted on samples of Belfast Upper Boulder Clay removed from a depth of 28 m below ground surface. Good agreement was found between predicted and measured yield loci. The study also examined the influence of subsequent isotropic compression on the yielding characteristics of the natural clay. The indications are that the anisotropy developed during deposition disappears when the sample is loaded to a stress level at least twice the stress generated during the original deposition process. The methods developed in the paper have also been applied to test results reported previously on Winnipeg clay, and good agreement was obtained.

Channelled Waves in Anisotropic Mesh Metamaterials

Channelled waves in 2-D periodic anisotropic L-C mesh metamaterials have been investigated. Circuit simulation and the newly developed analytical model of a unit cell have demonstrated full qualitative agreement for both lossless and lossy cases. Isofrequencies for a lattice unit cell and the circuit simulations of finite meshes have shown that propagating waves are channelled from a point source as pencil beams which can travel only along specific trajectories. The beam direction varies with frequency, and at the resonance frequency, the phase and group velocities of the travelling wave are orthogonal. The effect of losses was explored, and it was shown that losses cause qualitative changes of the channelled wave type. It was proven that the channelled waves do not follow the laws of geometrical optics (Snell's law, specular reflection, etc.) at the interfaces of L-C meshes but are governed by the conditions of phase synchronism and impedance matching. Only in the special case of dual L-C and C-L meshes with the interface parallel to the axis of rectangular grid excited at the resonance frequency (X=1) do the channels follow the trajectories of optical rays. A planar mesh test cell has been designed and used for retrieving the unit cell L-C parameters from the S-parameter measurements.

Anisotropic optical properties of arrays of gold nanorods embedded in alumina

A series of thin films comprising gold nanorods embedded in an alumina matrix have been fabricated with lengths ranging from 75 to 330 nm. Their optical properties, expressed in terms of extinction - In(T), where T is optical transmittance, have been measured as a function of wavelength, rod length, angle of incidence, and incident polarization state. The results are compared to a Maxwell-Garnett based theory modified to take into account the strongly anisotropic nature of the medium. Transverse and longitudinal plasmon resonances are observed. The interaction between the nanorods leads to the splitting of the longitudinal resonance with the longer-wavelength resonance being forbidden for direct optical observations. The shorter-wavelength resonance related to the symmetric coupling between longitudinal plasma excitations in the nanorods depends on rod length, polarization state, and angle of incidence of the probing light. The impact of electron confinement on the optical properties of the gold rods is also seen and may be incorporated into the Maxwell-Garnett theory by restricting the mean free path of the conduction electrons to produce excellent agreement between observations and the complete theory. Annealing experiments that modify the physical structure of the gold confirm this conclusion.

WAVES IN 2D ANISOTROPIC L-C LATTICE METAMATERIALS:PHENOMENOLOGY AND PROPERTIES

Anisotropic metamaterials composed of 2D periodic infi- nite and finite periodic lattices of lumped inductor (L) and capacitor (C) circuits have been explored. The unique features of wave channeling on such anisotropic lattices and scattering at their interfaces and edges are reviewed and illustrated by the examples of the specific arrangements. The lattice unit cells composed of inductors and capacitors (basic mesh) as well as of assemblies comprised of double series, double parallel, and mixed parallel-series L-C circuits are discussed.

Mechanical behaviour of unsaturated kaolin (with isotropic and anisotropic stress history). Part 1: wetting and compression behaviour

Over the last 40 years considerable progress has been made in understanding the complex behaviour of unsaturated soils. Research using constitutive modelling has extended the critical state framework and the concept of yielding in saturated soils to encompass unsaturated soils experiencing suction. However, validation testing of the framework for unsaturated soils has shown disagreement with the basic propositions. The main reason for this disparity is the anisotropic properties of the soil specimens tested as a result of preparation using one-dimensional compaction. The paper describes the detailed testing carried out to justify this statement. As part of the work presented, samples of unsaturated kaolin were prepared using isotropic compression. The suctions in these samples were reduced to predefined values by wetting under low isotropic loading. The pore size distributions, the pressure–volume relationships and yielding under subsequent isotropic loading are compared with tests on samples prepared by statically compressing kaolin into a one-dimensional compaction mould. The anisotropically compressed samples had initial water contents and specific volumes similar to those of the isotropically prepared samples and were also tested under reducing suctions; they exhibited distinctly different behaviour when tested under similar conditions. The results obtained from the isotropically prepared and tested samples have shown, probably for the first time, the existence of a unique normal compression surface that is not dependent on the initial conditions of the samples. The shape of the loading–collapse (LC) yield locus is shown to be different from the generally accepted form.

Mechanical behaviour of unsaturated kaolin (with isotropic and anisotropic stress history). Part 2: performance under shear loading

Validation of a framework for unsaturated soil behaviour has frequently resulted in disagreement with basic propositions. A primary reason for this disparity is considered to be attributable to the anisotropic properties of the soil specimens tested as a result of preparation using one-dimensional compaction. As part of the work presented, comparison is made between tests on samples of unsaturated kaolin prepared at identical specific volumes and specific water volumes using isotropic compression and one-dimensional compression. The suctions in the samples were reduced to predefined values by wetting under low isotropic loading in a triaxial cell. The samples were then taken through various stress paths to failure, defined as the critical state strength, while the suctions were held constant. Stress path tests were also performed on samples without reducing the suction to predefined values. In the latter, constant water mass tests, the suctions were allowed to vary and were measured using a psychrometer. The results of the tests at critical state are compared with the propositions of Wheeler and Sivakumar. The shear strengths of samples with isotropic previous history are shown to be significantly greater than those of samples with one-dimensional stress history when plotted against the mean net stress. The normal compression lines, critical state lines and yield characteristics are also shown to be significantly influenced by the previous stress history and are shown to be different for isotropically and one-dimensionally prepared samples.

Estimation of in situ stresses using anisotropic elasticity and suction measurements

Sampling and specimen preparation produce changes in mean effective stresses and pore water pressures, even with ‘perfect sampling’. The paper takes an existing simplified three-parameter cross-anisotropic elastic model and uses it to model these changes. The required ratio of cross-anisotropic parameters J/3G* can be obtained from standard CIU triaxial tests. If measurements are also made of suctions in unloaded specimens in the laboratory, then a combination of J/3G*, the measured suction, and the effective overburden pressure permits an estimation of the horizontal effective pressure and the K 0 ‘at rest’ coefficient. This can be helpful in numerical modelling that needs to start from in situ conditions, and in planning pressure levels for reconsolidation of clay specimens in the laboratory. Tests were done on Belfast Upper Boulder Clay from a depth of 28 m. Values of horizontal in situ effective stress estimated from these measurements compare favorably with conventional estimates of the ‘at rest’ coefficient K 0 and the overconsolidation ratio. Estimates of horizontal stress in London Clay were made using published data and the results compared with actual measurements. Again reasonable agreement was obtained.

Anisotropic seismic inversion using a multigrid Monte Carlo approach

We propose a new approach for the inversion of anisotropic P-wave data based on Monte Carlo methods combined with a multigrid approach. Simulated annealing facilitates objective minimization of the functional characterizing the misfit between observed and predicted traveltimes, as controlled by the Thomsen anisotropy parameters (epsilon, delta). Cycling between finer and coarser grids enhances the computational efficiency of the inversion process, thus accelerating the convergence of the solution while acting as a regularization technique of the inverse problem. Multigrid perturbation samples the probability density function without the requirements for the user to adjust tuning parameters. This increases the probability that the preferred global, rather than a poor local, minimum is attained. Undertaking multigrid refinement and Monte Carlo search in parallel produces more robust convergence than does the initially more intuitive approach of completing them sequentially. We demonstrate the usefulness of the new multigrid Monte Carlo (MGMC) scheme by applying it to (a) synthetic, noise-contaminated data reflecting an isotropic subsurface of constant slowness, horizontally layered geologic media and discrete subsurface anomalies; and (b) a crosshole seismic data set acquired by previous authors at the Reskajeage test site in Cornwall, UK. Inverted distributions of slowness (s) and the Thomson anisotropy parameters (epsilon, delta) compare favourably with those obtained previously using a popular matrix-based method. Reconstruction of the Thomsen epsilon parameter is particularly robust compared to that of slowness and the Thomsen delta parameter, even in the face of complex subsurface anomalies. The Thomsen epsilon and delta parameters have enhanced sensitivities to bulk-fabric and fracture-based anisotropies in the TI medium at Reskajeage. Because reconstruction of slowness (s) is intimately linked to that epsilon and delta in the MGMC scheme, inverted images of phase velocity reflect the integrated effects of these two modes of anisotropy. The new MGMC technique thus promises to facilitate rapid inversion of crosshole P-wave data for seismic slownesses and the Thomsen anisotropy parameters, with minimal user input in the inversion process.

Structural change of vortex patterns in anisotropic Bose-Einstein condensates

We study the changes in the spatial distribution of vortices in a rotating Bose-Einstein condensate due to an increasing eccentricity of the trapping potential. By breaking the rotational symmetry, the vortex system undergoes a rich variety of structural changes, including the formation of zigzag and linear configurations. These spatial rearrangements are well signaled by the change in the behavior of the vortex-pattern eigenmodes against the eccentricity parameter. This behavior allows to actively control the distribution of vorticity in many-body systems and opens the possibility of studying interactions between quantum vortices over a large range of parameters.