Solving the three-dimensional equations of the linear theory of elasticity

The three-dimensional Lamé equations are solved using Cartesian and curvilinear orthogonal coordinates. It is proved that the solution includes only three independent harmonic functions. The general solution of equations of elasticity for stresses is found. The stress tensor is expressed in both coordinate systems in terms of three harmonic functions. The general solution of the problem of elasticity in cylindrical coordinates is presented as an example. The three-dimensional stress–strain state of an elastic cylinder subjected, on the lateral surface, to arbitrary forces represented by a series of eigenfunctions is determined. An axisymmetric problem for a finite cylinder is solved numerically     

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.     

Mixed boundary value problem of elasticity for a quarter space

We consider the static elasticity problem for a quarter space with zero displacements on one of its surfaces and with given stresses on the other. The method for solving this problem is based on the use of newunknown functions in the formof a linear combination of the desired displacements, which reduces the system of three Lamé equations to two equations to be solved simultaneously and one equation to be solved separately. The exact solution of this problem was obtained earlier by the same method [1]. But it was shown in [2] that such a solution is exact only under certain restrictions on the given functions. In the present paper, the solution of this problem is constructed without restrictions on the given functions, which necessitates solving a one-dimensional integro-differential equation; this can be done approximately by the orthogonal polynomial method. We present numerical results obtained on the basis of our solution.     

A Mode- II Griffith Crack in Decagonal Quasicrystals

By using the method of stress functions, the problem of mode-II Griffith crack in decagonal quasicrystals was solved. First, the crack problem of two-dimensional quasicrystals was decomposed into a plane strain state problem superposed on anti-plane state problem and secondly, by introducing stress functions, the18 basic elasticity equations on coupling phonon-phason field of decagonal quasicrystals were reduced to a single higher-order partial differential equations. The solution of this equation under mixed boundary conditions of mode-II Griffith crack was obtained in terms of Fourier transform and dual integral equations methods. All components of stresses and displacements can be expressed by elemental functions and the stress intensity factor and the strain energy release rate were determined. Biography: GUO Yu-cui (1962-), Associate professor, Doctor     

Group foliation of the Lamé equations of the classical dynamical theory of elasticity

We perform the group foliation of the system of Lamé equations of the classical dynamical theory of elasticity for an infinite subgroup contained in a normal divisor of the main group. The resolving system of this foliation includes the following two classical systems of mathematical physics: the system of equations of vortex-free acoustics and the system of Maxwell equations, which allows one to use wider groups to obtain exact solutions of the Lamé equations. We obtain a first-order conformal-invariant system, which describes shear waves in a three-dimensional elastic medium. We also give examples of partially invariant solutions.     

General solutions of weakened equations in terms of stresses in the theory of elasticity

The general solutions of some weakened systems of equations expressed in terms of stresses in the isotropic theory of elasticity are analyzed. These systems are not equivalent to the classical one and involve the equilibrium equations and only three of the six equations of compatibility (either diagonal or off-diagonal ones). In the framework of elasticity theory, an equivalence of the formulations of quasistatic boundary value problems based on such systems and expressed in terms of stresses is discussed.     

Recursive-operator method in vibration problems for rod systems

Using linear differential equations with constant coefficients describing one-dimensional dynamical processes as an example, we show that the solutions of these equations and systems are related to the solution of the corresponding numerical recursion relations and one does not have to compute the roots of the corresponding characteristic equations. The arbitrary functions occurring in the general solution of the homogeneous equations are determined by the initial and boundary conditions or are chosen from various classes of analytic functions. The solutions of the inhomogeneous equations are constructed in the form of integro-differential series acting on the right-hand side of the equation, and the coefficients of the series are determined from the same recursion relations. The convergence of formal solutions as series of a more general recursive-operator construction was proved in [1]. In the special case where the solutions of the equation can be represented in separated variables, the power series can be effectively summed, i.e., expressed in terms of elementary functions, and coincide with the known solutions. In this case, to determine the natural vibration frequencies, one obtains algebraic rather than transcendental equations, which permits exactly determining the imaginary and complex roots of these equations without using the graphic method [2, pp. 448–449]. The correctness of the obtained formulas (differentiation formulas, explicit expressions for the series coefficients, etc.) can be verified directly by appropriate substitutions; therefore, we do not prove them here.     

Determination of a general solution of three-dimensional Lamé equations of elasticity theory

We integrate the Lamé equation and find new solutions in the case of three-dimensional elasticity theory, which are expressed in terms of harmonic functions. We prove that the solution obtained involves only three independent functions. In a curvilinear orthogonal coordinate system, a general solution of the Lamé equation is expressed in terms of three harmonic functions. __________ Translated from Neliniini Kolyvannya, Vol. 9, No. 1, pp. 109–116, January–March, 2006.     

Variational formulation of a simplified strain gradient elasticity theory and its application to a pressurized thick-walled cylinder problem

A detailed variational formulation is provided for a simplified strain gradient elasticity theory by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the complete boundary conditions of the theory for the first time. To supplement the stress-based formulation, the coordinate-invariant displacement form of the simplified strain gradient elasticity theory is also derived anew. In view of the lack of a consistent and complete formulation, derivation details are included for the tutorial purpose. It is shown that both the stress and displacement forms of the simplified strain gradient elasticity theory obtained reduce to their counterparts in classical elasticity when the strain gradient effect (a measure of the underlying material microstructure) is not considered. As a direct application of the newly obtained displacement form of the theory, the problem of a pressurized thick-walled cylinder is analytically solved. The solution contains a material length scale parameter and can account for microstructural effects, which is qualitatively different from Lamé’s solution in classical elasticity. In the absence of the strain gradient effect, this strain gradient elasticity solution reduces to Lamé’s solution. The numerical results reveal that microstructural effects can be large and Lamé’s solution may not be accurate for materials exhibiting significant microstructure dependence.     

The stress field near the front of an arbitrarily shaped crack in a three-dimensional elastic body

The problem stated in the title is investigated with special emphasis on the first three terms of the stress expansion, proportional to r ^-1/2, r ⁰=1 and r ^1/2 respectively, where r denotes the distance to the crack front. The particular case of a plane crack with a straight front and of stresses independent of the distance along the latter is studied first. It is shown that the classical plane strain and antiplane solutions must be supplemented by a few additional particular solutions to obtain the full stress expansion. The general case is then considered. The stress expansion is studied by writing the field equations (equilibrium, strain compatibility and boundary conditions) in a system of suitable curvilinear coordinates. It is shown that the number of independent constants in the stress expansion is the same as in the particular case considered previously but that the curvatures of the crack and its front and the non-uniformity of the stresses along the latter induce the appearance of corrective terms in this expansion.     

Automated Stress Separation Along Stress Trajectories

A procedure for the separation of principal stresses in automated photoelasticity is presented. It is based on the integration of indefinite equations of equilibrium along stress trajectories, also known as Lamè–Maxwell equations. A new algorithm for precise and reliable stress trajectory calculation, which is an essential feature of the procedure, has also been developed. Automated stress separation is carried out along stress trajectories starting from free boundaries. Experimental tests were performed on a disc in diametral compression and on a ring with internally applied pressure. Full-field principal stress values were obtained and results were compared with those from the theory of elasticity and with those obtained from the classical shear difference method. It was shown that the proposed method is more accurate and less affected by the presence of residual stresses or experimental errors at the boundaries than the shear difference method. In addition, the method requires little human interaction and is therefore well-suited for automated photoelasticity.     

A refined theory of elastic thick plates for extensional deformation

Based on elasticity theory, various two-dimensional (2D) equations and solutions for extensional deformation have been deduced systematically and directly from the three-dimensional (3D) theory of thick rectangular plates by using the Papkovich–Neuber solution and the Lur’e method without ad hoc assumptions. These equations and solutions can be used to construct a refined theory of thick plates for extensional deformation. It is shown that the displacements and stresses of the plate can be represented by the displacements and transverse normal strain of the midplane. In the case of homogeneous boundary conditions, the exact solutions for the plate are derived, and the exact equations consist of three governing differential equations: the biharmonic equation, the shear equation, and the transcendental equation. With the present theory a solution of these can satisfy all the fundamental equations of 3D elasticity. Moreover, the refined theory of thick plate for bending deformation constructed by Cheng is improved, and some physical or mathematical explanations and proof are provided to support our justification. It is important to note that the refined theory is consistent with the decomposition theorem by Gregory. In the case of nonhomogeneous boundary conditions, the approximate governing differential equations and solutions for the plate are accurate up to the second-order terms with respect to plate thickness. The correctness of the stress assumptions in the classic plane-stress problems is revised. In an example it is shown that the exact or accurate solutions may be obtained by applying the refined theory deduced herein.     

Interaction between elliptical and ellipsoidal inclusions under bending stress fields

Summary  This paper deals with interaction problems of elliptical and ellipsoidal inclusions under bending, using singular integral equations of the body force method. The problems are formulated as a system of singular integral equations with Cauchy-type or logarithmic-type singularities, where unknown functions are densities of body forces distributed in the x,y and r,θ,z directions in infinite bodies having the same elastic constants as those of the matrix and inclusions. In order to satisfy the boundary conditions along the elliptical and the ellipsoidal boundaries, the unknown functions are approximated by a linear combination of fundamental density functions and polynomials. The present method is found to yield the exact solutions for a single elliptical or spherical inclusion under a bending stress field. It yields rapidly converging numerical results for interface stresses in the interaction of inclusions. Received 9 September 1999; accepted for publication 15 January 2000     

变边界粘弹性轴对称问题的复变函数法

对边界几何形状、位置随时间变化的变边界结构,给出了用复变函数求解粘弹问题的解析方法。文中用拉普拉斯变换结合平面弹性复变方法,对内外边界变化时粘弹性轴对称问题进行求解。引入两个与时间、空间相关的解析函数,给出了变边界情况下应力、位移以及边界条件与解析函数的关系。当解析函数形式部分确定,则可用边界条件求解其中与时间相关的待定函数。求解待定函数的方程一般情况下为一系列积分方程,特殊情况可求得解析解。对轴对称问题中应力边值问题、位移边值问题以及混合边值问题,分别利用边界条件求得相关系数,从而得到了应力与位移的解析表达。当取Boltzmann粘弹模型时,进行不同边值问题的分析。分析显示,应力、位移的形态与大小均与边界变化过程相关,与固定边界粘弹性问题有较大不同。本文解答可用于粘弹性轴对称问题内外边界任意变化及各种边值问题的力学分析。此外,该法可进一步进行荷载非对称、复杂孔型变边界问题的求解。     

The derivation of a thick and thin plate formulation without ad hoc assumptions

For the plate formulation considered in this paper, appropriate three-dimensional elasticity solution representations for isotropic materials are constructed. No a priori assumptions for stress or displacement distributions over the thickness of the plate are made. The strategy used in the derivation is to separate functions of the thickness variable z from functions of the coordinates x and y lying in the midplane of the plate. Real and complex 3-dimensional elasticity solution representations are used to obtain three types of functions of the coordinates x, y and the corresponding differential equations. The separation of the functions of the thickness coordinate can be done by separately considering homogeneous and nonhomogeneous boundary conditions on the upper and lower faces of the plate. One type of the plate solutions derived involves polynomials of the thickness coordinate z. The other two solution forms contain trigonometric and hyperbolic functions of z, respectively. Both bending and stretching (or in-plane) solutions are included in the derivation.     

An elasticity solution of a nonhomogeneous half-plane problem

ＡＮＥＬＡＳＴＩＣＩＴＹＳＯＬＵＴＩＯＮＯＦＡＮＯＮＨＯＭＯＧＥＮＥＯＵＳＨＡＬＦ－ＰＬＡＮＥＰＲＯＢＬＥＭＷｏＧｕｏ－ｗｅｉ（沃国纬）（ＳｈａｎｇｈａｉＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ．Ｓｈａｎｇｈａｉ）（ＲｅｃｅｉｖｅｄＪａｎ．５．１９９４：...     

A finite element implementation of the stress gradient theory

In this contribution, a finite element implementation of the stress gradient theory is proposed. The implementation relies on a reformulation of the governing set of partial differential equations in terms of one primary tensor-valued field variable of third order, the so-called generalised displacement field. Whereas the volumetric part of the generalised displacement field is closely related to the classic displacement field, the deviatoric part can be interpreted in terms of micro-displacements. The associated weak formulation moreover stipulates boundary conditions in terms of the normal projection of the generalised displacement field or of the (complete) stress tensor. A detailed study of representative boundary value problems of stress gradient elasticity shows the applicability of the proposed formulation. In particular, the finite element implementation is validated based on the analytical solutions for a cylindrical bar under tension and torsion derived by means of Bessel functions. In both tension and torsion cases, a smaller is softer size effect is evidenced in striking contrast to the corresponding strain gradient elasticity solutions.

     

Numerical solution of singular integral equations in stress concentration problems under longitudinal shear loading

Summary The paper deals with numerical solutions of singular integral equations in stress concentration problems for longitudinal shear loading. The body force method is used to formulate the problem as a system of singular integral equations with Cauchy-type singularities, where unknown functions are densities of body forces distributed in the longitudinal direction of an infinite body. First, four kinds of fundamental density functions are introduced to satisfy completely the boundary conditions for an elliptical boundary in the range 0≤φ _k ≤2π. To explain the idea of the fundamental densities, four kinds of equivalent auxiliary body force densities are defined in the range 0≤φ _k ≤π/2, and necessary conditions that the densities must satisfy are described. Then, four kinds of fundamental density functions are explained as sample functions to satisfy the necessary conditions. Next, the unknown functions of the body force densities are approximated by a linear combination of the fundamental density functions and weight functions, which are unknown. Calculations are carried out for several arrangements of elliptical holes. It is found that the present method yields rapidly converging numerical results. The body force densities and stress distributions along the boundaries are shown in figures to demonstrate the accuracy of the present solutions. Received 26 May 1998; accepted for publication 27 November 1998     

On the nature of the relaxation time,the Maxwell–Cattaneo and Fourier law in the thermodynamics of a continuous medium,and the scale effects in thermal conductivity

We present the theory of space–time elasticity and demonstrate that it is the extended reversible thermodynamics and gives the coupled model of thermoelasticity and heat conductivity and involves traditional thermoelasticity. We formulate the generally covariant variational model’s dynamic thermoelasticity and heat conductivity in which the basic kinematic and static variables are unified tensor objects (subject, matter). Variation statement defines the whole set of the initial-boundary problems for the 4D vector governing equation (Euler equation), the spatial projections of which define motion equations and the time projection gives the heat conductivity equation. We show that space–time elasticity directly implies the Fourier and the Maxwell–Cattaneo laws of heat conduction. However, space–time elasticity is richer than classical thermoelasticity, and it advocates its own equations of motion for coupled thermoelasticity. Moreover, we establish that the Maxwell–Cattaneo law and Fourier law can be defined for the reversible processes as compatibility equations without introducing dissipation. We argue that the present framework of space–time elasticity should prove adequate to describe the thermoelastic phenomena at low temperatures for interpreting the results of molecular simulations of heat conduction in solids and for the optimal heat and stress management in the microelectronic components and the thermoelectric devices.