黏弹黏塑性有限元法应力更新隐式算法-计算力学学报.PDF

卷 期 棳 棾椀 椀 计算 力 学 学 报 斨旓旍棶棾椀 斘旓棶椀 年 月 旇旈旑斿旙斿斒旓旛旘旑斸旍旓旀 旓旐旛旚斸旚旈旓旑斸旍斖斿斻旇斸旑旈斻旙 斚斻旚旓斺斿旘棽棸棻椄 棽棸棻椄 棻棸 旔 椇 闫富有 斈斚斏棻棸棶椃椀棻棻 旙旍旞棽棸棻椃棸椆棽椆棸棸棻 旉 黏弹 黏塑性有限 法应力更新隐式算法 闫富有灣 崔 昊 张晓婉 刘忠玉 棬 棳 棭 郑州大学 土木工程学院 郑州 棿椀棸棸棸棻 椇 暍 棳 摘 要 黏弹 黏塑性耦合模型的黏弹性部分由弹簧 黏壶和 链串联而成 黏塑性部分为双曲线型 灢 斔斿旍旜旈旑 斈旘旛斻旊斿旘灢 暍 暎 棳 斝旘斸斿旘屈服函数 各向同性硬化和斝斿旘旡旑斸黏塑性流动模型 基于黏弹性蠕变柔度 通过定义与弹性问题相对应 旂 旟 棳 暎 棳 的与时间增量相关的黏弹性剪切模量和体积模量 导出增量递推形式的本构方程 为保证算法的收敛和稳定性 棳 把斝斿旘旡旑斸黏塑性流动方程转化为与弹塑性相似的一致性条件 建立黏塑性增量因子单侧逼近其收敛值的 斘灢斠 旟 暎 棳 暎 暍 迭代算法 最后 给出应力更新完全隐式算法和最终计算公式 分别采用黏弹性 黏弹 塑性和黏弹 黏塑性本构 灢 灢 棳 棳 暎 关系对一地基蠕变模型进行三 有限元分析和比较 结果表明 本文算法具有较高的计算效率和稳定性 椇 椈 椈 椈 椈 关键词 黏弹性 斝斿旘旡旑斸黏塑性 应力更新算法 一致切线算子 有限元法 旟 中图分类号椇 椈 文献标志码椇 文章编号椇 棬 棭 斣斦棿椃棸斚棾棿椀 斄 棻棸棸椃灢棿椃棸椄棽棸棻椄棸椀灢棸椂棻棻灢棸椄 棳 黏壶串联而成的黏弹性模型 其松弛模量并非显 棻 引 言 棳 棳 式 复杂模型松弛模量的推演也较难 导致计算程