Crack growth by grain boundary cavitation in the transient and extensive creep regimes

Crack growth due to cavity growth and coalescence along grain boundaries is analyzed under transient and extensive creep conditions in a compact tension specimen. Account is taken of the finite geometry changes accompanying crack tip blunting. The material is characterized as an elastic-power law creeping solid with an additional contribution to the creep rate arising from a given density of cavitating grain boundary facets. All voids are assumed present from the outset and distributed on a given density of cavitating grain boundary facets. The evolution of the stress fields with crack growth under three load histories is described in some detail for a relatively ductile material. The full-field plane strain finite element calculations show the competing effects of stress relaxation due to constrained creep, diffusion and crack tip blunting. and of stress increase due to the instantaneous elastic response to crack growth. At very high crack growth rates the Hui-Riedel fields dominate the crack tip region. However. the high growth rates are not sustained for any length of time in the compact tension geometry analyzed. The region of dominance of the Hui-Riedel field shrinks rapidly so that the near-tip fields are controlled by the HRR-type field shortly after the onset of crack growth. Crack growth rates under various conditions of loading and spanning the range of times from small scale creep to extensive creep are obtained. We show that there is a strong similarity between crack growth history and the behaviour of the C(t) and C(t) parameters. so that crack growth rates correlate rather well with C(t) and C(t). A relatively brittle material is also considered that has a very different near-tip stress field and crack growth history.

On dilatational plastic constitutive equation of ductile materials and plastic loading paths at bifurcation

The dilatational plastic constitutive equation presented in this paper is proved to be in a form of generality. Based on this equation, the constitutive behaviour of materials at the moment of bifurcation is demonstrated to follow a loading path with the response as "soft" as possible.

Influence of secondary void damage in the matrix material around voids

The mechanism of ductile damage caused by secondary void damage in the matrix around primary voids is studied by large strain, finite element analysis. A cylinder embedding an initially spherical void, a plane stress cell with a circular void and plane strain cell with a cylindrical or a flat void are analysed under different loading conditions. Secondary voids of smaller scale size nucleate in the strain hardening matrix, according to the requirements of some stress/strain criteria. Their growth and coalescence, handled by the empty element technique, demonstrate distinct mechanisms of damage as circumstances change. The macroscopic stress-strain curves are decomposed and illustrated in the form of the deviatoric and the volumetric parts. Concerning the stress response and the void growth prediction, comparisons are made between the present numerical results and those of previous authors. It is shown that loading condition, void growth history and void shape effect incorporated with the interaction between two generations of voids should be accounted for besides the void volume fraction.

激光辐照热障涂层中平面应变问题的热弹性变分分析

针对激光辐照热障涂层材料的平面应变问题,提出热障涂层热弹性分析的基本方程,对定常温度场给出级数形式解析解,并用最小余能原理和变分法分析了结构的热弹性应力场,研究了最大应力和界面应力的分布特征,并就一些物理参数的影响进行了讨论.结果表明,热障涂层的主要破坏因素为表面拉伸应力,界面应力相对较小,但在自由边界有集中现象,剥落应力大于剪切应力,是导致涂层破坏的重要原因.涂层厚度增加会改变厚度方向上的应力分布,界面应力向中心集中.

Catastrophic Rupture Induced Damage Coalescence in Heterogeneous Brittle Media

In heterogeneous brittle media, the evolution of damage is strongly influenced by the multiscale coupling effect. To better understand this effect, we perform a detailed investigation of the damage evolution, with particular attention focused on the catastrophe transition. We use an adaptive multiscale finite-element model (MFEM) to simulate the damage evolution and the catastrophic failure of heterogeneous brittle media. Both plane stress and plane strain cases are investigated for a heterogeneous medium whose initial shear strength follows the Weibull distribution. Damage is induced through the application of the Coulomb failure criterion to each element, and the element mesh is refined where the failure criterion is met. We found that as damage accumulates, there is a stronger and stronger nonlinear increase in stress and the stress redistribution distance. The coupling of the dynamic stress redistribution and the heterogeneity at different scales result in an inverse cascade of damage cluster size, which represents rapid coalescence of damage at the catastrophe transition.

An extension of the two dimensional JKR theory to the case with a large contact width

In the Hertz and JKR theories, parabolic assumptions for the rounded profiles of the sphere or cylinder are adopted under the condition that the contact radius (width) should be very small compared to the radius of the sphere or cylinder. However, a large contact radius (width) is often found in experiments even under a zero external loading. We aim at extending the plane strain JKR theory to the case with a large contact width. The relation between the external loading and the contact width is given. Solutions for the Hertz, JKR and rounded-profile cases are compared and analyzed. It is found that when the ratio of a/R is approximately larger than about 0.4, the parabolic assumptions in the Hertz and JKR theories are no longer valid and the exact rounded profile function should be used.

A discrete dislocation analysis of fatigue crack growth

Cyclic loading of a plane strain mode I crack under small scale yielding is analyzed using discrete dislocation dynamics. The dislocations are all of edge character, and are modeled as line singularities in an elastic solid. At each stage of loading, superposition is used to represent the solution in terms of solutions for edge dislocations in a half-space and a non-singular complementary solution that enforces the boundary conditions, which is obtained from a linear elastic, finite element solution. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation are incorporated into the formulation through a set of constitutive rules. An irreversible relation between the opening traction and the displacement jump across a cohesive surface ahead of the initial crack tip is also specified, which permits crack growth to emerge naturally. It is found that crack growth can occur under cyclic loading conditions even when the peak stress intensity factor is smaller than the stress intensity required for crack growth under monotonic loading conditions; however below a certain threshold value of ΔKI no crack growth was seen.

Discrete dislocation plasticity analysis of crack-tip fields in polycrystalline materials

Small scale yielding around a mode I crack is analysed using polycrystalline discrete dislocation plasticity. Plane strain analyses are carried out with the dislocations all of edge character and modelled as line singularities in a linear elastic material. The lattice resistance to dislocation motion, nucleation, interaction with obstacles and annihilation are incorporated through a set of constitutive rules. Grain boundaries are modelled as impenetrable to dislocations. The polycrystalline material is taken to consist of two types of square grains, one of which has a bcc-like orientation and the other an fcc-like orientation. For both orientations there are three active slip systems. Alternating rows, alternating columns and a checker-board-like arrangement of the grains is used to construct the polycrystalline materials. Consistent with the increasing yield strength of the polycrystalline material with decreasing grain size, the calculations predict a decrease in both the plastic zone size and the crack-tip opening displacement for a given applied mode I stress intensity factor. Furthermore, slip-band and kink-band formation is inhibited by all grain arrangements and, with decreasing grain size, the stress and strain distributions more closely resemble the HRR fields with the crack-tip opening approximately inversely proportional to the yield strength of the polycrystalline materials. The calculations predict a reduction in fracture toughness with decreasing grain size associated with the grain boundaries acting as effective barriers to dislocation motion.

Modelling of waterfront cantilever sheet pile walls subjected to earthquake loading

The seismic performance of waterfront cantilever sheet pile retaining walls is of continuing interest to geotechnical engineers as these structures suffer severe damage and even complete failure during earthquakes. This is often precipitated by liquefaction of the surrounding soil, either in the backfill or in front of the wall. This paper presents results from a series of small-scale plane strain models that were tested on a 1-g shaking table and recorded using a high-speed, high-resolution digital camera. The technique of Particle Image Velocimetry (PIV) was applied in order to allow the failure mechanisms to be visualised. It is shown that using PIV analyses it is possible to obtain failure mechanisms for a cantilever wall in liquefiable soil. These failure mechanisms are compared with those obtained for a cantilever wall in dry soil, previously carried out at a similar scale. It was observed that seismic liquefaction causes significant displacement in much larger zones of soil near the retaining wall compared to an equivalent dry case. The failure mechanism for a cantilever wall with liquefiable backfill, but with a remediated zone designed not to liquefy, is also presented and compared to the unremediated case.

Seismic behaviour of anchored quay walls with dry backfill

The seismic behaviour of anchored sheet pile walls is a complex soil-structure interaction problem. Damaged sheet pile walls are very expensive to repair and their seismic behaviour needs to be investigated in order to understand their possible mechanisms of failure. The research described in this paper involves both centrifuge testing and Finite Element (FE) analyses aimed at investigating the seismic behaviour of an anchored sheet pile wall in dry sand. The model wall is tied to the backfill with two tie rods connected to an anchor beam. The accelerations of the sheet pile wall, the anchor beam and the soil around the wall were measured using miniature piezoelectric accelerometers. The displacement at the tip of the wall was also measured. Stain gauges at five different locations on the wall were used to measure the bending moments induced in the the wall. The anchor forces in the tie rods were also measured using load cells. The results from the centrifuge tests were compared with 2-D, plane strain FE analyses conducted using DIANA-SWANDYNE II and the observed seismic behaviour was explained in the light of these findings. © 2011 Taylor & Francis.

Investigation of the poro-elastic response of seabed to solitary

A numerical model is established and validated to study the behavior of porous seabed under solitary wave propagation. Using Biot's poro-elastic theory, the problem is formulated as a two dimensional plane strain problem, and it is modelled using the Finite Element Method. The responses due to the solitary wave are compared with those of linear waves of the same height. It is found that regardless of the wave period, stresses due to solitary waves are generally larger. This indicates a higher potential for shear failure at the seabed under solitary waves. Implications on liquefaction need further investigation. Copyright © 2012 by the International Society of Offshore and Polar Engineers (ISOPE).

Failure mechanisms in brittle laminates

Steady-state tunneling and plane-strain delamination of an H-shape crack are examined for elastic, isotropic multi layers. Both tunneling and delamination are analysed by employing linear elastic fracture mechanics within a 2D finite element framework. Failure maps are produced to reveal the sensitivity of cracking path to the relative toughness of layer and interface, and to the stiffness mismatch of layers.

Insight into differences in nanoindentation properties of bone.

Nanoindentation provides the ideal framework to determine mechanical properties of bone at the tissue scale without being affected by the size, shape, and porosity of the bone. However, the values of tissue level mechanical properties vary significantly between studies. Since the differences in the bone sample, hydration state, and test parameters complicate direct comparisons across the various studies, these discrepancies in values cannot be compared directly. The objective of the current study is to evaluate and compare mechanical properties of the same bones using a broad range of testing parameters. Wild type C56BL6 mice tibiae were embedded following different processes and tested in dry and rehydrated conditions. Spherical and Berkovich indenter probes were used, and data analysis was considered within the elasto-plastic (Oliver-Pharr), viscoelastic and visco-elastic-plastic frameworks. The mean values of plane strain modulus varied significantly depending on the hydration state, probe geometry and analysis method. Indentations in dry bone analyzed using a visco-elastic-plastic approach gave values of 34 GPa. After rehydrating the same bones and indenting them with a spherical tip and utilizing a viscoelastic analysis, the mean modulus value was 4 GPa, nearly an order of magnitude smaller. Results suggest that the hydration state, probe geometry and the limitations and assumptions of each analysis method influence significantly the measured mechanical properties. This is the first time that such a systematic study has been carried out and it has been concluded that the discrepancies in the mechanical properties of bone measured by nanoindentation found in the literature should not be attributed only to the differences between the bones themselves, but also to the testing and analysis protocols.

An efficient model for the dynamic behaviour of a single pile in viscoelastic soil

This paper presents a novel, three-dimensional, single-pile model, formulated in the wavenumber domain and adapted to account for boundary conditions using the superposition of loading cases. The pile is modelled as a column in axial vibration, and a Euler-Bernoulli beam in lateral vibration. The surrounding soil is treated as a viscoelastic continuum. The response of the pile is presented in terms of the stiffness and damping coefficients, and also the magnitude and phase of the pile-head frequency-response function. Comparison with existing models shows that excellent agreement is observed between this model, a boundary-element formulation, and an elastic-continuum-type formulation. This three-dimensional model has an accuracy equivalent to a 3D boundary-element model, and a runtime similar to a 2D plane-strain analytical model. Analysis of the response of the single pile illustrates a difference in axial and lateral vibration behaviour; the displacement along the pile is relatively invariant under axial loads, but in lateral vibration the pile exhibits localised deformations. This implies that a plane-strain assumption is valid for axial loadings and only at higher frequencies for lateral loadings. © 2013 Elsevier Ltd.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.