Efficient numerical models for investigating the statistics of random dynamic systems

The study of random dynamic systems usually requires the definition of an ensemble of structures and the solution of the eigenproblem for each member of the ensemble. If the process is carried out using a conventional numerical approach, the computational cost becomes prohibitive for complex systems. In this work, an alternative numerical method is proposed. The results for the response statistics are compared with values obtained from a detailed stochastic FE analysis of plates. The proposed method seems to capture the statistical behaviour of the response with a reduced computational cost.

Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. © 2011 John Wiley & Sons, Ltd.

The effect of inclined soil layers on surface vibration from underground railways using a semi-analytical approach

Ground vibration due to underground railways is a significant source of disturbance for people living or working near the subways. The numerical models used to predict vibration levels have inherent uncertainty which must be understood to give confidence in the predictions. A semi-analytical approach is developed herein to investigate the effect of soil layering on the surface vibration of a halfspace where both soil properties and layer inclination angles are varied. The study suggests that both material properties and inclination angle of the layers have significant effect ( ± 10dB) on the surface vibration response. © 2009 IOP Publishing Ltd.

Flexoelectro-optic properties of chiral nematic liquid crystals in the uniform standing helix configuration.

The flexoelectro-optic effect describes the rotation of the optic axis of a short-pitch chiral nematic liquid crystal under the application of an electric field. We investigate the effect in the uniform standing helix, or "Grandjean" configuration. An in-plane electric field is applied. The director profile is determined numerically using a static one-dimensional continuum model with strong surface anchoring. The Berreman method is used to solve for plane-wave solutions to Maxwell's equations, and predict the optical properties of the resulting structure in general cases. By using a chiral nematic with short pitch between crossed polarizers an optical switch may be generated. With no applied field the configuration is nontransmissive at normal incidence, but becomes transmissive with an applied field. For this case, numerical results using the Berreman method are supplemented with an analytic theory and found to be in good agreement. The transmitted intensity as a function of tilt, the contrast ratio, and the tilt required for full intensity modulation are presented. The angular dependence of the transmission is calculated and the isocontrast curves are plotted. For typical material and cell parameters a switching speed of 0.017 ms and contrast ratio of 1500:1 at normal incidence are predicted, at a switch-on tilt of 41.5 degrees. Experimental verification of the analytic and numerical models is provided.

An experimentally based "bubble' model for clay

In order to extend the Cam clay model an inner kinematically hardening bubble is added in order to match experimentally observed hysteretic response. The new model is found to be quite successful in reproducing patterns of response observed in slow cyclic oedometer and triaxial tests. -from Authors

Energy scattering in weakly non-linear systems.

The Chinese Tam-Tam exhibits non-linear behavior in its vibro-acoustic response. The frequency content of the response during free, unforced vibration smoothly changes, with energy being progressively smeared out over a greater bandwidth with time. This is used as a motivating case for the general study of the phenomenon of energy cascading through weak nonlinearity. Numerical models based upon the Fermi-Pasta-Ulam system of non-linearly coupled oscillators, modified with the addition of damping, have been developed. These were used to study the response of ensembles of systems with randomized natural frequencies. Results from simulations will be presented here. For un-damped systems, individual ensemble members exhibit cyclical energy exchange between linear modes, but the ensemble average displays a steady state. For the ensemble response of damped systems, lightly damped modes can exhibit an effective damping which is higher than predicated by linear theory. The presence of a non-linearity provides a path for energy flow to other modes, increasing the apparent damping spectrum at some frequencies and reducing it at others. The target of this work is a model revealing the governing parameters of a generic system of this type and leading to predictions of the ensemble response.

Modelling of liquefaction-induced instability in pile

Centrifuge testing has been undertaken to investigate instability failure of pile groups during seismic liquefaction, with specific reference to the 'top-down' propagation of liquefaction during the earthquake and to account for initial imperfections in pile geometry. The results of these tests were used to validate numerical models within the finite element program ABAQUS, based on the popular p-y analysis method. Pseudostatic classical and post-buckling analyses were conducted to examine the collapse behaviour of the pile groups and were found to give reasonable predictions of collapse load and conservative predictions of the associated deflection conditions. This numerical model was compared to currently published methods which were found to over-predict collapse loads. The resulting insights into the collapse of axially loaded pile groups revealed that the failure load is strongly dependent on both the depth of liquefaction propagation and initial imperfections, which reduce the collapse load.

Predicting surface vibration from underground railways through inhomogeneous soil

Noise and vibration from underground railways is a major source of disturbance to inhabitants near subways. To help designers meet noise and vibration limits, numerical models are used to understand vibration propagation from these underground railways. However, the models commonly assume the ground is homogeneous and neglect to include local variability in the soil properties. Such simplifying assumptions add a level of uncertainty to the predictions which is not well understood. The goal of the current paper is to quantify the effect of soil inhomogeneity on surface vibration. The thin-layer method (TLM) is suggested as an efficient and accurate means of simulating vibration from underground railways in arbitrarily layered half-spaces. Stochastic variability of the soils elastic modulus is introduced using a KL expansion; the modulus is assumed to have a log-normal distribution and a modified exponential covariance kernel. The effect of horizontal soil variability is investigated by comparing the stochastic results for soils varied only in the vertical direction to soils with 2D variability. Results suggest that local soil inhomogeneity can significantly affect surface velocity predictions; 90 percent confidence intervals showing 8 dB averages and peak values up to 12 dB are computed. This is a significant source of uncertainty and should be considered when using predictions from models assuming homogeneous soil properties. Furthermore, the effect of horizontal variability of the elastic modulus on the confidence interval appears to be negligible. This suggests that only vertical variation needs to be taken into account when modelling ground vibration from underground railways. © 2012 Elsevier Ltd. All rights reserved.

Inhomogeneous soils and their effect on ground vibration due to underground railways

Ground vibration due to underground railways is a significant source of disturbance for people living or working near subways. Numerical models are commonly used to predict vibration levels; however, uncertainty inherent to these simulations must be understood to give confidence in the predictions. A semi-analytical approach is developed herein to investigate the effect of uncertainty in soil material properties on the surface vibration of layered halfspaces excited by an underground railway. The half-space is simulated using the thin-layer method coupled with the pipe-in-pipe (PiP) method for determining the load on the buried tunnel. The K-L expansion method is employed to smoothly vary the material properties throughout the soil by up to 10%. The simulation predicts a surface rms velocity variation of 5-10dB compared to a homogeneous, layered halfspace. These results suggest it may be prudent to include a 5dB error band on predicted vibration levels when simulating areas of varied material properties.

Modification of the trapped field in bulk high-temperature superconductors as a result of the drilling of a pattern of artificial columnar holes

The trapped magnetic field is examined in bulk high-temperature superconductors that are artificially drilled along their c-axis. The influence of the hole pattern on the magnetization is studied and compared by means of numerical models and Hall probe mapping techniques. To this aim, we consider two bulk YBCO samples with a rectangular cross-section that are drilled each by six holes arranged either on a rectangular lattice (sample I) or on a centered rectangular lattice (sample II). For the numerical analysis, three different models are considered for calculating the trapped flux: (i), a two-dimensional (2D) Bean model neglecting demagnetizing effects and flux creep, (ii), a 2D finite-element model neglecting demagnetizing effects but incorporating magnetic relaxation in the form of an E-J power law, and, (iii), a 3D finite element analysis that takes into account both the finite height of the sample and flux creep effects. For the experimental analysis, the trapped magnetic flux density is measured above the sample surface by Hall probe mapping performed before and after the drilling process. The maximum trapped flux density in the drilled samples is found to be smaller than that in the plain samples. The smallest magnetization drop is found for sample II, with the centered rectangular lattice. This result is confirmed by the numerical models. In each sample, the relative drops that are calculated independently with the three different models are in good agreement. As observed experimentally, the magnetization drop calculated in the sample II is the smallest one and its relative value is comparable to the measured one. By contrast, the measured magnetization drop in sample (1) is much larger than that predicted by the simulations, most likely because of a change of the microstructure during the drilling process.