In-plane free vibration of a single-crystal silicon ring

In this paper the natural frequencies and the associated mode shapes of in-plane free vibration of a single-crystal silicon ring are analyzed. It is found that the Si(1 1 1) ring is two-dimensionally isotropic in the (1 1 1) plane for elastic constants but three-dimensionally anisotropic, while the Si(1 0 0) ring is fully anisotropic. Hamilton’s principle is used to derive the equations of vibration, which is a set of partial differential equations with coefficients being periodic in polar variable. Expressing the radial and tangential displacements in sinusoidal form with non-predetermined amplitudes, and through the integration with respect to the circumferential variable, the original governing equations in partial differential form can be converted into the amplitude equations in ordinary differential form. The exact expressions for frequencies and mode shapes are obtained. It is found that for Si(1 0 0) rings the frequencies of a pair of modes, which are equal for an isotropic ring, split due to the anisotropic effect only for the second in-plane vibration mode. The phenomena of frequency splitting and degenerate modes can be proved either based on the conservation of averaged mechanical energy or by the concept of crystallographic symmetry groups. When the single-crystal silicon is replaced by the polycrystalline silicon, which is isotropic in elastic constants, the derived equations for frequencies correctly predict the vanishing of the phenomenon of frequency splitting.     

Hyperbolicity of Velocity-Stress Equations for Waves in Anisotropic Elastic Solids

This paper reports mathematical properties of the three-dimensional, first-order, velocity-stress equations for propagating waves in anisotropic, linear elastic solids. The velocity-stress equations are useful for numerical solution. The original equations include the equation of motion and the elasticity relation differentiated by time. The result is a set of nine, first-order partial differential equations (PDEs) of which the velocity and stress components are the unknowns. Cast into a vector-matrix form, the equations can be characterized by three Jacobian matrices. Hyperbolicity of the equations is formally proved by analyzing (i) the spectrum of a linear combination of the three Jacobian matrices, and (ii) the eigenvector matrix for diagonalizing the linearly combined Jacobian matrices. In the three-dimensional space, linearly combined Jacobian matrices are shown to be connected to the classic Christoffel matrix, leading to a simpler derivation for the eigenvalues and eigenvectors. The results in the present paper provide critical information for applying modern numerical methods, originally developed for solving conservation laws, to elastodynamics.     

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.     

Propagation and localization of two-dimensional in-plane elastic waves in randomly disordered layered piezoelectric phononic crystals

Two-dimensional in-plane wave propagation and localization in the disordered layered piezoelectric phononic crystals with material 6 mm are investigated taking the electromechanical coupling into account. The electric field is approximated as quasi-static. The analytical solutions of elastic waves are obtained. The 6 × 6 transfer matrix between two consecutive unit cells is obtained by means of the mechanical and electrical continuity conditions. The expressions of the localization factor and localization length in the disordered periodic structures are presented by regarding the variables of the mechanical and electrical fields as the elements of the state vector. The numerical results of the localization factors and localization lengths are presented for two kinds of disordered piezoelectric phononic crystals, i.e. ZnO–PZT–5H and PVDF–PZT–5H piezocomposites. It is seen from the results that the incident angle of elastic waves and the thickness of the piezoelectric ceramics have significant effects on the wave localization characteristics. For different piezoelectric phononic crystals, the effects of the incident angle are very different. Moreover, with the increase of the disorder degree, the localization phenomenon is strengthened.     

Analysis of non-axisymmetric wave propagation in a homogeneous piezoelectric solid circular cylinder of transversely isotropic material

A study concerning the propagation of free non-axisymmetric waves in a homogeneous piezoelectric cylinder of transversely isotropic material with axial polarization is carried out on the basis of the linear theory of elasticity and linear electro-mechanical coupling. The solution of the three dimensional equations of motion and quasi-electrostatic equation is given in terms of seven mechanical and three electric potentials. The characteristic equations are obtained by the application of the mechanical and two types of electric boundary conditions at the surface of the piezoelectric cylinder. A novel method of displaying dispersion curves is described in the paper and the resulting dispersion curves are presented for propagating and evanescent waves for PZT-4 and PZT-7A piezoelectric ceramics for circumferential wave numbers m = 1, 2, and 3. It is observed that the dispersion curves are sensitive to the type of the imposed boundary conditions as well as to the measure of the electro-mechanical coupling of the material.     

Surface electro-elastic Love waves in a layered structure with a piezoelectric substrate and a dielectric layer

The existence and behaviour of electro-elastic surface Love waves in a structure consisting of a piezoelectric substrate of crystal classe 6, 4, 6 mm, 4 mm, 622 or 422, an elastic layer and a dielectric medium is considered. The mathematical model obtained includes all the above crystal classes, i.e. the surface wave problems related to all these classes are presented in a single mathematical model. The dispersion equation for the existence of Love surface waves with respect to phase velocity is obtained. A detailed investigation of the electromechanical coupling coefficient is carried out depending on the dielectric and piezoelectric parameters of the problem. Geometrical investigation of the solutions of the dispersion equation is carried out.     

Surface electro-elastic Love waves in a layered structure with a piezoelectric substrate and two isotropic layers

Propagation of electro-elastic surface Love waves in a structure consisting of a piezoelectric half-space substrate of crystal class 6, 4, 6 mm or 4 mm and two layers, one of which (adjacent to the substrate) is a conducting material and the second is either a conducting or a dielectric material, is considered. The mathematical model obtained includes all the above crystal classes i.e. the surface wave problems related to all these classes are presented in a single mathematical model. The dispersion equation for the existence of Love surface waves with respect to phase velocity is obtained. Numerical calculations are carried out for three different layered structures. The effect of the second layer on the propagation behaviour of the surface Love wave in the structure is revealed.     

The dispersion of shear wave in multilayered magnetoelastic self-reinforced media

In this paper, we study the propagation of shear waves in a magnetoelastic self-reinforced medium using finite difference technique. Dispersion equation has been deduced for the case when (n ? 1) layers lie over a half space. It is observed that the obtained dispersion equation is in assertion with the classical Love wave equation for both the cases when a single and double layer lies over a half space. The stability condition for the used finite difference scheme and the expression for the phase and group velocity have been derived. The dispersion curve for different values of magnetoelastic coupling parameter, phase and group velocity variation for different values of stability ratio has been depicted by means of graphs.     

Failure in notched tension bars due to high-temperature creep: Interaction between nucleation controlled cavity growth and continuum cavity growth

For axi-symmetrically notched tension bars [Dyson, B.F., Loveday, M.S., 1981, Creep Fracture in Nimonic 80A under Tri-axial Tensile Stressing, In: Ponter A.R.S., Hayhurst, D.R. (Eds.), Creep in Structures, Springer-Verlag, Berlin, pp. 406–421] show two types of damage propagation are shown: for low stress, failure propagates from the outside notch surface to the centre-line; and for high stress, failure propagates from the centre-line to the outside notch surface. The objectives of the paper are to: identify the physics of the processes controlling global failure modes; and, describe the global behaviour using physics-based constitutive equations.Two sets of constitutive equations are used to model the softening which takes place in tertiary creep of Nimonic 80A at 750 °C. Softening by multiplication of mobile dislocations is firstly combined, for low stress, with softening due to nucleation controlled creep constrained cavity growth; and secondly combined, for high stress, with softening due to continuum void growth. The Continuum Damage Mechanics, CDM, Finite Element Solver DAMAGE XX has been used to study notch creep fracture. Low stress notch behaviour is accurately predicted provided that the constitutive equations take account of the effect of stress level on creep ductility. High stress notch behaviour is accurately predicted from a normalized inverse cavity spacing d/2? = 6, and an initial normalized cavity radius r_hi/? = 3.16 × 10^?3, where 2? is the cavity spacing, and d is the grain size; however, the constants in the strain rate equation required recalibration against high stress notch data. A void nucleation mechanism is postulated for high stress behaviour which involves decohesion where slip bands intersect second phase grain boundary particles. Both equation sets accurately predict experimentally observed global failure modes.     

Simulation of a valveless pump with an elastic tube

A valveless pump consisting of a pumping chamber with an elastic tube was simulated using an immersed boundary (IB) method. The interaction between the motion of the elastic tube and the pumping chamber generated a net flow toward the outlet throughout a full cycle of the pump. The net flow rate of the valveless pump was examined by varying the stretching coefficient (ϕ), bending coefficient (γ), the aspect ratio (l/d) of the elastic tube, and the frequency (f) of the pumping chamber. As the stretching and bending coefficients of the elastic tube increased, the net flow through the valveless pump decreased. Elastic tubes with aspect ratios in the range of 2 ⩽ l/d ⩽ 3 generated a higher flow rate than that generated for tubes with aspect rations of l/d = 1 or 4. As the frequency of the pumping chamber increased, the net flow rate of the pump for l/d = 2 increased. However, the net flow rate for l/d = 3 was nonlinearly related to the pumping frequency due to the complexity of the wave motions. Snapshots of the fluid velocity vectors and the wave motions of the elastic tube were examined over one cycle of the pump to gain a better understanding of the mechanism underlying the valveless pump. The relationship between the average gap in the elastic tube and the average flow rate of the pump was analyzed. A smaller gap in the elastic tube during the expansion mode and a wider gap in the elastic tube during the contraction mode played a dominant role in generating a high average flow rate in the pump, regardless of the stretching coefficient (ϕ), the aspect ratio (l/d) of the elastic tube, or the pumping frequency of the pumping chamber (f).     

Flow of mono-dispersed particles through horizontal bend

Pipeline slurry flow of mono-dispersed particles through horizontal bend is numerically simulated by implementing Eulerian two-phase model in FLUENT software. A hexagonal shape and Cooper type non-uniform three-dimensional grid is chosen to discretize the entire computational domain, and a control volume finite difference method is used to solve the governing equations. The modeling results are compared with the experimental data collected in 53.0 mm diameter horizontal bend with radius of 148.4 mm for concentration profiles and pressure drops. Experiments are performed on narrow-sized silica sand with mean diameter of 450 μm and for flow velocity up to 3.56 m/s (namely, 1.78, 2.67 and 3.56 m/s) and four efflux concentrations up to 16.28% (namely, 0%, 3.94%, 8.82% and 16.28%) by volume for each velocity. Eulerian model gives fairly accurate predictions for both the pressure drop and concentration profiles at all efflux concentrations and flow velocities.     

Surface responses induced by point load or uniform traction moving steadily on an anisotropic half-plane

Surface responses induced by point load or uniform traction moving steadily with subsonic speed on an anisotropic half-plane boundary are investigated. It is found that the effects of the material constant on surface displacements are through matrices L^?1(v) and S(v)L^?1(v), while those on surface stress components are through matrices Ω(v) and Γ(v). Explicit expressions for the elements of these four matrices are expressed in terms of elastic stiffness for general anisotropic materials. The special cases of monoclinic materials with symmetry plane at x₁ = 0, x₂ = 0 and x₃ = 0, and the case for orthotropic materials are all deduced. Results for isotropic material may be recovered from present results. For monoclinic materials with a plane of symmetry at x₃ = 0, two of the elements of matrix Ω(v) are found to be independent of subsonic speed.     

Cracked infinite cylinder with two rigid inclusions under axisymmetric tension

This paper considers the problem of an axisymmetric infinite cylinder with a ring shaped crack at z = 0 and two ring-shaped rigid inclusions with negligible thickness at z = ±L. The cylinder is under the action of uniformly distributed axial tension applied at infinity and its lateral surface is free of traction. It is assumed that the material of the cylinder is linearly elastic and isotropic. Crack surfaces are free and the constant displacements are continuous along the rigid inclusions while the stresses have jumps. Formulation of the mixed boundary value problem under consideration is reduced to three singular integral equations in terms of the derivative of the crack surface displacement and the stress jumps on the rigid inclusions. These equations, together with the single-valuedness condition for the displacements around the crack and the equilibrium equations along the inclusions, are converted to a system of linear algebraic equations, which is solved numerically. Stress intensity factors are calculated and presented in graphical form.     

Numerical simulation of the fountain flow instability in injection molding

For the first time, the viscoelastic flow front instability is studied in the full non-linear regime by numerical simulation. A two-component viscoelastic numerical model is developed which can predict fountain flow behavior in a two-dimensional cavity. The eXtended Pom-Pom (XPP) viscoelastic model is used. The levelset method is used for modeling the two-component flow of polymer and gas. The difficulties arising from the three-phase contact point modeling are addressed, and solved by treating the wall as an interface and the gas as a compressible fluid with a low viscosity. The resulting set of equations is solved in a decoupled way using a finite element formulation. Since the model for the polymer does not contain a solvent viscosity, the time discretized evolution equation for the conformation tensor is substituted into the momentum balance in order to obtain a Stokes like equation for computing the velocity and pressure at the new time level. Weissenberg numbers range from 0.1 to 10. The simulations reveal a symmetric fountain flow for Wi = 0.1–5. For Wi = 10 however, an oscillating motion of the fountain flow is found with a spatial period of three times the channel height, which corresponds to experimental observations.     

Laminar flow of power-law fluids past a rotating cylinder

In this work, the continuity and momentum equations have been solved numerically to investigate the flow of power-law fluids over a rotating cylinder. In particular, consideration has been given to the prediction of drag and lift coefficients as functions of the pertinent governing dimensionless parameters, namely, power-law index (1 ≥ n ≥ 0.2), dimensionless rotational velocity (0 ≤ α ≤ 6) and the Reynolds number (0.1 ≤ Re ≤ 40). Over the range of Reynolds number, the flow is known to be steady. Detailed streamline and vorticity contours adjacent to the rotating cylinder and surface pressure profiles provide further insights into the nature of flow. Finally, the paper is concluded by comparing the present numerical results with the scant experimental data on velocity profiles in the vicinity of a rotating cylinder available in the literature. The correspondence is seen to be excellent for Newtonian and inelastic fluids.     

Numerical analysis of turbulent convection heat transfer from an array of perforated fins

Numerical investigation is made for three-dimensional fluid flow and convective heat transfer from an array of solid and perforated fins that are mounted on a flat plate. Incompressible air as working fluid is modeled using Navier–Stokes equations and RNG based k ? ? turbulent model is used to predict turbulent flow parameters. Temperature field inside the fins is obtained by solving Fourier’s conduction equation. The conjugate differential equations for both solid and gas phase are solved simultaneously by finite volume procedure using SIMPLE algorithm. Perforations such as small channels with square cross section are arranged streamwise along the fin’s length and their numbers varied from 1 to 3. Flow and heat transfer characteristics are presented for Reynolds numbers from 2 × 10⁴ to 4 × 10⁴ based on the fin length and Prandtl number is taken Pr = 0.71. Numerical computations are validated with experimental studies of the previous investigators and good agreements were observed. Results show that fins with longitudinal pores, have remarkable heat transfer enhancement in addition to the considerable reduction in weight by comparison with solid fins.     

Adiabatic shear banding in plane strain tensile deformations of 11 thermoelastoviscoplastic materials with finite thermal wave speed

We study the initiation and propagation of adiabatic shear bands (ASBs) in 11 homogeneous materials each modeled as microporous, isotropic and thermoelastoviscoplastic, and deformed in plane strain tension. The heat conduction in each material is assumed to be governed by a hyperbolic heat equation; thus thermal and mechanical waves propagate with finite speeds. The decrease in the thermophysical parameters due to the increase in porosity is considered. An ASB is assumed to initiate at a material point when the maximum shear stress there has dropped to 80% of its peak value for that material point and it is deforming plastically. An approximate solution of the coupled nonlinear partial differential equations subject to suitable initial and boundary conditions is found by the finite element method (FEM). In contrast to the Considerè and the Hart criterion, it is found that an ASB initiates when the axial load drops rapidly and not when it peaks. The refinement of the 40 × 40 uniform FE mesh to 120 × 120 uniform elements decreased the ASB initiation time by 2.1% while increasing the CPU time by a factor of ∼26. By locating points where the ASB has initiated we find its current length, width and speed. The 11 materials are ranked according to the time of initiation of an ASB under otherwise identical geometric and loading conditions with the same initial nonuniform porosity distribution. This ranking of materials is found to differ somewhat from that ascertained by Batra and Kim (1992) who studied simple shearing deformations, and by Batra et al. (1995) who analyzed three-dimensional torsional deformations of thin-walled tubular specimens. The average axial strain determined from the maximum axial load condition differs noticeably from that when an ASB initiates.     

Wave propagation in a piezoelectric rod of 6 mm symmetry

The wave propagation in a piezoelectric rod of 6 mm symmetry is investigated by applying a 3-D piezoelectric elastic model. A self-adjoint method is introduced to solve this problem, this method avoids calculating the generalized eigenvalue equation, it completely draws the dispersion curves in the forms of Quasi-P wave, Quasi-SV wave and Quasi-SH wave under the self-adjoint boundary condition, and it can evaluate the dispersion curves of all kinds of boundary conditions. As an example, the dispersion curves of PLT-5H are completely drawn, we also found the Quasi-SV wave has standing wave phenomenon in the PLT-5H rod. In addition the relation of dispersion curves among different boundary conditions is discussed, and an experiment method is introduce to decide the dispersion curves for another boundary conditions.     

Simulation of initial formation and growth of martensitic band in NiTi micro-tube under tension

Previous experiments have shown that the distinct features of macro-martensitic band nucleation and propagation in micro-tube under tension are in three stages: the initiation and propagation of a single helical band → self-merging → propagation of the cylindrical band. In this paper, the martensitic formation and helical band propagation in the tube at different temperatures are modeled. The free energy function of the tube is formulated by introducing an equivalent method to calculate the stress and strain disturbances in the helical martensitic domain, and the phase transformation criterion is derived based on thermodynamics. The simulations successfully capture the main features of nucleation, pattern evolution and variation of front velocity of the helical martensitic band in the tube. The analytical results and the comparison with experiments are also discussed in this paper.     

Experimental and computational studies of shock wave-to-bubbly water momentum transfer

Momentum transfer from shock waves (SWs) of various intensity (from 0.05 MPa to 0.5 MPa in amplitude) to water containing air bubbles 2.5 to 4 mm of mean diameter is studied both experimentally and by means of numerical simulation. Experiments are performed in a vertical shock tube of a 50 × 100 mm² rectangular cross section consisting of a 495-mm long high-pressure section (HPS), 495-mm long low-pressure section (LPS), and 990 mm long test section (TS) equipped with an air bubbler and filled with water. Experiments have shown that as the initial gas volume fraction in water increases from 0 to 0.3 the momentum imparted in bubbly water by SWs increases monotonically, gradually levelling off at an air volume fraction of about 0.30. The experimental data are confirmed by two-dimensional (2D) simulation of SW propagation in bubbly water in terms of the SW velocity versus the air content, pressure profiles, as well as liquid and gas velocity behind the shock front.