Free nonaxisymmetric vibrations of radially polarized hollow piezoceramic cylinders of finite length

The three-dimensional problem of free nonaxisymmetric vibrations of hollow piezoceramic cylinders with axial polarization is considered. An efficient numerical analytic method to solve boundary-value problems is proposed. The original three-dimensional problem of electroelasticity is reduced to a two-dimensional problem by representing the displacement components as standing circumferential waves. Spline collocation with respect to the axial coordinate is used to reduce this two-dimensional problem to an eigenvalue boundary-value problem with respect to the radial coordinate. This problem is solved by the stable discrete-orthogonalization and incremental-search methods. Numerical results are presented and the natural frequencies of the cylinders are analyzed in a wide range of their geometric characteristics     

Free vibrations of axially polarized piezoceramic hollow cylinders of finite length

The problem of the free axisymmetric vibrations of longitudinally polarized piezoceramic hollow cylinders is solved by a numerical analytic method. The spline-collocation method with respect to the longitudinal coordinate is used to reduce the original problem of electroelasticity to an eigenvalue boundary-value problem for ordinary differential equations with respect to the radial coordinate. This problem is solved by the stable discrete-orthogonalization and incremental search methods. Numerical results are presented and the natural frequencies of the cylinders are analyzed for a wide range of their geometric characteristics     

Nonaxisymmetric electroelastic vibrations of a hollow sphere made of functionally gradient piezoelectric material

The nonaxisymmetric problem of natural vibrations of a hollow sphere made of functionally gradient piezoelectric material is solved based on 3D electroelasticity. The properties of the material change continuously along a radial coordinate according to an exponential law. The external surface of the sphere is free of tractions and either insulated or short-circuited by electrodes. After separation of variables and representation of the components of the displacements and of the stress tensor in terms of spherical functions, the initially three-dimensional problem is reduced to a boundary-value problem for the eigenvalues expressed by ordinary differential equations. This problem is solved by a stable discrete-orthogonalization technique in combination with a step-by-step search method with respect to the radial coordinate. Moreover, a numerical investigation is performed based on the algorithm used for solving the problem. In particular, we investigate the influence of the geometric and electric parameters on the frequency spectrum at the nonaxisymmetry of natural vibrations of an inhomogeneous piezoceramic thick-walled sphere.     

Free axisymmetric vibrations of solid cylinders: numerical problem solving

The problem of free vibrations of a solid cylinder with different boundary conditions is solved using the three-dimensional theory of elasticity and a numerical analytic approach. The spline-approximation and collocation methods are used to reduce the partial differential equations of elasticity to systems of ordinary differential equations of high order with respect to the radial coordinate. These equations are solved by stable numerical discrete orthogonalization and incremental search. Calculated results are presented for transversely isotropic and inhomogeneous materials of the cylinder and for several types of boundary conditions at its ends     

Using spline-approximation to solve problems of axisymmetric free vibration of thick-walled orthotropic cylinders

The three-dimensional theory of elasticity is used to study the free vibrations of an anisotropic hollow cylinder with different boundary conditions at the ends. The relevant problem is solved by a numerical-and-analytic method. Spline approximation and collocation is used to reduce the partial differential equations of elasticity to a boundary-value problem for a system of ordinary differential equations of high order for the radial coordinate, which is solved using the stable discrete-orthogonalization and incremental-search methods. The calculated results for an orthotropic inhomogeneous cylinder with boundary conditions of several types are presented Translated from Prikladnaya Mekhanika, Vol. 44, No. 10, pp. 74–85, October 2008.     

厚壳三维分析的半解析解

本文对座标系三维弹性力学问题采用周向与径向解析,轴向离散的半解析数值方法。通过引入部分解析函数,将三维问题归结为一维离散化方程。这种方法能适应于一大类复杂的弹性力学问题,方法简单,计算工作量少。本文用这方法来分析厚壳的三维变形与应力规律,研究大厚跨比的强厚壳的三维弹性理论解,为建立可靠的强厚壳理论提供依据。     

Conditions for Strong Ellipticity of Anisotropic Elastic Materials

We consider the problem of finding the transversely isotropic elasticity tensor closest to a given elasticity tensor with respect to the Frobenius norm. A similar problem was considered by other authors and solved analytically assuming a fixed orientation of the natural coordinate system of the transversely isotropic tensor. In this paper we formulate a method for finding the optimal orientation of the coordinate system—the one that produces the shortest distance. The optimization problem reduces to finding the absolute maximum of a homogeneous eighth-degree polynomial on a two-dimensional sphere. This formulation allows us a convenient visualization of local extrema, and enables us to find the closest transversely isotropic tensor numerically.      

Stability of viscous flow between rotating and stationary disks

The linear problem of the stability of viscous flow between rotating and stationary parallel disks is solved in the locally homogeneous formulation using the method of normal modes. The main flow is assumed to be selfsimilar with respect to the radial coordinate. The system of sixth-order equations, derived for the amplitude functions of the disturbances, is integrated by a finite difference method. The stability characteristics with respect to disturbances of four types are calculated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 79–87, November–December, 1991.     

Effect of weak rotation on natural convection in a horizontal porous cylinder

Natural convection in a fluid saturated porous medium confined in a horizontal circular cylinder and rotating about its axis, with isothermal boundary conditions and uniform internal heat sink, is studied by both numerical and perturbation methods. No symmetry with respect to the vertical diameter is expected for the flow and temperature fields and the whole region must be involved in the computation. Only the weak rotation regime, for which the centrifugal force is negligible compared to gravity, is considered. Governing equations for the two-dimensional flow field are solved in both rotating and non-rotating coordinate systems. Results indicate that rotation significantly decreases the radial amplitude of fluid particle trajectories in the radial direction and thus reduces the overall heat transfer.     

丝束轴向变角度复合材料梁的弹性分析

丝束变角度复合材料具有变刚度的特点,因此其结构分析具有相当难度．本文采用状态空间法和微分求积法联合的半解析数值方法对丝束轴向变角度复合材料梁的弯曲问题进行研究．假设纤维方向角沿梁的轴向按照任意连续函数变化,选取位移和位移的一阶导数作为状态变量,建立了丝束轴向变角度复合材料梁弹性分析的状态空间方程,将状态变量对轴向坐标的导数采用微分求积法进行求解,进而可得问题的半解析数值解．通过与现有文献及ABAQUS计算结果的比较,验证了本文方法的正确性,并对微分求积法求解本问题的收敛性进行了分析．通过数值算例研究了纤维方向角沿梁轴向的变化对丝束轴向变角度复合材料梁的位移及应力分布的影响,研究结果可为该种结构的设计提供一定的参考．     

The stress state of anisotropic conic shells with thickness varying in two directions

The problem on the stress state of thin layered conic shells consisting of rigidly joint layers is considered. The material of the layers has one plane of elastic symmetry. The thickness of each layer in the meridional direction varies linearly in such a manner that it is proportional to the distance from the axis of rotation to the coordinate surface and varies arbitrarily in the circumferential direction. The desired functions are represented by the product of a power function of the radial coordinate and the unknown function of the central angle in the cross section. This makes it possible to separate out the radial coordinate and to derive the resolving system of ordinary differential equations, which is solved numerically. An example of the solution of a specific problem is given. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev. Translated from Prikladnaya Mekhanika, Vol. 36, No. 5, pp. 81–88, May, 2000.     

Thermoelastoplastic Nonaxisymmetric Stress-Strain Analysis of Laminated Shells of Revolution

A technique for nonaxisymmetric thermoelastoplastic stress–strain analysis of laminated shells of revolution is developed. It is assumed that there is no slippage and the layers are not separated. The problem is solved using the geometrically linear theory of shells based on the Kirchhoff–Love hypotheses. The thermoplastic relations are written down in the form of the method of elastic solutions. The order of the system of partial differential equations obtained is reduced by means of trigonometric series in the circumferential coordinate. The systems of ordinary differential equations thus obtained are solved by Godunov's discrete-orthogonalization method. The nonaxisymmetric thermoelastoplastic stress–strain state of a two-layered shell is analyzed as an example     

Scattering by an anisotropic circle

The scattering by a circle is considered when the outside medium is isotropic and the inside medium is anisotropic (orthotropic). The problem is a scalar one and is phrased as a scattering problem for elastic waves with polarization out of the plane of the circle (SH wave), but the solution is with minor modifications valid also for scattering of electromagnetic waves. The equation inside the circle is first transformed to polar coordinates and it then explicitly contains the azimuthal angle through trigonometric functions. Making an expansion in a trigonometric series in the azimuthal coordinate then gives a coupled system of ordinary differential equations in the radial coordinate that is solved by power series expansions. With the solution inside the circle complete the scattering problem is solved essentially as in the classical case. Some numerical examples are given showing the influence of anisotropy, and it is noted that the effects of anisotropy are generally strong except at low frequencies where the dominating scattering only depends on the mean stiffness and not on the degree of anisotropy.     

Stability of the boundary layer above the surface of a wave travelling over a plate

The two-dimensional problem of the stability of the flow of an incompressible fluid over a rigid surface perturbed by a wave travelling in the propagation direction of the flow is discussed in the linear approximation. The problem is solved in the coordinate system at rest with respect to the travelling wave. The parameters of this wave are not eigenvalues of the corresponding linear problem of the stability. The solution is sought in the form of a series in powers of the wave amplitude with an accuracy out to the quadratic term inclusively. Calculations are made of the dependence of the neutral stability curve on the amplitude, wavelength, and phase velocity.Translated from Zhurnal Prikladnoi Mekhaniki Tekhnicheskoi Fiziki, No. 5, pp. 49–52, September–October, 1979.     

Temperature field of heat sources during fluid injection in an anisotropic inhomogeneous reservoir

The problem of the temperature field produced by sources whose position does not depend on the vertical coordinate and which are concentrated in a horizontal permeable layer surrounded by a heat-conducting medium with radial steady-state fluid flow. The problem is solved using an averagely accurate asymptotic method. Analytical expressions for the zero-order approximation and the first coefficient of the expansion. A condition is determined under which the averaged problem for the remainder term has a trivial solution.     

Hydromagnetic flow in a viscoelastic fluid due to the oscillatory stretching surface

An analysis is carried out to study the unsteady magnetohydrodynamic (MHD) two-dimensional boundary layer flow of a second grade viscoelastic fluid over an oscillatory stretching surface. The flow is induced due to an infinite elastic sheet which is stretched back and forth in its own plane. For the investigated problem, the governing equations are reduced to a non-linear partial differential equation by means of similarity transformations. This equation is solved both by a newly developed analytic technique, namely homotopy analysis method (HAM) and by a numerical method employing the finite difference scheme, in which a coordinate transformation is employed to transform the semi-infinite physical space to a bounded computational domain. The results obtained by means of both methods are then compared and show an excellent agreement. The effects of various parameters like visco-elastic parameter, the Hartman number and the relative frequency amplitude of the oscillatory sheet to the stretching rate on the velocity field are graphically illustrated and analysed. The values of wall shear stress for these parameters are also tabulated and discussed.     

Effective properties of a composite beam subjected to torsion

The problem of finding the effective characteristics of a rectilinear beam under pure torsion is considered. The problem can be reduced to determining the torsional stress function from the solution of a boundary-value problem in a cross section of the beam for a partial differential equation with variable coefficients. Two special boundary-value problems are formulated to find the effective characteristics. It is shown that the effective coefficients are reciprocal in the case of torsion of a layer with nonuniform thickness. In the two-dimensional case, the problem is solved by a finite element method. The cases of a square beam with single and multiple inclusions are discussed. The dependence of the effective characteristics on the inclusion volume fraction is analyzed.     

Non-axisymmetrical vibration of elastic circular plate on layered transversely isotropic saturated ground

The non-axisymmetrical vibration of elastic circular plate resting on a layered transversely isotropic saturated ground was studied.First,the 3-d dynamic equations in cylindrical coordinate for transversely isotropic saturated soils were transformed into a group of governing differential equations with 1-order by the technique of Fourier ex- panding with respect to azimuth,and the state equation is established by Hankel integral transform method,furthermore the transfer matrixes within layered media are derived based on the solutions of the state equation.Secondly,by the transfer matrixes,the general solutions of dynamic response for layered transversely isotropic saturated ground excited by an arbitrary harmonic force were established under the boundary conditions, drainage conditions on the surface of.ground as well as the contact conditions.Thirdly, the problem was led to a pair of dual integral equations describing the mixed boundary- value problem which can be reduced to the Fredholm integral equations of the second kind solved by numerical procedure easily.At the end of this paper,a numerical result concerning vertical and radical displacements both the surface of saturated ground and plate is evaluated.     

Compressive buckling analysis of plates in unilateral contact

The unilateral contact buckling behaviour of delaminated plates in a composite member is studied in this paper, where the two-dimensional mechanical problem is simplified to a one-dimensional mathematical model following the assumption of a buckling mode function in terms of the lateral coordinates. After the governing differential equations for the plate parts in the contact and non-contact regions have been solved, the problem is reduced to just two nonlinear algebraic equations allowing the buckling coefficient to be obtained through an iteration method. A simplified method is also deduced based on an assumed zero length of contact. The numerical results of the buckling coefficients based on both methods show good agreement.     

Non-linear,non-steady burning of a solid propellent in three coordinate systems

The non-linear partial differential equation controlling the temperature distribution in a burning solid propellent in a rectangular, cylindrical or spherical coordinate system is transformed from a fixed coordinate system to a moving Lagrangian coordinate system, and then, using a similarity variable is further transformed to a non-linear, second order, ordinary differential equation. The burning wall temperature is time dependent, and the burning rate is an explicit function of the wall surface temperature.The boundary conditions for the resulting non-linear differential equation are given, and the necessary form of the burning rate is determined. Additionally, the linear partial differential equation for this burning solid propellent problem is also treated in all three coordinate systems. A non-linear example is given and solved, in all three coordinate systems, to illustrate the method.