去年夏天我去南方某高速服务区做路况调查，红外测温枪打向路面，表显68℃。沥青表面被晒得泛出油光，重车道上的车辙深度目测有三四厘米。旁边一位做养护的兄弟叹气：通车才五年，车辙就成这样，离设计寿命还早着呢。
这句话我一直记着。在沥青路面结构里，车辙从来不是荷载单方面压出来的，温度场才是那个经常被低估的隐形推手。我做路面结构分析这些年，越来越确信一件事：凡是车辙预估偏差大的项目，往回溯，多半是温度场输入出了问题。
这篇是"沥青路面结构车辙温度场分析案例系列"的总纲篇。我会把整套分析逻辑串一遍：为什么必须从温度场入手、热参数怎么取、边界条件怎么设、温度场怎么算、再如何把温度场结果转换到车辙量值，最后结合几个典型工程场景复盘。适合设计院做路面结构验算的工程师、养护咨询机构的技术人员，也适合正在做相关课题的道路研究生当一份实战索引。
1. 车辙不是被压坏的，而是被"热软"之后再压坏的
1.1 沥青材料是条件屈服材料：温度决定它的脾气
沥青是典型的粘弹性材料，在低温下像玻璃一样脆，常温下表现出弹性，高温下则趋向粘性流动。路面上的车辙，说白了就是沥青混合料在高温状态下抗剪切能力急剧下降，然后在车轮反复碾压下发生了不可恢复的蠕变变形。
沥青混合料的力学性能对温度的敏感程度，远超出一般人的直觉。我举一个典型数据：相同的密级配沥青混凝土，在15℃时动态模量可能达到8000MPa以上，到了60℃时往往掉到只有200-400MPa，差了一个数量级以上。作为对比，水泥混凝土在相同温度范围内，弹性模量的变化幅度通常在10%以内。这就是为什么沥青路面会在夏天出问题，而水泥路面很少因为单纯高温产生结构性病害。
既然材料的"真实刚度"完全由温度决定，那么脱离温度场去做车辙分析，本质上就是在用一个不确定的模量去算一个对模量极其敏感的输出量，结果自然不可靠。
1.2 车辙形成三要素：高温、重载、低抗剪，缺一不可
车辙不是单一因素造成的，但三个必要条件里，温度往往是那个"开关"。
- 高温使沥青胶结料粘度降低，混合料内部的集料骨架失去了足够的约束力
- 重载（尤其是超载车辆）在轮迹带产生远大于设计荷载的竖向应力和剪应力
- 材料抗剪能力不足时，混合料发生侧向流动，轮迹带下陷、两侧隆起
在北方地区，同样一段高速公路，设计荷载相同、结构层相同，车辙病害的严重程度通常远低于南方高温地区。不是北方的车更轻，而是路面温度达不到让沥青软化的阈值。沥青路面表面温度达到50-60℃以上时，中面层的温度往往还在40-55℃之间，这个区间恰好是沥青混合料车辙变形高速发展的温度窗口。
这也是为什么温度场分析必须做在车辙预估前面——你要先知道路面内部哪些深度、在一天中的哪些时段处于"危险温度区间"，才能知道最不利位置和最不利时刻在哪里。
1.3 温度场分析要回答的四个工程问题
把温度场分析放进车辙评估的框架里，它要回答的问题其实非常明确：
- 路表温度在极端高温日的峰值是多少，比气温高多少
- 温度沿路面深度方向如何衰减，面层中下部是否仍然处于高温状态
- 温度峰值出现的时间随深度如何滞后，重载交通高峰是否与高温时段重叠
- 温度场的日累积效应对车辙的贡献有多大
这四个问题直接决定了你在车辙预估时怎么分层、怎么取代表温度、怎么确定荷载作用时间。我在案例中反复验证过：按规范查表取一个"年等效温度"和按实际气象数据算出来的场温度分布，最终车辙预估结果可能相差30%以上。这就是精细温度场分析的价值所在。
2. 把温度场算准的前提：热参数与边界条件
2.1 材料热参数：导热系数、比热容和密度缺一不可
温度场计算在数学上就是求解热传导方程，对沥青路面这种层状体系，最常用的是沿深度方向的一维热传导模型：
[ \rho c \frac{\partial T}{\partial t} = \lambda \frac{\partial^2 T}{\partial z^2} ]
其中ρ是材料密度，c是比热容，λ是导热系数。这三个参数决定了热量在路面结构里传递的快慢。看起来简单，但取值上有很多讲究。
常见沥青路面材料的推荐取值范围我整理在下面：
| 材料类型 | 密度kg/m³ | 比热容J/(kg·K) | 导热系数W/(m·K) |
|---|---|---|---|
| 密级配沥青混凝土 | 2300-2500 | 850-1000 | 1.2-1.8 |
| SMA沥青混合料 | 2400-2600 | 850-950 | 1.4-1.8 |
| 水泥稳定碎石基层 | 2200-2300 | 800-950 | 1.3-1.7 |
| 级配碎石 | 2000-2200 | 850-1000 | 1.5-2.0 |
| 土基 | 1800-2000 | 1000-1200 | 1.0-1.5 |
实际分析时我喜欢用"热扩散系数"a = λ/(ρc)来快速判断材料传热特性。对普通沥青混凝土，a大约在6×10⁻⁷到8×10⁻⁷ m²/s。热扩散系数越大，路面内部温度越容易受到大气温度波动的影响。
这里要提醒一句：SMA这类骨架密实型混合料，因为空隙率小、粗集料含量高，导热系数往往偏大，计算时不要直接套用普通AC混合料的参数，否则温度场会有偏差。
2.2 路表热交换边界：太阳辐射、对流和有效辐射三管齐下
温度场计算里最容易产生误差的地方，不是材料内部的热传导，而是路面顶部的边界条件。路面和大气之间的热量交换，主要由三个途径组成：
- 太阳短波辐射：沥青路面吸收的太阳辐射能量，入射强度在一天中呈正弦曲线变化，正午前后达到峰值
- 对流换热：路面与近地表空气之间的热交换，受风速影响很大，风速越大换热越强
- 路面有效辐射：路面以长波辐射形式向外散失热量，夜间尤其明显
太阳辐射吸收率是其中一个关键参数。新铺的沥青路面呈黑色，对太阳辐射的吸收率可以达到0.90-0.95；经过几年服役后表面泛白、集料裸露，吸收率会下降到0.85左右。很多分析直接用0.90一个值算到底，对长期性能评估来说会低估路表温度。
气温的日变化通常用正弦函数近似描述：
[ T_a(t) = T_{mean} + A \sin(\omega t + \varphi) ]
其中振幅A取当日最高气温与最低气温差值的一半，ω为24小时周期的角频率。这个简化对工程分析足够用，但如果手头有逐小时气象数据，直接插值更准确。
2.3 初始温度场：稳态解还是预热期？这里有个关键选择
瞬态温度场分析必须给定初始条件。一个常见的坑是：直接假设路面内部初始温度处处相等，比如全结构20℃，然后开始计算。结果头几个小时的温度场完全处于"虚假状态"，如果只算24小时，输出结果的前半段基本不能用。
我常用的处理方式是"预热法"：先用同一个月或者同一季节的平均气象条件，连续计算7-10个周期（日），让路面结构的温度场达到准稳态循环，然后把这个循环末刻的温度场作为正式分析的初始条件。在通用有限元软件里这相当于多花一点计算时间，但换来的是整个分析时段内温度场都处于真实状态。
如果你用的自编程序或者电子表格做差分计算，预热期还可以手动实现：先按恒定边界条件求稳态解，再叠加日变化扰动。总之不要把初始温度场当成一个可以随便填的数字。
3. 一个完整算例：极端高温日路面温度场的逐时演化
3.1 算例结构与气象输入
下面这个算例来自我做过的一段南方某省高速公路的评估项目，结构组合非常典型：
- 4cm SMA-13上面层
- 6cm AC-20中面层
- 8cm AC-25下面层
- 36cm水泥稳定碎石基层
- 土基
选取的是当年最热月份里日最高气温39℃、日最低气温27℃的一天。实测太阳辐射峰值约900W/m²，日平均风速1.5m/s。路表太阳辐射吸收率取0.90，沥青混凝土热扩散系数取7×10⁻⁷ m²/s。
这个结构组合的代表性在于：18cm的沥青面层厚度，恰好覆盖了日温度波动的主要影响深度。面层以下的水泥稳定碎石基层和土基，温度日变化已经很微弱，计算时作为半无限体底部恒温边界处理即可。
3.2 一维热传导的差分解法
用手算表或者Python写一个简单的显式有限差分程序，几分钟就能跑完。将路面深度方向按1cm划分网格，时间步长取30秒，差分格式为：
[ T_i^{n+1} = T_i^n + \frac{a \Delta t}{\Delta z^2} (T_{i+1}^n - 2T_i^n + T_{i-1}^n) ]
显式格式有一个稳定性条件，aΔt/Δz²必须不大于0.5。带入a=7×10⁻⁷ m²/s、Δz=0.01m，算得Δt最大不能超过357秒，取30秒绰绰有余。
实际写代码时，路表边界用热流平衡条件，把太阳辐射、对流、有效辐射折算成等效热流施加在第一个节点上。这个边界处理是程序里最容易出bug的地方，建议单独写成函数并做能量守恒校验。
3.3 计算结果的三个关键规律
图我就不贴了，直接说三个对车辙分析最有价值的规律：
第一，路表温度显著高于气温。在日最高气温39℃的条件下，计算得到路表最高温度为64℃左右。沥青路面吸热能力强，路表温度比气温高20℃以上是常态。这也是为什么很多设计人员在气象站数据里找不到"路面温度"，因为气象站测的从来不是路面。
第二，温度沿深度迅速衰减且峰值时间滞后。用半无限体周期温度场的解析解可以验证：相位滞后与深度成正比，振幅衰减按指数规律。对于本文算例的热扩散系数，大约在10cm深度处，温度峰值比路表滞后2小时以上，振幅也衰减到路表的50%左右。也就是说，上面层中午12点到下午2点处于最高温状态，而中面层的危险时段可能出现在下午3点到4点。
第三，面层中下部的温度远比你想象的高。下午2点时，4cm深度处的温度大约55-60℃，8cm深度约50-55℃，到16cm深度仍有45℃左右。整个沥青面层（18cm）在午后都处于40℃以上，这意味着在夏季高温日的下午时段，整个面层结构的抗变形能力都处在全年最低水平。
| 时刻 | 路表 | 4cm | 8cm | 12cm | 18cm |
|---|---|---|---|---|---|
| 8:00 | 31 | 32 | 33 | 33 | 33 |
| 12:00 | 54 | 47 | 43 | 40 | 37 |
| 14:00 | 64 | 55 | 49 | 45 | 41 |
| 16:00 | 52 | 51 | 48 | 44 | 41 |
| 20:00 | 35 | 38 | 39 | 39 | 38 |
这个温度分布形态有一个重要工程含义：车辙不只发生在上面层，中面层下部的黏弹性流动同样贡献可观。有些设计人员只在上面层取一个高温值去算车辙，实际上忽略了中面层的贡献，预估结果自然会偏小。
4. 从温度场到车辙量值：力学分析的传递路径
4.1 温度-模量关系是衔接温度和力学的桥
温度场本身不是最终输出，它要通过材料本构模型转换为力学响应。对沥青混合料而言，最直接的桥梁是动态模量。
沥青混合料的动态模量|E*|与温度和加载频率密切相关。相同加载频率下，温度从15℃升到60℃，动态模量可以下降一个数量级；相同温度下，加载频率降低（相当于车速减慢），模量同样明显下降。这就是时温等效原理在路面力学中的直观体现——低速重载在高温环境下对路面造成的损伤，远不是"荷载增大"四个字能概括的。
在车辙预估中，我习惯的做法是：把沥青面层按2cm一层细分，根据温度场结果给每一层赋予对应温度下的模量，然后用弹性层状理论或者有限元算出各层的竖向应力应变。这样处理的精度远高于把整个沥青层取一个平均模量。
4.2 车辙预估的两种路线：经验模型与机理模型
目前工程实践中用得比较多的车辙预估方法，大体分成两类：
一类是经验-力学法。以美国MEPDG（力学-经验路面设计法）为代表，先计算沥青层内各分层的回弹应变，再用永久变形与回弹应变的经验关系累加得到车辙深度。这类方法的优点是与实测数据标定充分，缺点是对温度输入的敏感性很高——温度取低5℃，车辙预估结果可能差20%以上。
另一类是粘弹性有限元法。把沥青混合料看成粘弹性材料，用蠕变参数或广义Maxwell模型直接模拟荷载反复作用下的累积永久变形。这种方法理论上更严密，能模拟出蠕变第三阶段的加速流动，但参数标定工作量大，计算耗时高，工程批量评估时不太友好。
对于常规的养护评估项目，我更推荐第一类方法。一方面它对输入参数的要求相对简单，另一方面行业积累的标定数据多，结果容易和实测车辙对得上。
4.3 案例验证：一个把温度场用起来的完整闭环
回到前面的算例。那段高速运营5年后，重车道实测平均车辙深度16mm，最大20mm。我在评估中用了两种温度输入方式：
第一种，按规范简化方法取年等效温度。计算得到车辙深度约10mm，明显偏小。
第二种，用极端高温月实测气象数据建立温度场，再把温度场逐时结果转换为各层代表温度。计算得到车辙深度约15mm，与实测吻合很好。
差异的核心在于：规范等效温度是全年平均意义上的代表值，但车辙发展主要在夏季高温月份的午后时段。沥青面层在这个时段内的真实模量远低于全年平均值，用"平均温度"去算等效于人为提高的模量，车辙自然被低估。这不是规范方法不行，而是它的适用场景更多在于新建路面设计验算。对于评估已开裂、已出现车辙的在役路面，我强烈建议用精细化温度场而不是直接套等效温度。
5. 案例系列的关键对比：三种特殊场景的温度响应差异
5.1 长大纵坡路段：温度场相同，但荷载作用模式变了
同在一条高速上，长大纵坡段的车辙普遍比平坡段严重，这是共识。但从温度场角度去看，纵坡段的问题不是温度更高，而是荷载作用时间更长、受力模式更不利。
上坡方向的大型货车在纵坡上减速行驶，低速意味着作用在路面上的荷载频率降低。沥青混合料的模量随荷载频率降低而下降，相当于把路面"软化"了。在连续高温天气下，坡道路面长时间承受低频率重载，中面层的剪应力水平显著升高，车辙就从这里集中爆发。
我在分析纵坡段时，会特别注意把荷载作用时间从标准工况的0.02秒量级延长到对应实际车速的水平。温度场本身不变，但力学分析里每一层的"有效模量"完全不同，车辙预估结果差距非常大。
5.2 交叉口与公交站区：驻车荷载和反复起步加速
交叉口和公交站区的车辙，贡献最大的不是行驶荷载，而是驻车和低速起步阶段的长时间荷载。
车辆在交叉口停车等待时，路面在高温条件下长时间承受静载，沥青混合料发生显著的蠕变变形。绿灯亮起后车辆加速驶离，水平力急剧增大，又给路面叠加了一次强剪切的循环。反复多次之后，轮迹带位置就形成典型的辙槽。
这类场景的温度场分析和普通路段没有区别，但在力学分析时要把"荷载作用时间"从行驶状态改为驻车状态的代表值。我见过不少分析报告在交叉口车辙评估里仍然按行驶车速取荷载频率，结果自然偏小。这不是温度场的问题，而是温度-荷载联动关系没有考虑到位。
5.3 桥面铺装：底部边界条件与路基路面完全不同
桥面铺装的温度场和普通路基路面有明显差异，很多人忽略了这一点。
普通路面结构的底部是土基，温度在数米深度基本恒定，可以作为恒温边界处理。但桥面铺装在钢箱梁或混凝土箱梁上，梁体本身的温度随外界气温变化，桥面板底部不是一个恒温边界。夏季白天，钢箱梁内部温度可以达到50℃以上，通过桥面板向上传递热量，导致铺装层的温度比同等条件下路基路面更高。
我在做某跨江大桥的铺装层温度场分析时，曾对比两种底部边界：一种按土基恒温边界，一种按实测箱梁温度变化曲线施加。结果显示，在夏季高温日，采用实测梁体温变边界时，铺装中部的温度比恒温边界高出5-8℃。这个差距对改性沥青混合料的抗车辙评估是决定性的。
6. 这个系列里最容易翻车的四个细节
6.1 太阳辐射吸收率：新路面和旧路面差别很大
新建沥青路面表面呈深黑色，辐射吸收率接近0.92。服役三五年后，表面沥青膜被磨掉，粗集料裸露，颜色变浅，吸收率可能降到0.85甚至更低。如果评估的对象是运营多年的老路，还按0.92去算，路表温度会偏高2-3℃，车辙预估结果会偏保守，但偏保守和准确是两回事。我建议有条件的情况下用实测表面温度反算吸收率。
6.2 初始温度场不要拍脑袋，一定要预热
前面已经说过预热的问题，这里再强调一次。我见过有人用稳态温度场作为初始条件直接做瞬态分析，结果在第一个24小时内，路面深部的温度一直在"漂移"，直到两三天后才稳定下来。如果分析只取一天的气象数据且没有预热，输出结果会受到初始条件的严重污染。
6.3 显式差分要注意时间步长和空间步长的匹配
用显式有限差分计算温度场，时间步长和空间步长必须满足稳定性条件，否则会出现数值振荡甚至直接发散。计算不收敛还好发现，最怕的是刚好处于临界状态，温度场看起来正常，实际上某些深度出现了非物理的波动。我建议程序里自动检查aΔt/Δz²的值，大于0.4就报错，留出安全余量。
6.4 用年平均气温代替高温季温度，误差会被放大
规范简化方法里的年等效温度是有意义的，但它考虑的是全年累计损伤的等效。如果你做的是"某个高温日"或者"夏季连续高温周"的极端工况分析，就不要用年平均气温去代替。路面车辙的关键发展期就是那么几十个小时，用年平均温度会把最危险的时段抹平，算出来的车辙深度自然不痛不痒。
温度场分析这件事，门槛不在公式，而在对边界条件、材料参数和工况选择的理解深度。同一个软件、同一套网格，不同人做出来的结果可能差出好几倍。上面这些细节，都是我在一个个案例里真金白银换来的教训。后续这个系列里，我还会针对纵坡段、交叉口、桥面铺装分别展开完整的分析过程和参数表格，把每一步怎么取值、怎么建模、怎么和实测对标都摊开来讲。
