A study of the micro-mechanics of granular media

The development of more realistic constitutive models for granular media, such as sand, requires ingredients which take into account the internal micro-mechanical response to deformation. Unfortunately, at present, very little is known about these mechanisms and therefore it is instructive to find out more about the internal nature of granular samples by conducting suitable tests. In contrast to physical testing the method of investigation used in this study employs the Distinct Element Method. This is a computer based, iterative, time-dependent technique that allows the deformation of granular assemblies to be numerically simulated. By making assumptions regarding contact stiffnesses each individual contact force can be measured and by resolution particle centroid forces can be calculated. Then by dividing particle forces by their respective mass, particle centroid velocities and displacements are obtained by numerical integration. The Distinct Element Method is incorporated into a computer program 'Ball'. This program is effectively a numerical apparatus which forms a logical housing for this method and allows data input and output, and also provides testing control. By using this numerical apparatus tests have been carried out on disc assemblies and many new interesting observations regarding the micromechanical behaviour are revealed. In order to relate the observed microscopic mechanisms of deformation to the flow of the granular system two separate approaches have been used. Firstly a constitutive model has been developed which describes the yield function, flow rule and translation rule for regular assemblies of spheres and discs when subjected to coaxial deformation. Secondly statistical analyses have been carried out using data which was extracted from the simulation tests. These analyses define and quantify granular structure and then show how the force and velocity distributions use the structure to produce the corresponding stress and strain-rate tensors.

The mechanical properties and formability of dual-phase steel

In recent years dual phase steels comprising of 5-20% martensite in a ferrite matrix have come into the limelight of high strength cold formable steels because of their potential for vehicle weight saving. They show the following features: no yield point; relatively low initial flow stress; high initial workhardening rate; well sustained work hardening. As a consequence of these characteristics, dual phase steels exhibit a better combination of strength and elongation than other HSLA steels. In this thesis, a broad view of the factors which influence their properties is presented. Mechanical properties and forming ability of a commercially available dual phase steel and an AL-Si killed steel processed to dual phase form are investigated to ascertain the effect of their microstructure on their properties. It is found that the yield phenomena are masked by the transformation induced stresses present during processing and so yield point could be recovered under suitable ageing treatment; that apart from giving the above properties dual phasing gives rise to very low strain-rate sensitivity and a low R value ~ 1; that the mechanical response under rolling conditions is not different from those under tension; that there is a danger of damage to tooling during forming operations of these steels if fracture should precede instability as a result of grain size dependent strength found for these steels. It is also found that very little deformation of the martensite islands took place during deformation except at high strains. The work-hardening and the strength levels can be controlled by either decreasing the grain size or increasing the martensite volume fraction, but it is found that increasing martensite has a detrimental effect on ductility and the ductility and fracture strength can be controlled better by refining the grain size. A remarkable effect found in the dual phase steel tested is that the compressive strength is higher than the tensile strength. The reason for this observation is not yet clear but it is suggested that it might be due to the introduction of emissary type dislocations into the ferrite lattice as a result of twins formed in the martensite during transformation from austenite. The twins are envisaged to be {111} <112> in character.

Effects of dynamic strain ageing on fatigue crack propagation and threshold behaviour in a Ni-base superalloy

In the temperature range 200-400 degree C the Ni-base superalloy, N901, develops marked dynamic strain ageing effects in its tensile behavior. These include inverse strain rate sensitivity, especially in UTS values, strongly serrated stress-strain curves and a heavily sheared failure mode at the higher test-temperatures. As for steels these effects seem to be due to interactions between the dislocations and the interstitial carbon atoms present. The results of tensile and fatigue threshold tests carried out between 20 degree C and 420 degree C are reported and the fatigue behavior is discussed in terms of the effects of surface roughness induced closure, temperature and strain aging interactions.

Microstructure characterization of 316L deformed at high strain rates using EBSD

Specimens from split Hopkinson pressure bar experiments, at strain rates between ~ 1000–9000 s− 1 at room temperature and 500 °C, have been studied using electron backscatter diffraction. No significant differences in the microstructures were observed at different strain rates, but were observed for different strains and temperatures. Size distribution for subgrains with boundary misorientations > 2° can be described as a bimodal lognormal area distribution. The distributions were found to change due to deformation. Part of the distribution describing the large subgrains decreased while the distribution for the small subgrains increased. This is in accordance with deformation being heterogeneous and successively spreading into the undeformed part of individual grains. The variation of the average size for the small subgrain distribution varies with strain but not with strain rate in the tested interval. The mean free distance for dislocation slip, interpreted here as the average size of the distribution of small subgrains, displays a variation with plastic strain which is in accordance with the different stages in the stress-strain curves. The rate of deformation hardening in the linear hardening range is accurately calculated using the variation of the small subgrain size with strain.

Impact on HDPE and PVC plates - Experimental tests and numerical simulations

This paper presents first material tests on HDPE and PVC, and subsequently impact tests on plates made of the same materials. Finally, numerical simulations of the plate impact tests are compared with the experimental results. A rather comprehensive series of mechanical material tests were performed to disclose the behaviour of PVC and HDPE in tension and compression. Quasi-static tests were carried out at three rates in compression and two in tension. Digital image correlation. DIC, was used to measure the in-plane strains, revealing true stress-strain curves and allowing to analyze strain-rate sensitivity and isotropy of Poisson`s ratio. In addition, dynamic compression tests were carried out in a split-Hopkinson pressure bar. Quasi-static and dynamic tests were also performed on clamped plates made of the same PVC and HDPE materials, using an optical technique to measure the full-field out-of-plane deformations. These tests, together with the material data, were used for comparative purposes of a finite element analysis. A reasonable agreement between experimental and numerical results was achieved. (C) 2010 Elsevier Ltd. All rights reserved.

Experiments on stress-triaxiality dependence of material behavior of aluminum alloys

The paper presents and discusses experimental procedures, visual observations and test results considered important to obtain data that can be used in validation of constitutive relations and failure criteria. The aim is to investigate the combined effects of stress intensity, stress-triaxiality and Lode parameter on the material response and failure behavior of aluminum alloys. Smooth and pre-notched tensile and shear specimens were manufactured from both very thin sheets and thicker plates to cover a wide range of stress triaxialities and Lode parameters. In addition, modified Arcan specimens were designed allowing investigation of the effect of sudden changes in stress states and deformation modes on the material behavior. (C) 2009 Elsevier Ltd. All rights reserved.

Scaling of structures subject to impact loads when using a power law constitutive equation

It is well known that structures subjected to dynamic loads do not follow the usual similarity laws when the material is strain rate sensitive. As a consequence, it is not possible to use a scaled model to predict the prototype behaviour. In the present study, this problem is overcome by changing the impact velocity so that the model behaves exactly as the prototype. This exact solution is generated thanks to the use of an exponential constitutive law to infer the dynamic flow stress. Furthermore, it is shown that the adopted procedure does not rely on any previous knowledge of the structure response. Three analytical models are used to analyze the performance of the technique. It is shown that perfect similarity is achieved, regardless of the magnitude of the scaling factor. For the class of material used, the solution outlined has long been sought, inasmuch as it allows perfect similarity for strain rate sensitive structures subject to impact loads. (C) 2009 Elsevier Ltd. All rights reserved.

A numerical investigation on the visco-plastic response of structures to different pulse loading shapes

The time varying intensity character of a load applied to a structure poses many difficulties in analysis. A remedy to this situation is to substitute a complex pulse shape by a rectangular equivalent one. It has been shown by others that this procedure works well for perfectly plastic elementary structures. This paper applies the concept of equivalent pulse to more complex structures. Special attention is given to the material behavior, which is allowed to be strain rate and strain hardening sensitive. Thanks to the explicit finite element solution, it is shown in this article that blast loads applied to complex structures made of real materials can be substituted by equivalent rectangular loads with both responses being practically the same. (c) 2007 Elsevier Ltd. All rights reserved.

Linearly-increasing-stress testing of carbon steel in 4 N NaNO3 and in Bayer liquor

This paper reports the application of linearly increasing stress testing (LIST) to the study of stress corrosion cracking (SCC) of carbon steel in 4 N NaNO3 and in Bayer liquor. LIST is similar to the constant extension-rate testing (CERT) methodology with the essential difference that the LIST is load controlled whereas the CERT is displacement controlled. The main conclusion is that LIST is suitable for the study of the SCC of carbon steels in 4 N NaNO3 and in Bayer liquor. The low crack velocity in Bayer liquor and a measured maximum stress close to that of the reference specimen in air both indicate that a low applied stress rate is required to study SCC in this system. (C) 1998 Chapman & Hall.

On the mechanism of microbial cell disruption in high-pressure homogenisation

The tensions produced in the wall of a rigid, thin-walled, liquid-filled sphere as it moves with an axisymmetric straining flow are examined. This problem has not been previously addressed. A generalised correlation for the maximum wall tension, expressed in dimensionless form as a Weber number (We), is developed in terms of the acceleration number (Ac) and Reynolds number (Re) of the straining flow. At low Reynolds number We is dominated by viscous forces, while inertial forces due to internal pressure gradients caused by sphere acceleration dominate at higher Re. The generalised correlation has been used to examine the case of a typical yeast cell (a thin-walled, liquid-filled sphere) passing through a typical high-pressure homogeniser (a straining-flow device). At 56 MPa homogenising pressure, a 6 mu m yeast cell experiences tensions in the inertially dominated regime (Re = 100). The correlation gives We = 0.206, corresponding to a maximum wall tension of 8 Nm(-1). This is equivalent to an applied compressive force of 150 mu N and compares favourably with the force required to break yeast cells under compressive micromanipulation (40-90 mu N). Inertial forces may therefore be an important and previously unrecognised. mechanism of microbial cell disruption during high-pressure homogenisation. Further work is required to examine the likelihood of cell deformation in the high-strain-rate short-residence-time environment of the homogeniser, and the effect that such deformation may have on the contribution of inertial forces to disruption. (C) 1998 Published by Elsevier Science Ltd. All rights reserved.

Evolution of three-dimensional folds for a non-Newtonian plate in a viscous medium

We derive analytical solutions for the three-dimensional time-dependent buckling of a non-Newtonian viscous plate in a less viscous medium. For the plate we assume a power-law rheology. The principal, axes of the stretching D-ij in the homogeneously deformed ground state are parallel and orthogonal to the bounding surfaces of the plate in the flat state. In the model formulation the action of the less viscous medium is replaced by equivalent reaction forces. The reaction forces are assumed to be parallel to the normal vector of the deformed plate surfaces. As a consequence, the buckling process is driven by the differences between the in-plane stresses and out of plane stress, and not by the in-plane stresses alone as assumed in previous models. The governing differential equation is essentially an orthotropic plate equation for rate dependent material, under biaxial pre-stress, supported by a viscous medium. The differential problem is solved by means of Fourier transformation and largest growth coefficients and corresponding wavenumbers are evaluated. We discuss in detail fold evolutions for isotropic in-plane stretching (D-11 = D-22), uniaxial plane straining (D-22 = 0) and in-plane flattening (D-11 = -2D(22)). Three-dimensional plots illustrate the stages of fold evolution for random initial perturbations or initial embryonic folds with axes non-parallel to the maximum compression axis. For all situations, one dominant set of folds develops normal to D-11, although the dominant wavelength differs from the Biot dominant wavelength except when the plate has a purely Newtonian viscosity. However, in the direction parallel to D-22, there exist infinitely many modes in the vicinity of the dominant wavelength which grow only marginally slower than the one corresponding to the dominant wavelength. This means that, except for very special initial conditions, the appearance of a three-dimensional fold will always be governed by at least two wavelengths. The wavelength in the direction parallel to D-11 is the dominant wavelength, and the wavelength(s) in the direction parallel to D-22 is determined essentially by the statistics of the initial state. A comparable sensitivity to the initial geometry does not exist in the classic two-dimensional folding models. In conformity with tradition we have applied Kirchhoff's hypothesis to constrain the cross-sectional rotations of the plate. We investigate the validity of this hypothesis within the framework of Reissner's plate theory. We also include a discussion of the effects of adding elasticity into the constitutive relations and show that there exist critical ratios of the relaxation times of the plate and the embedding medium for which two dominant wavelengths develop, one at ca. 2.5 of the classical Biot dominant wavelength and the other at ca. 0.45 of this wavelength. We propose that herein lies the origin of parasitic folds well known in natural examples.

Toxicity of biorational insecticides to Helicoverpa spp. (Lepidoptera : Noctuidae) and predators in cotton field

Field trials on upland cotton (Gossypium hirstum L.) during its reproductive phase were used to assess the toxicity of several biorational pesticides and chemicals to Helicoverpa armigera (Hubner) and H. puntigera Wallengren, as well as major predators at Dalby, Queensland, Australia. Moderate rate-dependent control was obtained in plots treated with neem (Azadirachta indica A. Juss) seed extract-azadirachtin (Aza) at rates of 30, 60 and 90 g/ha. Plots treated with Talstar EC (bifenthrin) applications achieved the best results, followed by treatment with alternation of chemicals (methomyl, bifenthrin, thiodicarb and endosulfan) and biorational insecticides (neem oil, azadirachtin and Bacillus thuringiensis kurstaki var. Berliner). Predators, including lady beetles, lacewings, spiders and predatory bugs, were insensitive to Aza, tooseendanin (Tsdn) and BT applications. In contrast, chemicals were very destructive of predators. All treatments provided some protection from infestation of H. armigera and H. puntigera. The effect of Aza on Helicoverpa spp. was reflected in a relatively higher yield of seed cotton harvested from Aza-treated plots compared with the control, but chemical control achieved significantly higher yields than any other treatment.

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.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.