功能梯度矩形板问题三维理论的半解析解

针对功能梯度材料矩形板问题,基于三维弹性理论,将位移和应力分量作为基本变量,通过双三角级数将其控制微分方程转化为常微分方程组的边值问题。采用插值矩阵法直接对常微分方程组边值问题进行求解,得到了功能梯度材料矩形板三维位移、应力场的半解析解。通过算例给出了材料参数按指数形式和幂函数形式变化情况下的功能梯度板的弯曲问题。对比有限元法和状态空间法,结果表明：本文提出的状态空间与插值矩阵法相结合的半解析法能有效地分析材料参数按任意形式连续变化的功能梯度矩形板问题,且具有良好的精度,精度可达10-4量级,能够满足工程需要;与其他方法相比,本文方法具有实施便捷、计算量小等优点,根据其力学场分析结果可设计出满足各种不同需求的功能梯度材料。     

圆板在物体撞击下的非线性动力响应

本文在Von Kármán大位移的意义上,利用虚位移原理伽辽金方法建立了圆板在物体撞击下的非线性动力响应的控制微分方程,在研究响应问题时,考虑了冲击载荷与圆板位移响应之间的耦合影响,文中使用时间增量法和奇异摄动理论求解问题的控制方程,获得了固支圆板非线性动力响应的近似解,并且求解了具体算例,绘出了圆板位移、应力响应曲线以及冲击力随时间的变化曲线。     

圆板在物体撞击下的非线性动力响应

本文在Von Kármán大位移的意义上,利用虚位移原理伽辽金方法建立了圆板在物体撞击下的非线性动力响应的控制微分方程,在研究响应问题时,考虑了冲击载荷与圆板位移响应之间的耦合影响,文中使用时间增量法和奇异摄动理论求解问题的控制方程,获得了固支圆板非线性动力响应的近似解,并且求解了具体算例,绘出了圆板位移、应力响应曲线以及冲击力随时间的变化曲线。     

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

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

各向异性复合材料尖劈和接头的奇性应力指数研究

提出了一个新的、基于位移的、求解三维尖劈端部奇性应力指数问题的非协调元特征分析法。该方法假定尖劈端部邻域内的位移场没有采用奇异变换技术，导出虚功方程的出发点不同于过去原有求解裂纹尖端近似场的有限元特征分析法，在有限元离散时采用的单元形式为非协调元。文中运用该方法给出了若干求解各向异性复合材料尖劈／接头端部奇性应力指数的算例。所有的计算结果表明，本文方法能够求解复杂尖劈／接头的全部奇性应力指数，使用的单元少而且精度高。     

板问题混合变量等参Hamiltonian元的半解析解

本文给出板问题混合变量方程的一种半离散半解析方法－混合变量等参Ｈａｍｉｌｔｏｎｉａｎ元。该方法沿板厚方向未作任何有关位移和应力的人为假设，而是采用控制论中方法给出真解，所以可以求解任意厚度板问题。     

无限粘弹性平面中椭圆孔口变边界过程的复变函数解答

工程中存在一类几何边界随时间变化的变边界结构,例如土木工程中处于施工阶段的结构。本文以粘弹性岩体中隧道开挖为背景,尝试用变边界问题对应关系和平面弹性复变方法求取无限平面中椭圆孔口自相似变边界情况下的解析解答。首先建立了复变函数法求解变边界粘弹性问题的基本步骤和公式。然后通过建立逆映射函数将已知?平面复位势转至z平面,从而解耦参与拉普拉斯变换的时间与孔口映射函数所带来的时间,从而导出了粘弹性类材料的应力与位移的统一表达。作为一个例子,本文选择Boltzmann粘弹性模型,代入模型参数后得到积分形式的位移、应力解析解,通过与数值解的比较验证了该解答的可靠性,并通过一个算例分析了变边界过程对位移、应力的影响。分析结果显示,采用不同变边界过程的位移、应力变化形态和数值均有差别。本文解答可用于进行地下椭圆孔型隧道在开挖过程中的力学分析,为实际工程提供初步设计的手段。此外,本文给出的方法可用于推导任意形状孔型变边界问题的解答。     

横观各向同性体三维裂纹的瞬态扩展

讨论一对集中力作用下横观各向同性体三维裂纹的瞬态扩展问题,其解答构成三维裂纹瞬态扩展问题的基本解。求解方法是基于积分变换技术,将混合边值问题化为Wiener-Hopf型积分方程,求得了裂纹所在平面应力和位移的封闭形式解。进一步利用Abel定理和Cagniard-de Hoop方法,求得了动态应力强度因子的精确解。最后通过数值结果揭示了横观各向同性材料三维扩展裂纹尖端场的动态特性。     

平板裂纹前缘三维效应区的弹性力学分析

本文利用Kane和Mindlin关于弹性乎板面内问题位移的基本假设及Fourier变换求解了无限大板的I型裂纹问题,得到了裂纹尖端应力位移场的渐近形式、应力强度因子随板厚的变化规律以及板厚对裂纹前缘三维效应区的影响.研究表明,对线性硬化材料,在塑性切线胰量不太小的情况下,线弹性分析的结果可近似适用于弹塑性材料.     

厚圆板轴对称自由振动的三维理论解

本文用三维弹性动力学理论研究厚圆板在柱面约束ｕ＝τｒｚ＝０下的轴对称自由振动问题，给出位移和应力振型的解析表达式，求得固有频率的精确值。研究结果表明：在板的厚度方向存在着不同半波数的反对称形式和对称形式的位移振型，就是说介质之间在厚度方向具有不同的运动形式，这是各种厚板理论都无法反映的。但是对于圆板的轴对称基本振动来说，Ｈｅｎｃｋｙ、Ｒｅｉｓｓｎｅｒ和Ｍｉｎｄｌｉｎ理论求得的固有频率与本文的弹性动     

Nonlinear three-dimension analysis for axially symmetrical circular plates and multilayered plates

Analytic nonlinear three-dimension solutions are presented for axially symmetrical homogeneous isotropic circular plates and multilayered plates with rigidly clamped boundary conditions and under transverse load.The geometric nonlinearily from a moderately large deflection is considered.A developmental perturbation method is used to solve the complicated nonlinear three-dimension differential equations of equilibrium.The basic idea of this perturbation method is using the two-dimension solutions as a basic form of the corresponding three-dimension solutions,and then processing the perturbation procedure to obtain the three-dimension perturbation solutions.The nonlinear three-dimension results in analytic expressions and in numerical forms for ordinary plates and multilayered plates are presented.All of the plate stresses are shown in figures.The results show that this perturbation method used to analyse nonlinear three-dimension problems of plates is effective.     

NONLINEAR THREE-DIMENSIONAL ANALYSIS OF COMPOSITELAMINATED PLATES

ＮＯＮＬＩＮＥＡＲＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＡＮＡＬＹＳＩＳＯＦＣＯＭＰＯＳＩＴＥＬＡＭＩＮＡＴＥＤＰＬＡＴＥＳ￥(江晓禹，张相周，陈百屏)ＪｉａｎｇＸｉａｏｙｕ；（ＳｏｕｔｈｗｅｓｔｅｒｎＪｉａｏｔｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｃｈｅｎｇｄｕ６...     

Static equilibrium of a thermoelastic half-plane: Green’s functions and solutions in integrals

This article presents in a closed form new influence functions of a unit point heat source on the displacements for three boundary value problems of thermoelasticity for a half-plane. We also obtain the corresponding new integral formulas of Green’s and Poisson’s types that directly determine the thermoelastic displacements and stresses in the form of integrals of the products of specified internal heat sources or prescribed boundary temperature and constructed already thermoelastic influence functions (kernels). All these results are presented in terms of elementary functions in the form of three theorems. Based on these theorems and on derived early by author the general Green-type integral formula, we obtain in elementary functions new solutions to two particular boundary value problems of thermoelasticity for half-plane. The graphical presentation of the temperature and thermal stresses of one concrete boundary value problems of thermoelasticity for half-plane also is included. The proposed method of constructing thermoelastic Green’s functions and integral formulas is applicable not only for a half-plane, but also for many other two- and three-dimensional canonical domains of different orthogonal coordinate systems.     

The periodic cracks of an infinite anisotropic media for plane skew-symmetric loadings

This paper attempts to solve the periodic crack problems of infinitive anisotropic media for plane skew-symmetric loadings by means of the method of complex function. The problems are now reduced to the determination of two complex functions that must satisfy certain boundary conditions. In this paper, the stresses, the displacements and the boundary conditions are assumed to be periodic, and further, the stresses are assumed to be bounded at infinity. The solutions are expressed in closed forms.     

加权残数配点法解正交各向异性板的积分方程

本文推导了一般各向异性板弯曲的积分方程,运用加权残数配点法求解了正交各向异性板弯曲的积发方程,本文将部分配点取在边界上,另一部分配点取在域外,只用关于找度的基本积分方程,而不用关于转角的补充积分方程,简化了方程求解和计算程序,由于正交各向异性板没有争析形式的、实用的基本解,本文提出了两种新的近似基本解;加权双三角级数;广义各向同性板解析形式的基本解和加权双三角级数的叠加,算例表明,本文提出的解法和近似基本解适用于各类边界条件的正交各向异性板,具有简单、可靠、精度高等优点。     

弹性薄板分析的条形传递函数方法

提出一种用于矩形弹性薄板变形分析的条形传递函数方法．一个矩形区域首先沿某一个方向被剖分成若干个条形子域,分割这些子域的直线称为结线,在结线上定义位移函数,它是结线坐标的一维函数,结线的两个端点称为结点．为适应复杂边界条件,在边界结线上定义若干结点,该结线的位移函数用结点位移参数插值表示．每个条形子域的变形用结线位移函数和适当的插值函数（形函数）表示．结线位移函数和结点位移参数满足的平衡微分方程及代数方程由变分原理给出     

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.     

Interactions of Penny-shaped Cracks in Three-dimensional Solids

The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is developed in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.The English text was polished by Keren Wang     

Thermoelastostatic equilibrium of a spatial quadrant: Green’s function and solution in integrals

In this study, a new Green??s function and a new Green-type integral formula for a 3D boundary value problem (BVP) in thermoelastostatics for a quarter-space are derived in closed form. On the boundary half-planes, twice mixed homogeneous mechanical boundary conditions are given. One boundary half-plane is free of loadings and the normal displacements and the tangential stresses are zero on the other one. The thermoelastic displacements are subjected by a heat source applied in the inner points of the quarter-space and by mixed non-homogeneous boundary heat conditions. On one of the boundary half-plane, the temperature is prescribed and the heat flux is given on the other one. When the thermoelastic Green??s function is derived, the thermoelastic displacements are generated by an inner unit point heat source, described by ??-Dirac??s function. All results are obtained in elementary functions that are formulated in a special theorem. As a particular case, when one of the boundary half-plane of the quarter-space is placed at infinity, we obtain the respective results for half-space. Exact solutions in elementary functions for two particular BVPs for a thermoelastic quarter-space and their graphical presentations are included. They demonstrate how to apply the obtained Green-type integral formula as well as the derived influence functions of an inner unit point body force on volume dilatation to solve particular BVPs of thermoelasticity. In addition, advantages of the obtained results and possibilities of the proposed method to derive new Green??s functions and new Green-type integral formulae not for quarter-space only, but also for any canonical Cartesian domain are also discussed.     

Application of the perturbation method in mechanics of deformable solids

Although the solutions of the classical problems of continuum mechanics have been studied sufficiently well, the smallest deviations, for example, of the body boundary or of the material characteristics from the traditional values prevent one from obtaining exact solutions of these problems. In this case, one has to use approximate methods, the most common of which is the perturbation method. The problems studied in [1–6] belong to classical problems in which the perturbation method is used to study the behavior of deformable bodies. A wide survey of studies analyzing the perturbations of the body boundary shape caused by variations in its stress-strain state is given in [5, 6]. In numerous studies, it was noted that the problem on the convergence of approximate solutions and hence the studies of the continuous dependence of the solution of the original problem on the characteristics of perturbations (“imperfections”) play an important role. In the present paper, we analyze the forms of mathematical models of deformable bodies by studying whether the solution of the original problem continuously depends on the characteristics of the perturbed shape of the body boundary on which the boundary conditions are posed in terms of stresses and on the characteristics of the material properties. We use the results of this analysis to conclude that, when using the perturbation method, one should state the boundary conditions in terms of stresses on the boundary of the real body in stressed state.