Anisotropic damage mechanics as a novel approach to improve pre-and post-failure borehole stability analysis

Anisotropic damage distribution and evolution have a profound effect on borehole stress concentrations. Damage evolution is an irreversible process that is not adequately described within classical equilibrium thermodynamics. Therefore, we propose a constitutive model, based on non-equilibrium thermodynamics, that accounts for anisotropic damage distribution, anisotropic damage threshold and anisotropic damage evolution. We implemented this constitutive model numerically, using the finite element method, to calculate stress–strain curves and borehole stresses. The resulting stress–strain curves are distinctively different from linear elastic-brittle and linear elastic-ideal plastic constitutive models and realistically model experimental responses of brittle rocks. We show that the onset of damage evolution leads to an inhomogeneous redistribution of material properties and stresses along the borehole wall. The classical linear elastic-brittle approach to borehole stability analysis systematically overestimates the stress concentrations on the borehole wall, because dissipative strain-softening is underestimated. The proposed damage mechanics approach explicitly models dissipative behaviour and leads to non-conservative mud window estimations. Furthermore, anisotropic rocks with preferential planes of failure, like shales, can be addressed with our model.

Comparison of moped, scooter and motorcycle crash risk and crash severity

The increased popularity of mopeds and motor scooters in Australia and elsewhere in the last decade has contributed substantially to the greater use of powered two-wheelers (PTWs) as a whole. As the exposure of mopeds and scooters has increased, so too has the number of reported crashes involving those PTW types, but there is currently little research comparing the safety of mopeds and, particularly, larger scooters with motorcycles. This study compared the crash risk and crash severity of motorcycles, mopeds and larger scooters in Queensland, Australia. Comprehensive data cleansing was undertaken to separate motorcycles, mopeds and larger scooters in police-reported crash data covering the five years to 30 June 2008. The crash rates of motorcycles (including larger scooters) and mopeds in terms of registered vehicles were similar over this period, although the moped crash rate showed a stronger downward trend. However, the crash rates in terms of distance travelled were nearly four times higher for mopeds than for motorcycles (including larger scooters). More comprehensive distance travelled data is needed to confirm these findings. The overall severity of moped and scooter crashes was significantly lower than motorcycle crashes but an ordered probit regression model showed that crash severity outcomes related to differences in crash characteristics and circumstances, rather than differences between PTW types per se. Greater motorcycle crash severity was associated with higher (>80 km/h) speed zones, horizontal curves, weekend, single vehicle and nighttime crashes. Moped crashes were more severe at night and in speed zones of 90 km/h or more. Larger scooter crashes were more severe in 70 km/h zones (than 60 km/h zones) but not in higher speed zones, and less severe on weekends than on weekdays. The findings can be used to inform potential crash and injury countermeasures tailored to users of different PTW types.

Cold-formed steel columns subject to local buckling at elevated temperatures

Fire safety design of building structures has received greater attention in recent times due to continuing losses of properties and lives in fires. However, the structural behaviour of thin-walled cold-formed steel columns under fire conditions is not well understood despite the increasing use of light gauge steels in building construction. Cold-formed steel columns are often subject to local buckling effects. Therefore a series of laboratory tests of lipped and unlipped channel columns made of varying steel thicknesses and grades was undertaken at uniform elevated temperatures up to 700°C under steady state conditions. Finite element models of the tested columns were also developed, and their elastic buckling and nonlinear analysis results were compared with test results at elevated temperatures. Effects of the degradation of mechanical properties of steel with temperature were included in the finite element analyses. The use of accurately measured yield stress, elasticity modulus and stress-strain curves at elevated temperatures provided a good comparison of the ultimate loads and load-deflection curves from tests and finite element analyses. The commonly used effective width design rules and the direct strength method at ambient temperature were then used to predict the ultimate loads at elevated temperatures by using the reduced mechanical properties. By comparing these predicted ultimate loads with those from tests and finite element analyses, the accuracy of using this design approach was evaluated.

Development of improved fire design rules for cold-formed steel wall systems

Traditionally the fire resistance rating of LSF wall systems is based on approximate prescriptive methods developed using limited fire tests. Therefore a detailed research study into the performance of load bearing LSF wall systems under standard fire conditions was undertaken to develop improved fire design rules. It used the extensive fire performance results of eight different LSF wall systems from a series of full scale fire tests and numerical studies for this purpose. The use of previous fire design rules developed for LSF walls subjected to non-uniform elevated temperature distributions based on AISI design manual and Eurocode3 Parts 1.2 and 1.3 was investigated first. New simplified fire design rules based on AS/NZS 4600, North American Specification and Eurocode 3 Part 1.3 were then proposed in this study with suitable allowances for the interaction effects of compression and bending actions. The importance of considering thermal bowing, magnified thermal bowing and neutral axis shift in the fire design was also investigated. A spread sheet based design tool was developed based on the new design rules to predict the failure load ratio versus time and temperature curves for varying LSF wall configurations. The accuracy of the proposed design rules was verified using the test and FEA results for different wall configurations, steel grades, thicknesses and load ratios. This paper presents the details and results of this study including the improved fire design rules for predicting the load capacity of LSF wall studs and the failure times of LSF walls under standard fire conditions.

Mathematical modelling of LiFePO4 cathodes

LiFePO4 is a commercially available battery material with good theoretical discharge capacity, excellent cycle life and increased safety compared with competing Li-ion chemistries. It has been the focus of considerable experimental and theoretical scrutiny in the past decade, resulting in LiFePO4 cathodes that perform well at high discharge rates. This scrutiny has raised several questions about the behaviour of LiFePO4 material during charge and discharge. In contrast to many other battery chemistries that intercalate homogeneously, LiFePO4 can phase-separate into highly and lowly lithiated phases, with intercalation proceeding by advancing an interface between these two phases. The main objective of this thesis is to construct mathematical models of LiFePO4 cathodes that can be validated against experimental discharge curves. This is in an attempt to understand some of the multi-scale dynamics of LiFePO4 cathodes that can be difficult to determine experimentally. The first section of this thesis constructs a three-scale mathematical model of LiFePO4 cathodes that uses a simple Stefan problem (which has been used previously in the literature) to describe the assumed phase-change. LiFePO4 crystals have been observed agglomerating in cathodes to form a porous collection of crystals and this morphology motivates the use of three size-scales in the model. The multi-scale model developed validates well against experimental data and this validated model is then used to examine the role of manufacturing parameters (including the agglomerate radius) on battery performance. The remainder of the thesis is concerned with investigating phase-field models as a replacement for the aforementioned Stefan problem. Phase-field models have recently been used in LiFePO4 and are a far more accurate representation of experimentally observed crystal-scale behaviour. They are based around the Cahn-Hilliard-reaction (CHR) IBVP, a fourth-order PDE with electrochemical (flux) boundary conditions that is very stiff and possesses multiple time and space scales. Numerical solutions to the CHR IBVP can be difficult to compute and hence a least-squares based Finite Volume Method (FVM) is developed for discretising both the full CHR IBVP and the more traditional Cahn-Hilliard IBVP. Phase-field models are subject to two main physicality constraints and the numerical scheme presented performs well under these constraints. This least-squares based FVM is then used to simulate the discharge of individual crystals of LiFePO4 in two dimensions. This discharge is subject to isotropic Li+ diffusion, based on experimental evidence that suggests the normally orthotropic transport of Li+ in LiFePO4 may become more isotropic in the presence of lattice defects. Numerical investigation shows that two-dimensional Li+ transport results in crystals that phase-separate, even at very high discharge rates. This is very different from results shown in the literature, where phase-separation in LiFePO4 crystals is suppressed during discharge with orthotropic Li+ transport. Finally, the three-scale cathodic model used at the beginning of the thesis is modified to simulate modern, high-rate LiFePO4 cathodes. High-rate cathodes typically do not contain (large) agglomerates and therefore a two-scale model is developed. The Stefan problem used previously is also replaced with the phase-field models examined in earlier chapters. The results from this model are then compared with experimental data and fit poorly, though a significant parameter regime could not be investigated numerically. Many-particle effects however, are evident in the simulated discharges, which match the conclusions of recent literature. These effects result in crystals that are subject to local currents very different from the discharge rate applied to the cathode, which impacts the phase-separating behaviour of the crystals and raises questions about the validity of using cathodic-scale experimental measurements in order to determine crystal-scale behaviour.

Drying kinetics of two grape varieties of Italy

Two varieties of grapes, white grape and red grape grown in the Campania region of Italy were selected for the study of drying characteristics. Comparisons were made with treated and untreated grapes under constant drying condition of 50o C in a conventional drying system. This temperature was selected to represent farm drying conditions. Grapes were purchased from a local market from the same supplier to maintain the same size of grapes and same properties. An abrasive physical treatment was used as pretreatment. The drying curves were constructed and drying kinetics was calculated using several commonly available models. It was found that treated samples show better drying characteristics than untreated samples. The objective of this study is to obtain drying kinetics which can be used to optimize the drying operations in grape drying.

Segmental torso masses and gravity-induced coronal plane joint moments in adolescent idiopathic scoliosis

Introduction: Calculating segmental (vertebral level-by-level) torso masses in Adolescent Idiopathic Scoliosis (AIS) patients allows the gravitational loading on the scoliotic spine during relaxed standing to be estimated. This study used supine CT scans of AIS patients to measure segmental torso masses and explored the joint moments in the coronal plane, particularly at the apex of a scoliotic major curve. Methods: Existing low dose CT data from the Paediatric Spine Research Group was used to calculate vertebral level-by-level torso masses and joint moments occurring in the spine for a group of 20 female AIS patients with right sided thoracic curves. The mean age was 15.0 ± 2.7 years and all curves were classified Lenke Type 1 with a mean Cobb angle 52 ± 5.9°. Image processing software, ImageJ (v1.45 NIH USA) was used to create reformatted coronal plane images, reconstruct vertebral level-by-level torso segments and subsequently measure the torso volume corresponding to each vertebral level. Segment mass was then determined by assuming a tissue density of 1.04x103 kg/m3. Body segment masses for the head, neck and arms were taken from published anthropometric data (Winter 2009). Intervertebral joint moments in the coronal plane at each vertebral level were found from the position of the centroid of the segment masses relative to the joint centres with the segmental body mass data. Results and Discussion: The magnitude of the torso masses from T1-L5 increased inferiorly, with a 150% increase in mean segmental torso mass from 0.6kg at T1 to 1.5kg at L5. The magnitudes of the calculated coronal plane joint moments during relaxed standing were typically 5-7 Nm at the apex of the curve, with the highest apex joint torque of 7Nm. The CT scans were performed in the supine position and curve magnitudes are known to be 7-10° smaller than those measured in standing, due to the absence of gravity acting on the spine. Hence, it can be expected that the moments produced by gravity in the standing individual will be greater than those calculated here.

Simultaneous determination of 6,7-dimethylesculetin and geniposide in rat plasma and its application to pharmacokinetic studies of Yin Chen Hao Tang preparation

A method for the rapid and simultaneous determination of 6,7-dimethylesculetin (CAS 120-08-1) and geniposide (CAS 24512-63-8) in rat plasma has been developed, using validated high performance liquid chromatography (HPLC) with solid phase extraction (SPE). The HPLC analysis was performed on a commercially available column (200 mm x 4.6 mm, 5 microm) with acetonitrile-methanol-0.1% aqueous formic acid as mobile phase and the UV detection at 343 nm and 238 nm for 6,7-dimethylesculetin and geniposide, respectively. The calibration curves for 6,7-dimethylesculetin and geniposide were linear over the range 0.4-25.6 microg/mL and 1.12-71.68 microg/mL, respectively. The lower limits of quantitation were 0.40 microg/ mL and 1.12 microg/mL, and the lower limits of detection were 0.06 microg/mL and 0.09 microg/ mL, respectively. The intra-day and inter-day precision for 6,7-dimethylesculetin and geniposide were < 5%, whereas the absolute recovery percentages were > 74%. A successful application of the developed HPLC analysis was demonstrated for the pharmacokinetic study of a Traditional Chinese Medicine formula of Yin Chen Hao Tang preparation.

SHOX gene is expressed in vertebral body growth plates in idiopathic and congenital scoliosis : implications for the etiology of scoliosis in Turner Syndrome

Reduced SHOX gene expression has been demonstrated to be associated with all skeletal abnormalities in Turner syndrome, other than scoliosis (and kyphosis). There is evidence to suggest that Turner syndrome scoliosis is clinically and radiologically similar to idiopathic scoliosis, although the phenotypes are dissimilar. This pilot gene expression study used relative quantitative real-time PCR (qRT-PCR) of the SHOX (short stature on X) gene to determine whether it is expressed in vertebral body growth plates in idiopathic and congenital scoliosis. After vertebral growth plate dissection, tissue was examined histologically and RNA was extracted and its integrity was assessed using a Bio-Spec Mini, NanoDrop ND-1000 spectrophotometer and standard denaturing gel electrophoresis. Following cDNA synthesis, gene-specific optimization in a Corbett RotorGene 6000 real-time cycler was followed by qRT-PCR of vertebral tissue. Histological examination of vertebral samples confirmed that only growth plate was analyzed for gene expression. Cycling and melt curves were resolved in triplicate for all samples. SHOX abundance was demonstrated in congenital and idiopathic scoliosis vertebral body growth plates. SHOX expression was 11-fold greater in idiopathic compared to congenital (n = 3) scoliosis (p = 0.027). This study confirmed that SHOX was expressed in vertebral body growth plates, which implies that its expression may also be associated with the scoliosis (and kyphosis) of Turner syndrome. SHOX expression is reduced in Turner syndrome (short stature). In this study, increased SHOX expression was demonstrated in idiopathic scoliosis (tall stature) and congenital scoliosis.

The measurement of adenosine and estrogen receptor expression in rat brains following ovariectomy using quantitative PCR analysis

In our laboratory we have developed a quantitative-polymerase chain reaction (Q-PCR) strategy to examine the differential expression of adenosine receptor (ADOR), A(1), A(2A), A(2B) and A(3), and estrogen receptors (ER) alpha and beta. Brain and uterine mRNA were first used to optimise specific amplification conditions prior to SYBR Green I real time analysis of receptor subtype expression. SYBR Green I provided a convenient and sensitive means of examining specific PCR amplification product in real time, and allowed the generation of standard curves from which relative receptor abundance could be determined. Real time Q-PCR analysis was then performed, to examine changes in receptor expression levels in brains of adult female Wistar rats 3-month post ovariectomy. Comparison with sham-operated age-matched control rats demonstrated both comparative and absolute-copy number changes in receptor levels. Evaluation of both analytical methods investigated 18S rRNA as an internal reference for comparative gene expression analysis in the brain. The results of this study revealed preferential repression of ADORA(2A) (>4-fold down) and consistent (>2-fold) down-regulation of ADORA(1), ADORA(3), and ER-beta, following ovariectomy. No change was found in ADORA(2B) or ER-alpha. Analysis of absolute copy number in this study revealed a correlation between receptor expression in response to ovariectomy, and relative receptor subtype abundance in the brain.

Ultimate load behaviour of box-columns under combined loading of axial compression and torsion

This paper describes a study of the theoretical and experimental behaviour of box-columns of varying b/t ratios under loadings of axial compression and torsion and their combinations. Details of the testing rigs and the testing methods, the results obtained such as the load-deflection curves and the interaction diagrams, and experimental observations regarding the behaviour of box-models and the types of local plastic mechanisms associated with each type of loading are presented. A simplified rigid-plastic analysis is carried out to study the collapse behaviour of box-columns under these loadings, based on the observed plastic mechanisms, and the results are compared with those of experiments.

Locked nucleic acid (LNA) single nucleotide polymorphism (SNP) genotype analysis and validation using real-time PCR

With an increased emphasis on genotyping of single nucleotide polymorphisms (SNPs) in disease association studies, the genotyping platform of choice is constantly evolving. In addition, the development of more specific SNP assays and appropriate genotype validation applications is becoming increasingly critical to elucidate ambiguous genotypes. In this study, we have used SNP specific Locked Nucleic Acid (LNA) hybridization probes on a real-time PCR platform to genotype an association cohort and propose three criteria to address ambiguous genotypes. Based on the kinetic properties of PCR amplification, the three criteria address PCR amplification efficiency, the net fluorescent difference between maximal and minimal fluorescent signals and the beginning of the exponential growth phase of the reaction. Initially observed SNP allelic discrimination curves were confirmed by DNA sequencing (n = 50) and application of our three genotype criteria corroborated both sequencing and observed real-time PCR results. In addition, the tested Caucasian association cohort was in Hardy-Weinberg equilibrium and observed allele frequencies were very similar to two independently tested Caucasian association cohorts for the same tested SNP. We present here a novel approach to effectively determine ambiguous genotypes generated from a real-time PCR platform. Application of our three novel criteria provides an easy to use semi-automated genotype confirmation protocol.

Aligned poly(L-lactic-co-ε-caprolactone) electrospun microfibers and knitted structure: a novel composite scaffold for ligament tissue engineering

We developed a novel technique involving knitting and electrospinning to fabricate a composite scaffold for ligament tissue engineering. Knitted structures were coated with poly(L-lactic-co-e-caprolactone) (PLCL) and then placed onto a rotating cylinder and a PLCL solution was electrospun onto the structure. Highly aligned 2-μm-diameter microfibers covered the space between the stitches and adhered to the knitted scaffolds. The stress–strain tensile curves exhibited an initial toe region similar to the tensile behavior of ligaments. Composite scaffolds had an elastic modulus (150 ± 14 MPa) similar to the modulus of human ligaments. Biological evaluation showed that cells proliferated on the composite scaffolds and they spontaneously orientated along the direction of microfiber alignment. The microfiber architecture also induced a high level of extracellular matrix secretion, which was characterized by immunostaining. We found that cells produced collagen type I and type III, two main components found in ligaments. After 14 days of culture, collagen type III started to form a fibrous network. We fabricated a composite scaffold having the mechanical properties of the knitted structure and the morphological properties of the aligned microfibers. It is difficult to seed a highly macroporous structure with cells, however the technique we developed enabled an easy cell seeding due to presence of the microfiber layer. Therefore, these scaffolds presented attractive properties for a future use in bioreactors for ligament tissue engineering.

Understanding mild acid pretreatment of sugarcane bagasse through particle scale modeling

Sugarcane bagasse is an abundant and sustainable resource, generated as a by-product of sugarcane milling. The cellulosic material within bagasse can be broken down into glucose molecules and fermented to produce ethanol, making it a promising feedstock for biofuel production. Mild acid pretreatment hydrolyses the hemicellulosic component of biomass, thus allowing enzymes greater access to the cellulosic substrate during saccharification. A particle-scale mathematical model describing the mild acid pretreatment of sugarcane bagasse has been developed, using a volume averaged framework. Discrete population-balance equations are used to characterise the polymer degradation kinetics, and diffusive effects account for mass transport within the cell wall of the bagasse. As the fibrous material hydrolyses over time, variations in the porosity of the cell wall and the downstream effects on the reaction kinetics are accounted for using conservation of volume arguments. Non-dimensionalization of the model equations reduces the number of parameters in the system to a set of four dimensionless ratios that compare the timescales of different reaction and diffusion events. Theoretical yield curves are compared to macroscopic experimental observations from the literature and inferences are made as to constraints on these “unknown” parameters. These results enable connections to be made between experimental data and the underlying thermodynamics of acid pretreatment. Consequently, the results suggest that data-fitting techniques used to obtain kinetic parameters should be carefully applied, with prudent consideration given to the chemical and physiological processes being modeled.

Real - time prediction of rainfall induced instability of residual soil slopes associated with deep cracks

The early warning based on real-time prediction of rain-induced instability of natural residual slopes helps to minimise human casualties due to such slope failures. Slope instability prediction is complicated, as it is influenced by many factors, including soil properties, soil behaviour, slope geometry, and the location and size of deep cracks in the slope. These deep cracks can facilitate rainwater infiltration into the deep soil layers and reduce the unsaturated shear strength of residual soil. Subsequently, it can form a slip surface, triggering a landslide even in partially saturated soil slopes. Although past research has shown the effects of surface-cracks on soil stability, research examining the influence of deep-cracks on soil stability is very limited. This study aimed to develop methodologies for predicting the real-time rain-induced instability of natural residual soil slopes with deep cracks. The results can be used to warn against potential rain-induced slope failures. The literature review conducted on rain induced slope instability of unsaturated residual soil associated with soil crack, reveals that only limited studies have been done in the following areas related to this topic: - Methods for detecting deep cracks in residual soil slopes. - Practical application of unsaturated soil theory in slope stability analysis. - Mechanistic methods for real-time prediction of rain induced residual soil slope instability in critical slopes with deep cracks. Two natural residual soil slopes at Jombok Village, Ngantang City, Indonesia, which are located near a residential area, were investigated to obtain the parameters required for the stability analysis of the slope. A survey first identified all related field geometrical information including slope, roads, rivers, buildings, and boundaries of the slope. Second, the electrical resistivity tomography (ERT) method was used on the slope to identify the location and geometrical characteristics of deep cracks. The two ERT array models employed in this research are: Dipole-dipole and Azimuthal. Next, bore-hole tests were conducted at different locations in the slope to identify soil layers and to collect undisturbed soil samples for laboratory measurement of the soil parameters required for the stability analysis. At the same bore hole locations, Standard Penetration Test (SPT) was undertaken. Undisturbed soil samples taken from the bore-holes were tested in a laboratory to determine the variation of the following soil properties with the depth: - Classification and physical properties such as grain size distribution, atterberg limits, water content, dry density and specific gravity. - Saturated and unsaturated shear strength properties using direct shear apparatus. - Soil water characteristic curves (SWCC) using filter paper method. - Saturated hydraulic conductivity. The following three methods were used to detect and simulate the location and orientation of cracks in the investigated slope: (1) The electrical resistivity distribution of sub-soil obtained from ERT. (2) The profile of classification and physical properties of the soil, based on laboratory testing of soil samples collected from bore-holes and visual observations of the cracks on the slope surface. (3) The results of stress distribution obtained from 2D dynamic analysis of the slope using QUAKE/W software, together with the laboratory measured soil parameters and earthquake records of the area. It was assumed that the deep crack in the slope under investigation was generated by earthquakes. A good agreement was obtained when comparing the location and the orientation of the cracks detected by Method-1 and Method-2. However, the simulated cracks in Method-3 were not in good agreement with the output of Method-1 and Method-2. This may have been due to the material properties used and the assumptions made, for the analysis. From Method-1 and Method-2, it can be concluded that the ERT method can be used to detect the location and orientation of a crack in a soil slope, when the ERT is conducted in very dry or very wet soil conditions. In this study, the cracks detected by the ERT were used for stability analysis of the slope. The stability of the slope was determined using the factor of safety (FOS) of a critical slip surface obtained by SLOPE/W using the limit equilibrium method. Pore-water pressure values for the stability analysis were obtained by coupling the transient seepage analysis of the slope using finite element based software, called SEEP/W. A parametric study conducted on the stability of an investigated slope revealed that the existence of deep cracks and their location in the soil slope are critical for its stability. The following two steps are proposed to predict the rain-induced instability of a residual soil slope with cracks. (a) Step-1: The transient stability analysis of the slope is conducted from the date of the investigation (initial conditions are based on the investigation) to the preferred date (current date), using measured rainfall data. Then, the stability analyses are continued for the next 12 months using the predicted annual rainfall that will be based on the previous five years rainfall data for the area. (b) Step-2: The stability of the slope is calculated in real-time using real-time measured rainfall. In this calculation, rainfall is predicted for the next hour or 24 hours and the stability of the slope is calculated one hour or 24 hours in advance using real time rainfall data. If Step-1 analysis shows critical stability for the forthcoming year, it is recommended that Step-2 be used for more accurate warning against the future failure of the slope. In this research, the results of the application of the Step-1 on an investigated slope (Slope-1) showed that its stability was not approaching a critical value for year 2012 (until 31st December 2012) and therefore, the application of Step-2 was not necessary for the year 2012. A case study (Slope-2) was used to verify the applicability of the complete proposed predictive method. A landslide event at Slope-2 occurred on 31st October 2010. The transient seepage and stability analyses of the slope using data obtained from field tests such as Bore-hole, SPT, ERT and Laboratory tests, were conducted on 12th June 2010 following the Step-1 and found that the slope in critical condition on that current date. It was then showing that the application of the Step-2 could have predicted this failure by giving sufficient warning time.