连续交通流模型及数值模拟.doc

连续交通流模型及数值模拟
[摘要]本文对现有的交通流宏观模型进行了研究，总结了各种模型的思想、优缺点以及适用条件，在此基础上，选取了Payne模型离散格式进行数值模拟，选取了某段高速公路的交通流作为模拟对象，展现了Payne模型模拟交通流的可行性。
[关键字]连续交通流；离散格式；数值模拟
0引言
交通流理论研究加深了人们对复杂多体系统远离平衡态时演变规律的认识，促进了统计物理、非线性动力学、应用数学、流体力学、交通工程学等学科的交叉和发展等多学科的交叉渗透和相互发展。交通流理论研究的对象是离散态物质，是一个复杂的非线性体系，对这类物质运动规律的描述，尚无成熟的理论。
在宏观的连续流模型中，交通流被比拟为连续的流体介质，即将流量、速度和密度等集聚变量视为时间和空间的连续函数。模型包含时间和空间的状态方程，考虑了车辆的加速度、惯性和可压缩性，能够合理准确描述交通流的动态特性，相比微观模型有更大的优势。连续流交通流模型通常用密度（k）、速度（u）、流量（q）三个变量来描述[1]。
1连续交通流模型
1.1LWR模型
1955年，Lighthill&Whitham提出了第一个交通流的流体力学模型——流体运动学模型[2]，随后P.I.Richards独立地提出了类似的交通流理论。LWR模型用k(x,t)和u(x,t)表示t时刻位于x处的交通流密度和平均速度，他们满足流体力学的连续方程：
（1-1）
此方程反映了车辆数守恒，其中g(x,t)是流量产生率，对没有进出匝道的公路，g(x,t)=0，
对进口匝道，g(x,t)>0，对出口匝道，g(x,t)=0。k为交通密度，也称为交通流量；x，t
分别为空间测度和时间测度。设u为空间平均速度，则存在以下关系：
(1-2)
对于平均速度u(x,t)，假设平衡速度——密度关系：
(1-3)
以上3个方程构成了完整的一阶连续交通流模型，LWR模型的优点是简单明了，可以采用流体力学和应用数学中的成熟工具进行分析，而且可以描述诸如交通阻塞形成和消散之类的交通现象，但是，由于该模型的速度是由平衡速度密度关系决定，并且没有考虑加速度和惯性影响，因此不适用于描述本质上处于非平衡态的交通现象，例如车辆上、下匝道的交通、“幽灵式”交通阻塞、交通迟滞、时走时停的交通等。于是，后来的学者们引进了高阶连续介质模型，考虑了加速度和惯性影响，将动量方程代替方程(1-3)。
1.2Payne模型
Pipes于1953年提出交通流加速度的一般表达式：
（1-4）
1971年，Payne根据LWR模型的思想，假设交通流速度是动态变化的，在引用连续性方程时，引进运动方程，导出高阶连续模型[3]。Payne从车辆跟驰理论的概念提出平均速度u与密度k存在以下关系：
（1-5）
并取，——车辆跟驰理论在的延滞时间，对上式分别作关于和的Taylor展开，得到：
（1-6）
上式中，：预期指数，将平衡速度简写成。
从而建立了由如下三个方程构成的Payne模型：
（1-6）
（1-7）
（1-2）
式(1—6)的右边第一项为期望项，反映驾驶员对前方交通状态改变的反应过程；第二项式弛豫项，描述车流速度在时间内向平衡速度的调整，最优速度函数和其他参数通过道路实测和参数辨识确定。
1979年，Payne编制了著名的FREFLO软件，有史以来第一次将交通流仿真模型应用于工程实践。但Payne模型并未充分考虑整个弛豫过程，而只是将其定为一个常数的弛豫时间，即使处于平衡状态时，弛豫时间变为零，在实际应用中，出现了一些问题。Rathi等人指出，使用Payne模型，车流速度到平衡态速度的调节过程过于缓慢。Ross也发现，当道路拓扑特性和交通量在短时间内突变时，由于车流速度调整到平衡速度过程缓慢，无法捕捉到实际交通流动态特性。Castillo等对Payne模型进行了线性稳定性分析，发现车辆总是在稳定的范围内行驶，这与实际不符。Payne本人也发现在高密度情况下，模型可能会遇到稳定性问题，车流密度可能出现大大偏离实际的高密度问题。后来的研究者在payne模型的基础上，不断加入新的项，构成了各自的模型。
1.3Kuhne模型
1984年，R.D.Khune引入交通流的粘性影响，基于Navier-Stokes方程建立如下方程：
（1-8）
式中:为直接与车辆跟驰的弹性有关的声速；为粘性系数。线性稳定性分析表明，当
时交通状态是稳定的。而
（1-9）
式中：称为临界密度。当超过时交通完全瘫痪。该模型可用于超拥挤状态的交通分析但仍需确定平衡状态下的速度——密度关系，因而并未根本解决Payne模型中的致命问题。
1.4吴正模型
我国学者吴正针对我国低速混合交通提出了～维管流模型，引用了交通压力、交通指数等新概念，并通过数值模拟分析了交通堵塞的形成和疏导过程，与实测相符[4