Material model of lung parenchyma based on living precision-cut lung slice testing.

We describe a novel constitutive model of lung parenchyma, which can be used for continuum mechanics based predictive simulations. To develop this model, we experimentally determined the nonlinear material behavior of rat lung parenchyma. This was achieved via uni-axial tension tests on living precision-cut rat lung slices. The resulting force-displacement curves were then used as inputs for an inverse analysis. The Levenberg-Marquardt algorithm was utilized to optimize the material parameters of combinations and recombinations of established strain-energy density functions (SEFs). Comparing the best-fits of the tested SEFs we found Wpar = 4.1 kPa(I1-3)2 + 20.7 kPa(I1 - 3)3 + 4.1 kPa(-2 ln J + J2 - 1) to be the optimal constitutive model. This SEF consists of three summands: the first can be interpreted as the contribution of the elastin fibers and the ground substance, the second as the contribution of the collagen fibers while the third controls the volumetric change. The presented approach will help to model the behavior of the pulmonary parenchyma and to quantify the strains and stresses during ventilation.

A zero-stiffness elastic shell structure

A remarkable shell structure is described that, due to a particular combination of geometry and initial stress, has zero stiffness for any finite deformation along a twisting path; the shell is in a neutrally stable state of equilibrium. Initially the shell is straight in a longitudinal direction, but has a constant, nonzero curvature in the transverse direction. If residual stresses are induced in the shell by, for example, plastic deformation, to leave a particular resultant bending moment, then an analytical inextensional model of the shell shows it to have no change in energy along a path of twisted configurations. Real shells become closer to the inextensional idealization as their thickness is decreased; experimental thin-shell models have confirmed the neutrally stable configurations predicted by the inextensional theory. A simple model is described that shows that the resultant bending moment that leads to zero stiffness gives the shell a hidden symmetry, which explains this remarkable property.

Mechanisms of erosion of unfilled elastomers by solid particle impact

The mechanisms of material removal were studied during the erosion of two unfilled elastomers (natural rubber and epoxidised natural rubber). The effects of impact velocity and of lubrication by silicone oil were investigated. The development of surface features due to single impacts and during the early stages of erosion was followed by scanning electron microscopy. The basic material removal mechanism at impact angles of both 30° and 90° involves the formation and growth of fine fatigue cracks under the tensile surface stresses caused by impact. No damage was observed after single impacts; it was found that many successive impacts are necessary for material removal. It was found that the erosion rate has a very strong dependance on impact velocity above about 50 ms-1.

Solid particle erosion of boronized steels

Boronizing is a thermochemical diffusion-based process for producing iron boride layers in the surface of steel components. The boride layer is wear resistant and is very hard. Large residual stresses are found to exist in the surface layers, which are a function of substrate steel composition and heat treatment. By slow cooling from the boronizing temperature (900°C), a large compressive stress is developed in the boride layer. Hardening the steel by rapid cooling, either directly from the boronizing treatment or after subsequent austenitizing, develops tension in the coating which causes it to fracture. Tempering of the martensite produces compression in the coating, closing but not welding the cracks. The results of solid particle erosion experiments using silicon carbide, quartz, and glass bead erodents on boronized steels are presented.

Visualization of turbulent wedges under favourable pressure gradients using both shear-sensitive and temperature-sensitive liquid crystals

Turbulent wedges induced by a 3D surface roughness placed in a laminar boundary layer over a flat plate were visualised for the first time using both shear-sensitive and temperature-sensitive liquid crystals. The experiments were carried out at three different levels of favourable pressure gradients. The purpose of this investigation was to examine the spreading angles of the turbulent wedges indicated by their associated surface shear stresses and heat transfer characteristics and hence obtain further insight about the difference in the behaviour of transitional momentum and thermal boundary layers when a streamwise pressure gradient exists. It was shown that under a zero pressure gradient the spreading angles indicated by the two types of liquid crystals are the same, but the difference increases as the level of favourable pressure gradient increases. The result from the present study could have an important implication to the transition modelling of thermal boundary layers over gas turbine blades.

The unsteady development of a turbulent wake through a downstream low-pressure turbine blade passage

This paper presents two-dimensional LDA measurements of the convection of a wake through a low-pressure (LP) turbine cascade. Previous studies have shown the wake convection to be kinematic but have not provided details of the turbulent field. The spatial resolution of these measurements has facilitated the calculation of the production of turbulent kinetic energy and this has revealed a mechanism for turbulence production as the wake converts through the bladerow. The measured ensemble-averaged velocity field confirmed the previously reported kinematics of wake convection while the measurements of the turbulence quantities showed the wake fluid to be characterised by elevated levels of turbulent kinetic energy (TKE) and to have an anisotropic structure. Based on the measured mean and turbulence quantities, the production of turbulent kinetic energy was calculated. This highlighted a TKE production mechanism that resulted in increased levels of turbulence over the rear suction surface where boundary layer transition occurs. The turbulence production mechanism within the bladerow was also observed to produce more nearly isotropic turbulence. Production occurs when the principal stresses within the wake are aligned with the mean strains. This coincides with the maximum distortion of the wake within the blade passage and provides a mechanism for the production of turbulence outside of the boundary layer.

LES and hybrid LES/RANS simulations for conjugate heat transfer over a matrix of cubes

Large Eddy Simulation (LES) and a novel k -l based hybrid LES/RANS approach have been applied to simulate a conjugate heat transfer problem involving flow over a matrix of surface mounted cubes. In order to assess the capability and reliability of the newly developed k -l based hybrid LES/RANS, numerical results are compared with new LES and existing RANS results. Comparisons include mean velocity profiles, Reynolds stresses and conjugate heat transfer. As well as for hybrid LES/RANS validation purposes, the LES results are used to gain insights into the complex flow physics and heat transfer mechanisms. Numerical simulations show that the hybrid LES/RANS approach is effective. Mean and instantaneous fluid temperatures adjacent to the cube surface are found to strongly correlate with flow structure. Although the LES captures more mean velocity field complexities, broadly time averaged wake temperature fields are found similar for the LES and hybrid LES/RANS. Copyright © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

A model for nonlinear viscoelastic mechanical responses of collagenous soft tissues

Experimental observations of the time-dependent mechanical responses of collagenous tissues have demonstrated behavior that deviates from standard treatments of linear or quasi-linear viscoelasticity. In particular, time-dependent deformation can be strongly coupled to strain level, and strain-rate independence can be observed under monotonic loading, even for a tissue with dramatic stress relaxation. It was postulated that this nonlinearity is fundamentally associated with gradual recruitment of individual collagen fibrils during applied mechanical loading. Based on previously observed experimental results for the time-dependent response of collagenous soft tissues, a model is developed to describe the mechanical behavior of these tissues under uniaxial loading. Tissue stresses, under applied strain-controlled loading, are assumed to be a sum of elastic and viscoelastic stress contributions. The relative contributions of elastic and viscoelastic stresses is assumed to vary with strain level, leading to strain- and time-dependent mechanical behavior. The model formulation is examined under conditions of monotonic loading at varying constant strain rates and stress-relaxation at different applied strain levels. The model is compared with experimental data for a membranous biological soft tissue, the amniotic sac, and is found to agree well with experimental results. The limiting behavior of the novel model, at large strains relative to the collagen recruitment, is consistent with the quasi-linear viscoelastic approach. © 2006 Materials Research Society.

Analysis of case histories on deep excavations in marine clay

Two case histories on deep excavation of marine clay are used to study the use of a decision-making tool based on a new deign method called the Mobilized Strength Design (MSD) method which allows the designer to use a simple method of predicting ground displacements during deep excavation. This application can approximately satisfy both safety and serviceability requirements by predicting stresses and displacements under working conditions by introducing the concept of "Mobilizable soil strength". The new method accommodates a number of features which are important to design of underground construction between retaining walls, including different deformation mechanism in different stages of excavation. The influence of wall depth, wall flexibility and stratified ground are the major focus of this paper. These developments should make it possible for a design engineer to take informed decisions on the influence of wall stiffness, or on the need for a jet-grouted base slab, for example, without having to conduct project-specific Finite Element Analysis.

Assessment of CFD models for shock boundary-layer interaction

A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and in general it was difficult to discern clear trends in the data. For the Reynolds Averaged Navier-Stokes methods the choice of turbulence model appeared to be the largest factor in solution accuracy. Large-eddy simulation methods produced error levels similar to RANS methods but provided superior predictions of normal stresses.

Centrifuge modelling of pile-soil interaction in liquefiable slopes

Piles passing through sloping liquefiable deposits are prone to lateral loading if these deposits liquefy and flow during earthquakes. These lateral loads caused by the relative soil-pile movement will induce bending in the piles and may result in failure of the piles or excessive pile-head displacement. Whilst the weak nature of the flowing liquefied soil would suggest that only small loads would be exerted on the piles, it is known from case histories that piles do fail owing to the influence of laterally spreading soils. It will be shown, based on dynamic centrifuge test data, that dilatant behaviour of soil close to the pile is the major cause of these considerable transient lateral loads which are transferred to the pile. This paper reports the results of geotechnical centrifuge tests in which models of gently sloping liquefiable sand with pile foundations passing through them were subjected to earthquake excitation. The soil close to the pile was instrumented with pore-pressure transducers and contact stress cells in order to monitor the interaction between soil and pile and to track the soil stress state both upslope and downslope of the pile. The presence of instrumentation measuring pore-pressure and lateral stress close to the pile in the research described in this paper gives the opportunity to better study the soil stress state close to the pile and to compare the loads measured as being applied to the piles by the laterally spreading soils with those suggested by the JRA design code. This test data shows that lateral stresses much greater than one might expect from calculations based on the residual strength of liquefied soil may be applied to piles in flowing liquefied slopes owing to the dilative behaviour of the liquefied soil. It is shown at least for the particular geometry studied that the current JRA design code can be un-conservative by a factor of three for these dilation-affected transient lateral loads.

Cell adhesion to plasma electrolytic oxidation (PEO) titania coatings, assessed using a centrifuging technique.

The adhesion of bovine chondrocytes and human osteoblasts to three titania-based coatings, formed by plasma electrolytic oxidation (PEO), was compared to that on uncoated Ti-6Al-4V substrates, and some comparisons were also made with plasma sprayed hydroxyapatite (HA) coatings. This was done using a centrifuge, with accelerations of up to 160,000 g, so as to induce buoyancy forces that created normal or shear stresses at the interface. It is shown that, on all surfaces, it was easier to remove cells under normal loading than under shear loading. Cell adhesion to the PEO coatings was stronger than that on Ti-6Al-4V and similar to that on HA. Cell proliferation rates were relatively high on one of the PEO coatings, which was virtually free of aluminium, but low on the other two, which contained significant levels of aluminium. It is concluded that the Al-free PEO coating offers promise for application to prosthetic implants.

The effects of a second tunnel on the propagation of ground-borne vibration from an underground railway

Accurate predictions of ground-borne vibration levels in the vicinity of an underground railway are greatly sought in modern urban centers. Yet the complexity involved in simulating the underground environment means that it is necessary to make simplifying assumptions about this environment. One such commonly-made assumption is to model the railway as a single tunnel, despite many underground railway lines consisting of twin-bored tunnels. A unique model for two tunnels embedded in a homogeneous, elastic full space is developed. The vibration response of this two-tunnel system is calculated using the superposition of two displacement fields: one resulting from the forces acting on the invert of a single tunnel, and the other resulting from the interaction between the tunnels. By partitioning of the stresses into symmetric and anti-symmetric mode number components using Fourier decomposition, these two displacement fields can by calculated with minimal computational requirements. The significance of the interactions between twin-tunnels is quantified by calculating the insertion gains that result from the existence of a second tunnel. The insertion-gain results are shown to be localized and highly dependent on frequency, tunnel orientation and tunnel thickness. At some locations, the magnitude of these insertion gains is greater than 20dB. This demonstrates that a high degree of inaccuracy exists in any surface vibration-prediction model that includes only one of the two tunnels. © 2012 Springer.

Moderate energy impact analysis combining phenomenological contact law with localised damage and integral equation method

A computational impact analysis methodology has been developed, based on modal analysis and a local contact force-deflection model. The contact law is based on Hertz contact theory while contact stresses are elastic, defines a modified contact theory to take account of local permanent indentation, and considers elastic recovery during unloading. The model was validated experimentally through impact testing of glass-carbon hybrid braided composite panels. Specimens were mounted in a support frame and the contact force was inferred from the deceleration of the impactor, measured by high-speed photography. A Finite Element analysis of the panel and support frame assembly was performed to compute the modal responses. The new contact model performed well in predicting the peak forces and impact durations for moderate energy impacts (15 J), where contact stresses locally exceed the linear elastic limit and damage may be deemed to have occurred. C-scan measurements revealed substantial damage for impact energies in the range of 30-50 J. For this regime the new model predictions might be improved by characterisation of the contact law hysteresis during the unloading phase, and a modification of the elastic vibration response in line with damage levels acquired during the impact. © 2011 Elsevier Ltd. All rights reserved.

Design methodology for composite structures: A small low air-speed wind turbine blade case study

A small low air-speed wind turbine blade case study is used to demonstrate the effectiveness of a materials and design selection methodology described by Monroy Aceves et al. (2008) [24] for composite structures. The blade structure comprises a shell of uniform thickness and a unidirectional reinforcement. The shell outer geometry is fixed by aerodynamic considerations. A wide range of lay-ups are considered for the shell and reinforcement. Structural analysis is undertaken using the finite element method. Results are incorporated into a database for analysis using material selection software. A graphical selection stage is used to identify the lightest blade meeting appropriate design constraints. The proposed solution satisfies the design requirements and improves on the prototype benchmark by reducing the mass by almost 50%. The flexibility of the selection software in allowing identification of trends in the results and modifications to the selection criteria is demonstrated. Introducing a safety factor of two on the material failure stresses increases the mass by only 11%. The case study demonstrates that the proposed design methodology is useful in preliminary design where a very wide range of cases should be considered using relatively simple analysis. © 2011 Elsevier Ltd.