基于离散单元法的碎石集料动力变形特性研究

碎石集料作为铁路道砟以及路基填料等在铁路工程中获得了大量的应用，承受着列车荷载的往复作用，所引起的沉降变形严重影响着行车舒适性、安全性与长期服役性[1]。碎石集料的动变形特性可以通过永久变形和回弹变形两个方面进行描述。其中，回弹变形特性可以采用回弹模量MR评价，其定义为循环荷载作用下，偏应力与可恢复的弹性应变的比值[2]为

(1)

Selig等[3]进行了大型动三轴试验，研究了循环荷载作用下的道砟集料的变形特性，结果表明，轴向应变及回弹模量随着荷载作用次数的增加而增加，并最终趋于稳定，集料在密实的同时持续硬化。Stewart等[4]基于大型动三轴试验，研究了荷载形式对集料永久变形的影响，结果表明：永久变形随幅值的增加明显增加，第一次加载所引起的永久变形最大。Akke等[5]基于室内试验分析了荷载幅值对动变形的影响，认为在动荷载幅值较低时，集料几乎不发生塑性变形。Shenton[6]研究了加载频率对于塑性变形的影响，认为频率对于塑性应变几乎没有影响。Indraratna等[7]基于大型三轴试验的结果认为道砟的累积塑性应变随频率的增加而增加。

除外部荷载外，集料的密实程度和孔隙状态也对动变形影响显著。Lekarp等[8]研究认为密实状态的微小降低都会引起集料塑性应变的明显增加。

采用完全随机排列，不设重复。试验共设10个小区，小区面积108平方米，长30米，宽3.6米，每个小区3个厢面，每个厢面种2行共6行，行距宽行0.7米，窄行0.5米，株距0.27米，亩栽4100株，单株留苗。四周设保护行，中间不留过道。成熟时实收测产。

在针对碎石集料的研究方法方面，除了大型室内试验以外，基于离散单元法、非连续变形分析法等数值分析手段也获得了广泛的应用，特别是以颗粒流软件PFC3D(Particle Flow Code)为代表的数值分析软件，在碎石集料力学性质的分析中体现出了较高的效率与精度[9-10]。

本文基于颗粒流计算方法建立了碎石集料的动三轴剪切试验模型，通过对比室内试验结果对模型进行了标定与验证，进一步地分析了动荷载形式、集料密实状态对于整体动变形特性的影响。

预应力混凝土结构，是在结构承受荷载之前，预先对其施加压力，使其在外荷载作用时的受压区产生拉应力，用以抵消或减小外荷载产生的压应力，因此，桥梁梁部采用预应力混凝土结构，可以提高梁部的承载力，从而大大增强整个工程项目结构的稳定性。

1 颗粒流模型的建立与验证

1.1 颗粒流试样生成

颗粒流模型中，可用无摩擦的刚性圆柱墙体单元模拟三轴侧边界，用颗粒单元模拟碎石块体。上下界面同样采用刚性墙体模拟加载板。计算时上下界面同时加载，通过伺服程序监测颗粒与侧边界的接触力并缩放侧边界墙体的半径保证恒定围压。

颗粒流模型中的基本单元为球体单元，采用clump命令生成不规则块体模拟碎石颗粒[11]。本文分别针对三种粒径级配的碎石集料进行了分析，级配A、B、C分别为1.27～2.54、2.54～3.81、3.81～5.08 cm。直接按照粒径大小生成颗粒后再替代为clump块体，完成试样制备。

图1 碎石试样与接触力分布

生成的三轴试验的颗粒流模型见图1(a)，试样高0.6 m，底面直径0.3 m，不同颜色单元代表不同形状的块体单元。三轴剪切前试样内部的接触力分布见图1(b)，力链分布较均匀，试样内部块体颗粒接触良好。

通过对比试验结果并进行调整，选择线刚性接触模型描述颗粒间的接触，最终确定的模型参数见表1。

表1 三轴试验细观参数

单元参数取值颗粒密度/(kg·m-3)2 650摩擦系数0.5法向刚度/(MN·m-1)100切向刚度/(MN·m-1)10墙体摩擦系数0侧墙法向刚度/(MN·m-1)10上下墙法向刚度/(MN·m-1)100摩擦系数0

1.2 室内动三轴试验

为了验证颗粒流模型，进行了大型室内动三轴剪切试验。室内试验所采用的三轴试验系统见图2，试样高度0.6 m，试样底面直径0.3 m。试验过程中所用的集料及筛分设备见图3，铁路道砟经筛分后形成相应粒径级配的集料。

图2 室内试验装备

图3 试样材料与筛分设备

数值计算及室内试验均施加10 000次的循环荷载。在颗粒流计算过程中，不能直接对墙体施加荷载，本研究编制了伺服程序对上下墙体单元的目标应力进行控制，通过调整上下墙体的位置实现动力加载。根据相关研究，将列车荷载简化为半正弦荷载[12-13]，见图4。

本文利用SACS软件建立风机模型，边界条件利用弹簧模拟桩土作用。由于疲劳分析过程中载荷均属于中等载荷，此时桩基变形为小变形，故基础弹簧使用线性弹簧进行模拟，弹簧刚度系数根据土体参数与API规范进行确定。疲劳分析的总体模型如图1所示，模型采用一阶梁单元模拟整个结构的刚度，在模型不同高度位置设置集中质量，模拟实际风机结构质量的分布情况，建立整体模型的目的是获得整个风机基础结构在风浪作用下的结构动力响应。风机基础疲劳计算的局部有限元模型加密位置如图2所示，加密位置网格大小为t×t，t为板厚，主要用于风机基础安装电缆位置的结构开孔的疲劳强度计算。

图4 动荷载形式

进一步地，可以将图4所示的动荷载写为

(2)

式中：σreq为墙体的目标控制应力;T为动荷载的周期。

1.3 模型验证分析

室内试验与颗粒流模型计算得到的，循环荷载条件下，碎石集料的轴向累积塑性应变εap及回弹模量MR的变化情况，如图5、图6所示。

图5 轴向累计应变与加载次数的关系

图6 回弹模量随加载次数的变化

动荷载作用下，εap和MR的变化规律基本一致，在动荷载作用初期，增加速率较快，并随着加载次数的增加而持续增加，但速率不断减小，并最终趋于平稳。

与室内试验结果相比，数值计算得到的永久变形偏大，主要是受颗粒棱角度和模型侧边界的影响，数值模型的颗粒间缺乏足够的嵌挤与摩擦，因此产生了较大的塑性变形[14]。数值计算结果可以较好地反应不同粒径集料的永久变形情况，并展开进一步计算分析。

2 动荷载形式的影响

本部分研究了荷载条件对集料动变形特性的影响，颗粒级配为3.81～5.08 cm。计算方案见表2，分别选择四种围压水平、幅值和频率值，利用围压、幅值与加载频率的数值组合来描述某一计算工况。相应的动荷载形式见图7。

表2 动荷载影响计算方案

围压/kPa荷载幅值/kPa加载频率/Hz30、55.1、100、150165.41、5、10、1530、50、100、165.41

图7 不同动荷载形式

数值计算得到的荷载频率为1 Hz时，不同围压以及荷载幅值条件下，集料的累积轴向应变情况，见图8。集料的累积轴向变形随幅值的增加而持续增加，动荷载作用初期最为明显，随着荷载次数的增加，趋于稳定。

图8 荷载幅值对轴向变形的影响

数值计算得到的荷载幅值为165.4 kPa时，加载频率对于永久变形的影响，见图9。随着加载频率的增加，集料的永久变形也明显增加，其增加的幅度也在增加。

防晒化妆品的功效评价主要包括三个方面的指标：防晒化妆品防晒指数 （Sun Protection Factor,SPF值）测定、防水性能测试以及长波紫外线防护指数（Protection Factor of UVA,PFA值）的测定。其中SPF值是指引起被防晒化妆品防护的皮肤产生红斑所需要的MED与未被防护的皮肤产生红斑所需的MED之比。MED（Minimal Erythema Dose）即最小红斑量，是指引起皮肤清晰可见的红斑，起范围达到照射点大部分区域所需要的紫外线照射最低剂量（J/m2）或最短时间（秒）。

图9 加载频率对轴向变形的影响

碎石集料的回弹模量MR随荷载作用次数变化较大，随着集料持续密实，最终趋于稳定，此时，回弹模量只与集料所受到的应力水平有关。因此可针对加载次数为N=10 000次时的MR展开分析，多种动荷载形式下的级配C集料的MR与应力不变量θ的关系见图10，采用k-θ模型进行描述，模型拟合结果列于图中，随着应力不变量θ的增加，回弹模量持续非线性增加。

第一种方案适宜渠道为总干渠一闸至四闸之间渠段与位于河套灌区上、中游乌兰布和灌溉区域与解放闸灌溉区域干渠，上述渠道均为国管渠道、多有防渗措施、渠道水利用系数比较高，水量损失小，渠道断面尺寸较大，能够在停水期内储存滴灌用水，在渠道两岸农田可发展直接引黄滴灌。通过对河套灌区各级渠道调研及遥感分析，第二种方案基本是在乌兰布和灌溉区域内沈乌干渠与东风渠两侧，有部分闲置荒地，可满足建设调蓄水池条件。

图10 回弹模量与应力不变量的关系

3 孔隙状态的影响

针对细粒含量对砂土力学性质的影响研究表明，相比于整体密实状态，大粒径骨架颗粒之间的相互咬合、嵌挤与摩擦的影响更加显著[15-16]。为了评价孔隙状态对含细颗粒的间断级配碎石集料整体力学性质的影响，提出了骨架孔隙比es的概念，其定义为碎石-细砂石间断级配混合料中孔隙体积与大粒径骨架碎石颗粒的体积比。

集料的整体孔隙比e为孔隙的体积与固体颗粒体积之比为

(3)

式中：Vgs和Vrs 分别为小粒径的砂石颗粒以及大粒径的碎石颗粒的体积；Vv为孔隙的体积。

面向枪声定位应用的无线传感器网络中，多个传感器节点构成一个新的麦克风阵列，有效地扩大系统检测口径。枪声定位系统网络如图1所示。

在前面提到数据层主要对后续应用层功能表现提供数据支持，即数据层是系统的基础环节，该层通过对底层机载系统采集传输的数据进行解析和存储，并在存储之前做一定程度的预处理，使存储的数据更具有结构性。本文将数据层实现分为两个部分：Redis实时缓存与MySQL持久化，以分别应对实时监测与大数据量分析运用的功能实现。

小粒径砂石含量m为砂石的质量与混合料总质量的百分比为

(4)

式中：ρgs和ρgs分别为砂石颗粒和碎石颗粒的密度，本研究中分别取为1 950 kg/m3和2 650 kg/m3。

7.2 防治方法：①农业防治：清除杂草及枯枝落叶，减少虫源。②药剂防治：加强虫情检查，控制在点片发生阶段，用药喷雾防治：1.8%阿维菌素乳油1000倍液；73%克螨特乳油1200倍液。

进一步推导出式(5)所示的关系，其中w表示砂石与碎石颗粒的体积比为

(5)

则可以推导出式(6)和式(7)所示的结果为

(6)

(7)

则根据骨架孔隙比的定义可以推到出其计算公式为

(8)

为分析孔隙状态对动变形的影响，进行了颗粒流三轴剪切计算分析，动荷载见图4，试验围压55.1 kPa。

混合料中大颗粒粒径为3.81～5.08 cm，细颗粒粒径为1～10 mm，计算方案如表3所示。在相同的整体孔隙比条件下，加入5种含量的细颗粒配置混合料进行研究，并计算出相应的骨架孔隙比。

表3 孔隙状态影响计算方案

碎石级配/cm整体孔隙比细粒含量/%骨架孔隙比3.81～5.080.460.540.680.7400.4600.5400.68100.74200.81300.89401.01501.1800.74

计算得到的不同整体孔隙比级配C集料的动应力-应变曲线见图11，荷载作用次数为100。整体而言，随着孔隙比的增加，集料在第一次动荷载作用下的轴向塑性应变在不断增加，在100次动荷载作用后，塑性应变随孔隙比的增加趋势更加明显，对应于孔隙比为0.46、0.54、0.68、0.74时的塑性应变分别为0.47、068、1.09、1.71。随着轴向应变的增加，动应力应变的斜率在不断增加，反应的是回弹模量的不断增加，而随着孔隙比的增加，曲线的斜率也在增加，说明孔隙较大时，在荷载作用下孔隙被不断压缩的同时以塑性变形为主，回弹变形则相对较小。

抚顺干馏工艺产生的气态产物经冷凝回收系统净化分离后，其中的油、水蒸汽生成了页岩油和干馏污水，气体部分则生成干馏瓦斯。干馏瓦斯按用途可分成三部分：一部分送蓄热式加热炉加热至500～750 ℃，携带大量显热回到干馏炉，作为油页岩干馏的循环热载体；一部分作为蓄热式加热炉的燃料气用于加热循环热载体；剩余部分则供燃气锅炉、电站等作燃料气使用，最终得到蒸汽、电及热能。

图11 加载次数N=100集料的应力应变曲线

孔隙状态对于永久变形的影响见图12(a)，并进一步地将不同孔隙比试样在荷载作用次数分别为100、1 000和10 000次时的塑性应变见图12(b)，随着孔隙比的增加，塑性变形增加较明显。

教书育人是一个漫长而艰巨的工程。作为高中地理教师，笔者深深感到地理教学的任重道远，它需要我们不断研究学生的特点，发现教材的变化，探索出符合学生和课改要求的新教法，提高地理课堂效率，减轻学生的课业负担，符合日新月异的现代教育发展的需求。这个过程中，既有付出的艰辛，也有收获的喜悦，也深深体会到地理教学的乐趣，同时发现教学中总有挖掘不完的东西，时时刻刻都有超越自我的快乐和成就感。

图12 孔隙状态对永久变形的影响

骨架孔隙比对混合料轴向塑性应变的影响，分别统计荷载作用次数为100、1 000和10 000次时的塑性应变见图13。随着细粒的增加，骨架孔隙比不断增加，轴向塑性应变几乎随之线性增加，细粒的增加改变了颗粒间的接触方式，大粒径的骨架效应在不断减弱，从而在相同的荷载条件下产生更大的塑性变形。

教师在了解了课程标准要求、教材内容、学生基本情况之后，就需要选择正确的方法与策略，将设计的教学目标准确表达出来。只有清晰的表达才能保证在教学过程中能够得到有效落实与执行。因此，教师应该根据不同的教学内容与教学目标性质，选择合适的方法与策略。

图13 骨架孔隙比对永久变形的影响

如图14所示为各集料在荷载作用次数为10 000次时的回弹模量MR随集料整体孔隙比e以及掺入细粒后的骨架孔隙比es的变化。

图14 孔隙状态对回弹模量的影响

随着e和es的增加均呈现出快速衰减的趋势，拟合曲线列于图中。本质上讲，无论是孔隙比还是骨架孔隙比的增加，其实质都是大颗粒的含量的降低，因此可以认为大颗粒含量越低，集料的回弹模量越小。

4 结论

本文基于室内动三轴试验与离散元数值计算方法研究了动荷载作用下碎石集料的永久变形与回弹变形特性，分析了集料密实程度、动荷载幅值及频率对动变形的影响，并提出了骨架孔隙比的概念。主要结论如下：

(1)动荷载作用下，碎石集料的永久变形与回弹变形规律相似，在动荷载作用初期，增加较快，随着动荷载作用次数的增加，持续缓慢增加，并最终趋于相对平稳。

(2)荷载条件下对集料的变形特性影响较大，随着幅值的增加，集料的累积轴向变形持续增加，动荷载作用初期增加的最为明显，随后缓慢增加并最终趋于稳定；随着加载频率的增加，集料的永久变形也明显增加。

(3)循环荷载作用下碎石集料最终的回弹模量只与集料所受到的应力水平有关(应力不变量θ)，可采用k-θ模型描述回弹模量与应力水平的关系，随着θ的增加，集料的回弹模量非线性地增加。

(4)骨架孔隙比的概念可评价含细小颗粒的间断级配碎石集料的孔隙状态，其定义为孔隙体积与骨架大颗粒体积比，随骨架孔隙比的增加，集料的塑性变形越大，回弹模量越小。

参考文献：

[1]李海锋．有碴轨道结构的变形设计及其几何状态行为理论的研究[D]．上海：同济大学，2006.

[2]SEED H B，CHAN C K，LEE C E.Resilience Characteristics Subgrade Soils and Their Relation to Fatigue Failures[C]//International Conference on the Structural Design of Asphalt Pavements.Washington D.C.:TRIS，1962：611-636.

[3]SELIG E T，WATERS J M. Track Geotechnology and Substructure Management[M].London:Thomas Telford Publications，1994.

[4]STEWART H E. Permanent Strains from Cyclic Variable Amplitude Loadings[J].Journal of Geotechnical Engineering，ASCE，1986，112(6)：646-660.

[5]AKKE S J，SUIKER，ERNEST T S. Static and Cyclic Triaxial Testing of Ballast and Subballast [J].Journal of Geotechnical and Geoenvironmental，2005，131(6)：771-782.

[6]SHENTON M J.Deformation of Railway Ballast under Repeated Loading Conditions. Railroad Track Mechanics and Technology，1978：405-425.

[7]INDRARATNA B，THAKUR P K，VINOD J S.Experimental and Numerical Study of Railway Ballast Behavior under Cyclic Loading[J].International Journal of Geomechanics，ASCE，2010，10(4)：136-144.

[8]LEKARP F，ISACSSON U，DAWSO A.State of the Art.II：Permanent Strain Response of Unbound Aggregates[J].Journal of Transportation Engineering，ASCE，2000，126(1)：76-83.

[9]王鹏程，刘建坤，李旭，等.双轴压缩条件下碎石集料静力剪切特性研究[J].铁道学报，2014，36(6)：87-92.

WANG Pengcheng，LIU Jiankun，LI Xu，et al.Static Shear Behavior of Crushed Rock Aggregate Subjected to Biaxial Compression[J].Chinese Journal of Journal of the China Railway Society，2014，36(6)：87-92.

[10]王鹏程，刘建坤，李旭，等.直剪条件下碎石集料力学特性研究[J].铁道学报，2014，37(2)：103-108.

WANG Pengcheng，LIU Jiankun，LI Xu，et al.Mechanical Behavior of Crushed Stone Aggregate Subjected to Direct Shear Conditions[J].Chinese Journal of Journal of the China Railway Society，2014，37(2)：103-108.

[11]LIU J K，WANG P C，LIU J Y．Macro-and Micro-mechanical Characteristics of Crushed Rock Aggregate Subjected to Direct Shearing [J]．Transportation Geotechnics，2015，2：10-19.

[12]QIAN Y.Integrated Computational and Experimental Framework for the Assessment of Railroad Ballast Life-cycle Behavior[D].Urbana-Champaign: University of Illinois at Urbana-Champaign，2015.

[13]QIAN Y，SEUNG J L，EROL T，et al.Simulating Ballast Shear Strength from Large-scale Triaxial Tests：Discrete Element Method[J].Journal of the Transportation Research Board，2013,2374：126-135.

[14]李旭，王鹏程，敖国栋，等.多边形道砟颗粒堆积和压缩特性的DDA数值模拟研究[J].铁道学报，2014，36(4)：71-75.

LI Xu，WANG Pengcheng，AO Guodong，et al.Stacking and Compaction Process of Polygon Ballast Particles Simulated By DDA Method[J].Journal of the Transportation Research Board，2014，36(4)：71-75.

[15]BELKHATIR M，ARAB A，DELLA N，et al.Experimental Study of Undrained Shear Strength of Silty Sand: Effect of Fines and Gradation[J].Geotechnical and geological Engineering，2012，30(5)：1103-1118.

[16]RAO C E，TUTUMLUER，KIM I T. Quantification of Coarse Aggregate Angularity Based on Image Analysis[J].Journal of the Transportation Research Board，2002，1787：117-124.