一种改进修正剑桥模型计算精度算法.docVIP

一种改进修正剑桥模型计算精度算法 　　[摘 要] Roscoe和Burland将剑桥模型的屈服面设定为椭圆，得出修正剑桥模型，修正剑桥模型充分利用等向压缩试验曲线，同时只利用了三轴试验中临界状态线，对三轴实验成果利用较少。本文通过运用修正剑桥模型从理论上计算土体的体积应变和剪切应变，运用CU三轴试验测出相应土体的体积应变和剪切应变，利用最小二乘法理论拟合理论计算和试验结果的比值，对修正剑桥模型进行改进，得到的弹塑性矩阵改善模型的计算精度。 　　[关键词] 修正剑桥模型；三轴试验；曲线拟合；计算精度； 　　中图分类号：TU43 文献标识码：A 　　引言 　　修正剑桥模型是英国剑桥大学的Roscoe和Burland[1]根据正常固结粘性土和弱超固结粘性土的三轴试验，采用状态边界面的概念，由塑性理论的流动法则和塑性位势理论，利用简单曲线配合法，建立塑性与硬化定律的函数。提出临界状态线、状态边界面、弹性墙等一系列物理概念，开创了土力学临界状态理论。 　　经历半个多世纪的发展，国内外许多专家学者对修正剑桥模型进行了一系列改进，如基于经典苏醒理论框架的修正、非关联流动法则、次塑性理论与临界状态理论、有限应变、时间有关、各向异性及结构性以及循环荷载下的临界状态[2]等。 　　在实际工作中，人们发现修正剑桥模型充分利用等向压缩试验曲线，如：土性参数λ、K以及压缩指数Ce等[3]，而对三轴试验，则只利用了临界状态线[4]，对三轴试验成果利用较少。 　　本文通过用修正剑桥模型计算土体的体积应变和剪应变，根据对应三轴试验结果，拟合两者比值得出修正系数，代入土的本构模型增量关系式中形成新的弹塑性矩阵（刚度矩阵或柔度矩阵），从而达到改进模型计算精度的目的，既验证了剑桥模型的适用性，又充分利用了三轴试验的实验成果。 　　1 修正剑桥本构模型基本理论 　　对修正剑桥模型有 　　2 修正剑桥模型的改进 　　2.1 理论计算土样参数 　　修正剑桥模型是Roscoe和Schofield等于1968年在三轴试验的基础上提出的用于正常固结土或弱超固结土的模型，该模型于1974年得到MIT路堤试验结果的验证[7]，至今仍广泛应用。 　　本文基于太沙基理论运用修正剑桥模型计算出三轴试验（CU）中土体的体应变和剪应变，在CU试验中，由于土样已经固结完全，在三轴剪切过程中土样已属正常固结土，其应力路径如图1。 　　在土力学中，土介质的不均匀性、各向异性导致其应力路径也是不均匀和各向异性的[8]，即一种模型在不同加卸荷载作用下，其路径是不同的。故理论计算时采用“微分”的观念，对土中的有效应力计算按时步划分，确保计算精度，确定每一计算时步有效应力的计算增量为： 　　至此：可得到塑性体积应变增量、塑性剪应变增量计算式。 　　2.2 三轴试验计算土样参数 　　可求出广义剪应变，配合每个时步的应力―应变关系的变化，求出体积应变增量以及剪应力增量[11]。 　　2.3 修正剑桥模型改进方法 　　2.4 修正系数αβ的曲线拟合 　　由于利用最小二乘法只是对数值上进行拟合，故不会影响土的本构模型数学基础，故由此得到（α，β）关系式代入式（10），得到新的弹塑性刚度矩阵，可改善现有修正剑桥模型的计算精度。 　　3 算例及分析 　　图2给出了在给定应力路径情况下改进后模型与改进前模型的q-ε1结果比较，从图中可以看出：修正后的模型与修正前都能比较好的反映这一应力路径下的应力―应变规律，但是修正后模型计算的轴向应变比修正前计算结果小，即修正后的模型有效的提高了计算精度。 　　4 结论 　　（1）本文基于太沙基三维固结理论，利用修正剑桥模型在理论上推导CU三轴试验中土体的塑性体积应变与塑性剪应变。 　　（2）利用CU三轴试验直接测出试验中土体的体积应变，并通过广义剪应变计算式得出土体的剪应变。 　　（3）利用最小二乘法对理论推导与实测体积应变和剪应变进行曲线拟合，得出修正剑桥模型的修正系数关系式可充分利用实验数据对修正剑桥模型进行修正。 　　（4）根据拟合结果得到的弹塑性矩阵可提高已有模型的计算精度或者由此构造新的修正形式。 　　[参考文献] 　　[1] 谢定义，姚仰平，党发宁.高等土力学[M].北京：高等教育出版社，2008.109-164. 　　[2] 徐舜华，徐光黎，程瑶. 土的剑桥模型发展综述[J].长江科技学院院报，2007，24（3）：27-33. 　　[3] 张云.修正剑桥模型参数对计算结果的影响[J].岩土力学，2006，27（3）：441-445（Zhang yun. Effects of parameters of modified Cambridge model on computed results[J]. Rock and Soid Mechanics，2006，27（