板弯曲求解新体系及其应用1)

建立平面弹性与板弯曲的相似性理论,给出了板弯曲经典理论的另一套基本方程与求解方法,然后进入哈密顿体系用直接法研究板弯曲问题.新方法论应用分离变量、本征函数展开方法给出了条形板问题的分析解,突破了传统半逆解法的限制.结果表明新方法论有广阔的应用前景.     

功能梯度矩形板的三维弹性分析

将功能梯度三维矩形板的位移变量按双三角级数展开，以弹性力学的平衡方程为基础．导出位移形式的平衡方程。引入状态空间方法，以三个位移分量及位移分量的一阶导数为状态变量，建立状态方程。考虑四边简支的边界条件，由状态方程得到了功能梯度三维矩形板的静力弯曲问题和自由振动问题的精确解。由给出的均匀矩形板自由振动问题的计算结果表明．与已有的理论解以及有限元方法的计算结果相吻合。假设功能梯度三维矩形板的材料常数沿板的厚度方向按照指数函数的规律变化．进一步给出了功能梯度三维矩形板的自由振动问题和静力弯曲问题的算例分析，并讨论了材料性质的梯度变化对板的动力响应和静力响应的影响。     

功能梯度与均匀圆板弯曲解的线性转换关系

基于一阶剪切理论, 研究了功能梯度材料圆板与均匀圆板轴对称弯曲解之间 的线性转换关系. 通过理论分析和比较 功能梯度材料圆板和均匀圆板在一阶剪切理论下的位 移形式的轴对称弯曲控制方程, 发现了功能梯度材料圆板的转角与均匀圆板的转角之间 的相似转换关系. 在假设材料性质沿板厚连续变化的情况下, 给出了相似转换系数的解析表 达式. 在此基础上, 进一步导出了一阶剪切理论下功能梯度圆板的挠度与经典理论 下, 均匀圆板的挠度之间的线性关系. 从而, 可将功能梯度材料圆板在一阶剪切理 论下的弯曲问题求解, 转化为相应均匀薄圆板在经典理论下的弯曲问题求解, 以及 转换系数的计算问题. 这一方法为功能梯度非均匀中厚度圆板的求解提供了简捷途 径, 而且更便于工程应用. 作为例子, 采用上述方法分别求得了周边简支和夹紧条 件下, 梯度圆板在均布横向载荷作用下的弯曲解析解, 该解答与Reddy得到的结果 完全吻合.     

含椭圆孔有限大薄板弯曲应力分析

本文利用各向异性体弹性平面理论中的复势方法，以Faber级数为工具，对含椭圆孔有限大各向异性板弯曲问题进行应力分析，得出有限大含椭圆孔各向异性板弯曲的级数解形式，分析了有限大含椭圆孔板在受到弯曲载荷时孔边的应力分布，并讨论了各种参数对应力分布的影响，给出了有益的结论     

宽板塑性弯曲成形过程中的板厚变化规律

根据宽板弯曲过程中的变形与应力分布特征，提出了一种计算弯曲过程中板料厚度随弯曲程度变化的新的近似解答方法．该方法基于塑性增量理论，应用塑性成形过程的体积不变假设和弯曲过程的平面假设．作为算例，得到了理想刚塑性材料的板料厚度、变薄系数以及应变中性层内移系数随板料弯曲内表面曲率半径变化的规律，并与实验数据进行了比较，两吻合良好．     

应用混合变量余能原理求解不同载荷作用下四边简支矩形板的弯曲

由功的互等定理导出了弯曲矩形板混合变量的总余能,应用混合变量余能原理求解四边简支矩形板弯曲问题,给出了四边简支矩形板的混合变量余能表达式,求解了不同载荷作用下四边简支矩形板的弯曲问题,并将计算的不同载荷作用下的挠度和弯矩值与ANSYS有限元的计算值比较,给出了一种求解四边简支矩形板的新方法。     

基于平面弹性—板弯曲模拟关系的薄板有限元列式——从膜单元到薄板单元

本文全面讨论了基于平面弹性－－板弯曲模拟关系的薄板有限单元的理论和方法，由于直接对弯矩函数进行插值，c1连续性的要求得以自然避免，薄板单元可以直接在c0连续的层面上加以构造，无需借用Reissner-Mindlin的中厚板理论，由之引发的闭锁问题也得以避免，本文系统地阐明了平面弹性膜单元与薄板弯曲单元的对应关系，及由平面弹性膜单元的向薄板弯曲单元转换的一整套方法。为薄板单元的构造提供了一条新的有余的途径，文中给出了对应于平面弹性膜单元CST，LST，Q4，Q8的薄板单元，我们称之为MPS板单元，MPS板元以挠度和转角为自由度，便于实际应用，和其它板单元相比具有非常高的精度。     

求解弯曲矩形板受迫振动的混合变量法

应用混合变量最小作用量原理,求解了两邻边固定、另两邻边自由弯曲矩形板在均匀谐载作用下的受迫振动问题,得到其受迫振动的稳态解．所获得的计算结果,跟现有文献进行比较,证明其给出的方法是正确的,因此不仅具有重要的理论价值,而且可以为工程实际直接采用．     

固定式厚管板的弯曲问题

本文应用E.Reissner的厚板理论对高压固定式热交换器管板的弯曲问题进行了研究,给出了厚管板在固定和简支两种边界条件下的应力和挠度的解析公式。由实例说明,理论值与试验值相当符合,本文公式可供厚管板设计时参考。     

热/机械载荷下功能梯度材料矩形厚板的弯曲行为

采用Reddy高阶剪切板理论，考虑材料物性参数随坐标和温度变化的特性，研究在均匀变化的温度场内功能梯度材料矩形板在面内与横向载荷共同作用下的横向弯曲问题，基于一维DQ法和Galerkin技术，给出了一对边固支，另对边任意约束时板弯曲问题的半解析解，以Si3N4/SUS304板为例考察了材料组份，温度场，面内载荷及边界约束条件等对功能梯度材料板弯曲行为的影响。     

New solution system for circular sector plate bending and its application

Instead of the biharmonic type equation, a set of new governing equations and solving method for circular sector plate bending is presented based on the analogy between plate bending and plane elasticity problems. So the Hamiltonian system can also be applied to plate bending problems by introducing bending moment functions. The new method presents the analytical solutions for the circular sector plate. The results show that the new method is effective. Project supported by National Natural Science Foundation (No. 19732020) and the Doctoral Research Foundation of China.     

板弯曲与平面弹性问题的多类变量变分原理

进一步完善板弯曲与平面弹性问题的多类变量变分原理，给出了相关边界积分项的具体表达式．多类交量变分原理涵盖了平衡、应力函数、应力、位移一应变、协调和物性共五大类基本方程和所有边界条件，是一个具有更加广泛意义的变分原理．     

薄板理论的正交关系及其变分原理

利用平面弹性与板弯曲的相似性理论，将弹性力学新正交关系中构造对偶向量的思路推广到 各向同性薄板弹性弯曲问题，由混合变量求解法直接得到对偶微分方程并推导了对应的变分 原理. 所导出的对偶微分矩阵具有主对角子矩阵为零矩阵的特点. 发现了两个独立的、对称 的正交关系，利用薄板弹性弯曲理论的积分形式证明了这种正交关系的成立. 在恰当选择对 偶向量后，弹性力学的新正交关系可以推广到各向同性薄板弹性弯曲理论.     

On new symplectic elasticity approach for exact bending solutions of rectangular thin plates with two opposite sides simply supported

This paper presents a bridging research between a modeling methodology in quantum mechanics/relativity and elasticity. Using the symplectic method commonly applied in quantum mechanics and relativity, a new symplectic elasticity approach is developed for deriving exact analytical solutions to some basic problems in solid mechanics and elasticity which have long been bottlenecks in the history of elasticity. In specific, it is applied to bending of rectangular thin plates where exact solutions are hitherto unavailable. It employs the Hamiltonian principle with Legendre’s transformation. Analytical bending solutions could be obtained by eigenvalue analysis and expansion of eigenfunctions. Here, bending analysis requires the solving of an eigenvalue equation unlike in classical mechanics where eigenvalue analysis is only required in vibration and buckling problems. Furthermore, unlike the semi-inverse approaches in classical plate analysis employed by Timoshenko and others such as Navier’s solution, Levy’s solution, Rayleigh–Ritz method, etc. where a trial deflection function is pre-determined, this new symplectic plate analysis is completely rational without any guess functions and yet it renders exact solutions beyond the scope of applicability of the semi-inverse approaches. In short, the symplectic plate analysis developed in this paper presents a breakthrough in analytical mechanics in which an area previously unaccountable by Timoshenko’s plate theory and the likes has been trespassed. Here, examples for plates with selected boundary conditions are solved and the exact solutions discussed. Comparison with the classical solutions shows excellent agreement. As the derivation of this new approach is fundamental, further research can be conducted not only on other types of boundary conditions, but also for thick plates as well as vibration, buckling, wave propagation, etc.     

构造极坐标中Airy应力函数的观察法

通过引入Airy应力函数，平面问题可以归结为在给定的边界条件下求解一个双调和方程．因此对双调和函数性质的研究将有利于平面问题的求解．首先给出一个有关双调和函数的引理，并分别从复变和微分两种角度提供该引理的证明．借助这个引理，提出了一种构造极坐标中Airy应力函数的观察法．最后，举例说明了该观察法在几个经典平面问题中的应用．这些例子说明，利用本的观察法可以将某些平面问题应力函数构造的过程简单化。     

加强板的弯矩函数列式

本文首先谇薄板弯曲问题矩函数的物理意义，据此，将弯矩函数列式推广到具有加强条的薄板弯曲问题，给出了与平面弹性问题完全对应的余能原理。     

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.     

The use of complex valued functions for the solution of three-dimensional elasticity problems

For the treatment of plane elasticity problems the use of complex functions has turned out to be an elegant and effective method. The complex formulation of stresses and displacements resulted from the introduction of a real stress function which has to satisfy the 2-dimensional biharmonic equation. It can be expressed therefore with the aid of complex functions. In this paper the fundamental idea of characterizing the elasticity problem in the case of zero body forces by a biharmonic stress function represented by complex valued functions is extended to 3-dimensional problems. The complex formulas are derived in such a way that the Muskhelishvili formulation for plane strain is included as a special case. As in the plane case, arbitrary complex valued functions can be used to ensure the satisfaction of the governing equations. Within the solution of an analytical example some advantages of the presented method are illustrated.     

Elasticity solutions for plane anisotropic functionally graded beams

This paper considers the plane stress problem of generally anisotropic beams with elastic compliance parameters being arbitrary functions of the thickness coordinate. Firstly, the partial differential equation, which is satisfied by the Airy stress function for the plane problem of anisotropic functionally graded materials and involves the effect of body force, is derived. Secondly, a unified method is developed to obtain the stress function. The analytical expressions of axial force, bending moment, shear force and displacements are then deduced through integration. Thirdly, the stress function is employed to solve problems of anisotropic functionally graded plane beams, with the integral constants completely determined from boundary conditions. A series of elasticity solutions are thus obtained, including the solution for beams under tension and pure bending, the solution for cantilever beams subjected to shear force applied at the free end, the solution for cantilever beams or simply supported beams subjected to uniform load, the solution for fixed–fixed beams subjected to uniform load, and the one for beams subjected to body force, etc. These solutions can be easily degenerated into the elasticity solutions for homogeneous beams. Some of them are absolutely new to literature, and some coincide with the available solutions. It is also found that there are certain errors in several available solutions. A numerical example is finally presented to show the effect of material inhomogeneity on the elastic field in a functionally graded anisotropic cantilever beam.     

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