Solvability on boundary-value problems of elasticity of three-dimensional quasicrystals

Weak solution (or generalized solution) for the boundary-value problems of partial differential equations of elasticity of 3D (three-dimensional) quasicrystals is given, in which the matrix expression is used. In terms of Korn inequality and theory of function space,we prove the uniqueness of the weak solution.This gives an extension of existence theorem of solution for classical elasticity to that of quasicrystals,and develops the weak solution theory of elasticity of 2D quasicrystals given by the second author of the paper and his students.     

准晶数学弹性力学和缺陷力学

对准晶数学弹性理论的基本概念和基本框架作了介绍,在此基础上分别针对目前已经发现的几类一维准晶、二维准晶和三维准晶讨论了其数学弹性的理论体系．为了求解准晶弹性的边值问题或初值一边值问题,还必须发展相应的方法论．物理工作者在研究准晶位错弹性问题中发展了Ｇｒｅｅｎ函数方法．针对一维与二维准晶弹性中几类问题提出了分解与叠加程序,这一程序的使用,使极其复杂的准晶弹性问题得到简化,进而引进位移函数或应力函数,把数目。庞大的准晶弹性基本方程化成一个或少数几个高阶偏微分方程,进一步使求解步骤大为简化．对三维立方准晶弹性也采用了类似步骤使求解过程大为简化．在以上化简的基础上,发展了准晶弹性的边值问题或初值一边值问题的复交函数方法和 Ｆｏｕｒｉｅｒ分析方法,求得了一系列准晶位错问题和裂纹问题的分析解（古典解）．在研究准晶弹性的边值问题古典解的同时,也讨论了同这些边值问题相对应的变分问题和广义解（弱解）以及这种弱解的数值方法──有限元法．在物理学家工作基础上开展的这些工作可以看作对经典数学弹性理论和方法、经典Ｖｏｌｔｅｒｒａ位错理论、普通结构材料断裂力学和经典有限元的某些发展．此外,还把一维六方准晶弹性动力学的结果与统计物理的某些     

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

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

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).     

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     

十二次对称二维准晶中的无摩擦接触问题

利用积分变换的方法讨论了在一个刚性压头作用下十二次对称二维准晶的无摩擦接触问题. 通过引入位移势函数,将数量巨大而复杂的偏微分方程转化为两个独立的双调和方程,应用Fourier分析与对偶积分方程理论解决了十二次对称二维准晶材料的无摩擦接触问题,得到了相应的接触应力解析表达式,结果表明：如果接触位移是一常数,则接触应力在接触区域边缘具有-1/2阶奇异性;反之,如果接触应力在接触区域边缘具有-1/2阶的奇异性,则接触位移一定为一常数,这为准晶材料的接触变形提供了重要的力学参数.     

A Theory of General Solutions of Plane Problems in Two-Dimensional Octagonal Quasicrystals

A theory of general solutions of plane problems is developed for the coupled equations in plane elasticity of two-dimensional octagonal quasicrystals. In virtue of the operator method, the general solutions of the antiplane and inplane problems are given constructively with two displacement functions. The introduced displacement functions have to satisfy higher order partial differential equations, and therefore it is difficult to obtain rigorous analytic solutions directly and is not applicable in most cases. In this case, a decomposition and superposition procedure is employed to replace the higher order displacement functions with some lower order displacement functions, and accordingly the general solutions are further simplified in terms of these functions. In consideration of different cases of characteristic roots, the general solution of the antiplane problem involves two cases, and the general solution of the inplane problem takes three cases, but all are in simple forms that are convenient to be applied. Furthermore, it is noted that the general solutions obtained here are complete in x ₃-convex domains.      

MOVING SCREW DISLOCATION IN CUBIC QUASICRYSTAL

The elasticity theory of the dislocation of cubic quasicrystals is developed. The governing equations of anti-plane elasticity dynamics problem of the quasicrystals were reduced to a solution of wave equations by introducing displacement functions, and the analytical expressions of displacements, stresses and energies induced by a moving screw dislocation in the cubic quasicrystalline and the velocity limit of the dislocation were obtained. These provide important information for studying the plastic deformation of the new solid material.     

Analysis of plate in second strain gradient elasticity

The bending analysis of a thin rectangular plate is carried out in the framework of the second gradient elasticity. In contrast to the classical plate theory, the gradient elasticity can capture the size effects by introducing internal length. In second gradient elasticity model, two internal lengths are present, and the potential energy function is assumed to be quadratic function in terms of strain, first- and second-order gradient strain. Second gradient theory captures the size effects of a structure with high strain gradients more effectively rather than first strain gradient elasticity. Adopting the Kirchhoff’s theory of plate, the plane stress dimension reduction is applied to the stress field, and the governing equation and possible boundary conditions are derived in a variational approach. The governing partial differential equation can be simplified to the first gradient or classical elasticity by setting first or both internal lengths equal to zero, respectively. The clamped and simply supported boundary conditions are derived from the variational equations. As an example, static, stability and free vibration analyses of a simply supported rectangular plate are presented analytically.     

General plane elasticity solution for non-homogeneous materials

In this study, we consider the general plane elasticity problem for non-homogeneous materials subjected to mechanical and thermal loads. The general equations are derived in x–y coordinate system. It is shown that, if the non-homogeneity varies exponentially, then the governing equations are reduced to partial differential equations with constant coefficients. To illustrate the validity and usefulness of the formulation, problems involving both mechanical and thermal loads with the non-homogeneity varying in both directions are considered. As examples, infinite non-homogeneous beams subjected to axial load and bending moment are considered and the solutions obtained in closed form. Then an example illustrating the effect of temperature change is presented. Again a closed form solution is obtained for the stress distribution. Also for limiting values of the non-homogeneity parameter, extensive asymptotic analyses are presented. In addition two problems where the Young’s Modulus or the coefficient of thermal expansion coefficient varies in both directions are presented. From the examples solved, it is found that non-homogeneity affects the stress distribution significantly, and leads to counterintuitive results with serious implication to the study of fracture problems in functionally graded materials.     

非局部平面应变和平面应力问题界定及其精确性讨论

在非局部弹性理论框架下对平面应变和平面应力状态重新界定.首先,分别在其相应简化假设下推导控制方程,并与经典局部情况进行比较.然后,引入变形协调条件对两类非局部平面问题的精确性进行讨论.其中,对于非局部平面应力状态,通过应变协调方程的Fourier变换形式来进行研究,使问题得以简化.通过以上分析,最终得到一些有价值的结论.     

On the fundamental equations of piezoelasticity of quasicrystal media

The three-dimensional fundamental equations of elasticity of quasicrystals with extension to quasi-static electric effect are expresses in both differential and variational invariant forms for a regular region of quasicrystal material. The principle of conservation of energy is stated for the regular region and the constitutive relations are obtained for the piezoelasticity of material. A theorem is proved for the uniqueness in solutions of the fundamental equations by means of the energy argument. The sufficient boundary and initial conditions are enumerated for the uniqueness. Hamilton’s principle is stated for the regular region and a three-field variational principle is obtained under some constraint conditions. The constraint conditions, which are generally undesirable in computation, are removed by applying an involutory transformation. Then, a unified variational principle is obtained for the regular region, with one or more fixed internal surface of discontinuity. The variational principle operating on all the field variables generates all the fundamental equations of piezoelasticity of quasicrystals under the symmetry conditions of the phonon stress tensor and the initial conditions. The resulting equations, which are expressible in any system of coordinates and may be used through simultaneous approximation upon all the field variables in a direct method of solutions, pave the way to the study of important dislocation, fracture and interface problems of both elasticity and piezoelasticity of quasicrystal materials.     

Nonlinear stability and dynamics of composite skew plates under nonuniform loadings using differential quadrature method

The buckling, postbuckling and postbuckled vibration behaviour of composite skew plates subjected to nonuniform inplane loadings are presented here. The skew plate is modelled using first order shear deformation theory accounting for von-Kármán geometric nonlinearity and initial geometric imperfections. The different types of nonuniform loads that have been considered in this study are concentrated load, partial load and parabolic load. The explicit analytical expressions for prebuckling stress distributions within composite skew plate subjected to three different types of nonuniform inplane loadings are developed by solving plane elasticity problem using Airy's stress function approach. It is observed that the inplane normal stress distributions within the skew plate due to above nonuniform loadings do not become uniform even at mid-section. The generalized differential quadrature (GDQ) method is used to solve the nonlinear governing partial differential equations. It is observed that the postbuckled load carrying capacity of skew plate under concentrated loading is the lowest compared to other nonuniform and uniform loadings.     

Plane Deformations in Incompressible Nonlinear Elasticity

The substantially general class of plane deformation fields, whose only restriction requires that the angular deformation not vary radially, is considered in the context of isotropic incompressible nonlinear elasticity. Analysis to determine the types of deformations possible, that is, solutions of the governing systems of nonlinear partial differential equations and constraint of incompressibility, is developed in general. The Mooney-Rivlin material model is then considered as an example and all possible solutions to the equations of equilibrium are determined. One of these is interpreted in the context of nonradially symmetric cavitation, i.e., deformation of an intact cylinder to one with a double-cylindrical cavity. Results for general incompressible hyperelastic materials are then discussed. The novel approach taken here requires the derivation and use of a material formulation of the governing equations; the traditional approach employing a spatial formulation in which the governing equations hold on an unknown region of space is not conducive to the study of deformation fields containing more than one independent variable. The derivation of the cylindrical polar coordinate form of the equilibrium equations for the nominal stress tensor (material formulation) for a general hyperelastic solid and a fully arbitrary cylindrical deformation field is also given. This revised version was published online in August 2006 with corrections to the Cover Date.     

Bending analyses of 1D orthorhombic quasicrystal plates

The meshless Petrov–Galerkin method (MLPG) is applied to plate bending analysis in 1D orthorhombic quasicrystals (QCs) under static and transient dynamic loads. The Bak and elasto-hydrodynamic models are applied for phason governing equation in the elastodynamic case. The phason displacement for the orthorhombic QC in the first-order shear deformation plate theory depends only on the in-plane coordinates on the mean plate surface. Nodal points are randomly distributed over the mean surface of the considered plate. Each node is the center of a circle surrounding this node. The coupled governing partial differential equations are satisfied in a weak-form on small fictitious subdomains. The spatial variations of the phonon and phason displacements are approximated by the moving least-squares (MLS) scheme. After performing the spatial MLS approximation, a system of ordinary differential equations (ODEs) for nodal unknowns is obtained. The system of the ODEs of the second order is solved by the Houbolt finite-difference scheme. Our numerical examples demonstrate clearly the effect of the coupling parameter on both static and dynamic phonon/phason deflections.     

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

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

弹性基础上对称叠层厚板的基本解

利用平面波分解法将弹性地基础上计入剪切变形的对称叠层厚板的偏微分方程组转化常微分方程组，然后利用Ｈ¨ｏｒｍａｎｄｅｒ算子法将常微分方程组转化为一个六阶的常微分方程，以定积分的形式提出了弹性基础上对称叠层厚板的基本解．     

用路径守恒积分计算平面准晶裂纹扩展的能量释放率

将准晶（包括点群５,５ｍ,８ｍｍ,１０,１０ｍｍ,１２ｍｍ）弹性边值问题化为广义能量泛函的变分问题,建立了求解准晶弹性变形的有限元法,并提出了含裂纹准晶的路径守恒Ｅ－积分,指出Ｅ－积分在数值上等于准晶裂纹扩展的能量释放率,作为实例,计算了平面五次、八次、十二次对称准晶裂纹扩展的能量释放率,数值计算表明准晶中路径积分具有较好的守恒性,同时数值结果表明准晶中相位子场的出现使准晶裂纹扩展的能量释放率增大     

Superposition of Generalized Plane Strain on Anti-Plane Shear Deformations in Isotropic Incompressible Hyperelastic Materials

The purpose of this research is to investigate the basic issues that arise when generalized plane strain deformations are superimposed on anti-plane shear deformations in isotropic incompressible hyperelastic materials. Attention is confined to a subclass of such materials for which the strain-energy density depends only on the first invariant of the strain tensor. The governing equations of equilibrium are a coupled system of three nonlinear partial differential equations for three displacement fields. It is shown that, for general plane domains, this system decouples the plane and anti-plane displacements only for the case of a neo-Hookean material. Even in this case, the stress field involves coupling of both deformations. For generalized neo-Hookean materials, universal relations may be used in some situations to uncouple the governing equations. It is shown that some of the results are also valid for inhomogeneous materials and for elastodynamics.     

无网格方法在平面粘弹性力学问题中的应用

本文综合应用无网格方法(EFGM)、线性粘弹性与弹性力学之间的对应原理,Laplace变换和逆变换等方法求解了拟静态平面弹性和粘弹性力学问题。首先,利用Laplace变换和逆变换推导了平面问题的粘弹性本构关系,建立了拟静态粘弹性平面问题的边值问题;其次,利用粘弹性与弹性力学之间的对应原理得到了Laplace变换域中平面问题的基本方程,在Laplace变换域中建立了相应的泛函,并得到了用无网格方法离散的控制方程;同时,求解了几个拟静态弹性和粘弹性平面问题,给出了它们的表达式和数值结果;最后,采用Laplace逆变换和数值逆变换,得到了粘弹性力学平面问题在物理空间中的解,并比较了由解析解和无网格数值方法所得到的数值结果,可以看到它们是非常吻合的。说明本文方法的正确性和有效性。