A refined state space formalism for piezothermoelasticity

The state space formalism for piezothermoelasticity [Tarn, J.Q., 2002c. A state space formalism for piezothermoelasticity. International Journal of Solids and Structures 39, 5173–5184.] is refined by introducing the generalized displacement vector and generalized stress vectors as the fundamental variables in which appropriate electrical variables are included. The basic equations of piezoelectricity with temperature change are formulated neatly into a state equation and an output equation in terms of the generalized displacement vector and generalized stress vectors. The formalism bears a remarkable resemblance to its elastic counterpart. Various problems of piezothermoelasticity can be solved by simple extension of the corresponding solutions of anisotropic elasticity. For illustration, some fundamental problems are studied within the context and exact solutions are obtained in a systematic and self-contained manner.     

Solution of contact problems on the basis of a refined theory of plates and shells

Solutions of contact mixed boundary-value problems for a plate and for a cylindrical shell are given. These solutions are obtained with the use of equations for shells constructed by expanding solutions of elasticity theory equations with respect to the Legendre polynomials. Results of numerical simulations of the stress state in the vicinity of points with changing conditions on the frontal faces of the shell are presented. The results obtained are compared with analytical solutions of elasticity theory problems and with solutions obtained on the basis of the classical equations of the shell theory. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 169–176, September–October, 2008.     

Three-Dimensional Elasticity Solutions for Rectangular Orthotropic Plates

A state space approach for three-dimensional analysis of rectangular orthotropic elastic plates subjected to external loads on the top and bottom faces is developed. Through Hamiltonian variational formulation via Legendre’s transformation, the basic equations of elasticity are formulated into the state space framework in which the state equation exhibits Hamiltonian characteristics and the associated eigensystem possesses symplectic orthogonality. By means of separation of variables and eigenfunction expansion, three-dimensional elasticity solutions for orthotropic rectangular plates with two opposite edges simply supported and the other two arbitrary—which can be any combinations of simply-supported, clamped, and traction-free edges—are determined in a systematic way. The existing elasticity solution for the fully simply-supported plate is recovered. The through-thickness variations of the displacements and stresses are evaluated within the context.     

A state space formalism for anisotropic elasticity. Part I: Rectilinear anisotropy

A state space formalism for anisotropic elasticity including the thermal effect is developed. A salient feature of the formalism is that it bridges the compliance-based and stiffness-based formalisms in a natural way. The displacement and stress components and the thermoelastic constants of a general anisotropic elastic material appear explicitly in the formulation, yet it is simple and clear. This is achieved by using the matrix notation to express the basic equations and grouping the stress in such a way that it enables us to cast neatly the three-dimensional equations of anisotropic elasticity into a compact state equation and an output equation. The homogeneous solution to the state equation for the generalized plane problem leads naturally to the eigen relation and the sextic equation of Stroh. Extension, twisting, bending, temperature change and body forces are accounted for through the particular solution. Based on the formalism the general solution for generalized plane strain and generalized torsion of an anisotropic elastic body are determined in an elegant manner.     

On displacement methods in planar anisotropic elasticity

Two displacement formulation methods are presented for problems of planar anisotropic elasticity. The first displacement method is based on solving the two governing partial differential equations simultaneously/ This method is a recapitulation of the orignal work of Eshelby, Read and Shockley [7] on generalized plane deformations of anisotropic elastic materials in the context of planar anisotropic elasticity.The second displacement method is based on solving the two governing equations separately. This formulation introduces a displacement function, which satisfies a fourth-order partial differential equation that is identical in the form to the one given by Lekhnitskii [6] for monoclinic materials using a stress function. Moreover, this method parallels the traditional Airy stress function method and thus the Lekhnitskii method for pure plane problems. Both the new approach and the Airy stress function method start with the equilibrium equations and use the same extended version of Green's theorem (Chou and Pagano [13], p. 114; Gao [11]) to derive the expressions for stress or displacement components in terms of a potential (stress or displacement) function (see also Gao [10, 11]). It is therefore anticipated that the displacement function involved in this new method could also be evaluated from measured data, as was done by Lin and Rowlands [17] to determine the Airy stress function experimentally.The two different displacement methods lead to two general solutions for problems of planar anisotropic elasticity. Although the two solutions differ in expressions, both of the depend on the complex roots of the same characteristic equation. Furthermore, this characteristic equation is identical to that obtained by Lekhnitskii [6] using a stress formulation. It is therefore concluded that the two displacement methods and Lekhnitskii's stress method are all equivalent for problems of planar anisotropic elasticity (see Gao and Rowlands [8] for detailed discussions).     

Explicit frequency equations of free vibration of a nonlocal Timoshenko beam with surface effects

Considerations of nonlocal elasticity and surface effects in micro-and nanoscale beams are both important for the accurate prediction of natural frequency. In this study, the governing equation of a nonlocal Timoshenko beam with surface effects is established by taking into account three types of boundary conditions: hinged–hinged, clamped–clamped and clamped–hinged ends. For a hinged–hinged beam, an exact and explicit natural frequency equation is obtained. However, for clamped–clamped and clamped–hinged beams, the solutions of corresponding frequency equations must be determined numerically due to their transcendental nature. Hence, the Fredholm integral equation approach coupled with a curve fitting method is employed to derive the approximate fundamental frequency equations, which can predict the frequency values with high accuracy. In short,explicit frequency equations of the Timoshenko beam for three types of boundary conditions are proposed to exhibit directly the dependence of the natural frequency on the nonlocal elasticity, surface elasticity, residual surface stress, shear deformation and rotatory inertia, avoiding the complicated numerical computation.     

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.     

圆柱型各向异性弹性力学平面问题

本文对圆柱型各向异性弹性力学平面问题的基本方程进行了改写。在此基础上，导出了应力函数Ｇ和位移函数φ，它们满足相同的控制方程，比文〔１〕的应力函数Ｆ的控制方程要简单，便于求得特解，并有Ｆ＝ｒＧ的关系。还对若干经典问题进行了求解。     

Equilibrium of an internal transverse crack in a semiinfinite elastic body with thin coating

Static elasticity problems for a half-plane and a strip weakened by a rectilinear transverse crack are studied. In each case, the upper boundary of the body is reinforced by a flexible patch. Various versions of conditions on the lower boundary are considered in the case of the strip. The crack is maintained in the open state by distributed normal forces. The method of generalized integral transforms reduces solving the problem for the equations of equilibriumto solving a singular integral equation of the first kind with the Cauchy kernel with respect to the derivative of the crack opening function. The solutions of the integral equation are constructed by the small parameter and collocation methods for various combinations of the geometric and physical parameters of the problem, and the structure of the solutions is analyzed. The values of the stress intensity factor (SIF) near the crack vertex are obtained.     

Research on an orthogonal relationship for orthotropic elasticity

IntroductionAsymplecticsystematicmethodology[1- 3]forelasticitywasestablishedbyZhongWan_xie .Hepresentedcreativelythedualvectorsandthesymplecticorthogonalrelationshipandopenedaworkplatformparalleledtothetraditionalelasticity[4 - 9].AnewdualvectorandanewdualdifferentialmatrixLwerepresentedforasymplecticsystematicmethodologyfortwo_dimensionalelasticityandaneworthogonalrelationshipwasdiscoveredforisotropicplaneproblems[4 ]byLuoJian_hui.Theneworthogonalrelationshipisgeneralizedfororthotropicelas…     

压电压磁圆柱壳的自由振动分析

本文基于三维弹性理论,结合状态空间理论和离散奇异卷积算法分析了压电压磁圆柱壳的自由振动问题．圆柱壳的厚度方向被作为状态空间理论的传递方向,同时应用离散奇异卷积算法对面内域进行离散．因此,初始的偏微分运动方程被转化为由一阶常微分方程构成的状态方程．离散奇异卷积算法的引入使得本方法可以处理不同边界条件,从而扩展了常规状态空间方法的应用范围．本文对数值算例的计算验证了此方法的有效性和精确度．     

球面各向同性弹性力学的位移解法

本文引入三个位移函数（ｗ，Ｇ，ψ），将球面各向同性弹性力学运动方程，简化为关于ψ的二阶偏微分方程，和关于Ｗ和Ｇ的联立方程。在静力学问题中，联立方程可进一步简化，ｗ和Ｇ可用另一位移函数Ｆ表示，而Ｆ满足一个四阶偏微分方程。在球壳固有振动问题中，则简化为一个独立的二阶常微分方程，和另两个二阶的联立的常微分方程，证明了在多层球壳中，它们分别对应独立的两类振动。改进了常微分方程的解法，并计算了一个二层球壳的频率。     

Characteristic equation solution strategy for deriving fundamental analytical solutions of 3D isotropic elasticity

A simple characteristic equation solution strategy for deriving the fundamental analytical solutions of 3D isotropic elasticity is proposed. By calculating the determinant of the differential operator matrix obtained from the governing equations of 3D elasticity, the characteristic equation which the characteristic general solution vectors must satisfy is established. Then, by substitution of the characteristic general solution vectors, which satisfy various reduced characteristic equations, into various reduced adjoint matrices of the differential operator matrix, the corresponding fundamental analytical solutions for isotropic 3D elasticity, including Boussinesq-Galerkin (B-G) solutions, modified Papkovich-Neuber solutions proposed by Min-zhong WANG (P-N-W), and quasi HU Hai-chang solutions, can be obtained. Furthermore, the independence characters of various fundamental solutions in polynomial form are also discussed in detail. These works provide a basis for constructing complete and independent analytical trial functions used in numerical methods.     

On the asymptotic behavior of a phase-field model for elastic phase transitions

We study the asymptotic behavior of a one-dimensional, dynamical model of solid-solid elastic transitions in which the phase is determined by an order parameter. The system is composed of two coupled evolution equations, the mechanical equation of elasticity which is hyperbolic and a parabolic equation in the order parameter. Due to the strong coupling and the lack of smoothing in the hyperbolic equation, the asymptotic behavior of solutions is difficult to determine using standard methods of gradient-like systems. However, we show that under suitable assumptions all solutions approach the equilibrium set weakly, while the phase field stabilizes strongly.     

On solutions of the problem in stresses with the use of Maxwell stress functions

The spatial problems of elasticity are mainly solved in displacements [1, 2], i.e., the Lamé equations are taken as the initial equations. This is related to the lack of general solutions for the system of basic equations of elasticity expressed in stresses. In this connection, a new variational statement of the problem in stresses was developed in [3, 4]; this statement consists in solving six generalized equations of compatibility for six independent components of the stress tensor, while the three equilibrium equations are transferred to the set of boundary conditions. This method is more convenient for the numerical solution of problems in stresses and has been tested when solving various boundary value problems. In the present paper, analyzing the completeness of the Saint-Venant identities and using the Maxwell stress functions, we obtain a new resolving system of three differential equations of strain compatibility for the three desired stress functions φ, ξ, and ψ. This system is an alternative to the three Lamé equilibrium equations for three desired displacement components u, v, w and is simpler in structure. Moreover, both of these systems of resolving equations can be solved by the new recursive-operator method [5, 6]. In contrast to well-known methods for constructing general solutions of linear differential equations and their systems, the solutions obtained by the recursive-operator method are constructed as operator-power series acting on arbitrary analytic functions of real variables (not necessarily harmonic), and the series coefficients are determined from recursive relations (matrix in the case of systems of equations). The arbitrary functions contained in the general solution can be determined directly either from the boundary conditions (the obtained system of inhomogeneous equations with a right-hand side can also be solved by the recursive-operator method [6]) or by choosing them from various classes of analytic functions (elementary, special); a complete set of particular solutions can be obtained in the same function classes, and the coefficients of linear combinations of particular solutions can be determined by the Trefftz method, the least-squares method, and the collocation method.     

The solution of weak formulation for axisymmetric problem of orthotropic cantilever cylindrical shell

Weak formulation of equilibrium equations including boundary conditions of laminated cylindrical shell is presented, and thermal stresses mixed state equation for axisymmetric problem of closed cantilever cylindrical shell is established. A unified approach and weak solutions are obtained for closed laminated cantilever cylindrical shell of arbitrary thickness under thermal and mechanical loadings. The equation and boundary conditions proposed in this paper are weakened, the method of this paper would be easy to popularize in dynamics analysis of elasticity. Foundation item: the National Natural Science Foundation of China (19392305)     

A penny shaped crack in a transversely isotropic material under non-axisymmetric impact loads

A solution is given for problems involving non-axisymmetric dynamic impact loading of a penny shaped crack in a transversely isotropic medium. Laplace and Hankel transforms are used to reduce the equations of elasticity to integral equations, and solutions are obtained for the three modes of fracture. The stress intensity factors are determined for a penny shaped crack loaded by concentrated normal impact forces and concentrated radial shear impact forces. The integral equations are solved by numerical methods, and the results are plotted showing how the dynamic stress intensity factors are influenced by the asymmetric loading.     

具有固支边的强厚度叠层板的精确解

抛弃任何有关位移或应力模式的人为假设,在文献[1]、[2]的基础上,引入δ-函数,对具有固支边的强厚度叠层矩形板在任意荷载作用下建立其状态方程。给出静力、动力和稳定问题的精确解,其解满足弹性力学所有方程,并计及了所有弹性常数。     

Lattice with long-range interaction of power-law type for fractional non-local elasticity

Lattice models with long-range interactions of power-law type are suggested as a new type of microscopic model for fractional non-local elasticity. Using the transform operation, we map the lattice equations into continuum equation with Riesz derivatives of non-integer orders. The continuum equations that are obtained from the lattice model describe fractional generalization of non-local elasticity models. Particular solutions and correspondent asymptotic of the fractional differential equations for displacement fields are suggested for the static case.