Elastic Half-space Research Articles

Partial Slipping of Flat Punch in Thermomechanical Contact with Elastic Half-Space

Partial Slipping of Flat Punch in Thermomechanical Contact with Elastic Half-Space

Scattering of SH wave in an elastic half-space by a semi-elliptical crater with surface elasticity

Cratering occurs inevitably on the surface of materials at all scales due to comprehensive factors. The presence of craters may have a significant influence on the mechanical performance of the surface of an elastic body. In this paper, we study the SH wave propagation in an elastic semi-infinite medium whose surface has a semi-elliptical crater. We explore the scattering of the SH wave by the crater when surface elasticity is incorporated in the model of deformation to account for possible small-scale effects of the crater. We obtain a series solution for the scattered wave function by employing the Mathieu functions in the corresponding elliptical coordinate system. We discuss in some numerical examples the dynamic stress concentration around the crater relative to the wave frequency and few geometric parameters.

Computationally Efficient Emulation of Spheroidal Elastic Deformation Sources Using Machine Learning Models: A Gaussian‐Process‐Based Approach

AbstractElastic continuum mechanical models are widely used to compute deformations due to pressure changes in buried cavities, such as magma reservoirs. In general, analytical models are fast but can be inaccurate as they do not correctly satisfy boundary conditions for many geometries, while numerical models are slow and may require specialized expertise and software. To overcome these limitations, we trained supervised machine learning emulators (model surrogates) based on parallel partial Gaussian processes which predict the output of a finite element numerical model with high fidelity but >1,000× greater computational efficiency. The emulators are based on generalized nondimensional forms of governing equations for finite non‐dipping spheroidal cavities in elastic halfspaces. Either cavity volume change or uniform pressure change boundary conditions can be specified, and the models predict both surface displacements and cavity (pore) compressibility. Because of their computational efficiency, using the emulators as numerical model surrogates can greatly accelerate data inversion algorithms such as those employing Bayesian Markov chain Monte Carlo sampling. The emulators also permit a comprehensive evaluation of how displacements and cavity compressibility vary with geometry and material properties, revealing the limitations of analytical models. Our open‐source emulator code can be utilized without finite element software, is suitable for a wide range of cavity geometries and depths, includes an estimate of uncertainties associated with emulation, and can be used to train new emulators for different source geometries.

Effect of variable thermal conductivity in semiconducting medium underlying an elastic half-space

Effect of variable thermal conductivity in semiconducting medium underlying an elastic half-space

Semi-analytical method for calculating ground vibrations from a tunnel in a homogeneous half-space with an irregular surface

This paper proposes a novel two-dimensional (2D) semi-analytical method for calculating ground vibrations from a tunnel situated in a homogeneous half-space with an irregular surface. The circular tunnel is conceptualized as an elastic solid, while the soil is modelled as an elastic, isotropic, and homogeneous half-space with an irregular surface. A virtual horizontal interface is introduced to divide the soil domain into an irregular region with an arbitrary-shaped surface and a half-space with a circular hollow. The wavefield scattered by the irregular surface is simulated by the boundary integral equation. Through the application of the discrete wavenumber approach and expansion of wavefield on the irregular surface, the boundary integral equation is transformed into a set of simultaneous matrix equations containing unknown expansion coefficients. By utilizing the transformation between cylindrical and plane waves, the boundary conditions on soil-tunnel interface are satisfied, yielding a solution for a harmonic point load acting on the tunnel in a half-space with an irregular surface. The proposed method can simulate topographies with arbitrary shapes and only requires the discretization of the irregular part of the ground surface, which provides a highly efficient tool to investigate the propagation characteristics of train-induced vibrations in complex topographies. The proposed method is verified through comparisons with the existing analytical methods and the finite element method. The dynamic responses of a tunnel embedded in a half-space with a Gaussian-shaped surface are meticulously investigated. A case study involving a site for the Shenzhen metro line 5 is performed. The numerical results demonstrate that the presence of irregular topography alters the distribution of ground vibrations, leading to significant differences in responses within the irregular region and certain disparities in the regular region. The effect of the location, size, and undulating direction of the irregular topography on ground vibrations is frequency-dependent.

Coseismic deformation and seismogenic structure of the 2024 Hualien Earthquake measured by InSAR and GNSS

On April 3, 2024, an M7.3 earthquake occurred in the offshore area of Hualien County, Taiwan, China. The seismogenic structure at the epicentral location was highly complex, and studying this earthquake is paramount for understanding regional fault activity. In this study, we employed ascending and descending orbit Sentinel-1 Synthetic Aperture Radar (SAR) data and utilized differential interferometry (InSAR) technique to obtain the co-seismic deformation field of this event. The line-of-sight deformation field revealed that the main deformation caused by this earthquake was predominantly uplift, with maximum uplift values of approximately 38.8 ​cm and 46.1 ​cm for the ascending and descending orbits, respectively. By integrating the three-dimensional GNSS co-seismic deformation field, we identified the seismogenic fault located in the offshore thrust zone east of Hualien, trending towards the northwest. The fault geometry parameters, obtained through the inversion of an elastic half-space homogeneous model, indicated an optimal fault strike of 196°, a dip angle of 30.9°, and an average strike-slip of 0.4 ​m and dip-slip of −2.6 ​m. This suggests that the predominant motion along the seismogenic fault is thrusting. The distribution of post-seismic Coulomb stress changes revealed that aftershocks mainly occurred in stress-loaded regions. However, stress loading was observed along the northern segment of the Longitudinal Valley Fault, with fewer aftershocks. This highlights the importance of closely monitoring the seismic hazard associated with this fault segment.

SH waves in orthotropic piezomaterials considered surface effects

Nowadays, orthotropic piezomaterials have enormous potential in producing electronic devices owing to their inherent merits, e.g., intrinsically robust piezoelectricity. Nevertheless, the previous theoretical studies merely report surface effects on the physical properties of SH waves in a layered transversely isotropic piezoelectric structure. Filling such a blank is a special necessity for optimizing Surface Acoustic Wave (SAW) sensors. Thus, the present paper aims to determine how surface effects and interfacial imperfection affect the attributes of SH waves in a layered orthotropic piezoelectric semi-space by means of surface piezoelectricity theory and a linear spring model. The structure mainly consists of the piezoelectric film's bottom and top surface layers and the imperfect bonding between the piezoelectric film and elastic half-space. To explicitly figure out surface and imperfect interface effects on the mechanical behavior of SH waves in a layered piezoelectric half-space, we extensively derive the phase velocity equations of SH waves in an orthotropic piezoelectric plate and an orthotropic piezoelectric half-space. Then, a comparison study among SH waves in a layered piezoelectric semi-space, a piezoelectric plate, and a piezoelectric half-space is conducted. As a result, it is revealed that the attributes of SH waves dramatically vary as the piezoelectric films’ thickness changes from microscale to nanoscale. In addition, phase velocity increases or decreases with an increasing surface elastic constant or surface density due to the enhancement of stiffness or density of the whole structure. When the ratio of piezoelectric film thickness to wavelength is small, imperfect parameters substantially reduce the phase velocity owing to a decreasing interface elastic stiffness. The comparison study shows that the imperfect parameter does not influence the waves’ velocity, while the ratio of film thickness to wavelength is relatively considerable. Anisotropic piezoelectric constant essentially weakens surface effects on the waves’ velocity. Intriguingly, a linear function is acquired and is used to measure surface density via the slope of the linear function. The summaries supply unique guidelines and instructions to open an avenue for designing and manipulating surface acoustic wave nanosensors in the future.

Conical indentation over a transversely isotropic and layered elastic half-space

We propose a novel method for solving the static response of a conical indenter on a transversely isotropic and layered elastic half-space. The newly developed Fourier-Bessel series (FBS) system of vector functions, along with the unconditionally stable dual-variable and position method, is employed to derive the Green's function in the transversely isotropic and layered elastic half-space under a vertical ring load on the surface. To calculate the response at different field points on the surface, we apply discrete love numbers within the FBS vector system. The load densities in the discretized rings within the contact radius of the conical indenter are determined using the integral least-square method, along with a self-adaptive algorithm developed in this study. Finally, the relationship between the indentation depth (vertical displacement) and the applied load is obtained through force balance between the external load and the summed contact traction. The developed scheme is validated using existing exact solutions for the reduced homogeneous half-space case. Selected numerical results clearly demonstrate the effect of anisotropic material and layering on the indentation response. It is observed that, regardless of whether the structure is a stratified half-space or a layered structure with a rigid substrate, the material properties in the top layer have the most significant influence on the indentation behavior. In the case of a layered structure with an underlying elastic half-space, the material properties in the interlayer and bottom layer could also affect the indentation behaviors.

Active shortening and aseismic slip along the Cephalonia Plate Boundary (Paliki Peninsula, Greece): Evidence from InSAR and GNSS data

We present a comprehensive analysis of geodetic data, including InSAR and GNSS, to assess the interseismic deformation of the Paliki peninsula in western Cephalonia, Greece. The region is prone to frequent earthquakes, due to its proximity to the Cephalonia Transform Fault (CTF), a 140 km long dextral strike-slip fault (striking NNE-SSW) that accommodates the relative motion between the Apulian and Aegean lithospheric plates. Our analysis covers the period from 2016 to 2022 and leverages LiCSBAS, an open-source package, for InSAR time series analysis with the N-SBAS method. The results of the InSAR analysis demonstrate deformation rates between 2 and 5 mm/yr in the line-of-sight (LOS) direction of the satellite. We constructed E-W velocity profiles that provided a velocity gradient terminating against the outcropping trace of the seismic fault of the 3 February 2014 earthquake (M6.1) near Atheras. This geodetic evidence indicates a ∼ 1 mm/yr minimum aseismic slip (fault creep) motion along the February 2014 seismic fault, during 2016–2022. Positive LOS values in both satellite imaging geometries show that the coastal town of Lixouri experiences uplift of 1 mm/yr. East-West velocity cross-sections across western Cephalonia including Gulf of Argostoli reveal several velocity discontinuities, possibly bounded by active faults and/or landslides. The E-W shortening rate between Lixouri and Argostoli areas amounts to 1.5 mm/yr corresponding to −187 ns/yr (nanostrain/yr) of tectonic strain. Our results suggest a complex deformation pattern on the Paliki peninsula with strain accumulation along strike-slip and reverse-slip faults. We also inverted GNSS velocities from Italy and Greece, across the Cephalonia segment of CTF (assuming elastic half-space) and obtained a locking depth of 13 km, and a slip rate of 17.3 ± 0.8 mm/yr.

Static deformation of a two-phase medium consisting of a rigid boundary elastic layer and an isotropic elastic half-space induced by a very long tensile fault

An analytical model is proposed to determine the static response of a multi-layered medium comprising of an underlying isotropic (ISO) elastic half-space that is in welded contact with an isotropic (ISO) elastic layer of uniform thickness H with a rigid boundary, induced by a tensile dislocation embedded in the layer. The integral expressions for stress fields are obtained by using the airy stress function. These integrals are calculated approximately, using Sneddon’s approach, by substituting the terms under the integral sign with a finite sum of exponential expressions. It is studied numerically and graphically how stresses vary with the change in the epicentral distance for different source locations embedded in the layer. The effect of the rigidity ratio is also analysed on the stress field. To examine the impact of the internal boundary, the stress fields in the layered and uniform half-spaces are compared graphically.

Non-adhesive and adhesive contacts of an elastic quarter-or eighth-space with freely sliding sides

The contact of an elastic quarter- or eighth-space is studied under the condition that the movement of the side surface of the quarter-space is constrained: It can slide freely along the plane of the side surface but its normal movement is blocked (for example, by a rigid wall). The solution of this contact problem can be easily achieved by additionally applying a mirrored load to an elastic half-space. Non-adhesive contact and the Johnson-Kendall-Roberts (JKR)-type adhesive contact between a rigid sphere and an elastic quarter-space under such a boundary condition is numerically simulated using the fast Fourier transform (FFT)-assisted boundary element method (BEM). Contacts of an elastic eighth-space are investigated using the same idea. Depending on the position of the sphere relative to the side edge, different contact behavior is observed. In the case of adhesive contact, the force of adhesion first increases with increasing the distance from the edge of the quarter-space, achieves a maximum, and decreases further to the JKR-value in large distance from the edge. The enhancement of the force of adhesion compared to the half-space-contact is associated with the pinning of the contact area at the edge. We provide the maps of the force of adhesion and their analytical approximations, as well as pressure distributions in the contact plane and inside the quarter-/eighth-space.

Development of governing partial differential equations of reinforcing thin films

In this paper, a novel interface model for elastic isotropic thin films perfectly bonded to the surrounding elastic media is introduced. Considering an average stress and displacement through the thin film thickness, the shear stress discontinuity across the interface is expressed in terms of spatial and time derivatives of displacement components. The derived interfacial equations are presented in a general form. As such, they are applicable to asymmetric elastodynamic problems of anisotropic solids with thin film interfaces. Employing the introduced interface model and Hankel transform technique, static Green's functions of an elastic isotropic half-space reinforced by a buried thin film are obtained. Results of some special cases, including axisymmetric loading, inextensible membrane, unreinforced half-space, and surface-coated half-space, are also presented. The robustness and effectiveness of the proposed interface model are shown by two solved numerical examples, and some elastic responses are plotted to show the impact of thin film stiffness. According to the numerical results, modeling elastic thin films as inextensible membranes can significantly alter the elastic responses.

Deformation of solid earth by surface pressure: equivalence between Ben-Menahem and Singh’s formula and Sorrells’ formula

SUMMARY Atmospheric pressure changes on Earth’s surface can deform the solid Earth. Sorrells derived analytical formulae for displacement in a homogeneous, elastic half-space, generated by a moving surface pressure source with speed $c$. Ben-Menahem and Singh derived formulae when an atmospheric P wave impinges on Earth’s surface. For a P wave with an incident angle close to the grazing angle, which essentially meant a slow apparent velocity $c_a$ in comparison to P- ($\alpha ^{\prime }$) and S-wave velocities ($\beta ^{\prime }$) in the Earth ($c_a \ll \beta ^{\prime } \lt \alpha ^{\prime }$), they showed that their formulae for solid-Earth deformations become identical with Sorrells’ formulae if $c_a$ is replaced by $c$. But this agreement was only for the asymptotic cases ($c_a \ll \beta ^{\prime }$). The first point of this paper is that the agreement of the two solutions extends to non-asymptotic cases, or when $c_a /\beta ^{\prime }$ is not small. The second point is that the angle of incidence in Ben-Menahem and Singh’s problem does not have to be the grazing angle. As long as the incident angle exceeds the critical angle of refraction from the P wave in the atmosphere to the S wave in the solid Earth, the formulae for Ben-Menahem and Singh’s solution become identical to Sorrell’s formulae. The third point is that this solution has two different domains depending on the speed $c$ (or $c_a$) on the surface. When $c/\beta ^{\prime }$ is small, deformations consist of the evanescent waves. When $c$ approaches Rayleigh-wave phase velocity, the driven oscillation in the solid Earth turns into a free oscillation due to resonance and dominates the wavefield. The non-asymptotic analytical solutions may be useful for the initial modelling of seismic deformations by fast-moving sources, such as those generated by shock waves from meteoroids and volcanic eruptions because the condition $c / \beta ^{\prime } \ll 1$ may be violated for such fast-moving sources.

A Note on Exact Closed-Form Solutions for Lamb’s Problem with an Arbitrary Poisson’s Ratio

Abstract Lamb’s Problem, a cornerstone in theoretical seismology, investigates the displacement of an elastic half-space due to surficial or buried point sources. Although the exact closed-form solution for surface displacements from surficial sources is well established, the solutions for buried sources in the 3D case have remained limited, either confined to integral forms or applicable only for certain Poisson’s ratios. This study addresses this gap by generalizing the existing closed-form solutions to accommodate Poisson’s ratio for buried sources. Our solutions are expressed in terms of elementary algebraic expressions, elliptic integrals, and Gauss hypergeometric functions. The derived solution is validated by numerical tests. With the newly derived solution, we also explore the properties of the leaking mode, which stems from the complex conjugate roots of the Rayleigh equation.

Surface effect on the partial-slip contact of a nano-sized flat indenter

The partial-slip contact at nano-scale is always influenced by a surface effect. A new model based on the Gurtin-Murdoch (G-M) surface elasticity theory is proposed in this paper to study the partial-slip contact behavior of a rigid flat nano-indenter on an elastic half-space. The surface effect in the partial-slip contact is characterized via a non-classical boundary condition involving a surface-induced normal traction related to the residual surface stress and a surface-induced tangential traction related to the surface elasticity, the influences of which on stress and displacement fields and the stick-slip state at the contact surface are investigated. It is found that, with the increase of residual surface stress, both the normal pressure and vertical and horizontal displacements in the contact zone decrease, the relative slip between indenter and substrate reduces and the stick region is enlarged. Such effects are basically attributed to the action of the surface-induced normal traction opposite to the externally applied compression. However, the increase of surface elastic constants is conductive to the relative slip and results in a shrinkage of the stick region, which is attributed to the action of the surface-induced tangential traction opposite to the frictional stress. An interesting phenomenon is further unveiled that, when the frictional coefficient increases, the dominant role in affecting the stick-slip state changes from the surface-induced tangential traction to the normal one, thus inspiring a feasible route to manipulate the surface effect by tuning the frictional coefficient of substrate. The present research enables one to better understand the partial-slip contact behavior of nano-indenters, which is of guiding value for anti-wear designs in nano-mechanical devices.

On the depth of cylindrical indentation of an elastic half-space for different types of displacement boundary conditions

Nonlinear relationships between the penetration depth δ and applied force F for cylindrical indentation of a half-space are derived corresponding to two types of displacement boundary conditions: (i) zero vertical displacement imposed at points x=±b of the free-surface of a half-space, and (ii) zero vertical displacement imposed at a point z=h below the indenter. The requirement that the work done by the indentation force must be equal to the work done by the contact pressure defines the minimum values of b and h for which the indentation is geometrically and physically possible. The minimum value of b is found to be about 10 times greater than the semi-width of the contact zone a=(4FR/πE∗)1/2, with the corresponding minimum value of the indentation depth δmin≈3.5δ0, where δ0=2F/πE∗ is the height of the contact zone under applied force F, and E∗=E/(1−ν2) is the plane strain elastic modulus. The minimum value of h is equal to about 21a in the case ν=1/3, and about 27a in the case ν=1/2. The relationship between h and b, corresponding to the same indentation depth δ≥δmin, is derived and shown to be nearly linear for large values of h/a and b/a. The obtained results are used to estimate the indentation depth of finite-size elastic bodies under different boundary constraints. Comparison with FEM prediction is made.

Approximate Near-Front Ray Solutions in the Axisymmetric Strain Dynamics of a Linear Elastic Half-Space

Approximate Near-Front Ray Solutions in the Axisymmetric Strain Dynamics of a Linear Elastic Half-Space

Minimal Conditioned Stiffness Matrices with Frequency-Dependent Path Following for Arbitrary Elastic Layers over Half-Spaces

This paper introduces an efficient computational procedure for analyzing the propagation of harmonic waves in layered elastic media. This offers several advantages, including the ability to handle arbitrary frequencies, depths, and the number of layers above an elastic half-space, and efforts to follow dispersion curves and flag up possible singularities are investigated. While there are inherent limitations in terms of computational accuracy and capacity, this methodology is straightforward to implement for studying free or forced vibrations and obtaining relevant response data. We present computations of wavenumber dispersion diagrams, phase velocity plots, and response data in both the frequency and time domains. These computational results are provided for two example cases: plane strain and axisymmetry. Our methodology is grounded in a well-conditioned dynamic stiffness approach specifically tailored for deep-layered strata analysis. We introduce an innovative method for efficiently computing wavenumber dispersion curves. By tracking the slope of these curves, users can effectively manage continuation parameters. We illustrate this technique through numerical evidence of a layer resonance in a real-life case study characterized by a fold in the dispersion curves. Furthermore, this framework is particularly advantageous for engineers addressing problems related to ground-borne vibrations. It enables the analysis of phenomena such as zero group velocity (ZGV), where a singularity occurs, both in the frequency and time domains, shedding light on the unique characteristics of such cases. Given the reduced dimension of the problem, this formulation can considerably aid geophysicists and engineers in areas such as MASW or SASW techniques.

Analytical prediction and numerical verification of stress concentration profile around an in-situ tow break in resin-impregnated filament-wound composites

A new empirical analytical approach is developed for predicting the stress concentration profile around an in-situ tow break in filament-wound composites. A shear-lag analysis is firstly performed to solve for the perturbational axial displacement of the broken tow and resultant debonding lengths. Solution of stress field caused by tangential load on the surface of an elastic half space is then utilized in combination with superposition principle to obtain the overload magnitudes in the neighboring tows. Subsequently, high-fidelity finite element analysis on a representative uni-directional laminate model under different stress states is performed, and excellent overall agreement is observed between analytical and numerical predictions. The proposed method takes into account essential aspects such as transversely isotropic material properties, in-situ stress states and their effect on the interfacial frictional forces in the debonded interfaces, and thus provides a convenient way to evaluate the stress concentration factors in damaged filament-wound composites. In addition, this approach can be applied to yield auxiliary failure evaluation criteria for statistical strength prediction or finite element modeling of filament-wound composites or similar structures.

Elastic and Elastoplastic Contact Mechanics of Concentrated Coated Contacts

Machines operate under increasingly harsher contact conditions, causing significant wear and contact fatigue. Sub-surface stresses are responsible for the premature contact fatigue of rolling element bearings, meshing gears, and cam–follower pairs. Surface protection measures include hard, wear-resistant coatings. Traditionally, contact integrity has been predicted using classical Hertzian contact mechanics. However, the theory is only applicable when the contact between a pair of ellipsoidal solids of revolution may be considered as a rigid indenter penetrating a semi-infinite elastic half-space. Many coatings act as thin bonded elastic layers that undergo considerably higher pressures than those predicted by the classical theory. Furthermore, inelastic deformation of bonded solids can cause plastic flow, work-hardening, and elastoplastic behaviour. This paper presents a comprehensive, integrated contact mechanics analysis that includes induced sub-surface stresses in concentrated counterformal finite line contacts for all the aforementioned cases. Generated pressures and deformation are predicted for hard coated surfaces, for which there is a dearth of relevant analysis. The contact characteristics, which are of particular practical significance, of many hard, wear-resistant advanced coatings are also studied. The paper clearly demonstrates the importance of using efficient semi-analytical, detailed holistic contact mechanics rather than the classical idealised methods or empirical numerical ones such as FEA. The novel approach presented for the finite line contact of thin-layered bonded solids has not hitherto been reported in the open literature.