Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Adsorption in mesopores: A molecular-continuum model with application to MCM-41

The classical model of surface layering followed by capillary condensation during adsorption in mesopores, is modified here by consideration of the adsorbate solid interaction potential. The new theory accurately predicts the capillary coexistence curve as well as pore criticality, matching that predicted by density functional theory. The model also satisfactorily predicts the isotherm for nitrogen adsorption at 77.4 K on MCM-41 material of various pore sizes, synthesized and characterized in our laboratory, including the multilayer region, using only data on the variation of condensation pressures with pore diameter. The results indicate a minimum mesopore diameter for the surface layering model to hold as 14.1 Å, below which size micropore filling must occur, and a minimum pore diameter for mechanical stability of the hemispherical meniscus during desorption as 34.2 Å. For pores in-between these two sizes reversible condensation is predicted to occur, in accord with the experimental data for nitrogen adsorption on MCM-41 at 77.4 K.

Thermal effects in the evolution of initially layered mantle material

A simplified model for anisotropic mantle convection based on a novel class of rheologies, originally developed for folding instabilities in multilayered rock (MUHLHAUS et al., 2002), is extended ¨ through the introduction of a thermal anisotropy dependent on the local layering. To examine the effect of the thermal anisotropy on the evolution of mantle material, a parallel implementation of this model was undertaken using the Escript modelling toolkit and the Finley finite-element computational kernel (DAVIES et al., 2004). For the cases studied, there appears too little if any effect. For comparative purposes, the effects of anisotropic shear viscosity and the introduced thermal anisotropy are also presented. These results contribute to the characterization of viscous anisotropic mantle convection subject to variation in thermal conductivities and shear viscosities.

On uniquely determining cell-wall material properties with the compression experiment

A finite element model (FEM) of the cell-compression experiment has been developed in dimensionless form to extract the fundamental cell-wall-material properties (i.e. the constitutive equation and its parameters) from experiment force-displacement data. The FEM simulates the compression of a thin-walled, liquid-filled sphere between two flat surfaces. The cell-wall was taken to be permeable and the FEM therefore accounts for volume loss during compression. Previous models assume an impermeable wall and hence a conserved cell volume during compression. A parametric study was conducted for structural parameters representative of yeast. It was shown that the common approach of assuming reasonable values for unmeasured parameters (e.g. cell-wall thickness, initial radial stretch) can give rise to nonunique solutions for both the form and constants in the cell-wall constitutive relationship. Similarly, measurement errors can also lead to an incorrectly defined cell-wall constitutive relationship. Unique determination of the fundamental wall properties by cell compression requires accurate and precise measurement of a minimum set of parameters (initial cell radius, initial cell-wall thickness, and the volume loss during compression). In the absence of such measurements the derived constitutive relationship may be in considerable error, and should be evaluated against its ability to predict the outcome of other mechanical experiments. (C) 1998 Elsevier Science Ltd. All rights reserved.

Wall material properties of yeast cells. Part II. Analysis

In the preceding paper (Part I) force-deformation data were measured with the compression experiment in conjunction with the initial radial stretch ratio and the initial wall-thickness to cell-radius ratio for baker's yeast (Saccharomyces cerevisiae). In this paper, these data have been analysed with the mechanical model of Smith et al. (Smith, Moxham & Middelberg (1998) Chemical Engineering Science, 53, 3913-3922) with the wall constitutive behaviour defined a priori as incompressible and linear-elastic. This analysis determined the mean Young's modulus ((E) over bar), mean maximum von Mises stress-at-failure (<(sigma)over bar>(VM,f)) and mean maximum von Mises strain-at failure (<(epsilon)over bar>(VM,f)) to be (E) over bar = 150 +/- 15 MPa, <(sigma)over bar>(VM,f) = 70 +/- 4 MPa and <(epsilon)over bar>(VM,f) = 0.75 +/- 0.08, respectively. The mean Young's modulus was not dependent (P greater than or equal to 0.05) on external osmotic pressure (0-0.8 MPa) nor compression rate (1.03-7.68 mu m/s) suggesting the incompressible linear-elastic relationship is representative of the actual cell-wall constitutive behaviour. Hydraulic conductivities were also determined and were comparable to other similar cell types (0-2.5 mu m/MPa s). The hydraulic conductivity distribution was not dependent on external osmotic pressure (0-0.8 MPa) nor compression rate (1.03-7.68 mu m/s) suggesting inclusion of cell-wall permeability in the mechanical model is justified. <(epsilon)over bar>(VM,f) was independent of cell diameter and to a first-approximation unaffected (P greater than or equal to 0.01) by external osmotic pressure and compression rate, thus providing a reasonable failure criterion. This criterion states that the cell-wall material will break when the strain exceeds <(epsilon)over bar>(VM,f) = 0.75 +/- 0.08. Variability in overall cell strength during compression was shown to be primarily due to biological variability in the maximum von Mises strain-at-failure. These data represent the first estimates of cell-wall material properties for yeast and the first fundamental analysis of cell-compression data. They are essential for describing cell-disruption at the fundamental level of fluid-cell interactions in general bioprocesses. They also provide valuable new measurements for yeast-cell physiologists. (C) 2000 Elsevier Science Ltd. All rights reserved.

The role of mobile belts for the longevity of deep cratonic lithosphere: The crumple zone model

Many Archean cratons are surrounded by Proterozoic mobile belts that have experienced episodes of tectonic re-activation over their lifetimes. This suggests that mobile belt lithosphere may be associated with long lived, inherited weakness. It is proposed that the proximity of this weakness can increase the longevity of deep Archean lithosphere by buffering Archean cratons from mantle derived stresses. The physical plausibility of this idea is explored through numerical simulations of mantle convection that include continents and allow for material rheologies that model the combined brittle and ductile behavior of the lithosphere. Within the simulations, the longevity of deep cratonic lithosphere does increase if it is buffered by mobile belts that can fail at relatively low stress levels.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Briefing: Factored material properties and limit state loads-unlikely extreme or impossible pretense

In the limit state design (LSD) method each design criterion is formally stated and assessed using a performance function. The performance function defines the relationship between the design parameters and the design criterion. In practice, LSD involves factoring up loads and factoring down calculated strengths and material parameters. This provides a convenient way to carry out routine probabilistic-based design. The factors are statistically calculated to produce a design with an acceptably low probability of failure. Hence the ultimate load and the design material properties are mathematical concepts that have no physical interpretation. They may be physically impossible. Similarly, the appropriate analysis model is also defined by the performance function and may not describe the real behaviour at the perceived physical equivalent limit condition. These points must be understood to avoid confusion in the discussion and application of partial factor LSD methods.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Kinetic model for selective cultivation of microfungi in a microscreen process for food processing wastewater treatment and biomass production

The research was aimed at developing a technology to combine the production of useful microfungi with the treatment of wastewater from food processing. A recycle bioreactor equipped with a micro-screen was developed as a wastewater treatment system on a laboratory scale to contain a Rhizopus culture and maintain its dominance under non-aseptic conditions. Competitive growth of bacteria was observed, but this was minimised by manipulation of the solids retention time and the hydraulic retention time. Removal of about 90% of the waste organic material (as BOD) from the wastewater was achieved simultaneously. Since essentially all fungi are retained behind the 100 mum aperture screen, the solids retention time could be controlled by the rate of harvesting. The hydraulic retention time was employed to control the bacterial growth as the bacteria were washed through the screen at a short HRT. A steady state model was developed to determine these two parameters. This model predicts the effluent quality. Experimental work is still needed to determine the growth characteristics of the selected fungal species under optimum conditions (pH and temperature).

HPGR model verification and scale-up

This paper summarises test results that were used to validate a model and scale-up procedure of the high pressure grinding roll (HPGR) which was developed at the JKMRC by Morrell et al. [Morrell, Lim, Tondo, David,1996. Modelling the high pressure grinding rolls. In: Mining Technology Conference, pp. 169-176.]. Verification of the model is based on results from four data sets that describe the performance of three industrial scale units fitted with both studded and smooth roll surfaces. The industrial units are currently in operation within the diamond mining industry and are represented by De Beers, BHP Billiton and Rio Tinto. Ore samples from the De Beers and BHP Billiton operations were sent to the JKMRC for ore characterisation and HPGR laboratory-scale tests. Rio Tinto contributed an historical data set of tests completed during a previous research project. The results conclude that the modelling of the HPGR process has matured to a point where the model may be used to evaluate new and to optimise existing comminution circuits. The model prediction of product size distribution is good and has been found to be strongly dependent of the characteristics of the material being tested. The prediction of throughput and corresponding power draw (based on throughput) is sensitive to inconsistent gap/diameter ratios observed between laboratory-scale tests and full-scale operations. (C) 2004 Elsevier Ltd. All rights reserved.

The effect of adsorption, filter material and point of dilution on antibiotic elimination by haemofiltration - An in vitro study of levofloxacin

We studied an in vitro model of continuous venous-venous haemofiltration (CVVH), into which levofloxacin 100 mg was infused, to determine levofloxacin adsorption and to determine the effect of filter material and point of dilution (pre- or post-filter) on sieving coefficient. Mean (standard deviation; S.D.) adsorption was 18.7 (5.3) mg for the polyamide filter and 40.2 (2.0) mg for the polyacrylonitrile (PAN) filter (P < 0.001). Post-dilution resulted in a minor, but statistically significant, decrease in sieving coefficient (pre-dilution 0.96 (S.D. 0.10), post-dilution 0.88 (S.D. 0.11) with the PAN filter. These data indicate that the variability in published values for levofloxacin sieving coefficient are not due to variation in point of dilution or membrane type (PAN or polyamide). Significant adsorption of levofloxacin onto PAN filters occurs. (C) 2004 Elsevier B.V. and the International Society of Chemotherapy. All rights reserved.

The rellationship between specific features of aphasia-friendly written material and comprehension of written material for people with aphasia: An exploratory study

Background: Written material is often inaccessible fro people with aphasia. The format of written material needs to be adapted to enable people with aphasia to read with understanding. Aims: This study aimed to further explore some issues raised in Rose, Worrall, and MacKenna (2003) concerning the effects of aphasia-friendly formats on the reading comprehension of people with aphasia. It was hypothesised that people with aphasia would comprehend significantly more paragraphs that were formatted in an aphasia-friendly manner than control paragraphs. This study also aimed to investigate if each single aspect of aphasia-friendly formatting (i.e., simplified vocabulary and syntax, large print, increased white spacem and pictures) used in isolation would result in increased comprehension compared to control paragraphs. Other aims were to compare the effect of aphasia-friendly fromatting with the effects of each single adaptation, and to investigate if the effects of aphasia-friendly formates were related to aphasia severity. Methods & Procedures: Participants with mild to moderately severe aphasia (N = 9) read a battery of 90 paragraphs and selected the best word of phrase from a choice of four to complete each paragraph. A linear mixed model (p < .05) was used to analyse the differences in reading comprehension with each paragraph fromat across three reading grade levels. Outcomes & Results: People with aphasia comprehended significantly more aphasia-friendly paragraphs than control paragraphs. They also comprehended significantly more paragraphs with each of the following single adaptations: simplified vocabulary and syntax, large ptint, and increased white spaces. Although people with aphasia tended to comprehend more paragraphs with pictures added than control paragraphs, this difference was not significant. No significant correlation between aphasia severity and the effect of aphasia-friendly formatting was found. Conclusion: This study supports the idea that aphasia-friendly formats increase the reading comprehension of people with aphasia. It suggests that adding pictures, particularly Clip Art pictures, may not significantly improve the reading the reading comprehension of people with aphasia. These findings have implications for all written communication with people with aphasia, both in the clinical setting and in the wider community. Applying these findings may enable people with aphasia to have equal access to written information and to participate in society.

Potential to increase yield in lucerne (Medicago sativa subsp sativa) through introgression of Medicago sativa subsp falcata into Australian adapted material

The effect of interspecific heterosis in crosses between Medicago sativa subsp. sativa and M. sativa subsp. falcata was assessed. Three sativa and 3 falcata plants were crossed in a diallel design. Progeny dry matter yield and natural plant height were assessed in a replicated field experiment at Gatton, Queensland. Yield data were analysed using the method of residual maximum likelihood (REML) and Griffing's model 1. There were significant differences between the reciprocal, general combining ability (GCA), and specific combining ability (SCA) effects. As expected, S-1 populations were lower yielding than their respective intraspecific cross and falcata x falcata crosses were significantly lower yielding than sativa x sativa crosses. Some of the interspecific crosses indicated substantial SCA effects, yielding at least as well as the best sativa x sativa crosses. We have demonstrated the potential usefulness of unselected M. sativa subsp. falcata as a heterotic group in the improvement of yield in northern Australian adapted lucerne material, and discuss how it could be incorporated into future breeding to overcome the yield stagnation currently being experienced in Australian programs.