Cohesive zone modelling of mineralized collagen fibril arrays in bending

Hard biological materials such as bone possess superior material properties of high stiffness and toughness. Two unique characteristics of bone microstructure are a large aspect ratio of mineralized collagen fibrils (MCF), and an extremely thin and large area of extrafibrillar protein matrix located between the MCF. The objective of this study is to investigate the effects of: (1) MCF aspect ratio, and (2) energy dissipation in extrafibrillar protein matrix on the mechanical behaviour of MCF arrays. In this study, notched specimens of MCF arrays in extrafibrillar protein matrix are subjected to bending. Cohesive zone model was implemented to simulate the failure of extrafibrillar protein matrix. The study reveals that the MCF array with a higher MCF aspect ratio and the MCF array with a higher protein energy dissipation in the interface direction are able to sustain a higher bending force and dissipate higher energy.

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.

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.

Modelling sea water intrusion in coastal aquifers using heterogeneous computing

The objective of this PhD research program is to investigate numerical methods for simulating variably-saturated flow and sea water intrusion in coastal aquifers in a high-performance computing environment. The work is divided into three overlapping tasks: to develop an accurate and stable finite volume discretisation and numerical solution strategy for the variably-saturated flow and salt transport equations; to implement the chosen approach in a high performance computing environment that may have multiple GPUs or CPU cores; and to verify and test the implementation. The geological description of aquifers is often complex, with porous materials possessing highly variable properties, that are best described using unstructured meshes. The finite volume method is a popular method for the solution of the conservation laws that describe sea water intrusion, and is well-suited to unstructured meshes. In this work we apply a control volume-finite element (CV-FE) method to an extension of a recently proposed formulation (Kees and Miller, 2002) for variably saturated groundwater flow. The CV-FE method evaluates fluxes at points where material properties and gradients in pressure and concentration are consistently defined, making it both suitable for heterogeneous media and mass conservative. Using the method of lines, the CV-FE discretisation gives a set of differential algebraic equations (DAEs) amenable to solution using higher-order implicit solvers. Heterogeneous computer systems that use a combination of computational hardware such as CPUs and GPUs, are attractive for scientific computing due to the potential advantages offered by GPUs for accelerating data-parallel operations. We present a C++ library that implements data-parallel methods on both CPU and GPUs. The finite volume discretisation is expressed in terms of these data-parallel operations, which gives an efficient implementation of the nonlinear residual function. This makes the implicit solution of the DAE system possible on the GPU, because the inexact Newton-Krylov method used by the implicit time stepping scheme can approximate the action of a matrix on a vector using residual evaluations. We also propose preconditioning strategies that are amenable to GPU implementation, so that all computationally-intensive aspects of the implicit time stepping scheme are implemented on the GPU. Results are presented that demonstrate the efficiency and accuracy of the proposed numeric methods and formulation. The formulation offers excellent conservation of mass, and higher-order temporal integration increases both numeric efficiency and accuracy of the solutions. Flux limiting produces accurate, oscillation-free solutions on coarse meshes, where much finer meshes are required to obtain solutions with equivalent accuracy using upstream weighting. The computational efficiency of the software is investigated using CPUs and GPUs on a high-performance workstation. The GPU version offers considerable speedup over the CPU version, with one GPU giving speedup factor of 3 over the eight-core CPU implementation.

Cell migration and proliferation on homogeneous and non-homogeneous domains : modelling on the scale of individuals and populations

Cell migration is a behaviour critical to many key biological effects, including wound healing, cancerous cell invasion and morphogenesis, the development of an organism from an embryo. However, given that each of these situations is distinctly different and cells are extremely complicated biological objects, interest lies in more basic experiments which seek to remove conflating factors and present a less complex environment within which cell migration can be experimentally examined. These include in vitro studies like the scratch assay or circle migration assay, and ex vivo studies like the colonisation of the hindgut by neural crest cells. The reduced complexity of these experiments also makes them much more enticing as problems to mathematically model, like done here. The primary goal of the mathematical models used in this thesis is to shed light on which cellular behaviours work to generate the travelling waves of invasion observed in these experiments, and to explore how variations in these behaviours can potentially predict differences in this invasive pattern which are experimentally observed when cell types or chemical environment are changed. Relevant literature has already identified the difficulty of distinguishing between these behaviours when using traditional mathematical biology techniques operating on a macroscopic scale, and so here a sophisticated individual-cell-level model, an extension of the Cellular Potts Model (CPM), is been constructed and used to model a scratch assay experiment. This model includes a novel mechanism for dealing with cell proliferations that allowed for the differing properties of quiescent and proliferative cells to be implemented into their behaviour. This model is considered both for its predictive power and used to make comparisons with the travelling waves which result in more traditional macroscopic simulations. These comparisons demonstrate a surprising amount of agreement between the two modelling frameworks, and suggest further novel modifications to the CPM that would allow it to better model cell migration. Considerations of the model’s behaviour are used to argue that the dominant effect governing cell migration (random motility or signal-driven taxis) likely depends on the sort of invasion demonstrated by cells, as easily seen by microscopic photography. Additionally, a scratch assay simulated on a non-homogeneous domain consisting of a ’fast’ and ’slow’ region is also used to further differentiate between these different potential cell motility behaviours. A heterogeneous domain is a novel situation which has not been considered mathematically in this context, nor has it been constructed experimentally to the best of the candidate’s knowledge. Thus this problem serves as a thought experiment used to test the conclusions arising from the simulations on homogeneous domains, and to suggest what might be observed should this non-homogeneous assay situation be experimentally realised. Non-intuitive cell invasion patterns are predicted for diffusely-invading cells which respond to a cell-consumed signal or nutrient, contrasted with rather expected behaviour in the case of random-motility-driven invasion. The potential experimental observation of these behaviours is demonstrated by the individual-cell-level model used in this thesis, which does agree with the PDE model in predicting these unexpected invasion patterns. In the interest of examining such a case of a non-homogeneous domain experimentally, some brief suggestion is made as to how this could be achieved.

Development of optimal designs of insulated rail joints

Proper functioning of Insulated Rail Joints (IRJs) is essential for the safe operation of the railway signalling systems and broken rail identification circuitries. The Conventional IRJ (CIRJ) resembles structural butt joints consisting of two pieces of rails connected together through two joint bars on either side of their web and the assembly is held together through pre-tensioned bolts. As the IRJs should maintain electrical insulation between the two rails, a gap between the rail ends must be retained at all times and all metal contacting surfaces should be electrically isolated from each other using non-conductive material. At the gap, the rail ends lose longitudinal continuity and hence the vertical sections of the rail ends are often severely damaged, especially at the railhead, due to the passage of wheels compared to other continuously welded rail sections. Fundamentally, the reason for the severe damage can be related to the singularities of the wheel-rail contact pressure and the railhead stress. No new generation designs that have emerged in the market to date have focussed on this fundamental; they only have provided attention to either the higher strength materials or the thickness of the sections of various components of the IRJs. In this thesis a novel method of shape optimisation of the railhead is developed to eliminate the pressure and stress singularities through changes to the original sharp corner shaped railhead into an arc profile in the longitudinal direction. The optimal shape of the longitudinal railhead profile has been determined using three nongradient methods in search of accuracy and efficiency: (1) Grid Search Method; (2) Genetic Algorithm Method and (3) Hybrid Genetic Algorithm Method. All these methods have been coupled with a parametric finite element formulation for the evaluation of the objective function for each iteration or generation depending on the search algorithm employed. The optimal shape derived from these optimisation methods is termed as Stress Minimised Railhead (SMRH) in this thesis. This optimal SMRH design has exhibited significantly reduced stress concentration that remains well below the yield strength of the head hardened rail steels and has shifted the stress concentration location away from the critical zone of the railhead end. The reduction in the magnitude and the relocation of the stress concentration in the SMRH design has been validated through a full scale wheel – railhead interaction test rig; Railhead strains under the loaded wheels have been recorded using a non-contact digital image correlation method. Experimental study has confirmed the accuracy of the numerical predications. Although the SMRH shaped IRJs eliminate stress singularities, they can still fail due to joint bar or bolt hole cracking; therefore, another conceptual design, termed as Embedded IRJ (EIRJ) in this thesis, with no joint bars and pre-tensioned bolts has been developed using a multi-objective optimisation formulation based on the coupled genetic algorithm – parametric finite element method. To achieve the required structural stiffness for the safe passage of the loaded wheels, the rails were embedded into the concrete of the post tensioned sleepers; the optimal solutions for the design of the EIRJ is shown to simplify the design through the elimination of the complex interactions and failure modes of the various structural components of the CIRJ. The practical applicability of the optimal shapes SMRH and EIRJ is demonstrated through two illustrative examples, termed as improved designs (IMD1 & IMD2) in this thesis; IMD1 is a combination of the CIRJ and the SMRH designs, whilst IMD2 is a combination of the EIRJ and SMRH designs. These two improved designs have been simulated for two key operating (speed and wagon load) and design (wheel diameter) parameters that affect the wheel-rail contact; the effect of these parameters has been found to be negligible to the performance of the two improved designs and the improved designs are in turn found far superior to the current designs of the CIRJs in terms of stress singularities and deformation under the passage of the loaded wheels. Therefore, these improved designs are expected to provide longer service life in relation to the CIRJs.

Modelling the failure of thin layered mortar joints in masonry

This paper deals with the failure of high adhesive, low compressive strength, thin layered polymer mortar joints in masonry through a contact modelling in finite element framework. Failure due to combined shear, tensile and compressive stresses are considered through a constitutive damaging contact model that incorporates traction–separation as a function of displacement discontinuity. The modelling method is verified using single and multiple contact analyses of thin mortar layered masonry specimens under shear, tensile and compressive stresses and their combinations. Using this verified method, the failure of thin mortar layered masonry under a range of shear to tension ratios and shear to compression ratios has been examined. Finally, this model is applied to thin bed masonry wallettes for their behaviour under biaxial tension–tension and compression–tension loadings perpendicular and parallel to the bed joints.

How to increase the accuracy of analysis and reduce the computational time in ANSYS in the case of deformation study of orthopedic bone plates

Currently, finite element analyses are usually done by means of commercial software tools. Accuracy of analysis and computational time are two important factors in efficiency of these tools. This paper studies the effective parameters in computational time and accuracy of finite element analyses performed by ANSYS and provides the guidelines for the users of this software whenever they us this software for study on deformation of orthopedic bone plates or study on similar cases. It is not a fundamental scientific study and only shares the findings of the authors about structural analysis by means of ANSYS workbench. It gives an idea to the readers about improving the performance of the software and avoiding the traps. The solutions provided in this paper are not the only possible solutions of the problems and in similar cases there are other solutions which are not given in this paper. The parameters of solution method, material model, geometric model, mesh configuration, number of the analysis steps, program controlled parameters and computer settings are discussed through thoroughly in this paper.

Heat and mass transfer modeling of the osmo-convective drying of yacon roots (Smallanthus sonchifolius)

Osmotic treatments are often applied prior to convective drying of foods to impart sensory appeal aspects. During this process a multicomponent mass flow, composed mainly of water and osmotic agent, takes place. In this work, a heat and mass transfer model for the osmo-convective drying of yacon was developed and solved by the Finite Element Method using COMSOL Multiphysics®, considering a 2-D axisymmetric geometry and moisture dependent thermophysical properties. Yacon slices were osmotically dehydrated for 2 hours in a solution of sucralose and then dried in a tray dryer for 3 hours. The model was validated by experimental data of temperature, moisture content and sucralose uptake (R²> 0.90).

Behaviour and design of crest-fixed profiled steel roof claddings under wind uplift

Profiled steel roof claddings in Australia are commonly made of very thin high tensile steel and are crest-fixed with screw fasteners. At present the design of these claddings is entirely based on testing. In order to improve the understanding of the behaviour of these claddings under wind uplift, and thus the design methods, a detailed investigation consisting of a finite element analysis and laboratory experiments was carried out on two-span roofing assemblies of three common roofing profiles. It was found that the failure of the roof cladding system was due to a local failure (dimpling of crests/pull-through) at the fasteners. This paper presents the details of the investigation, the results and then proposes a design method based on the strength of the screwed connections, for which testing of small-scale roofing models and/or using a simple design formula is recommended.

Effect of overload cycles on thin steel roof claddings during cyclonic winds

During an investigation on thin steel roof claddings under simulated cyclonic wind loading, it was found that trapezoidal roof claddings behaved quite differently to corrugated (arc and tangent type) roof claddings due to the presence of overload cycles. The overload cycles caused a reduction in fatigue life for corrugated roofing whereas the reverse occurred for trapezoidal roofing. This contrasting behavior of the two crest-fixed roof claddings was investigated using small scale roofing models instead of the commonly used large scale two-span roof claddings. It was found that overload cycles formed a weaker locally dimpled mechanism around the fastener holes of corrugated roofing and thus accelerated the fatigue-caused pull-through failure. In contrast, a stronger deformed shape was formed in trapezoidal roofing which delayed the pull-through failure. Both laboratory testing and finite element analysis of small scale models were used to study the contrasting behavior of roof claddings.

Buckling experiments on hollow flange beams with web stiffeners

A new cold-formed and resistance welded section known as the Hollow Flange Beam (HFB) has been developed recently in Australia. In contrast to the common lateral torsional buckling mode of I-beams, this unique section comprising two stiff triangular flanges and a slender web is susceptible to a lateral distortional buckling mode of failure involving lateral deflection, twist and cross-section change due to web distortion. This lateral distortional buckling behaviour has been shown to cause significant reduction of the available flexural strength of HFBs. An investigation using finite element analyses and large scale experiments was carried out into the use of transverse web plate stiffeners to improve the lateral buckling capacity of HFBs. This paper presents the details of the experimental investigation, the results, and the final stiffener arrangement whereas the details of the finite element analyses are presented in a companion paper at this conference.

Local plastic mechanisms in thin steel plates under in-plane compression

Thin-walled steel plates subjected to in-plane compression develop two types of local plastic mechanism, namely the roof-shaped mechanism and the so-called flip-disc mechanism, but the intriguing question of why two mechanisms should develop was not answered until recently. It was considered that the location of first yield point shifted from the centre of the plate to the midpoint of the longitudinal edge depending on the b/t ratio, imperfection level, and yield stress of steel, which then decided the type of mechanism. This paper has verified this hypothesis using analysis and laboratory experiments. An elastic analysis using Galerkin's method to solve Marguerre's equations was first used to determine the first yield point, based on which the local plastic mechanism/imperfection tolerance tables have been developed which give the type of mechanism as a function of b/t ratio, imperfection level and yield stress of steel. Laboratory experiments of thin-walled columns verified the imperfection tolerance tables and thus indirectly the hypothesis. Elastic and rigid-plastic curves were them used to predict the effect on the ultimate load due to the change of mechanism. A finite element analysis of selected cases also confirmed the results from simple analyses and experiments.

Buckling experiments on hollow flange beams with web stiffeners

The hollow flange beam (HFB) is a new cold-formed and resistance-welded section developed in Australia. Due to its unique geometry comprising two stiff triangular flanges and a slender web, the HFB is susceptible to a lateral-distortional buckling mode of failure involving web distortion. Investigation using finite-element analyses showed that the use of transverse web plate stiffeners effectively eliminated lateral-distortional buckling of HFBs and thus any associated reduction in flexural capacity. A detailed experimental investigation was then carried out to validate the results from the finite-element analysis and to improve the stiffener configuration further. This led to the development of a special stiffener that is screw-fastened to the flanges on alternate sides of the web. This paper presents the details of the experimental investigations, the results, and the final stiffener arrangement whereas the details of the finite-element analyses are presented in a companion paper.

Flexural capacity of hollow flange beams

The hollow flange beam (HFB) is a unique cold-formed steel section developed in Australia for use as a flexural member. Research has identified that the HFB section's flexural capacity for intermediate span members is limited by lateral distortional buckling, which is characterized by simultaneous lateral deflection, twist, and web distortion. This buckling behaviour is mainly due to the unique geometry of the section, comprising two torsionally stiff triangular flanges connected by a slender web. This paper presents a finite element analytical model suitable for non-linear analysis of HFB flexural members. The model includes all significant effects that may influence the ultimate capacity of such members, including material inelasticity, local buckling, member instability, web distortion, residual stresses, and geometric imperfections. It was found to accurately predict both the elastic lateral distortional buckling moments and the ultimate capacities of HFB flexural members, and was therefore used in the development of design curves and suitable design procedures.