Effect of open joint on stress wave propagation

Stress wave characteristics are drastically altered by joints and other inhomogenities. This paper addresses the effect of an open joint on stress wave transmission. An elastodynamic analysis is developed to supplement and explain some recent observations by Fourney and Dick(1995) on open as well as filled joints. The analytical model developed here assuming spherical symmetry can be extended to filled joints between dissimilar media, but results are presented only for open joints separating identical materials. As a special case, stress wave transmission across a joint with no gap is also addressed.

Investigation of Dynamic Fracture Toughness under Compressive-Shear Combinede Stress Wave Loading

提出亚微秒单脉冲应力波载荷作用下II型裂纹的平板冲击实验技术。加载率为dK/dt-10~8MPa·m~{”/d}·s~{-1}。实验中由锰铜应力片和弹性波理论分别测定和计算了压应力；通过微观分析确定了动态裂纹的平均扩展长度；引进等效应力强度因子,用动态断裂理论确定了60号钢的动态断裂韧性K_{Id}和K_{IId}；建立了亚微秒冲击载荷作用下确定材料动态断裂韧性的方法。

Stress Wave Propagation in a Gradient Elastic Medium

The gradient elastic constitutive equation incorporating the second gradient of the strains is used to determine the monochromatic elastic plane wave propagation in a gradient infinite medium and thin rod. The equation of motion, together with the internal material length, has been derived. Various dispersion relations have been determined. We present explicit expressions for the relationship between various wave speeds, wavenumber and internal material length.

Thermal stress wave and spallation induced by an electron beam

Thermal stress wave and spallation in aluminium alloy exposed to a high fluency and low energy electron beams are studied theoretically. A simple model for the study of energy deposition of electrons in materials is presented on the basis of some empirical formulae. Under the stress wave induced by energy deposition, microcracks and/or microvoids may appear in target materials, and in this case, the inelastic volume deformation should not vanish. The viscoplastic model proposed by Bodner and Partom with corresponding Gurson's yield function requires modification for this situation. The new constitutive model contains a scalar field variable description of the material damage which is taken as the void volume fraction of the polycrystalline material. Incorporation of the damage parameter permits description of rate-dependent, compressible, inelastic deformation and ductile fracture. The melting phenomenon has been observed in the experiment, therefore one needs to take into account the melting process in the intermediate energy deposition range. A three-phase equation of state used in the paper provides a more detailed and thermodynamical description of metals, particularly, in the melting region. The computational results based on the suggested model are compared with the experimental test for aluminium alloy, which is subjected to a pulsed electron beam with high fluency and low energy. (C) 1997 Elsevier Science Ltd.

An asymptotic analysis for one dimensional stress wave propagation in damaged viscoelastic media

A regular perturbation technique is suggested to deal with the problem of one dimensional stress wave propagation in viscoelastic media with damage. Based upon the first order asymptotic solution obtained, the characteristics of wave attenuation are studied. In fact, there exist three different time-dependent phenomena featuring the dynamic response of the materials, the first expressing the characteristics of wave propagation, the second indicating the innate effect of visco-elastic matrix and the third coming from the time dependent damage. The comparision of first order asymptotic solution with the numerical results calculated by a finite difference procedure shows that the perturbation expansion technique may offer a useful approach to the problem concerned.

Analysis and optimization of stress wave propagation in two-dimensional granular crystals with defects

Granular crystals are compact periodic assemblies of elastic particles in Hertzian contact whose dynamic response can be tuned from strongly nonlinear to linear by the addition of a static precompression force. This unique feature allows for a wide range of studies that include the investigation of new fundamental nonlinear phenomena in discrete systems such as solitary waves, shock waves, discrete breathers and other defect modes. In the absence of precompression, a particularly interesting property of these systems is their ability to support the formation and propagation of spatially localized soliton-like waves with highly tunable properties. The wealth of parameters one can modify (particle size, geometry and material properties, periodicity of the crystal, presence of a static force, type of excitation, etc.) makes them ideal candidates for the design of new materials for practical applications. This thesis describes several ways to optimally control and tailor the propagation of stress waves in granular crystals through the use of heterogeneities (interstitial defect particles and material heterogeneities) in otherwise perfectly ordered systems. We focus on uncompressed two-dimensional granular crystals with interstitial spherical intruders and composite hexagonal packings and study their dynamic response using a combination of experimental, numerical and analytical techniques. We first investigate the interaction of defect particles with a solitary wave and utilize this fundamental knowledge in the optimal design of novel composite wave guides, shock or vibration absorbers obtained using gradient-based optimization methods.

Study of Brazilian commercial Oriented Strand Board panels using stress wave

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Behaviour of stress wave propagation in utility timber pole

Non-destructive testing has been used for many years to evaluate the in situ condition of timber piles. Longitudinal impact is usually applied on the top of piles to induce longitudinal wave to detect faults in piles due to the fact that the longitudinalwave has less dispersive nature at lowfrequency. On the other hand,when it comes to evaluation of poles in situ, it is different as poles are partly embedded in soil and it is more practical to produce bending waves, as the top of the pole is not easily accessible. However, bending wave is known for its highly dispersive nature; especially in the low frequency range which is usually induced in low strain integrity testing. As bending wave can be considered as a hybrid of longitudinal and shear waves, it will be helpful, if it could detect the component of these twowaves separately.To do so, components of displacements or accelerations along radial and longitudinal directions need to be determined. By applying Fast Fourier Transform (FFT) on the signals, the dominant frequencies can be obtained. It has been found that, the longitudinal component decreases along radial direction which indicates the presence of bending wave component and this finding allows to the application of ContinuousWavelet Transform (CWT) on the longitudinal component of wave signals in order to obtain phase velocity. Phase velocities at different frequencies are then determined to draw the dispersive curve and compare with analytical phase velocity curve. The dispersion curve matched well with the analytical curve. © 2013 Taylor & Francis Group.

Determinations of stress wave velocity in a timber pole using wavelet transform

This paper presents an application of Wavelet Transfonn (WT) for determination of stress wave velocity for Non-destructive Testing of timber utility poles in service. For surface Non-destructive Testing (NDT), the hammer impact, which produces generally broadband frequency excitation, is used to generate stress wave. Moreover, due to practicality the impact location for field testing of a utility pole is on the side of the pole and 1.5 m above ground level. And the geometry of utility pole could not guarantee non-dispersive longitudinal wave. All of these issues have resulted in lack of accuracy and reliability of results from surface NDT in field testing. In recognition of such problem, this research explores methods to reliably calculate desired wave velocity by isolating wave mode and studying dispersive nature of utility pole. Fast Fourier Transfonn (FFT) is firstly conducted to determine the suitable frequency from a stress wave data. Then WT is applied on the wave data mentioned to perfonn time-frequency analysis. Velocity can be detennined by time history data of desired frequency from WT results which will be compared with the available analytical solution for longitudinal wave velocity. The results of the investigation showed that wavelet transfonn analysis can be a reliable signal processing tool for non-destructive testing in tenns of velocity detennination, which in tum also helps to detennine the embedded length of the timber pole.

Simulation of AE wave propagation in thin plate for source identification

In most materials, short stress waves are generated during the process of plastic deformation, phase transformation, crack formation and crack growth. These phenomena are applied in acoustic emission (AE) for the detection of material defects in wide spectrum areas, ranging from non-destructive testing for the detection of materials defects to monitoring of microeismical activity. AE technique is also used for defect source identification and for failure detection. AE waves consist of P waves (primary/longitudinal waves), S waves (shear/transverse waves) and Rayleight (surface) waves as well as reflected and diffracted waves. The propagation of AE waves in various modes has made the determination of source location difficult. In order to use the acoustic emission technique for accurate identification of source location, an understanding of wave propagation of the AE signals at various locations in a plate structure is essential. Furthermore, an understanding of wave propagation can also assist in sensor location for optimum detection of AE signals. In real life, as the AE signals radiate from the source it will result in stress waves. Unless the type of stress wave is known, it is very difficult to locate the source when using the classical propagation velocity equations. This paper describes the simulation of AE waves to identify the source location in steel plate as well as the wave modes. The finite element analysis (FEA) is used for the numerical simulation of wave propagation in thin plate. By knowing the type of wave generated, it is possible to apply the appropriate wave equations to determine the location of the source. For a single plate structure, the results show that the simulation algorithm is effective to simulate different stress waves.

I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

Controlling wave propagation through nonlinear engineered granular systems

We study the fundamental dynamic behavior of a special class of ordered granular systems in order to design new, structured materials with unique physical properties. The dynamic properties of granular systems are dictated by the nonlinear, Hertzian, potential in compression and zero tensile strength resulting from the discrete material structure. Engineering the underlying particle arrangement of granular systems allows for unique dynamic properties, not observed in natural, disordered granular media. While extensive studies on 1D granular crystals have suggested their usefulness for a variety of engineering applications, considerably less attention has been given to higher-dimensional systems. The extension of these studies in higher dimensions could enable the discovery of richer physical phenomena not possible in 1D, such as spatial redirection and anisotropic energy trapping. We present experiments, numerical simulation (based on a discrete particle model), and in some cases theoretical predictions for several engineered granular systems, studying the effects of particle arrangement on the highly nonlinear transient wave propagation to develop means for controlling the wave propagation pathways. The first component of this thesis studies the stress wave propagation resulting from a localized impulsive loading for three different 2D particle lattice structures: square, centered square, and hexagonal granular crystals. By varying the lattice structure, we observe a wide range of properties for the propagating stress waves: quasi-1D solitary wave propagation, fully 2D wave propagation with tunable wave front shapes, and 2D pulsed wave propagation. Additionally the effects of weak disorder, inevitably present in real granular systems, are investigated. The second half of this thesis studies the solitary wave propagation through 2D and 3D ordered networks of granular chains, reducing the effective density compared to granular crystals by selectively placing wave guiding chains to control the acoustic wave transmission. The rapid wave front amplitude decay exhibited by these granular networks makes them highly attractive for impact mitigation applications. The agreement between experiments, numerical simulations, and applicable theoretical predictions validates the wave guiding capabilities of these engineered granular crystals and networks and opens a wide range of possibilities for the realization of increasingly complex granular material design.

Spectral element modelling of wave propagation with boundary and structural discontinuity reflections

Spectral element method is very efficient in modelling high-frequency stress wave propagation because it works in the frequency domain. It does not need to use very fine meshes in order to capture high frequency wave energy as the time domain methods do, such as finite element method. However, the conventional spectral element method requires a throw-off element to be added to the structural boundaries to act as a conduit for energy to transmit out of the system. This makes the method difficult to model wave reflection at boundaries. To overcome this limitation, imaginary spectral elements are proposed in this study, which are combined with the real structural elements to model wave reflections at structural boundaries. The efficiency and accuracy of this proposed approach is verified by comparing the numerical simulation results with measured results of one dimensional stress wave propagation in a steel bar. The method is also applied to model wave propagation in a steel bar with not only boundary reflection, but also reflections from single and multiple cracks. The reflection and transmission coefficients, which are obtained from the discrete spring model, are adopted to quantify the discontinuities. Experimental tests of wave propagation in a steel bar with one crack of different depths are also carried out. Numerical simulations and experimental results show that the proposed method is effective and reliable in modelling wave propagation in one-dimensional waveguides with reflections from boundary and structural discontinuities. The proposed method can be applied to effectively model stress wave propagation for structural damage detection.

Wave-induced changes in seaweed toughness entail plastic modifications in snail traits maintaining consumption efficacy, link to supplementary material

1. Environmental stress can influence species traits and performance considerably. Using a seaweed-snail system from NW (Nova Scotia) and NE (Helgoland) Atlantic rocky shores, we examined how physical stress (wave exposure) modulates traits in the seaweed Fucus vesiculosus and indirectly in its main consumer, the periwinkle Littorina obtusata. 2. In both regions, algal tissue toughness increased with wave exposure. Reciprocal-transplant experiments showed that tissue toughness adjusts plastically to the prevailing level of wave exposure. 3. Choice experiments tested the feeding preference of snails from sheltered, exposed, and very exposed habitats for algae from such wave exposures. Snails from exposed and very exposed habitats consumed algal tissues at similar rates irrespective of the exposure of origin of the algae. However, snails from sheltered habitats consumed less algal tissues from very exposed habitats than tissues from sheltered and exposed habitats. Choice assays using reconstituted algal food (triturated during preparation) identified high thallus toughness as the explanation for the low preference of snails from sheltered habitats for algae from very exposed habitats. 4. Ultrastructural analyses of radulae indicated that rachidian teeth were longest and the number of cusps in lateral teeth (grazing-relevant traits) was highest in snails from very exposed habitats, suggesting that radulae are best suited to rupture tough algal tissues in such snails. 5. No-choice feeding experiments revealed that these radular traits are also phenotypically plastic, as they adjust to the toughness of the algal food. 6. Synthesis. This study indicates that the observed plasticity in the feeding ability of snails is mediated by wave exposure through phenotypic plasticity in the tissue toughness of algae. Thus, plasticity in consumers and their resource species may reduce the potential effects of physical stress on their interaction.