Analogy solutions of axisymmetrical twist problems by axisymmetrical finite element program

It used to be considered that an axisymmetrical problem and a twist problem of an axisymmetrical body cannot be simulated by each other, because the number of unknown variables in an axisymmetrical problem is greater than that in a twist problem, and the governing equations are not the same. This paper proposes a degenerated analogy method, by which the twist problems of axisymmetrical bodies can be simulated by axisymmetrical problems with finite element programs.An ordinary structural analysis method can be used to analyze an axisymmetrical problem, but a twist problem of axisymmetrical bodies is treated as a 3-dimensional problem usually. According to the method proposed in this paper, the analysis of a twist problem can be simulated by the analysis of an axisymmetrical body with a structural analysis problem. The example of analysis computation is also given. Thecomputed result is in agreement with the theoretical result.In this paper, the constitutive relation of the degenerated analogy problem is given.The authors suggest that a twist problem of a body made of any materials is simulated by an axisymmetrical problem of a body made of orthotropic material. If you have to use some program for the axisymmetrical problem to be limited to isotropic materials the penalty coefficient method can be used to solve the problem.     

Approximation of thermoelasticity contact problem with nonmonotone friction

The paper presents the formulation and approximation of a static thermoelasticity problem that describes bilateral frictional contact between a deformable body and a rigid foundation. The friction is in the form of a nonmonotone and multivalued law. The coupling effect of the problem is neglected. Therefore, the thermic part of the problem is considered independently on the elasticity problem. For the displacement vector, we formulate one substationary problem for a non-convex, locally Lipschitz continuous functional representing the total potential energy of the body. All problems formulated in the paper are approximated with the finite element method.     

Homogenization of the equations of the Cosserat theory of elasticity of inhomogeneous bodies

The paper deals with the homogenization of a boundary value problem for an inhomogeneous body with Cosserat properties, which is referred to as the original problem. The homogenization process is understood as a method for representing the solution of the original problem in terms of the solution of precisely the same problem for a body with homogeneous properties. The problem for a body with homogeneous properties is called the accompanying problem, and the body itself, the accompanying homogeneous body. As a rule, a constructive homogenization procedure includes the following three stages: at the first stage, the properties of the inhomogeneous body are used to find the properties of the accompanying homogeneous body (efficient properties); at the second stage, the boundary value problem is solved for the accompanying body; at the third stage, the solution of the accompanying problem is used to find the solution of the original problem. This approach was implemented in mechanics of composite materials constructed of numerous representative elements. A significant contribution to the development of mechanics of composites is due to Rabotnov [1–3] and his students. Recently, the homogenization method has been widely used to solve problems for composites of regular structure by expanding the solution of the original problem in a power series in a small geometric parameter equal to the ratio of the characteristic dimension of the periodicity cell to the characteristic dimension of the entire body. The papers by Bakhvalov [4–6] and Pobedrya [7] were the first in the field. At present, there are numerous monographs partially or completely dealing with the method of a small geometric parameter [8–14]. Isolated problems for inhomogeneous bodies with nonperiodic dependence of their properties on the coordinates were considered by many authors. Most of such papers published before 1973 are collected in two vast bibliographic indices [15, 16]. General methods were considered, and many specific problems of the theory of elasticity of continuously inhomogeneous bodies were solved in Lomakin’s papers and his monograph [17]. The theory of torsion of inhomogeneous anisotropic rods was considered in [18]. In 1991, in his Doctoral dissertation, one of the authors of this paper proposed a version of the homogenization method based on an integral formula representing the solution of the original static problem of inhomogeneous elasticity via the solution of the accompanying problem [19, 20]. An integral formula for the dynamic problem of elasticity was published somewhat later [21]. This integral formula was used to develop a constructive method for the homogenization of the dynamic problem of inhomogeneous elasticity, which can be used in the case of both periodic and nonperiodic inhomogeneity of the properties [22]. The integral formula in the case of the Cosserat theory of elasticity was published in [23]. The present paper briefly presents constructive methods for homogenizing the problems of the Cosserat theory of elasticity based on the integral formula.     

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     

Geometric shapes of the interface surface of bicomponent flows between two concentric rotating cylinders

In this paper, the shape problem of interface of bicomponent flows between two concentric rotating cylinders is investigated. With tensor analysis, the problem is reduced to an energy functional isoperimetric problem when neglecting the effects of the dissipative energy caused by viscosity. We derive the associated Euler-Lagrangian equation, which is a nonlinear elliptic boundary value problem of the second order. Moreover, by considering the effects of the dissipative energy, we propose another total energy functional to characterize the geometric shape of the interface, and obtain the corresponding Euler-Lagrangian equation, which is also a nonlinear elliptic boundary value problem of the second order. Thus, the problem of the geometric shape is converted into a nonlinear boundary value problem of the second order in both cases.     

Nonaxisymmetric vibrations of radially polarized hollow cylinders made of functionally gradient piezoelectric materials

The nonaxisymmetric problem of natural vibrations of radially polarized hollow cylinders made of functionally gradient piezoelectric materials is solved. The properties of the material change continuously along a radial coordinate according to an exponential law. The lateral surfaces of the cylinder are free of external tractions and short circuited by electrodes. After separation of variables and representation of the components of the displacement vector in the form of standing circumferential waves, the initially three-dimensional problem is reduced to a two-dimensional partial differential equation problem. By using the method of spline-collocations with respect to a longitudinal coordinate, this two-dimensional problem is reduced to a boundary-value problem for the eigenvalues expressed in terms of ordinary differential equations. This problem is solved by the stable discrete-orthogonalization technique in combination with a step-by-step search method with respect to the radial coordinate. Results were obtained numerically and subsequently analyzed in this paper.     

On the existence of non-polarized azimuthal-symmetric electromagnetic waves in circular dielectric waveguide filled with nonlinear isotropic homogeneous medium

The propagation of monochromatic nonlinear symmetric hybrid waves in a cylindrical nonlinear dielectric waveguide is considered. The physical problem is reduced to solving a transmission eigenvalue problem for a system of ordinary differential equations. Spectral parameters of the problem are propagation constants of the waveguide. The problem is reduced to the new type of nonlinear eigenvalue problem. The analytical method of solving this problem is presented. New propagation regime is found.     

Dynamics of vapor bubbles

The nonlinear problem of the thermal, mass, and dynamic interaction of a single vapor bubble with the surrounding liquid is discussed. This problem has ramifications in research on flows of vapor-liquid mixtures with a bubble-matrix structure, in particular, the propagation of shock waves in such media. Results are given from a numerical solution of the problem of the radial motion imparted to a bubble by a sudden change of pressure in the liquid; this problem corresponds, in particular, to the behavior of bubbles behind a shock front when the latter enters a bubble curtain.     

A numerical approximation of the thermal coupling of fluids and solids

In this article we analyze the problem of the thermal coupling of fluids and solids through a common interface. We state the global thermal problem in the whole domain, including the fluid part and the solid part. This global thermal problem presents discontinuous physical properties that depend on the solution of auxiliary problems on each part of the domain (a fluid flow problem and a solid state problem). We present a domain decomposition strategy to iteratively solve problems posed in both subdomains and discuss some implementation aspects of the algorithm. This domain decomposition framework is also used to revisit the use of wall function approaches used in this context. Copyright © 2008 John Wiley & Sons, Ltd.     

Symmetric shapes of contact of a jet with the surface of a heavy liquid

The plane linearized problem of oblique impingement of a weightless jet of an ideal incompressible fluid on the surface of a heavy fluid is considered. Flows are sought with symmetric forms of the contact region. Mathematically we arrive at the problem of the eigenvalues and eigenfunctions of an integral equation; solving this problem we obtain various contact forms. The fundamental result for the infinitesimally thin jet of finite intensity is derived by passing to the limit, yielding a result analagous with the forms of free vibrations of a string. Some results are presented for the problem under consideration in the nonlinear formulation.The two-dimensional problem on (vertical) impingement of a jet on a liquid was solved by Olmstead and Raynor [1]. Some results for oblique impingement of a sufficiently thin, slightly curved jet are presented by Frolov [2], Information on other studies, primarily experimental, is presented in [3].This problem is related to the model of a jet curtain of an air-cushion vehicle; in this regard we note the study of Stepanov [4] in which, in particular, a result is obtained for an infinitesimally thin jet curtain.     

New Statement of the elasticity problem for a layer

The problem on the equilibrium of an inhomogeneous anisotropic elastic layer is considered. The classical statement of the problem in displacements consists of three partial differential equations with variable coefficients for the three displacements and of three boundary conditions posed at each point of the boundary surface. Sometimes, instead of the statement in displacements, it is convenient to use the classical statement of the problem in stresses [1] or the new statement of the problem in stresses proposed by B. E. Pobedrya [2]. In the case of the problem in stresses, it is necessary to find six components of the stress tensor, which are functions of three coordinates. The choice of the statement of the problem depends on the researcher and, of course, on the specific problem. The fact that there are several statements of the problem makes for a wider choice of the method for solving the problem. In the present paper, for a layer with plane boundary surfaces, we propose a new statement of the problem, which, in contrast to the other two statements indicated above, can be called a mixed statement. The problem for a layer in the new statement consists of a system of three partial differential equations for the three components of the displacement vector of the midplane points. The system is coupled with three integro-differential equations for the three longitudinal components of the stress tensor. Thus, in the new statement, just as in the other statements in stresses, one should find six functions. In the new statement, three of these functions (the displacements of the midplane points) are functions of two coordinates, and the other three functions (the longitudinal components of the stress tensor) are functions of three coordinates. It is shown that all equations in the new statement are the Euler equations for the Reissner functional with additional constraints. After the problem is solved in the new statement, three components of the displacement vector and three transverse components of the stress tensor are determined at each point of the layer. The new statement of the problem can be used to construct various engineering theories of plates made of composite materials.     

On one class of differential equations of fractional order

We consider the Darboux problem for a differential equation of fractional order that contains a regularized mixed derivative. Sufficient conditions for the existence and uniqueness of a solution of this problem are obtained in the class of continuous functions. We also propose a method for finding an approximate solution of this problem and prove the convergence of this method.     

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.     

Control of homogeneous shear flow of multimode Maxwell fluids

We consider a homogeneous parallel shear flow of a multimode Maxwell fluid. This problem results in a set of ordinary differential equations for the stresses. In this system, we view the shear rate as a control and consider the problem of steering the system to a given state of stress. The objective is to steer the system from given initial stresses to a final state of stress, allowing the shear rate to vary in an arbitrary fashion. We show that this problem is related to a calculus of variations problem. For the case of two modes, we obtain a characterization of the set of achievable streses.     

Solution of multiple crack problem in a finite plate using coupled integral equations

This paper studies a numerical solution of multiple crack problem in a finite plate using coupled integral equations. After using the principle of superposition, the multiple crack problem in a finite plate can be converted into two problems: (a) the multiple crack problem in an infinite plate and (b) a usual boundary value problem for the finite plate. For the former problem, the Fredholm integral equation is used. For the latter problem, a BIE based on complex variable is suggested in which a Cauchy singular kernel exists. For the proposed BIE, after using the inverse matrix technique, the dependence of the traction at a domain point from the boundary tractions is formulated indirectly. This is a particular advantage of the present study. Several numerical examples are provided and the computed results for stress intensity factor and T-stress at crack tips are given.     

Calculation of neutral monotonic instability curves in the problem of concentration convection in a mixture with an infinite number of components

The Oberbeck-Boussinesq approximation for concentration convection in a mixture with an infinite number of components is constructed. The features of the formulation of the problem are described in detail. The large-parameter asymptotics are constructed for the linear problem of hydrodynamic stability. The problem is substantially simplified and equations not previously encountered in hydrodynamic stability theory are obtained. In the case of the non self-adjoint problem the asymptotics of the eigenvalues and eigenfunctions are obtained. Numerical results which, in particular, show that the spectrum of the boundary value problem is not connected are presented. The critical values obtained make it possible to solve the important practical problem of improving the process of mixture separation by the isoelectric focusing method.Rostov-on-Don. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, pp. 11–20, September–October, 1995.     

Stress analysis near the tips of a transverse crack in an elastic semi-strip

The plane elastic problem for a semi-strip with a transverse crack is investigated. The initial problem is reduced to a one-dimensional continuous problem by use of an integral transformation method with a generalized scheme. The one-dimensional problem is first formulated as a vector boundary problem, and then reduced to a system of three singular integral equations(SIEs). The system is solved by use of an orthogonal polynomial method and a special generalized method. The contribution of this work is the consideration of kernel fixed singularities in solving the system. The crack length and its location relative to the semi-strip's lateral sides are investigated to simplify the problem's statement. This simplification reduces the initial problem to a system of two SIEs.     

Optimal Microstructures Drag Reducing Mechanism of Riblets

We consider an optimal shape design problem of periodically distributed three-dimensional microstructures on surfaces of swimming bodies in order to reduce their drag. Our model is restricted to the flow in the viscous sublayer of the boundary layer of a turbulent flow. The costs for the optimization problem are very high because the three-dimensional flow equations have to be solved several times. We avoid this problem by approximations: the microscopic optimization problem is reduced applying homogenization. Considering a special geometry (riblets) the resulting so-called macroscopic optimization problem can be additionally reduced to a two-dimensional problem. We analyze the drag reducing mechanism of riblets which are believed to be optimal structures. Therefore we perform direct simulations on the total rough channel for different shapes of microstructures: riblets and fully three-dimensional structures.     

Transient piezothemoelasticity of a hexagonal plate of class 6 mm

Summary Piezoelastic materials entail piezoelectric polarization that is directly proportional to the applied strain. Therefore, piezoelastic materials have recently attracted considerable attention because of their potential use in intelligent structural systems. In this paper, we treat a transient problem of piezothermoelasticity in a hexagonal plate of class 6 mm subjected to transient temperature change. We analyze the transient problem by use of the potential function method, which can separate the coupling between the thermoelastic and electric fields. Numerical calculations for the transient problem are carried out for a cadmium selenide solid, and the results are presented graphically in comparison with those for a steady problem.     

Hertz接触问题的解的唯一性条件

讨论了接触面为圆面的Ｈｅｒｔｚ接触问题。若压力分布是轴对称的，则该接触问题的解必是唯一的。且在上述条件下，该接触问题的积分方程可化为两个推广的Ａｂｅｌ积分方程组，此方程组的解便给出此接触问题的解。