Derivation Of A Nonlinear Reynolds Stress Model Using Renormalization Group Analysis And Two-Scale Expansion Technique

Adopting Yoshizawa's two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

Interface Layer Effect On The Stress Distribution Of A Wafer-Bonded Bilayer Structure

The interface layer plays an important role in stress transfer in composite structures. However, many interface layer properties such as the modulus, thickness, and uniformity are difficult to determine. The model developed in this article links the influence of the interface layer on the normal stress distribution along the layer thickness with the layer surface morphology before bonding. By doing so, a new method of determining the interfacial parameter(s) is suggested. The effects of the layer thickness and the surface roughness before bonding on the normal stress distribution and its depth profile are also discussed. For ideal interface case with no interfacial shear stress, the normal stress distribution pattern can only be monotonically decreased from the interface. Due to the presence of interfacial shear stress, the normal stress distribution is much more complex, and varies dramatically with changes in the properties of the interface layer, or the dimensions of the bonding layers. The consequence of this dramatic stress field change, such as the shift of the maximum stress from the interface is also addressed. The size-dependent stress distribution in the thickness direction due to the interface layer effect is presented. When the interfacial shear stress is reduced to zero, the model presented in this article is also demonstrated to have the same normal stress distribution as obtained by the previous model, which does not consider the interface layer effect.

Stress intensity factor histories of a steadily propagating crack under 3-D combined mode traction

Elastodynamic stress intensity factor histories of an unbounded solid containing a semi-infinite plane crack that propagates at a constant velocity under 3-D time-independent combined mode loading are considered. The fundamental solution, which is the response of point loading, is obtained. Then, stress intensity factor histories of a general loading system are written out in terms of superposition integrals. The methods used here are the Laplace transform methods in conjunction with the Wiener-Hopf technique.

A second-order dynamic subgrid-scale stress model

A second-order dynamic model based on the general relation between the subgrid-scale stress and the velocity gradient tensors was proposed. A priori test of the second-order model was made using moderate resolution direct numerical simulation date at high Reynolds number ( Taylor microscale Reynolds number R-lambda = 102 similar to 216) for homogeneous, isotropic forced flow, decaying flow, and homogeneous rotating flow. Numerical testing shows that the second-order dynamic model significantly improves the correlation coefficient when compared to the first-order dynamic models.

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.

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.

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.

Elastodynamic stress-intensity factors for a semiinfinite crack under 3-d combined mode loading

The dynamic stress intensity factor histories for a half plane crack in an otherwise unbounded elastic body are analyzed. The crack is subjected to a traction distribution consisting of two pairs of suddenly-applied shear point loads, at a distance L away from the crack tip. The exact expression for the combined mode stress intensity factors as the function of time and position along the crack edge is obtained. The method of solution is based on the direct application of integral transforms together with the Wiener-Hopf technique and the Cagniard-de Hoop method, which were previously believed to be inappropriate. Some features of solutions are discussed and the results are displayed in several figures.

Elastodynamic stress intensity factor history for a semi-infinite crack under three-dimensional transient loading

The dynamic stress intensity factor history for a semi-infinite crack in an otherwise unbounded elastic body is analyzed. The crack is subjected to a pair of suddenly-applied point loadings on its faces at a distance L away from the crack tip. The exact expression for the mode I stress intensity factor as a function of time is obtained. The method of solution is based on the direct application of integral transforms, the Wiener-Hopf technique and the Cagniard-de Hoop method. Due to the existence of the characteristic length in loading this problem was long believed a knotty problem. Some features of the solutions are discussed and graphical result for numerical computation is presented.