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.     

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.     

The general and basic boundary value problems of plane steady jet flows and the corresponding nonlinear equations

The general boundary value problem, including known plane steady jet flows of an ideal incompressible fluid, is formulated. The simplest problem retaining all the specific features of the general problem, known as the basic problem, is separated from the general problem. The solution of the basic problem is reduced to solving a non-linear integro-differential equation and also to solving nonlinear integral equations. Examples of flows whose determination is reduced to' solving the basic problem are cited.     

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.     

轴对称体扭转问题的轴对称有限元模拟解

轴对称体的轴对称问题与扭转问题一向被认为是两个互相不能模拟的问题.前者的未知量与方程多于后者,形式也不相同.本文提出一种退化模拟方法.能够把扭转问题模拟为轴对称问题的一类特殊情况来解.一般的结构分析程序都能够分析轴对称问题,但轴对称体的扭转问题通常作为三维问题处理.按本文提出的方法,可用结构分析程序的轴对称分析功能模拟扭转分析.本文还给出模拟计算的算例.计算结果表明与理论解完全一致.本文对退化模拟的材料本构关系进行了研究,建议在数值计算时以各向异性材料的轴对称问题模拟任何材料的扭转问题.当限定用各向同性材料的轴对称问题来模拟时,采用了罚系数法.     

Apollonius coordinates,theN-body problem,and continuation of periodic solutions

This paper treats theN-body problem and its relation to various restricted problems. For each solution of the Kepler problem a generalization of the pulsating coordinates used to express the Hamiltonian of the elliptic restricted three-body problem is given. These coordinates are called Apollonius coordinates. The method of symplectic scaling is used to give a precise derivation of the elliptic restricted problem showing the precise asymptotic relationship between the restricted problem and the full three-body problem. This derivation obviates the proof of the fact that a nondegenerate periodic solution of the elliptic restricted three-body problem can be continued into the full three-body problem under mild nonresonance assumptions. Also, the method of symplectic scaling is used to give a precise derivation of the elliptic Hill lunar equation showing the precise relationship between the elliptic Hill lunar equation and the full three-body problem. A similar continuation theorem is established.     

Meshless methods for the inverse problem related to the determination of elastoplastic properties from the torsional experiment

The problem of determining the elastoplastic properties of a prismatic bar from the given experimental relation between the torsional moment M and the angle of twist per unit length of the rod’s length θ is investigated as an inverse problem. The proposed method to solve the inverse problem is based on the solution of some sequences of the direct problem by applying the Levenberg-Marquardt iteration method. In the direct problem, these properties are known, and the torsional moment is calculated as a function of the angle of twist from the solution of a non-linear boundary value problem. This non-linear problem results from the Saint-Venant displacement assumption, the Ramberg–Osgood constitutive equation, and the deformation theory of plasticity for the stress–strain relation. To solve the direct problem in each iteration step, the Kansa method is used for the circular cross section of the rod, or the method of fundamental solutions (MFS) and the method of particular solutions (MPS) are used for the prismatic cross section of the rod. The non-linear torsion problem in the plastic region is solved using the Picard iteration.     

A qualitative analysis of the brachistochrone problem with dry friction and maximizing the horizontal range

The problem of maximizing the horizontal coordinate of a point moving in a vertical plane under the action of gravity and dry friction and the corresponding brachistochrone problem are considered. The optimal control problem is reduced to a boundary value problem for a system of two nonlinear differential equations. A qualitative analysis of the trajectories of this system is carried out, their typical features are found and illustrated by numerical solving of the boundary value problem. It is shown that the normal component of the support reaction should be positive when moving along the optimal curve. The optimality of the found extremals is discussed.     

On One Approach to the Solution of Stress Problems for Noncircular Hollow Cylinders

An approach is proposed to solve three-dimensional stress problems for noncircular hollow cylinders. The end conditions are such that the problem can be reduced to a two-dimensional problem. This problem is reduced to a one-dimensional problem by introducing additional functions into the resolvable system of equations. These functions are determined using discrete Fourier series. The one-dimensional problem is solved by a stable numerical method. As an example, the stress state of cylinders with an elliptic cross section is analyzed depending on their thickness and degree of ellipticity.     

Flux difference splitting for open-channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.     

Identification of the thermal and thermostressed states of a two-layer cylinder from surface displacements

The problem of identifying the law of time variation in the temperature of one boundary surface of a two-layer cylinder and its thermal and thermostressed state from the temperature and radial displacement of the other surface is formulated and solved. The inverse problem of thermoelasticity to which the problem posed is reduced is analyzed for well-posedness. The solution of the direct problem of thermoelasticity is used to numerically test the technique of solving the inverse problem __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 1, pp. 40–47, January 2008.     

Symmetry groups and similarity solutions for a free convective boundary-layer problem

The free convective boundary-layer problem due to the motion of an elastic surface into an electrically conducting fluid is studied with group-theoretical methods. The symmetry groups admitted by the corresponding boundary value problem are obtained. Particular attention is paid on the group of scaling which provides the similarity solution of the problem. Also, the admissible form of the data, in order to be conformed to the obtained symmetries, is provided. Finally, with the use of the entailed similarity solution the problem is transformed into a boundary value problem of ODEs and is solved numerically.     

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     

Optimal synthesis in the time minimization problem under phase constraints

The time minimization problem in a linear system with a scalar control under phase constraints is considered. This problem is reduced to the problem of maximum deviations. Necessary conditions of control optimality are obtained. The results are applied to the problem of uniaxial stabilization of a satellite with consideration of a gravitational torque. The stabilization is performed by a flywheel engine with a limited angular momentum.     

One-dimensional non-Darcy flow in a semi-infinite porous media: a multiphase implicit Stefan problem with phases divided by hydraulic gradients

One-dimensional non-Darcy flow in a semi-infinite porous media is investigated. We indicate that the non-Darcy relation which is usually determined from exper-imental results can always be described by a piecewise linear function,and the problem can be equivalently transformed to a multiphase implicit Stefan problem.The novel feature of this Stefan problem is that the phases of the porous media are divided by hydraulic gradients,not the excess pore water pressures.Using the similarity transformation technique,an exact solution for the situation that the external load increases in proportion to the square root of time is developed. The study on the existence and uniqueness of the solution leads to the requirement of a group of inequalities.A similar Ste-fan problem considering constant surface seepage velocity is also investigated, and the solution, which we indicate to be uniquely existent under all conditions,is established. Meanwhile,the relation between our Stefan problem and the traditional multiphase Stefan problem is demonstrated.In the end,computational examples of the solution are presented and discussed.The solution provides a useful benchmark for verifying the accuracy of general approximate algorithms of Stefan problems, and it is also attractive in the context of inverse problem analysis.     

波流与结构物相互作用的数值模拟

采用高阶边界元法建立波流与任意形状结构物相互作用的时域数学模型。在小流速假设下,将速度势进行摄动展开,边界条件分解为双物体假设下的零阶波陡稳定问题和一阶波陡下的不稳定波浪问题。物面速度势及自由水面速度势的法向导数等未知量通过求解高阶边界元积分方程得到,而积分方程的求解则通过一个数值程序来实现。在每个时间步上采用四阶Runge-Kutta法更新下一时刻自由水面波面和速度势,自由水面上采用人工阻尼层消除散射波。通过波流与直立圆柱相互作用的数值计算研究了一阶激振力和二阶慢漂力随波数的变化关系以及圆柱周围波幅的分布曲线,并与已有频域结果和时域结果对比验证了所建模型的准确性,进而应用本文模型研究了波流与实际工程结构物相互作用的问题。     

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.     

MHD Falkner-Skan flow of Maxwell fluid by rational Chebyshev collocation method

The magnetohydrodynamics （MHD） Falkner-Skan flow of the Maxwell fluid is studied. Suitable transform reduces the partial differential equation into a nonlinear three order boundary value problem over a semi-infinite interval. An efficient approach based on the rational Chebyshev collocation method is performed to find the solution to the proposed boundary value problem. The rational Chebyshev collocation method is equipped with the orthogonal rational Chebyshev function which solves the problem on the semi-infinite domain without truncating it to a finite domain. The obtained results are presented through the illustrative graphs and tables which demonstrate the affectivity, stability, and convergence of the rational Chebyshev collocation method. To check the accuracy of the obtained results, a numerical method is applied for solving the problem. The variations of various embedded parameters into the problem are examined.     

Inplane deformation of a circular inhomogeneity with imperfect interface

The plane elastic problem of a circular inhomogeneity with an imperfect interface of spring-constant-type is reduced to the solution of a Somigliana dislocation problem, when the solution for the corresponding problem with a perfect interface is known. The Burger's vector of the Somigliana dislocation is determined so that its components satisfy two interfacial conditions involving the traction components of the corresponding problem with a perfect interface. Employing complex variables, a two-phase potential solution to the Somigliana dislocation inhomogeneity problem is developed for a general form of the Burger's vector. Detailed results are reported for a uniform eigenstrain in the inhomogeneity, and for a remote uniform heat flow in the matrix. In the latter case, the inhomogeneity behaves as a void, when it begins to slide.     

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.