亚塑性本构模型隐式数值积分算法.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 亚塑性本构模型的隐式数值积分算法 岑威钧 河海大学水利水电工程学院 江苏南京 摘要 针对亚塑性本构方程率型隐式的特点 推导了求解亚塑性本构方程的隐式数值积分算法 并对该算法做了进 一步改进 引入自适应时间步长技术 在保证精度和收敛速度的前提下 通过合理选择时间步长 降低了计算时间 保证了该计算程序的可靠性 通过采用该算法所编制的 亚塑性模型有限元程序 对已知算例进行验 证计算结果表明 本文算法比较精确 最后对一堆石大堤进行了应力变形的亚塑性数值求解 所得结果合理 反映 了该大堤堆石料的主要应力应变特性 表明本文算法求解实际工程问题时可行 关键词 亚塑性 本构模型 率型 模型 隐式数值积分 算法 中图分类号 文献标识码 亚塑性 本构理论 是 世纪中后期由德国学者 等人在有理力学框架体 系内发展起来的主要用于描述散粒型土体应力应变关系的一种全新本构理论 该理论以张量运算符号 为工具 直接建立起联系应力率空间和应变率空间的各向同性非线性张量函数表达式 与经典的弹塑 性理论不同的是 亚塑性理论扬弃了总应变分解为弹性应变和塑性应变 屈服准则 硬化规律 塑性势及 流动法则等基本概念 但这些概念又能在亚塑性本构表达式中得到相应的解释 亚塑性本构理论起初 用于描述砂 粉粒等无黏性散粒型土体的非弹性和非线性等主要力学特性 近年来不断向粗颗粒土 和 黏性土 方向拓广发展 近期发展起来的一些亚塑性模型较以往模型在理论上更加完善 更能全面合 理地反映土体的非弹性 非线性 压硬性 剪胀剪缩性等主要力学特性 并对实际工程问题进行了尝试应 用 由于亚塑性本构方程采用的是率型表达式 且方程的两边同时包含着当前应力状态 不能将应力单 独分离出来进行显式表示 为此 本文推导了一种亚塑性本构方程求解的隐式数值积分算法 并对算法 做了两点改进 编制了相应的亚塑性有限元数值计算程序 用 个算例验证了本算法的合理性和可靠 性 亚塑性本构模型的率型表达式 亚塑性本构模型采用张量符号表示的率型表达式 如 　 式中 为 应力率张量 　 为 应力张量 　为 的时间导数 为旋转 张量 ? 为空间速率梯度张量 为一关于　除　 连续可微的各向同性张量函 数　为应变率张量 为状态变量 是标量 如孔隙比 等 可以没有 一个或多个 收稿日期 基金项目 科技部中奥政府间国际科技合作项目 作者简介 岑威钧 男 浙江慈溪人 博士 讲师 主要从事土石料本构建模分析 土石坝结构计算及水工建筑物渗流分析等 研究 根据各向同性张量函数的一般表述定理 和亚塑性理论率无关性的要求 不失一般性 可 分解成线性部分和非线性部分的组合 即 　 　 式中 为 阶线性算子张量 　 为非线性算子张量 为 范数运算符 如果式 中第二项非线性部分为零 则亚塑性本构模型式 退化为亚弹性模型 从这个意义上 讲 亚塑性本构模型是亚弹性本构模型的推广 从理论上讲 对于不同的土体 散粒体 在满足一些基本本构定律要求和假设的前提下 只要构造 合适的非线性项和线性项函数表达式 就能很好地刻画这种材料的应力应变特性 早期提出的一些亚 塑性模型不涉及任何状态变量 仅为 和　的函数 后来逐渐发展成引入孔隙比 作为状态变量 的改进形式 可考虑应力状态和孔隙比 互相影响 使得模型更能反映无黏性散粒型土体的剪胀剪缩性 等重要力学特性 模型本构参