Corneal biomechanics: a review

Biomechanics is often defined as ‘mechanics applied to biology’. Due to the variety and complexity of the behaviour of biological structures and materials, biomechanics is better defined as the development, extension and application of mechanics for a better understanding of physiology and physiopathology and consequently for a better diagnosis and treatment of disease and injury. Different methods for the characterisation of corneal biomechanics are reviewed in detail, including those that are currently commercially available (Ocular Response Analyzer and CorVis ST). The clinical applicability of the parameters provided by these devices are discussed, especially in the fields of glaucoma, detection of ectatic disorders and orthokeratology. Likewise, other methods are also reviewed, such as Brillouin microscopy or dynamic optical coherence tomography and others with potential application to clinical practice but not validated for in vivo measurements, such as ultrasonic elastography. Advantages and disadvantages of all these techniques are described. Finally, the concept of biomechanical modelling is revised as well as the requirements for developing biomechanical models, with special emphasis on finite element modelling.

A simple method to predict elastic settlements in foundations resting on two soils of differing deformability

This paper shows the analysis results obtained from more than 200 finite element method (FEM) models used to calculate the settlement of a foundation resting on two soils of differing deformability. The analysis considers such different parameters as the foundation geometry, the percentage of each soil in contact with the foundation base and the ratio of the soils’ elastic moduli. From the described analysis, it is concluded that the maximum settlement of the foundation, calculated by assuming that the foundation is completely resting on the most deformable soil, can be correlated with the settlement calculated by FEM models through a correction coefficient named “settlement reduction factor” (α). As a consequence, a novel expression is proposed for calculating the real settlement of a foundation resting on two soils of different deformability with maximum errors lower than 1.57%, as demonstrated by the statistical analysis carried out. A guide for the application of the proposed simple method is also explained in the paper. Finally, the proposed methodology has been validated using settlement data from an instrumented foundation, indicating that this is a simple, reliable and quick method which allows the computation of the maximum elastic settlement of a raft foundation, evaluates its suitability and optimises its selection process.

Seismically-Induced Rock-Slope Failure: Numerical Investigation using the Bonded Particle Model

Thesis (Ph.D.)--University of Washington, 2016-06

Prediction of the growth of wear-type rail corrugation

The growth behaviour of the vibrational wear phenomenon known as rail corrugation is investigated analytically and numerically using mathematical models. A simplified feedback model for wear-type rail corrugation that includes a wheel pass time delay is developed with an aim to analytically distil the most critical interaction occurring between the wheel/rail structural dynamics, rolling contact mechanics and rail wear. To this end, a stability analysis on the complete system is performed to determine the growth of wear-type rail corrugations over multiple wheelset passages. This analysis indicates that although the dynamical behaviour of the system is stable for each wheel passage, over multiple wheelset passages, the growth of wear-type corrugations is shown to be the result of instability due to feedback interaction between the three primary components of the model. The corrugations are shown analytically to grow for all realistic railway parameters. From this analysis an analytical expression for the exponential growth rate of corrugations in terms of known parameters is developed. This convenient expression is used to perform a sensitivity analysis to identify critical parameters that most affect corrugation growth. The analytical predictions are shown to compare well with results from a benchmarked time-domain finite element model. (C) 2004 Elsevier B.V. All rights reserved.

Effect of granule morphology on breakage behaviour during compression

In the area of dry particle breakage, Discrete Element Method (DEM) simulations have been widely used to analyse the sensitivity of various physical parameters to the behaviour of agglomerates during breakage. This paper looks at the effect of agglomerate shape and structure on the mechanisms and extent of breakage of dry agglomerates under compressive load using DEM simulations. In the simulations, a spherical-shaped agglomerate produced within the DEM code is compared with an irregularly shaped agglomerate, whose structure is that of an actual granule that was characterised with X-ray microtomography (muCT). Both agglomerates have identical particle size distribution, coordination number and surface energy values, with only the agglomerate shape and structure differing between the two. The work here details the breakage behaviour with a number of traditional DEM output parameters (i.e., contact/cluster distributions) with showing vastly different behaviour between the two agglomerates. (C) 2004 Elsevier B.V. All rights reserved.

The igneous charnockite—high-K alkali-calcic I-type granite—incipient charnockite association in Trivandrum Block, southern India

The Pan-African (640 Ma) Chengannoor granite intrudes the NW margin of the Neoproterozoic high-grade metamorphic terrain of the Trivandrum Block (TB), southern India, and is spatially associated with the Cardamom hills igneous charnockite massif (CM). Geochemical features characterize the Chengannoor granite as high-K alkali-calcic I-type granite. Within the constraints imposed by the high temperature, anhydrous, K-rich nature of the magmas, comparison with recent experimental studies on various granitold source compositions, and trace- and rare-earth-element modelling, the distinctive features of the Chengannoor granite reflect a source rock of igneous charnockitic nature. A petrogenetic model is proposed whereby there was a period of basaltic underplating; the partial melting of this basaltic lower crust formed the CM charnockites. The Chengannoor granite was produced by the partial melting of the charnoenderbites from the CM, with subsequent fractionation dominated by feldspars. In a regional context, the Chengannoor I-type granite is considered as a possible heat source for the near-UHT nature of metamorphism in the northern part of the TB. This is different from previous studies, which favoured CM charnockite as the major heat source. The Occurrence of incipient charnockites (both large scale as well as small scale) adjacent to the granite as well as pegmatites (which contain CO2, CO2-H2O, F and other volatiles), suggests that the fluids expelled from the alkaline magma upon solidification generated incipient charnockites through fluid-induced lowering of water activity. Thus the granite and associated alkaline pegmatites acted as conduits for the transfer of heat and volatiles in the Achankovil Shear Zone area, causing pervasive as well as patchy charnockite formation. The transport Of CO2 by felsic melts through the southern Indian middle crust is suggested to be part of a crustal-scale fluid system that linked mantle heat and CO2 input with upward migration of crustally derived felsic melts and incipient charnockite formation, resulting in an igneous charnockite - I-type granite - incipient charnockite association.

A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics

The Lattice Solid Model has been used successfully as a virtual laboratory to simulate fracturing of rocks, the dynamics of faults, earthquakes and gouge processes. However, results from those simulations show that in order to make the next step towards more realistic experiments it will be necessary to use models containing a significantly larger number of particles than current models. Thus, those simulations will require a greatly increased amount of computational resources. Whereas the computing power provided by single processors can be expected to increase according to Moore's law, i.e., to double every 18-24 months, parallel computers can provide significantly larger computing power today. In order to make this computing power available for the simulation of the microphysics of earthquakes, a parallel version of the Lattice Solid Model has been implemented. Benchmarks using large models with several millions of particles have shown that the parallel implementation of the Lattice Solid Model can achieve a high parallel-efficiency of about 80% for large numbers of processors on different computer architectures.

Volcanic resurfacing and the early terrestrial crust: Zircon U-Pb and REE constraints from the Isua Greenstone Belt, southern West Greenland

The oldest known bona fide succession of elastic metasediments Occurs in the Isua Greenstone Belt. SW Greenland and consists of a variety of mica schists and rare metaconglomerates. The metasediments are in direct contact with a felsic metavolcanic lithology that has previously been dated to 3.71 Ga. Based on trace element geochemical data for 30 metasediments, we selected the six samples with highest Zr concentrations for zircon extraction. These samples all yielded very few or no zircon, Those extracted from mica schists yielded ion probe U/Pb ages between 3.70 and 3,71 Ga. One metaconglomerate sample yielded just a single zircon of 3.74 Ga age. The mica schist hosted zircons have U/Pb ages. Th/U ratios, REE patterns and Eu anomalies indistinguishable from zircon in the adjacent 3.71 Ga felsic metavolcanic unit. Trace element modelling requires the bulk of material in the metasediments to be derived from variably weathered mafic lithologies but some metasediments contain substantial contribution from more evolved source lithologies. The paucity of zircon in the mica schists is thus explained by incorporation of material from largely zircon-free volcanic lithologies. The absence of older zircon in the mica schists and the preponderance of mafic source material imply intense, mainly basaltic resurfacing of the early Earth. The implications of this process are discussed, Thermal considerations suggest that horizontal growth of Hadean crust by addition of mafic ultramafic lavas must have triggered self-reorganisation of the protocrust by remelting. Reworking oft Hadean crust may have been aided by burial of hydrated (weathered) metabasalt due to semi-continuous addition of new voluminous basalt Outpouring,;, This process Causes a bias towards eruption of Zr-saturated partial melts at the surface with O-isotope corn posit ion,, potentially different from the mantle. The oldest zircons hosted in sediments would have been buried to substantial depth or formed in plutons that crystallised at some depth from which it took hundreds of millions of years for them to be exhumed and incorporated into much younger sediments. (C) 2005 Elsevier B.V.All rights reserved.

Parallel 3D Simulation of a Fault Gouge using the Lattice Solid Model

Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained.

The mechanical behaviour of cylindrical GRP panels

Glass reinforced plastic (GRP) is now an established material for the fabrication of sonar windows. Its good mechanical strength, light weight, resistance to corrosion and acoustic transparency, are all properties which fit it for this application. This thesis describes a study, undertaken at the Royal Naval Engineering College, Plymouth, into the mechanical behaviour of a circular cylindrical sonar panel. This particular type of panel would be used to cover a flank array sonar in a ship or submarine. The case considered is that of a panel with all of its edges mechanically clamped and subject to pressure loading on its convex surface. A comprehensive program of testing, to determine the orthotropic elastic properties of the laminated composite panel material is described, together with a series of pressure tests on 1:5 scale sonar panels. These pressure tests were carried out in a purpose designed test rig, using air pressure to provide simulated hydrostatic and hydrodynamic loading. Details of all instrumentation used in the experimental work are given in the thesis. The experimental results from the panel testing are compared with predictions of panel behaviour obtained from both the Galerkin solution of Flugge's cylindrical shell equations (orthotropic case), and finite element modelling of the panels using PAFEC. A variety of appropriate panel boundary conditions are considered in each case. A parametric study, intended to be of use as a preliminary design tool, and based on the above Galerkin solution, is also presented. This parametric study considers cases of boundary conditions, material properties, and panel geometry, outside of those investigated in the experimental work Final conclusions are drawn and recommendations made regarding possible improvements to the procedures for design, manufacture and fixing of sonar panels in the Royal Navy.

Micromechanics of agglomerate damage processes

This thesis reports a detailed investigation of the micromechanics of agglomerate behaviour under free-fall impact, double (punch) impact and diametrical compression tests using the simulation software TRUBAL. The software is based on the discrete element method (DEM) which incorporates the Newtonian equations of motion and contact mechanics theory to model the interparticle interactions. Four agglomerates have been used: three dense (differing in interface energy and contact density) and one loose. Although the simulated agglomerates are relatively coarse-grained, the results obtained are in good agreement with laboratory test results reported in the literature. The computer simulation results show that, in all three types of test, the loose agglomerate cannot fracture as it is unable to store sufficient elastic energy. Instead, it becomes flattened for low loading-rates and shattered or crushed at higher loading-rates. In impact tests, the dense agglomerates experience only local damage at low impact velocities. Semi-brittle fracture and fragmentation are produced over a range of higher impact velocities and at very high impact velocities shattering occurs. The dense agglomerates fracture in two or three large fragments in the diametrical compression tests. Local damage at the agglomerate-platen interface always occurs prior to fracture and consists of local bond breakage (microcrack formation) and local dislocations (compaction). The fracture process is dynamic and much more complex than that suggested by continuum fracture mechanics theory. Cracks are always initiated from the contact zones and propagate towards the agglomerate centre. Fracture occurs a short time after the start of unloading when a fracture crack "selection" process takes place. The detailed investigation of the agglomerate damage processes includes an examination of the evolution of the fracture surface. Detailed comparisons of the behaviour of the same agglomerate in all three types of test are presented. The particle size distribution curves of the debris are also examined, for both free-fall and double impact tests.

Eulerian-Lagrangian two phase debris flow model

The main objective of this work is to develop a quasi three-dimensional numerical model to simulate stony debris flows, considering a continuum fluid phase, composed by water and fine sediments, and a non-continuum phase including large particles, such as pebbles and boulders. Large particles are treated in a Lagrangian frame of reference using the Discrete Element Method, the fluid phase is based on the Eulerian approach, using the Finite Element Method to solve the depth-averaged Navier-Stokes equations in two horizontal dimensions. The particle’s equations of motion are in three dimensions. The model simulates particle-particle collisions and wall-particle collisions, taking into account that particles are immersed in a fluid. Bingham and Cross rheological models are used for the continuum phase. Both formulations provide very stable results, even in the range of very low shear rates. Bingham formulation is better able to simulate the stopping stage of the fluid when applied shear stresses are low. Results of numerical simulations have been compared with data from laboratory experiments on a flume-fan prototype. Results show that the model is capable of simulating the motion of big particles moving in the fluid flow, handling dense particulate flows and avoiding overlap among particles. An application to simulate debris flow events that occurred in Northern Venezuela in 1999 shows that the model could replicate the main boulder accumulation areas that were surveyed by the USGS. Uniqueness of this research is the integration of mud flow and stony debris movement in a single modeling tool that can be used for planning and management of debris flow prone areas.

Eulerian-Lagrangian Two Phase Debris Flow Model

The main objective of this work is to develop a quasi three-dimensional numerical model to simulate stony debris flows, considering a continuum fluid phase, composed by water and fine sediments, and a non-continuum phase including large particles, such as pebbles and boulders. Large particles are treated in a Lagrangian frame of reference using the Discrete Element Method, the fluid phase is based on the Eulerian approach, using the Finite Element Method to solve the depth-averaged Navier–Stokes equations in two horizontal dimensions. The particle’s equations of motion are in three dimensions. The model simulates particle-particle collisions and wall-particle collisions, taking into account that particles are immersed in a fluid. Bingham and Cross rheological models are used for the continuum phase. Both formulations provide very stable results, even in the range of very low shear rates. Bingham formulation is better able to simulate the stopping stage of the fluid when applied shear stresses are low. Results of numerical simulations have been compared with data from laboratory experiments on a flume-fan prototype. Results show that the model is capable of simulating the motion of big particles moving in the fluid flow, handling dense particulate flows and avoiding overlap among particles. An application to simulate debris flow events that occurred in Northern Venezuela in 1999 shows that the model could replicate the main boulder accumulation areas that were surveyed by the USGS. Uniqueness of this research is the integration of mud flow and stony debris movement in a single modeling tool that can be used for planning and management of debris flow prone areas.

Piezoelectric energy harvesting using blood-flow-like excitation for implantable devices

The goal of this research is to produce a system for powering medical implants to increase the lifetime of the implanted devices and reduce the battery size. The system consists of a number of elements – the piezoelectric material for generating power, the device design, the circuit for rectification and energy storage. The piezoelectric material is analysed and a process for producing a repeatable high quality piezoelectric material is described. A full width half maximum (FWHM) of the rocking curve X-Ray diffraction (XRD) scan of between ~1.5° to ~1.7° for test wafers was achieved. This is state of the art for AlN on silicon and means devices with good piezoelectric constants can be fabricated. Finite element modelling FEM) was used to design the structures for energy harvesting. The models developed in this work were established to have an accuracy better than 5% in terms of the difference between measured and modelled results. Devices made from this material were analysed for power harvesting ability as well as the effect that they have on the flow of liquid which is an important consideration for implantable devices. The FEM results are compared to experimental results from laser Doppler vibrometry (LDV), magnetic shaker and perfusion machine tests. The rectifying circuitry for the energy harvester was also investigated. The final solution uses multiple devices to provide the power to augment the battery and so this was a key feature to be considered. Many circuits were examined and a solution based on a fully autonomous circuit was advanced. This circuit was analysed for use with multiple low power inputs similar to the results from previous investigations into the energy harvesting devices. Polymer materials were also studied for use as a substitute for the piezoelectric material as well as the substrate because silicon is more brittle.

Strength of the Iberian intraplate lithosphere: Cenozoic deformations and seismicity

Previous studies about the strength of the lithosphere in the Iberia centre fail to resolve the depth of earthquakes because of the rheological uncertainties. Therefore, new contributions are considered (the crustal structure from a density model) and several parameters (tectonic regime, mantle rheology, strain rate) are checked in this paper to properly examine the role of lithospheric strength in the intraplate seismicity and the Cenozoic evolution. The strength distribution with depth, the integrated strength, the effective elastic thickness and the seismogenic thickness have been calculated by a finite element modelling of the lithosphere across the Central System mountain range and the bordering Duero and Madrid sedimentary basins. Only a dry mantle under strike-slip/extension and a strain rate of 10-15 s-1, or under extension and 10-16 s-1, causes a strong lithosphere. The integrated strength and the elastic thickness are lower in the mountain chain than in the basins. These anisotropies have been maintained since the Cenozoic and determine the mountain uplift and the biharmonic folding of the Iberian lithosphere during the Alpine deformations. The seismogenic thickness bounds the seismic activity in the upper–middle crust, and the decreasing crustal strength from the Duero Basin towards the Madrid Basin is related to a parallel increase in Plio–Quaternary deformations and seismicity. However, elasto–plastic modelling shows that current African–Eurasian convergence is resolved elastically or ductilely, which accounts for the low seismicity recorded in this region.