Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Journal of Mechanics and Physics of Solids42, 813–843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis et al. (Journal of the Mechanics and Physics of Solids43, 319–341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice et al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity.     

Wave propagation in transversely isotropic porous piezoelectric materials

Wave propagation in porous piezoelectric material (PPM), having crystal symmetry 6 mm, is studied analytically. Christoffel equation is derived for the propagation of plane harmonic waves in such a medium. The roots of this equation give four complex wave velocities which can propagate in such materials. The phase velocities of propagation and the attenuation quality factors of all these waves are described in terms of complex wave velocities. Phase velocities and attenuation of the waves in PPM depend on the phase direction. Numerical results are computed for the PPM BaTiO₃. The variation of phase velocity and attenuation quality factor with phase direction, porosity and the wave frequency is studied. The effects of anisotropy and piezoelectric coupling are also studied. The phase velocities of two quasi dilatational waves and one quasi shear waves get affected due to piezoelectric coupling while that of type 2 quasi shear wave remain unaffected. The phase velocities of all the four waves show non-dispersive behavior after certain critical high frequency. The phase velocity of all waves decreases with porosity while attenuation of respective waves increases with porosity of the medium. The characteristic curves, including slowness curves, velocity curves, and the attenuation curves, are also studied in this paper.     

Dynamic interaction of various beams with the underlying soil – finite and infinite, half-space and Winkler models

Various beams lying on the elastic half-space and subjected to a harmonic load are analyzed by a double numerical integration in wavenumber domain. The compliances of the beam–soil systems are presented for a wide frequency range and for a number of realistic parameter sets. Generally, the soil stiffness G has a strong influence on the low-frequency beam compliance whereas the beam parameters EI and m^′ are more important for the high-frequency compliance. An important parameter is the elastic length l=(EI/G)^1/4 of the beam–soil system. Around the corresponding frequency ω_l=v_S/l, the wave velocity of the combined beam–soil system changes from the Rayleigh wave v_R≈v_S to the bending wave velocity v_B and the combined beam–soil wave has typically a strong damping. The interaction frequency ω_l is found not far from the characteristic frequency ω₀=(G/m^′)^1/2 where an amplification compared to the static compliance is observed for special parameter constellations. In contrast, real foundation beams show no resonance effects as they are highly damped by the radiation into the soil. At medium and high frequencies, asymptotes for the compliance of the beam–soil system are found, u/P(ρv_Paiω)^−3/4 in case of the dominating damping and u/P(−m^′ω²)^−3/4 for high frequencies. The low-frequency compliance of the coupled beam–soil system can be approximated by u/P1/Gl, but it also depends weakly on the width a of the foundation. All numerical results of different beam–soil systems are evaluated to yield a unique relation u/P₀=f(a/l). The integral transform method is also applied to ballasted and slab tracks of railway lines, showing the influence of train speed on the deformation of the track beam. The presented results of infinite beams on half-space are compared with results of finite beams and with infinite beams on a Winkler support. Approximating Winkler parameters are given for realistic foundation-soil systems which are useful when vehicle-track interaction is analyzed for the prediction of railway induced vibration.     

Lock-on characteristics of a cavity shear layer

This investigation focuses on defining the lock-on regions of a cavity shear layer subject to local periodic excitations. A circular cylinder of small diameter (d=4 mm), located very close to the upstream edge of cavity, is used to generate the local periodic excitations in the form of oscillatory rotation about its center with various angular amplitudes (Δθ) and frequencies (f_e). All the experiments were conducted in a recirculating water channel at three different Reynolds numbers that are based on the momentum thickness at the upstream edge of cavity (Re_θ0=152, 216 and 278). The LDV system and the laser-sheet technique are employed to perform the quantitative velocity measurements and the qualitative flow visualization, respectively. For cavity flows at three Reynolds numbers studied, the resonant lock-on is found to be the primary lock-on region within the range of frequency ratio (f_e/f₀=0.28–2.0). Here f₀ denotes the natural instability frequency of an unexcited cavity shear layer. The frequency bandwidth of resonant lock-on region does increase with increasing excitation amplitudes (Δθ). While the excitation amplitudes are smaller than 5° (Δθ5°), the resonant lock-on region, at Reynolds numbers 216 and 278, distributes asymmetrically about f_e/f₀=1.0 and biases to the high frequency (or large f_e/f₀) side. However, the sidewise expansion of resonant lock-on region is enlarged and the degree of asymmetric distribution is alleviated at large excitation amplitudes (Δθ>5°). The amount of sidewise expansion of the resonant lock-on region biased toward the high-frequency side is more significant at the lowest Reynolds number (152) than those at two higher Reynolds numbers (216 and 278). Besides, there exists a sub-harmonic lock-on region only at the lowest Reynolds number 152. The existence of a sub-harmonic lock-on region clearly reveals that the differential equation governing the self-excited oscillation within a cavity contains the quadratic nonlinear term. Further, at the lowest Reynolds number (152), the sidewise expansion of the sub-harmonic lock-on region is much narrower than that of the resonant lock-on region.     

On the propagation waves in the theory of thermoelasticity with microtemperatures

The present paper studies the propagation of plane time harmonic waves in an infinite space filled by a thermoelastic material with microtemperatures. It is found that there are seven basic waves traveling with distinct speeds: (a) two transverse elastic waves uncoupled, undamped in time and traveling independently with the speed that is unaffected by the thermal effects; (b) two transverse thermal standing waves decaying exponentially to zero when time tends to infinity and they are unaffected by the elastic deformations; (c) three dilatational waves that are coupled due to the presence of thermal properties of the material. The set of dilatational waves consists of a quasi-elastic longitudinal wave and two quasi-thermal standing waves. The two transverse elastic waves are not subjected to the dispersion, while the other two transverse thermal standing waves and the dilatational waves present the dispersive character. Explicit expressions for all these seven waves are presented. The Rayleigh surface wave propagation problem is addressed and the secular equation is obtained in an explicit form. Numerical computations are performed for a specific model, and the results obtained are depicted graphically.     

Oblique stagnation-point flow of a viscoelastic fluid with heat transfer

The two-dimensional forced convection stagnation-point flow and heat transfer of a viscoelastic second grade fluid obliquely impinging on an infinite plane wall is considered as an exact solution of the full partial differential equations. This oblique flow consists of an orthogonal stagnation-point flow to which a shear flow whose vorticity is fixed at infinity is added. The relative importance of these flows is measured by a parameter γ. The viscoelastic problem is reduced to two ordinary differential equations governed by the Weissenberg number W_e, two parameters α and β, the later being a free parameter β, introduced by Tooke and Blyth [A note on oblique stagnation-point flow, Physics of Fluids 20 (2008) 033101-1–3], and the Prandtl number Pr. The two cases when α=β and α≠β are, respectively, considered. Physically the free parameter may be viewed as altering the structure of the shear flow component by varying the magnitude of the pressure gradient. It is found that the location of the separation point x_s of the boundary layer moves continuously from the left to the right of the origin of the axes (x_s<0).     

Crack repair using an elastic filler

The effect of repairing a crack in an elastic body using an elastic filler is examined in terms of the stress intensity levels generated at the crack tip. The effect of the filler is to change the stress field singularity from order 1/r^1/2 to 1/r^(1-λ) where r is the distance from the crack tip, and λ is the solution to a simple transcendental equation. The singularity power (1-λ) varies from (the unfilled crack limit) to 1 (the fully repaired crack), depending primarily on the scaled shear modulus ratio γ_r defined by G₂/G₁=γ_rε, where 2πε is the (small) crack angle, and the indices (1, 2) refer to base and filler material properties, respectively. The fully repaired limit is effectively reached for γ_r≈10, so that fillers with surprisingly small shear modulus ratios can be effectively used to repair cracks. This fits in with observations in the mining industry, where materials with G₂/G₁ of the order of 10^-3 have been found to be effective for stabilizing the walls of tunnels. The results are also relevant for the repair of cracks in thin elastic sheets.     

Piezomagnetic and piezoelectric poling effects on mode I and II crack initiation behavior of magnetoelectroelastic materials

The ferrite and ferroelectric phase of magnetoelectroelastic (MEE) material can be selected and processed to control the macroscopic behavior of electron devices using continuum mechanics models. Once macro- and/or microdefects appear, the highly intensified magnetic and electric energy localization could alter the response significantly to change the design performance. Alignment of poling directions of piezomagnetic and piezoelectric materials can add to the complexity of the MEE material behavior to which this study will be concerned with.Appropriate balance of distortional and dilatational energy density is no longer obvious when a material possesses anisotropy and/or nonhomogeneity. An excess of the former could result in unwanted geometric change while the latter may lead to unexpected fracture initiation. Such information can be evaluated quantitatively from the stationary values of the energy density function dW/dV. The maxima and minima have been known to coincide, respectively, with possible locations of permanent shape change and crack initiation regardless of material and loading type. The direction of poling with respect to a line crack and the material microstructure described by the constitutive coefficients will be specified explicitly with reference to the applied magnetic field, electric field and mechanical stress, both normal and shear. The crack initiation load and direction could be predicted by finding the direction for which the volume change is the largest. In contrast to intuition, change in poling directions can influence the cracking behavior of MEE dramatically. This will be demonstrated by the numerical results for the BaTiO₃–CoFe₂O₄ composite having different volume fractions where BaTiO₃ and CoFe₂O₄ are, respectively, the inclusion and matrix.To be emphasized is that mode I and II crack behavior will not have the same definition as that in classical fracture mechanics where load and crack extension symmetry would coincide. A striking result is found for a mode II crack. By keeping the magnetic poling fixed, a reversal of electric poling changed the crack initiation angle from θ₀=+80° to θ₀=−80° using the line extending ahead of the crack as the reference. This effect is also sensitive to the distance from the crack tip. Displayed and discussed are results for r/a=10⁻⁴ and 10⁻¹. Because the theory of magnetoelectroelasticity used in the analysis is based on the assumption of equilibrium where the influence of material microstructure is homogenized, the local space and temporal effects must be interpreted accordingly. Among them are the maximum values of (dW/dV)_max and (dW/dV)_min which refer to as possible sites of yielding and fracture. Since time and size are homogenized, it is implicitly understood that there is more time for yielding as compared to fracture being a more sudden process. This renders a higher dW/dV in contrast to that for fracture. Put it differently, a lower dW/dV with a shorter time for release could be more detrimental.     

Crack size and speed interaction characteristics at micro-, meso- and macro-scale

The motivation to examine physical events at even smaller size scale arises from the development of use-specific materials where information transfer from one micro- or macro-element to another could be pre-assigned. There is the growing belief that the cumulated macroscopic experiences could be related to those at the lower size scales. Otherwise, there serves little purpose to examine material behavior at the different scale levels. Size scale, however, is intimately associated with time, not to mention temperature. As the size and time scales are shifted, different physical events may be identified. Dislocations with the movements of atoms, shear and rotation of clusters of molecules with inhomogeneity of polycrystals; and yielding/fracture with bulk properties of continuum specimens. Piecemeal results at the different scale levels are vulnerable to the possibility that they may be incompatible. The attention should therefore be focused on a single formulation that has the characteristics of multiscaling in size and time. The fact that the task may be overwhelmingly difficult cannot be used as an excuse for ignoring the fundamental aspects of the problem.Local nonlinearity is smeared into a small zone ahead of the crack. A “restrain stress” is introduced to also account for cracking at the meso-scale.The major emphasis is placed on developing a model that could exhibit the evolution characteristics of change in cracking behavior due to size and speed. Material inhomogeneity is assumed to favor self-similar crack growth although this may not always be the case. For relatively high restrain stress, the possible nucleation of micro-, meso- and macro-crack can be distinguished near the crack tip region. This distinction quickly disappears after a small distance after which scaling is no longer possible. This character prevails for Mode I and II cracking at different speeds. Special efforts are made to confine discussions within the framework of assumed conditions. To be kept in mind are the words of Isaac Newton in the Fourth Regula Philosophandi:

“Men are often led into error by the love of simplicity which disposes us to reduce things to few principles, and to conceive a greater simplicity in nature than there really is … We may learn something of the way in which nature operates from fact and observation; but if we conclude that it operates in such a manner, only because to our understanding that operates to be the best and simplest manner, we shall always go wrong.”––Isaac Newton

Non-linear oscillation of circular cylindrical shells

The method of multiple scales is used to analyze the non-linear forced response of circular cylindrical shells in the presence of a two-to-one internal (autoparametric) resonance to a harmonic excitation having the frequency Ω. If ω_r and a_r denote the frequency and amplitude of a flexural mode and ω_b and a_b denote the frequency and amplitude of the breathing mode, the steady-state response exhibits a saturation phenomenon when ω_b ≈ 2ω_r, if the excitation frequency Ω is near ω_b. As the amplitude ƒ of the excitation increases from zero, a_b increases linearly whereas a_r remains zero until a threshold is reached. This threshold is a function of the damping coefficients and ω_b−2ω_r. Beyond this threshold a_b remains constant (i.e. the breathing mode saturates) and the extra energy spills over into the flexural mode. In other words, although the breathing mode is directly excited by the load, it absorbs a small amount of the input energy (responds with a small amplitude) and passes the rest of the input energy into the flexural mode (responds with a large amplitude). For small damping coefficients and depending on the detunings of the internal resonance and the excitation, the response exhibits a Hopf bifurcation and consequently there are no steadystate periodic responses. Instead, the responses are amplitude- and phase-modulated motions. When Ω ≈ ω_r, there is no saturation phenomenon and at close to perfect resonance, the response exhibits a Hopf bifurcation, leading again to amplitude- and phase-modulated or chaotic motions.     

Rayleigh waves obtained by the indeterminate couple-stress theory

Rayleigh waves in a linear elastic couple-stress medium are investigated; the constitutive equations involve a length parameter l that characterizes the microstructure of the material. With , c_T=conventional transversal speed and q=wave number, an explicit expression is derived for the relation between , lq and Poisson's ratio ν. The Rayleigh speed turns out to be dispersive and always larger than the conventional Rayleigh speed. It is of interest that when lq=1 and ν≥0, it always holds that . The displacement field is investigated and it is shown that no Rayleigh wave motions exist when lq→∞ and when lq=1, ν≥0. Moreover, a principal change of the displacement field occurs when lq passes unity. The peculiarity that no Rayleigh wave motions exist when lq=1, ν≥0 may support the criticism by Eringen (1968) against the couple-stress theory adopted here as well as in much recent literature.     

Mechanics of adhesive contact on a power-law graded elastic half-space

We consider adhesive contact between a rigid sphere of radius R and a graded elastic half-space with Young's modulus varying with depth according to a power law E=E₀(z/c₀)^k (0<k<1) while Poisson's ratio ν remaining a constant. Closed-form analytical solutions are established for the critical force, the critical radius of contact area and the critical interfacial stress at pull-off. We highlight that the pull-off force has a simple solution of P_cr=−(k+3)πRΔγ/2 where Δγ is the work of adhesion and make further discussions with respect to three interesting limits: the classical JKR solution when k=0, the Gibson solid when k→1 and ν=0.5, and the strength limit in which the interfacial stress reaches the theoretical strength of adhesion at pull-off.     

Axisymmetric resonant disturbances in a homogeneous wave guide

The initial boundary-value linear stability problem for small localised axisymmetric disturbances in a homogeneous elastic wave guide, with the free upper surface and the lower surface being rigidly attached to a half-space, is formally solved by applying the Laplace transform in time and the Hankel transforms of zero and first orders in space. An asymptotic evaluation of the solution, expressed as a sum of inverse Laplace-Hankel integrals, is carried out by using the approach of the mathematical formalism of absolute and convective instabilities. It is shown that the dispersion-relation function of the problem D₀ (κ, ω), where the Hankel parameter κ is substituted by a wave number (and the Fourier parameter) κ, coincides with the dispersion-relation function D₀ (k, ω) for two-dimensional (2-D) disturbances in a homogeneous wave guide, where ω is the frequency (and the Laplace parameter) in both cases. An analysis for localised 2-D disturbances in a homogeneous wave guide is then applied. We obtain asymptotic expressions for wave packets, triggered by axisymmetric perturbations localised in space and finite in time, as well as for responses to axisymmetric sources localised in space, with the time dependence satisfying e^−iω₀t + O(e^−εt) for t → ∞, where Im ω₀ = 0, ε > 0, and t denotes time, i.e. for signalling with frequency ω₀. We demonstrate that, for certain combinations of physical parameters, axisymmetric wave packets with an algebraic temporal decay and axisymmetric signalling with an algebraic temporal growth, as √t, i.e., axisymmetric temporal resonances, are present in a neutrally stable homogeneous wave guide. The set of physically relevant wave guides having axisymmetric resonances is shown to be fairly wide. Furthermore, since an axisymmetric part of any source is L²-orthogonal to its non-axisymmetric part, a 3-D signalling with a non-vanishing axisymmetric component at an axisymmetric resonant frequency will generally grow algebraically in time. These results support our hypothesis concerning a possible resonant triggering mechanism of certain earthquakes, see Brevdo, 1998, J. Elasticity, 49, 201–237.     

Properties of Elastic Waves in a Non-Newtonian (Maxwell) Fluid-Saturated Porous Medium

The present study investigates novelties brought into the classic Biot's theory of propagation of elastic waves in a fluid-saturated porous solid by inclusion of non-Newtonian effects that are important, for example, for hydrocarbons. Based on our previous results (Tsiklauri and Beresnev, 2001), we investigated the propagation of rotational and dilatational elastic waves by calculating their phase velocities and attenuation coefficients as a function of frequency. We found that the replacement of an ordinary Newtonian fluid by a Maxwell fluid in the fluid-saturated porous solid results in: (a) an overall increase of the phase velocities of both the rotational and dilatational waves. With the increase of frequency these quantities tend to a fixed, higher level, as compared to the Newtonian limiting case, which does not change with the decrease of the Deborah number . (b) The overall decrease of the attenuation coefficients of both the rotational and dilatational waves. With the increase of frequency these quantities tend to a progressively lower level, as compared to the Newtonian limiting case, as decreases. (c) Appearance of oscillations in all physical quantities in the deeply non-Newtonian regime.     

Effect of rigid boundary on propagation of torsional surface waves in porous elastic layer

The paper presents the effect of a rigid boundary on the propagation of torsional surface waves in a porous elastic layer over a porous elastic half-space using the mechanics of the medium derived by Cowin and Nunziato (Cowin, S. C. and Nunziato, J. W. Linear elastic materials with voids. Journal of Elasticity, 13(2), 125–147 (1983)). The velocity equation is derived, and the results are discussed. It is observed that there may be two torsional surface wave fronts in the medium whereas three wave fronts of torsional surface waves in the absence of the rigid boundary plane given by Dey et al. (Dey, S., Gupta, S., Gupta, A. K., Kar, S. K., and De, P. K. Propagation of torsional surface waves in an elastic layer with void pores over an elastic half-space with void pores. Tamkang Journal of Science and Engineering, 6(4), 241–249 (2003)). The results also reveal that in the porous layer, the Love wave is also available along with the torsional surface waves. It is remarkable that the phase speed of the Love wave in a porous layer with a rigid surface is different from that in a porous layer with a free surface. The torsional waves are observed to be dispersive in nature, and the velocity decreases as the oscillation frequency increases.     

Two dimensional wave problems in rotating elastic media

The rotation of an elastic medium makes it act anisotropically and dispersively. The eigenvectors for plane wave propagation are in general complex and thus the waves are elliptically polarized. In general the waves are neither pure shear nor pure compressional waves, and their speeds depend on the ratio of rotational frequency of the medium and the angular frequency of the wave.The class of problems discussed here involves waves propagating perpendicularly to the axis of rotation and in particular we discuss plane strain modes. The reflection and refraction of plane waves is considered.The plane waves are used to construct a general solution in cylindrical coordinates. The solution is given in terms of Bessel functions. The cylindrical solution is applied to scattering by circular cylinders. The problem of free oscillations is mentioned briefly.     

Uniformly valid solution of limit cycle of the Duffing–van der Pol equation

The limit cycle of the Duffing–van der Pol equation is studied. By considering the product of the frequency ω of the limit cycle and the coefficient ε as an independent parameter μ=εω, an equivalent equation is obtained and then solved by Liao’s homotopy analysis method. The frequency ω is deduced as a function of μ and δ. This function provides us with an algebraic equation for ω, according to which we have an analytical approximation for the frequency. Numerical examples show that the attained approximation is very accurate. More importantly, the results are uniformly valid for all positive values of ε.     

Similarities Between Rotation and Stratification Effects on Homogeneous Shear Flow

Linear theory is applied to examine rotation and buoyancy effects on homogeneous turbulent shear flows with given vertical velocity shear, S=d/dx ₃. In the rotating shear case (where the rotation vector is perpendicular to the plane of the mean flow, Ω _i =Ωδ _{i 2}), general solutions for the Fourier components of the fluctuating velocity are proposed. These solutions are compared with those proposed in the literature for the Fourier components of the fluctuating velocity and density in the case of a homogeneous stratified shear flow with vertical density gradient, S _ρ=d/dx ₃. It is shown that from the normal mode stability stand point the Bradshaw parameter B=2Ω/S(1+2Ω/S) (in the rotating shear case) and the Richardson number R _i (in the statified shear case) play similar roles in identifying the stability for all the wave components except in the case where Ω·k=0, for which rotation has no effects on the flow. Analysis of the long-time behavior of the non-dimensional spectral density of energy, e ^g , is carried out. In the stable case, e ^g has decaying oscillations or undergoes a power law decay in time. Analytical solutions for the streamwise two-dimensional energy ℰ _ii ¹/2 (i.e. the limit at k ₁=0 of the one-dimensional energy spectra) are proposed. At large time, ℰ _ii ¹(t)/ℰ _ii ¹(0) oscillates around the value (3R _i +1)/(4R _i ) except at R _i =1 it stays constant in time. Similar behavior for ℰ _ii ¹(t)/ℰ _ii ¹(0) is also observed in the rotating shear case (ℰ _ii ¹(t)/ℰ _ii ¹(0) oscillates around the value (1+4B)/(4B)). Due to the behavior of the dimensionless spectral density of energy in both flow cases, the turbulent kinetic energy, /2, the production rate, ?, and the rate due to the buoyancy forces, ℬ, are split into two parts, , ?=?₁+?₂, ℬ=ℬ₁+ℬ₂ (in the stratified shear case, both ?₁ and ℬ₁ vanish when R _i >?, while in the rotating shear case one has ℬ=0). It is shown that when rotation is “cyclonic” (i.e. Ω/S>0), part reaches maximum magnitudes at St ≈2, independent of the B value, and the first time to a zero crossing of ?₂ occurs at this particular value. When rotation is “anticyclonic” (i.e. Ω/S<0) one finds St ≈1.6 instead of St ≈2. In the stratified shear case, both ?₂ and ℬ₂ cross zero at Nt=St ≈2, and part reaches maximum magnitudes at this particular value. These results and in particular those for the turbulent kinetic energy are compared with previous direct numerical simulation (DNS) results in homogeneous stratified shear flows. Received 30 July 2001 and accepted 19 February 2002     

PIV–LES analysis of channel flow rotating about the streamwise axis

The flow field of a channel rotating about the streamwise axis is analyzed experimentally and numerically. The current investigations were carried out at a bulk velocity based Reynolds number of Re_m = 2850 and a friction velocity based Reynolds number of Re_τ = 180, respectively. Particle-image velocimetry (PIV) measurements are compared with large-eddy simulation data to show earlier direct numerical simulation findings to generate too large a reverse flow region in the center region of the spanwise flow. The development of the mean spanwise velocity distribution and the influence of the rotation on the turbulent properties, i.e., the Reynolds stresses and the two-point correlations of the flow, are confirmed in both investigations. The rotation primarily influences those components of the Reynolds shear stresses, which contain the spanwise velocity component. The size of the correlation areas and thus the length scales of the flow generally grow in all three coordinate directions leading to longer structures. Furthermore, experimental results of the same channel flow at a significantly lower bulk Reynolds number of Re_{m, l} = 665, i.e., a laminar flow in a non-rotating channel, are introduced. The experiments show the low Reynolds number flow to become turbulent under rotation and to develop the same characteristics as the high Reynolds number flow.