Constraint effects in adhesive joint fracture

Constraint effects in adhesive joint fracture are investigated by modelling the adherents as well as a finite thickness adhesive layer in which a single row of cohesive zone elements representing the fracture process is embedded. Both the adhesive and the adherents are elastic-plastic with strain hardening. The bond toughness Γ (work per unit area) is equal to Γ₀+Γ_p, where Γ₀ is the intrinsic work of fracture associated with the embedded cohesive zone response and Γ_p is the extra contribution to the bond toughness arising from plastic dissipation and stored elastic energy within the adhesive layer. The parameters of the model are identified from experiments on two different adhesives exhibiting very different fracture properties. Most of the tests were performed using the wedge-peel test method for a variety of adhesives, adherents and wedge thicknesses. The model captures the constraint effects resulting from the change in Γ_p: (i) the plastic dissipation increases with increasing bond line thickness in the fully plastic regime and then decreases to reach a constant value for very thick adhesive layers; (ii) the plastic dissipation in the fully plastic regime increases drastically as the thickness of the adherent decreases. Finally, this model is used to assess a simpler approach which consists of simulating the full adhesive layer as a single row of cohesive elements.     

Finite strain analysis of crack tip fields in incompressible hyperelastic solids loaded in plane stress

A new method is developed to determine the dominant asymptotic stress and deformation fields near the tip of a Mode-I traction free plane stress crack. The analysis is based on the fully nonlinear equilibrium theory of incompressible hyperelastic solids. We show that the dominant singularity of the near tip stress field is governed by the asymptotic solution of a linear second order ordinary differential equation. Our method is applicable to any hyperelastic material with a smooth work function that depends only on the trace of the Cauchy-Green tensor and is particularly useful for materials that exhibit severe strain hardening. We apply this method to study two types of soft materials: generalized neo-Hookean solids and a solid that hardens exponentially. For the generalized neo-Hookean solids, our method is able to resolve a difficulty in the previous work by Geubelle and Knauss (1994a). Our theoretical results are compared with finite element simulations.     

Asymptotic mechanical fields at the tip of a mode I crack in rubber-like solids

Asymptotic analyses of the mechanical fields in front of stationary and propagating cracks facilitate the understanding of the mechanical and physical state in front of crack tips, and they enable prediction of crack growth and failure. Furthermore, efficient modelling of arbitrary crack growth by use of XFEM (extended finite element method) requires accurate knowledge of the asymptotic crack tip fields. In the present work, we perform an asymptotic analysis of the mechanical fields in the vicinity of a propagating mode I crack in rubber. Plane deformation is assumed, and the material model is based on the Langevin function, which accounts for the finite extensibility of polymer chains. The Langevin function is approximated by a polynomial, and only the term of the highest order contributes to the asymptotic solution. The crack is predicted to adopt a wedge-like shape, i.e. the crack faces will be straight lines. The angle of the wedge and the order of the stress singularity depend on the hardening of the strain energy function. The present analysis shows that in materials with a significant hardening, the inertia term in the equations of motion becomes negligible in the asymptotic analysis. Hence, there is no upper theoretical limit to the crack speed.     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

Piezomagnetic and piezoelectric poling effects on mode I and II crack initiation behavior of magnetoelectroelastic materials

The ferrite and ferroelectric phase of magnetoelectroelastic (MEE) material can be selected and processed to control the macroscopic behavior of electron devices using continuum mechanics models. Once macro- and/or microdefects appear, the highly intensified magnetic and electric energy localization could alter the response significantly to change the design performance. Alignment of poling directions of piezomagnetic and piezoelectric materials can add to the complexity of the MEE material behavior to which this study will be concerned with.Appropriate balance of distortional and dilatational energy density is no longer obvious when a material possesses anisotropy and/or nonhomogeneity. An excess of the former could result in unwanted geometric change while the latter may lead to unexpected fracture initiation. Such information can be evaluated quantitatively from the stationary values of the energy density function dW/dV. The maxima and minima have been known to coincide, respectively, with possible locations of permanent shape change and crack initiation regardless of material and loading type. The direction of poling with respect to a line crack and the material microstructure described by the constitutive coefficients will be specified explicitly with reference to the applied magnetic field, electric field and mechanical stress, both normal and shear. The crack initiation load and direction could be predicted by finding the direction for which the volume change is the largest. In contrast to intuition, change in poling directions can influence the cracking behavior of MEE dramatically. This will be demonstrated by the numerical results for the BaTiO₃–CoFe₂O₄ composite having different volume fractions where BaTiO₃ and CoFe₂O₄ are, respectively, the inclusion and matrix.To be emphasized is that mode I and II crack behavior will not have the same definition as that in classical fracture mechanics where load and crack extension symmetry would coincide. A striking result is found for a mode II crack. By keeping the magnetic poling fixed, a reversal of electric poling changed the crack initiation angle from θ₀=+80° to θ₀=−80° using the line extending ahead of the crack as the reference. This effect is also sensitive to the distance from the crack tip. Displayed and discussed are results for r/a=10⁻⁴ and 10⁻¹. Because the theory of magnetoelectroelasticity used in the analysis is based on the assumption of equilibrium where the influence of material microstructure is homogenized, the local space and temporal effects must be interpreted accordingly. Among them are the maximum values of (dW/dV)_max and (dW/dV)_min which refer to as possible sites of yielding and fracture. Since time and size are homogenized, it is implicitly understood that there is more time for yielding as compared to fracture being a more sudden process. This renders a higher dW/dV in contrast to that for fracture. Put it differently, a lower dW/dV with a shorter time for release could be more detrimental.     

Investigation on critical triaxiality along unified yield loci at crack tip under mixed mode (I + II) fracture

Metallic components used in industries and day to day appliances often contain micro-cracks. In general, cracks occur in various orientations to the loading axis. The present paper discusses the criticality of stress triaxiality, a well-known ductile fracture parameter, on the yield loci at the crack tip. In the process, an old model of stress triaxiality has been generalized using unified strength theory to incorporate various convex and nonconvex failure criteria, including single shear, twin shear, etc. The new triaxiality model also reveals about the effect of intermediate principal stress at the crack tip for materials with and without strength difference. The crack initiation angles at the crack tip, obtained through the proposed model have been found to be in unison with those obtained through other fracture criteria.     

Interaction between crack tip advancement and fluid flow in fracturing saturated porous media

We address stepwise crack tip advancement and pressure fluctuations, which have been observed in the field and experimentally in fracturing saturated porous media. Both fracturing due to mechanical loading and pressure driven fracture are considered. After presenting the experimental evidence and the different explanations for the phenomena put forward and mentioning briefly what has been obtained so far by published numerical and analytical methods we propose our explanation based on Biot’s theory. A short presentation of three methods able to simulate the observed phenomena namely the Central Force Model, the Standard Galerkin Finite Element Method SGFEM and extended finite element method XFEM follows. With the Central Force Model it is evidenced that already dry geomaterials break in an intermittent fashion and that the presence of a fluid affects the behavior more or less depending on the loading and boundary conditions. Examples dealing both with hydraulic fracturing and mechanical loading are shown. The conditions needed to reproduce the observed phenomena with FE models at macroscopic level are evidenced. They appear to be the adoption of a crack tip advancement/time step algorithm which interferes the least possible with the three interacting velocities, namely the crack tip advancement velocity on one side, the seepage velocity of the fluid in the domain and from the crack (leak-off), and the fluid velocity within the crack on the other side. Further the crack tip advancement algorithm must allow for reproducing jumps observed in the experiments.     

复合材料Ⅰ、Ⅱ混合型分层试验研究

通过对T300／5405、1300／913、1300／HD03三种航空材料Ⅰ、Ⅱ混合型分层研究，进一步研究了新的Ⅰ、Ⅱ混合型分层断裂测试的JJ法，并和MMB法进行了比较，发现两者较为吻合。但JJ法耗费试件少，简单易行，操作方便，有利于建立混合型分层的失效判据，它可以适用于不同环境下的混合型分层实验的测试。     

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅰ)—THE ANALYSES OF CONSTRAINT PARAMETERS AND FRACTURE PARAMETERS

In the present paper, three dimensional analyses of some general constraint parameters and fracture parameters near the crack tip. of Mode I CT specimens in two different thicknesses are carried out by employing ADINA program. The results reveal that the constraints along the thickness direction are obviously separated into two parts: the keeping similar high constraint field (Z₁) and rapid reducing constraints one (Z₂). The two fields are experimentally confiremed to correspond to the smooth region and the shear lip on the fracture face respectively. So the three dimensional stress structure of Mode I specimens can be derived through discussing the two fields respectively. The distribution of the Crack Tip Opening Displacement (CTOD) along the thickness direction and the three dimensional distribution of the void growth ratio (V_g) near the crack tip are also obtained. The two fracture parameters are in similar trends along the thickness direction, and both of them can reflect the effect of thickness and that of the loading level to a certain degree.     

Coupled continuum and discrete analysis of random heterogeneous materials: Elasticity and fracture

Recent work has suggested that the heterogeneous distribution of mechanical properties in natural and synthetic materials induces a toughening mechanism that leads to a more robust structural response in the presence of cracks, defects or other types of flaws. Motivated by this, we model an elastic solid with a Young′s modulus distribution described by a Gaussian process. We study the pristine system using both a continuum and a discrete model to establish a link between the microscale and the macroscale in the presence of disorder. Furthermore, we analyze a flawed discrete particle system and investigate the influence of heterogeneity on the fracture mechanical properties of the solid. We vary the variability and correlation length of the Gaussian process, thereby gaining fundamental insights into the effect of heterogeneity and the essential length scales of heterogeneity critical to enhanced fracture properties. As previously shown for composites with complex hierarchical architectures, we find that materials with disordered elastic fields toughen by a ‘distribution-of-weakness’ mechanism inducing crack arrest and stress delocalization. In our systems, the toughness modulus can increase by up to 30% due to an increase in variability in the elastic field. Our work presents a foundation for stochastic modeling in a particle-based micromechanical environment that can find broad applications within natural and synthetic materials.     

On strain and damage interactions during tearing: 3D in situ measurements and simulations for a ductile alloy (AA2139-T3)

Strain and damage interactions during tearing of a ductile Al-alloy with high work hardening are assessed in situ and in 3D combining two recently developed experimental techniques, namely, synchrotron laminography and digital volume correlation. Digital volume correlation consists of registering 3D laminography images. Via simultaneous assessments of 3D strain and damage at a distance of 1-mm ahead of a notch root of a thin Compact Tension-like specimen, it is found that parallel crossing slant strained bands are active from the beginning of loading in a region where the crack will be slanted. These bands have an intermittent activity but are stable in space. Even at late stages of deformation strained bands can stop their activity highlighting the importance of plasticity on the failure process rather than damage softening. One void is followed over the loading history and seen to grow and orient along the slant strained band at very late stages of deformation. Void growth and strain are quantified. Gurson–Tvergaard–Needleman-type simulations using damage nucleation for shear, which is based on the Lode parameter, are performed and capture slant fracture but not the initial strain fields and in particular the experimentally found slant bands. The band formation and strain distribution inside and outside the bands are discussed further using plane strain simulations accounting for plastic material heterogeneity in soft zones.     

Development of fracture facets from a crack loaded in mode I+III: Solution and application of a model 2D problem

It is experimentally well-known that a crack loaded in mode I+III propagates through formation of discrete fracture facets inclined at a certain tilt angle on the original crack plane, depending on the ratio of the mode III to mode I initial stress intensity factors. Pollard et al. (1982) have proposed to calculate this angle by considering the tractions on all possible future infinitesimal facets and assuming shear tractions to be zero on that which will actually develop. In this paper we consider the opposite case of well-developed facets; the stress field near the lateral fronts of such facets becomes independent of the initial crack and essentially 2D in a plane perpendicular to the main direction of crack propagation.To determine this stress field, we solve the model 2D problem of an infinite plate containing an infinite periodic array of cracks inclined at some angle on a straight line, and loaded through uniform stresses at infinity. This is done first analytically, for small values of this angle, by combining Muskhelishvili's (1953) formalism and a first-order perturbation procedure. The formulae found for the 2D stress intensity factors are then extended in an approximate way to larger angles by using another reference solution, and finally assessed through comparison with some finite element results.To finally illustrate the possible future application of these formulae to the prediction of the stationary tilt angle, we introduce the tentative assumption that the 2D mode II stress intensity factor is zero on the lateral fronts of the facets. An approximate formula providing the tilt angle as a function of the ratio of the mode III to mode I stress intensity factors of the initial crack is deduced from there. This formula, which slightly depends on the type of loading imposed, predicts somewhat smaller angles than that of Pollard et al. (1982).     

Comments on “Explicit transient solutions for a mode III crack subjected to dynamic concentrated loading in a piezoelectric material” by Yi-Shyong Ing, Mau-Jung Wang [Int. J. Solids Struct. 41 (2004) 3849–3864]

In this comment it is pointed out that the analysis of the dynamic stress intensity factor, dynamic electric displacement intensity factor and dynamic energy release rate conducted by Ing and Wang [Ing, Y.S., Wang, M.J., 2004. Explicit transient solutions for a mode III crack subjected to dynamic concentrated loading in a piezoelectric material. International Journal of Solids and Structures 41, 3849–3864] is incorrect. The correct analysis and corresponding correct plots are presented.     

Extension of the Coherent Gradient Sensor (CGS) to the Combined Measurement of In-Plane and Out-of-Plane Displacement Field Gradients

The Coherent Gradient Sensor (CGS) is extended to the optical differentiation of specular, diffracted wave fronts leading to the combined measurement of in- and out-of-plane displacement field gradients. A derivation of the underlying optical interference principles is presented along with an analysis of the effective instrument sensitivity. In order to demonstrate the capabilities of the technique, experimental measurements of crack-tip deformation fields were conducted under various loading conditions corresponding to mode-I, mode-II, and mixed mode near-tip crack fields. The experimental procedures and results of these tests are presented as validation of the technique.     

Comments on ‘‘Dynamic crack propagation in piezoelectric materials—Part I. Electrode solution’’ by Shaofan Li, Peter A. Mataga [J. Mech. Phys. Solids 44 (1996) 1799-1830]

In this comment, it is pointed out that the paper [Li and Mataga, 1996. J. Mech. Phys. Solids 44, 1799-1830], which presents original and valid solution strategy for an important problem of dynamic crack propagation in piezoelectric materials, contains ultimate quantitatively and qualitatively incorrect expressions, conclusions and plots due calculation errors. The correct calculations and corresponding correct conclusions and plots are presented.     

TOPOLOGY AND VORTEX STRUCTURES OF A CURVING TURBINE CASCADE WITH TIP CLEARANCE (Ⅱ)- TOPOLOGICAL FLOW PATTERN AND VORTEX STRUCTURE IN THE TRANSVERSE SECTION OF A BLADE CASCADE

By means of ink trace visualization of the flows in conventional straight, positively curved and negatively curved cascades with tip clearance, and measurement of the aerodynamic parameters in transverse section, and by appling topology theory, the topological structures and vortex structure in the transverse section of a blade cascade were analyzed. Compared with conventional straight cascade, blade positive curving eliminates the separation line of the upper passage vortex, and leads the secondary vortex to change from close separation to open separation, while blade negative curving effects merely the positions of singular points and the intensities and scales of vortex. Foundation items: 973 Project of China; the Doctoral Foundation of Education Ministry of China (EDAF24403003) Biography: YANG Qing-hai (1969−)     

Numerical simulation for 2D double-diffusive convection (DDC) in rectangular enclosures based on a high resolution upwind compact streamfunction model I: numerical method and code validation

Applied Mathematics and Mechanics - A high resolution upwind compact streamfunction numerical algorithm for two-dimensional (2D) double-diffusive convection (DDC) is developed. The unsteady...