A spectral method for numerical elastodynamic fracture analysis without spatial replication of the rupture event

Abstract

Perrin et al. (1995) and Geubelle and Rice (1995) have introduced a spectral method for numerical solution of two-and three-dimensional elastodynamic fracture problems. The method applies for ruptures confined to a plane separating homogeneous elastic half spaces. In this method, the physical variables, such as the traction components of stress and displacement discontinuity on the rupture plane, are represented as Fourier series in space with time-dependent coefficients. An analytical solution is found for each Fourier mode, in that each Fourier coefficient for stress is expressed by the time convolution of the corresponding coefficient for displacement with a convolution kernel specific to the rupture mode. Once the 2D formulation of the method is known, the method is readily generalizable to 3D problems in that it involves only linear combinations of the convolution kernels found for each rupture mode in 2D. This conceptual simplicity has, however, a major drawback : due to the Fourier series representations of the physical variables, the problem solved is in fact an infinite and periodic replication of rupture events on the fracture plane. So, in order to study the evolution of a single rupture, one has to use a spatial period large enough in order that the waves coming from the replication cracks do not enter the zone of interest during the time duration studied, or provide negligible stress alteration when they do arrive. We show here how to rigorously offset this defect while retaining the modal independence.Once expressed in the spatial domain, the method amounts to truncating in space the space-time convolution kernels, in a manner that provides an exact evaluation for all positions within the rupture domain (where the constitutive law between stress and displacement discontinuity is to be imposed), but not outside. In order for the method to be identical in structure to the method of Perrin et al. (1995) and Geubelle and Rice (1995). the oeriod of the Fourier series is requested to be only twice as large as the rupture domain of interest. The only difference, then, to the original spectral method is that the convolution kernels in the Fourier domain require more elaborate calculations to be established, but this has to be done only once to allow simulations on a given domain. '$3 1997

Dynamic modelling of the flat 2-D crack by a semi-analytic BIEM scheme

“…More recent modelling works, based on the displacement discontinuity method, the regularization technique and the time-domain representation, dealt with the at two-dimensional (2-D) anti-plane crack [3; 4], at 3-D crack [5][6][7], non-planar 2-D anti-plane and=or in-plane cracks [8][9][10][11][12][13], and the non-planar 3-D crack [14][15][16]. Meanwhile, an alternative formulation, based on a Fourier transform in space, was also developed [17][18][19].…”

A coupled FEM–Bernstein approach for computing the J k integrals

“…Geubelle and Rice [30] and Cochard and Rice [31] developed a formulation based on a Fourier representation in spatial coordinates along a fracture plane of the tractions and relative displacements. Legendre polynomials were introduced by Woo and Lee [32] yielding a p-FEM implementation, while Rahulkumar et al [33] developed a p-version of the singular quarter-point element.…”

The inverse problem of seismic fault determination using part time measurements

“…In simple cases listed above, some forward models in infinite elastic medium can give good prediction of the acceleration signals to be compared to measurements [Campillo 1983;Cochard and Rice 1997;Vallée 2003;Peyrat et al 2004;Lapusta et al 2001]. True inverse problems are divided into two kinds.…”