Numerical modelling of the scattering of elastic waves in plates

This paper describes the application of finite difference methods to the calculation of the scattering of elastic waves. The emphasis is on cracklike defects in plates, and it is shown that a common numerical technique can span a range of wavelengths from Lamb waves to ultrasonic waves with many reflections from the surfaces of the plate. Quantitative results are given for the scattering of Lamb waves and ultrasonic shear waves from surface-breaking cracks.     

Automated dynamic fracture procedure for modelling mixed-mode crack propagation using explicit time integration brick finite elements

Several fracture codes have been developed in recent years to perform analyses of dynamic crack propagation in arbitrary directions. However, general-purpose, commercial finite-element software which have capabilities to do fracture analyses are still limited in their use to stationary cracks and crack propagation along trajectories known a priori . In this paper, we present an automated fracture procedure implemented in the large-scale, nonlinear, explicit, finite-element code DYNA3D which can be used to simulate dynamic crack propagation in arbitrary directions. The model can be used to perform both generation- and application-phase simulations of self-similar as well as non-self-similar dynamic crack propagation in linear elastic structures without user intervention. It is developed based on dynamic fracture mechanics concepts and implemented for three-dimensional solid elements. Energy approach is used in the model to check for crack initiation/propagation. Dynamic energy release rate and stress intensity factors are determined from far-field finite-element field solutions using finite-domain integrals. Fracture toughness is input as a function of crack-tip velocity, and when the criterion for crack growth is satisfied, an element deletion-and-replacement re-meshing procedure is used along with a gradual nodal release technique to update the crack geometry and model the crack propagation. Direction of crack propagation is determined using the maximum circumferential stress criterion. Numerical simulations of experiments involving non-self-similar crack propagation are performed, and results are presented as verification examples.     

Boundary integral equation analysis for radiation and scattering of elastic waves in three dimensions

Abstract

In this paper, the boundary integral equation (BIE) method is employed to investigate the radiation and scattering of time‐harmonic elastic waves by obstacles of arbitrary shape embedded in an infinite medium. Based on the vector BIE, entirely free of Cauchy principal value integrals, an efficient numerical scheme using quadratic isoparametric boundary elements is proposed. Furthermore, the difficulty of non‐uniquess of a solution inherent with BIE formulations for exterior elastodynamic problems is studied numerically and analytically. The counterparts of the combined Helmholtz integral formulation method for elastodynamics together with the least‐square or Lagrange‐multiplier technique are derived and applied to overcome this difficulty successfully. In addition, the elastic‐wave fields radiated or scattered by either a spherical cavity or a rigid sphere in an infinite medium are calculated and the results are compared with the analytical solutions to demonstrate the accuracy and versatility of the proposed numerical scheme.     

Numerical modelling of fatigue crack initiation and growth of martensitic steels

This paper presents a numerical simulation of micro‐crack initiation that is based on Tanaka‐Mura micro‐crack nucleation model. Three improvements were added to this model. First, multiple slip bands where micro‐cracks may occur are used in each grain. Second improvement deals with micro‐crack coalescence by extending existing micro‐cracks along grain boundaries and connecting them into a macro‐crack. The third improvement handles segmented micro‐crack generation, where a micro‐crack is not nucleated in one step like in Tanaka‐Mura model, but is instead generated in multiple steps. High cycle fatigue testing was also performed and showed reasonably good correlation of proposed model to experimental results. Because numerical model was directed at simulating fatigue properties of thermally cut steel, edge properties of test specimens were additionally inspected in terms of surface roughness and micro‐structural properties.     

Numerical modelling of micro-pitting of gear teeth flanks

A two‐dimensional computational model for simulation of contact fatigue of gear teeth flanks is presented. In the model, it is assumed that the initial crack of length 0.015 mm is initiated at the surface due to previous mechanical or heat treatment of the material as well as a consequence of the running in process. The discretized model with the initial crack is then subjected to normal contact pressure, which takes into account the elasto‐hydro‐dynamic (EHD) lubrication conditions, and tangential loading due to friction between contacting surfaces. The model also considers the moving contact of gear flanks, fluid trapped in the crack and residual stresses due to heat treatment of the material on crack propagation. The virtual crack extension (VCE) method, implemented in the finite element method, is then used for simulating the fatigue crack growth from the initial crack up to the formation of the surface pit. The computational results show that the initial surface crack of length 15 μm and the considered boundary conditions lead to the appearance of very small surface pits, which can be termed as micro‐pitting on gear teeth flanks. The numerical results correspond well with available experimental data.     

Numerical resolution of the boundary integral equations for elastic scattering by a plane crack

The problem of wave scattering by a plane crack is solved, either in the case of acoustic waves or in the case of elastic waves incidence using the boundary integral equation method. A collocation method is often used to solve that equation, but here we will use a variational method, first writing the problem of Fourier variables, and then writing the associated integrals in the sesquilinear form with weak singularity kernels. This representation is used in the numerical approach, made with a finite element method in the surface of the crack. Numerical tests were made with circular and elliptical cracks, but this method can be extended to other shapes, with the same convergence profiles. Extensive results are given concerning the crack opening displacement, the scattering cross-section, the back-scattered amplitude and far-field patterns.     

A generalized approach to the mechanics of three-dimensional nonlinear ocean waves

The paper suggests a generalized approach to the analysis of second-order three-dimensional random wave fields that can represent either ocean waves in an undisturbed field or standing ocean waves. A nonlinear stochastic family is introduced that has a linear stationary Gaussian random contribution and a second-order correction exact in a Stokes expansion. This family can describe the nonlinear mechanics of a random wave field. Subsequently, in the considered random ocean field, the nonlinear mechanics of wave groups can be investigated conditional to the occurrence of a very high wave. The free-surface displacement and velocity potential may be obtained, as well as the nonlinear wave velocity and acceleration, for both the stochastic family and the very high wave groups. The results are of interest for the analysis of nonlinear wave kinematics when extreme surface wave events occur and thus for Morison′s wave forces on slender structures.     

Elastoplastic modelling of subsurface crack growth in rail/wheel contact problems

Propagation of small subsurface cracks subjected to shear under repeated rolling contact load is studied. An analytical crack model (Dugdale) with plastic strips at the two crack tips is employed. Compressive stresses promoting crack closure and friction between crack faces are considered. The triaxial stress state is used in the yield criterion. A damage criterion is suggested based on experimental LCF data. In a numerical study, critical crack lengths are found below which propagation of an existing crack should be effectively suppressed.     

Linearized inversion methods for three-dimensional electromagnetic imaging in the multiple scattering regime

In optical or microwave computational tomography, the sample permittivity is reconstructed numerically from the measurements of its scattered field for various illuminations. When the light sample interaction involves multiple scattering, the relationship between the scattered field and the permittivity is non-linear and a direct reconstruction is not possible. Using a simple physical approach, adapted to the three-dimensional vectorial electromagnetic framework, we derive an iterative inversion technique, based on the linearization of the scattering operator, for imaging (possibly anisotropic) targets in the multiple scattering regime. We investigate the performances of different approximations of this operator accounting for more or less multiple scattering. Our method is applied to the reconstruction of targets in the microwave domain using experimental data.     

含裂纹悬臂梁的振动与疲劳耦合分析

基于Paris方程和同步分析方法考虑振动与疲劳裂纹扩展耦合之影响,提出一种含裂纹梁的振动疲劳寿命分析思路.振动分析过程中,利用线性弹簧等效裂纹段,复弹性模量引入阻尼损耗因子,得到考虑裂纹扩展、激励频率和阻尼等因素影响的动应力响应.结果表明:裂纹扩展、激励频率和阻尼等因素对疲劳寿命具有重要的影响.通过振动分析与疲劳裂纹扩展寿命估算同步进行,可进一步提高疲劳寿命估算精度.     

SV波入射下多重散射问题的改进算法以及双排非连续刚性屏障的隔振研究

摘 要：相比起Aviles等学者的单排非连续刚性屏障隔振理论,采用更为完整的Fourier-Bessel函数系的级数展开表达式来研究弹性波的多重散射问题,结合位移连续的边界条件,并引入更具一般意义上的Graf加法定理,推导了无限均匀介质中任意分布的多个刚性圆柱体对于平面SV波入射时散射系数的理论解析解。随后,取所有的圆柱体位于相平行的两条直线上,则问题演变成为双排刚性桩屏障对于平面SV波的隔离问题。重点讨论了两排桩之间的间距h对屏障整体隔振效果的影响及其他的一些隔振性状。特别地,当h=0时,该问题退化为常见的单排非连续屏障的隔振问题。     

Numerical modelling of elastic wave scattering in frequency domain by the partition of unity finite element method

In this paper, we investigate a numerical approach based on the partition of unity finite element method, for the time‐harmonic elastic wave equations. The aim of the proposed work is to accurately model two‐dimensional elastic wave problems with fewer elements, capable of containing many wavelengths per nodal spacing, and without refining the mesh at each frequency. The approximation of the displacement field is performed via the standard finite element shape functions, enriched by superimposing pressure and shear plane wave basis, which incorporate knowledge of the wave propagation. A variational framework able to handle mixed boundary conditions is described. Numerical examples dealing with the radiation and the scattering of elastic waves by a circular body are presented. The results show the performance of the proposed method in both accuracy and efficiency. Copyright © 2008 John Wiley & Sons, Ltd.     

Dual boundary element method and finite element method for mixed‐mode crack propagation simulations in a cracked hollow shaft

Three‐dimensional mixed‐mode crack propagation simulations were performed by means of the dual boundary element method code BEASY and 2 finite element method‐based crack propagation codes: ZENCRACK (ZC) and CRACKTRACER3D (CT3D). The stress intensity factors (SIFs) along the front of an initial semielliptical crack, initiated from the external surface of a shaft, were calculated for 4 different load cases: bending, press fit, shear, and torsion. The methods used for the SIF assessment along the crack front were the J‐integral for BEASY and ZC and the quarter point element stress method for CT3D. Subsequently, crack propagation simulations were performed, with the crack growth rate evaluated by using Paris' law, calibrated for the material at stake (American Society for Testing and Materials A469 steel). The kink angles were evaluated by using the minimum strain energy density and maximum tangential stress criteria for BEASY, the maximum energy release rate and maximum tangential stress for ZC, and the maximum principal asymptotic stress for CT3D. The results obtained in terms of SIFs and crack propagation life show very good agreement among the 3 codes. Also, the shape of the propagated crack, which is significantly out‐of‐plane for the shear and torsion loading, matched very well.     

Measurement and modelling of residual stress effects on cracks

Defects that form by mechanisms such as fatigue and stress corrosion cracking are influenced both by external loads on engineering structures and internal, residual stresses that are generated during the manufacture and operation. This paper describes a programme of experimental and analytical work on a high‐strength, low‐toughness aluminium alloy (AL2024‐T351) to assess the influence of residual stress on crack opening displacement (COD) and crack‐driving force (CDF) for a range of fatigue crack lengths in compact tension (CT) specimens containing a mechanically induced residual stress field. Comparison of experimentally measured and numerically predicted CODs, at the mid‐plane and surface of CT specimens, show generally good agreement for cracks introduced into the finite‐element model in a progressive, element‐by‐element manner. Cracks introduced in a simultaneous manner give larger than observed CODs. The CDFs for the progressively introduced crack are always smaller than for simultaneously introduced. These results have implications for the assessment of initiation for slowly growing cracks.     

Dynamic stress intensity factors due to scattering of harmonic waves by a crack in an infinite medium

An analytical method for calculating dynamic stress intensity factors in the mixed mode (combination of opening and sliding modes) using complex functions theory is presented. The crack is in infinite medium and subjected to the plane harmonic waves. The basis of the method is grounded on solving the two‐dimensional wave equations in the frequency domain and complex plane using mapping technique. In this domain, solution of the resulting partial differential equations is found in the series of the Hankel functions with unknown coefficients. Applying the boundary conditions of the crack, these coefficients are calculated. After solving the wave equations, the stress and displacement fields, also the J‐integrals are obtained. Finally using the J‐integrals, dynamic stress intensity factors are calculated. Numerical results including the values of dynamic stress intensity factors for a crack in an infinite medium subjected to the dilatation and shear harmonic waves are presented.     

Predicting the crack opening displacements of a partially debonded pipeline in rock mass under P waves

A crack open may lead to the failure because of rupture or leak in the pipeline. In this paper, the dynamic crack opening displacements (CODs) of a partially debonded pipeline subjected to P and SV waves are predicted, and a mathematical model is proposed. The debonded region around the pipeline is modelled as an interface crack with non‐contacting faces. The wave fields in the bonded and debonded areas are expressed by wave function expansion method, and the expanded coefficients are determined by satisfying the boundary conditions with consideration of cracks. The CODs in the debonding area are represented by Chebyshev polynomials. Some examples are illustrated to analyse the CODs under different thicknesses of pipeline and different wave frequencies. In the region of high frequency, the value of CODs is more sensitive to the crack size. The difference between the CODs in the compression direction and shear directions is also examined. Comparison with existing results is given to validate this dynamic model.     

Finite element modelling of lamb waves propagation in 2D plates and thin sheets for damage detection

Structural health monitoring is an evolving technology applied to identify, locate and quantify severity of damages in structures before failure. Lamb waves have become a keen interest for inspection since they can be used to monitor a large area from one single location. The objective of this research is to simulate lamb wave response using finite element method and its application to crack detection and identification in thin metallic engineering structures. Two types of specimen i. e. two‐dimensional aluminium plate and thin aluminium sheets are simulated using commercially available finite element package ABAQUS. Initially phase velocity and group velocity dispersion curve are plotted for aluminium material. Thereafter simulation of individual specimens with cracks and without cracks is performed. Simulation results were compared and validated with actual results and were found to be in reasonably good agreement. This is certainly done by employing group velocity and time of flight for the distance travelled between the monitoring point and crack position. Assessment of lamb waves sensitivity to various sizes and shapes of cracks like rectangular and circular are also investigated and its effect on the structure is discussed in detail.     

An analytical method for mixed-mode propagation of pressurized fractures in remotely compressed rocks

An analytical method for mixed-mode (mode I and mode II) propagation of pressurized fractures in remotely compressed rocks is presented in this paper. Stress intensity factors for such fractured rocks subjected to two-dimensional stress system are formulated approximately. A sequential crack tip propagation algorithm is developed in conjunction with the maximum tensile stress criterion for crack extension. For updating stress intensity factors during crack tip propagation, a dynamic fictitious fracture plane is used. Based on the displacement correlation technique, which is usually used in boundary element/finite element analyses, for computing stress intensity factors in terms of nodal displacements, further simplification in the estimation of crack opening and sliding displacements is suggested. The proposed method is verified comparing results (stress intensity factors, propagation paths and crack opening and sliding displacements) with that obtained from a boundary element based program and available in literatures. Results are found in good agreements for all the verification examples, while the proposed method requires a trivial computing time.     

Simulation of crack growth in composite material shell structures

This paper implements a domain integral energy method for modelling crack growth in composite material shell structures using the finite element method. Volume integral expressions to evaluate the dynamic energy release rate in a through‐thickness three‐dimensional crack are derived. Using the domain integral, the energy release rate computation is implemented in the DYNA3D explicit non‐linear dynamic finite element analysis program wherein crack propagation is modelled by releasing the constraints between initially constrained node pairs. The implementation enables the program to either determine the energy resistance response for the material (provided experimental data is available) or predict the rate of crack propagation in shell structures. The numerical implementation was verified by simulating mode I and mode III slow crack growth problems in semi‐infinite transversely isotropic media, for which analytic solutions are available. Oscillations of energy following the release of nodal constraints as the crack propagates in discrete increments were suppressed using light mass proportional damping and a moving averaging scheme. Copyright © 2004 John Wiley & Sons, Ltd.