水泥道面传力杆裹附混凝土应力响应模拟分析

袁 捷， 鲁 行， 黄崇伟， 孙 陈， 殷 巍

(1. 民航飞行区设施耐久与运行安全重点实验室， 上海 201804；2. 同济大学 道路与交通工程教育部重点实验室， 上海 201804； 3. 上海理工大学 交通运输工程系， 上海 200093)

接缝是水泥道面的重要组成部分，通常分为胀缝、缩缝、施工缝三种，主要用以释放温度造成的道面内应力.我国机场道面设计规范规定在胀缝、部分缩缝和缩缝处的施工缝宜设置传力杆[1].传力杆的设置可以起到增加传荷能力、减小面板最大主应力、缓解路面板脱空和错台的发展速度等作用[2-5].

在传力杆设计理论体系中，采用Friberg提出的理论公式[6]，以Timoshenko和Lessels提出的理论模型为基础[7].Timoshenko将传力杆视为弹性介质中的半无限长悬臂梁，传力杆受到荷载作用后在弹性介质中产生波浪形的弯曲，混凝土对传力杆的支承应力与弹性介质的压应变或者传力杆的挠度成比例.在传力杆的设计方法中，主要考虑的问题是混凝土对传力杆的支承应力是否在混凝土的承受范围内，但是此设计方法对传力杆系统是否有效和耐久工作的考虑过于简单，因为由于传力杆和混凝土的材料性质的差异，以及接缝结构的特异性，在传力杆的裹附混凝土中极易出现应力集中，从而导致结构损伤，进而使传力杆系统的荷载传递能力逐渐失效[8]，所以不能仅仅笼统地考虑混凝土的支承应力.传力杆裹附混凝土的应力集中，混凝土的塑性变形，疲劳损伤都是亟待研究的问题.

目前国内外对于道面传力杆和混凝土相互作用的力学行为方面有一定的研究.由于试验观测混凝土的损伤规律比较困难，使用有限元模拟传力杆和混凝土相互作用成为一些研究者的主要研究手段.诸多学者如蒋应军和戴经梁[9-11]、Kuo[12]、Maitra[13]、刘旭峰[14]等虽使用3D有限元对传力杆和混凝土界面接触应力进行了分析，但对于界面接触应力的影响因素并没有进行深入研究，例如温度翘曲的影响，未考虑到道面温度和荷载叠加作用效果.

本文针对传力杆裹附混凝土细部应力集中这一问题进行深入研究，重点考察已有研究未涉及的温度这一影响因素，分析正负温度梯度和三种位置荷载的叠加作用效果，进而得到温度、三种位置荷载单独作用产生的应力集中分布规律，以及温度荷载叠加后的应力集中的变化规律.

1 有限元模型

1.1 模型及参数设置

有限元模型的尺寸根据实际机场道面尺寸确定，如图1所示，模型中材料参数为：水泥混凝土为热-弹性，弹性模量为36 000 MPa，泊松比为0.15，密度为2 400 kg·m-3，热膨胀系数为10-5°C-1；基层为纯弹性，弹性模量为1 400 MPa，泊松比为0.2，密度为2 300 kg·m-3；传力杆为纯弹性，弹性模量为2×105MPa，泊松比为0.3.模型在混凝土板和传力杆接触的部分使用更密的网格，如图1所示.在传力杆孔口周围划分规则细密的网格，如图2所示.

图1 有限元模型尺寸 (单位：cm)

图2 传力杆周围混凝土网格划分

实际情况中，传力杆一端会做防粘处理以保证其在混凝土中可自由滑动，故模型中传力杆依据实际情况一端使用tie连接，另一端设置摩擦系数0.8，且交错布置.共设置14根传力杆，其间距为0.35 m，按照图3中所示编号1～14，黑色区域为加载位置.其中奇数号的传力杆在受荷板中的一端为滑动端，偶数号的传力杆在受荷板中的一端为固定端，在非受荷板中则相反.飞机荷载采用B-737-800机型，主起落架构型为单轴双轮，轮印长宽为0.435 m×0.3 m，轮胎压力为1.47 MPa，横向轮距为0.86m.加载位置考虑如图3所示三种位置，分别为板边中心，板角传力杆滑动端，以及板角传力杆固定端.为使后文叙述简洁，三种荷载位置使用代码L-Mid(load-middle)，L-CUB(load-corner-unbonded)，L-CB(load-corner-bonded).

温度与荷载的叠加方法：为了使模型和实际情况更加吻合，分两步进行，道面板在温度作用下先产生变形，其后荷载作用在变形后的道面板上.

设板的上部比板的下部温度低时为“负温度梯度”，板的上部比板的下部温度高时为“正温度梯度”.根据实际温度测量结果，沿板竖向的温度梯度约为20°C·m-1，因此上下温度差如果为负温度梯度设为-8°C，如果为正温度梯度则设为8°C.

图3 传力杆布设和编号及加载位置

Fig.3Distribution and numbering of dowel and loading area

1.2 模型准确性验证

有限元模型的准确性通过将计算结果与Westergaard的解析解对比实现.解析解假设板为半无限板条和无限长板，板底充分支撑，温克勒地基，无板重.Westergaard解析解沿板横向中心线的板竖向位移值δt为[15]：

(1)

(2)

(3)

(4)

式中:E为板的弹性模量；v为泊松比；α为混凝土的热膨胀系数；t为板上下温度差(板上-板下)；b为板的宽度；y为y轴上的距离(与交通方向垂直)；k为地基反应模量；h为板厚.

计算负温度梯度为-8°C时的板竖向位移，与解析解进行对比.由于解析解未考虑基层的传荷，而使用综合地基反应模量，且不考虑板重，板长无限，故本研究除了计算实际情形(板长L=5 m，混凝土密度ρs=2 400 kg·m-3)以外，还计算了几乎没有板重和基层重的情形(为避免模型计算不收敛，混凝土密度ρs和基层密度ρb取1 kg·m-3)，板长很长即L=20 m的情形，以及无基层板的情形.计算得到位移如图4所示(图中NB代表无基层).

排尿性晕厥又称小便猝倒，是指在排尿开始、排尿过程中或排完尿离开厕所时突然晕倒。主要是由于血管舒张和收缩障碍造成低血压，引起大脑一时供血不足所致。晕厥持续的时间，少则数秒钟，多则半小时，一般休息几个小时就会恢复正常，不会留下后遗症。排尿性晕厥多见于中老年男性，一般好发在夜间，常常突然发生，之前多无先兆。患有肺结核、神经衰弱和气血两虚的病人易发生此种晕厥。此外，病后体虚、过度疲劳以及饮酒等也可诱发这种现象。

图4 板的竖向位移计算值与理论值比较

Fig.4 Comparison of computed results and theoretical results on vertical displacement of slab

从图4中可知，当无基层，板长很长且密度为1 kg·m-3时，板的竖向位移曲线和理论解几乎完全重合，充分说明了有限元计算的准确性.有基层时(L=20 m，ρs=ρb=1 kg·m-3)板的挠曲线在板中位置略有不同.从图4中挠曲线①(L=20 m，ρs=ρb=1 kg·m-3)过渡到④(实际情形)，先后反映了板长、板密度、基层密度的影响.

2 单独作用下的传力杆裹附混凝土应力分析

2.1 温度梯度的单独作用效果

在正负温度梯度作用下，板会产生温度翘曲，以负温度梯度为例，板的形态见图5，板边向上翘曲.从图6可看出，设置传力杆的板在传力杆的约束作用下挠度降低约15%，传力杆对翘曲的约束作用是温度梯度下传力杆周围混凝土的应力集中的主要原因.

图5 混凝土板在负温度梯度-8°C影响下的翘曲(放大2 000倍后的效果)

图6 传力杆对温度翘曲的约束作用(以负温度为例)

由于传力杆的模量和强度均比混凝土大很多，故在传力杆和裹附混凝土的相互作用过程中，混凝土必定先出现破坏或损伤，且最易在孔口处出现应力集中，见图7b.为了进一步分析传力杆裹附混凝土的应力状态，在传力杆圆周上在圆周上划分角度作为混凝土位置的标定(见图7c)，角度从孔口底部逆时针方向增大，根据右手定则判定的法向方向指向混凝土外部.按照此角度规定，提取温度梯度作用下最大主应力在传力杆圆周的应力.正负温度梯度引起的最大主应力的最大值出现在6号传力杆固定端的圆周上，正温度梯度下出现在13号杆固定端，分布如图8所示.根据最大主应力屈服准则，最大主应力超过单轴拉伸屈服强度时，产生屈服.负温度梯度产生的最大主应力未出现压应力且绝对值较大，正温度梯度产生的最大主应力有压应力，并且绝对值较小.从而可知负温度梯度更容易引起混凝土的破坏.

图7 分析传力杆周围应力集中的角度规定

图8 温度梯度作用下传力杆裹附混凝土最大主应力

2.2 荷载的单独作用效果

在荷载单独作用下，传力杆裹附混凝土的应力如图9所示.从图中可见，在三种荷载单独作用下混凝土会出现明显的应力集中.可以看到以下规律：作用在板角的荷载引起的混凝土应力集中大于板边荷载引起的应力集中，作用在传力杆滑动端上的荷载引起的混凝土应力集中大于作用在传力杆固定端上的荷载引起的混凝土应力集中，且混凝土的应力集中突出反映于传力杆滑动端裹附的混凝土(表1).

图9 荷载作用下传力杆裹附混凝土最大主应力

荷载作用位置板边中心板角滑动端板角固定端最大主应力位置9号传力杆滑动端1号传力杆滑动端14号传力杆滑动端

3 温度、荷载叠加作用下的传力杆裹附混凝土应力分析

3.1 温度梯度和板边中心荷载叠加

从图10中可以看到温度梯度和板边中心荷载叠加作用对各应力的影响.可以看到负温度梯度使最大主应力提高约26.9%，正温度梯度使最大主应力的降低约37.6%.

图10 传力杆裹附混凝土的最大主应力(L-Mid和ΔT叠加)

Fig.10 The maximum principal stress of concrete around dowel (under superposition of L-Mid and ΔT)

3.2 温度和板角荷载叠加(荷载于传力杆滑动端)

从图11中可以看到温度梯度和板角滑动端荷载叠加作用对最大主应力的影响.负温度梯度使最大主应力升高约15.2%，正温度梯度使最大主应力降低约15.3%.

Fig.11 The maximum principal stress of concrete around dowel (under superposition of L-CUB and ΔT)

3.3 温度和板角荷载叠加(荷载于传力杆固定端)

从图12中可以看到温度梯度和板角固定端荷载叠加作用时最大主应力的变化.负温度梯度对混凝土的应力集中不利，最大主应力升高18.5%，正温度梯度使最大主应力降低约21.0%.

图12 传力杆裹附混凝土的最大主应力(L-CB和ΔT叠加)

Fig.12 The maximum principal stress of concrete around dowel (under superposition of L-CB and ΔT)

4 结论

本研究使用数值模拟分析了温度和荷载单独作用和叠加作用条件下的机场道面接缝处传力杆裹附混凝土的应力响应，主要结论如下：

(2) 在温度单独作用下，传力杆对温度翘曲有约15%的挠度约束作用.

(3) 对于三种荷载位置，作用在板角传力杆滑动端上方时，混凝土内产生应力最大.混凝土的应力集中突出反映于传力杆滑动端裹附的混凝土.

(4) 温度和荷载叠加时，不论荷载作用在哪种位置，负温度梯度均不利于传力杆裹覆混凝土的应力集中，正温度梯度对其有利，因此在温度骤降时行车更容易造成传力杆系统的损坏.