《干砂强夯动力颗粒定向细观颗粒流探讨》由会员分享，可在线阅读，更多相关《干砂强夯动力颗粒定向细观颗粒流探讨（17页珍藏版）》请在金锄头文库上搜索。
1、从本学科出发，应着重选对国民经济具有一定实用价值和理论意义的课题。课题具有先进性，便于研究生提出新见解，特别是博士生必须有创新性的成果干砂强夯动力颗粒定向细观颗粒流探讨摘要:通过引进和开发二维颗粒流程序,基于相似理论建立了可以模拟砂土地基强夯加固的细观颗粒流模型,结合小比尺室内细观模型试验,从颗粒细观力学角度入手对干砂在强夯冲击加固过程中的颗粒定向特性进行了数值模拟。结果表明,数值模拟结果与试验结果具有较好的一致性,实现从细观角度揭示干砂强夯动力微观加固机理,研究工作为今后砂土强夯加固宏细观机理研究提供了一条新的思路。关键词:强夯;干砂;颗粒流;颗粒定向1引言强夯法已成为软弱地基特别是砂性土地
2、基最为常用的加固措施之一。强夯加固地基就是利用巨大的冲击能量,在地基中产生极大的冲击波和动应力,使地基土体产生密实或动力固结,从而达到改善地基工程特性的目的1-7。但至今强夯加固法的研究仍然处于半理论、半经验状态之中,理论严重滞后于实践是不容质疑的事实。目前,学术界普遍认为强夯加固存在三种理论:(一)动力固结理论;(二)波动理论;(三)微观结构理论。到目前为止,传统解释强夯加固机理的“微观理论观点”仍然是人们基于一系列理论方法推导假设而来,在学术界和工程界普遍被人们认为是正确的,但至今也没有很好的方法从微观的角度上来阐述和证明它,依据很多参考文献和工程实例,我们知道在冲击荷载作用下,位于不同位
3、置不同区域的土体颗粒长轴定向性规律可能会有所不同。但就目前已有的室内模型试验装置和数码拍照摄像技术的水平,还不可能完全知晓模型箱内任意一个区域的颗粒长轴定向性的分布规律。然PFC颗粒流(Particle Flow Code in Dimension)8却给我们提供了一个新的研究手段,通过PFC颗粒流程序的二次开发,利用FISH语言统计每次冲击后不同区域的颗粒数目和颗粒长轴的定向性分布规律,我们能够研究数值模型中任意一个区域的颗粒长轴定向性分布规律。因此,PFC颗粒流在研究这个问题上具有其它数值模拟方法不可比拟的优越性,研究的成果有助于我们进一步加深对“微观理论观点”理解和认识。本文通过引进和开
4、发颗粒流程序,采用PFC2D内置的Fishtank函数库和FISH语言,基于散体颗粒力学特性建立了可以模拟砂土地基强夯加固过程的细观颗粒流模型,对冲击荷载作用下锤底及周边影响范围内的颗粒定向规律进行了系统的模拟和验证,研究结果有助于正确理解冲击过程中砂土颗粒细观力学响应和宏观密实间的关系。2颗粒流数值模型的建立颗粒接触模型在建立砂土强夯细观颗粒流模型时,砂土颗粒的接触模型对模拟结果影响较大。为了有效模拟砂土在冲击荷载作用下的动力反映特性,笔者比较了3种接触模型(线弹性、Hertz和滞回阻尼)模拟结果的优缺点,最终选用了滞回阻尼接触模型(Hysteretic Damping Model)。图1为
5、滞回阻尼模型法向和切向本构模型示意图,从图中可以看出,在受到外荷载作用时颗粒法向刚度在加载阶段和卸载阶段不同,卸载刚度明显大于加载刚度。该模型这个特性能够很好地反映砂土在强夯过程中的刚度变化和瞬间冲击密实现象。模型箱结合实际强夯的对称性和影响范围,为加快计算速度,在数值建模时采用半模数值模型(见图2)。模型箱尺寸为60 cm60 cm,是根据实际强夯影响范围和相似性原理按照一定比例缩小后得到的。模型箱由墙体形成,共定义左、右和下部3道墙体。为了能够使半模的结果与全模一致,通过多次调整参数试算,确定左侧墙体的法向刚度1101N/m,切向刚度为0,墙壁光滑,其余两边的墙体参数见表1。颗粒及其细观力
6、学参数为减少计算量,在数值建模时采用分块建模方式,如图2所示。在冲击荷载主要影响区域1采用较小的颗粒粒径,在b(区域1宽度)h(区域1长度)=的范围内生成20362个颗粒,而在影响较弱的2、3区域采用较大的颗粒粒径,共生成18814个颗粒。数值砂颗粒的细观力学参数是基于双轴试验结果通过多次试算和调整后得到的9,具体见表1。夯锤及其细观力学参数通过将84个直径为1cm的颗粒球叠合在一起构成CLUMP块来模拟实际夯锤,夯锤尺寸为bh=4cm,如图3所示。夯锤细观力学参数见表1。颗粒形状的选择关于颗粒形状对砂土力学性质的影响一直是土力学试验研究的一项重要内容。PFC2D中虽然不能直接生成椭圆颗粒或多
7、边形颗粒,但是提供了“Clump”方法(“团颗粒”方法)以供非圆颗粒单元的二次开发,该方法将若干个团聚在一起的圆盘颗粒视为单一团颗粒,根据团聚形式的不同可以构造任意形状的团颗粒,圆盘颗粒之间可相互重叠也可不重叠。根据面积和质量等效原则,开发了近似椭圆团颗粒的构造,见图4所示。3数值模拟及结果分析数值模拟步骤与室内模型试验一样,数值模拟也可分制样和冲击加荷两个过程,强夯加固的颗粒流数值模拟可分以下几步:建立模型箱和在预定范围内生成相应大小的砂颗粒,赋予其合适的细观力学参数,详见表1,施加重力场,让其在自重的作用下逐渐沉积,直到系统达到平衡稳定状态,待纯圆颗粒模型生成并在自重作用下沉积并得到平衡后
8、,利用面积和质量等效原则,用长细比Se=的近似椭圆颗粒替换所有的纯圆颗粒,得到椭圆颗粒数值试样,在自重作用下平衡,如图6所示。然后在颗粒间设置滞回阻尼接触模型并运行足够多的时步,直至系统再次达到平衡状态。运用Clump Logic生成半个重锤,重锤共有84个直径为1 cm的颗粒球组成,如图3所示。为了节省计算时间,在保证动能和动量相等的前提下,重锤自由下落阶段通过在试样表面赋予重锤一个初始速度 代替。当落高为H=40cm时,重锤直接在试样表面生成后赋予初始速度为/s。第1击冲击模拟开始后时刻监测重锤的速度和整个系统的动能,当重锤的速度基本保持0不变和系统动能最小时,可以停止计算,认为第一击完成
9、。第1击完成后,删除重锤,让整个系统在自重的作用下再次达到静力平衡状态后,重新在第一击的同一位置生成重锤,赋予相同初始速度再次冲击土体,当重锤的速度基本保持0不变和系统动能最小时,可以停止计算,认为第2击完成。重复操作、,逐步完成所需要的N击。 3.颗粒长轴定向性分布规律 在进行颗粒长轴定向性分布规律分析时,根据强夯的一般影响范围对冲击荷载影响较为强烈的区域进行统计,为了能够与室内模型试验结果进行对比,首先对锤底1D2D(D为锤径,下同)范围内的颗粒长轴与水平方向的夹角进行统计,其次,我们对锤底2D3D,以及3D4D范围内的颗粒也进行统计,共设置8个分析区域,统计区域的具体位置见图6,图7为局
10、部颗粒放大图。 区域1统计结果分析图8为不同击数下统计区域1内的颗粒数目,可以看出经过6次冲击后,落入该区域的颗粒数目从468个增加到550个,颗粒数目增加了82个,表明该区域(深度范围在1D2D)之间的土体在冲击荷载的作用下随着击数的增加不断地被密实,颗粒之间的孔隙明显减小,土体得到了明显的加固。图9分别为冲击前和冲击六次后的颗粒长轴定向性演化玫瑰花图,可以看出,冲击前位于010和170180这两个区间的颗粒数目之和为37个,一击后位于这两个区间的颗粒数目之和为96个,落入90160之间的颗粒数目在第一次冲击后颗粒明显减小。随着击数的增加,颗粒长轴落入010和170180的颗粒数目不断增加,
11、到6击后,这两个区间的颗粒数目之和达到了133个,比冲击前增加了96个,而落入70160之间的颗粒数目仍然在不断减少。表明该区域的颗粒在冲击荷载的作用下颗粒长轴与水平方向的夹角逐渐趋于零,即颗粒的长轴趋于水平向分布,颗粒的长轴定向性分布规律相当明显。区域2统计结果分析图10为不同击数下统计区域2内的颗粒数目,经过6次冲击后,落入该区域的颗粒数目从462个增加到507个,颗粒数目增加了45个,表明该区域(深度范围在2D3D之间)的土体在冲击荷载的作用下随着击数的增加也在不断地被压密,颗粒之间的孔隙明显减小,土体得到了明显的加固,但是相对区域1来说,该区域在6击后颗粒增加的数目小于区域1的,即该区
12、域的加固程度要稍弱于区域1的。下面的定向性规律更能说明这一点。图11为冲击前和冲击六次后的颗粒长轴定向性演化玫瑰花图,可以看出,冲击前位于010和170180这两个区间的颗粒数目之和为50个,一击后位于这两个区间的颗粒数目之和为73个,这一增加量的数目是小于区域1的,随着冲击次数的增加,落入010和170180这两个区间的颗粒数目也在不断增加,落入90160区间的颗粒数目在不断减少,到第六击后,落入010和170180区间的颗粒总数目为105个,相比冲击前增加了55个(这一数目也小于区域1),表明该区域(2D3D)土体在冲击荷载作用前后颗粒的长轴仍然出现了明显的定向分布规律,但定向性程度较区域
13、1要稍弱一些。区域3统计结果分析图12为不同击数下区域3内的颗粒数目变化情况,可以看出在冲击荷载的作用下颗粒数目发生了明显的增加,其原因是该区域位于3D4D之间,该区域跨越粒径差别较大的两种颗粒,详见图所示,所以该区域内的初始状态的颗粒数目明显小于区域1和区域2的。在冲击荷载的作用下,位于该区域内的颗粒会向下运动,导致该区域下部的大粒径颗粒不断移出该区域,其上面的小粒径颗粒不断从该区域的上部移入,这一过程的发生造成了该区域的颗粒数目发生了明显的增加。图13为冲击前和冲击六次后的颗粒长轴定向性演化玫瑰花图,可以看出,冲击前位于010和170180这两个区间的颗粒数目之和为35个,一击后位于这两个
14、区间的颗粒数目之和为32个,反而有所减小,但随着冲击次数的增加,到第六击后,落入010和170180区间的颗粒总数目达到57个,相比冲击前增加了22个,表明该区域(3D4D)的颗粒长轴仍然还存在定向性分布规律,但是在冲击荷载作用前后玫瑰花图的大体形状变化不是很大,颗粒长轴与水平方向夹角的分布规律也基本上没有太明显的变化。与区域2相比,该区域的砂土颗粒的长轴定向性排列程度又要弱一些。区域4统计结果分析当竖向深度达到4D5D的区域4后,在冲击荷载的作用下,落入该区域的颗粒数目已经基本上不再增加了,经过6次冲击后,该区域的颗粒数目只增加了5个,因此基本上可以认为冲击前后该区域颗粒数目基本保持恒定,这
15、一规律表明:冲击荷载作用到达该深度后对土体的影响已经不明显。图15分别为冲击前和冲击六次后的颗粒长轴定向性演化玫瑰花图,可以看出,随着冲击次数的增加,颗粒长轴落入任意一角度区间的颗粒数目基本保持不变。表明该深度区域受到冲击荷载的影响很弱。区域5统计结果分析图16为不同击数下区域5内的颗粒数目变化情况,可以看出在冲击荷载的作用下颗粒数目发生有少许增加,第一击后落入该区域数目仅增加了5个,到六击后颗粒数目相对与冲击前也只是增加了23个,随着冲击次数增加颗粒增加的数目有略微的波动。图17分别为冲击前和冲击六次后的颗粒长轴定向性演化玫瑰花图,可以看出,随着击数的增加,颗粒长轴方向落入010和170180的颗粒数目基本保持不变,长轴定向性玫瑰花图的在不同冲击次数下的形状也大体相同,通过6次冲击后,颗粒长轴方向落入010和170180的颗粒总数目为124个,这个数目比冲击前的102个增加了22个,颗粒长轴方向落入010和170180这两个区间的颗粒总数目时而增加时而减少,呈现出不稳定状态,基本上不存在较为明显的定向
