统一硬化参数模型在有限元中的应用及其验证.pdfVIP

70 岩土力学与工程新进展 统一硬化参数模型在有限元中的应用及验证 甄文战，孙德安 (上海大学土木工程系，上海200072) 摘耍针对能反映砂土剪胀及粘土变形特性的统一硬化参数本构模型，采用弹性预测及塑性修正两阶段的回映 算法，给出了该模型的应力更新表达式及其对应的刚度矩阵，根据其算法，编制了非线性有限元软件子程序，实现 了该模型在有限元软件中的应用，并对该模型进行了数值测试，模拟了三轴压缩、三轴伸长及平面应变等单元试验 问题，通过数值分析结果和试验数据的对比，验证了有限元子程序的正确性，从而为岩土工程中采用该模型进行有 限元分析提供了可能。 关键调砂土，硬化参数；回映算法}有限元子程序 1 引 言 剑桥模型‘13(Cam—clay)是Roscoe等人在对正 (2) 常固结土和弱超固结土进行大量试验的基础上，以 临界状态土力学为理论基础，而提出的一个经典弹 塑性本构模型。剑桥模型以塑性体应变为硬化参 式中，M为体积应变由剪缩发展为剪胀时的应力比，即 数，不能反映砂土这类岩土材料所具有的剪胀特性， 变相时的应力比，M为峰值时的应力比，如图1所示。 为了克服此缺点，姚仰平等人提出了以H作为硬化 等向压缩状态 参数，可以考虑砂土具有的剪胀特性，并给出了该参 (1)当刁=0时， H亍I出： 数的具体表达式[2]，后来通过试验验证了该硬化参 d￡：0 (2)当0r／M时，剪缩状态 数的合理性及应力路径无关等特性r3H]。 (3)当刀=M时，变向状态 d￡：=0 本文通过对H为统一硬化参数模型的分析，采 de：0 (4)当M刁Mf时，剪胀状态 用弹性预测及塑性修正两阶段的回映算法，给出了 (5)当刁=Mf时，临界状态 该模型的应力更新表达式及其对应的刚度矩阵，并 由式(3)可知，当M=Mf时 编制了非线性有限元软件子程序(UMAT)，成功实 现了该模型在有限元软件中的应用，通过对三轴压 H—IdH=Id￡： (4) 缩、三轴伸长及平面应变等单元问题的数值分析结 果和试验数据的对比，验证了子程序的正确性。 则此时硬化参数H就变为以塑性体应变为硬化参 数，即式(1)是式(2)的一种特例。当M≠Mf时，硬 2统一硬化参数 化参数H可以反映砂土的剪胀特性，所以H可作 剑桥模型以塑性体应变为硬化参数，而且采用 为粘土及砂土的硬化参数，故称为统一硬化参数。 相关联流动法则，屈服函数表达式为 f2驴=lAh-,,IFln丢+ln(1+舌)卜：一0 (1) 式中，≯为塑性势函数；A为压缩指数}Ⅳ为回弹指 数；￡：为塑性体应变；可为应力l七(q／p)；加及P分别 为初始平均应力及平均应力。 为了考虑砂具有的剪胀特性，采用H为硬化参 数E利，采用相关联流动法则，屈服函数表达式为