DOI:10.3969/j.issn.1674 ̄0696.2019.10.08
重庆交通大学学报(自然科学版)
JOURNALOFCHONGQINGJIAOTONGUNIVERSITY(NATURALSCIENCE)
Vol.38 No.10
Oct.2019
沥青混凝土路面装配式基层结构力学特性研究
孙建诚ꎬ鲍克端
(河北工业大学土木与交通学院ꎬ天津300401)
摘要 沥青混凝土路面装配式基层结构相比传统的沥青路面结构具有环保 快速和经济等优点 为了深入研究沥青混凝土路面装配式基层结构的力学特性 运用 研究了预制板块的边长 厚度 弹性模量对基 有限元软件 层底部力学行为的影响 结果表明 当预制板块边长在 范围内 约为 时 装配式沥青混凝土路面基层预制板块边长对基底应力行为的影响程度大于预制板块厚度和弹性模量的影响 关 键 词 道路工程 装配式基层结构 有限元 力学分析 尺寸效应中图分类号
文献标志码
文章编号
结构基层底部最大竖向应变值最小 在一定范围内 预制板块的厚度和模量对基底最大竖向应变值的影响并不大
旜
斚
旜 旜 旜 旜 旜 旜 旜 旜 旜
旜 旜 旜 旜
旜
0 引 言
我国国民经济的迅速发展带来了交通建设的大发展ꎬ而公路建设作为交通建设的重要组成部分也随之迅速发展ꎮ十三五期间国家将大力推动节能低碳发展ꎬ鼓励标准化设计及工厂预制ꎬ综合利用废旧路面、疏浚土、钢轨、轮胎和沥青等材料以及无害化处理后的工业废料、建筑垃圾ꎮ笔者所研究的预制板块可以使用建筑工业产生的废料ꎬ对于保护环境有着积极的作用ꎮ
传统的公路施工工艺虽然已经相当成熟ꎬ但还是有很多人为因素影响着施工质量ꎬ例如施工质量把控不严ꎬ偷工减料等ꎮ而装配式道路可以减少人为因素对公路施工质量的影响ꎬ对建设合格的公路有着积极的意义ꎮK.TANG等[1]研究发现ꎬ在混凝土铺路块生产中使用建筑和拆除废料(C&DW)具有可能性ꎻC.RODRÍGUEZ等[2]研究表明ꎬ使用再生混合料(RMA)制备预制非结构混凝土具有可能性ꎻ王对路表弯沉进行了数值模拟计算ꎻ马健生等[4]针对火明等[3]通过建立有限元路面结构力学分析模型ꎬ
收稿日期:2018 ̄03 ̄19ꎻ修订日期:2018 ̄11 ̄05
第一作者:孙建诚(1969—)ꎬ男ꎬ河北冀州人ꎬ副教授ꎬ博士ꎬ主要从事道路方面的研究ꎮE ̄mail:sunjiancheng2000@126.comꎮ
第10期孙建诚ꎬ等:沥青混凝土路面装配式基层结构力学特性研究
51
新型的装配式道路基层结构研制了其填缝所用的自流平水泥基砂浆材料ꎬ并对其力学、耐久和疲劳性能进行了试验研究ꎻ严秋荣等
[5]
通过数值模拟手段ꎬ
对装配式水泥混凝土路面圆形企口缝的力学行为进行了研究ꎬ为装配式水泥混凝土路面应用提供了一定的理论参考ꎮ
现有的装配式道路研究大多集中于装配式路面面层的力学特性研究ꎬ而沥青混凝土路面装配式基层结构是将预制板块用于路面基层的修筑之中ꎮ笔者使用ABAQUS有限元软件和正交试验法研究预制板块边长、厚度及弹性模量对基层底部力学行为的影响ꎮ
1 沥青混凝土路面装配式基层结构
传统的装配式路面将预制混凝土板用于道路路面的铺装ꎬ不仅对预制板的力学性能有一定要求ꎬ而且对预制板的耐久性(如表面的抗裂性)及平整度有着较高的要求ꎬ故对所用预制板的材料要求较高ꎮ而装配式基层结构是将预制板块用于基层的铺装之中ꎬ能降低对预制板块平整度的要求ꎬ只需满足传统的道路基层所提供的力学性能即可ꎬ因此可以使用建筑施工中的废弃材料ꎮ将废弃的材料粉碎后进行分筛作为装配式基层的原材料ꎮ沥青混凝土路面装配式基层结构如图1ꎮ由于所研究的装配式基层板块与板块之间存在缝隙ꎬ故利用结构中的级配碎石层吸收应力来防止底部的裂缝扩散到面层导致路面损坏ꎮ
图1 装配式基层结构
Fig.1 Assemblytypebasestructure
2 沥青混凝土路面装配式基层结构力
学模型建立
2.1 有限元力学参数的确定
模型所选的材料结构层弹性模量、泊松比、厚度如表1ꎮ
表1 材料参数
Table1 Materialparameters
材料弹性模量E/MPa
泊松比μ厚度d/cm
沥青混凝土
1200
0.30
10
(续表1)
材料弹性模量E/MPa
泊松比μ厚度d/cm
水泥稳定碎石15000.2525级配碎石5000.3515预制板块31000
0.1520压实土
50
0.40
—
2.2 荷载加载方式
许多研究选择集中荷载作为道路的加载方式ꎬ这与道路实际的荷载作用方式不同ꎮ如果采用单轴双轮组的荷载模式更符合预制板块的实际工作状态ꎬ因此笔者决定采用单轴双轮组的荷载模式进行加载ꎮ车辆轴载采用标准轴载单轴双轮组BZZ ̄利以及考虑到车轮的实际作用情况100ꎬ胎压为0.7MPaꎮ为了简化计算ꎬ、将车轮与路面施加荷载的便
的接触形式简化为矩形ꎬ如图2[6]ꎮ
图2 车轮与路面接触简化(单位:cm)Fig.2 Simplificationofwheel ̄roadcontact
2.3 边界条件
沿着行车方向两端行车方向的两个边界定义好约束ꎬ约束类型为位移/转角ꎬ选择约束其法向和竖直方向ꎬ约束土基底部边界的3个方向(竖直、法向和切向)[7]2.4 土基的计算深度
ꎮ
在ABAQUS有限元软件中ꎬ构件模型的尺寸必须是有限的ꎬ模型中的土基计算深度越大则越符合实际情况ꎬ但是计算深度越大ꎬ软件的计算量也就越大ꎮ行车荷载在路基中引起的附加应力会随着土基深度的增加而不断降低ꎮ当达到一定深度后ꎬ由于荷载附加应力与路基自重应力之比很小ꎬ便可忽略车辆荷载的影响ꎮ这一深度即为由车辆荷载所引起的路基附加应力分布范围ꎬ并以此深度作为路基的应力计算深度[8]轴重作用下的路基应力计算深度ꎮ通过有限元模型分析不同轴型ꎬ得出路基应力计、算深度应不小于3.0m[8]中ꎬ将土基计算深度设定为ꎮ5.0在笔者所研究的模型mꎮ2.5 层间接触条件
为简化计算ꎬ在参考文献[9]的基础上ꎬ将模型的上部部分结构层(包括沥青混凝土、水泥稳定碎石
52 重庆交通大学学报(自然科学版)
第38卷
和级配碎石层)设为完全连续ꎮ同时将上部部分和装配式基层之间的接触连接应用罚函数公式进行接触模拟ꎬ其中接触的摩擦系数为0.6[9 ̄10]层和土基之间的接触连接同样采用罚函数进行接触ꎬ装配式基模拟ꎬ摩擦系数为0.5[9 ̄10]2.6 单元类型及网格划分ꎮ
在网格模块部分单元类型的选择中ꎬ通常隐式方法(Standard)能够求解线性和非线性问题ꎬ包括静态分析、动态分析ꎻ而显式方法(Explicit)适用于求解复杂非线性动力学问题和静态问题ꎬ特别是模拟短暂、瞬时的动态事件[9]使用的模型ꎬ使用隐式方法ꎮ参考文献(Standard)[9]求解模型并集合本次所ꎮ
模型采用六面体单元进行计算ꎬ可选用的单元基本性质有C3D8R(八节点六面体单元ꎬ减缩积分)和C3D20R(二十节点二次六面体单元ꎬ减缩积分)ꎮ由于剪力自锁ꎬC3D8R的计算精度低于C3D20Rꎬ但使用C3D20R的计算成本较高ꎮ为进一步说明模型网格划分的可靠性ꎬ以板块边长为1.5mꎬ厚度为
0.2表2ꎮ
m为例进行计算说明ꎮ网格收敛性计算结果如表2 网格收敛性计算分析结果
Table2 Gridconvergencecalculationanalysisresults网格大小/m
σ/MPaC3D8R
ε/mmσ/MPaC3D20R
ε/mm0.050.100.2810.3640.150.2190.6360.3590.6340.6360.20
0.0450.6440.038
0.6820.688
0.3570.355
0.6400.640
注:由表σ为基底最大拉应力2可以看出ꎬꎬε网格大小和单元类型对于应
为基底最大竖向应变
变计算的影响较小ꎬ其中最大的计算应变值仅比最小的计算应变值大8.5%ꎬ但是网格大小和单元类型对于应力计算结果影响较大ꎬ特别是对于C3D8R单元ꎬ其最大计算结果比最小计算结果大86.5%ꎬ而对于C3D20R单元类型ꎬ其最大应力值仅比最小应变值大2.53%ꎮ因此选择C3D20R单元类型应用于应力应变求解区域范围内ꎬ将C3D8R单元类型应用于不影响计算结果的计算区域ꎮ在参考相关文献的基础上[9]算精度较高的ꎬ对于模型上部部分和预制板块部分采用计C3D20Rꎬ对于土基部分采用精度较低的C3D8Rꎮ
依据黄仰贤二维有限元计算结果:单元边长比
的变化对挠度影响不大ꎬ但对应力计算结果的影响很大ꎮ因而ꎬ在单元划分时ꎬ应尽量保持单元的边长比为lꎬ且单元边长比不宜超过2[11]体单元ꎬ评价单元形状需要a/b和b/cꎮ两个指标对于三维实(a、b、c分别为单元的长度、宽度和高度)ꎮ取a/b=1ꎮ上部部分(包括沥青混凝土、水泥稳定碎石、级配碎石层)在竖直方向上采取3种不同的网格划分密度ꎬ其中b/c分别为1、0.80、1.33ꎻ在预制板部分和土基部分ꎬ取b/c=1ꎬ均能满足上述要求ꎮ
综上所述ꎬ指派上部部分和预制板部分网格属性为结构网格Standardꎬ分ꎻ指派土基部分网格尺寸为几何阶次为二次ꎬ网格大小为ꎬ属性为0.1mꎬ指派单元类型为0.4mꎬC3D20Rꎬ网格属性为结减缩积构网格ꎬ指派单元类型为Standardꎬ几何阶次为线性的ꎬ属性为C3D8Rꎮ计算模型如图3ꎮ
图3 计算模型Fig.3 Calculationmodel
3 沥青混凝土路面装配式基层结构力
学模型结果分析
3.1 预制板块尺寸对应力行为的影响
将静荷载分别作用于模型中心靠近板缝处和模6、型边缘靠近板缝处图7ꎮ
ꎬ如图4、图5ꎬ力学计算结果如图
图4 荷载作用于模型中间示意(单位:mm)Fig.4 Loadactingonthemiddleofthemodel
第10期孙建诚ꎬ等:沥青混凝土路面装配式基层结构力学特性研究
53
图5 荷载作用于模型边缘加载示意(单位:mm)Fig.5 Loadactingontheedgeofthemodel
图6 荷载作用于模型中间应力应变
Fig.6 Stressandstrainoftheloadactingonthemiddleofthemodel
图7 荷载作用于模型边缘应力应变
Fig.7 Stressandstrainoftheloadactingonthemodeledge
由图6可得ꎬ当荷载作用于模型中间时ꎬ预制板块的边长从1m增加到2mꎬ基底的最大拉应力增加了1.78倍ꎬ由此可得ꎬ预制板块的边长越大ꎬ基底的最大拉应力越大ꎮ基底的最大竖向应变随着预制板块边长先增加后减小ꎬ在预制板块边长为1.4m时出现了突变ꎬ其值在边长为1.4m时最小ꎮ当预制板块边长从1m增加到1.4m时ꎬ基底最大竖向
2应变减少了m时ꎬ基底最大竖向应变增加了9.0%ꎻ预制板块边长从15.3%ꎮ
1.4m增加到
从图7可以看出ꎬ当荷载作用于模型边缘时ꎬ基
底的最大应力随着预制板块边长的增长呈上升趋势ꎬ当预制板块的边长从1m增加到2m时ꎬ基底的1.4最大拉应力增加了83.9%ꎮ当预制板块的边长为
长为m1时mꎬ基底的最大应变处于最低值时和边长为1.4m时的基底最大竖向应变ꎮ预制板块边相差20.8%ꎻ预制板块边长为1.4m时和边长为2m
时的基底最大竖向应变相差17.0%ꎮ
综上所述ꎬ板块边长对基底的最大竖向应变值有一定的影响ꎮ在路面结构中ꎬ应变过大会导致道路的破坏ꎬ而在此次的模型分析中ꎬ基底的最大应力并不算大ꎬ无论荷载作用于模型中间还是模型边缘ꎬ最大的基底拉应力也仅为0.618MPaꎬ小于基底的承受能力ꎮ因此笔者将基底的最大应变选为主要控制1.4因素ꎮ从图6、图7可以看出ꎬ当预制板块边长为
的长期使用m时ꎬ基底的最大应变处于较低值ꎮ再依据实际施工的模数要求ꎬ有利于路面ꎬ预制板块的边长为1.5m比较符合实际施工要求ꎮ3.2 预制板块厚度对应力行为的影响
将荷载作用于模型中心处ꎬ改变预制件板块的1.5厚度mꎬꎬ得出基底最大应力和应变模量为31000MPa为例)ꎬ(以预制板块边长为
如图8ꎮ
图8 预制板块厚度对应力应变的影响
Fig.8 Influenceofprecastplatethicknessonstressandstrain
在荷载作用于模型中心处ꎬ随着预制板块厚度的增加ꎬ基底的最大拉应力值和最大竖向应变值均随之减小ꎮ由图8可得ꎬ基底最大拉应力和板块厚度-1大.892致86呈x+线0.718性关38ꎬ系Rꎮ线性拟合函数为y=2度从0.2m增加到了0.3m=0.952时ꎬ基底最大应力减小了14ꎮ当预制板块厚
54.32%ꎬ由此可得ꎬ预制板块的厚度越大ꎬ基底的最
54 重庆交通大学学报(自然科学版)
第38卷
大拉应力越小ꎬ预制板块的厚度对应力影响较大ꎮ基底最大竖向应变和板块厚度大致呈线性关系ꎬ线性拟0.980合69ꎮ函数预制板块的厚度越大为y=-0.48571x+ꎬ基底的最大竖向应
0.72976ꎬR2=变值越小ꎮ当预制板块的厚度从0.2m增加到0.3m
时ꎬ基底的最大应变仅减小了7.86%ꎮ由此可得ꎬ虽然增加预制板块的厚度可以减少基底的最大竖向应变ꎬ但是预制板块的厚度对基底的最大竖向应变值影响较小ꎮ
3.3 预制板块模量对应力行为的影响
将荷载作用于模型中心处ꎬ改变预制件板块(板块边长为1.5mꎬ厚度为0.2m)的模量ꎬ得出基底最大拉应力值和最大竖向应变值ꎬ如图9ꎮ图9 预制板块弹性模量对应力应变的影响
Fig.9 Influenceofprecastplateelasticmodulusonstressandstrain
从图9可以看出ꎬ基底最大应力随着预制板块模量的增大而增大ꎬ二者大致呈线性关系ꎬ线性拟合函数为y=0.00000115x+0.32265ꎬR2但是当预制板块模量从MPa25000MPa增加到=0.9723384ꎮ000制板块模量对于基底的最大应力影响不显著时ꎬ基底的最大应力仅仅增加了2.56%ꎬ可见预ꎮ基底最大竖向应变随着预制板块模量的增加而减小ꎬ并且呈线性关系0.667ꎬRꎬ线性拟合函数为y=-0.000001x+2到33000=MPa1ꎮ当预制板块模量从时ꎬ基底的最大应25变000仅仅MPa减少增加
1.24%ꎬ了应变影响不显著可以看出ꎮꎬ预制板块模量对基底的最大竖向综上所述ꎬ预制板块模量对基底的最大应力和最大竖向应变影响均不显著ꎮ
3.4 预制板块的正交分析
正交试验设计是研究多因素多水平的一种设计方法ꎬ能够选取有代表性的实验ꎬ使统计分析变得简单而有条理ꎮ考虑分析板块边长、厚度和模量等3种因素ꎬ每个因素选取4个水平进行比较ꎮ将荷载作用于模型中间得出计算结果ꎬ如表3ꎮ
表3 正交试验结果Table3 Orthogonaltestresults
试验板块边长m
/
板块厚度组数/mm
1/m
板块模量/基底最大基底最大拉应力/MPa应变21.00.2233MPa
31.00.24310000.16241.0290000.1450.64251.00.26270000.1270.63361.20.28310000.1050.62171.20.22330000.2420.61781.20.241.20.26270000.2240.616290000.1930.612109
1.50.28290000.1650.609111.50.22270000.3020.602121.50.24330000.2510.623131.50.26310000.2070.614142.00.28270000.1890.601152.00.22290000.4690.594162.00.242.00.260.28310000.4320.677330000000.4010.6680.377
0.6590.644
4、 表 根据表5ꎬ其中3K的模拟计算结果建立极差分析ꎬ如表i为对应该列因素第ꎬi个水平的试验数
据之和ꎮ例如ꎬ表4中第一列因素“板块边长”中K对应第一水平“1.0m”4个试验的基底最大拉应力
1之和ꎬR为各因素的极差ꎮ
由表4中的R值可知ꎬ3个因素对基底最大拉应力影响的主次因素为:板块边长>板块厚度>板块模量ꎮ此外ꎬ对于此次正交试验而言ꎬ通过比较各列因素Ki的大小ꎬ可以确定本次试验各因素的优水平ꎬ即在板块边长为33000MPa时ꎬ板底最大拉应力最小1.0mꎬ厚度为0.28ꎮ
m和模量为
表4 应力分析结果
Table4 Stressanalysisresults
MPa
分析指标
影响因素
板块边长板块厚度板块模量KK10.539K20.8241.175K30.9491.0521.0184R
1.6790.9281.0261.140
0.8360.9770.339
0.9700.056
向应变影响的主次因素为 由表5中的R值可知:ꎬ3板块边长个因素对基底最大竖>板块厚度>板
第10期孙建诚ꎬ等:沥青混凝土路面装配式基层结构力学特性研究
55
块模量ꎮ针对此次正交试验ꎬ通过比较各列因素Ki的大小ꎬ可以确定本次试验的优水平ꎬ即在板块边长为1.5mꎬ厚度为0.28m和模量为33000MPa时ꎬ板底最大竖向应变最小ꎮ
表5 应变分析结果Table5 Strainanalysisresults
mm
分析指标影响因素
板块边长板块厚度板块模量KK12.5132.558K22.4392.5272.5172.514K32.4324R
2.6482.4900.216
2.4572.5020.101
2.4990.018
4
结 论
1.41)当预制板块的边长在1~2m范围内ꎬ边长为
的最大应力值和最大竖向应变值最小m时ꎬ沥青混凝土路面装配式基层结构基层底部ꎬ综合考虑实际应用情况效减少沥青混凝土路面装配式基层结构基层底部的2)在一定范围内ꎬ可以取预制板块的边长为ꎬ增加预制板块的厚度可以有1.5mꎮ
最大应力值ꎬ但是对基底最大竖向应变值的影响并不大ꎮ考虑到几组基底的最大应力值均处于较低水平ꎬ因此可以将基底最大竖向位移作为主要控制因素ꎬ并不需要过度增加预制板块的厚度大拉应力及最大竖向应变影响并不太大3)在一定范围内ꎬ预制板块的模量对基底的最ꎮ
大于预制板块的厚度和弹性模量的影响4)预制板块边长对基底应力行为的影响程度要
ꎮꎮ
参考文献(References):
[1]TANGmadeMaterialswithKꎬrecycledSOUTSOSandTechnologies ̄InternationaldemolitionMNꎬMILLARDaggregatesS[C]G.Conference/Concrete/SustainablepavingonConstructionproducts
Sustainable
ConstructionMaterialsandTechnologies.Rotterdam:TaylorandFrancis ̄[2]Balkemaꎬ2007:77 ̄84.
constructionRODRÍGUEZnon ̄structuralandCꎬdemolitionPARRACꎬwastesCASADOasrecycledGꎬetmixedal.Incorporationaggregatesof
inProductionꎬ2016ꎬ127:152 ̄161.
concreteprecastpieces[J].JournalofCleaner[3]王火明ꎬ赵军.预制块体路面结构力学特性的有限元计算分析
[J].WANG公路交通技术analysisHuomingꎬꎬ2011(6):16 ̄20.
formechanicalZHAOcharacteristicsJun.Finiteofelementprecastcomputationblockpavementand
20.
structure[J].TechnologyofHighwayandTransportꎬ2011(6):16 ̄[4]马健生计及性能分析ꎬ孙大伟[J].ꎬ余地公路ꎬ等ꎬ2017(10):17 ̄21.
.装配式道路基层结构填缝材料配比设
MAperformanceJianshengꎬanalysisSUNDaweiꎬofcaulkingYUDiꎬmaterialsetal.Proportioningforfabricateddesignroadbaseand[5]structure[J].严秋荣Highwayꎬ2017(10):17 ̄21.
[J].YAN重庆交通大学学报ꎬ彭齐腾.装配式水泥混凝土路面圆企口应力应变响应
(自然科学版)ꎬcementQiurongꎬChongqingconcretePENGJiaotongwithUniversitycircularQiteng.tongueStress ̄strain2017ꎬ36(11):(NaturalandgrooveresponseSciencejoints[J].ofprefabricated
28 ̄32.)ꎬ2017ꎬ36(Journal11):of[6]28 ̄32.
余文成[J]ꎬYU中国水运ꎬ谭波.汽车荷载作用下水泥混凝土路面力学性能分析
concreteWenchengꎬꎬpavementTAN2008ꎬ8underBo.(9):Mechanical239 ̄243.
vehicleload[performanceJ].ChinaanalysisWaterTransportofcement2008ꎬ8(9):239 ̄243.
ꎬ
[7]李娣东南大学.预制混凝土板快速修复水泥路面结构设计研究ꎬ2015:18 ̄28.
[D].南京:
LIRepairDi.StudywithonPrecastConcreteConcretePavementSlabStructural[D].DesignNanjing:MethodSoutheastofQuick[8]林小平Universityꎬ2015:18 ̄28.
[J].ꎬLIN同济大学学报凌建明ꎬ官盛飞stressXiaopingꎬlevelofsubgradeLING(自然科学版ꎬ等.水泥混凝土路面路基应力水平分析
JianmingꎬundercementGUAN)ꎬ2010ꎬconcreteShengfeiꎬ38(4):545 ̄551.pavement[J].etal.AnalysisJournalonofTongjiUniversity(NaturalScience)ꎬ2010ꎬ38(4):545 ̄551.[9][M].严明星杭州ꎬ王金昌:浙江大学出版社.ABAQUS有限元软件在路面结构分析中的应用
YANElementMingxingꎬSoftwareWANGinPavementJinchang.ꎬ2016.StructuralApplicationAnalysisof[M].ABAQUSHangzhou:Finite[10]Zhejiang廖公云UniversityPressꎬ2016.
[M].LIAO南京ꎬ黄晓明.ABAQUS有限元软件在道路工程中的应用
ElementGongyunꎬ:东南大学出版社SoftwareHUANGinRoadXiaoming.ꎬ2016.
EngineeringApplication[M].Nanjing:ofABAQUSSoutheastFinite[11]姚祖康UniversityPressꎬ2016.社ꎬ1999.
.水泥混凝土路面设计[M].合肥:安徽科学技术出版
YAOAnhuiZukang.ScienceDesignandTechnologyofCementPublishingConcreteHouseꎬ1999.Pavement[M].Hefei:(责任编辑:罗素琴)
因篇幅问题不能全部显示,请点此查看更多更全内容
Copyright © 2019- vipyiyao.com 版权所有 湘ICP备2023022495号-8
违法及侵权请联系:TEL:199 18 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务