Simulation of crack induced nonlinear elasticity using the combined finite-discrete element method
Introduction
Non-Destructive Testing (NDT) for detection and quantification of defects in solids (e.g. crack, delamination, debonding, pore and inter-granular contact) is of significant industrial and academic importance in many areas [1], [2], [3], [4], [5]. Probing and imaging using ultrasonic waves is a leading tool for such applications and there is a need for continued development of robust ultrasonic techniques in NDT [2], [3]. Among the many ultrasonic wave applications in NDT, a major class of them employ the principles of wave reflection, transmission or scattering. While extremely useful, these linear ultrasonic techniques are less capable of accurately detecting contact-type defects, and are also less sensitive to micro or closed cracks [2], [3], [4]. Since defects can behave nonlinearly under sufficient excitation, the Nonlinear Elastic Wave Spectroscopy (NEWS) methods have shown remarkable potential for defect detection and characterization [2], [3], [4], [6], [7], [8], [9], [10]. The nonlinear methods often involve exciting a solid with an ultrasonic signal of certain frequency and generating an output frequency spectrum consisting of harmonics and subharmonics of the exciting frequency. These effects, often referred to as Contact Acoustic Nonlinearity (CAN) [10], are mainly induced by clapping and frictional contacts of the defects [2], [3], [11]. Because of their high sensitivity, nonlinear ultrasonic approaches have seen an increasing interest in the NDT community over the past decades [1], [2], [4], [9], [12], [13].
A large number of experimental studies have been conducted using NEWS techniques for various types of defect detection in solid materials such as composite plates, metals, concretes and rocks, and have demonstrated successes in the field of NDT [6], [7], [8], [14], [15], [16], [17], [18], [19], [20], [21]. However, to date, the underlying microscopic mechanism of nonlinear techniques for defect detection is still poorly understood [4]. Numerical simulations, which are capable of providing detailed analyses of the nonlinear behavior at a level of spatial and temporal resolution not accessible experimentally, are therefore necessary for decoding the complicated mechanism associated with the nonlinear ultrasonic techniques. In addition to the ease of implementation, the numerical approaches are indispensable also because they can provide dedicated comparisons of nonlinear indicators with experimental results, and thus link measured macroscopic events to defect internal parameters (both physical and geometric) and in this way, a complete characterization of the defects can be achieved [2].
Simulation of ultrasonic wave propagation in solids with defects by considering the nonlinearity introduced by them is challenging, and has been the object of study for two decades [2], [3]. The numerical approaches used in this regard mainly include Finite-Difference Time-Domain (FDTD) and Finite Element Method (FEM). For example, Sarens et al. [22] implemented a three-dimensional finite difference, staggered grid simulation to model the contact nonlinear acoustic generation in a composite plate containing an artificial defect; Marhenke et al. [23] used FDTD simulations to cross-validate the simulated interference effects resulting from multiple ultrasonic reflections within the delamination layers with the laboratory experiments. In general, the FDTD is easy to implement; however, it has many restrictions, e.g. the defects in the FDTD simulations are usually restricted to rectangular shapes [2]. As a more flexible alternative, FEM is widely used in crack-wave interaction simulations. In particular, Kawashima et al. [24] used a FEM model to study CAN in which Rayleigh waves were employed to detect surface cracks; Blanloeuil et al. [12], [25] studied the nonlinear scattering of ultrasonic waves by closed cracks subjected to CAN; to investigate the clapping and friction induced nonlinearity in solids containing cracks, Van Den Abeele and colleagues [2], [3], [11], [26] implemented a series of comprehensive normal and tangential constitutive models into FEM to control the nonlinear behavior of crack surfaces. Among the many numerical simulations, some of them use hypothetical defects in which artificial nonlinear stress-strain relations are introduced into special elements to represent defects [24]; others employ physical defects by splitting the computational nodes along defects, and then the normal and tangential contact stresses, which are calculated based on the relative distances between the corresponding Gauss points located on the defect surfaces, are applied to the same Gauss point pairs as boundary conditions for the bulk material simulation [2], [3], [4], [5], [11], [12], [22], [23], [25], [26], [27], [28], [29]. The former only captures the defect behavior in an approximate manner; the latter may be difficult to explicitly realize complicated scenarios such as defects with irregular shapes and especially, the interactions between many defects of different types.
From a computational mechanics point of view, a solid with defects is essentially a combination of continua (bulk material) and discontinua (interaction between defect surfaces). Considering this, a numerical tool that has the capability of handling continua and discontinua simultaneously would be helpful. Fortunately, a recently developed numerical method – the combined finite–discrete element method (FDEM) [30], [31], [32], [33], which merges finite element-based analysis of continua with discrete element-based transient dynamics, contact detection and contact interaction solutions of discontinua, provides a natural solution for such simulation. To date, a systematic application of FDEM in NDT is not available in the literature. The goal of this paper is to introduce FDEM to the NDT community, and to demonstrate its power on simple problems.
A simple FDEM realization of a solid plate with a crack is presented in Fig. 1 where the solid (excluding the crack since it is considered as void) is discretized into finite elements to capture the motion and deformation of the bulk material, and the finite elements along the two sides of the crack also behave as discrete elements to track the normal and tangential interactions between crack surfaces. By employing FDEM, the system can be explicitly described and particularly, the contacts along the sides of the defects can be uniformly processed using well-developed discrete element method (DEM)-based algorithms.
The focus of the current work is to demonstrate the applicability of FDEM to the simulation of crack induced nonlinear elasticity in solids, and to introduce the approach as another alternative for NDT numerical based analysis. A comparison of the simulated results with laboratory experiments is beyond the scope of the present paper and thus it is left for future work. In the following sections, we first provide a brief introduction to the theories of FDEM. Then we illustrate the numerical model setup and present how the normal (clapping) and tangential (friction) contact may influence the nonlinear behavior of a cracked solid. The applicability of FDEM for NDT simulation is demonstrated, and the corresponding conclusions are drawn.
Section snippets
The combined finite-discrete element method (FDEM)
FDEM was originally developed by Munjiza in the early 1990s to simulate the material transition from continuum to discontinuum [33]. The essence of this method is to merge the algorithmic advantages of DEM with those of the FEM. The theory of the FDEM can be broken down into the following parts: governing equations, finite strain-based formulation for deformation description, contact detection, and contact interaction algorithms [34], [35], [36].
FDEM simulation examples
As an introductory illustrative application of the use of FDEM for crack induced nonlinearity simulations, two-dimensional rectangle models with a single crack are excited using compressive sinusoidal waves. We first demonstrate the model setup, then a set of simulations using different combinations of crack aperture and excitation amplitudes are presented for the model with a single horizontal crack to examine the applicability of FDEM for this type of applications. Finally, a more detailed
Summary
In this paper we have introduced FDEM by describing the application to an intact solid and two simple crack models – one normal to the forcing wave and one with crack at a 30° angle to the wave. In the FDEM model, the solid is discretized into finite elements to capture the wave propagation in the bulk material, and the finite elements along the two sides of the crack also behave as discrete elements to track the normal and tangential interactions between crack surfaces. By employing FDEM, the
Acknowledgments
The Department of Energy Office of Basic Energy Research, Geoscience, supported this work. Technical support and computational resources from the Los Alamos National Laboratory Institutional Computing Program are highly appreciated.
References (59)
- et al.
Modelling of nonlinear crack–wave interactions for damage detection based on ultrasound—A review
J. Sound Vib.
(2014) - et al.
Two dimensional modeling of elastic wave propagation in solids containing cracks with rough surfaces and friction – Part II: Numerical implementation
Ultrasonics
(2018) - et al.
Two dimensional modeling of elastic wave propagation in solids containing cracks with rough surfaces and friction – Part I: Theoretical background
Ultrasonics
(2018) - et al.
Finite element modeling of nonlinear acoustics/ultrasonics for the detection of closed delaminations in composites
Ultrasonics
(2017) - et al.
Numerical study of nonlinear interaction between a crack and elastic waves under an oblique incidence
Wave Motion
(2014) - et al.
Micro-damage diagnostics using nonlinear elastic wave spectroscopy (NEWS)
NDT&E Int.
(2001) - et al.
CAN: an example of nonclassical acoustic nonlinearity in solids
Ultrasonics
(2002) - et al.
Three-dimensional finite element simulation of closed delaminations in composite materials
Ultrasonics
(2012) - et al.
Analytical extension of Finite Element solution for computing the nonlinear far field of ultrasonic waves scattered by a closed crack
Wave Motion
(2016) - et al.
Impact-based nonlinear acoustic testing for characterizing distributed damage in concrete
J. Nondestruct. Eval.
(2017)
Modeling of delamination detection utilizing air-coupled ultrasound in wood-based composites
NDT&E Int.
Nonlinear acoustic response through minute surface cracks: FEM simulation and experimentation
Ultrasonics
2D finite element modeling of the non-collinear mixing method for detection and characterization of closed cracks
NDT&E Int.
Method of memory diagrams for mechanical frictional contacts subject to arbitrary 2D loading
Int. J. Solids Struct.
Numerical simulation of Rayleigh wave interaction with surface closed cracks under external pressure
Wave Motion
Influence of boundary constraints on stress heterogeneity modelling
Comput. Geotech.
Modelling stress-dependent permeability in fractured rock including effects of propagating and bending fractures
Int. J. Rock Mech. Min. Sci.
Role of natural fractures in damage evolution around tunnel excavation in fractured rocks
Eng. Geol.
Continuum–discontinuum analysis of failure mechanisms around unsupported circular excavations in anisotropic clay shales
Int. J. Rock Mech. Min. Sci.
Assessment of crack initiation and propagation in rock from explosion-induced stress waves and gas expansion by cross-hole seismometry and FEM–DEM method
Int. J. Rock Mech. Min. Sci.
Large scale simulation of red blood cell aggregation in shear flows
J. Biomech.
A calibration procedure for two-dimensional laboratory-scale hybrid finite–discrete element simulations
Int. J. Rock Mech. Min. Sci.
Nonlinear elastic wave spectroscopy (NEWS) techniques to discern material damage, Part I: nonlinear wave modulation spectroscopy (NWMS)
Res. Nondestruct. Eval.
Nonlinear elastic wave spectroscopy (NEWS) techniques to discern material damage, Part II: single-mode nonlinear resonance acoustic spectroscopy
Res. Nondestruct. Eval.
Nonlinear Mesoscopic Elasticity: the Complex Behaviour of Rocks, Soil, Concrete
Nonlinear ultrasonic techniques for nondestructive assessment of micro damage in material: A review
Int. J. Precis. Eng. Man.
Resonance and elastic nonlinear phenomena in rock
J. Geophys. Res. Solid Earth
A local defect resonance to enhance acoustic wave-defect interaction in ultrasonic nondestructive evaluation
Appl. Phys. Lett.
Crack detection using nonlinear acoustics and piezoceramic transducers—instantaneous amplitude and frequency analysis
Smart Mater Struct.
Cited by (20)
Effect of intermittent joint distribution on the mechanical and acoustic behavior of rock masses
2023, Journal of Rock Mechanics and Geotechnical EngineeringImpact of fracture normal and shear stiffnesses on the scattering attenuation of P and S waves in a naturally fractured rock
2022, Journal of Applied GeophysicsFourier spectral element for simulation of vibro-acoustic modulation caused by contact nonlinearity in the beam
2022, Thin-Walled StructuresCitation Excerpt :Modeling the CAN is the essential step in VAM simulation. The CAN was modeled using a variety of numerical approaches, including the time-domain finite difference method [34,35], finite discrete element method [36], local interaction simulation approach [37], time-domain spectral finite method [38], and FEM [39–41]. Zhang et al. investigated the effect of microscale surface contact on the dynamic responses of a composite joint using fractal contact theory and the elastic–viscoelastic correspondence principle in a multi-scale numerical model [42].
Characterization of delamination features of orthotropic CFRP laminates using higher harmonic generation technique: Experimental and numerical studies
2022, Composite StructuresCitation Excerpt :Azusa Sugawara et al. [24] further extended the FDTD model with damped double nodes to actual testing conditions and compared the simulation results with experimental data. Although the FDTD method is widely used, it has some limitations, such as the inflexibility in the model geometry [21,23]. As a more flexible alternative, FEM is widely applied in CAN interaction simulations.
A simplified analytical model for the investigation of contact acoustic nonlinearity in pipe structures
2021, International Journal of Mechanical SciencesCitation Excerpt :Therefore, analysis of the generation of second harmonic waves in a pipe structure and a proper nonlinear index are both in high demand for the measurement of CAN. Apart from the analytical model analysis of CAN, different numerical methods were also developed to study the interaction between guided waves and close crack in a structure, including finite discrete element method (FDEM) [32], local interaction simulation approach (LISA) [33–35], finite difference time domain method (FDTD) [36,37], time-domain spectral finite method [38] and finite element method (FEM) [39,40]. The FEM model was adopted widely for the development of CAN model in both steel and composite components [41,42].