Free vibration analysis of a finite-length isotropic solid elliptic cylinder using exact three dimensional elasticity theory

A three-dimensional analytical model based on Navier’s displacement equation of motion is developed to describe the free vibrations of a simply supported elastic isotropic solid elliptical cylinder of finite length. Helmholtz decomposition theorem is invoked and the classical technique of separation of variables in elliptical coordinates is employed to solve the subsequent uncoupled governing differential equations as products of ordinary and modified Mathieu functions. Inclusive numerical data are provided in a systematic manner for the first several symmetric/anti-symmetric eigen-frequencies as a function of cross sectional eccentricity for selected cylinder length ratios. The incidences of “mode crossing” and “frequency veering” among different modes of the same symmetry group as well as swapping of the corresponding mode shapes in the veering region are examined. Also, selected three dimensional deformation mode shapes are illustrated. The precision of calculations is checked by performing proper convergence studies, and the correctness of results is demonstrated using a commercial finite element software as well as by comparison with the information in the accessible literature. The collection of presented data is believed to be the first meticulous attempt to determine the free vibration frequencies of finite-length simply-supported elastic elliptical cylinders in an exact manner.     

Axisymmetric deformation of cylinders made of homogeneous and fiber-reinforced elastic materials in butt-end torsion

The equations of axisymmetric deformation of a circular cylinder made of homogeneous and weakly fiber-reinforced materials in a macroscopically two-dimensional problem statement are presented. The fibers are placed in the axial, circumferential, and radial directions and along spirals in the cylinder wall. They also can lie along opposite spirals in adjacent cylindrical layers, forming an angle-ply reinforcement scheme. At a low content of fibers, it is assumed that they are in the uniaxial stress state, determined by the axial tension or constraint compression in the deformed cylinder. Based on this mathematical model, the deformation of homogeneous and fiber-reinforced elastic cylinders of various length, with free outer surfaces, in butt-end torsion is investigated. The effect of length of the cylinders on their deformation character in torsion is evaluated. The torsion of the cylinders under the conditions of sliding fit along the inner surface and under the action of pressure on the surface, when the initial and deformed configurations are practically congruent, is considered.     

Effects of functionally graded interlayers on dispersion relations of shear horizontal waves in layered piezoelectric/piezomagnetic cylinders

The effects of functionally graded interlayers on dispersion relations of shear horizontal waves in layered piezoelectric/piezomagnetic cylinders are studied. First, the basic physical quantities of elastic waves in piezoelectric cylinder are derived by assuming that the SH waves propagate along the circumferential direction steadily. Then the transfer matrices of the functional graded interlayer and outer piezomagnetic cylinder are obtained by solving the state transfer equations with spatial-varying coefficients. Furthermore, making use of the electro-magnetic surface conditions of the outer cylinder, the dispersion relations for the shear horizontal waves in layered piezoelectric/piezomagnetic cylinders are obtained and the numerical results are shown graphically. Seven kinds of functionally graded interlayers and four kinds of electro-magnetic surface conditions are considered. It is found that the functionally graded interlayers have evident influences on the dispersion relations of shear horizontal waves in layered piezoelectric/piezomagnetic cylinders. The high order modes are more sensitive to the gradient interlayers while the low order modes are more sensitive to the electro-magnetic surface conditions.     

Thermal shock problem in a homogeneous isotropic hollow cylinder with energy dissipation

In this work, we constructed the equations of generalized thermoelasticity of a homogeneous isotropic hollow cylinder. The formulation is applied in the context of the Green and Naghdi theory of types II and III. The material of the cylinder is assumed to be homogeneous isotropic both mechanically and thermally. The problem has been solved numerically using a finite-element method. Numerical results for the temperature distribution, displacement, radial stress, and hoop stress are represented graphically. Comparisons are made with the results predicted by the types II and III. The results obtained in this paper can be used to design various homogeneous thermoelastic elements under thermal load to meet special engineering requirements.     

Solution of an axisymmetric problem of free vibrations of piezoceramic hollow cylinders of finite length by the spline collocation method

An axisymmetric problem of longitudinal free vibrations of piezoceramic hollow cylinders for some types of boundary conditions on the end faces is considered. The lateral surfaces of the cylinder are covered with thin short-circuited electrodes. The polarization of piezoceramics is directed along the axis of the cylinder. Using the method of separation of variables and method of spline collocation along the length of the cylinder, we reduce the problem to a system of ordinary differential equations, which is solved by the method of step-by-step search. We present results of calculations for a PZT 4 piezoceramic cylinder.     

Buckling stability of multilayer plates and shells with interfacial structural defects under axial compression

Based on the discrete-structural theory of thin plates and shells, a variant of the equations of buckling stability, containing a parameter of critical loading, is put forward for the thin-walled elements of a layered structure with a weakened interfacial contact. It is assumed that the transverse shear and compression stresses are equal on the interfaces. Elastic slippage is allowed over the interfaces between adjacent layers. The stability equations include the components of geometrically nonlinear moment subcritical buckling conditions for the compressed thin-walled elements. The buckling of two-layer transversely isotropic plates and cylinders under axial compression is investigated numerically and experimentally. It is found that variations in the kinematic and static contact conditions on the interfaces of layered thin-walled structural members greatly affect the magnitude of critical stresses. In solving test problems, a comparative analysis of the results of stability calculations for anisotropic plates and shells is performed with account of both perfect and weakened contacts between adjacent layers. It is found that the model variant suggested adequately reflects the behavior of layered thin-walled structural elements in calculating their buckling stability. __________ Translated from Mekhanika Kompozitnykh Materialov, Vol. 43, No. 4, pp. 513–530, July–August, 2007.     

Optimizing the packing of cylinders into a rectangular container: A nonlinear approach

The container loading problem has important industrial and commercial applications. An increase in the number of items in a container leads to a decrease in cost. For this reason the related optimization problem is of economic importance. In this work, a procedure based on a nonlinear decision problem to solve the cylinder packing problem with identical diameters is presented. This formulation is based on the fact that the centers of the cylinders have to be inside the rectangular box defined by the base of the container (a radius far from the frontier) and far from each other at least one diameter. With this basic premise the procedure tries to find the maximum number of cylinder centers that satisfy these restrictions. The continuous nature of the problem is one of the reasons that motivated this study. A comparative study with other methods of the literature is presented and better results are achieved.     

Transient analysis of thermo-elastic waves in thick hollow cylinders using a stochastic hybrid numerical method,considering Gaussian mechanical properties

This article attempts to study the stochastic coupled thermo-elasticity of thick hollow cylinders subjected to thermal shock loading considering uncertainty in mechanical properties. The thermo-elastic governing equations based on Green–Naghdi theory (without energy dissipation) are stochastically solved using a hybrid numerical method (combined Galerkin finite element and Newmark finite difference methods). The mechanical properties are considered as random variables with Gaussian distribution, which are generated using Monte Carlo simulation method with various coefficients of variations (COVs). The effects of uncertainty in mechanical properties with various coefficients of variations on thermo-elastic wave propagation are studied in detail. Also, the maximum, mean and variance of temperature, displacement and stresses are illustrated across thickness of cylinder in various times.     

Calculation of the Thermoelastoplastic Nonaxisymmetric Stress-Strain State of Layered Orthotropic Shells of Revolution

The methods for determining the nonaxisymmetric thermoelastoplastic stress-strain state of layered orthotropic shells of revolution are developed. It is assumed that the layered package deforms without mutual slippage or separation of layers. The problem is solved using the geometrically nonlinear theory of shells based on the Kirchhoff-Love hypotheses. In the isotropic layers, plastic deformations may appear, whereas the orthotropic layers deform in the elastic region. It is assumed that the mechanical properties of the materials are temperature-dependent. The thermoplasticity equations are presented in a form corresponding to the method of additional deformations. The order of the system of partial differential equations obtained is reduced with the help of trigonometric series in the circumferential coordinate. The resulting systems of ordinary differential equations are solved by the Godunov technique of discrete orthogonalization. The nonaxisymmetric thermoelastoplastic stress-strain states of layered shells of revolution are considered as examples.     

Calculation of contact stresses in axisymmetric underwater explosions

The problem of an underwater electrical explosion in a long layered elastoplastic cylinder is solved. The solution is obtained numerically with use of artificial viscosity. The law for energy produced by the explosion is determined experimentally and residual stresses are determined by the establishment method. Contact stress is presented as a function of explosion energy for various cylinders.Translated from Teoreticheskaya i Prikladnaya Mekhanika, No. 20, pp. 106–109, 1989.     

Exact solutions of problems of the three-dimensional flow of ideal viscous incompressible fluids near cylindrical surfaces

Three-dimensional flows of an incompressible fluid, the parameters of which depend on two coordinates and time, are considered. The stream surfaces of such flows are cylindrical. The equations of continuity and the Navier-Stokes equations can be transformed to relations, one of which is the equation for the stream function the other is the integral of the equations relating the pressure and the stream function, and the third is a linear equation for the projection of the velocity vector onto the axis parallel to the generatrix of the cylindrical surfaces. The problems of modelling the flows are considered on the basis of the exact solutions of the Navier-Stokes equations and Euler's equations using examples. Relations for the distribution of the flow parameters in the channel created by hyperbolical cylinders are derived for the case of unsteady inviscid flow. The streamlines of these flows are situated on the side surfaces of the hyperbolical cylinders and intercept the generatrices of the cylinders at certain indirect angles. The flow around a circular cylinder and the flow of fluid inside an elliptic cylinder are considered in the case of steady inviscid flow. The streamlines on the circular cylinder are arranged transverse to the cylinder (the projection of the velocity vector onto the coordinate axis, parallel to the generatrix of the cylinder, is equal to zero). Far from the cylinder the streamlines are also situated on a cylindrical surfaces, but not transverse to the cylinder, making certain indirect angles with the generatrix. Viscous three-dimensional flows, possessing a certain symmetry, are considered. In the case of radial symmetry the streamlines are helical lines. The non-planar Couette flow between parallel moving planes is characterized by the fact that the velocity vectors, being situated in the same plane, are collinear, while the velocity vectors in parallel planes are not collinear. Relations for viscous steady three-dimensional flows, using well-known relations, obtained for the stream function of two-dimensional flows, are given.     

Numerical solution of the free‐surface viscous flow on a horizontal rotating elliptical cylinder

The numerical solution of the free‐surface fluid flow on a rotating elliptical cylinder is presented. Up to the present, research has concentrated on the circular cylinder for which steady solutions are the main interest. However, for noncircular cylinders, such as the ellipse, steady solutions are no longer possible, but there will be periodic solutions in which the solution is repeated after one full revolution of the cylinder. It is this new aspect that makes the investigation of noncircular cylinders novel. Here we consider both the time‐dependent and periodic solutions for zero Reynolds number fluid flow. The numerical solution is expedited by first mapping the fluid film domain onto a rectangle such that the position of the free‐surface is determined as part of the solution. For the time‐dependent case a simple time‐marching method of lines approach is adopted. For the periodic solution the discretised nonlinear equations have to be solved simultaneously over a time period. The resulting large system of equations is solved using Newton's method in which the form of the Jacobian enables a straightforward decomposition to be implemented, which makes matrix inversion manageable. In the periodic case all derivatives have been approximated pseudospectrally with the time derivative approximated by a differentiation matrix which has been specially derived so that the weight of fluid is algebraically conserved. Of interest is the solution for which the weight of fluid is at its maximum possible value, and this has been obtained by increasing the weight until a consistency break‐down occurs. Time‐dependent solutions do not produce the periodic solution after a long time‐scale but have protuberances which are constantly appearing and disappearing. Periodic solutions exhibit spectral accuracy solutions and maximum supportable weight solutions have been obtained for ranges of eccentricity and angular velocity. The maximum weights are less than and approximately proportional to those obtained for the circular case. The shapes of maximum weight solutions is distinctly different from sub‐maximum weight solutions. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Numerical solution of the problem of the stress-strain state in hollow cylinders using spline approximations

The three-dimensional theory of elasticity is used for a study of the stress-strain state in a hollow cylinder with varying stiffness. The corresponding problem is solved by a method that is partly analytical and partly numerical in nature: Spline approximations and collocation are used to reduce the partial differential equations of elasticity to a boundary-value problem for a system of ordinary differential equations of higher order for the radial coordinate, which is then solved using the method of stable discrete orthogonalization. Results for an inhomogeneous cylinder for various types of stiffness are presented.     

通过有漏孔管道时的层流解析解

研究了通过有漏孔管道时的层流,并解析地求解了动量方程和能量守恒方程．由Hagen-Poiseuille的速度分布,同时定义轴向和径向速度分量的未知函数,得到了压力和质量输运方程,并根据不同的参数,画出其分布图．结果表明,管道中的轴向速度、径向速度、质量输运参数和压力,随着流体沿管道的流动而减小．     

Computer simulation of a twisted nanotube buckling

We develop procedures for solving the problems of dynamic nanostructure deformation and buckling numerically. The procedures are based on discretization with respect to time of the nonlinear equations of molecular mechanics whose matrices and vectors are determined using the Morse potential for the central forces of interaction between atoms and fictitious truss elements accounting for the variations of the angle between atomic bonds. To determine the critical values of deformation parameters and the shapes of buckling nanostructures we use a stability loss criterion for solutions to nonlinear ordinary differential equations on a finite time interval. We implemented our procedures in the PIONER code, using which we solve the problem of a twisted nanotube buckling in the conditions of a quasistatic deformation. To determine the postcritical equilibrium modes we solve the same problem in a dynamic formulation. We show that the modes of equilibrium configurations of the nanotube in the initial postcritical deformation correspond to a buckling mode obtained both at the bifurcation point of quasistatic solutions and at the quasibifurcation point of dynamic solutions.     

Guided ultrasonic waves in composite cylinders

The ultrasonic nondestructive evaluation of composite cylinders requires a thorough understanding of the propagation of waves in these materials. In this paper, the propagation of flexural and longitudinal guided waves in fiber-reinforced composite (FRC) rods with transversely isotropic symmetry properties is studied. The frequency equations obtained for free cylinders and the effect of increased fiber volume fraction (increased anisotropy) on the dispersion characteristics of the rod are explored. The numerical results reveal a number of previously unnoticed characteristics of dispersion curves for composite cylinders. The mode shapes of longitudinal waves propagating in glass/epoxy cylinders are also plotted. These plots can be used to choose an appropriate strategy for inspecting composite cylinders by ultrasonic nondestructive evaluation techniques. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 43, No. 3, pp. 411–426, May–June, 2007.     

Modelling two-dimensional flow past arbitrary cylindrical bodies using boundary element formulations

This paper presents an efficient method of solving Queen's linearized equations for steady plane flow of an incompressible, viscous Newtonian fluid past a cylindrical body of arbitrary cross-section. The numerical solution technique is the well known direct boundary element method. Use of a fundamental solution of Oseen's equations, the ‘Oseenlet’, allows the problem to be reduced to boundary integrals and numerical solution then only requires boundary discretization. The formulation and solution method are validated by computing the net forces acting on a single circular cylinder, two equal but separated circular cylinders and a single elliptic cylinder, and comparing these with other published results. A boundary element representation of the full Navier-Stokes equations is also used to evaluate the drag acting on a single circular cylinder by matching with the numerical Oseen solution in the far field.     

Solution of problems of free torsional vibrations of thick-walled orthotropic inhomogeneous cylinders

With the use of the 3D theory of elasticity, we investigate the problem of free torsional vibrations of an anisotropic hollow cylinder with different boundary conditions at its end faces. We have proposed a numerical-analytic approach for the solution of this problem. The original partial differential equations of the theory of elasticity with the use of spline approximation and collocation are reduced to an eigenvalue problem for a system of ordinary differential equations of high order in the radial coordinate. This system is solved by the stable numerical method of discrete orthogonalization together with the method of step-by-step search. We also present numerical results for the case of orthotropic and inhomogeneous material of the cylinder for some kinds of boundary conditions.     

Second-approximation shear theory for shallow layered shells and plates

Analysis of a second-approximation refined shear model for shallow layered composite shells and plates with a substantially inhomogeneous structure over the thickness is presented. The tangential displacements and corresponding normal stresses are expressed in the form of a polynomial of the fith degree in the transverse coordinate and contain squared rigidity characteristics. In this way, the accuracy of results and practical coincidence with the 3D solutions is ensured. Based on the refined model, a theory of shallow layered shells is developed. A system of resolving equations of sixteenth power together with appropriate boundary conditions was obtained and solved analytically. It is shown that the area of application of the formed model is extended as compared with the model of the first approximation. The model proposed allows us to examine the stress-strain state of layered composite structures of substantially different thickness and physical-mechanical characteristics of the layers, including the possibility of simulating relatively large shear deformations of rigid layers separated by a low-modulus thin interlayer pliable to transverse shear.Presented at the 10th International Conference on the Mechanics of Composite Materials (Riga, April 20–23, 1998).Ukrainian Transport University, Kiev, Ukraine. Translated from Mekhanika Kompozitnykh Materialov, Vol. 34, No. 3, pp. 363–370, May–June, 1998.     

The stress state of the finite cylinder with the circular cracks

The axysimmetrical torsion problem for a finite cylinder with an arbitrary quantity of the cylindrical layers is solved. The cylinder is weaked by the parallel circular cracks in the first internal layer. The stated boundary valued problem problem is reduced to the system of the integro-duferential equations solving with the help of the orthogonal polinomials method. The stress intensity factors (SIF) are obtained. The dependences of SIF values from the cracks' sizes, their location, and ratio of the layers' shear moduluses are concretized for the case of the two layers. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)