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.

Differential perturbation of neuronal and glial glutamate transport systems in retinal ischaemia

Glutamate is the major excitatory neurotransmitter in the retina and is removed from the extracellular space by an energy-dependent process involving neuronal and glial cell transporters. The radial glial Muller cells express the glutamate transporter, GLAST, and preferentially accumulate glutamate. However, during an ischaemic episode, extracellular glutamate concentrations may rise to excitotoxic levels. Is this catastrophic rise in extracellular glutamate due to a failure of GLAST? Using immunocytochemistry, we monitored the transport of the glutamate transporter substrate, D-aspartate, in the retina under normal and ischaemic conditions. Two models of compromised retinal perfusion were compared: (1) Anaesthetised rats had their carotid arteries occluded for 7 days to produce a chronic reduction in retinal blood flow. Retinal function was assessed by electroretinography. D-aspartate was injected into the eye for 45 min, Following euthanasia, the retina was processed for D-aspartate. GLAST and glutamate immunocytochemistry. Although reduced retinal perfusion suppresses the electroretinogram b-wave, neither retinal histology, GLAST expression, nor the ability of Muller cells to uptake D-aspartate is affected. As this insult does not appear to cause excitotoxic neuronal damage, these data suggest that GLAST function and glutamate clearance are maintained during periods of reduced retinal perfusion. (2) Occlusion of the central retinal artery for 60 min abolishes retinal perfusion, inducing histological damage and electroretinogram suppression. Although GLAST expression appears to be normal. its ability to transport D-aspartate into Muller cells is greatly reduced. Interestingly, D-aspartate is transported into neuronal cells, i.e. photoreceptors, bipolar and ganglion cells. This suggests that while GLAST is vitally important for the clearance of excess extracellular glutamate, its capability to sustain inward transport is particularly susceptible to an acute ischaemic attack. Manipulation of GLAST function could alleviate the degeneration and blindness that result from ischaemic retinal disease. (C) 2001 Elsevier Science Ltd, All rights reserved.

Assessment of regional long-axis function during dobutamine echocardiography

Echocardiographic analysis of regional left ventricular function is based upon the assessment of radial motion. Long-axis motion is an important contributor to overall function. but has been difficult to evaluate clinically until the recent development of tissue Doppler techniques. We sought to compare the standard visual assessment of radial motion with quantitative tissue Doppler measurement of peak systolic velocity. timing and strain rate (SRI) in 104 patients with known or suspected coronary artery disease undergoing dobutamine stress echocardiography (DbE). A standard DbE protocol was used with colour tissue Doppler images acquired in digital cine-loop format. peak systolic velocity (PSV), time to peak velocity (TPV) and SRI were assessed off-line by an independent operator. Wall motion was assessed by an experienced reader. Mean PSV, TPV and SRI values were compared with wall motion and the presence of coronary artery disease by angiography. A further analysis included assessing the extent of jeopardized myocardium by comparing average values of PSV, TPV and SRI against the previously validated angiographic score. Segments identified as having normal and abnormal radial wall motion showed significant differences in mean PSV (7.9 +/- 3.8 and 5.9 +/- 3.3 cm/s respectively; P < 0.001), TPV (84 40 and 95 +/- 48 ms respectively; P = 0.005) and SRI (- 1.45 +/- 0.5 and - 1.1 +/- 0.9 s(-1) respectively; P < 0.001). The presence of a stenosed subtending coronary artery was also associated with significant differences from normally perfused segments for mean PSV (8.1 3.4 compared with 5.7 +/- 3.7 cm/s; P < 0.001), TPV (78 50 compared with 92 +/- 45 ms; P < 0.001) and SRI (- 1.35 0.5 compared with - 1.20 +/- 0.4 s(-1); P = 0.05). PSV, TPV and SRI also varied significantly according to the extent of jeopardized myocardium within a vascular territory. These results suggest that peak systolic velocity, timing of contraction and SRI reflect the underlying physiological characteristics of the regional myocardium during DbE, and may potentially allow objective analysis of wall motion.

Solution techniques for transport problems involving steep concentration gradients: application to noncatalytic fluid solid reactions

Some efficient solution techniques for solving models of noncatalytic gas-solid and fluid-solid reactions are presented. These models include those with non-constant diffusivities for which the formulation reduces to that of a convection-diffusion problem. A singular perturbation problem results for such models in the presence of a large Thiele modulus, for which the classical numerical methods can present difficulties. For the convection-diffusion like case, the time-dependent partial differential equations are transformed by a semi-discrete Petrov-Galerkin finite element method into a system of ordinary differential equations of the initial-value type that can be readily solved. In the presence of a constant diffusivity, in slab geometry the convection-like terms are absent, and the combination of a fitted mesh finite difference method with a predictor-corrector method is used to solve the problem. Both the methods are found to converge, and general reaction rate forms can be treated. These methods are simple and highly efficient for arbitrary particle geometry and parameters, including a large Thiele modulus. (C) 2001 Elsevier Science Ltd. All rights reserved.

A beam-forming network for a circular switched-beam phased array antenna

This paper reports on the design and development of a dividing/phasing network for a compact switched-beam array antenna for Land-vehicle mobile satellite communications, The device is formed by a switched radial divider/combiner and 1-bit phase shifters and generates a sufficient number of beams for the proper satellite tracking.

Hyaline tissue of thermally unaltered conodont elements and the enamel of vertebrates

Comparison of the ultrastructure of the hyaline tissue of conodont elements and the enamel of vertebrates provides little support for a close phylogenetic relationship between conodonts and vertebrates. Transmission and scanning electron microscopy shows that the mineralised component of the hyaline tissue of Panderodus and of Cordylodus elements consists of large, flat, oblong crystals, arranged in layers that run parallel to the long axis of the conodont. Enamel in the dentition of a living vertebrate, the lungfish Neoceratodus forsteri, has crystals of calcium hydroxyapatite, arranged in layers, and extending in groups from the dentine-enamel junction; the crystals are slender, elongate spicules perpendicular to the surface of the tooth plate, Similar crystal arrangements to those of lungfish are found in other vertebrates, but none resembles the organisation of the hyaline tissue of conodont elements, The crystals of hydroxyapatite in conodont hyaline tissue are exceptionally large, perpendicular or parallel to the surface of the element, with no trace of prisms, unlike the protoprismatic radial crystallite enamel of fish teeth and scales, or the highly organised prismatic enamel of mammals.

Computation of the linear elastic properties of random porous materials with a wide variety of microstructure

A finite-element method is used to study the elastic properties of random three-dimensional porous materials with highly interconnected pores. We show that Young's modulus, E, is practically independent of Poisson's ratio of the solid phase, nu(s), over the entire solid fraction range, and Poisson's ratio, nu, becomes independent of nu(s) as the percolation threshold is approached. We represent this behaviour of nu in a flow diagram. This interesting but approximate behaviour is very similar to the exactly known behaviour in two-dimensional porous materials. In addition, the behaviour of nu versus nu(s) appears to imply that information in the dilute porosity limit can affect behaviour in the percolation threshold limit. We summarize the finite-element results in terms of simple structure-property relations, instead of tables of data, to make it easier to apply the computational results. Without using accurate numerical computations, one is limited to various effective medium theories and rigorous approximations like bounds and expansions. The accuracy of these equations is unknown for general porous media. To verify a particular theory it is important to check that it predicts both isotropic elastic moduli, i.e. prediction of Young's modulus alone is necessary but not sufficient. The subtleties of Poisson's ratio behaviour actually provide a very effective method for showing differences between the theories and demonstrating their ranges of validity. We find that for moderate- to high-porosity materials, none of the analytical theories is accurate and, at present, numerical techniques must be relied upon.

Viscoelasticity of radiation-formed PVA/PVP hydrogel

The dynamic theological behaviour of gamma-irradiated 12.8 wt% poly(vinyl alcohol) (PVA), 12.8 wt% poly(vinyl pyrrolidone) (PVP), and a blend of 8 wt% PVA and 4.8 wt% PVP aqueous solutions have been studied pre- and post-gelation. The non-irradiated solutions displayed theological behaviour typical of dilute to semi-dilute polymer solutions, with the complex viscosity being independent of the frequency and shear rate (i.e. Newtonian behaviour) over the range of frequencies tested and the loss modulus G(omega) and storage modulus G(omega) being nearly proportional to omega and omega(2) respectively. After a set of doses of gamma-radiation, the magnitudes of the dynamic moduli G'(omega) and G(omega) increased as the absorbed dose increased, with notable differences between the two homopolymers and the blend. The stages of gelation were effectively monitored by means of dynamic theological measurements, allowing the possible mechanisms of network formation to be elucidated. The doses required for gelation of the PVA, PVP, and blend samples, determined on the basis of the Winter and Chambon criteria for gelation, were found to be 12 kGy for the 12.8 wt% PVA, 4 kGy for the 12.8 wt% PVP, and 5 kGy for the 8 wt% PVA/4.8 wt% PVP solutions. The unexpected lower gelation dose demonstrated by the blend sample, compared with predictions based on the blend composition, and the associated gelation mechanism are also discussed.

Drop penetration into porous powder beds

The kinetics of drop penetration were studied by filming single drops of several different fluids (water, PEG200, PEG600, and HPC solutions) as they penetrated into loosely packed beds of glass ballotini, lactose, zinc oxide, and titanium dioxide powders. Measured times ranged from 0.45 to 126 s and depended on the powder particle size,viscosity, surface tensions, and contact angle. The experimental drop penetration times were compared to existing theoretical predictions by M. Denesuk et al. (J. Colloid Interface Sci. 158, 114, 1993) and S. Middleman (Modeling Axisymmetric Flows: Dynamics of Films, Jets, and Drops, Academic Press, San Diego, 1995) but did not agree. Loosely packed powder beds tend to have a heterogeneous bed structure containing large macrovoids which do not participate in liquid flow but are included implicitly in the existing approach to estimating powder pore size. A new two-phase model was proposed where the total volume of the macrovoids was assumed to be the difference between the bed porosity and the tap porosity. A new parameter, the effective porosity (epsilon)eff, was defined as the tap porosity multiplied by the fraction of pores that terminate at a macrovoid and are effectively blocked pores. The improved drop penetration model was much more successful at estimating the drop penetration time on all powders and the predicted times were generally within an order of magnitude of the experimental results. (C) 2002 Elsevier Science (USA).

Investigation of ultrasonic properties of PAG and MAGIC polymer gel dosimeters

Ultrasonic speed of propagation and attenuation were investigated as a function of absorbed radiation dose in PAG and MAGIC polymer gel dosimeters. Both PAG and MAGIC gel dosimeters displayed a dependence of ultrasonic parameters on absorbed dose with attenuation displaying significant changes in the dose range investigated. The ultrasonic attenuation dose sensitivity at 4 MHz in MAGIC gels was determined to be 4.7 +/- 0.3 dB m(-1) Gy(-1) and for PAG 3.9 +/- 0.3 dB m(-1) Gy(-1). Ultrasonic speed dose sensitivities were 0.178 +/- 0.006 m s(-1) Gy(-1) for MAGIC gel and -0.44 +/- 0.02 m s(-1) Gy(-1) for PAG. Density and compressional elastic modulus were investigated to explain the different sensitivities of ultrasonic speed to radiation for PAG and MAGIC gels. The different sensitivities were found to be due to differences in the compressional elastic modulus as a function of dose for the two formulations. To understand the physical phenomena underlying the increase in ultrasonic attenuation with dose, the viscoelastic properties of the gels were studied. Results suggest that at ultrasonic frequencies, attenuation in polymer gel dosimeters is primarily due to volume viscosity. It is concluded that ultrasonic attenuation significantly increases with absorbed dose. Also, the ultrasonic speed in polymer gel dosimeters is affected by changes in dosimeter elastic modulus that are likely to be a result of polymerization. It is suggested that ultrasound is a sufficiently sensitive technique for polymer gel dosimetry.

Application of the Williams-Landel-Ferry model to the viscosity-temperature relationship of Australian honeys

The rheological behaviour of nine unprocessed Australian honeys was investigated for the applicability of the Williams-Landel-Ferry (WLF) model. The viscosity of the honeys was obtained over a range of shear rates (0.01-40 s(-1)) from 2degrees to 40 degreesC, and all the honeys exhibited Newtonian behaviour with viscosity reducing as the temperature was increased. The honeys with high moisture were of lower viscosity, The glass transition temperatures of the honeys, as measured with a differential scanning calorimeter (DSC), ranged from -40degrees to -46 degreesC, and four models (WLF. Arrhenius, Vogel-Tammann-Fulcher (VTF), and power-law) were investigated to describe the temperature dependence of the viscosity. The WLF was the most suitable and the correlation coefficient averaged 0.999 +/- 0.0013 as against 0.996 +/- 0.0042 for the Arrhenius model while the mean relative deviation modulus was 0-12% for the WLF model and 10-40% for the Arrhenius one. With the universal values for the WLF constants, the temperature dependence of the viscosity was badly predicted. From non-linear regression analysis, the constants of the WLF models for the honeys were obtained (C-1 = 13.7-21.1: C-2 = 55.9-118.7) and are different from the universal values. These WLF constants will be valuable for adequate modeling of the rheology of the honeys, and they can be used to assess the temperature sensitivity of the honeys. (C) 2002 Elsevier Science Ltd. All rights reserved.

Testing landscape metrics as indicators of habitat loss and fragmentation in continuous eucalypt forests (Queensland, Australia)

Landscape metrics are widely applied in landscape ecology to quantify landscape structure. However, many are poorly tested and require rigorous validation if they are to serve as reliable indicators of habitat loss and fragmentation, such as Montreal Process Indicator 1.1e. We apply a landscape ecology theory, supported by exploratory and confirmatory statistical techniques, to empirically test landscape metrics for reporting Montreal Process Indicator 1.1e in continuous dry eucalypt forests of sub-tropical Queensland, Australia. Target biota examined included: the Yellow-bellied Glider (Petaurus australis); the diversity of nectar and sap feeding glider species including P. australis, the Sugar Glider P. breviceps, the Squirrel Glider P. norfolcensis, and the Feathertail Glider Acrobates pygmaeus; six diurnal forest birds species; total diurnal bird species diversity; and the density of nectar-feeding diurnal bird species. Two scales of influence were considered: the stand-scale (2 ha), and a series of radial landscape extents (500 m - 2 km; 78 - 1250 ha) surrounding each fauna transect. For all biota, stand-scale structural and compositional attributes were found to be more influential than landscape metrics. For the Yellow-bellied Glider, the proportion of trace habitats with a residual element of old spotted-gum/ironbark eucalypt trees was a significant landscape metric at the 2 km landscape extent. This is a measure of habitat loss rather than habitat fragmentation. For the diversity of nectar and sap feeding glider species, the proportion of trace habitats with a high coefficient of variation in patch size at the 750 m extent was a significant landscape metric. None of the landscape metrics tested was important for diurnal forest birds. We conclude that no single landscape metric adequately captures the response of the region's forest biota per se. This poses a major challenge to regional reporting of Montreal Process Indicator 1.1e, fragmentation of forest types.

Direction selectivity in the retina

The neuronal circuitry underlying the generation of direction selectivity in the retina has remained elusive for almost 40 years. Recent studies indicate that direction selectivity may be established within the radial dendrites of 'starburst' amacrine cells and that retinal ganglion cells may acquire their direction selectivity by the appropriate weighting of excitatory and inhibitory inputs from starburst dendrites pointing in different directions. If so, this would require unexpected complexity and subtlety in the synaptic connectivity of these CNS neurons.

Leaf water relations of Eucalyptus cloeziana and Eucalyptus argophloia in response to water deficit

Leaf water relations responses to limited water supply were determined in 7-month-old plants of a dry inland provenance of Eucalyptus argophloia Blakely and in a humid coastal provenance (Gympie) and a dry inland provenance (Hungry Hills) of Eucalyptus cloeziana F. Muell. Each provenance of E. cloeziana exhibited a lower relative water content at the turgor loss point, a lower apoplastic water content, a smaller ratio of dry mass to turgid mass and a lower bulk modulus of elasticity than the single provenance of E. argophloia. Osmotic potential at full turgor and water potential at the turgor loss point were significantly lower in E. argophloia and the inland provenance of E. cloeziana than in the coastal provenance of E. cloeziana. There was limited osmotic adjustment in response to soil drying in E. cloeziana, but not in E. argophloia. Between-species differences in water relations parameters were larger than those between the E. cloeziana provenances. Both E. cloeziana provenances maintained turgor under moderate water stress through a combination of osmotic and elastic adjustments. Eucalyptus argophloia had more rigid cell walls and reached lower water potentials with less reduction in relative water content than either of the E. cloeziana provenances, thereby enabling it to extract water from dryer soils.

Modelling the size of the crushed zone around a blasthole

A new model to predict the extent of crushing around a blasthole is presented. The model is based on the back-analysis of a comprehensive experimental program that included the direct measurement of the zone of crushing from 92 blasting tests on concrete blocks using two commercial explosives. The concrete blocks varied from low, medium to high strength and measured 1.5 in in length, 1.0 m in width and 1.1 m in height. A dimensionless parameter called the crushing zone index (CZI) is introduced. This index measures the crushing potential of a charged blasthole and is a function of the borehole pressure, the unconfined compressive strength of the rock material, dynamic Young's modulus and Poisson's ratio. It is shown that the radius of crushing is a function of the CZI and the blasthole radius. A good correlation between the new model and measured results was obtained. A number of previously proposed models could not approximate the conditions measured in the experimental work and there are noted discrepancies between the different approaches reviewed, particularly for smaller diameter holes and low strength rock conditions. The new model has been verified with full scale tests reported in the literature. Results from this validation and model evaluations show its applicability to production blasting. (C) 2003 Elsevier Science Ltd. All rights reserved.