利用matlab对老路纵断面线形拟合

《利用matlab对老路纵断面线形拟合》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、万方数据4中外公路第30卷第5期2O1O年1O月文章编号：1671—2579【2010)05一o004一031前言利用Matlab对老路纵断面线形拟合施金君。程建川(东南大学交通学院，江苏南京210096)摘要：在进行老路改扩建过程中，常因道路修建较早或者有些路段已经经过多次改建扩建，原有设计文件往往不易得到或已与实际情况不符。现阶段研究表明：对老路的纵断面线形进行测量并拟合分析测量数据能有效地获取老路竖曲线参数。目前在中国，一般采用人工拉线法，由于操作起来主观性较强，不但费工时，还不能有效地控制误差。笔者介绍一个模型，

2、二次开发Matlab提供的工具箱，半自动地拟合老路纵断面线形，目标函数为每个测量点与拟合曲线的绝对误差总和最小。最后通过实例验证，获得了理想的结果。关键词：纵断面线形；竖曲线；老路；曲线拟合；Matlab老路改扩建工程中，常需要对已有纵断面线形参数进行准确恢复，以利公路改造(老路基利用、路面洗刨、加铺等)设计更好地满足高程约束。道路纵断面线形的恢复工作常用做法是：测量老路中心线桩号和高程，拟合出原道路线形设计资料，并通过技术手段将是误差控制在允许范围内，为新的道路线形设计提供客观依据。目前中国的一般做法是采用人工拉线法，

3、首先使用全站仪现场测量出一系列道路中线点坐标，再将这些坐标导入AutoCAD软件中，人工穿线定出交点。人工拉线主观性较强，不但工作量大，而且不能有效地控制误差。国外的一般做法为采用函数优化法，把每个测量点的横纵坐标视为决策变量，在LONGO中编程，函数优化目标为拟合出来的曲线与测量点的高程绝对差总和最小。该方法的缺点：每个测量点都视为鼍鬻黼囊鬻崇鬻鬻睾素鬻荣鬻鬻鬻豢鬻鬻鬻{*麝荣鬻豢鬻黼鬻鬻寨鬻豢鬻鬻素豢幕豢睾黼蕞蕞荣幸料荣鬻繁荣豢囊豢豢豢鬻黼黼豢睾料豢鬻黼豢囊料崇荣黼崇等黼豢鬻_琦}靠亲善豢弗睾鬻鬻鬻亲蕞豢鬻黼鬻鬻(3

4、)锚索框架。抗滑桩上部坡面设置锚索框架防减少对原有生态的破坏，保护现有生态植被资源。护，锚索为6索，采用无粘结钢绞线。锚索直接在框架运营期间需防止初期路面水进入水体，沿线设置上张拉，待框架达到设计强度的90％以上后再张拉油水分离池对污水进行处理后再排放，避免污染。锚索。(4)排水措施。坡面外侧设置截水沟，挡墙后设置平台排水沟，汇水排入路基边沟。排水沟采用M10水泥砂浆抹面。3．3错台路基设计。降低土石方开挖及占地对于横坡较陡，地形较顺适路段，采用错台路基设计，尽量减少路基开挖土石方量。3．4对水环境的保护设计重点加强赤水

5、河生态环境保护，减少对沿河线造成水土流失以及对水环境的破坏影响；强化沿河线的生态绿化，促进旅游业发展。合理选择路线走廊带，顺应地形、地貌特点布设路线，注意避让不良地质，使公路立体线形贴近自然，尽量减少工程规模，最大限度地收稿日期：20lO一02—25作者简介：施金君，男，硕士研究生．E—mail：jinbi007@163．com结语沿赤水河路线设计过程中，充分考虑了地形、地质、路线与赤水河的关系以及环保等，不断优化调整路线平、纵、横面设计以及环保设计，使路线建成后为赤水河畔增添一道壮美风景。参考文献：[1]交通部公路司．

6、新理念公路设计指南[M]．北京：人民交通出版社，2005．[2]交通部公路司．降低造价公路设计指南[M]．北京：人民交通出版社，2005．万方数据2010年第5期施金君，等：利用Matlab对老路纵断面线形拟合5变量，编写的程序语言复杂，计算时间长，可视性差。笔者二次开发Matlab中提供的拟合工具箱，编写出一个模型半自动的拟合公路的纵断面测量数据，能较快而又准确地获取纵断面线形参数，包括竖曲线半径、竖曲线的起终点位置和高程等。模型简单，充分利用Matlab提供的可视化平台。下文先介绍拟合的方法和模型的结构，最后给出案例

7、分析和结论。2在Matlab中拟合方法考查典型凸形纵断面。通过现场测量，测得老路中心线纵断面的一系列坐标数据：(z；，y；)，i=1，2，3，⋯，咒。z，为测量点的桩号，弘为高程。假设竖曲点(SQD)和竖直点(SZD)的坐标分别为(z。，y。)、(z2，yz)。模型拟合目标为：当未知的z。和z。在切线上移动时，寻找总误差最小的最佳竖曲线。如图1所示。莉毛勉图l测量点和拟合曲线示意图．对于任意测量点i，拟合估计值y和测量值y；的误差d；为：df=y一弘(1)yn=i1工+口(3)3f12一i2z+6(4)得到两条切线的坡度

8、：i。和i。。然后把z。、z。值分别代入式(3)、(4)，获得了SQD和SZD的坐标为：(z。，z：)、(z2，了2)。于是竖曲线(抛物线)的长度L为：L=z2一zl(5)竖曲线的半径R表达为：R=L／(f2一il)(6)又根据竖曲线抛物线的特征，可推算得抛物线的方程：y=y。+i。(z—z。)+躺(z—z-)2(7