832吉林大学学报(地球科学版第38卷0引言随着我国经济建设和交通运输事业的发展,车流量、行车速度和载重量不断增加,特别是随着高速铁路建设,如京沪高速等在我国的大规模开展,循环施加的交通荷载[1刮对路基受力变形特性的影响也越来越大[3刮;而我国东部沿海地区广泛分布着高含水量的粘性土,路基的沉降病害问题十分突出,如高速公路路面因路基下沉而致开裂。因此,对粘性土受循环荷载作用下的动力特性问题进行研究具有十分重要的意义。土体动三轴试验可以很好地模拟动荷载的施加过程及土体现场实际的三向受力状态,因此试验结果更能反映土体实际变形状况。但动三轴试验本身从制件到荷载试验,需耗费大量物力和时间,故本文尝试采用有限元数值模拟来替代部分实际试验,以节省研究的时间和费用,推动土动力学研究以更快的速度发展[6。7J。l循环荷载下粘性土弹塑性模型迄今为止,研究土体动荷载作用下变形特性的本构模型有三大类,即边界面模型、套叠屈服面模型及基于广义塑性理论的弹塑性模型"剖。本文在有限元分析中使用基于Mohr—Coulomb屈服准则的弹塑性模型。2动三轴试验2.1试验土样实验采用长春市卫星路取土场的中液限粉质粘土,取样深度2.5m。其基本物理力学参数如下:密度为1.96g/cm3,最佳含水量为21%,液限为34%,塑性指数为11.6,内聚力(f为47kPa,内摩擦角(垆为34。。2.2动三轴试验.’试验土样为直径70mm、高度140mm的圆柱体击实成型试件(图1。试样真空抽气饱和后等向固结,在不排水条件下先施加大小为0.5qa静荷载,然后施加正弦波荷载,振幅0.5qa,频率为1Hz。其中q。为试验中的竖向荷载,单位kPa。仪器为日本DTC一306型多功能电液伺服动态三轴仪,采用应力控制方式。一台WX2400X—Y绘图仪用于记录试验数据,选用应变5%为破坏标准。图1动三轴试验现场及土样Fig.1Sceneandsoilspecimenofthedynamictri-axialtest2.3滞回曲线图2为实验所得的不同动应力和围压条件下部分动应力和位移关系图。从图中可以看出,土在往复荷载作用下会因土粒相互滑移形成新的排列而产生不可恢复的永久塑性变形。3有限元计算分析3.1三维有限元模型的建立用ABAQUS软件建立的三维有限元模型L9。10]如图3所示。模型中定义的土材料属性:密度为1.96g/cm3,杨氏模量为18MPa,泊松比为0.35,渗透系数为1.0×10~mm/s,孔隙率为0.3,内摩擦角34。,膨胀角0。,固结不排水。根据土的参数及ABAQUS软件有限元模型适用条件,采用Mohr—Coulomb模型进行数值模拟。所施加荷载为正弦波形且考虑土样自重,频率为1Hz,即q(￡一0.5qd-t-O.5qdsin2nt。其中:q(￡为任意t时刻的轴向动倚载值,见图4。3.2模拟和试验塑性变形的比较与分析由于实验刚开始阶段不稳定,加上有非线性静载段的影响,处理数据复杂而且不准确,所以将第二次加载作为处理和分析数据的起点,每一次的塑性变形从第二次加载开始进行累加,将实验的数值和有限元模拟的数值进行对比。本文取出了4组数据进行分析比较,如图5。从图5可以看出,随着循环次数的增加,累积塑性变形逐渐增加,实验和模拟的每一组数据中,塑性变形随循环次数n的变化趋势非常一致,而且每一第5期戴文亭,陈星,张弘强:粘性土的动力特性实验及数值模拟833(aol-480kPa,o3=100kPa,9周破坏(bol=580kPa,o3=200kPa,15周破坏(c01—640kPa,03=300kPa,23周破坏(doI=680kPa,o,=300kPa,7周破坏图2不同动应力和围压条件下动应力和位移关系图Fig.2Dynamicstcess-displacementgraphunderdifferentdynamicloadamplitudeandconfiningpressure图3三维有限元模型Fig.3Thethree-dimensionalfiniteelementmodelqd量萼和/2O0.250.500.751.001.25t/s图4循环荷载谱Fig.4Cyclicloadingspectrum次循环的塑性变形和累计塑性变形都吻合很好。经过统计分析,累积塑性变形模拟值与实验值问的误差最大为16.30%,在有限元允许误差20%以内,最小的仅为0.19%,很好地满足了精度要求。这些都充分说明所建立的有限元模型是正确的,模型中选取的参数是合理的,也表明用模拟的方法来替代动三轴实验是可行的,该方法将极大地节省研究的时间和各种费用,在土动力学分析研究中有推广应用的价值。实验和模拟之间的误差主要有3个方面的原因:一是仪器原因,如记录数据采用的是WX2400X—Y绘...