Influence of threshold stress on the estimation of the Weibull statistics

The influence of threshold stress on the estimation of the Weibull statistics is discussed in terms of the Akaike information criterion. Numerical simulations show that, if sample data are limited in number and threshold stress is not too large, the two-parameter Weibull distribution is still a preferred choice. For example, the fit of strength data of glass and ceramics to the two- and three-parameter Weibull distributions is compared.

Electrothermal stress in conductive body with collinear cracks

The temperature and stress field in a thin plate with collinear cracks interrupting an electric current field are determined. This is accomplished by using a complex function method that allows a direct means of finding the distribution of the electric current, the temperature and stress field. Temperature dependency for the heat-transfer coefficient, coefficient of linear expansion and the elastic modulus are considered. As an example, temperature distribution is calculated for an alloy (No. GH2132) plate with two collinear cracks under high temperature. Relationships between the stress, temperature, electric density and crack length are obtained. Crack trajectories emanating from existing crack are predicted by application of the strain energy density criterion which can also be used for finding the load carrying capacity of the cracked plate. (C) 2003 Elsevier Ltd. All rights reserved.

Semi-weight function method on computation of stress intensity factors in dissimilar materials

Semi-weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi-weight functions were obtained as virtual displacement and stress fields with eigenvalue-lambda. Integral expression of fracture parameters, K-I and K-II, were obtained from reciprocal work theorem with semi-weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi-weight function method is a simple, convenient and high precision calculation method.

Singularity of dynamic stress fields around an interface crack between viscoelastic bodies

The singular nature of the dynamic stress fields around an interface crack located between two dissimilar isotropic linearly viscoelastic bodies is studied. A harmonic load is imposed on the surfaces of the interface crack. The dynamic stress fields around the crack are obtained by solving a set of simultaneous singular integral equations in terms of the normal and tangent crack dislocation densities. The singularity of the dynamic stress fields near the crack tips is embodied in the fundamental solutions of the singular integral equations. The investigation of the fundamental solutions indicates that the singularity and oscillation indices of the stress fields are both dependent upon the material constants and the frequency of the harmonic load. This observation is different from the well-known -1/2 oscillating singularity for elastic bi-materials. The explanation for the differences between viscoelastic and elastic bi-materials can be given by the additional viscosity mismatch in the case of viscoelastic bi-materials. As an example, the standard linear solid model of a viscoelastic material is used. The effects of the frequency and the material constants (short-term modulus, long-term modulus and relaxation time) on the singularity and the oscillation indices are studied numerically.

Coupling effects of heterogeneity and stress fluctuation on rupture

The coupling of mesoscopic strength distribution and stress fluctuation on damage evolution and rupture are examined. The numerical simulations show that there is only weak stress fluctuation at the initial damage stage when the mean field approximation is in effect. As the damage fraction becomes larger than the threshold value, the fluctuation is amplified significantly, and damage localization appears. The coupling between stress fluctuation, disordered heterogeneity and the damage localization may play an essential role in catastrophic rupture. (C) 2003 Elsevier Ltd. All rights reserved.

Validity ranges of interfacial wave theories in a two-layer fluid system

In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness d(1), and lower layer thickness d(2), instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehautes plot for free surface waves if water depth ratio r = d(1)/d(2) approaches to infinity and the upper layer water density rho(1) to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of sigma = (rho(2) - rho(1))/rho(2) -> 1.0 and r > 1.0. In the end, several figures of the validity ranges for various interfacial wave theories in the two-layer fluid are given and compared with the results for surface waves.

Finite element analysis on stresses field of normalized layer thickness within ceramic coating on aluminized steel

Multilayer ceramic coatings were fabricated on steel substrate using a combined technique of hot dipping aluminum(HDA) and plasma electrolytic oxidation(PEO). A triangle of normalized layer thickness was created for describing thickness ratios of HDA/PEO coatings. Then, the effect of thickness ratio on stresses field of HDA/PEO coatings subjected to uniform normal contact load was investigated by finite element method. Results show that the surface tensile stress is mainly affected by the thickness ratio of Al layer when the total thickness of coating is unchanged. With the increase of A] layer thickness, the surface tensile stress rises quickly. When Al2O3 layer thickness increases, surface tensile stress is diminished. 'Meanwhile, the maximum shear stress moves rapidly towards internal part of HDA/PEO coatings. Shear stress at the Al2O3/Al interface is minimal when Al2O3 layer and Al layer have the same thickness.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.

Atomistic investigation of the effects of temperature and surface roughness on diffusion bonding between Cu and Al

Molecular dynamics (MD) simulations are carried out to analyze the diffusion bonding at Cu/Al interfaces. The results indicate that the thickness of the interfacial layer is temperature-dependent, with higher temperatures yielding larger thicknesses. At temperatures below 750 K, the interface thickness is found to increase in a stepwise manner as a function of time. At temperatures above 750 K, the thickness increases rapidly and smoothly. When surface roughness is present, the bonding process consists of three stages. In the first stage, surfaces deform under stress, resulting in increased contact areas. The second stage involves significant plastic deformation at the interface as temperature increases, resulting in the disappearance of interstices and full contact of the surface pair. The last stage entails the diffusion of atoms under constant temperature. The bonded specimens show tensile strengths reaching 88% of the ideal Cu/Al contact strength. (c) 2007 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Deflections and curvatures of a film-substrate structure with the presence of gradient stress in MEMS applications

The close form solutions of deflections and curvatures for a film-substrate composite structure with the presence of gradient stress are derived. With the definition of more precise kinematic assumption, the effect of axial loading due to residual gradient stress is incorporated in the governing equation. The curvature of film-substrate with the presence of gradient stress is shown to be nonuniform when the axial loading is nonzero. When the axial loading is zero, the curvature expressions of some structures derived in this paper recover the previous ones which assume the uniform curvature. Because residual gradient stress results in both moment and axial loading inside the film-substrate composite structure, measuring both the deflection and curvature is proposed as a safe way to uniquely determine the residual stress state inside a film-substrate composite structure with the presence of gradient stress.

Thermal stress wave and spallation induced by an electron beam

Thermal stress wave and spallation in aluminium alloy exposed to a high fluency and low energy electron beams are studied theoretically. A simple model for the study of energy deposition of electrons in materials is presented on the basis of some empirical formulae. Under the stress wave induced by energy deposition, microcracks and/or microvoids may appear in target materials, and in this case, the inelastic volume deformation should not vanish. The viscoplastic model proposed by Bodner and Partom with corresponding Gurson's yield function requires modification for this situation. The new constitutive model contains a scalar field variable description of the material damage which is taken as the void volume fraction of the polycrystalline material. Incorporation of the damage parameter permits description of rate-dependent, compressible, inelastic deformation and ductile fracture. The melting phenomenon has been observed in the experiment, therefore one needs to take into account the melting process in the intermediate energy deposition range. A three-phase equation of state used in the paper provides a more detailed and thermodynamical description of metals, particularly, in the melting region. The computational results based on the suggested model are compared with the experimental test for aluminium alloy, which is subjected to a pulsed electron beam with high fluency and low energy. (C) 1997 Elsevier Science Ltd.

Micromechanics analysis of particulate-reinforced composites and their failure mechanisms

Stress fields and failure mechanisms have been investigated in composites with particles either surface treated or untreated under uniaxial tension. Previous experimental observation of failure mechanisms in a composite with untreated particles showed that tensile cracks occurred mostly at the polar region of the particle and grew into interfacial debonding. In a composite with surface-treated particles, however, shear yielding and shear cracking proceeded along the interphase-matrix interface at the polar area of the matrix and thus may improve the mechanical behaviour of the material. The finite element calculations showed that octahedral shear stress at the polar and longitudinal areas of the particle treated by coupling agents is much larger than that of materials with untreated particles, and the shear stress distribution around the interface is sensitive to the interphase property. The results suggest that a th ree-phase model can describe the composites with surface-treated fillers.

Micromechanics analyses of interface crack

A new mechanics model based on Peierls concept is presented in this paper, which can clearly characterize the intrinsic features near a tip of an interfacial crack. The stress and displacement fields are calculated under general combined tensile and shear loadings. The near tip stress fields show some oscillatory behaviors but without any singularity and the crack faces open completely without any overlapping when remote tensile loading is comparable with remote shear loading. A fracture criterion for predicting interface toughness has been also proposed, which takes into account for the shielding effects of emitted dislocations. The theoretical toughness curve gives excellent prediction, as compared with the existing experiment data.

Three-dimensional finite element analysis of composites with coated spherical inclusions

A three-dimensional finite element analysis has been used to determine the internal stresses in a three-phase composite. The stresses have been determined for a variety of interphase properties, the thicknesses of the interphase and the volume fractions of particles. Young's modulus has been calculated from a knowledge of these stresses and the applied deformation. The calculations show that stress distributions in the matrix and the mechanical properties are sensitive to the interphase property in the three-phase composites. The interfacial stresses in the three-dimensional analysis are in agreement with results obtained by an axisymmetric analysis. The predicted bulk modulus in three-dimensional analysis agrees well with the theoretical solution obtained by Qui and Weng, but it presents a great divergence from that in axisymmetric analyses. An investigation indicates that this divergence may be caused by the difference in the unit cell structure between two models. A comparison of the numerically predicted bulk and shear modulus for two-phase composites with the theoretical results indicates that the three-dimensional analysis gives quite satisfactory results.