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

第 30 卷 第4 期 应 用 力 学 学 报 Vol.30 No.4 2013 年 8 月 CHINESE JOURNAL OF APPLIED MECHANICS Aug. 2013 文章编号：1000- 4939(2013) 04-0581-06 功能梯度矩形板问题三维理论的半解析解 杨智勇 1,2 牛忠荣 1 程长征 1 伍章健 3 (1 合肥工业大学工程力学系 230009 合肥；2 铜陵学院机械工程系 244000 铜陵； 3 曼彻斯特大学机械、航空、土木工程学院 M13 9PL 曼彻斯特) 摘要：针对功能梯度材料矩形板问题，基于三维弹性理论，将位移和应力分量作为基本变量，通 过双三角级数将其控制微分方程转化为常微分方程组的边值问题。采用插值矩阵法直接对常微分 方程组边值问题进行求解，得到了功能梯度材料矩形板三维位移、应力场的半解析解。通过算例 给出了材料参数按指数形式和幂函数形式变化情况下的功能梯度板的弯曲问题。对比有限元法和 状态空间法，结果表明：本文提出的状态空间与插值矩阵法相结合的半解析法能有效地分析材料 -4 参数按任意形式连续变化的功能梯度矩形板问题，且具有良好的精度，精度可达 10 量级，能够 满足工程需要；与其他方法相比，本文方法具有实施便捷、计算量小等优点，根据其力学场分析 结果可设计出满足各种不同需求的功能梯度材料。 关键词：矩形板；功能梯度材料(FGM)；三维；弹性理论；插值矩阵法 中图分类号：TB12 文献标识码：A DOI ：10.11776/cjam.30.04.C033 条件、各种载荷，缺点是进行三维求解时的计算量 1 引 言 较大。文献[2]研究了功能梯度圆柱壳在均匀内压作 用下的线弹性弯曲问题。文献[3]采用直接位移法分 多材料组合结构在结合界面边角处会产生应 析了功能梯度圆板的轴对称自由振动问题。文献[4] 力集中、奇异应力，并常常成为结构破坏的首要薄 采用无网格Petrov-Galerkin法对FGM矩形板进行了 弱处，特别是弹性模量差异较大的结合材料。构造 三维静力分析。文献[5]采用高阶剪切应变理论分析 材料参数，设计功能梯度变化的特性，以连续变化 了功能梯度板的静态响应及自由振动。文献[6]采用 材料参数避免结合结构应力集中，这是优化材料设 打靶法对FGM矩形板进行了三维静动态响应分析。 计的有效途径。因此非均质材料结构的力学场是一 在现有解法中，状态空间法的计算效率较高， 项有待于研究的新课题。 在许多板问题求解中得到应用。文献[7]针对三维弹 文献[1]采用有限单元法分析了横向加载下的 性板壳提出了状态空间法，对矩形板形状采用双三 功能梯度板问题，材料属性是沿厚度方向分层定义 角级数展开；通过求出状态方程的特征值、特征向 的；为保障计算精度，划分了较多的层数，随着层 量，给出了状态转移矩阵，从而得到了三维弹性板 数的增加，对长宽方向也作了进一步细分，使得单 各位移、应力分量的解析解。然而对于功能梯度矩 元数目剧增。有限元方法的优点在于适宜各种边界 形板，所得的是变系数常微分方程组，其状态转移 基金项目：国家自然科学基金(11