Elastic wave propagation in anisotropic spherical curved plates

Spherical plate-like structures are used in pressure vessels, spherical domes of power plants, and in many other industrial applications. For non-destructive evaluation of such spherical structures, the mechanics of elastic wave propagation in spherical curved plates must be understood. The current literature shows some valuable studies on Rayleigh surface wave propagation in isotropic solids with spherical boundaries. However, the guided wave propagation problem in an anisotropic spherical curved plate, which has not been studied before, is solved for the first time in this paper.The wave propagation, in both isotropic and anisotropic spherical curved plates, is investigated. The differential equations of motion and the stress-free boundary conditions on the inner and outer surfaces of a hollow sphere are approximately solved by a general solution technique. This solution technique was successfully utilized by the authors for solving the wave propagation problem in cylindrical plates, in their earlier works. Dispersion curves for spherical plates made of isotropic aluminum, steel, and anisotropic composite material are presented as well.     

Affine transformations of the equations of the linear theory of elasticity

This paper deals with the general formulas of affine transformations that preserve invariance of the static equations of the linear theory of elasticity in the case of arbitrary anisotropic materials. The invariance of the equations with respect to affine transformations allows one to model a given anisotropic material by another material. All anisotropic materials are divided into classes of mutually congruent materials. The congruency conditions are obtained for orthotropic and isotropic materials and for orthotropic and transversely isotropic materials. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 4, pp. 124–134, July–August, 2006.     

Stress-strain state of an anisotropic plate with curved cracks and thin rigid inclusions

We solve the bending problem for an anisotropic plate with flaws like smooth curved nonoverlapping through cracks and rigid inclusions. The problem is solved by the method of Lekhnitskii complex potentials specified as Cauchy type integrals over the flaw contours with an unknown integrand density function. We use the Sokhotskii—Plemelj formulas to reduce the boundary-value problem to a system of singular integral equations with the additional conditions that the displacements in the plate are single-valued when going around the cut contours and the equilibrium conditions for stress-free rigid inclusions. After the singular integrals are approximated by the Gauss-Chebyshev quadrature formulas, the problem is reduced to solving a system of linear algebraic equations. We study the local stress distribution near flaw tips. We analyze the mutual influence of flaws on the stress distribution character near their vertices and compare the well-known solutions for isotropic plates with the solutions obtained by passing to the limit in the anisotropy parameters (“weakly anisotropic material”) and by using the method proposed here.     

Perturbation Formulas for Polarization Ratio and Phase Shift of Rayleigh Waves in Prestressed Anisotropic Media

This paper completes an earlier study (Tanuma and Man, Journal of Elasticity, 85, 21–37, 2006) where we derive a first-order perturbation formula for the phase velocity of Rayleigh waves that propagate along the free surface of a macroscopically homogeneous, anisotropic, prestressed half-space. We adopt the formulation of linear elasticity with initial stress and assume that the deviation of the prestressed anisotropic medium from a suitably-chosen, comparative, unstressed and isotropic state be small. No assumption, however, is made on the material anisotropy of the incremental elasticity tensor. With the help of the Stroh formalism, here we derive first-order perturbation formulas for the changes in polarization ratio and phase shift of Rayleigh waves from their respective comparative isotropic value. Examples are given, which show that the perturbation formulas for phase velocity and polarization ratio can serve as a starting point for investigations on the possible advantages of using Rayleigh-wave polarization, as compared with using wave speed, for acoustoelastic measurement of stress.      

The degeneration of image singularities from anisotropic materials to isotropic materials for an elastic half-plane

The degeneration of image singularities from an anisotropic material to an isotropic material for a half-plane is discussed in this study. The Green’s functions for anisotropic and isotropic half-planes with traction free boundary subjected to concentrated forces and dislocations have been obtained by many authors. It was commonly accepted that the solution of isotropic problem cannot be derived from anisotropic solutions. However, we believe that this possibility exists as we will demonstrate in this paper. Anisotropic materials include only image singularities of order O(1/r) (i.e., forces and dislocations) existing on image points. There are many image points for anisotropic materials and the locations of these image points depend on the material constants. However, isotropic materials have only one image point with higher order image singularities (O(1/r²), O(1/r³)). From the analysis provided in this study, it is found that the higher order image singularities for an isotropic half-plane are generated by combining the concentrated forces and dislocations when an anisotropic material degenerates to an isotropic material. The solutions of higher order image singularities for isotropic material are dependent. Therefore, these image singularities can be combined to form only three or four simpler image singularities acting on an image point of the isotropic material.     

An inclusion problem

An inclusion is a special region in a material, and this region experiences a transformation of the following nature. If the inclusion were free, then it would acquire a certain deformation with no stress arising in it; but since the inclusion is “pasted” into the material, this prevents free deformations and causes stresses arising in the inclusion itself and in the environment. Three systems of equations describing the problem are derived. For a space with a homogeneous isotropic matrix, an equivalent system of integral equations is obtained whose solution, for a homogeneous anisotropic ellipsoidal inclusion, is reduced to a system of linear algebraic equations. For the case where the moduli of elasticity in the inclusion and the homogeneous matrix coincide, an explicit solution for an inclusion of arbitrary shape is obtained.     

Constitutive behaviour of anisotropic materials under shock loading

Thermodynamically and mathematically consistent constitutive equations suitable for shock wave propagation in an anisotropic material are presented in this paper. Two fundamental tensors α_ij and β_ij which represent anisotropic material properties are defined and can be considered as generalisations of the Kronecker delta symbol, which plays the main role in the theory of isotropic materials. Using two fundamental tensors α_ij and β_ij, the concept of total generalised “pressure” and pressure corresponding to the thermodynamic (equation of state) response are redefined. The equation of state represents mathematical and physical generalisation of the classical Mie–Grüneisen equation of state for isotropic material and reduces to the Mie–Grüneisen equation of state in the limit of isotropy. Based on the generalised decomposition of the stress tensor, the modified equation of state for anisotropic materials, and the modified Hill criteria, combined with the associated flow rule, a system of constitutive equations suitable for shock wave propagation is formulated. The behaviour of aluminium alloy 7010-T6 under shock loading conditions is considered. A comparison of numerical simulations with existing experimental data shows good agreement of the general pulse shape, Hugoniot Elastic Limits (HELs), and Hugoniot stress levels, and suggests that the constitutive equations are performing satisfactorily. The results are presented and discussed, and future studies are outlined.     

An approach for obtaining approximate formulas for the Rayleigh wave velocity

In this paper, we introduce an approach for finding analytical approximate formulas for the Rayleigh wave velocity for isotropic elastic solids and anisotropic elastic media as well. The approach is based on the least-square principle. To demonstrate its application, we applied it in order to obtain an explanation for Bergmann’s approximation, the earliest known approximation of the Rayleigh wave velocity for isotropic elastic solids, and used it to establish a new approximation. By employing this approach, the best approximate polynomials of the second order of the cubic power and the quartic power in the interval [0, 1] were found. By using the best approximate polynomial of the second order of the cubic power, we derived an approximate formula for the Rayleigh wave speed in isotropic elastic solids which is slightly better than the one given recently by Rahman and Michelitsch by employing Lanczos’s approximation. Also by using this second order polynomial, analytical approximate expressions for orthotropic, incompressible and compressible elastic solids were found. For incompressible case, it is shown that the approximation is comparable with Rahman and Michelitsch’s approximation, while for the compressible case, it is shown that our approximate formulas are more accurate than Mozhaev’s ones. Remarkably, by using the best approximate polynomials of the second order of the cubic power and the quartic power in the interval [0, 1], we derived an approximate formula of the Rayleigh wave velocity in incompressible monoclinic materials, where the explicit exact formulas of the Rayleigh wave velocity so far are not available.     

Fractal effect and anisotropic constitutive model for concrete

An elasto-anisotropic damage constitutive model for concrete is developed in this work. Disregarding the coupling between the isotropic and the anisotropic damage, the isotropic damage variables are defined as functions of the microcrack fractal dimension, and the anisotropic parts are expressed by the lengths of cracks in concrete which various in different directions. The Helmholtz free energy is decomposed into the elastic deforming, damage and irreversible deforming components, with the last component used to replace the plastic deformation. Therefore the damage constitutive formulas for concrete are derived based on continuum damage mechanics. Evolution laws for both isotropic and anisotropic damage variables are derived, in which the anisotropic parts are obtained by modifying an empirical model. The critical fracture stress and the fracture toughness are investigated for materials with a single fractal crack based on the fractal geometry and the Griffith fracture criterion. Numerical computation is conducted for concrete under the uniaxial and the biaxial compression. The results indicate that the material stiffness degradation can be well addressed when the anisotropic damage is incorporated; the irreversible deformation is greatly related to the behavior of the descending branch beyond the peak load. The validation of the presented model is proofed by comparing results with the experimental data. This model provides an approach to link the macro properties of a material with its micro-structure change.     

Reflection and transmission of electromagnetic waves from a nonlinear anisotropic slab between two linear isotropic media

Summary In this paper, the author has studied the reflection from and transmission through a homogenous nonlinear anisotropic slab (anisotropy being due to an external magnetic field) bounded by two linear isotropic media. Nonlinear equations describing the growth of the two modes of propagation of an electromagnetic wave in the direction of the magnetic field, in an anisotropic nonlinear medium, have been set up and solved; the solutions have been used to obtain expressions for the reflected and transmitted components of the incident wave. The simpler problem of the reflection and transmission from an isotropic nonlinear slab has also been discussed as a special case.     

Solvability of Static Contact Problems with Coulomb Friction for Orthotropic Material

We calculate an exact upper bound for the magnitude of the coefficient of friction that ensures the existence of a solution to a static contact problem with Coulomb friction. The approach is based on a general existence result that is valid under the assumption that the coefficient of friction is bounded by a certain constant depending on the constants in two special trace type estimates for a half space domain. We calculate these constants for orthotropic material and two space dimensions with the help of a representation for a partial Fourier transform of the solution to the corresponding system of elasticity equations. The result is compared to the formula for general anisotropic material. The new bound for orthotropic material is significantly larger than the old one for general material, if the material is close to an isotropic material with Poisson ration greater than zero. For some cases the new bound can be even larger than one.      

Strain energy density bounds for linear anisotropic elastic materials

Upper and lower bounds are presented for the magnitude of the strain energy density in linear anisotropic elastic materials. One set of bounds is given in terms of the magnitude of the stress field, another in terms of the magnitude of the strain field. Explicit algebraic formulas are given for the bounds in the case of cubic, transversely isotropic, hexagonal and tetragonal symmetry. In the case of orthotropic symmetry the explicit bounds depend upon the solution of a cubic equation, and in the case of the monoclinic and triclinic symmetries, on the solution of sixth order equations.     

A Complete Solution of the Wave Equations for Transversely Isotropic Media

A transversely isotropic material in the sense of Green is considered. A complete solution in terms of retarded potential functions for the wave equations in transversely isotropic media is presented. In this paper we reduce the number of potential functions to only one, and we discuss the required conditions. As a special case, the torsionless and rotationally symmetric configuration with respect to the axis of symmetry of the material is discussed. The limiting case of elastostatics is cited, where the solution is reduced to the Lekhnitskii–Hu–Nowacki solution. The solution is simplified for the special case of isotropy. In this way, a new series of potential functions (to the best knowledge of the author) for the elastodynamics problem of isotropic materials is presented This solution is reduced to a special case of the Cauchy–Kovalevski–Somigliana solution, if the displacements satisfy specific conditions. Finally, Boggio's Theorem is generalized for transversely isotropic media which may be of interest to the reader beyond the present application. Dedicated to Morton E. Gurtin     

Analytical Solution for Isothermal Flow in a Shock Tube Containing Rigid Granular Material

Analytical solution of shock wave propagation in pure gas in a shock tube is usually addressed in gas dynamics. However, such a solution for granular media is complex due to the inclusion of parameters relating to particles configuration within the medium, which affect the balance equations. In this article, an analytical solution for isothermal shock wave propagation in an isotropic homogenous rigid granular material is presented, and a closed-form solution is obtained for the case of weak shock waves. Fluid mass and momentum equations are first written in wave and (mathematical) non-conservation forms. Afterwards by redefining the sound speed of the gas flowing inside the pores, an analytical solution is obtained using the classical method of characteristics, followed by Taylor’s series expansion based on the assumption of weak flow which finally led to explicit functions for velocity, density and pressure. The solution enables plotting gas velocity, density and pressure variations in the porous medium, which is of high interest in the design of granular shock isolators.     

Stress singularities in an anisotropic body of revolution

To fill the gap in the literature on the application of three-dimensional elasticity theory to geometrically induced stress singularities, this work develops asymptotic solutions for Williams-type stress singularities in bodies of revolution that are made of rectilinearly anisotropic materials. The Cartesian coordinate system used to describe the material properties differs from the coordinate system used to describe the geometry of a body of revolution, so the problems under consideration are very complicated. The eigenfunction expansion approach is combined with a power series solution technique to find the asymptotic solutions by directly solving the three-dimensional equilibrium equations in terms of the displacement components. The correctness of the proposed solution is verified by convergence studies and by comparisons with results obtained using closed-form characteristic equations for an isotropic body of revolution and using the commercial finite element program ABAQUS for orthotropic bodies of revolution. Thereafter, the solution is employed to comprehensively examine the singularities of bodies of revolution with different geometries, made of a single material or bi-materials, under different boundary conditions.     

WAVES IN PRE-STRETCHED INCOMPRESSIBLE SOFT ELECTROACTIVE CYLINDERS:EXACT SOLUTION

This paper studies wave propagation in a soft electroactive cylinder with an underlying finite deformation in the presence of an electric biasing field.Based on a recently proposed nonlinear framework for electroelasticity and the associated linear incremental theory,the basic equations governing the axisymmetric wave motion in the cylinder,which is subjected to homogeneous pre-stretches and pre-existing axial electric displacement,are presented when the electroactive material is isotropic and incompressible.Exact wave solution is then derived in terms of(modified) Bessel functions.For a prototype model of nonlinear electroactive material,illustrative numerical results are given.It is shown that the effect of pre-stretch and electric biasing field could be significant on the wave propagation characteristics.     

Topological derivative for an anisotropic and heterogeneous heat diffusion problem

The topological derivative measures the sensitivity of a given shape functional with respect to an infinitesimal singular domain perturbation. According to the literature, the topological derivative has been fully developed for a wide range of physical phenomenon modeled by partial differential equations, considering homogeneous and isotropic constitutive behavior. In fact, only a few works dealing with heterogeneous and anisotropic material behavior can be found in the literature, and, in general, the derived formulas are given in an abstract form. In this work, we derive the topological derivative in its closed form for the total potential energy associated to an anisotropic and heterogeneous heat diffusion problem, when a small circular inclusion of the same nature of the bulk phase is introduced at an arbitrary point of the domain. In addition, we provide a full mathematical justification for the derived formula and develop precise estimates for the remainders of the topological asymptotic expansion. Finally, the influence of the heterogeneity and anisotropy are shown through some numerical examples of heat conductor topology optimization.     

Plane strain dynamics of elastic solids with intrinsic boundary elasticity, with application to surface wave propagation

In this paper, in a development of the static theory derived by Steigmann and Ogden (Proc. Roy. Soc. London A 453 (1997) 853), we establish the equations of motion for a non-linearly elastic body in plane strain with an elastic surface coating on part or all of its boundary. The equations of (linearized) incremental motions superposed on a finite static deformation are then obtained and applied to the problem of (time-harmonic) surface wave propagation on a pre-stressed incompressible isotropic elastic half-space with a thin coating on its plane boundary. The secular equation for (dispersive) wave speeds is then obtained in respect of a general form of incompressible isotropic elastic strain-energy function for the bulk material and a general energy function for the coating material. Specialization of the form of strain-energy function enables the secular equation to be cast as a quartic equation and we therefore focus on this for illustrative purposes. An explicit form for the secular equation is thereby obtained. This involves a number of material parameters, including residual stress and moment in the properties of the coating. It is shown how this equation relates to previous work on waves in a half-space with an overlying thin layer set in the classical theory of isotropic elasticity and, in particular, the significant effect of omission of the rotatory inertia term, even at small wave numbers, is emphasized. Corresponding results for a membrane-type coating, for which the bending moment, inertia and residual moment terms are absent, are also obtained. Asymptotic formulas for the wave speed at large wave number (high frequency) are derived and it is shown how these results influence the character of the wave speed throughout the range of wave number values. A bifurcation criterion is obtained from the secular equation by setting the wave speed to zero, thereby generalizing the bifurcation results of Steigmann and Ogden (Proc. Roy. Soc. London A 453 (1997) 853) to the situation in which residual stress and moment are present in the coating. Numerical results which show the dependence of the wave speed on the various material parameters and the finite deformation are then described graphically. In particular, features which differ from those arising in the classical theory are highlighted.     

Finite Inhomogeneous Shearing Deformations of a Transversely Isotropic Incompressible Material

The present paper deals with finite inhomogeneous shearing deformations of a slab of a special anisotropic solid. Two cases according to the directions of the anisotropic director of the medium are examined. In one case the solution reduces to a quadrature and gives an exact deformation field for particular values of the material constants. In the other case an exact solution is obtained. All such solutions reduce to the same existing solution for the corresponding isotropic elastic material. This revised version was published online in August 2006 with corrections to the Cover Date.     

Theory of averaging of elastic fields in media with a monoclinic structure

Conclusions The proposed relations of averaging theory, together with complex Kolosov-Muskhelishvili potentials for isotropic matrices and Lekhnitskii potentials for rectilinearly anisotropic matrices with prismatic fillers, constitute a closed system of equations in the problem of determining the internal fields and the complete set of effective elastic constants of composite media with uniform external static stresses.By combining relations of the averaging theory and well-known solutions of boundary-value problems on the stress-state of an infinite medium with an individual inclusion, we can directly construct the solution of the problem of determining the macroscopic parameters of a composite system with an arbitrary structure.Conformal mapping of the external boundary of the determining element onto a unit circle is an efficient method of calculating contour integrals in averaging theory with a high degree of accuracy.When the initial terms are retained in an expansion of the complex potentials in degrees of inclusion interaction, it is possible to obtain approximate analytic formulas for all of the effective constants. In special cases, these formulas coincide with the asymptotic formulas found from the exact solutions.Institute of Mechanics, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Prikladnaya Mekhanika, Vol. 23, No. 1, pp. 3–18, January, 1987.