A novel singular finite element of mixed-mode crack problems with arbitrary crack tractions

A novel singular finite element is presented to study cracked plates with arbitrary traction acting on crack surfaces. Firstly, the analytical solution around crack tips is determined using the symplectic dual approach. Subsequently, the solution is used to develop a novel singular finite element, which depicts accurately the characteristic of singular stresses field near crack tips. And the novel element can be applied to solve cracked plates, and both Mode I and Mode II stress intensity factors can be determined directly and accurately. Lastly, two numerical examples are given to illustrate the present method.     

Using M-integral for multi-cracked problems subjected to nonconservative and nonuniform crack surface tractions

In this paper, an energy parameter based on the concept of the M-integral is proposed for describing the fracture behavior of a multi-cracked solid subjected to nonconservative and nonuniform crack surface tractions. By using the M-integral with a suitably chosen closed contour, one can evaluate the ‘surface creation energy’ (SCE) required for creation of the stressed cracks. Also, it is demonstrated that the property of path-independence holds even under the action of crack surface tractions. Therefore, the singular stress field in the near-tip areas is not directly involved in the calculation so that a complicated finite element model around the crack tips is not required in evaluation of the M-integral.     

Stress intensity factors of an oblique edge crack subjected to normal and shear tractions

Stress intensity factors (SIFs) were obtained for an oblique crack under normal and shear traction and remote extension loads. The oblique crack was modeled as the pseudodislocation. The stress field due to tractions was solved by the Flamant solution. The SIR of Mode I and Mode II (K_Iand K_II) were then obtained. Finite element analysis was performed with ABAQUS and compared with the analytical solutions. The analytical solutions were in good agreement with the results of FEM. From investigating SIFs and their ranges, the following results were obtained. The growth rate of an oblique edge crack decreased due to the reduction in the SIF ranges. The crack driving force depended on the obliquity, the normal traction and the ratio of crack to traction length. The peak value of shear traction was found as a key parameter to accelerate the crack growth.     

Analysis of multiple cracks in an infinite plate under arbitrary crack surface tractions

Summary Based on the analytical solution derived for a crack under arbitrary crack surface tractions an alternating method is presented to analyze two-dimensional mixed-mode fracture problems with multiple cracks which are arbitrarily distributed in an infinite plate. During the iterative solution process the reversed residual stresses at the location of fictitious cracks induced by other cracks are repeatedly evaluated and treated as new applied loading until the residual stresses on each crack surface are negligible. The influence of interaction effects among cracks on the determination of stresses and deformations is studied in detail. Several numerical example are solved to demonstrate the validity of the method. Good agreement between the computed stress intensity factors and referenced solutions shows the accuracy and efficiency of the this work.

---

Analyse mehrfacher risse in einer unendlichen platte mit beliebigen randspannungen am riß

Übersicht Eine alternierende Methode, die auf dem Ergebnis für einen Riß unter beliebigen Spannungen an der Rißoberfläche basiert, wird hier zur Analyse zweidimensionaler mixed-mode Bruchprobleme mit mehrfachen Rissen präsentiert. Die mehrfachen Risse können in der Platte beliebig verteilt sein. Mit einem Iterationsprozeß wird die umgekehrte Restspannung an der Stelle angenommener Risse, die von anderen Rissen verursacht werden, solange berechnet und als neue Last betrachtet, bis die Restspannung in jeder Rißfläche sehr klein geworden ist. Bei der Bestimmung der Spannung und der Verformung wird der wechselseitige Einfluß zwischen den Rissen eingehend untersucht. Um die Genauigkeit der Methode zu prüfen, werden mehrere Beispiele gegeben. Die Ergebnisse zeigen, daß die gerechnete Spannungsintensitätsfaktoren mit bekannten Lösungen sehr gut übereinstimmen.

     

Stochastic framework for modeling the linear apparent behavior of complex materials: Application to random porous materials with interphases

This paper is concerned with the modeling of randomness in multiscale analysis of heterogeneous materials. More specifically, a framework dedicated to the stochastic modeling of random properties is first introduced. A probabilistic model for matrix-valued second-order random fields with symmetry propertries, recently proposed in the literature, is further reviewed. Algorithms adapted to the Monte Carlo simulation of the proposed representation are also provided. The derivations and calibration procedure are finally exemplified through the modeling of the apparent properties associated with an elastic porous microstructure containing stochastic interphases.     

Simulation of crack spacing using a reinforced concrete model with an internal length parameter

Summary A gradient-enhanced smeared crack model and bond-slip interface elements are utilized in finite element simulations of reinforced concrete. The crack model is rooted in an enhanced plasticity theory. It uses the Rankine failure surface dependent on an equivalent inelastic strain measure as well as on its Laplacian. As a result, finitely sized fracture process zones and realistic crack spacings are obtained. A reinforced concrete bar in uniaxial tension is analyzed to demonstrate the regularizing influence of the internal length parameter in the model and to evaluate the influence of the model parameters on the energy dissipation in multiple cracks. A comparison of numerical simulations with experimental results for a beam without shear reinforcement in four-point bending concludes the analysis. Received 4 November 1997; accepted for publication 23 April 1998     

Transient analysis of a propagating crack with finite length subjected to a horizontally polarized shear wave

In this study, the transient response of a finite crack subjected to an incident horizontally polarized shear wave and then propagated with a constant speed in an unbounded elastic solid is investigated. Initially, the finite crack with crack length l is stress-free and at rest. At time t = 0, an incident horizontally polarized shear wave strikes at one of the crack tips and will arrive at the other tip at a later time. Then, two crack tips propagate along the crack tip line with different velocities as the corresponding stress intensity factors reach their fracture toughness. The correspondent configuration is shown in Fig. 1

1. Download : Download high-res image (31KB)
2. Download : Download full-size image

Fig. 1. Configuration and coordinate systems of a finite crack in an unbounded medium.

. In analyzing this problem, diffracted waves generated by two propagating crack tips must be taken into account and it makes the analysis extremely difficult. In order to solve this problem, the transform formula in the Laplace transform domain between moving and stationary coordinates is first established. Complete solutions are determined by superposition of proposed fundamental solutions in the Laplace transform domain. The fundamental solutions to be used are from the problems of applying exponentially distributed traction and screw dislocation on crack faces and along the crack tip line, respectively. The exact transient solutions of dynamic stress intensity factor for the first few diffracted waves that arrive at two crack tips are obtained and expressed in compact formulations. Numerical calculations of dynamic stress intensity factors for both tips are evaluated and the results are discussed in detail.     

Fracture behavior of a bonded magneto-electro-elastic rectangular plate with an interface crack

In this paper the anti-plane problem for an interface crack between two dissimilar magneto-electro-elastic plates subjected to anti-plane mechanical and in-plane magneto-electrical loads is investigated. The interface crack is assumed to be either magneto-electrically impermeable or permeable, and the position of the interface crack is arbitrary. The finite Fourier transform method is employed to reduce the mixed boundary-value problem to triple trigonometric series equations. The dislocation density functions and proper replacement of the variables are introduced to reduce these series equations to a standard Cauchy singular integral equation of the first kind. The resulting integral equation together with the corresponding single-valued condition is approximated as a system of linear algebra equations which can be easily solved. Field intensity factors and energy release rates are determined numerically and discussed in detail. Numerical results show the effects of crack configuration and loading combination parameters on the fracture behaviors of crack tips according to energy release rate criterion. The study of this problem is expected to have applications to the investigation of dynamic fracture properties of magneto-electro-elastic materials with cracks.     

Stochastic stability of parametrically excited random systems

Multidegree-of-freedom dynamic systems subjected to parametric excitation are analyzed for stochastic stability. The variation of excitation intensity with time is described by the sum of a harmonic function and a stationary random process. The stability boundaries are determined by the stochastic averaging method. The effect of random parametric excitation on the stability of trivial solutions of systems of differential equations for the moments of phase variables is studied. It is assumed that the frequency of harmonic component falls within the region of combination resonances. Stability conditions for the first and second moments are obtained. It turns out that additional parametric excitation may have a stabilizing or destabilizing effect, depending on the values of certain parameters of random excitation. As an example, the stability of a beam in plane bending is analyzed.Published in Prikladnaya Mekhanika, Vol. 40, No. 10, pp. 135–144, October 2004.     

Stochastic non-linear flutter of a panel subjected to random in-plane forces

—An analysis of non-linear flutter of a simply-supported panel exposed to supersonic gas flow and random in-plane forces is presented for two- and three-mode interactions. A first order quasi-steady state aerodynamic piston theory is used to model the aerodynamic loading. The Fokker-Planck equation is used to derive a general moment equation for two- and three-mode interactions. For stability analysis the moment equation is consistent and the mean square stability boundaries of the equilibrium are obtained in terms of the system parameters. The stability boundaries reveal common features to those predicted by the deterministic theory of panel nutter. For the non-linear response the moment equation is found inconsistent and a cumulant-neglect closure is used by setting cumulants of fifth and sixth orders to zero. This first order non-Gaussian closure is carried out to solve for the response statistics in terms of the air-to-plate mass ratio, aerodynamic pressure, modal damping, and in-plane random force spectral density. It is found that the non-Gaussian solution yields higher levels for the response statistics than those obtained by the Gaussian solution. The inclusion of more modes results in a reduction of the response levels and expands the stability region.     

Stability,bifurcation and post-critical behavior of a homogeneously deformed incompressible isotropic elastic parallelepiped subject to dead-load surface tractions

We study the equilibrium homogeneous deformations of a homogeneous parallelepiped made of an arbitrary incompressible, isotropic elastic material and subject to a distribution of dead-load surface tractions corresponding to an equibiaxial tensile stress state accompanied by an orthogonal uniaxial compression of the same amount. We show that only two classes of homogeneous equilibrium solutions are possible, namely symmetric deformations, characterized by two equal principal stretches, and asymmetric deformations, with all different principal stretches. Following the classical energy-stability criterion, we then find necessary and sufficient conditions for both symmetric and asymmetric equilibrium deformations to be weak relative minimizers of the total potential energy. Finally, we analyze the mechanical response of a parallelepiped made of an incompressible Mooney–Rivlin material in a monotonic dead loading process starting from the unloaded state. As a major result, we model the actual occurrence of a bifurcation from a primary branch of locally stable symmetric deformations to a secondary, post-critical branch of locally stable asymmetric solutions.     

Elastic-plastic behavior of elastically homogeneous materials with a random field of inclusions

A multiphase material is considered, which consists of a homogeneous elastic-plastic matrix containing a homogeneous statistically uniform random set of ellipsoidal elastic-plastic inclusions. The elastic properties of the matrix and the inclusions are the same, but the so-called “stress-free strains”, i.e. the strain contributions due to temperature loading, phase transformations, and the plastic strains, fluctuate. A general theory of the yielding for arbitrary loading (by the stress and by the temperature) is employed. The realization of an incremental plasticity scheme is based on averaging over each component of the nonlinear yield criterion. Usually, averaged stresses are used inside each component for this purpose. In distinction from this usual practice physically consistent assumptions about the dependence of these functions on the component's values of the second stress moments are applied. The application of the proposed theory to the prediction of the thermomechanical deformation behavior of a model material is shown.     

Elastostatic crack tip behavior for a rubber-like material

Analyzed in this work is the elastostatic field near a crack tip in a rubber-like material. Asymptotic equations for a crack opened symmetrically about its plane are derived from assumed forms of the strain energy density and constitutive relation that applies to large and finite strain and remains valid even when the strain tends to infinity in the limit. Near field solutions are obtained in regions that decreases and increases in size as the crack tip is approached. Their singular character depends on the constitutive parameters and is evaluated numerically.     

Asymptotic behavior of the energy functional for a three-dimensional body with a rigid inclusion and a crack

A model of a three-dimensional elastic body containing a rigid inclusion and a crack located on the interface between the inclusion and the body is considered. Natural boundary conditions are imposed on the crack. A derivative of the energy functional with respect to the perturbation parameter is derived for an arbitrary, rather smooth perturbation of the domain, in particular, the Griffith formula is obtained.     

Stochastic thermoelastic problem of a functionally graded plate under random temperature load

This study attempts to derive the statistics of temperature and thermal stress in functionally graded material (FGM) plates exposed to random external temperatures. The thermomechanical properties of the FGM plates are assumed to vary arbitrarily only in the plate thickness direction. The external temperatures are expressed as random functions with respect to time. The transient temperature field in the FGM plate is determined by solving a nonhomogeneous heat conduction problem for a multilayered plate with linear nonhomogeneous thermal conductivity and different homogeneous heat capacity in each layer. The autocorrelations and power spectrum densities (PSDs) of temperature and thermal stress are derived analytically. These statistics for FGM plates composed of partially stabilised zirconia (PSZ) and austenitic stainless steel (SUS304) are computed under the condition that the fluctuation in the external temperature can be considered as white noise or a stationary Markov process.     

Stochastic model reduction for robust dynamical characterization of structures with random parameters

In this paper, we characterize random eigenspaces with a non-intrusive method based on the decoupling of random eigenvalues from their corresponding random eigenvectors. This method allows us to estimate the first statistical moments of the random eigenvalues of the system with a reduced number of deterministic finite element computations. The originality of this work is to adapt the method used to estimate each random eigenvalue depending on a global accuracy requirement. This allows us to ensure a minimal computational cost. The stochastic model of the structure is thus reduced by exploiting specific properties of random eigenvectors associated with the random eigenfrequencies being sought. An indicator with no additional computation cost is proposed to identify when the method needs to be enhanced. Finally, a simple three-beam frame and an industrial structure illustrate the proposed approach.     

Stochastic modelling of crack growth based on damage accumulation

The growth rate of a fatigue crack is modelled from a damage accumulation standpoint. The material ahead of the crack-tip is considered to be composed of assembly of uniaxial fatigue elements which accumulate damage per load cycle. Each element is subjected to increased levels of stress and strain ranges as the crack propagates. A linearly accumulated damage criterion is assumed, and failure of an element indicates a void initiation at its position. Both deterministic and stochastic analyses are included. The historical damage of the material before it reaches the crack tip vicinity is quantified and is shown to be significant for the first few elements. The predicted results agree fairly well with the experimental data.     

非线性呼吸行为对带横向裂纹Jeffcott转子稳定性的影响

在分析现有力学模型局限性的基础上,考虑转子在失稳转速附近复杂的涡动状态,构建了不依赖于重力占优假设的含横向裂纹的Jeffcott转子的通用运动方程和一种拟合的裂纹呼吸模型。根据Floquet稳定性理论,用打靶法分析了该裂纹转子模型运动的稳定性。结果表明,刚度相对变化量、系统阻尼、裂纹相位角等系统参数对模型的运动稳定性都有较大影响,特别是裂纹相位角对稳定性的影响很复杂;在临界转速附近,裂纹的非线性呼吸行为可以提高转子的稳定性,且改变了失稳发散的行为;在亚临界转速附近,裂纹的非线性呼吸行为会降低转子的稳定性,使周期运动发生结构失稳。     

Determination of the length of a short crack at a v-notch from a full field measurement

Determination of the length of a short crack at the root of a v-notch, from a full kinematic field measurement, is performed using a direct method. It is based on a matched asymptotic expansions procedure together with the theory of singularities. The first corrective term of the outer expansion can be straightforwardly expressed as a function of the crack length. Its extraction is achieved through the calculation of the associated generalized stress intensity factors for elastic homogeneous materials as well as bimaterials. Numerical simulations are carried out on a finite element solution disturbed by a random noise. In addition, the method used to compute the generalized stress intensity factors proved accurate and robust.     

Stochastic optimal active control of a half-car nonlinear suspension under random road excitation

The random response analysis and the stochastic optimal active control of a half-car model with nonlinear suspension stiffness and damping traversing a rough road are studied in this paper. The road roughness height is modeled as the output of a first-order linear filter to Gaussian white noise. Considering the hysteretic nonlinear stiffness and the square damping of the vehicle model, the response statistics of the nonlinear suspension with active control are obtained by using the equivalent linearization method. The performance indexes of the active suspension are evaluated and compared with those of the corresponding passive suspension. It is found that the nonlinear active suspension gives a better vehicle performances like ride comfort, suspension stroke and overall performance. Finally, the theoretical results are verified through Monte Carlo simulation.