Rayleigh and Love surface waves in isotropic media with negative Poisson’s ratio

The behavior of Rayleigh surface waves and the first mode of the Love waves in isotropic media with positive and negative Poisson’s ratio is compared. It is shown that the Rayleigh wave velocity increases with decreasing Poisson’s ratio, and it increases especially rapidly for negative Poisson’s ratios less than ?0.75. It is demonstrated that, for positive Poisson’s ratios, the vertical component of the Rayleigh wave displacements decays with depth after some initial increase, while for negative Poisson’s ratios, there is a monotone decrease. The Rayleigh waves are characterized by elliptic trajectories of the particle motion with the change of the rotation direction at critical depths and by the linear vertical polarization at these depths. It is found that the elliptic orbits are less elongated and the critical depths are greater for negative Poisson’s ratios. It is shown that the stress distribution in the Rayleighwaves varies nonmonotonically with the dimensionless depth as (positive or negative) Poisson’s ratio varies. The stresses increase strongly only as Poisson’s ratio tends to?1. It is shown that, in the case of an incompressible thin covering layer, the velocity of the first mode of the Love waves strongly increases for negative Poisson’s ratios of the half-space material. If the thickness of the incompressible layer is large, then the wave very weakly penetrates into the halfspace for any value of its Poisson’s ratio. For negative Poisson’s ratios, the Love wave in a layer and a half-space is mainly localized in the covering layer for any values of its thickness and weakly penetrates into the half-space. For the first mode of the Love waves, it was discovered that there is a strong increase in the maximum of one of the shear stresses on the interface between the covering layer and the half-space as Poisson’s ratios of both materials decrease. For the other shear stress, there is a stress jump on the interface and a more complicated dependence of the stress on Poisson’s ratio on both sides of the interface.     

FE analysis of the in-plane mechanical properties of a novel Voronoi-type lattice with positive and negative Poisson’s ratio configurations

This work presents a novel formulation for a Voronoi-type cellular material with in-plane anisotropic behaviour, showing global positive and negative Poisson’s ratio effects under uniaxial tensile loading. The effects of the cell geometry and relative density over the global stiffness, equivalent in-plane Poisson’s ratios and shear modulus of the Voronoi-type structure are evaluated with a parametric analysis. Empirical formulas are identified to reproduce the mechanical trends of the equivalent homogeneous orthotropic material representing the Voronoi-type structure and its geometry parameters.     

Parametric identification of crystals having a cubic lattice with negative Poisson’s ratios

A two-dimensional model of an anisotropic crystalline material with cubic symmetry is considered. This model consists of a square lattice of round rigid particles, each possessing two translational and one rotational degree of freedom. Differential equations that describe propagation of elastic and rotational waves in such a medium are derived. A relationship between three groups of parameters is found: second-order elastic constants, acoustic wave velocities, and microstructure parameters. Values of the microstructure parameters of the considered anisotropic material at which its Poisson’s ratios become negative are found.     

Determining the Mie potential parameters using Poisson’s ratio

The known value of Poisson’s ratio specifying the relation between the strains along the principal directions in the case of uniaxial strain is used to propose an approach to derive an equation relating this ratio to the exponents of the Mie pair potential. An example of determining one of these exponents is discussed when the other exponent is given.     

A new constitutive equation for solid propellant with the effects of aging and viscoelastic Poisson’s ratio

A new constitutive equation for solid propellant with the effects of aging and viscoelastic Poisson’s ratio is proposed. Effects of thermo-oxidative aging and viscoelastic Poisson’s ratio are considered in this comprehensive constitutive equation with two sets of reduced time system coping with the time and temperature dependence. In order to simulate the single and combined effects of aging and viscoelastic Poisson’s ratio, constitutive equation is rewritten into an incremental form and implemented in the user subroutine UMAT at the platform of finite element code ABAQUS. Detailed procedure for acquiring the parameters in constitutive equation is introduced and conducted for the subsequent applied analysis. Two typical loading cases during the service life of solid rocket motor and four sets of combined constitutive models are simulated. Von Mises strain and stress distribution and their changes versus time are utilized as the main analysis index. The results show that the effects of aging and viscoelastic Poisson’s ratio or their combinations will improve or decrease the level and change the distribution of Von Mises strain and stress in varying degrees.     

Using the simple compression test to determine Young’s modulus,Poisson’s ratio and the Coulomb friction coefficient

Analytical solutions are derived for the compression of cylinders with bonded surfaces and with Coulomb friction conditions at the interfaces. The bonded solution assumes that the radial displacement is linearly dependent on radius which leads to simple forms. These are compared with FE data and the apparent modulus is found to be within about 8% for the whole range of aspect ratios (10^?2–10³), and thus degrees of constraint for the cylinders. The apparent moduli are shown to be strong functions of both ν and μ and the solutions thus provide schemes for finding both parameters experimentally using inverse methods. This is demonstrated by using the FE results as such data to explore how many tests, and what aspect ratios, are needed. Some preliminary experimental results are also given.     

Anti-plane shear Green’s functions for an isotropic elastic half-space with a material surface

This paper presents analytical Green’s function solutions for an isotropic elastic half-space subject to anti-plane shear deformation. The boundary of the half-space is modeled as a material surface, for which the Gurtin–Murdoch theory for surface elasticity is employed. By using Fourier cosine transform, analytical solutions for a point force applied both in the interior or on the boundary of the half-space are derived in terms of two particular integrals. Through simple numerical examples, it is shown that the surface elasticity has an important influence on the elastic field in the half-space. The present Green’s functions can be used in boundary element method analysis of more complicated problems.     

Green’s functions of a surface-stiffened transversely isotropic half-space

Green’s functions of a transversely isotropic half-space overlaid by a thin coating layer are analytically obtained. The surface coating is modeled by a Kirchhoff thin plate perfectly bonded to the half-space. With the aid of superposition technique and employing appropriate displacement potential functions, the Green’s functions are expressed in two parts; (i) a closed-form part corresponding to the transversely isotropic half-space with surface kinematic constraints, and (ii) a numerically evaluated part reflecting the interaction between the half-space and the plate in the form of semi-infinite integrals. Some limiting cases of the problem such as surface-stiffened isotropic half-space, Boussinesq and Cerruti loadings, and extremely flexible and rigid plates are also studied. For the classical Cerruti problem in transversely isotropic materials, the effects of incompressibility are highlighted. Numerical results are provided to show the effects of material anisotropy, relative stiffness factor, and load buried depth. The obtained Green’s functions play a key role in treating further mixed-boundary-value problems in surface stiffened transversely isotropic half-spaces.     

Three-dimensional Green’s functions for transversely isotropic thermoelastic bimaterials

Green’s functions for transversely isotropic thermoelastic biomaterials are established in the paper. We first express the compact general solutions of transversely isotropic thermoelastic material in terms of harmonic functions and introduce six new harmonic functions. The three-dimensional Green’s function having a concentrated heat source in steady state is completely solved using these new harmonic functions. The analytical results show some new phenomena of temperature and stress distributions at the interface. The temperature contours are normal to the interface for the isotropic material but not for the orthotropic one. The normal stress contours are parallel to the interface at the boundary in the isotropic region only and shear failure is most likely at the heat source due to the highly degenerated direction of shear stress contours.     

Influence of the spatial variation of Poisson’s ratio upon the elastic field in nonhomogeneous axisymmetric bodies

The effect of a nonconstant Poisson’s ratio upon the elastic field in functionally graded axisymmetric solids is analyzed. Both of the elastic coefficients, i.e. Young’s modulus and Poisson’s ratio, are permitted to vary in the radial direction. These elastic coefficients are considered to be functions of composition and are related on this basis. This allows a closed form solution for the stress function to be obtained. Two cases are discussed in this investigation: first, both Young’s modulus and Poisson’s ratio are allowed to vary across the radius and the effect of spatial variation of Poisson’s ratio upon the maximum radial displacement is investigated; secondly, Young’s modulus is taken as constant and the change in the maximum hoop stress resulting from a variable Poisson’s ratio is calculated.     

Green’s functions of an exponentially graded transversely isotropic half-space

By virtue of a complete set of displacement potential functions and Hankel transform, the analytical expressions of Green’s function of an exponentially graded elastic transversely isotropic half-space is presented. The given solution is analytically in exact agreement with the existing solution for a homogeneous transversely isotropic half-space. Employing a robust asymptotic decomposition technique, the Green’s function is decomposed to the closed-form Green’s function corresponding to the homogeneous transversely isotropic half-space and grading term with strong decaying integrands. This representation is very useful for numerical methods which are based on boundary-integral formulations such as boundary-element method since the numerically evaluated part is not responsible for the singularity. The high accuracy of the proposed numerical scheme is confirmed by some numerical examples.     

Very Large Poisson’s Ratio with a Bounded Transverse Strain in Anisotropic Elastic Materials

In a recent paper by Ting and Chen [18] it was shown by examples that Poisson’s ratio can have no bounds for all anisotropic elastic materials. With the exception of cubic materials, the examples presented involve a very large transverse strain. We show here that a very large Poisson’s ratio with a bounded transverse strain exists for all anisotropic elastic materials. The large Poisson’s ratio with a bounded transverse strain occurs when the axial strain is in the direction very near or at the direction along which Young’s modulus is very large. In fact the transverse strain has to be very small for the material to be stable. If the non-dimensionalized Young’s modulus is of the order δ⁻¹, where δ is very small, the axial strain, the transverse strain and Poisson’s ratio are of the order δ, δ^1/2 and δ^−1/2, respectively. Mathematics Subject Classifications (2000) 74B05, 74E10.T.C.T. Ting: Professor Emeritus of University of Illinois at Chicago and Consulting Professor of Stanford University.     

Three-dimensional dynamic Green’s functions for a multilayered transversely isotropic half-space

With the aid of a method of displacement potentials, an efficient and accurate analytical derivation of the three-dimensional dynamic Green’s functions for a transversely isotropic multilayered half-space is presented. Constituted by proper algebraic factorizations, a set of generalized transmission–reflection matrices and internal source fields that are free of any numerically unstable exponential terms are proposed for effective computations of the potential solution. Three-dimensional point-load Green’s functions for stresses and displacements are given, for the first time, in the complex-plane line-integral representations. The present formulations and solutions are analytically in exact agreement with the existing solutions given by Pak and Guzina (2002) for the isotropic case. For the numerical computation of the integrals, a robust and effective methodology which gives the necessary account of the presence of singularities including branch points and poles on the path of integration is laid out. A comparison with the existing numerical solutions for multilayered isotropic half-space is made to confirm the accuracy of the numerical solutions.     

Extrema of Young’s modulus for cubic and transversely isotropic solids

For a homogeneous anisotropic and linearly elastic solid, the general expression of Young’s modulus E(n), embracing all classes that characterize the anisotropy, is given. A constrained extremum problem is then formulated for the evaluation of those directions n at which E(n) attains stationary values. Cubic and transversely isotropic symmetry classes are dealt with, and explicit solutions for such directions n are provided. For each case, relevant properties of these directions and corresponding values of the modulus are discussed as well. Results are shown in terms of suitable combinations of elements of the elastic tensor that embody the discrepancy from isotropy. On the basis of such material parameters, for cubic symmetry two classes of behavior can be distinguished and, in the case of transversely isotropic solids, the classes are found to be four. For both symmetries and for each class of behavior, some examples for real materials are shown and graphical representations of the dependence of Young’s modulus on direction n are given as well.     

Resonant-Based Identification of the Poisson’s Ratio of Orthotropic Materials

The resonant-based identification of the in-plane elastic properties of orthotropic materials implies the estimation of four principal elastic parameters: E ₁ , E ₂ , G ₁₂ , and ν ₁₂ . The two elastic moduli and the shear modulus can easily be derived from the resonant frequencies of the flexural and torsional vibration modes, respectively. The identification of the Poisson’s ratio, however, is much more challenging, since most frequencies are not sufficiently sensitive to it. The present work addresses this problem by determining the test specimen specifications that create the optimal conditions for the identification of the Poisson’s ratio. Two methods are suggested for the determination of the Poisson’s ratio of orthotropic materials: the first employs the resonant frequencies of a plate-shaped specimen, while the second uses the resonant frequencies of a set of beam-shaped specimens. Both methods are experimentally validated using a stainless steel sheet.     

Maxwell’s Equations with Vector Hysteresis

Electromagnetic processes in magnetic materials are described by Maxwells equations. In ferrimagnetic insulators, assuming that D = E, we have the equationIn ferromagnetic metals, neglecting displacement currents and assuming Ohms law, we instead getAlternatively, under quasi-stationary conditions, for either material we can also deal with the magnetostatic equations:(Here f_ext and J_ext are prescribed time-dependent fields.) In any of these settings, the dependence of M on H is represented by a constitutive law accounting for hysteresis: M= (H), being a vector extension of the relay model. This is characterized by a rectangular hysteresis loop in a prescribed x-dependent direction, and accounts for high anisotropy and nonhomogeneity. The discontinuity in this constitutive relation corresponds to the possible occurrence of free boundaries.Weak formulations are provided for Cauchy problems associated with the above equations; existence of a solution is proved via approximation by time-discretization, derivation of energy-type estimates, and passage to the limit. An analogous representation is given for hysteresis in the dependence of P on E in ferroelectric materials. A model accounting for coupled ferrimagnetic and ferroelectric hysteresis is considered, too.Acknowledgement This research was partly supported by the project Free boundary problems in applied sciences of Italian M.I.U.R.. I gratefully acknowledge the useful suggestions from the reviewers.     

On Saint-Venant’s Problem with Concentrated Loads

The classical Saint-Venant problem is to find a solution of the traction problem of elastostatics in a finite cylinder ?? loaded over its bases. We prove that the problem has a unique solution for equilibrated surface forces $\hat{ \boldsymbol { s}}\in W^{-1,q}(\partial\Omega)$ , with q??(2?? ₀,+??) for some positive ? ₀ depending on ??. Hence $\hat{ \boldsymbol { s}}$ can model force acting on ???, concentrated on sets of zero Lebesgue surface measure of ???. Moreover, if $\hat{ \boldsymbol { s}}$ is equilibrated on each basis, we give a simple proof of the Toupin estimate expressing Saint-Venant??s principle.     

Image singularities of Green’s functions for isotropic elastic bimaterials subjected to concentrated forces and dislocations

Green’s functions for isotropic materials in the two-dimensional problem for elastic bimaterials with perfectly bonded interface are reexamined in the present study. Although the Green’s function for an isotropic elastic bimaterial subjected to a line force or a line dislocation has been discussed by many authors, the physical meaning and the structure of the solution are not clear. In this investigation, the Green’s function for an elastic bimaterial is shown to consist of eight Green’s functions for a homogeneous infinite plane. One of the novel features is that Green’s functions for bimaterials can be expressed directly by knowing Green’s functions for the infinite plane. If the applied load is located in material 1, the solution for the half-plane of material 1 is constructed with the help of five Green’s functions corresponding to the infinite plane. However, the solution for the half-plane of material 2 only consists of three Green’s functions for the infinite plane. One of the five Green’s functions of material 1 and all the three Green’s functions of material 2 have their singularities located in the half-plane where the load is applied, and the other four image singularities of material 1 are located outside the half-plane at the same distance from the interface as that of the applied load. The nature and magnitude of the image singularities for both materials are presented explicitly from the principle of superposition, and classified according to different loads. It is known that for the problem of anisotropic bimaterials subjected to concentrated forces and dislocations, the image singularities are simply concentrated forces and dislocations with the stress singularity of order O(1/r). However, higher orders (O(1/r²) and O(1/r³)) of stress singularities are found to exist in this study for isotropic bimaterials. The highest order of the stress singularity is O(1/r³) for the image singularities of material 1, and is O(1/r²) for material 2. Using the present solution, Green’s functions associated with the problems of elastic half-plane with free and rigidly fixed boundaries, for homogeneous isotropic elastic solid, are obtained as special cases.     

Flow structures of turbulent Rayleigh–Bénard convection in annular cells with aspect ratio one and larger

Acta Mechanica Sinica - We present an experimental study of flow structures in turbulent Rayleigh–Bénard convection in annular cells of aspect ratios $$\varGamma =1$$ , 2 and 4, and...