The crack tip strain field of AISI 4340 - Part I - Measurement technique

This is the first paper in a study on the influence of the environment on the crack tip strain field for AISI 4340. A stressing stage for the environmental scanning electron microscope (ESEM) was constructed which was capable of applying loads up to 60 kN to fracture-mechanics samples. The measurement of the crack tip strain field required preparation (by electron lithography or chemical etching) of a system of reference points spaced at similar to 5 mu m intervals on the sample surface, loading the sample inside an electron microscope, image processing procedures to measure the displacement at each reference point and calculation of the strain field. Two algorithms to calculate strain were evaluated. Possible sources of errors were calculation errors due to the algorithm, errors inherent in the image processing procedure and errors due to the limited precision of the displacement measurements. Estimation of the contribution of each source of error was performed. The technique allows measurement of the crack tip strain field over an area of 50 x 40 mu m with a strain precision better than +/- 0.02 at distances larger than 5 mu m from the crack tip. (C) 1999 Kluwer Academic Publishers.

The crack tip strain field of AISI 4340 - Part II - Experimental measurements

Crack tip strain maps have been measured for AISI 4340 high strength steel. No significant creep was observed. The measured values of CTOD were greater than expected from the HRR model. Crack tip branching was observed in every experiment. The direction of crack branching was in the same direction as a major ridge'' of epsilon(yy) strain, which in turn was in the same direction as predicted by the HRR model. Furthermore, the measured magnitudes of the epsilon(y)y strain in this same direction were in general greater than the values predicted by the HRR model. This indicates more plasticity in the crack tip region than expected from the HRR model. This greater plasticity could be related to the larger than expected CTOD values. The following discrepancies between the measured strain fields for AISI 4340 and the HRR predictions are noteworthy: (1) The crack branching. (2) Values of CTOD significantly higher than predicted by HRR. (3) The major ridge'' of epsilon(yy) strain an angle of about 60 degrees with the direction of overall propagation of the fatigue precrack, in which the measured magnitudes of the epsilon(yy) strain were greater than the values predicted by the HRR model. (4) Asymmetric shape of the plastic zone as measured by the epsilon(yy) strain. (5) Values of shear strain gamma(xy) significantly higher than predicted by the HRR model. (C) 1999 Kluwer Academic Publishers.

The crack tip strain field of AISI 4340 - Part III - Hydrogen influence

This paper studied the influence of hydrogen and water vapour environments on the plastic behaviour in the vicinity of the crack tip for AISI 4340. Hydrogen and water vapour (at a pressure of 15 Torr) significantly increased the crack tip opening displacement. The crack tip strain distribution in 15 Torr hydrogen was significantly different to that measured in vacuum. In the presence of sufficient hydrogen, the plastic zone was larger, was elongated in the direction of crack propagation and moreover there was significant creep. These observations support the hydrogen enhanced localised plasticity model for hydrogen embrittlement in this steel. The strain distribution in the presence of water vapour also suggests that SCC in AISI 4340 occurs via the hydrogen enhanced localised plasticity mechanism. (C) 1999 Kluwer Academic Publishers.

NMR studies of the low-density lipoprotein receptor-binding peptide of apolipoprotein E bound to dodecylphosphocholine micelles

Circular dichroism and NMR spectroscopy have been used to determine the structure of the low-density lipoprotein (LDL) receptor-binding peptide, comprising residues 130-152, of the human apolipoprotein E. This peptide has little persistent three-dimensional structure in solution, but when bound to micelles of dodecylphosphocholine (DPC) it adopts a predominantly alpha-helical structure. The three-dimensional structure of the DPC-bound peptide has been determined by using H-1-NMR spectroscopy: the structure derived from NOE-based distance constraints and restrained molecular dynamics is largely helical. The derived phi and psi angle order parameters show that the helical structure is well defined but with some flexibility that causes the structures not to be superimposable over the full peptide length. Deuterium exchange experiments suggest that many peptide amide groups are readily accessible to the solvent, but those associated with hydrophobic residues exchange more slowly, and this helix is thus likely to be positioned on the surface of the DPC micelles. In this conformation the peptide has one hydrophobic face and two that are rich in basic amino acid side chains. The solvent-exposed face of the peptide contains residues previously shown to be involved in binding to the LDL receptor.

A hybrid, inverse approach to the design of magnetic resonance imaging magnets

This paper describes a hybrid numerical method of an inverse approach to the design of compact magnetic resonance imaging magnets. The problem is formulated as a field synthesis and the desired current density on the surface of a cylinder is first calculated by solving a Fredholm equation of the first, kind. Nonlinear optimization methods are then invoked to fit practical magnet coils to the desired current density. The field calculations are performed using a semi-analytical method. The emphasis of this work is on the optimal design of short MRI magnets. Details of the hybrid numerical model are presented, and the model is used to investigate compact, symmetric MRI magnets as well as asymmetric magnets. The results highlight that the method can be used to obtain a compact MRI magnet structure and a very homogeneous magnetic field over the central imaging volume in clinical systems of approximately 1 m in length, significantly shorter than current designs. Viable asymmetric magnet designs, in which the edge of the homogeneous region is very close to one end of the magnet system are also presented. Unshielded designs are the focus of this work. This method is flexible and may be applied to magnets of other geometries. (C) 2000 American Association of Physicists in Medicine. [S0094-2405(00)00303-5].

Effect of density on spatial distribution, pit formation and pit diameter of Myrmeleon acer Walker, (Neuroptera : Myrmeleontidae): patterns and processes

The effect of increasing population density on the formation of pits, their size and spatial distribution, and on levels of mortality was examined in the antlion Myrmeleon acer Walker. Antlions were kept at densities ranging from 0.4 to 12.8 individuals per 100 cm(2). The distribution of pits was regular or uniform across all densities, but antlions constructed proportionally fewer and smaller pits as density increased. Mortality through cannibalism was very low and only occurred at densities greater than five individuals per 100 cm(2). Antlions in artificially crowded situations frequently relocated their pits and when more space became available, individuals became more dispersed with time. Redistribution of this species results from active avoidance of other antlions and sand throwing associated with pit construction and maintenance, rather than any attempt to optimise prey capture per se.

Comparative analysis of dendritic cell density and total number in commonly transplanted organs: morphometric estimation in normal mice

Dendritic cells (DC) are considered to be the major cell type responsible for induction of primary immune responses. While they have been shown to play a critical role in eliciting allosensitization via the direct pathway, there is evidence that maturational and/or activational heterogeneity between DC in different donor organs may be crucial to allograft outcome. Despite such an important perceived role for DC, no accurate estimates of their number in commonly transplanted organs have been reported. Therefore, leukocytes and DC were visualized and enumerated in cryostat sections of normal mouse (C57BL/10, B10.BR, C3H) liver, heart, kidney and pancreas by immunohistochemistry (CD45 and MHC class II staining, respectively). Total immunopositive cell number and MHC class II+ cell density (C57BL/10 mice only) were estimated using established morphometric techniques - the fractionator and disector principles, respectively. Liver contained considerably more leukocytes (similar to 5-20 x 10(6)) and DC (similar to 1-3 x 10(6)) than the other organs examined (pancreas: similar to 0.6 x 10(6) and similar to 0.35 x 10(6): heart: similar to 0.8 x 10(6) and similar to 0.4 x 10(6); kidney similar to 1.2 x 10(6) and 0.65 x 10(6), respectively). In liver, DC comprised a lower proportion of all leukocytes (similar to 15-25%) than in the other parenchymal organs examined (similar to 40-60%). Comparatively, DC density in C57BL/10 mice was heart > kidney > pancreas much greater than liver (similar to 6.6 x 10(6), 5 x 10(6), 4.5 x 10(6) and 1.1 x 10(6) cells/cm(3), respectively). When compared to previously published data on allograft survival, the results indicate that the absolute number of MHC class II+ DC present in a donor organ is a poor predictor of graft outcome. Survival of solid organ allografts is more closely related to the density of the donor DC network within the graft. (C) 2000 Elsevier Science B.V. All rights reserved.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

Evolution of stress deficit and changing rates of seismicity in cellular automaton models of earthquake faults

We investigate the internal dynamics of two cellular automaton models with heterogeneous strength fields and differing nearest neighbour laws. One model is a crack-like automaton, transferring ail stress from a rupture zone to the surroundings. The other automaton is a partial stress drop automaton, transferring only a fraction of the stress within a rupture zone to the surroundings. To study evolution of stress, the mean spectral density. f(k(r)) of a stress deficit held is: examined prior to, and immediately following ruptures in both models. Both models display a power-law relationship between f(k(r)) and spatial wavenumber (k(r)) of the form f(k(r)) similar tok(r)(-beta). In the crack model, the evolution of stress deficit is consistent with cyclic approach to, and retreat from a critical state in which large events occur. The approach to criticality is driven by tectonic loading. Short-range stress transfer in the model does not affect the approach to criticality of broad regions in the model. The evolution of stress deficit in the partial stress drop model is consistent with small fluctuations about a mean state of high stress, behaviour indicative of a self-organised critical system. Despite statistics similar to natural earthquakes these simplified models lack a physical basis. physically motivated models of earthquakes also display dynamical complexity similar to that of a critical point system. Studies of dynamical complexity in physical models of earthquakes may lead to advancement towards a physical theory for earthquakes.

Race and sex effects on the association between muscle strength, soft tissue, and bone mineral density in healthy elders: The Health, Aging, and Body Composition Study

Two factors generally reported to influence bone density are body composition and muscle strength. However, it is unclear if these relationships are consistent across race and sex, especially in older persons. If differences do exist by race and/or sex, then strategies to maintain bone mass or minimize bone loss in older adults may need to be modified accordingly. Therefore, we examined the independent effects of bone mineral-free lean mass (LM), fat mass (FM), and muscle strength on regional and whole body bone mineral density (BMD) in a cohort of 2619 well-functioning older adults participating in the Health, Aging, and Body Composition (Health ABC) Study with complete measures. Participants included 738 white women, 599 black women, 827 white men, and 455 black men aged 70-79 years. BMD (g/cm(2)) of the femoral neck, whole body, upper and lower limb, and whole body and upper limb bone mineral-free LM and FM was assessed by dual-energy X-ray absorptiometry (DXA). Handgrip strength and knee extensor torque were determined by dynamometry. In analyses stratified by race and sex and adjusted for a number of confounders, LM was a significant (p < 0.001) determinant of BMD, except in white women for the lower limb and whole body. In women, FM also was an independent contributor to BMD at the femoral neck, and both PM and muscle strength contributed to limb BMD. The following were the respective Beta-weights (regression coefficients for standardized data, Std beta) and percent difference in BMD per unit (7.5 kg) LM: femoral neck, 0.202-0.386 and 4.7-6.9 %; lower limb,.0.209-0.357 and 2.9-3.5%; whole body, 0.239-0.484 and 3.0-4.7 %; and upper limb (unit = 0.5 kg), 0.231-0.407 and 3.1-3.4%. Adjusting for bone size (bone mineral apparent density [BMAD]) or body size BMD/height) diminished the importance of LM, and the contributory effect of FM became more pronounced. These results indicate that LM and FM were associated with bone mineral depending on the bone site and bone index used. Where differences did occur, they were primarily by sex not race. To preserve BMD, maintaining or increasing LM in the elderly would appear to be an appropriate strategy, regardless of race or sex.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Mantle convection modeling with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling

The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.

Aspects of the biology of the yabby Trypea australiensis (Dana) (Decapoda : Thalassinidea) and the potential of burrow counts as an indirect measure of population density

The thalassinidean shrimp Trypea australiensis (the yabby) commonly occurs on intertidal sandflats and subtidal regions of sheltered embayments and estuaries along the east coast of Australia and is harvested commercially and recreationally for use as bait by anglers. The potential for counts of burrow openings to provide a reliable indirect estimate of the abundance of yabbies was examined on intertidal sandflats on North Stradbroke Island (Queensland, Australia). The relationship between the number of burrow openings and the abundance of yabbies was generally poor and also varied significantly through time, casting doubt on previous estimates of abundance for this species based on unvalidated hole counts. Spatial and temporal variation in population density, the size at maturity and the reproductive period of the yabby were also assessed. Except for an initial peak in abundance as a result of recruitment, the density of yabbies was constant throughout the study but considerably less than that estimated from a previous study in the same area. Ovigerous females were recorded at 3 mm carapace length (CL) which is smaller than previously recorded for this species and thalassinideans in general. Small ovigerous females were found throughout the study, including the summer months, which is unusual for thalassinideans in the intertidal zone. It was hypothesised that in the intertidal zone, small female yabbies may be able to balance the metabolic demands of reproduction and respiration at higher temperatures than can larger females allowing them to reproduce in the warmer months.