A finite element analysis of quasistatic crack growth in a pressure sensitive constrained ductile layer

Polymeric adhesive layers are employed for bonding two components in a wide variety of technological applications, It has been observed that, unlike in metals, the yield behavior of polymers is affected by the state of hydrostatic stress. In this work, the effect of pressure sensitivity of yielding and layer thickness on quasistatic interfacial crack growth in a ductile adhesive layer is investigated. To this end, finite deformation, finite element analyses of a cracked sandwiched layer are carried out under plane strain, small-scale yielding conditions for a wide range of mode mixities. The Drucker-Prager constitutive equations are employed to represent the behavior of the layer. Crack propagation is simulated through a cohesive zone model, in which the interface is assumed to follow a prescribed traction-separation law. The results show that for a given mode mixity, the steady state Fracture toughness [K](ss) is enhanced as the degree of pressure sensitivity increases. Further, for a given level of pressure sensitivity, [K](ss) increases steeply as mode Il loading is approached. (C) 2000 Elsevier Science Ltd. All rights reserved.

Effect of interface friction on the mechanics of indentation of a finite layer resting on a rigid substrate

Copper strips of 2.5 mm thickness resting on stainless steel anvils were normally indented by wedges under nominal plane strain conditions. Inflections in the hardness-penetration characteristics were identified. Inflections separate stages where each stage has typical mechanics of deformation. These are arrived at by studying the distortion of 0.125 mm spaced grids inscribed on the deformation plane of the strip. The sensitivity of hardness and deformation mechanics to wedge angle and the interfacial friction between strip and anvil were investigated within the framework of existing slip line field models of indentation of semi-infinite and finite blocks.

A finite element analysis of dynamic fracture initiation by ductile failure mechanisms in a 4340 steel

In some recent dropweight impact experiments [5] with pre-notched bend specimens of 4340 steel, it was observed that considerable crack tunneling occurred in the interior of the specimen prior to gross fracture initiation on the free surfaces. The final failure of the side ligaments happened because of shear lip formation. The tunneled region is characterized by a flat, fibrous fracture surface. In this paper, the experiments of [5] (corresponding to 5 m/s impact speed) are analyzed using a plane strain, dynamic finite element procedure. The Gurson constitutive model that accounts for the ductile failure mechanisms of micro-void nucleation, growth and coalescence is employed. The time at which incipient failure was observed near the notch tip in this computation, and the value of the dynamic J-integral, J d, at this time, compare reasonably well with experiments. This investigation shows that J-controlled stress and deformation fields are established near the notch tip whenever J d , increases with time. Also, it is found that the evolution of micro-mechanical quantities near the notch root can be correlated with the time variation of J d .The strain rate and the adiabatic temperature rise experienced at the notch root are examined. Finally, spatial variations of stresses and deformations are analyzed in detail.

Dynamic growth of tensile cracks by ductile and brittle fracture mechanisms in a viscoplastic material

In this paper, a finite element analysis of steady-state dynamic crack growth under Mode I, plane strain, small-scale yielding conditions is performed in a rate dependent plastic material characterized by the over-stress model. The main objective of the paper is to obtain theoretically the dependence of dynamic fracture toughness on crack speed. Crack propagation due to a ductile (micro-void) mechanism or a brittle (cleavage) mechanism, as well as transition from one mode to another are considered. The conversion from ductile to brittle has been observed experimentally but has received very little attention using analytical methods. Local fracture criteria based on strains and stresses are used to describe ductile and brittle fracture mechanisms. The results obtained in this paper are in general agreement with micro-structural observations of mode conversion during fracture initiation. Finally, the particular roles played by material rate sensitivity and inertia are examined in some detail.

Influence of Inclination Angle and Machining Direction on Friction and Transfer Layer Formation

In the present investigation, unidirectional grinding marks were created on a set of steel plates. Sliding experiments were then conducted with the prepared steel plates using Al-Mg alloy pins and an inclined pin-on-plate sliding tester. The goals of the experiments were to ascertain the influence of inclination angle and grinding mark direction on friction and transfer layer formation during sliding contact. The inclination angle of the plate was held at 0.2 deg, 0.6 deg, 1 deg, 1.4 deg, 1.8 deg, 2.2 deg, and 2.6 deg in the tests. The pins were slid both perpendicular and parallel to the grinding marks direction. The experiments were conducted under both dry and lubricated conditions on each plate in an ambient environment. Results showed that the coefficient of friction and the formation of transfer layer depend on the grinding marks direction and inclination angle of the hard surfaces. For a given inclination angle, under both dry and lubricated conditions, the coefficient of friction and transfer layer formation were found to be greater when the pins slid perpendicular to the unidirectional grinding marks than when the pins slid parallel to the grinding marks. In addition, a stick-slip phenomenon was observed under lubricated conditions at the highest inclination angle for sliding perpendicular to the grinding marks direction. This phenomenon could be attributed to the extent of plane strain conditions taking place at the asperity level during sliding. DOI: 10.1115/1.4002604]

Role of Surface Texture, Roughness, and Hardness on Friction During Unidirectional Sliding

In the present investigation, experiments were conducted by unidirectional sliding of pins made of FCC metals (Pb, Al, and Cu) with significantly different hardness values against the steel plates of various surface textures and roughness using an inclined pin-on-plate sliding apparatus in ambient conditions under both the dry and lubricated conditions. For a given material pair, it was observed that transfer layer formation and the coefficient of friction along with its two components, namely adhesion and plowing, are controlled by the surface texture of the harder mating surfaces and are less dependent of surface roughness (R (a)) of the harder mating surfaces. The effect of surface texture on the friction was attributed to the variation of the plowing component of friction for different surfaces. It was also observed that the variation of plowing friction as a function of hardness depends on surface textures. More specifically, the plowing friction varies with hardness of the soft materials for a given type of surface texture and it is independent of hardness of soft materials for other type of surface texture. These variations could be attributed to the extent of plane strain conditions taking place at the asperity level during sliding. It was also observed that among the surface roughness parameters, the mean slope of the profile, Delta (a), correlated best with the friction. Furthermore, dimensionless quantifiable roughness parameters were formulated to describe the degree of plowing taking place at the asperity level.

Asymptotic and finite element analyses of mode III dynamic crack growth at a ductile-brittle interface

In this work, dynamic crack growth along a ductile-brittle interface under anti-plane strain conditions is studied. The ductile solid is taken to obey the J(2) flow theory of plasticity with linear isotropic strain hardening, while the substrate is assumed to exhibit linear elastic behavior. Firstly, the asymptotic near-tip stress and velocity fields are derived. These fields are assumed to be variable-separable with a power singularity in the radial coordinate centered at the crack tip. The effects of crack speed, strain hardening of the ductile phase and mismatch in elastic moduli of the two phases on the singularity exponent and the angular functions are studied. Secondly, full-field finite element analyses of the problem under small-scale yielding conditions are performed. The validity of the asymptotic fields and their range of dominance are determined by comparing them with the results of the full-field finite element analyses. Finally, theoretical predictions are made of the variations of the dynamic fracture toughness with crack velocity. The influence of the bi-material parameters on the above variation is investigated.

Finite element simulation of wedge indentation

A finite element simulation of frictionless wedge indentation of a copper strip has been carried out under plane strain conditions. The problem was first modelled using an one-pass contact algorithm. The difficulties associated with using this method to model wedge indentation problems are explained. An alternative procedure which alleviates some of the problems associated with the one-pass contact algorithm is proposed for modelling frictionless wedge indentation. Also, a re-meshing procedure which has to be carried out when the distortion of the elements around the indenter becomes significant, is discussed. A sample problem involving indentation of a 4 mm copper strip by a rigid wedge indenter has been modelled and the results are compared with experimental and theoretical results.

Finite element simulations of ductile rupture in a constrained metal foil

A numerical study of the ductile rupture in a metal foil constrained between two stiff ceramic blocks is performed. The finite element analysis is carried out under the conditions of mode I, plane strain, small-scale yielding. The rate-independent version of the Gurson model that accounts for the ductile failure mechanisms of microvoid nucleation, growth and coalescence is employed to represent the behavior of the metal foil. Different distributions of void nucleating sites in the metal foil are considered for triggering the initiation of discrete voids. The results clearly show that far-field triaxiality-induced cavitation is the dominant failure mode when the spacing of the void nucleating sites is large. On the contrary, void coalescence near the notch tip is found to be the operative failure mechanism when closely spaced void nucleating sites are considered.

Variability of sif and cod of stochastic structural systems

Stochastic structural systems having a stochastic distribution of material properties and stochastic external loadings in space are analysed when a crack of deterministic size is present. The material properties and external loadings are considered to constitute independent, two-dimensional, univariate, real, homogeneous stochastic fields. The stochastic fields are characterized by their means, variances, autocorrelation functions or the equivalent power spectral density functions, and scale fluctuations. The Young's modulus and Poisson's ratio are treated to be stochastic quantities. The external loading is treated to be a stochastic field in space. The energy release rate is derived using the method of virtual crack extension. The deterministic relationship is derived to represent the sensitivities of energy release rate with respect to both virtual crack extension and real system parameter fluctuations. Taylor series expansion is used and truncation is made to the first order. This leads to the determination of second-order properties of the output quantities to the first order. Using the linear perturbations about the mean values of the output quantities, the statistical information about the energy release rates, SIF and crack opening displacements are obtained. Both plane stress and plane strain cases are considered. The general expressions for the SIF in all the three fracture modes are derived and a more detailed analysis is conducted for a mode I situation. A numerical example is given.

Plate-shaped transformation products in zirconium-base alloys

Plate-shaped products resulting from martensitic, diffusional, and mixed mode transformations in zirconium-base alloys are compared. in the present study. These alloys are particularly suitable for the comparison in view of the fact that the lattice correspondence between the parent beta (bcc) and the product alpha (hcp) or gamma-hydride (fct) phases are remarkably similar for different types of transformations. Crystallographic features such as orientation relations, habit planes, and interface structures associated with these transformations have been compared:, with a view toward examining whether the transformation mechanisms have characteristic imprints on these experimental observables. Martensites exhibiting dislocated lath, internally twinned plate, and self-accommodating three-plate cluster morphologies have been encountered in Zr-2.5Nb alloy. Habit planes corresponding to all these morphologies have been found to be consistent with the predictions based on the invariant plane strain (IFS) criterion. Different morphologies have been found to reflect the manner in which the neighboring martensite variants are assembled. Lattice-invariant shears (LISs) for all these cases have been identified to be either {10 (1) over bar 1}(alpha) ((1) over bar 123)(alpha) slip or twinning on (10 (1) over bar 1)(alpha) planes. Widmanstatten alpha precipitates, forming in a step-quenching treatment, have been shown to have a lath morphology, the alpha/beta interface being decorated with a periodic array of (c + a) dislocations at a spacing of 8 to 10 nm. The line vectors of these dislocations are nearly parallel to the invariant lines. The alpha precipitates, forming in the retained beta phase on aging, exhibit an internally twinned structure with a zigzag habit plane. Average habit planes for the morphologies have been found to lie near the {103}(beta) - {113}(beta) poles, which are close to the specific variant of the {112}(beta) plane, which transforms into a prismatic plane of the type {1 (1) over bar 00}(alpha). The crystallography of the formation of the gamma-hydride phase (fct) from both the alpha and beta phases is seen to match the IFS predictions. While the beta-gamma transformation can be treated approximately as a simple shear on the basal plane involving a change in the stacking sequence, the alpha-gamma transformation call be conceptually broken into a alpha --> beta transformation following the Burgers correspondence and the simple beta-gamma shear process. The active eutectoid decomposition in the Zr-Cu system, beta --> alpha + beta', has been described in terms of cooperative growth of the alpha phase from the beta phase through the Burgers correspondence and of the partially ordered beta' (structurally similar to the equilibrium Zr2Cu phase) through an ordering process. Similarities and differences in crystallographic features of these transformations have been discussed. and the importance of the invariant line vector in deciding the geometry of the corresponding habit planes has been pointed out.

Constraint loss under dynamic loading in rate independent plastic solids

The objectives of this paper are to examine the loss of crack tip constraint in dynamically loaded fracture specimens and to assess whether it can lead to enhancement in the fracture toughness at high loading rates which has been observed in several experimental studies. To this end, 2-D plane strain finite element analyses of single edge notched (tension) specimen and three point bend specimen subjected to time varying loads are performed. The material is assumed to obey the small strain J(2) flow theory of plasticity with rate independent behaviour. The results demonstrate that a valid J-Q field exists under dynamic loading irrespective of the crack length and specimen geometry. Further, the constraint parameter Q becomes strongly negative at high loading rates, particularly in deeply cracked specimens. The variation of dynamic fracture toughness K-dc with stress intensity rate K for cleavage cracking is predicted using a simple critical stress criterion. It is found that inertia-driven constraint loss can substantially enhance K-dc for (K) over dot > 10(5) MPa rootm/s.

Symmetry-breaking transitions in equilibrium shapes of coherent precipitates: Effect of elastic anisotropy and inhomogeneity

We examine the symmetry-breaking transitions in equilibrium shapes of coherent precipitates in two-dimensional (2-D) systems under a plane-strain condition with the principal misfit strain components epsilon(xx)*. and epsilon(yy)*. For systems with cubic elastic moduli, we first show all the shape transitions associated with different values of t = epsilon(yy)*/epsilon(xx)*. We also characterize each of these transitions, by studying its dependence on elastic anisotropy and inhomogeneity. For systems with dilatational misfit (t = 1) and those with pure shear misfit (t = -1), the transition is from an equiaxed shape to an elongated shape, resulting in a break in rotational symmetry. For systems with nondilatational misfit (-1 < t < 1; t not equal 0), the transition involves a break in mirror symmetries normal to the x- and y-axes. The transition is continuous in all cases, except when 0 < t < 1. For systems which allow an invariant line (-1 less than or equal to t < 0), the critical size increases with an increase in the particle stiffness. However, for systems which do not allow an invariant line (0 < t less than or equal to 1), the critical size first decreases, reaches a minimum, and then starts increasing with increasing particle stiffness; moreover, the transition is also forbidden when the particle stiffness is greater than a critical value.

Effect of lattice orientation on crack tip constraint in ductile single crystals

In this work, the effect of lattice orientation on the fields prevailing near a notch tip is investigated pertaining to various constraint levels in FCC single crystals. A modified boundary layer formulation is employed and numerical solutions under mode I, plane strain conditions are generated by assuming an elastic-perfectly plastic FCC single crystal. The analysis is carried out corresponding to different lattice orientations with respect to the notch line. It is found that the near-tip deformation field, especially the development of kink or slip shear bands is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the notch tip are also strongly influenced by the level of T-stress. The present results clearly establish that ductile single crystal fracture geometries would progressively lose crack tip constraint as the T-stress becomes more negative irrespective of lattice orientation. Also, the near-tip field for a range of constraint levels can be characterized by two-parameters such as K-T or J-Q as in isotropic plastic solids.

Design of Geocell Reinforcement for Supporting Embankments on Soft Ground

The methods of design available for geocell-supported embankments are very few. Two of the earlier methods are considered in this paper and a third method is proposed and compared with them. The first method is the slip line method proposed by earlier researchers. The second method is based on slope stability analysis proposed by this author earlier and the new method proposed is based on the finite element analyses. In the first method, plastic bearing failure of the soil was assumed and the additional resistance due to geocell layer is calculated using a non-symmetric slip line field in the soft foundation soil. In the second method, generalpurpose slope stability program was used to design the geocell mattress of required strength for embankment using a composite model to represent the shear strength of geocell layer. In the third method proposed in this paper, geocell reinforcement is designed based on the plane strain finite element analysis of embankments. The geocell layer is modelled as an equivalent composite layer with modified strength and stiffness values. The strength and dimensions of geocell layer is estimated for the required bearing capacity or permissible deformations. These three design methods are compared through a design example.