Dynamic Stress Concentration Research Articles

The dynamic performance of semi-infinite piezoelectric strip containing irregular boundaries based on the boundary element method

In this article, the problem of a semi-infinite piezoelectric strip with irregular boundaries is studied by boundary element method, and SH guided wave and line source loads are considered as external loads acting on piezoelectric strip. The superiority of the boundary element method is induced by two different arithmetic examples. Firstly, when the irregular boundary is not subjected to line source loads, the dynamic characteristics are analyzed in the first example, employing the image method and Graf addition theorem. Then, the Green’s identities are introduced and the Green’s function in infinite three-dimensional space is solved. When irregular boundaries are subjected to line source loads, the dynamic characteristics are investigated using the boundary element method in the second example. Finally, results clarified the influence on the dynamic stress concentration factor and electric field intensity concentration factor under proper conditions. Besides, the analytical solutions are compared with the finite element solutions to verify the accuracy of the conclusions in this article.

Mechanical response and failure mechanism of circular inclusion embedded in brittle materials under dynamic impact

Inclusions are prevalent in both natural and synthetic materials. Gaining a comprehensive understanding of their dynamic response and interaction with the surrounding matrix is essential in fundamental mechanics and engineering applications. This study aims to achieve such understanding through theoretical analysis, experimental investigations, and numerical simulations, focusing on the dynamic destruction of inclusions and the underlying mechanisms. Firstly, the circumferential, radial, and shear stresses around elastic heterogeneous inclusions were derived by the wave function expansion and Duhamel integration methods. Subsequently, a set of experiments on sandstone specimens containing three types of inclusions (plaster, epoxy resin, cement) were conducted utilizing the Split Hopkinson Pressure Bar (SHPB) system in conjunction with high-speed photography and Digital Image Correlation (DIC) techniques to obtain the surface deformation of inclusion specimens. Eventually, a series of Finite Element Method (FEM) numerical simulations adopted suitable materials were also carried out to investigate the fracture process of inclusion and matrix under dynamic impact. The results demonstrate that the stress distributions and fracture mode of matrix and inclusion are highly dependent on the physical and mechanical properties of the inclusions and the surrounding matrix, specifically, density ρ, elastic modulus E, Poisson's ratio v, and strength. With an increase in the E and v, there is a reduction in the concentration of circumferential stress, while the radial and shear stresses experience an increase. The experimental and numerical results corroborated the theoretical findings and indicated that the localized dynamic stress concentrations induced by wave scattering around the inclusions directly dominated both local and overall specimen failure. Due to the dissimilarities in elastic parameters and strength between the inclusion and the matrix, the stress state of the inclusion and the interface is not homogeneous. Under dynamic loading, the weaker inclusion experiences tensile cracking at both ends of the loading, and as the mechanical properties (ρ, E, v, and strength) of the inclusion rise, a transition from tensile- hybrid-shear failure occurs within the inclusion. The findings of this study help us understand the dynamic failure mechanisms of dissimilar inclusions in brittle material.

Response of semicircular canyons and movable cylindrical cavities to SH waves in anisotropic half-space geology

The development of tunnels or the laying of underground pipelines are essential engineering projects in modern society, and in canyon tunnels and underground pipeline projects, the surface motion and cavity edge motion have been topics of concern in ground vibration problems. We investigate the wave scattering problem in an elastic half-space anisotropic medium containing a semicircular canyon and a subsurface movable cylindrical cavity by using the wave function expansion method, the complex function method, and the mirror method. By deriving the governing equation and transforming it into the standard form of the Helmholtz equation satisfying the zero-stress boundary condition, we solve the corresponding displacement functions. Introducing a position correction coefficient, the scattered wavefield in a half-space anisotropic medium is constructed by the mirror method, which improves the problem of scattered wave source singularity in anisotropic half-space media. Then, combining the free boundary conditions with a Fourier series expansion method, we solve the unknown coefficients in the equations. The correctness of the method is verified by degenerating it to a classical analytic solution. Finally, using frequency- and time-domain analysis, we investigate the effects of the relevant parameters on the surface motion [Formula: see text], the dynamic stress concentration factor, and the displacement amplitude [Formula: see text]. The results indicate that rock anisotropy and the presence of semicircular canyons have a significant effect on the dynamic response of subsurface structures. This not only provides a theoretical basis for practical unlined tunnels or pipeline projects but also can provide a basis for seismic design of underground structures.

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.

Blasting induced dynamic stress concentration and failure characteristics of deep-buried rock tunnel

Blasting induced dynamic stress concentration and failure characteristics of deep-buried rock tunnel

Analytical analysis of the interaction between cavities and a crack in bonded isotropic and anisotropic half spaces under SH wave

Interaction between two cavities and a linear crack subjected to elastic shear horizontally (SH) wave is investigated analytically in this work. Two cavities are buried in two bounded half spaces respectively and the cross-section of each cavity has a different shape. The media of two half spaces are isotropic and anisotropic respectively, and the complex function method is applied to first solve the wave equation in each medium. Then, with the help of image principle, Green’s functions and wave field expressions in each separated half space are constructed. Subsequently, utilizing “conjunction” technique, a series of Fredholm integral equations containing undetermined forces at the interface are obtained and reduced in order to be solved. Dynamic stress distribution is mainly considered to discuss the interaction between the cavities and the crack. Finally, a parametric study on dynamic stress concentration factor (DSCF) and dynamic stress intensity factor (DSIF) presents the interaction is affected by several dimensionless factors simultaneously.

Influence of Incident Orientation on the Dynamic Response of Deep U-Shaped Cavern Subjected to Transient Loading

In deep rock engineering, caverns are often disturbed by engineering loads from different directions. To investigate the dynamic response of deep U-shaped caverns under different incident orientations, a theoretical solution of the dynamic stress concentration factor along the cavern boundary was derived based on the wave function expansion and conformal mapping method, and the failure characteristics around the cavern were further investigated by PFC2D (Particle Flow Code in two dimensions). As the incident orientation increases from 0° to 90°, the dynamic compressive stress concentration area transforms from both the roof and the floor to the sidewalls, and the peak dynamic stress concentration factor of the roof decreases from 2.98 to −0.20. The failure of the floor converts from dynamic compression shear failure to dynamic tensile failure. Compared to a stress wave incident from the curved boundary, a stress wave incident from the flat boundary causes severer damage. When the stress wave is incident from the sidewall, the cavern with a larger height-to-width (h/w) ratio exhibits severer damage. Conversely, the cavern with a smaller h/w ratio tends to fail as the stress wave is incident from the floor. This paper provides a basic understanding of dynamic responses of the deep U-shaped cavern.

Scattering of SH waves by elliptical cavity and type-III crack in deep anisotropic geology

Most of the underground rock formations are anisotropic, and the study of the dynamic response of anisotropic geology containing internal defects is necessary for the development of underground tunnels or caverns. In this paper, the wave function expansion method, the complex function method, and the Green's function method are used to solve the problem of scattering of SH waves in anisotropic rock in the presence of elliptical cavity and type-III crack. Boundary conditions for solving the problem are given using the conformal transformation method, and crack in the medium is constructed using the "crack cutting" technique. The scattered wave field generated by the crack is constructed based on the Green function method, and the displacement and stress fields in the model are obtained by the wave field superposition principle in the case of SH wave incidence. The dynamic stress concentration factor (DSCF) on the boundary of the elliptical cavity and the dynamic stress intensity factor (DSIF) at the crack tip are derived. The correctness of the method is verified by degenerating it to a classical analytic solution. Finally, the effects of the anisotropic strength of the rock mass, the dimensions of the cavern, and the incidence angle and wave number of the SH wave on the DSCF and DSIF are analyzed by numerical calculations in both the frequency and time domains. The results show that the anisotropy of the rock mass, the shape of the cavern and the defects of the surrounding medium have a great influence on the dynamic response of the underground structures, which should be given enough attention in the development and design of deep underground caverns or tunnels.

Dynamic stress prediction for a Pump-Turbine in Low-Load Conditions: Experimental validation and phenomenological analysis

To meet the increasing demand for higher flexibility in hydropower, an operating range extension in an existing hydroelectric powerplant is considered. To this end, the dynamic flow and structural behavior of the low-head pump-turbine is investigated in turbine mode in partial load and deep partial load conditions, to evaluate the feasibility of this flexibilization due to the increased high cycle fatigue damage. In a first step, different numerical simulation approaches that combine fluid dynamics and structural dynamics are applied and compared to in-situ dynamic strain and pressure measurements on the impeller. As a result, the modelling requirements for a numerical workflow that leads to unmatched accuracy in predicting the broad-band dynamic stresses of deep partial load are identified. In a second step, the operating-point-dependent fatigue crack initiation spots are identified based on dynamic stresses. Herein, the emergence of an additional dominant critical spot near the blade leading edge is observed, which is unexpected for low-head pump turbines. Thanks to the detailed numerical flow and structural analyses, the root cause for the dynamic stress concentration is studied and the emergence of the new critical spot is understood. It turns out in deep partial load, that the main source of pressure oscillations on the impeller blade surfaces is the vortex-dominated flow detachment zone around the blade leading edge, which, in combination with the geometry of the investigated impeller, translates into the main stress oscillations. The results of this work were later used to extend the operating range of the hydroelectric powerplant from initially 25% to 100% power output to now 0% to 100%.

Dynamic inverse plane characteristics of circular holes in a strip-shaped medium under line source excitation

Based on the complex variable function method and the wave function expansion method, the dynamic anti-plane characteristics of circular holes in the strip domain under line source excitation were studied. The total wave field in the strip domain is constructed using the concept of infinite mirroring, which eliminates the mutual impact of the wave field under the action of the upper and lower limits of the strip domain. To determine the overall wave field, solve the wave field's unknown coefficients in accordance with the boundary conditions. The right number of mirror images is found using the dynamic stress distribution of a circular hole in the band-shaped domain as a specific calculation example. The effect of altering the incident frequency and band-shaped domain width on the dynamic stress concentration coefficient surrounding the circular hole is also explored. As a consequence of the findings, the dynamic stress concentration coefficient around the circular hole in the zone domain at low frequencies is substantially bigger than it is at medium and high frequencies. Under identical conditions, the dynamic stress concentration coefficient will exhibit periodic extreme values and converge as the distance increases between the lower boundary and the circular hole’s center.

Anti‐plane stress analysis for V‐notch in a piezomagnetic half space under SH wave

AbstractIn this paper, the anti‐plane stress analysis of a V‐notch with complex boundary conditions in a piezomagnetic half space is studied. Firstly, SH wave is considered as an external load acting on piezomagnetic half space, on the basis of repeated image superposition, the analytical expression of scattering wave is conducted, which satisfies the boundary conditions on the boundary of the half space. Then, the analytical expression of standing wave is established, which satisfies the stress free and magnetic insulation conditions on the boundaries of V‐notch by the fractional Bessel function expansion method and Graf addition theorem. Finally, Green's function method is applied, the half space is divided into two parts along the vertical interface, a pair of in‐plane magnetic field and out‐plane forces are applied on the vertical interface, and the first kind of Fredholm integral equations are set up and solved by applying orthogonal function expansion technique and effective truncation. Results clarified the influence on the dynamic stress concentration factor and magnetic field intensity concentration factor under proper conditions. Besides, the analytical solutions are compared with the finite element solutions to verify the accuracy of the conclusions in this article.

Dynamic anti-planar characteristics of circular holes in laminated structures under the action of line source

Based on the complex variable function method and the wave function expansion method, the dynamic anti-plane characteristics of the circular hole in the laminated structure under line source excitation were studied. To deal with the interface problem within the laminated structure, the idea of “fit” is adopted to construct continuity conditions that satisfy displacement and stress. The circular hole boundary stress-free condition is used, the wave field is constructed using the mirror method, and the unknown coefficients of the scattered wave field are solved through Fourier series expansion, and finally the total wave field in the laminated structure is obtained. Taking the dynamic stress distribution of a circular hole in a laminated structure as a specific calculation example, the effects of changes in line source position, layer thickness, and medium hard-soft ratio on the dynamic stress concentration coefficient around the circular hole are discussed. The results show that when the line source load is applied directly above the circular hole and the lower layer is narrow, the dynamic stress concentration coefficient will have a sudden change compared with other situations. When other conditions are the same, under low, medium, and high-frequency wave fields, the dynamic stress concentration caused by different soft-to-hard ratios shows completely different rules. These rules can be used to better detect defects in laminated structures.

Dynamic response of semi-cylindrical depression, cylindrical cavity and type-III crack to SH wave in half-space anisotropic media

In this study, the anti-plane dynamic response of an elastic half-space anisotropic medium containing surface semi-cylindrical depressions and internal cylindrical cavity and type-III crack is solved analytically. The wave function expansion method, the complex function method and the Green's function method can be used to effectively construct the free wave field equation and the scattered wave field equation in the case of SH wave incidence, in which the scattered wave field generated by the crack is constructed by the Green's function method and the "crack cut" method. For the scattered wave field generated by a circular cavity is constructed using the mirror method, where the positional offset coefficient in the mirror method of a half-space anisotropic medium is introduced, which well solves the problem of the asymmetry of the scattered wave field. Finally, the effects of anisotropic strength of the medium, SH wave incidence angle and dimensionless incident wave number on the surface displacement amplitude (|W|), dynamic stress concentration factor (DSCF) at the boundary of the circular cavity and dynamic stress intensity factor (DSIF) at the crack tip are investigated in both frequency and time domains under different working conditions, and conclusions different from those of isotropic medium are obtained, which provide some theoretical references for the earthquake-resistant design of the underground structure with canyon topography and the surface building.

Harmonic wave scattered by an inclusion in an elastic plane: The complete Gurtin-Murdoch model

We consider the propagation of a harmonic elastic wave in a composite inclusion–matrix structure subjected to plane deformation. The interface between the inclusion and matrix is described by the complete Gurtin-Murdoch model with non-vanishing interface tension and interface stretching rigidity. We consider an inclusion of general shape and formulate the corresponding boundary value problem for the wave functions in the inclusion and matrix when a harmonic compressional or shear wave is incident on the edge of the matrix. The problem is then solved by series expansion methods for the case of a circular inclusion embedded in an infinite matrix. The series solution obtained is validated by checking its convergence and via comparisons with existing static and dynamic solutions for certain reduced cases. Numerical examples are presented for the case of a small-sized circular hole embedded in a soft matrix demonstrating the influence of surface tension on the incident wave-induced dynamic stress concentration in the matrix. We find that the presence of surface tension relieves the peak stress around the circular hole when the frequency of the incident wave is below a certain critical value, while it tends to intensify the peak stress for a high-frequency incident wave.

Analysis on scattering characteristics of SH guided wave due to V-notch in a piezoelectric bi-material strip

In this paper, the problem of a V-notch with complex boundary conditions in a piezoelectric bi-material strip is studied. First, SH guided wave is considered as an external load acting on the piezoelectric bi-material strip; on the basis of repeated image superposition, the analytical expression of scattering wave is conducted, which satisfies the stress-free and electric insulation conditions on the upper and lower horizontal boundaries of the strip. Then, the analytical expression of standing wave is established, which satisfies the stress-free and electric insulation conditions on the boundaries of V-notch by the fractional Bessel function expansion method and Graf addition theorem. Finally, Green’s function method is applied, the bi-material strip is divided into two parts along the vertical interface, a pair of in-plane electric field and out-plane forces is applied on the vertical interface, and the first kind of Fredholm integral equations is set up and solved by applying the orthogonal function expansion technique and effective truncation. Results clarified the influence on the dynamic stress concentration factor and electric field intensity concentration factor under proper conditions. Besides, the analytical solutions are compared with the finite element solutions to verify the accuracy of the conclusions in this article.

Analysis of dynamic stress concentration in three different types of poro-viscoelastic rock medium

Analysis of dynamic stress concentration in three different types of poro-viscoelastic rock medium

Studying the influence of heat treatment modes on the properties of medium-carbon steel

Problem. Parts of modern machines and structures operate under high dynamic loads, stress concentrations, and low or high temperatures. All of this contributes to brittle fracture and reduces the reliability of machines. Therefore, structural steels must, in addition to high mechanical properties, have high structural strength, which is manifested in the conditions of its actual use in the form of parts, structures, etc. [1] A properly selected heat treatment mode allows to obtain an optimal combination of properties in a part, thereby ensuring its reliability and operational durability.
 Goal. The aim of the study is to determine the effect of heat treatment parameters on the hardness, impact strength, and structure of 40 and 40X medium-carbon steels. Method. In the course of the work, a set of mechanical studies was carried out on samples of 40 and 40X steels in the initial state and after various heat treatment modes. Results. Based on the analysis of mechanical studies, the influence of heat treatment parameters on the mechanical properties and structure of medium-carbon steels 40 and 40X was studied. On the basis of mechanical and microstructural studies, the relationship between the effect of economical alloying of medium-carbon steel and the main indicators of mechanical properties in the range of heating temperatures was established. The experiments conducted in this work allowed us to obtain a comparative characterization of the effect of various technological parameters of heat treatment on the structure and properties of medium-carbon steels 40 and 40X, which are being improved. The factors that may affect the quality of parts made from these steel grades were also identified. Summarizing the results of metallographic studies and mechanical properties, the following conclusions can be drawn: 1. Steel 40, which has low supercooled austenite stability, must be cooled sharply (cold water, 8-12% aqueous solutions of NaOH or NaCl, etc.) to obtain continuous hardenability of parts. 2. The addition of 1 % chromium increases the stability of supercooled austenite, i.e. increases the hardenability of steel, so parts made of economical 40X alloy steel can be cooled in oil. However, if the part is large, a higher cooling rate should be used – in water. 3. Normalizing can be recommended for 40 and 40X steels as a work hardening heat treatment that increases ductility and impact strength. However, if the hardness level of the normalized alloy steel is high, then full or isothermal annealing should be used. 4. To prevent decarburization of steel at high temperatures and long holding times, protective atmospheres must be used. In the case of processing parts in furnaces with an oxidizing atmosphere, the parts must be filled with graphite, because cast iron chips do not provide reliable protection against oxidation. 5. In the range of tempering temperatures of improving steels recommended by the reference literature (400-600 °C), a monotonic decrease in hardness of both steels occurs, but in carbon steels, the hardening is more intense than in alloy steels. 6. To prevent the development of tempering brittleness of the second kind, to which 40X steel is prone, cooling during tempering should be accelerated (in oil, water), especially if the tempering temperature exceeds 550 °C.

Metamaterial invisibility of a corner target with proportional integral active control

Wave invisibility in confined space has been proved practical significance in mechanical engineering and military fields. In this work, a flexural wave cloak of a corner-domain target in a metamaterial plate with proportional integral (PI) active control circuit is proposed. According to the double reflection principle and mechanical characteristic at the boundary, wave fields for a corner-domain target are defined as a collection of incident, scattering and three mirror domains. The unknown scattering mode coefficients are determined by the image method and boundary conditions. The dynamic stress concentration, scattering amplitude and scattering cross-section for the corner-domain target are obtained. The cloaking configuration is composed of multiple equally spaced concentric layers, which are filled with piezoelectric (PZT) patches. In addition, these PZT patches are connected to the active control circuit with proportional integral (PI) control module. The results indicate that the flexural wave cloak can reduce the dynamic stress concentration, scattering amplitude and scattering cross-section in corner domain. Compared to the original structure without PI active control, the effective frequency range of the active cloak can be extended by 20%.

Dynamic Responses of a Shallow Lined Tunnel with Imperfect Interface Under Transient P Wave

Shallow-buried cavities and underground structures are inevitably subjected to transient perturbation from earthquakes or engineering activities during their life circle, leading to the occurrence of dynamic failure. The investigation presented by this paper concentrates on the dynamic responses triggered by an aperiodic transient [Formula: see text] wave around a circular lined tunnel with an imperfect interface in a half-space and the spring model is employed to describe the interface between rock mass and liner. Resorting to the wave function expansion scheme, the large arc assumption, the Graf’s addition theorem and the numerical inversion algorithm of Laplace transformation, the dynamic responses around a circular shallow lined tunnel with an imperfect interface are evaluated theoretically. The good agreement of time–history curves and contour snapshots for stress between the theoretical solutions and the numerical counterparts validates the analytical scheme. The effect of disturbance direction, embedment depth, thickness of liner and interface coefficients on the dynamic responses around the shallow-buried lined tunnel is analyzed. The results indicate that the positions of peak hoop dynamic stress concentration factor (DSCF) are approximately perpendicular to the perturbation direction. The DSCF distributions converge to be stable when the ratio of embedment depth to excavation radius exceeds 10.0. The thickness of the liner mainly affects the dynamic responses around the interface. Increase of elastic coefficient induces a decreasing DSCF in rock mass and an increasing counterpart in the liner. The results may provide the theoretical basis for the transient responses and aseismic design of underground tunnels.

Dynamic Stress Concentration and Particle Velocity Response of a Composite Lining Tunnel with Imperfect Interfaces Subjected to Blasting P Waves

Dynamic Stress Concentration and Particle Velocity Response of a Composite Lining Tunnel with Imperfect Interfaces Subjected to Blasting P Waves