Wave propagation of a functionally graded plate via integral variables with a hyperbolic arcsine function

Several studies on functionally graded materials(FGMs)have been done by researchers,but few studies have dealt with the impact of the modification of the properties of materials with regard to the functional propagation of the waves in plates.This work aims to explore the effects of changing compositional characteristics and the volume fraction of the constituent of plate materials regarding the wave propagation response of thick plates of FGM.This model is based on a higher-order theory and a new displacement field with four unknowns that introduce indeterminate integral variables with a hyperbolic arcsine function.The FGM plate is assumed to consist of a mixture of metal and ceramic,and its properties change depending on the power functions of the thickness of the plate,such as linear,quadratic,cubic,and inverse quadratic.By utilizing Hamilton’s principle,general formulae of the wave propagation were obtained to establish wave modes and phase velocity curves of the wave propagation in a functionally graded plate,including the effects of changing compositional characteristics of materials.

Non-dimensional analysis on blast wave propagation in foam concrete:Minimum thickness to avoid stress enhancement

Foam concrete is a prospective material in defense engineering to protect structures due to its high energy absorption capability resulted from the long plateau stage.However,stress enhancement rather than stress mitigation may happen when foam concrete is used as sacrificial claddings placed in the path of an incoming blast load.To investigate this interesting phenomenon,a one-dimensional difference model for blast wave propagation in foam concrete is firstly proposed and numerically solved by improving the second-order Godunov method.The difference model and numerical algorithm are validated against experimental results including both the stress mitigation and the stress enhancement.The difference model is then used to numerically analyze the blast wave propagation and deformation of material in which the effects of blast loads,stress-strain relation and length of foam concrete are considered.In particular,the concept of minimum thickness of foam concrete to avoid stress enhancement is proposed.Finally,non-dimensional analysis on the minimum thickness is conducted and an empirical formula is proposed by curve-fitting the numerical data,which can provide a reference for the application of foam concrete in defense engineering.

Effect of porous surface layer on wave propagation in elastic cylinder immersed in fluid

To study the damage to an elastic cylinder immersed in fluid, a model of an elastic cylinder wrapped with a porous medium immersed in fluid is designed. This structure can both identify the properties of guided waves in a more practical model and address the relationship between the cylinder damage degree and the surface and surrounding medium. The principal motivation is to perform a detailed quantitative analysis of the longitudinal mode and flexural mode in an elastic cylinder wrapped with a porous medium immersed in fluid. The frequency equations for the propagation of waves are derived each for a pervious surface and an impervious surface by employing Biot theory. The influences of the various parameters of the porous medium wrapping layer on the phase velocity and attenuation are discussed. The results show that the influences of porosity on the dispersion curves of guided waves are much more significant than those of thickness,whereas the phase velocity is independent of the static permeability. There is an apparent “mode switching” between the two low-order modes. The characteristics of attenuation are in good agreement with the results from the dispersion curves.This work can support future studies for optimizing the theory on detecting the damage to cylinder or pipeline.

Wave propagation responses of porous bi-directional functionally graded magneto-electro-elastic nanoshells via nonlocal strain gradient theory

This study examines the wave propagation characteristics for a bi-directional functional grading of barium titanate(BaTiO_(3)) and cobalt ferrite(CoFe_(2)O_(4)) porous nanoshells,the porosity distribution of which is simulated by the honeycomb-shaped symmetrical and asymmetrical distribution functions.The nonlocal strain gradient theory(NSGT) and first-order shear deformation theory are used to determine the size effect and shear deformation,respectively.Nonlocal governing equations are derived for the nanoshells by Hamilton's principle.The resulting dimensionless differential equations are solved by means of an analytical solution of the combined exponential function after dimensionless treatment.Finally,extensive parametric surveys are conducted to investigate the influence of diverse parameters,such as dimensionless scale parameters,radiusto-thickness ratios,bi-directional functionally graded(FG) indices,porosity coefficients,and dimensionless electromagnetic potentials on the wave propagation characteristics.Based on the analysis results,the effect of the dimensionless scale parameters on the dispersion relationship is found to be related to the ratio of the scale parameters.The wave propagation characteristics of nanoshells in the presence of a magnetoelectric field depend on the bi-directional FG indices.

Mesh Size Effect in Numerical Simulation of Blast Wave Propagation and Interaction with Structures

Numerical method is popular in analysing the blast wave propagation and interaction with structures.However,because of the extremely short duration of blast wave and energy trans-mission between different grids,the numerical results are sensitive to the finite element mesh size.Previous numerical simulations show that a mesh size acceptable to one blast scenario might not be proper for another case,even though the difference between the two scenarios is very small,indicating a simple numerical mesh size convergence test might not be enough to guarantee accu-rate numerical results.Therefore,both coarse mesh and fine mesh were used in different blast scenarios to investigate the mesh size effect on numerical results of blast wave propagation and interaction with structures.Based on the numerical results and their comparison with field test re-sults and the design charts in TM5-1300,a numerical modification method was proposed to correct the influence of the mesh size on the simulated results.It can be easily used to improve the accu-racy of the numerical results of blast wave propagation and blast loads on structures.

Dynamic mechanical properties and wave propagation of composite rock-mortar specimens based on SHPB tests

Filled inclusions in rock discontinuities play a key role in the mechanical characteristics of the rock and thereby influence the stability of rock engineering. In this study, a series of impact tests were performed using a split Hopkinson pressure bar system with high-speed photography to investigate the effect of interlayer strength on the wave propagation and fracturing process in composite rock-mortar specimens.The results indicate that the transmission coefficient, nominal dynamic strength, interlayer closure, and specific normal stiffness generally increase linearly with increasing interlayer stiffness. The cement mortar layer can serve as a buffer during the deformation of composite specimens. The digital images show that tensile cracks are typically initiated at the rock-mortar interface, propagate along the loading direction, and eventually result in a tensile failure regardless of the interlayer properties. However, when a relatively weaker layer is sandwiched between the rock matrix, an increasing amount of cement mortar is violently ejected and slight slabbing occurs near the rock-mortar interface.

Experimental verification of effect of horizontal inhomogeneity of evaporation duct on electromagnetic wave propagation

The evaporation duct which forms above the ocean surface has a significant influence on electromagnetic wave propagation above 2 GHz over the ocean. The effects of horizontal inhomogeneity of evaporation duct on electromagnetic wave propagation are investigated, both in numerical simulation and experimental observation methods, in this paper. Firstly, the features of the horizontal inhomogeneity of the evaporation duct are discussed. Then, two typical inhomogeneous cases are simulated and compared with the homogeneous case. The result shows that path loss is significantly higher than that in the homogeneous case when the evaporation duct height （EDH） at the receiver is lower than that at the transmitter. It is also concluded that the horizontal inhomogeneity of the evaporation duct has a significant influence when the EDH is low or when the electromagnetic wave frequency is lower than 13 GHz. Finally, experimental data collected on a 149-km long propagation path in the South China Sea in 2013 are used to verify the conclusion. The experimental results are consis- tent with the simulation results. The horizontal inhomogeneity of evaporation duct should be considered when modeling electromagnetic wave propagation over the ocean.

Wave propagation control in periodic track structure through local resonance mechanism

Excessive vibration and noise radiation of the track structure can be caused by the operation of high speed trains.Though the track structure is characterized by obvious periodic properties and band gaps,the bandwidth is narrow and the elastic wave attenuation capability within the band gap is weak.In order to effectively control the vibration and noise of track structure,the local resonance mechanism is introduced to broaden the band gap and realize wave propagation control.The locally resonant units are attached periodically on the rail,forming a new locally resonant phononic crystal structure.Then the tuning of the elastic wave band gaps of track structure is discussed,and the formation mechanism of the band gap is explicated.The research results show that a new wide and adjustable locally resonant band gap is formed after the resonant units are introduced.The phenomenon of coupling and transition can be observed between the new locally resonant band gap and the original band gap of the periodic track structure with the band gap width reaching the maximum at the coupling position.The broader band gap can be applied for vibration and noise reduction in high speed railway track structure.

Effects of axial static stress on stress wave propagation in rock considering porosity compaction and damage evolution

A wave equation of rock under axial static stress is established using the equivalent medium method by modifying the Kelvin-Voigt model.The analytical formulas of longitudinal velocity,space and time attenuation coefficients and response frequency are obtained by solving the equation using the harmonic method.A series of experiments on stress wave propagation through rock under different axial static stresses have been conducted.The proposed models of stress wave propagation are then verified by comparing experimental results with theoretical solutions.Based on the verified theoretical models,the influences of axial static stress on longitudinal velocity,space and time attenuation coefficients and response frequency are investigated by detailed parametric studies.The results show that the proposed theoretical models can be used to effectively investigate the effects of axial static stress on the stress wave propagation in rock.The axial static stress influences stress wave propagation characteristics of porous rock by varying the level of rock porosity and damage.Moreover,the initial porosity,initial elastic modulus of the rock voids and skeleton,viscous coefficient and vibration frequency have significant effects on the P-wave velocity,attenuation characteristics and response frequency of the stress wave in porous rock under axial static stress.

An extended displacement discontinuity method for analysis of stress wave propagation in viscoelastic rock mass

An extended displacement discontinuity method (EDDM) is proposed to analyze the stress wave propagation in jointed viscoelastic rock mass (VRM).The discontinuities in a rock mass are divided into two groups.The primary group with an average geometrical size larger than or in the same order of magnitude of wavelength of a concerned stress wave is defined as 'macro-joints',while the secondary group with a high density and relatively small geometrical size compared to the wavelength is known as 'micro-defects'.The rock mass with micro-defects is modeled as an equivalent viscoelastic medium while the macro-joints in the rock mass are modeled explicitly as physical discontinuities.Viscoelastic properties of a micro-defected sedimentary rock are obtained by longitudinally impacting a cored long sedimentary rod with a pendulum.Wave propagation coefficient and dynamic viscoelastic modulus are measured.The EDDM is then successfully employed to analyze the wave propagation across macro-joint in VRM.The effect of the rock viscosity on the stress wave propagation is evaluated by comparing the results of VRM from the presented EDDM with those of an elastic rock mass (ERM) from the conventional displacement discontinuity method (CDDM).The CDDM is a special case of the EDDM under the condition that the rock viscosity is ignored.Comparison of the reflected and transmitted waves shows that the essential rock viscosity has a significant effect on stress wave attenuation.When a short propagation distance of a stress wave is considered,the results obtained from the CDDM approximate to the EDDM solutions,however,when the propagation distance is sufficiently long relative to the wavelength,the effect of rock viscosity on the stress wave propagation cannot be ignored.

Surface effect and non-local elasticity in wave propagation of functionally graded piezoelectric nano-rod excited to applied voltage

A non-local solution for a functionally graded piezoelectric nano-rod is pre- sented by accounting the surface effect. This solution is used to evaluate the charac- teristics of the wave propagation in the rod structure. The model is loaded under a two-dimensional （2D） electric potential and an initially applied voltage at the top of the rod. The mechanical and electrical properties are assumed to be variable along the thick- ness direction of the rod according to the power law. The Hamilton principle is used to derive the governing differential equations of the electromechanical system. The effects of some important parameters such as the applied voltage and gradation of the material properties on the wave characteristics of the rod are studied.

Nonlinear Effect of Wave Propagation in Shallow Water

In this paper, a nonlinear model is presented to describe wave transformation in shallow water with the zero- vorticity equation of wave- number vector and energy conservation equation. The nonlinear effect due to an empirical dispersion relation (by Hedges) is compared with that of Dalrymple's dispersion relation. The model is tested against the laboratory measurements for the case of a submerged elliptical shoal on a slope beach, where both refraction and diffraction are significant. The computation results, compared with those obtained through linear dispersion relation, show that the nonlinear effect of wave transformation in shallow water is important. And the empirical dispersion relation is suitable for researching the nonlinearity of wave in shallow water.

Applications of A Numerical Model to Wave Propagation on Mild Slopes

Based on the mild slope equation that has heen deeomposed inlo three equations related to wave phase function, wave amplitude and wave approach angle, a refraction-diffraction model is developed. The finite difference method has been selected as the solution method. The model results are compared with experimental results and the model is applied to coastal waters of the Fethiye Bay, whieh is located at the Mediterranean Sea of Turkey.

Vibration and wave propagation analysis of twisted micro-beam using strain gradient theory

In this research, vibration and wave propagation analysis of a twisted micro- beam on Pasternak foundation is investigated. The strain-displacement relations （kine-matic equations） are calculated by the displacement fields of the twisted micro-beam. The strain gradient theory （SGT） is used to implement the size dependent effect at micro-scale. Finally, using an energy method and Hamilton＇s principle, the governing equations of motion for the twisted micro-beam are derived. Natural frequencies and the wave prop- agation speed of the twisted micro-beam are calculated with an analytical method. Also, the natural frequency, the phase speed, the cut-off frequency, and the wave number of the twisted micro-beam are obtained by considering three material length scale parameters, the rate of twist angle, the thickness, the length of twisted micro-beam, and the elastic medium. The results of this work indicate that the phase speed in a twisted micro-beam increases with an increase in the rate of twist angle. Moreover, the wave number is in- versely related with the thickness of micro-beam. Meanwhile, it is directly related to the wave propagation frequency. Increasing the rate of twist angle causes the increase in the natural frequency especially with higher thickness. The effect of the twist angle rate on the group velocity is observed at a lower wave propagation frequency.

Numerical Perspective of Second-Harmonic Generation of Circumferential Guided Wave Propagation in a Circular Tube

The effect of second-harmonic generation （SHG） by primary （fundamental） circumferential guided wave （CGW） propagation is investigated from a numerical standpoint. To enable that the second harmonic of the primary CGW mode can accumulate along the circumferential direction, an appropriate mode pair of primary and double frequency CGWs is chosen. Finite element simulations and evaluations of nonlinear CGW propagation are analyzed for the selected CGW mode pair. The numerical simulations performed directly demonstrate that the response of SHG is completely generated by the desired primary CGW mode that satisfies the condition of phase velocity matching at a specific driving frequency, and that the second harmonic of the primary CGW mode does have a cumulative effect with circumferential angles. The numerical perspective obtained yields an insight into the complicated physical process of SHG of primary CGW propagation unavailable previously.

A Study of the Stability Properties in Simulation of Wave Propagation with SPH Method

When ordinary Smoothed Particle Hydrodynamics （SPH） method is used to simulate wave propagation in a wave tank, it is usually observed that the wave height decays and the wave length elongates along the direction of wave propagation. Accompanied with this phenomenon, the pressure under water decays either and shows a big oscillation simultaneously. The reason is the natural potential tensile instability of modeling water motion with ordinary SPH which is caused by particle negative stress in the computation. I＇o deal with the problems, a new sextic kernel function is proposed to reduce this instability. An appropriate smooth length is given and its computation criterion is also suggested. At the same time, a new kind dynamic boundary condition is introduced. Based on these improvements, the new SPH method named stability improved SPH （SISPH） can simulate the wave propagation well. Both the water surface and pressure can be well expressed and the oscillation of pressure is nearly eliminated. Compared with other improved methods, SISPH can truly reveal the physical reality without bringing some new problems in a simple way.

Accelerating the discontinuous Galerkin method for seismic wave propagation simulations using multiple GPUs with CUDA and MPI

We have successfully ported an arbitrary highorder discontinuous Galerkin method for solving the threedimensional isotropic elastic wave equation on unstructured tetrahedral meshes to multiple Graphic Processing Units （GPUs） using the Compute Unified Device Architecture （CUDA） of NVIDIA and Message Passing Interface （MPI） and obtained a speedup factor of about 28.3 for the single-precision version of our codes and a speedup factor of about 14.9 for the double-precision version. The GPU used in the comparisons is NVIDIA Tesla C2070 Fermi, and the CPU used is Intel Xeon W5660. To effectively overlap inter-process communication with computation, we separate the elements on each subdomain into inner and outer elements and complete the computation on outer elements and fill the MPI buffer first. While the MPI messages travel across the network, the GPU performs computation on inner elements, and all other calculations that do not use information of outer elements from neighboring subdomains. A significant portion of the speedup also comes from a customized matrix-matrix multiplication kernel, which is used extensively throughout our program. Preliminary performance analysis on our parallel GPU codes shows favorable strong and weak scalabilities.

EXACT ANALYSIS OF WAVE PROPAGATION IN AN INFINITE RECTANGULAR BEAM

The Fourier series method was extended for the exact analysis of wave propagation in an infinite rectangular beam.Initially,by solving the three-dimensional elastodynamic equations a general analytic solution was derived for wave motion within the beam.And then for the beam with stress-free boundaries,the propagation characteristics of elastic waves were presented.This accurate wave propagation model lays a solid foundation of simultaneous control of coupled waves in the beam.

PARTITION OF UNITY FINITE ELEMENT METHOD FOR SHORT WAVE PROPAGATION IN SOLIDS

A partition of unity finite element method for numerical simulation of short wave propagation in solids is presented. The finite element spaces were constructed by multiplying the standard isoparametric finite element shape functions, which form a partition of unity, with the local subspaces defined on the corresponding shape functions, which include a priori knowledge about the wave motion equation in trial spaces and approximately reproduce the highly oscillatory properties within a single element. Numerical examples demonstrate the performance of the proposed partition of unity finite element in both computational accuracy and efficiency.

A Refined Study on ECR Wave Propagation and Absorption in the Weakly Relativistic Plasma

A refined derivation of refraction and absorption of the pure O-mode and X-mode Electron Cyclotron Resonance (ECR) wave in tokamak plasma is carried out. The weakly- relativistic dielectric tensor elements are used and the results show that the refraction only changes a little, compared to that deduced from the cold-plasma dispersion relation even in the inner re- gion. Refined formulae of the wave damping rate are then obtained for both the O-mode and the X-mode fundamental waves.