应力主轴旋转试验离散元细观分析.docxVIP

应力主轴旋转试验离散元细观分析 1 应力路径对粗粒料力学性能的影响 由于粗粒体尺寸大，等级组合理论，抗弯性好，密度高，变形小，具有海绵技术的特点。随着我国西部大开发战略的实施以及高坝水库建设的增多,对于粗粒料力学特性的研究亟待深入,但由于试验条件的限制,对粗粒料应力–应变、变形和强度特性的研究远远落后于细料土。通常认为应力路径对粗粒料强度特性的影响很小,但是不同应力路径对粗粒料力学特性影响的研究大多基于室内三轴试验,而更复杂应力路径(如扭剪、应力主轴旋转等)的试验研究较少,其细观机制更加不明确。 空心圆柱扭剪仪(hollow cylinder apparatus,简称HCA)是目前可以实现包括应力主轴旋转在内的多种复杂应力路径的最先进的土工试验设备。本文利用同济大学GDS空心圆柱扭剪仪(TJ–1)对粗粒料在三轴压缩、扭剪、三轴–扭剪联合、应力主轴旋转等应力路径下的力学特性展开研究,总结不同应力路径对粗粒料强度和变形特性的影响规律。并利用三维离散元程序(particle flow code in three dimensions,简称PFC3D)对试验过程进行数值模拟,通过对颗粒运动过程的监测,进一步揭示应力路径对粗粒料力学特性影响的细观机制。 2 应力场及参数 同济大学GDS空心圆柱扭剪仪可制备试样的最大尺寸为60 mm(内径)×100 mm(外径)×200 mm(高),试样厚度为20 mm。根据相关研究成果以及试验规程,确定试验用料的允许最大粒径为4 mm,属于砾石范畴。鉴于目前国内外尚未研制成功大型空心圆柱扭剪仪,因此对于粗粒料在复杂应力路径下力学特性的研究也只能被限制在砾石土的范围。试验所用砾石料由上海市建筑科学研究院提供,试样颗分级配曲线见图1。试样按设备可制备的最大尺寸制样,最大粒径为4 mm,制样干密度为1.92g/cm3,相对密度为0.8。 HCA工作的基本原理是通过在压力室内安设空心圆柱状试样,对其施加轴力W,扭矩MT,内围压iP及外围压oP,图2给出了空心圆柱试样单元体应力状态。相关应力、应变以及应力主轴旋转角的公式推导详见D.W.Hight等的研究成果,在此仅列出加载参数与应力参数以及应力主轴转角之间的定量关系: 本文采用HCA对砾石土进行了三轴压缩、扭剪、三轴–扭剪联合、应力主轴旋转等不同应力路径的试验。为了对比不同应力路径的影响,反映砾石土在复杂应力状态下的受剪程度,需采用统一的强度标准。本文采用广义剪应力τg和广义剪应变γg作为判别砾石土强度的控制标准,广义剪应力和广义剪应变分别按下式计算: 式中:σ1,σ2,σ3和ε1,ε2,ε3分别为大、中、小主应力和大、中、小主应变。 3 应力主轴旋转的试验过程 由于是初次尝试利用HCA研究不同应力路径对砾石土力学特性的影响,试验过程中保持内、外围压不变且相等,维持在200 k Pa,并且本文仅研究静力加载的情况。 (1)三轴压缩试验:试验采用位移控制,待等向压力变形稳定之后施加轴向位移,轴向剪切速率为0.2 mm/min,试样受剪破坏或轴向应变达到15%时结束试验。 (2)扭剪试验:试验采用位移控制,待等向压力变形稳定之后施加环向扭剪位移,扭转剪切速率为0.003 rad/min,试样受剪破坏或剪应变达到10%时结束试验。 (3)三轴–扭剪联合试验:试验采用应力控制,待等向压力变形稳定之后施加轴向应力σz和环向应力τzθ,试验过程中保持扭剪应力τzθ与(σz-σθ)/2相等,即在τzθ-(σz-σθ)/2应力路径平面上加载曲线为45°向上的斜直线,如图3所示。扭剪应力τzθ的加载速率为1.5 k Pa/min,直至试样破坏。 (4)应力主轴旋转应力路径试验:笔者通过反复试验发现,目前利用HCA研究应力主轴连续旋转的最大限制在于只能在低应力水平下进行,而在高应力水平下则无法实现,如图4所示。在τzθ-(σz-σθ)/2平面上,水平坐标轴正向至曲线上的点与原点连线的夹角为应力主轴旋转角α的2倍(见式(8))。在高应力水平下,当应力主轴旋转至左半平面,会导致试样瞬间破坏,扭剪变形剧烈增大。其原因在于,σθ受内、外围压的制约(式(3)),其值只能在一定范围内变化,而如果要使(σz-σθ)/2足够小,试样所受轴向应力σz必须减小,甚至为负数(即试样轴向受拉),在过低的轴向力约束下试样承受扭剪荷载,必然发生瞬间破坏。 为了深化对应力主轴旋转应力路径的认识,对比高、低应力水平下应力主轴旋转对粗粒料力学特性的影响,本文设计了全新的加载路径,即应力主轴半周旋转分级加载路径,如图5所示。该应力路径在每一级荷载下可实现应力主轴在-45°~45°范围内往复旋转,待旋转一定的周次后,可按照某一预设路径加载到下一级荷载,再次进行主轴旋转,直至试样破坏。该路径虽然无法进行应力主轴在-90