双势积分算法在非关联材料中的应用-应用数学和力学.PDFVIP

　 ， ３９ １ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１８年１月１５ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３９，Ｎｏ．１，Ｊａｎ．１５，２０１８ 文章编号：１０００⁃０８８７（２０１８）０１⁃００１１⁃１８ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 双势积分算法在非关联材料中的应用∗ 周洋靖，　 冯志强，　 彭　 磊 （西南交通大学 力学与工程学院，成都 ６１００３１） （我刊编委冯志强来稿） 摘要：　 在双势理论的框架下，根据材料自由能形式，材料可以被划分为显式标准材料和隐式标准 材料．以经典的非关联Ｄ⁃Ｐ模型为例，对其本构锥体进行了描述，并引入了一对对偶锥体．证明了在 对偶锥体的描述下，不仅能满足非关联Ｄ⁃Ｐ 模型自身本构关系，其应力和塑性应变也能满足隐式 流动表达．结合双势理论和Ｄ⁃Ｐ模型自身的本构特点，推导出了非关联Ｄ⁃Ｐ模型率形式弹性状态 下、率形式塑性状态下、增量形式弹性状态下、增量形式塑性状态下和增量形式弹塑性状态下的双 势函数，从而得到了非关联Ｄ⁃Ｐ模型的双势积分算法．通过数值模拟算例验证了双势积分算法的 准确性和稳定性． 关　 键　 词：　 非关联材料；　 双势理论；　 正交性；　 非关联Ｄ⁃Ｐ模型；　 数值模拟 中图分类号：　 Ｏ３９；Ｏ３４１　 　 　 文献标志码：　 Ａ ＤＯＩ：１０．２１６５６／ １０００⁃０８８７．３８０１３９ 引　 　 言 ［１］ 探究材料本构关系一直是固体力学中一个重要的研究课题 ，能正确模拟材料的应力⁃应 ［２］ 变关系，即数值实现材料的本构积分，是材料非线性研究的核心 ．对于经典的关联性材料，如 ［３］ 大部分金属材料，其理论框架较为成熟，研究方法较为经典，研究成果也较为丰富 ．同样，大 ［４］ 量的研究人员也对非关联材料进行了深入的研究 ．但直到现在，如何在复杂的工程实践中， ［５］ 找到更有效、更准确的模拟方法仍然是一个难题 ． 根据不同的理论假设，不同的数值实现方法，非关联材料模型的本构积分方法存在一定的 差异，文献［６］对非关联模型本构进行了简单的讨论．其中，Ｄｒｕｃｋｅｒ⁃Ｐｒａｇｅｒ模型（以下简称Ｄ⁃Ｐ ［７］ 模型）是一种常用的非关联材料模型 ，常被用于模拟非关联土体材料．能高效准确地模拟简 ［８］ 单工况下非关联Ｄ⁃Ｐ模型的本构关系，是成功模拟复杂非关联材料的基础 ． 针对非关联Ｄ⁃Ｐ模型应力⁃应变关系，传统塑性力学采用显式本构积分算法（ｅｘｐｌｉｃｉｔ ｉｎｔｅ⁃ ［９］ ｇｒａｔｉｏｎ ａｌｇｏｒｉｔｈｍ） ，在每一步迭代计算中，只通过结构满足整体平衡来计算求解材料应力增 量．这种方法的优势在于模型本构积分只依赖于材料的屈服势函数Ｆ 与流动势函数Ｐ 的函数 形式．因此采用显式本构积分算法推导模型本构相对容易．但是由于此方法每一步计算其结构 ［１０］ 只满足整体平衡，可能会造成计算过程的不稳定，计算效率有一