A high-order theory for functionally graded axially symmetric cylindrical shells

A high-order theory for functionally graded axially symmetric cylindrical shell based on expansion of the axially symmetric equations of elasticity for functionally graded materials into Legendre polynomials series has been developed. The axially symmetric equations of elasticity have been expanded into Legendre polynomials series in terms of a thickness coordinate. In the same way, functions that describe functionally graded relations has been also expanded. Thereby, all equations of elasticity including Hook’s law have been transformed to corresponding equations for coefficients of Legendre polynomials expansion. Then system of differential equations in terms of displacements and boundary conditions for the coefficients of Legendre polynomials expansion coefficients has been obtained. Cases of the first and second approximations have been considered in more details. For obtained boundary-value problems’ solution, a finite element has been used and numerical calculations have been done with COMSOL MULTIPHYSICS and MATLAB.     

A problem for a micropolar rectangular region in the fifth approximation

An isotropic micropolar two-dimensional region is considered. Several equations of the fifth approximation for displacements and rotations are derived in terms of moments with respect to the Legendre polynomials. Based on these equations, the solutions obtained in the framework of the micropolar theory are compared with the solutions obtained in the framework of the classical theory of elasticity.     

On One Technique of Constructing the General Solution to Equilibrium Equations for Nonthin Plates

The mathematical theory of plates based on the expansion of functions into Fourier series in terms of Legendre polynomials is used to state a method for determining the general solution to a system of equilibrium equations describing the stress–strain state of nonthin transversally isotropic plates. The state is assumed symmetric about the median plane     

Some Nonlinear Relations of the Mathematical Theory of Nonthin Shells

The expansion of functions into Fourier series in terms of Legendre polynomials is used to state some relations of the geometrically nonlinear theory of nonthin anisotropic shells with a variable thickness. A system of equilibrium equations and corresponding boundary conditions are constructed     

Equations of an Elastic Anisotropic Layer

Differential equations of an elastic orthotropic layer are constructed on the basis of expansion of the solutions of the elasticity theory in terms of the Legendre polynomials. The order of the system of differential equations is independent of the form of the boundary conditions on the layer surfaces, which allows a correct formulation of conditions on contact surfaces.     

Representations of elastic fields of circular dislocation and disclination loops in terms of spherical harmonics and their application to various problems of the theory of defects

Elastic fields of circular dislocation and disclination loops are represented in explicit form in terms of spherical harmonics, i.e. via series with Legendre and associated Legendre polynomials. Representations are obtained by expanding Lipschitz-Hankel integrals with two Bessel functions into Legendre series. Found representations are then applied to the solutions of elasticity boundary-value problems of the theory of defects and to the calculation of elastic fields of segmented spherical inclusions. In the framework of virtual circular dislocation–disclination loops technique, a general scheme to solving axisymmetric elasticity problems with boundary conditions specified on a sphere is given. New solutions for elastic fields of a twist disclination loop in a spherical particle and near a spherical pore are demonstrated. The easy and straightforward way for calculations of elastic fields of segmented spherical inclusion with uniaxial eigenstrain is shown.     

Solution of contact problems on the basis of a refined theory of plates and shells

Solutions of contact mixed boundary-value problems for a plate and for a cylindrical shell are given. These solutions are obtained with the use of equations for shells constructed by expanding solutions of elasticity theory equations with respect to the Legendre polynomials. Results of numerical simulations of the stress state in the vicinity of points with changing conditions on the frontal faces of the shell are presented. The results obtained are compared with analytical solutions of elasticity theory problems and with solutions obtained on the basis of the classical equations of the shell theory. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 169–176, September–October, 2008.     

Stress-strain state of non-thin plates and shells. Generalized theory (survey)

Conclusions Thus, this part of our survey has presented the main approaches that have been taken to the construction of two-dimensional (in terms of the space coordinates) equations of a generalized theory of plates and shells. The solutions of these equations represent a certain approximation of the solution of the initial three-dimensional problem. They are based on expansion of the sought functions into Fourier series in Legendre polynomials of the thickness coordinate. Studies completed on the basis of the given variants of plate and shell theory were systematized and analyzed. In terms of the method of its construction, the theory involves a regular process of replacing the solution of the three-dimensional problem by the solution (or sequence of solutions) of two-dimensional boundary-value problems or initial-boundary-value problems. Numerical results illustrating the convergence of the successive approximation were presented. It should be noted that to make comparison with the results of classical or applied theories, several of the studies cited here presented solutions of problems for thin plates and shells with allowance only for the initial terms of expansions of the stress and displacement components into base functions (Legendre polynomials).S. P. Timoshenko Institute of Mechanics, Ukrainian Academy of Sciences, Kiev. Translated from Prikladnaya Mekhanika, Vol. 29, No. 11, pp. 3–34, November, 1993.     

Nonlinear equations of elastic deformation of plates

A method for constructing nonlinear equations of elastic deformation of plates with boundary conditions for stresses and displacements at the face surfaces in an arbitrary coordinate system is proposed. The initial three–dimensional problem of the nonlinear theory of elasticity is reduced to a one–parameter sequence of two–dimensional problems by approximating the unknown functions by truncated series in Legendre polynomials. The same unknowns are approximated by different truncated series. In each approximation, a linearized system of equations whose differential order does not depend on the boundary conditions at the face surfaces which can be formulated in terms of stresses or displacements is obtained.     

Approximate solutions for the problem of liquid film flow over an unsteady stretching sheet with thermal radiation and magnetic field

The proposed method is based on replacement of the unknown function by a truncated series of the shifted Legendre polynomial expansion. An approximate formula of the integer derivative is introduced. Special attention is given to study the convergence analysis and derive an upper bound of the error for the presented approximate formula. The introduced method converts the proposed equation by means of collocation points to a system of algebraic equations with shifted Legendre coefficients. Thus, after solving this system of equations, the shifted Legendre coefficients are obtained. This efficient numerical method is used to solve the system of ordinary differential equations which describe the thin film flow and heat transfer with the effects of the thermal radiation, magnetic field, and slip velocity.     

Flow of a viscous fluid around a sphere of finite radius

We consider the stead y flow of a viscous fluid around a sphere of finiteradius in non-linear formulation. The equations of motion are written in non-dimensional terms. We seek their solution as the expansion of the unknown stream function in a series of powers of the Reynolds number, the coefficients of which are polynomials in associated Legendre functions of the first kind. Recurrence relations are given for the sequential determination of all coefficients. The velocity and pressure fields are determined. The drag is calculated. Numerical calculations are carried out.     

Vibrations of multilayer plates under the effect ofimpulse loads. Three-dimensional theory

A new analytical–numerical approach to investigation of the response of multilayerplates to impulse loading is described in this paper. The plates behaviour is described by theequations of the three-dimensional elasticity theory. According to the approach being proposed,the sought for functions included in the system of equations and the boundary and initialconditions are presented as Fourier series expansions in the tangential directions. The derivativesof these functions in the transverse direction are replaced by their finite-difference presentations.As a result of such transforms, the problem of vibration of a multilayer plate is reduced tointegration of a system of ordinary differential equations with constant coefficients. Integration isperformed by expansion into the Taylors series. The possibilities of the approach proposed andthe validity of results obtained is illustrated by several examples of calculating vibration processesand the processes of propagation of elastic waves. A comparison of the results obtained on thebasis of other approaches has been performed.     

Parametric finite-volume micromechanics of periodic materials with elastoplastic phases

The recently incorporated parametric mapping capability into the finite-volume direct averaging micromechanics (FVDAM) theory has produced a paradigm shift in the theory’s development. The use of quadrilateral subvolumes made possible by the mapping facilitates efficient modeling of microstructures with arbitrarily shaped heterogeneities, and eliminates artificial stress concentrations produced by the rectangular subvolumes employed in the standard version. Herein, the parametric FVDAM theory is extended to the inelastic domain by implementing additional formulation required to accommodate plastic and thermal loading. Two different approaches of implementing plasticity have been investigated. The first approach is based on the treatment employed in previous versions of the theory wherein plastic strain fields are represented by a series expansion in Legendre polynomials. The second approach is based on direct surface-averaging of plastic strains calculated at a number of collocation points along the quadrilateral subvolumes’ surfaces, and offers substantial simplification in the parametric finite-volume theory’s elastic–plastic framework. Moreover, substantial reductions in execution times without loss of accuracy are realized due to the elimination of redundant plastic strain calculations in the subvolumes’ interiors employed in the evaluation of the Legendre polynomial coefficients. Numerical studies demonstrate the advantages of the parametric FVDAM theory relative to the standard version, together with new results that highlight its modeling capabilities vis-a-vis an emerging class of periodic lamellar materials with wavy microstructures and the thus-far undocumented architectural effects amplified by plasticity.     

Nonclassical models of the theory of plates and shells

Publications dealing with the study of methods of reducing a three-dimensional problem of the elasticity theory to a two-dimensional problem of the theory of plates and shells are reviewed. Two approaches are considered: the use of kinematic and force hypotheses and expansion of solutions of the three-dimensional elasticity theory in terms of the complete system of functions. Papers where a three-dimensional problem is reduced to a two-dimensional problem with the use of several approximations of each sought function (stresses and displacements) by segments of Legendre polynomials are also reviewed.     

On new symplectic elasticity approach for exact bending solutions of rectangular thin plates with two opposite sides simply supported

This paper presents a bridging research between a modeling methodology in quantum mechanics/relativity and elasticity. Using the symplectic method commonly applied in quantum mechanics and relativity, a new symplectic elasticity approach is developed for deriving exact analytical solutions to some basic problems in solid mechanics and elasticity which have long been bottlenecks in the history of elasticity. In specific, it is applied to bending of rectangular thin plates where exact solutions are hitherto unavailable. It employs the Hamiltonian principle with Legendre’s transformation. Analytical bending solutions could be obtained by eigenvalue analysis and expansion of eigenfunctions. Here, bending analysis requires the solving of an eigenvalue equation unlike in classical mechanics where eigenvalue analysis is only required in vibration and buckling problems. Furthermore, unlike the semi-inverse approaches in classical plate analysis employed by Timoshenko and others such as Navier’s solution, Levy’s solution, Rayleigh–Ritz method, etc. where a trial deflection function is pre-determined, this new symplectic plate analysis is completely rational without any guess functions and yet it renders exact solutions beyond the scope of applicability of the semi-inverse approaches. In short, the symplectic plate analysis developed in this paper presents a breakthrough in analytical mechanics in which an area previously unaccountable by Timoshenko’s plate theory and the likes has been trespassed. Here, examples for plates with selected boundary conditions are solved and the exact solutions discussed. Comparison with the classical solutions shows excellent agreement. As the derivation of this new approach is fundamental, further research can be conducted not only on other types of boundary conditions, but also for thick plates as well as vibration, buckling, wave propagation, etc.     

Polynomial approach for modeling a piezoelectric disc resonator partially covered with electrodes

The frequency spectrum of a partially metallized piezoelectric disc resonator was studied using Legendre polynomials. The formulation, based on three-dimensional equations of linear elasticity, takes into account the high contrast between the electroded and non-electroded regions. The mechanical displacement components and the electrical potential were expanded in a double series of orthonormal functions and were introduced into the equations governing wave propagation in piezoelectric media. The boundary and continuity conditions were automatically incorporated into the equations of motion by assuming position-dependent physical material constants or delta-functions. The incorporation of electrical sources is illustrated. Structure symmetry was used to reduce the number of unknowns. The vibration characteristics of the piezoelectric discs were analyzed using a three-dimensional modelling approach with modal and harmonic analyses. The numerical results are presented as resonance and anti-resonance frequencies, electric input admittance, electromechanical coupling coefficient and field profiles of fully and partially metallized PIC151 and PZT5A resonator discs. In order to validate our model, the results obtained were compared with those published previously and those obtained using an analytical approach.     

Stress concentrations in the particulate composite with transversely isotropic phases

The accurate series solution have been obtained of the elasticity theory problem for a transversely isotropic solid containing a finite or infinite periodic array of anisotropic spherical inclusions. The method of solution has been developed based on the multipole expansion technique. The basic idea of method consists in expansion the displacement vector into a series over the set of vectorial functions satisfying the governing equations of elastic equilibrium. The re-expansion formulae derived for these functions provide exact satisfaction of the interfacial boundary conditions. As a result, the primary spatial boundary-value problem is reduced to an infinite set of linear algebraic equations. The method has been applied systematically to solve for three models of composite, namely a single inclusion, a finite array of inclusions and an infinite periodic array of inclusions, respectively, embedded in a transversely isotropic solid. The numerical results are presented demonstrating that elastic properties mismatch, anisotropy degree, orientation of the anisotropy axes and interactions between the inclusions can produce significant local stress concentration and, thus, affect greatly the overall elastic behavior of composite.     

Equations of motion and boundary conditions of physical meaning of micropolar theory of thin bodies with two small cuts

Nowadays, microcontinuous mechanics (mechanics of media with microstructure) is being developed very intensively, which is testified by recently published papers [1–14] and by many others, as well as by the symposiumdedicated to the hundredth anniversary of the brothers Cosserat monograph [15], held inParis in 2009. A survey of foreign papers is given in [16], and a special place is occupied by earlier publications of Soviet scientists on micropolar theory of elasticity [17–24]. A brief survey of Cosserat theory of elasticity and an analysis and prospects of such theories in mechanics of rigid deformable bodies is given in [21]. It should be noted that, in a majority of cases, the structure strength calculations are based on the classical theory of elasticity. But there are materials such as animal bones, graphite, several polymers, polyurethane films, porous materials (pumice), various synthetic materials, and materials with inclusions which, under certain conditions, exhibit micropolar properties. There are effects which cannot be prescribed by the classical theory. In statics, nonclassical behavior can be observed in bending of thin films and cantilevers, in torsion of thin and thin-walled rods, and in the case of stress concentration near holes, corner points, cracks, and inclusions. For example, thin specimens are more rigid in bending and torsion as is prescribed by the classical theory [25–27]. The stress concentration near holes decreases, and the concentration factor depends on the radius [28]. The stress concentration near cracks also becomes lower. Conversely, the stress concentration near inclusions is higher than predicted by the classical theory [29–31]. If the material has no center of symmetry of elastic properties, then calculations according to the micropolar theory shows that the specimen is twisted in tension [32]. In dynamical problems, several phenomena also differ from the classical concepts. For example, shear waves propagate with dispersion, microrotation waves arise, and the vibration natural modes differ from the classical ones [2, 7, 11–13, 33]. All these phenomena are used to determine material constants of the micropolar theory of elasticity. There are many methods for determining such constants [2, 34]. Since thin bodies (one-, two-, three-, and multilayer structures) are widely used, it is necessary to create new refined microcontinual theories of thin bodies and advanced methods for their computations. In the present paper, various representations of the system of equations of motion are obtained in the micropolar theory of thin bodies with two small parameters in momenta with respect to a system of Legendre polynomials in the case where an arbitrary line is taken for the base. In this connection, a vector parametric equation of the region of a thin body is given for the parametrization under study, different families of bases (frames) are introduced, and expressions for components of the unit tensor of rank two (UTRT) are obtained. Representations of gradient, tensor divergence, equations of motion, and boundary conditions for the considered parametrization are given. Definitions of (m, n)th-order moment of a variable with respect to an arbitrary system of orthogonal polynomials and a system of Legendre polynomials is given. Expressions for themoments of partial derivatives and several expressions with respect to a system of Legendre polynomials and boundary conditions in moments are obtained.     

Two-dimensional elasticity solution for bending of functionally graded beams with variable thickness

This paper studies the stress and displacement distributions of functionally graded beam with continuously varying thickness, which is simply supported at two ends. The Young’s modulus is graded through the thickness following the exponential-law and the Poisson’s ratio keeps constant. On the basis of two-dimensional elasticity theory, the general expressions for the displacements and stresses of the beam under static loads, which exactly satisfy the governing differential equations and the simply supported boundary conditions at two ends, are analytically derived out. The unknown coefficients in the solutions are approximately determined by using the Fourier sinusoidal series expansions to the boundary conditions on the upper and lower surfaces of the beams. The effect of Young’s modulus varying rules on the displacements and stresses of functionally graded beams is investigated in detail. The two-dimensional elasticity solution obtained can be used to assess the validity of various approximate solutions and numerical methods for the aforementioned functionally graded beams.     

Representation of solutions to equations of hyperbolic type

The general solutions to hyperbolic equations of fourth and sixth order are obtained using Vekua’s method for the representation of the general solutions to elliptic equations of order 2n with the aid of n analytic functions. It is assumed that the right-hand sides of the hyperbolic equations can be expanded in time series of sines. The systems of equations of various approximations for a prismatic thin body in terms of moments with respect to the system of Legendre polynomials can be reduced to these equations and to some hyperbolic-type equations of higher order.