热环境中黏弹性功能梯度材料和其结构的蠕变.pdfVIP

第 44 卷第 2 期 力 学 学 报 Vol. 44, No. 2 2012 年 3 月 Chinese Journal of Theoretical and Applied Mechanics Mar., 2012 热环境中黏弹性功能梯度材料及其结构的蠕变1) 彭 凡 顾勇军 马庆镇 (湖南大学机械与运载工程学院, 长沙 410082) 摘要 基于经典的对应原理, 将Mori-Tanaka 方法等细观力学结果推广于定常温度环境下的黏弹性情形. 根据 泊松比与时间呈弱相关的特点, 给出Laplace 象空间中功能梯度材料的松弛模量和热膨胀系数, 并直接建立耦 合热应变的多维黏弹性本构关系. 在此基础上, 求解黏弹性功能梯度圆柱薄壳在热环境中的轴对称弯曲蠕变变 形问题. 考虑材料热物参数的温度相关性, 首先确定稳态温度场, 导出相空间中轴对称弯曲变形的解析解, 采用 数值反演得到蠕变变形. 算例表明, 蠕变初期, 热环境的影响明显, 随着时间增加, 热应力松弛, 影响逐渐消失. 当圆柱薄壳受轴压时, 相比于两端固支, 两端简支的端部变形更加明显. 通过圆柱薄壳的轴对称弯曲求解, 给出 体积含量呈任意分布的黏弹性功能梯度结构在热机载荷下的蠕变分析途径. 关键词 功能梯度材料, 热环境, 黏弹性, 圆柱薄壳, 蠕变, 细观力学 中图分类号: O034 文献标识码: A 文章编号: 0459-1879(2012)02-0308-09 引 言 征, 这是功能梯度材料及其结构的力学分析的主要 思想 . 梁军等 基于对应原理和Eshelby 等效夹 热环境中, 金属/陶瓷和聚合物/陶瓷等功能梯 杂方法预测黏弹性颗粒复合材料的松弛与蠕变特性. 度材料表现出时间相依的力学行为, 研究其蠕变变 刘书田等 给出了基于均匀化理论分析复合材料黏 形、蠕变失稳及延迟开裂机制等问题是十分必要的. 弹性的多尺度数值方法, 在相空间中进行有限元离 - 散求解单胞域上的细观均匀化问题以及宏观位移和 文献[1 2] 基于Norton 蠕变模型, 假设梯度材料的蠕 变系数为径向坐标的函数, 分别得到了受内压球罐 宏观应力, 结合曲线拟合和Laplace 逆变换得到材料 的球对称变形解和受内压厚壁圆筒的轴对称变形解. 的黏弹性参量与热应力松弛模量, 分析过程较复杂. 一些学者利用黏弹性力学模型描述功能梯度材料的 Muliana 等- 提出了黏弹性功能梯度结构的双尺 度分析方法, 采用有限元法求解宏观场, 在积分节点 热流变, Paulino 等 和Khazanovich 提出弹性-- 黏 处, 构造细观场并求解弹性、黏弹性应变, 给出宏细 弹性对应原理; Zhang 等 分析了功能梯度薄板的 观耦合求解算法, 这是一种纯数值求解途径. 目前, 热黏弹性弯曲问题; 李伟杰等 利用Paulino 的对应 原理研究了黏弹性功能梯度材料中垂直梯度方向的 关于功能梯度材料及其结构的许多热弹性问题得到 界面裂纹问题, 得到应力强度因子随时间变化的特 了解析或半解析解, 如何在考虑组分材料体积含量 征. 以上工作都是针对某类特殊的功能梯度材料展 呈任意分布以及热流变特性的前提下, 进一步给