《连续交通流模型及数值模拟》由会员分享，可在线阅读，更多相关《连续交通流模型及数值模拟（33页珍藏版）》请在金锄头文库上搜索。
1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-date连续交通流模型及数值模拟连续交通流模型及数值模拟连续交通流模型及数值模拟摘要本文对现有的交通流宏观模型进行了研究，总结了各种模型的思想、优缺点以及适用条件，在此基础上，选取了Payne模型离散格式进行数值模拟，选取了某段高速公路的交通流作为模拟对象，展现了Payne模型模拟交通流的可行性。 关键字 连续交通流；离散格式；数值模拟0 引言交通流理论研究加深了人们对复杂多
2、体系统远离平衡态时演变规律的认识，促进了统计物理、非线性动力学、应用数学、流体力学、交通工程学等学科的交叉和发展等多学科的交叉渗透和相互发展。交通流理论研究的对象是离散态物质，是一个复杂的非线性体系，对这类物质运动规律的描述，尚无成熟的理论。 在宏观的连续流模型中，交通流被比拟为连续的流体介质，即将流量、速度和密度等集聚变量视为时间和空间的连续函数。模型包含时间和空间的状态方程，考虑了车辆的加速度、惯性和可压缩性，能够合理准确描述交通流的动态特性，相比微观模型有更大的优势。连续流交通流模型通常用密度（k）、速度（u）、流量（q）三个变量来描述1。1 连续交通流模型1.1 LWR模型1955年，
3、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个方程构成
4、了完整的一阶连续交通流模型，LWR模型的优点是简单明了，可以采用流体力学和应用数学中的成熟工具进行分析，而且可以描述诸如交通阻塞形成和消散之类的交通现象，但是，由于该模型的速度是由平衡速度密度关系决定，并且没有考虑加速度和惯性影响，因此不适用于描述本质上处于非平衡态的交通现象，例如车辆上、下匝道的交通、“幽灵式”交通阻塞、交通迟滞、时走时停的交通等。于是，后来的学者们引进了高阶连续介质模型，考虑了加速度和惯性影响，将动量方程代替方程(1-3)。1.2 Payne模型Pipes于1953年提出交通流加速度的一般表达式： （1-4）1971年，Payne根据LWR模型的思想，假设交通流速度是动态变
5、化的，在引用连续性方程时，引进运动方程，导出高阶连续模型3。Payne从车辆跟驰理论的概念提出平均速度u与密度k存在以下关系： （1-5）并取，车辆跟驰理论在的延滞时间，对上式分别作关于和的Taylor展开，得到： （1-6）上式中，：预期指数，将平衡速度简写成。从而建立了由如下三个方程构成的Payne模型： （1-6） （1-7） （1-2）式(16)的右边第一项为期望项，反映驾驶员对前方交通状态改变的反应过程；第二项式弛豫项，描述车流速度在时间内向平衡速度的调整，最优速度函数和其他参数通过道路实测和参数辨识确定。1979年，Payne编制了著名的FREFLO软件，有史以来第一次将交通流仿真
6、模型应用于工程实践。但Payne模型并未充分考虑整个弛豫过程，而只是将其定为一个常数的弛豫时间，即使处于平衡状态时，弛豫时间变为零，在实际应用中，出现了一些问题。Rathi等人指出，使用Payne模型，车流速度到平衡态速度的调节过程过于缓慢。Ross也发现，当道路拓扑特性和交通量在短时间内突变时，由于车流速度调整到平衡速度过程缓慢，无法捕捉到实际交通流动态特性。Castillo等对Payne模型进行了线性稳定性分析，发现车辆总是在稳定的范围内行驶，这与实际不符。Payne本人也发现在高密度情况下，模型可能会遇到稳定性问题，车流密度可能出现大大偏离实际的高密度问题。后来的研究者在payne模型的
7、基础上，不断加入新的项，构成了各自的模型。1.3 Kuhne模型1984年，R.D.Khune引入交通流的粘性影响，基于Navier-Stokes方程建立如下方程： （1-8）式中: 为直接与车辆跟驰的弹性有关的声速；为粘性系数。线性稳定性分析表明，当时交通状态是稳定的。而 （1-9）式中：称为临界密度。当超过时交通完全瘫痪。该模型可用于超拥挤状态的交通分析但仍需确定平衡状态下的速度密度关系，因而并未根本解决Payne 模型中的致命问题。1.4 吴正模型我国学者吴正针对我国低速混合交通提出了维管流模型，引用了交通压力、交通指数等新概念，并通过数值模拟分析了交通堵塞的形成和疏导过程，与实测相符4
8、。一维管道内物质流的基本方程如下： （1-10） （1-11）其中变量k车流密度，即单位长度路段上的车辆数；A路段宽度或车道数；u车流速度； 车流经过单位面积路面时所受阻力；P交通压力。提出了交通流压力比拟的基本假定： （） （1-12）其中C和n是两个可调参数，改变它们的取值就能使模型适用于不同的交通情况。n称为交通状态指数。吴正模型可以较好地模拟实际交通过程，包括局部地方发生严重交通堵塞后的疏散过程。东明等人在吴正的运动方程中考虑正比于速度的阻尼项的作用，数值模拟分析了地面交通对高架路交通的影响。2 交通流离散模型数值模拟2. 1 交通流计算模型笔者采用的交通流计算模型是Payne型模型
9、其运动微分方程为： （2-1） （2-2） （2-3） 2. 2 模型的离散格式采用Euler积分对Payne运动微分方程进行离散，得到Payne模型的离散方程： （2-4） （2-5） （2-6）式中：n=1,2,3,为时间步；i=1,2,3,为空间节点位置；为第i段道路长度；为进出匝道的流率。本文采用的时间步长、空间步长随时间序列和空间节点位置两者的变化而变化的变步长。令 （2-7） （2-8）其中（）为畅行速度。3 算例及计算结果分析以广深高速公路某时段的实测交通数据为样本进行数值模拟。计算时边界条件均采用自由流边界条件：， （3-1） 算例的计算路段长度L= 6 km 设计速度 =12
10、0km/h 阻塞密度= 200pcu/ ( km车道)最大波速选自广深高速公路2车道的实测数据回归值 = =100km/h。计算中，将路段L等距分成12个区段 即= 0.5 km， L= 12，模拟100步时间步，初值条件如表3-1所列。表3-1 变步长计算初值条件区段123456789101112密度/pcu(km车道)-117.918.418.418.218.152.151.851.718.618.117.818.1速度/(kmh-1)87.787.687.687.687.744.945.145.287.687.787.787.7经过100步迭代计算后密度和速度的计算结果如图3-1、图3-
11、2 所示，部分结果如表3-2所示。图3-1 变步长密度关系图图3-2 变步长速度关系图表3-2 100步迭代部分结果区段计算步长97步98步99步100步密度速度密度速度密度速度密度速度117.9109.2117.9109.2117.9109.2117.9109.21217.9109.2117.9109.2117.9109.2117.9109.21317.9109.2117.9109.2117.9109.2117.9109.21417.9109.2117.9109.2117.9109.2117.9109.21518108.5918108.5918108.5918108.59625.01108.
12、1625.01108.1625.01108.1625.01108.16725.1107.7825.1107.7825.1107.7825.1107.78825.07107.925.07107.925.07107.925.07107.9918.56108.0418.56108.0418.56108.0418.56108.041018.54108.1618.54108.1618.54108.1618.54108.161118.51108.2918.52108.2918.51108.2918.52108.291218.25109.9618.22109.9618.25109.9618.22109.964 结论本文首先对现有的交通流宏观模型进行了研究，总结了各种模型的思想、优缺点以及适用条件，在此基础上，选取了Payne模型离散格式进行数值模拟，选取了某段高速公路的交通流作为模拟对象，展现了Payne模型模拟交通流的可行性。参考文献：1李进平,张乐文.交通流的流体动力学模型研究,武汉理工大学学报(交通科学与
