Dimensionless Stress Research Articles

Viscoplastic corner eddies

When viscous fluid in a corner is disturbed, eddies can form in the absence of inertia. Examples of flow configurations in which this motion occurs include flow through an abrupt contraction and over a cavity. Six decades ago, Moffatt (J. Fluid Mech., vol. 18, 1964, pp. 1–18) calculated the slow viscous flow of Newtonian fluids in sharp corners, detailing his eponymous ‘Moffatt eddies’. In this study we examine corner flows of viscoplastic materials, a class of non-Newtonian fluids which exhibit solid-like behaviour for stresses below a yield stress. Specifically, we consider a Bingham fluid, for which the material is perfectly rigid at stresses below the yield stress. While a static unyielded plug forms at the tip of the corner, eddies analogous to those found by Moffatt can also form. We examine these viscoplastic eddies numerically, by computing finite element solutions using the augmented-Lagrangian method, and analytically, by employing a viscoplastic boundary layer formulation and scaling arguments. We measure the depth of the static plug as a function of the Bingham number (dimensionless yield stress), show that the process of a new eddy forming as the Bingham number is decreased is driven by the pressure in the yielded fluid adjacent to the static plug, and provide a heuristic argument for the critical Bingham number at which this occurs.

Investigating the graded-index influence on elastic responses of axisymmetric pressurized and heated thick-walled functionally graded material of cylindrical vessel

This article presents the influence of a graded-index on the elastic responses of an axisymmetric pressurized and heated thick-walled functionally graded material (FGM) cylindrical vessel under plane strain conditions. It is assumed that the inside surface/outside surface of a cylinder is made up of ceramic /metal. The temperature-independent thermo-mechanical parameters of an FGM cylinder are described by using the rule-of-mixture with the function of radius. The profiles of thermo-mechanical material characteristics are assumed to vary continuously and gradually along the thickness direction of a cylinder. The numerical solutions of the thick-walled FGM cylinder governing differential equations are examined by using the non-linear shooting and fourth-order Runge-Kutta approaches. The thermo-elastic solutions of the thick-walled FGM cylinder over the radial direction are computed using Jupiter, Python open-access software, and graphically reported. The analytical solutions (AS) of a ceramic cylinder are compared with numerical solutions (NS) of an FGM cylinder (n = 0) to validate the proposed method. The results reveal that the graded-index n has a higher influence on the profiles of dimensionless hoop stress, dimensionless axial stress, dimensionless radial displacement, hoop strain, radial strain, and thermo-mechanical material properties in the dimensionless radial direction of an FGM cylinder. Nevertheless, the graded-index has an insignificant impact on dimensionless radial stress across the dimensionless radial direction of an FGM cylinder. Based on these findings, it is recommended that graded-index is an important parameter which can be used to regulate the responses of FGM structures under thermal and mechanical loading conditions.

Hazards in mining: A novel model for the prediction of run-out distances in tailings dams using CFD

Tailings dams are potential facilities capable of failure, with the potential to generate high-impact environmental, economic, and social catastrophes given that they do collapse. In this study, numerical simulations of tailings as non-Newtonian fluids based on the Herschel–Bulkley model have been carried out to characterize the runout distances under different physical conditions. Using scaling arguments, a mathematical expression formed by three dimensionless groups including dimensionless yield stress, a modified Reynolds number , and the angle of the inclination of the system is proposed to explain the runout distance in tailings dams. The subsequent results show a good correlation between the variables, allowing us to deduce a high influence of the yield stress and the slope, albeit reflected both in the modified Reynolds number and the yield coefficient. Finally, a simple mathematical model is obtained from multidimensional analysis and the possibility of an abacus for regular intervals is proposed, which in turn supports an additional complement to the existing models for estimating the runout distances. • Prediction of dam breaking in mining. • Runout distances using CFD. • Dense suspensions in mining.

Effect of Mechanical Loadings on Two Unequal Slanted Cracks Length in Bi-Materials Plate

Although a lot of crack problems in bi-materials plate were previously treated, few solutions are available under mechanical loadings, arbitrary crack lengths and material combinations. In this paper the dimensionless stress intensity factors (SIFs) of two slanted cracks in the upper plate of bi-materials are considered under mechanical loadings with varying the crack length and material combinations systematically. In order to calculate the dimensionless SIFs accurately, the hypersingular integral equations (HSIEs) was formulated by using the modified complex potentials (MCP) function. The details numerical results of the dimensionless SIFs are given in tabular form and graphical presentations. Comparisons with the existing exact solutions show that the numerical results in this paper have high accuracy. Our results are described with clarifying the effect of the mechanical loadings, bi-elastic constant ratio and element size of cracks on the dimensionless SIFs.

Modeling and simulation of interstitial fluid flow around an osteocyte in a lacuno-canalicular network

Experiments have shown that external mechanical loading plays an important role in bone development and remodeling. In fact, recent research has provided evidence that osteocytes can sense such loading and respond by releasing biochemical signals (mechanotransduction, MT) that initiate bone degradation or growth. Many aspects on MT remain unclear, especially at the cellular level. Because of the extreme hardness of the bone matrix and complexity of the microenvironment that an osteocyte lives in, in vivo studies are difficult; in contrast, modeling and simulation are viable approaches. Although many computational studies have been carried out, the complex geometry that can involve 60+ irregular canaliculi is often simplified to a select few straight tubes or channels. In addition, the pericellular matrix (PCM) is usually not considered. To better understand the effects of these frequently neglected aspects, we use the lattice Boltzmann equations to model the fluid flow over an osteocyte in a lacuno-canalicular network in two dimensions. We focus on the influences of the number/geometry of the canaliculi and the effects of the PCM on the fluid wall shear stress (WSS) and normal stress (WNS) on an osteocyte surface. We consider 16, 32, and 64 canaliculi using one randomly generated geometry for each of the 16 and 32 canaliculi cases and three geometries for the 64 canaliculi case. We also consider 0%, 5%, 10%, 20%, and 40% pericellular matrix density. Numerical results on the WSS and WNS distributions and on the velocity field are visualized, compared, and analyzed. Our major results are as follows: (1) the fluid flow generates significantly greater force on the surface of the osteocyte if the model includes the pericellular matrix (PCM); (2) in the absence of PCM, the average magnitudes of the stresses on the osteocyte surface are not significantly altered by the number and geometry of the canaliculi despite some quantitative influence of the latter on overall variation and distribution of those stresses; and (3) the dimensionless stress (stress after non-dimensionalization) on the osteocyte surface scales approximately as the reciprocal of the Reynolds number and increasing PCM density in the canaliculi reduces the range of Reynolds number values for which the scaling law holds.

Experimental study on the critical condition of river blockage by a viscous debris flow

A debris flow flowing into a river may cause river blockage and form a dam. When the dam breaks, large floods can occur in the areas downstream, which can lead to fatalities and property damage. In this study, the parameters of river blockage were classified as three types of factors: morphological factors, dynamic factors, and characteristics factors of debris flows. The confluent angle, the dimensionless volume, the unit discharge ratio, and the dimensionless yield stress are important parameters for river blockage. They are all independent and dimensionless parameters. The degree of river blockage was divided into 4 stages: minor blockage, partial blockage, major blockage, and complete blockage. A classification based on the blockage coefficient is proposed to assess the degree of blockage. An experimental study was conducted with a total of 61 experimental runs. Based on the blockage coefficient of each run of the tests, the relationships between the degree of blockage and the key parameters (i.e., the confluent angle, the dimensionless volume, the unit discharge ratio, and the dimensionless yield stress) were obtained. A prediction model of the degree of river blockage is proposed based on these relationships. The area total of 21 debris flow tributaries and their corresponding rivers in China were selected to validate the prediction model. In these 21 tributaries, 24 events involving a debris flow flowing into a river resulted in partial blockage, major blockage, and completely blockage of the rivers. The prediction model was successful validated using the 22 events in China involving partial blockage, major blockage, and complete blockage.

Re-profiling response of heavily squeezed tunnels via ground reaction curves

This paper presents a simple method for determining the re-profiling response of a heavily squeezed tunnel using ground reaction curves. The ground reaction curve of certain tunnel is unique since the mechanical behavior of the original surrounding rock of each re-profiled tunnel remains the same. The dimensionless stress distribution also remains the same for a same support pressure. According to the corresponding relation of material points in the origin and current coordinates, the excavation radius of the re-profiled tunnel is determined. The location of re-profiled tunnel on the ground reaction curve and evolution curves of softening and residual radii can be determined, and the tunnel response induced by the decreasing support pressure can be further predicted using the above three curves. The critical excavation radius for a desired clearance and its corresponding plastic and residual radii can be theoretically determined. As the increase of re-profiling times and/or critical plastic shear strain, the increments of tunnel wall displacement, plastic and residual radii gradually decrease for a constant re-profiling radius.

Residual stress indentation model based on material equivalence

A novel residual stress indentation model for conical indentation loading is proposed to describe the relationship between the residual stress, material constitutive parameters, load, and displacement for materials with a uniaxial constitutive relationship that obeys Hollomon’s power law (H-law). The novel model was established based on the principle that the equivalent material without residual stress corresponds to the original material with residual stress, conical indentation theoretical model based on energy density equivalence, and an assumed power-law relationship between the dimensionless residual stress and relative difference of the yield stresses of the equivalent material and original material. Sixty imaginary H-law materials with ten equibiaxial and ten uniaxial residual stresses were investigated by Finite Element Analysis (FEA). The residual stresses predicted by the novel model from the indentation load–displacement curves simulated for the imaginary materials are in close agreement with those applied by the FEA. Finally, indentation tests for Cr12MoV steel, 45 steel, and 6061-T6511 aluminum alloy were carried out on their specimens without residual stress and their bending specimens with equibiaxial and uniaxial residual stresses. The residual stresses predicted by the novel model according to the indentation load–displacement test curves are in good agreement with those applied by the tests.

Simulation of time-dependent radiative heat motion over a stretching/shrinking sheet of hybrid Nanofluid: Stability analysis for dual solutions

The heat transfer characteristics for the flow of a time-dependent hybrid nanofluid with thermal radiation and source/sink over a stretching/shrinking sheet are examined in the current investigation. We have transformed the governing equations of the presented study into the similarity equations utilizing similarity variables. However, a numerical solution is obtained by using in-build MATLAB code bvp5c. The mass and energy profiles for diverse values of thermophysical parameters are studied together with their physical quantities. It is observed that dual solutions exist, that is, one is upper, and the other is lower branch solution for a definite choice of the unsteadiness parameter. Also, stability analysis is executed to determine the long-term stability of dual solutions, indicating that out of the two, only one is stable and the other is unstable. It is revealed that comparatively, the first solution shows stability, while the second solution shows instability. There is a considerable influence of second-order slip on the problem’s respective flow and heat transfer characteristics. Further, major outcomes also show the dimensionless frictional stress and the magnitude of conventional heat transfer enhancement with growing suction parameter values.

STUDYING THE EXTRUSION PROCESS OF A SINTERED THIN-WALL PIPE LOADED WITH INTERNAL HIGH PRESSURE IN CASE OF A BOUNDARY CONDITION IN THE DIE INLET

An analytical study of the process of sintered thin-wall pipe extrusion loaded with internal high pressure, on its cross section in the conical die inlet by compressing meridio¬n¬a¬l dimensionless stress has been carried out. It is assumed that the external die pressure p_v and the internal high pressure p_r acting on the pipe are equal. In this case, the circumferential stress is determined and a simplified plasticity condition is obtained. The dimensionless values of meridional, circumferential and average stresses depending on the degree of pipe deformation were determined, which made it possible to study the processes of their compaction using the DTPPM equations at different values of coefficients of friction and initial porosities of the pipe material. A computer simulation of the process was carried out in the MS EXEL software environment, presetting different initial values of the compressive meridional dimensionl¬e¬s¬s stress. The pipe extrusion process will correspond to such an initial value of this stress, when the meridional stress of the pipe in the outlet of the die is equal to zero. In this case, the value of strain intensities is ε_i=0.258. The graphs of changes in 10, 20 and 30% of the initial porosities of the material show that the porosity becomes zero in the interval of value of deformation intensities ε_i from 0.116 to 0.132, and also in the case of large values of porosities and contact friction, the compaction of the material occurs more intensively.

A FFM analysis on mode III static and fatigue crack initiation from sharp V-notches

The brittle or quasi-brittle failure initiation at sharp V-notches under mode III loading conditions is investigated by means of the coupled Finite Fracture Mechanics (FFM) approach. The model is developed by assuming that failure is shear stress governed, and by imposing a condition of consistency of both energy and stress requirements. The expression for the notch fracture toughness, which results a function of the material torsional strength, mode III fracture toughness and notch opening angle, is firstly derived. This allows to achieve the dimensionless failure stress through the analytical shape functions presented in the literature. By considering, as a general rule, the transition of fracture mechanism from tensile to shear in terms of the material strength ratio, experimental data referring to cylindrical V-notched samples are then discussed in the static regime. Although FFM predictions on the failure stress are fairly accurate, thus revealing simple tools for engineering purposes, a linear elastic analysis results simplistic since it disregards non linear effects (friction/ductility) observed during tests. Finally, the approach is extended to the fatigue framework, where the basic assumptions are proved to be more reliable. The analysis includes a discussion on the material properties to be implemented, and a satisfactory comparison with a large variety of experimental data.

Experimental study on Mode-I fracture toughness using Chevron Straight-notched Semi-circular Bend (CSNSCB) method

Rock fracture toughness is an important parameter for the engineering geological design. In this present study, chevron straight-notched semi-circular bend (CSNSCB) specimens with wide-flat liner chevron ligament (Type-A) and long-narrow liner chevron ligament (Type-B) were designed and prepared using wire cutting method. Mode-I fracture toughness KIC values and dimensionless stress intensity factor Y* were calculated. Moreover, these two types rock failure properties and KIC values variation features were both analyzed. The obtained results indicate that the width of precast artificial crack is small which can ensure the flatness and symmetry of cutting plane and then the tested Stress Intensity Factor (SIF) which has smaller deviation is more accurate and reliable. The Y* values for each type of configuration were calculated using Finite Element Method (FEM) and Slice Synthesis Method (SSM) and they firstly had a decreased trend and then an increased trend with the increasing of normalized crack length α, which represent the transition from stable crack propagation (α > αm) to unstable crack propagation (α < αm). Compared to Type-A specimens, Type-B specimens had the calculated minimum dimensionless stress intensity factor Y*min with smaller deviation and average KIC values with smaller variation coefficients. Moreover, Type-A specimen has the average fractal dimension D value of 1.2209 which is 4.79% larger than that of 1.1651 for Type-B specimen, indicating its more complex crack generation and rougher fracture surfaces. That is to say, Type-B specimen has the advantages of small crack FPZ and crack propagation without deviation as well as smooth fracture surface. Since it can offset the harmful influence of rock anisotropy and good crack guided effect, its KIC value is much satisfactory and closer to the true fracture toughness of rock-like materials.

Influence of porosity distribution on static and buckling responses of porous functionally graded plates

This article studies the effect of porosity distribution on static and buckling response of a functionally graded (FG) porous plate with all its four edges simply supported and subjected to a transverse load. The plate’s displacement field is approximated based on an inverse hyperbolic shear deformation theory (IHSDT) involving five variables. In this theory, inplane displacements vary as an inverse hyperbolic function through the plate thickness, so account for non-linear variation of transverse shear stresses and satisfies traction free boundary conditions on the surfaces of plate without the need of any shear correction factors. The governing equations of motion are derived using the principle of virtual work, and thus, the expressions of its strain energy and work done are obtained. Effective material properties of porous plate are assumed with an additional term of porosity in the direction of thickness. Power-law is used for continuous variation of material properties along the thickness direction. In this study, five types of porosity distribution functions are considered. The effect of the porosity parameter, the power-law exponent, side-thickness ratio, and aspect ratio on the static and buckling responses of the porous FG plate is evaluated. In problem-solving, the Navier solution technique is used to obtain the exact solutions. The non-dimensional deflections and critical buckling loads of FG porous plate are thus obtained. The results are compared with available results in the literature, and good agreement is found. Moreover, the results are presented to study the effects of various parameters on the dimensionless deflections, stresses, and critical uni-axial and bi-axial buckling loads.

Numerical exploration of thermal transport in water-based nanoparticles: A computational strategy

This report examines the three-dimensional flow characteristics and heat transfer of TiO2–Cu/water nanofluid past over a bidirectional surface with the effect of MHD. Tiwari and Das model of nanofluid is used with the thermophysical properties of nanofluid. The derived physical laws have been simplified under boundary layer theory. The complex sets of coupled PDEs have developed in form of ODEs by utilizing the appropriate transformations. The transformed equations handled via shooting method. The impacts of different influential parameters have been sketched. The involvement of different parameters for dimensionless stress, rate of heat and mass are tabulated. Limiting case of present investigation is compared with those of published ones and an excellent settlement is noted. It is noted that the enhancement in fluid parameters behave contrary to fluid velocity as compared with thermal and concentration profile. Moreover, thermal profile depreciate against higher values of Eckert and Prandtl numbers.

Experimental analysis of the pseudo-compact tension (pCT) testing configuration using two alternative sample geometries

The pseudo-compact tension (pCT) method recently proposed by Muñoz-Ibáñez et al. (2020) is a satisfactory approach to measure mode I fracture toughness (KIC) in rocks and other materials using disc-shaped samples loaded under pure tensile conditions. In contrast to other methods, such as the semi-circular bend (SCB) suggested by the ISRM (2014), the pCT test provides with good control after peak load, making it possible to further characterize the processes involved in fracture propagation. In this work we assess the influence of the testing configuration at the onset of unstable crack propagation. In order to extend the pCT concept to complementary geometries with potential interest we studied an alternative to the SCB specimen, which we call pseudo-SCB (pSCB). To compute KIC in this configuration we have derived the corresponding dimensionless stress intensity factor function (Y’) based on the finite element method. The results show that the pSCB test provides with consistent values of KIC and it also allows to control the propagation of the crack beyond peak load, which reinforces the idea that the loading conditions may be a more determinant factor than the sample geometry in controlling post-peak behaviour. In addition, an expression of Y’ is presented for cubic samples tested using the pCT approach. This configuration may be useful for testing other materials amenable of moulding such as mortar, concrete, ceramics, etc.

V-notch crack toughness of PMMA by strain gages and its comparison with caustics results under impact loading

Strain gages and dynamic caustics were used in this experiment to measure the V-notch crack initiation toughness and propagation toughness of polymethyl methacrylate (PMMA) specimens. Furthermore, a method of pasting strain gages on V-notch PMMA specimens for three-point bending test was proposed. The specimens were divided into five groups according to the opening angle (ω), and each group corresponded to two pasting methods: acute-angle pasting method and obtuse-angle pasting method. The results of notch stress intensity factors (NSIFs) obtained by strain gage measurement were consistent with the results calculated using the dynamic caustics method. The concept of the dimensionless notch stress intensity factor (DNSIF) was introduced for a better analysis of the effect of opening angles on cracking. The experimental results revealed that the initiation of V-notch crack in the specimens were inhibited more when ω was larger than 90° than for small opening angles. Under this condition, the specimens displayed higher cracking initiation toughness. Furthermore, the time required for crack initiation was longer, and the cracking was instantaneous. An analysis of the vicinity of the V-notch crack tip showed that with the increase in ω, the stress component decreased, and the rate of decrease became higher. This also explains the phenomenon wherein the specimens transformed from being easy to crack to not being so as ω increased. The experimental results showed that the crack propagation velocity and propagation toughness were also affected by the opening angle of the V-notch crack.

Effect of frequency of minute-amplitude oscillatory shear loadings on ultimate yield stress of cohesive sediments

The ultimate yield stress of cohesive sediments changes with oscillatory loading or vibration. Given a strong enough vibration, the sediment may be fluidized. In this study, the effect of the frequency of minute-amplitude oscillatory shear loadings on the ultimate yield stress of cohesive sediments was experimentally and analytically investigated in a laboratory. The test samples were prepared with cohesive sediments with a median particle size of 31 μm and at four water content (31.8% to 37.03%). The sediments were subjected to continuous oscillatory shear loadings with a constant frequency (0–131.6 Hz). The findings of the study indicated that the ultimate yield stress of cohesive sediments under minute-amplitude oscillatory shear loadings was dependent on the oscillatory shear frequency, water content, and the ultimate yield stress of the sediments prior to the application of the oscillatory shear loading. This study revealed that there exists a critical shear loading frequency fcr. When the oscillatory shear frequency (f) was less than the critical shear frequency (fcr), the ultimate yield stress decreased as f increased. However, when f was greater than fcr, the ultimate yield stress was constant, irrespective of further increase of the oscillatory shear frequency. The former condition is identified as the non-equilibrium fluidization stage and the latter, the equilibrium fluidization stage. The fcr value decreased as the water content increased. In the equilibrium fluidization stage, the dimensionless yield stress Su/Su0 did not vary with the water content or oscillatory shear frequency. Two discrete formulae were proposed to estimate the ultimate yield stress of cohesive sediments subjected to minute-amplitude oscillatory shear loadings for the non-equilibrium and equilibrium fluidization stages. To explore the engineering applications of the formula for the equilibrium fluidization stage, a series of vertical pull out tests were performed using a dynamic torpedo anchor. This study provides an important method for distinguishing the degree of fluidization and sheds lights on potential applications in coastal and marine engineering, such as pile or anchor installation and extraction and scour protection in cohesive bed.

The influence of the rigidity of mining security structures on the stability of side rocks in the coal-rock massif

The object of research is the processes of controlling the state of side rocks to prevent the collapse of the stratified rock strata in the coal-rock massif containing the workings. The studies carried out made it possible to establish the influence of the rigidity of the guard structures of mine workings on the stability of side rocks in the coal-rock massif. It is proved that as a result of the action of an external compressive load on the supporting lateral rocks, the structure, in the form of a model of bunches made of wooden posts, increases its rigidity until the destruction of the security structure. After that, there is an increase in the convergence of side rocks and their destruction. The change in the stiffness of the crushed rock in the filling massif model, which is used to support the lateral rocks, occurs as a result of the compaction of the original material. This is due to repackaging of crushed rock fractions of different sizes and its additional grinding. As a result of this interaction of the side rocks with the filling massif, the integrity of the roof and soil is ensured and convergence is limited. To assess the stability of side rocks, it is proposed to use a dimensionless stress concentration factor k. This coefficient takes into account the rigidity of the guard structures of the mine workings and the flexural rigidity of the side rocks. It was found that when the values of the coefficient k are close to zero (k→0), there is a loss of stability of the guard structures of mine workings and the destruction of side rocks in the coal-rock massif. The preservation of the integrity of the side rocks and the stability of security structures is ensured at values of k&gt;0.1, which corresponds to the parameters of the pliable supporting structures. Most favorably on the condition of side rocks in the coal-rock massif is influenced by the method of backing up the mined-out space of crushed rock. The use of this method excludes the collapse of side rocks. When solving the problem of stability of mine workings at the stage of making technical decisions, it is necessary to predetermine the issues of rigidity of security structures with deformation characteristics of side rocks.

A New Approach to the Modeling of Magnetorheological Dampers and Application to Resonance Control

Abstract A predetermined flow pattern in a magnetorheological damper providing continuously variable resistance to flow is required for efficient damping of a given load. The required predetermined flow pattern rests on the a priori determination of the constitutive properties of the magnetorheological (MR) fluid determined to generate variable resistance to flow. The inverse problem of constructing the predetermined response of the damper with a specific displacement pattern of the piston in the damper for efficient damping of a given load is solved. The MR fluid in the damper is modeled as a Bingham phase change material with time-dependent yield stress offering continuously variable resistance to the flow in the piston to achieve the required specific displacement pattern. The governing equations are solved for any time history of the dimensionless yield stress of the fluid which in turn is determined from the imposed response of the damper. Analytical tools developed can be used in optimizing damper performance. The application of the method to resonance mitigation is illustrated.

The Effect of Depth‐Dependent Stress in Controlling Free‐Surface‐Induced Supershear Rupture on Strike‐Slip Faults

AbstractSupershear rupture is crucial in determining ground motion characteristics and its various transition mechanisms is in key connection with rupture dynamics, geometrical complexities and stress heterogeneities. Using 3D dynamic rupture models with numerous depth‐dependent stress regimes, we show that depth‐dependent stress plays a major role in controlling the occurrence of free‐surface‐induced (FSI) supershear ruptures. The effect of stress gradient, initial normal stress at the free surface, fault width and critical depth of depth‐dependent stress regime are examined. The 2D tangent rupture velocity distributions show that generation of a sustained FSI supershear rupture prefers a larger stress gradient and initial normal stress at the free surface. The truncation of fault width could transit the rupture to a sub‐Rayleigh rupture. We propose a dimensionless stress parameter representing a measure of the overall level of depth‐dependence in the initial normal stress, and a larger moment magnitude is more likely to trigger a sustained FSI supershear rupture. By comparing the near‐field ground motions, three distinct ground motion characteristics are presented using the sustainability of Mach fronts and displacement polarization. This work helps to understand the occurrence of FSI supershear ruptures and its relationship with underground depth‐dependent stresses.