1、 纵断面优化的大致思路基本相似，即纵断面优化的大致思路基本相似，即使得工程量（工程费用）、工程和运营总使得工程量（工程费用）、工程和运营总费用最小！在新线设计中，主要使土石方费用最小！在新线设计中，主要使土石方数量最小，在既有线中，主要使起落道量数量最小，在既有线中，主要使起落道量最小。从而达到工程费用最小的目的。最小。从而达到工程费用最小的目的。在本章中我们仅以既有线为例来说明在本章中我们仅以既有线为例来说明优化的思路和过程，优化的原理主要是怎优化的思路和过程，优化的原理主要是怎么样使得目标函数最小。所以我们在此仅么样使得目标函数最小。所以我们在此仅以最小二乘法结合黄金分割法来实现目标以最小
2、二乘法结合黄金分割法来实现目标函数的最小值，从而实现优化的目的。函数的最小值，从而实现优化的目的。既有线在运营过程中，个别路段的路基既有线在运营过程中，个别路段的路基会因沉陷、冻害而变形，在经常维修过程中，会因沉陷、冻害而变形，在经常维修过程中，由于更换道碴、起道、落道，也要引起轨面由于更换道碴、起道、落道，也要引起轨面标高的改变。所以既有线轨面的纵断面多与标高的改变。所以既有线轨面的纵断面多与原设计不同，而原设计标准又多偏低，不符原设计不同，而原设计标准又多偏低，不符合现行合现行线规线规标准；延长站线而需加长站标准；延长站线而需加长站坪长度时，引起站坪两端纵断面的改建；削坪长度时，引起站坪两
3、端纵断面的改建；削减超限坡时，需要抬高或降低路基标高；线减超限坡时，需要抬高或降低路基标高；线路受洪水威胁地段，则需加高路基。这一切路受洪水威胁地段，则需加高路基。这一切都要引起线路纵断面的改建。都要引起线路纵断面的改建。改建既有纵断面设计，以线路纵断面测量出的轨改建既有纵断面设计，以线路纵断面测量出的轨面高程为准。面高程为准。一般情况下，起道高程小于一般情况下，起道高程小于50cm时，用道碴起道；时，用道碴起道；等于等于50100cm时，用渗水土壤起道；大于时，用渗水土壤起道；大于100cm或落或落道后道床厚度小于规定标准时，需抬降路基面。为了方道后道床厚度小于规定标准时，需抬降路基面。为了
4、方便施工及减轻对运营的干扰，一般不采用挖切路基的办便施工及减轻对运营的干扰，一般不采用挖切路基的办法来降低轨面高程，仅在受建筑限界与结构物构造控制，法来降低轨面高程，仅在受建筑限界与结构物构造控制，以及为消除路基病害的地段方可采用。亦不宜降低既有以及为消除路基病害的地段方可采用。亦不宜降低既有线轨面高程，以免挖切道床影响正常运营，仅在个别地线轨面高程，以免挖切道床影响正常运营，仅在个别地段，为避免改建桥隧建筑物，避免挖切路基，或为了减段，为避免改建桥隧建筑物，避免挖切路基，或为了减少线路改建工程，才允许挖切道床以降低轨面高程。道少线路改建工程，才允许挖切道床以降低轨面高程。道床厚度仅允许较规定
5、标准减薄床厚度仅允许较规定标准减薄5cm以内，但最小道床厚以内，但最小道床厚度不得小于度不得小于25cm。既有线纵断面设计，要求细致准确，以保证充分利既有线纵断面设计，要求细致准确，以保证充分利用既有线建筑物，减少改建工程，故应采用距离为用既有线建筑物，减少改建工程，故应采用距离为1：10000，高程为，高程为1：100的放大纵断面图进行设计。放大的放大纵断面图进行设计。放大纵断面图的下部，自下而上设有既有线平面、百米标与纵断面图的下部，自下而上设有既有线平面、百米标与加标、地面高程、既有道床厚度、既有轨面高程、轨面加标、地面高程、既有道床厚度、既有轨面高程、轨面设计坡度、轨面设计高程、设计与
6、既有轨面高程差（抬设计坡度、轨面设计高程、设计与既有轨面高程差（抬降值）、病害、工程地质特征等栏。放大纵断面图的上降值）、病害、工程地质特征等栏。放大纵断面图的上半部绘有地面线、既有道床底面线、既有轨面线、设计半部绘有地面线、既有道床底面线、既有轨面线、设计轨面线、计算轨面线，并标明建筑物的特征，如车站、轨面线、计算轨面线，并标明建筑物的特征，如车站、道口的中心里程，隧道洞门里程及长度，以及桥涵类型、道口的中心里程，隧道洞门里程及长度，以及桥涵类型、孔径、中心里程及设计洪水位高程等。孔径、中心里程及设计洪水位高程等。（一）竖曲线（一）竖曲线 竖曲线的坡度代数差，一、二级铁路差值超过竖曲线的坡度
7、代数差，一、二级铁路差值超过3，三级铁路差值超过，三级铁路差值超过4时需设竖曲线，且竖曲时需设竖曲线，且竖曲线不应与缓和曲线重合，竖曲线不应设在明桥面上，线不应与缓和曲线重合，竖曲线不应设在明桥面上，竖曲线不应与道岔重叠。竖曲线不应与道岔重叠。（二）最大坡度的折减（二）最大坡度的折减 由于平面上出现曲线和遇到长于由于平面上出现曲线和遇到长于400米的隧道米的隧道时，因为附加阻力增到、粘着系数降低，而需将最时，因为附加阻力增到、粘着系数降低，而需将最大坡度值减缓。具体折减方法在资料中可查。大坡度值减缓。具体折减方法在资料中可查。（三）桥涵（三）桥涵 有碴桥涵梁上，一般应该按计算轨面高程设计纵有碴
8、桥涵梁上，一般应该按计算轨面高程设计纵断面。通常不允许落低既有线高程，以免因降低墩台顶断面。通常不允许落低既有线高程，以免因降低墩台顶面高程引起施工困难，此时需加高粱的边墙，以免道碴面高程引起施工困难，此时需加高粱的边墙，以免道碴溢出。轨面抬高值一般在溢出。轨面抬高值一般在1015cm,以免加厚道碴后影以免加厚道碴后影响桥梁的应力与稳定性。当抬高值较大时，则需要加高响桥梁的应力与稳定性。当抬高值较大时，则需要加高墩台顶面高程，施工困难。墩台顶面加高值大于墩台顶面高程，施工困难。墩台顶面加高值大于0.4m时时应进行强度和稳定性检算。应进行强度和稳定性检算。明桥面桥梁上，轨面高程的变动必将引起抬降
9、墩明桥面桥梁上，轨面高程的变动必将引起抬降墩台顶面高程等困难高程，因此应根据既有轨面高程设计台顶面高程等困难高程，因此应根据既有轨面高程设计纵断面。纵断面。涵洞处所允许适当抬降既有轨面高程。但抬降值涵洞处所允许适当抬降既有轨面高程。但抬降值过大时，往往需要改建涵洞的端墙与翼墙，甚至接长涵过大时，往往需要改建涵洞的端墙与翼墙，甚至接长涵洞。若大量降低既有轨面高程并挖低路基时，应保证涵洞。若大量降低既有轨面高程并挖低路基时，应保证涵洞顶部到道床底面的最小填土高程。洞顶部到道床底面的最小填土高程。（四）隧道（四）隧道 隧道内需要提高隧道净空或削减隧道内的坡度时，隧道内需要提高隧道净空或削减隧道内的坡
10、度时，一般采用落道方法，以免破坏隧道的拱圈，但降低值不一般采用落道方法，以免破坏隧道的拱圈，但降低值不宜大于宜大于0.4m，以免破坏隧道边墙基础。，以免破坏隧道边墙基础。（五）车站站坪（五）车站站坪 车站站坪内正线的纵断面一般不宜过多抬降，以免车站站坪内正线的纵断面一般不宜过多抬降，以免引起站内建筑物（如车站站线、咽喉区、站台、天桥信引起站内建筑物（如车站站线、咽喉区、站台、天桥信号与给水等设备）的改建。号与给水等设备）的改建。当减缓站坪坡度、延长站坪长度、增设车站、或消当减缓站坪坡度、延长站坪长度、增设车站、或消减限制坡度引起站坪纵断面改建时，应该全面考虑，使减限制坡度引起站坪纵断面改建时，
11、应该全面考虑，使整个改建工程量最小。站坪车道上坡度一般不应大于整个改建工程量最小。站坪车道上坡度一般不应大于2.5。站坪长度由到发线有效长度和线路形式决定。站坪长度由到发线有效长度和线路形式决定。车站正线抬降时，可用站线作为施工临时通车线车站正线抬降时，可用站线作为施工临时通车线。（六）路基（六）路基 挡土墙、护坡地段抬道时，应考虑加宽路基后，挡土墙、护坡地段抬道时，应考虑加宽路基后，不使其填土坡脚盖过挡土墙或护坡，必要时，可用干不使其填土坡脚盖过挡土墙或护坡，必要时，可用干砌片石加陡坡度。砌片石加陡坡度。路基病害地段，如沙害、雪害，以及因毛细水上路基病害地段，如沙害、雪害，以及因毛细水上升引
12、起的冻害或翻浆冒泥地段均应结合抬道整治病害。升引起的冻害或翻浆冒泥地段均应结合抬道整治病害。路基基床土质不良及道床排水不畅引起的道碴陷路基基床土质不良及道床排水不畅引起的道碴陷囊，一般可结合落道给予整治。囊，一般可结合落道给予整治。路堑地段落道时，应考虑施工扩堑对行车的干扰，路堑地段落道时，应考虑施工扩堑对行车的干扰，特别是石质路堑需要爆破施工的地段，干扰更为严重。特别是石质路堑需要爆破施工的地段，干扰更为严重。此外应考虑路堑边坡的稳定与地下水位的影响。如设此外应考虑路堑边坡的稳定与地下水位的影响。如设计的路肩高程低于地下水位，应采用降低地下水位的计的路肩高程低于地下水位，应采用降低地下水位的
13、措施。措施。（一）根据外业勘查资料，填写既有线平面、一）根据外业勘查资料，填写既有线平面、百米桩与加标、地面高程、既有道床厚度及既有轨百米桩与加标、地面高程、既有道床厚度及既有轨面高程各栏数据，并标明路基病害地段与工程地质面高程各栏数据，并标明路基病害地段与工程地质特征。特征。地面高程一般按线路前进方向左侧既有线路堑地面高程一般按线路前进方向左侧既有线路堑坡顶或路堤坡脚点的高程填写，桥涵处按实际沟底坡顶或路堤坡脚点的高程填写，桥涵处按实际沟底高程填写；隧道处按地形图高程填写。高程填写；隧道处按地形图高程填写。（二）道床底面高程和既有轨面高程的计算：（二）道床底面高程和既有轨面高程的计算：道床底
14、面高程道床底面高程=既有轨面高程既有轨面高程-既有轨道高度既有轨道高度 计算轨面高程计算轨面高程=道床底面高程道床底面高程+设计轨道高度设计轨道高度 轨道高度轨道高度=钢轨高度钢轨高度+垫板厚度垫板厚度+轨枕高度轨枕高度+道床厚度道床厚度 （三）根据既有轨面高程绘出既有轨面线，根（三）根据既有轨面高程绘出既有轨面线，根据道床底面高程和计算轨面高程绘出道床底面线和据道床底面高程和计算轨面高程绘出道床底面线和计算轨面线。标明车站、道口的中心里程与长度，计算轨面线。标明车站、道口的中心里程与长度，以及桥涵类型、孔径与中心里程。以及桥涵类型、孔径与中心里程。（四）设计轨面线应符合设计标准，并使其尽（四
15、）设计轨面线应符合设计标准，并使其尽量接近且不低于既有轨面线与计算轨面线，以宁抬量接近且不低于既有轨面线与计算轨面线，以宁抬勿降的原则进行坡度设计。若设计轨面线高出计算勿降的原则进行坡度设计。若设计轨面线高出计算轨面线过多，则垫铺的道碴过多，引起浪费。轨面线过多，则垫铺的道碴过多，引起浪费。（五）根据初步绘出的设计轨面线定出轨面设（五）根据初步绘出的设计轨面线定出轨面设计坡度（取计坡度（取0.10.1的倍数）的倍数）.坡段长应为坡段长应为10m10m的倍数。的倍数。根据轨面设计坡度，标出各百米标与设计轨面高程根据轨面设计坡度，标出各百米标与设计轨面高程（准确到厘米）。变坡点的坡度差（准确到厘米
16、）。变坡点的坡度差 i i33时，应时，应设置圆曲线形竖曲线。并计算竖曲线要素。设置圆曲线形竖曲线。并计算竖曲线要素。（六）计算既有轨面抬降值，即设计轨面高程（六）计算既有轨面抬降值，即设计轨面高程与既有轨面高程之差。竖曲线应按设置竖曲线后的与既有轨面高程之差。竖曲线应按设置竖曲线后的设计轨面高程来计算抬降值设计轨面高程来计算抬降值.路堤地段，特别是高路堤地段抬道时，应考路堤地段，特别是高路堤地段抬道时，应考虑加宽路基稳定的影响。虑加宽路基稳定的影响。在本设计中，放大纵断面图的设计，是先在在本设计中，放大纵断面图的设计，是先在EXCELEXCEL表中设计轨面坡度，然后汇总各项设计资料表中设计轨
17、面坡度，然后汇总各项设计资料绘制详细断面图。设计方法与步骤如下：绘制详细断面图。设计方法与步骤如下：(一一)根据外业勘查资料，在根据外业勘查资料，在EXCELEXCEL表中填写里表中填写里程、既有轨面标高、既有道床厚三列数据。程、既有轨面标高、既有道床厚三列数据。三、用三、用ExcelExcel表格进行设计坡度计算表格进行设计坡度计算 在设计坡度计算中以放大纵断面设计各项要求，在设计坡度计算中以放大纵断面设计各项要求，注意事项来设计。反复调整设计坡度使其充分满足注意事项来设计。反复调整设计坡度使其充分满足设计要求。在设计坡度时一般要考虑：站坪的位置设计要求。在设计坡度时一般要考虑：站坪的位置（
18、坡度（坡度2.52.5，坡段长度，坡段长度1300m1300m）；平曲线的位）；平曲线的位置（竖曲线不能合缓和曲线重合）；坡段长度置（竖曲线不能合缓和曲线重合）；坡段长度（200m200m）；桥涵的位置（抬降量的限制）等。）；桥涵的位置（抬降量的限制）等。图中道床厚度单位为厘米所以计算时要单位统一转化。图中道床厚度单位为厘米所以计算时要单位统一转化。0.3640.364为：既有钢轨高度为：既有钢轨高度+垫板厚度垫板厚度+轨枕高度值。轨枕高度值。图中图中0.8380.838为：设计钢轨高度为：设计钢轨高度+垫板厚度垫板厚度+轨枕高轨枕高度值度值+道床厚度道床厚度 外包线是以外包线是以宁抬勿降宁抬
19、勿降的原则来选择既有和设计轨的原则来选择既有和设计轨面标高而计算。面标高而计算。设计轨面高程是在给定一个设计坡度后而得到的计算设计轨面高程是在给定一个设计坡度后而得到的计算值。注意设计起点不变，即起点的设计轨面高等于既有值。注意设计起点不变，即起点的设计轨面高等于既有轨面高。图中轨面高。图中1877.761877.76也是既有轨面高。设计坡度时注也是既有轨面高。设计坡度时注意上面所提到的设计要求。意上面所提到的设计要求。路肩高程路肩高程=道床底面标高道床底面与路肩之间的道床底面标高道床底面与路肩之间的垂直距离（垂直距离（10.89cm10.89cm）道床底面与路肩之间的垂直距离可以通过道床底面
20、与路肩之间的垂直距离可以通过CADCAD作出作出设计线标准横断面图而量出此值。设计线标准横断面图而量出此值。1.1.计算坐标值计算坐标值 由于由于CADCAD中长度单位为毫米，图中单位为米，中长度单位为毫米，图中单位为米，横向（里程）采用横向（里程）采用1 1：1000010000的比例，竖向（高程）的比例，竖向（高程）采用采用1 1：200200的比例，所以只需给表格中的里程数据的比例，所以只需给表格中的里程数据除以除以1010即为即为1 1：1000010000的比例，高程数据乘以的比例，高程数据乘以5 5即可即可为为1 1：200200的比例，中间的的比例，中间的&,&,&符号可以将以上
21、两符号可以将以上两数输在同一个单元格中。这样可分别列出既有线、数输在同一个单元格中。这样可分别列出既有线、设计线和计算线的数据。设计线和计算线的数据。图中的公式中：图中的公式中：-168800-168800，-7000-7000的原因是，坐的原因是，坐标值太大。同时减去一个定值是坐标值变小，但整标值太大。同时减去一个定值是坐标值变小，但整体线与线的关系不变。在体线与线的关系不变。在CADCAD中选择多段线，再将中选择多段线，再将坐标数据复制粘贴到指令栏中即可。设计线坐标只坐标数据复制粘贴到指令栏中即可。设计线坐标只取变坡点坐标来画出设计线路。取变坡点坐标来画出设计线路。2 2所需数据的调整：由
22、于此图从右开始向左所需数据的调整：由于此图从右开始向左绘制，数据的方向需要在图中向左，所以必须将数绘制，数据的方向需要在图中向左，所以必须将数据颠倒。据颠倒。3 3抬降量分列：在图中抬降量各有一行，而抬降量分列：在图中抬降量各有一行，而表格中只有一列，必须将抬降量分开。在表格中用表格中只有一列，必须将抬降量分开。在表格中用判断语句可将其分开。判断语句可将其分开。4.CAD4.CAD中数据的调整中数据的调整 粘贴到图中的数据是一个整体，可以通过调整粘贴到图中的数据是一个整体，可以通过调整其列间距（其列间距（1010），角度（），角度（9090），字高（），字高（2.52.5），），宽度（宽度（1
23、010），对齐方式（中下）既可以准确地放入），对齐方式（中下）既可以准确地放入数据。数据。在图中数据有些是保留两位小数，为了数据的准在图中数据有些是保留两位小数，为了数据的准确必须在确必须在ExcelExcel表格中也将相对应的数据计算保留两位表格中也将相对应的数据计算保留两位小数。小数。注意：如上图所示有竖曲线的地方修正标高注意：如上图所示有竖曲线的地方修正标高1824.67=1824.67=设计标高（设计标高（1824.591824.59）+EO+EO（0.080.08）在图中需要画的在图中需要画的4 4条线是：既有线，设计线，条线是：既有线，设计线，计算线，道床底面线。计算线，道床底面线
24、。既有轨面线：一般直线既有轨面线：一般直线 设计轨面线：宽度为设计轨面线：宽度为0.7mm0.7mm 计算轨面线：点划线（比例计算轨面线：点划线（比例0.50.5）道床底面线：虚线（比例道床底面线：虚线（比例0.70.7）一、纵断面程序设计过程一、纵断面程序设计过程 （1 1）做数据文件，读入数据，包括既有轨面标高、）做数据文件，读入数据，包括既有轨面标高、道床厚度、平曲线要素、限制坡度等。道床厚度、平曲线要素、限制坡度等。（2 2）判断变坡点位置（二阶差分法）。）判断变坡点位置（二阶差分法）。（3 3）求单个坡段最优坡度值（黄金分割结合最小）求单个坡段最优坡度值（黄金分割结合最小二乘）。二乘
25、）。（4 4）变坡点处的修正（位置，标高）。）变坡点处的修正（位置，标高）。（5 5）曲线折减检核。）曲线折减检核。（6 6）控制点检核。）控制点检核。二、不考虑约束问题时坡度拟合二、不考虑约束问题时坡度拟合如下图所示：如下图所示：未拟合坡度线既有坡度线拟合范围HHikiIN2N1OXY图4-5 坡度拟合图对上式求导，并使对上式求导，并使 0Isdd0)(2iiiisLDLIdd0112iiiiiiiLDIL121iiiiiiiLLDImin2CBi（二）最优初始点的选择（用黄金分割法）（二）最优初始点的选择（用黄金分割法）铁路规范：降量铁路规范：降量5cm5cm，抬升量，抬升量50cm50c
26、m在在-0.05-0.05，+0.5+0.5范围内时，可以计算出最优值范围内时，可以计算出最优值I I。由上分析易知，初始点的标高不同，最后所得由上分析易知，初始点的标高不同，最后所得的坡度值也不相同，所以我们要采用优化方法来找的坡度值也不相同，所以我们要采用优化方法来找出最优的初始可行点，在这里，我们采用黄金分割出最优的初始可行点，在这里，我们采用黄金分割法来完成。如图所示：法来完成。如图所示：1.选定初始点（按实测轨面标高）定分选定初始点（按实测轨面标高）定分割区间割区间（D=0.2）baHH,05.0iaHH50.0ibHH 建立（建立（0.618）法单谷函数曲线（四个点）法单谷函数曲线
27、（四个点):bbaaSSSS,HaHHby(H)xOS0.618 DHaHbHaHbSaSbSaSbH0.382 DHaHHby(H)xOS0.618 DHaHbHaHbSaSbSaSbH0.382 D 取消目标函数取消目标函数S S最大点，重新定出最大点，重新定出单谷函数曲线单谷函数曲线 判断判断 ，定出最优初始点高程定出最优初始点高程baSS3.3.变坡点判断变坡点判断 直线范围内：直线范围内：kxy 求导后求导后 0y故可用二阶导数判断变坡点，利用差分法求二阶导数：故可用二阶导数判断变坡点，利用差分法求二阶导数：SSyyyynnnn112式中：式中：NN点号（点号（5050米整倍米整倍数
28、）数）SS相邻点的里程差相邻点的里程差判断如果判断如果 112.12.1nnn则则n n点为变坡点点为变坡点 可用可用坡度差进行计算坡度差进行计算 4.4.变坡点修正变坡点修正 利用上述方法，分别对各个坡段进行优化，优利用上述方法，分别对各个坡段进行优化，优化后在坡段的交点处标高不一定相同，需要对其进化后在坡段的交点处标高不一定相同，需要对其进行整饰。我们采用以下方法：行整饰。我们采用以下方法：5.5.设计坡度的计算设计坡度的计算11iiiiLHHI三、三、考虑约束时算法考虑约束时算法 1.1.高程约束高程约束 限制起道量（隧道）和落道量（桥梁）的地段作为限制起道量（隧道）和落道量（桥梁）的地
29、段作为高程约束点，以无约束最优初始点和最优坡度为中值，对高程约束点，以无约束最优初始点和最优坡度为中值，对称平移或旋转坡度线，直到达到要求为止。称平移或旋转坡度线，直到达到要求为止。2.2.限坡约束限坡约束iix3.3.平面约束平面约束 竖缓、竖岔不能重合，重合时移动变坡点的位竖缓、竖岔不能重合，重合时移动变坡点的位置，移动方向以改变纵断面较小为准。所以，必须置，移动方向以改变纵断面较小为准。所以，必须计算平面曲线的坡度折减值。计算平面曲线的坡度折减值。平面曲线的坡度折减在直缓、缓直之间进行。平面曲线的坡度折减在直缓、缓直之间进行。开始开始输入数据输入数据 角度转化角度转化判断变坡点判断变坡点
30、变坡点高程优化变坡点高程优化各坡段坡度优化各坡段坡度优化变坡点修正变坡点修正圆曲线坡度折减计算圆曲线坡度折减计算实际坡度的计算实际坡度的计算 输出数据输出数据结束结束设计坡度的计算设计坡度的计算抬降量的计算抬降量的计算一、输入数据文件 读入测点数，高程，里程读入测点数，高程，里程 read(1,read(1,*)N)N do 5,i=1,N do 5,i=1,N read(1,read(1,*)H(i),LC(i)H(i),LC(i)continue continue二、计算每两个测点之间的坡度值二、计算每两个测点之间的坡度值 do 10,i=1,Ndo 10,i=1,N PD(I)=(H(I
31、+1)-H(I)/50.0 PD(I)=(H(I+1)-H(I)/50.010 continue10 continue 三、坡度差计算三、坡度差计算 do 20,i=2,Ndo 20,i=2,N PDC(I)=PD(I)-PD(I-1)PDC(I)=PD(I)-PD(I-1)20 continue20 continue四、计算变坡点序号四、计算变坡点序号 J=2 do 30,i=5,N-4 IF(PDC(I).GT.1.2*PDC(I-1).AND.PDC(I).GT.1.2*PDC(I+1)THEN BPD(J)=I J=J+1 END IF30 continue BPD(1)=1 BPD(
32、J)=N五、选择最优初始点高程（五、选择最优初始点高程（H值确定，采用值确定，采用黄金分割法）黄金分割法）DO 40 I=1,J-1DO 40 I=1,J-1 JB=BPD(I)JB=BPD(I)第一变坡点，即初始点第一变坡点，即初始点 JA=BPD(I+1)JA=BPD(I+1)第二变坡点第二变坡点 LOW=H(JB)-0.05LOW=H(JB)-0.05 HIGH=H(JB)+0.5 HIGH=H(JB)+0.5 H1=LOW+0.618H1=LOW+0.618*(HIGH-LOW)(HIGH-LOW)H2=HIGH-0.618 H2=HIGH-0.618*(HIGH-LOW)(HIGH-
33、LOW)LBC(I)=LC(JA)-LC(JB)LBC(I)=LC(JA)-LC(JB)坡度长坡度长度（第一二变坡点之间的）度（第一二变坡点之间的）IF(HIGH-LOW.GT.1E-3)THENIF(HIGH-LOW.GT.1E-3)THEN DL1=0 DL1=0 进行累加前准备进行累加前准备 DL2=0DL2=0 LL=0 LL=0 S1=0 S1=0 S2=0 S2=0 DO 50,II=JB,JADO 50,II=JB,JA LCC=LC(II)-LC(JB)LCC=LC(II)-LC(JB)任意一测点到任意一测点到第一个变坡点之间的距离第一个变坡点之间的距离 DD1=H(II)-H
34、1 DD1=H(II)-H1 相当于相当于 值值 DD2=H(II)-H2DD2=H(II)-H2 DL1=DL1+DD1 DL1=DL1+DD1*LCC dLCC d和和l l的乘积，的乘积，DL2=DL2+DD2DL2=DL2+DD2*LCCLCC LL=LL+LCC LL=LL+LCC*2 l2 l的平方，计算的平方，计算函数函数s s须首先计算坡度须首先计算坡度i i 50 continue50 continue AI1=DL1/LL AI1=DL1/LL 坡度坡度1 1 AI2=DL2/LL AI2=DL2/LL 坡度坡度2 2 do 60,III=JB,JA do 60,III=J
35、B,JA LCC=LC(III)-LC(JB)LCC=LC(III)-LC(JB)S1=S1+(LCC S1=S1+(LCC*AI1-(H(III)-H1)AI1-(H(III)-H1)*2 2 S2=S2+(LCC S2=S2+(LCC*AI2-(H(III)-H2)AI2-(H(III)-H2)*2 2continuecontinue 利用黄金分割法改变利用黄金分割法改变H值，寻找最优坡度值，寻找最优坡度时的时的H值值 IF(S1.GT.S2)THENIF(S1.GT.S2)THEN HIGH=H1 HIGH=H1 H1=H2 H1=H2 H2=HIGH-0.618 H2=HIGH-0.6
36、18*(HIGH-LOW)(HIGH-LOW)ELSE ELSE LOW=H2 LOW=H2 H2=H1 H2=H1 H1=LOW+0.618 H1=LOW+0.618*(HIGH-LOW)(HIGH-LOW)END IF END IF 这只是这只是1 1次逼近，还需进一步逼近，这里次逼近，还需进一步逼近，这里需要用到一个循环，逐步逼近。留作大家思需要用到一个循环，逐步逼近。留作大家思考！考！以最优坡度时的以最优坡度时的H H计算最小的计算最小的S S IF(S1.GT.S2)THEN IF(S1.GT.S2)THEN AIL(I)=AI2 AIL(I)=AI2 HBX(I)=H2 HBX(I
37、)=H2 ELSE ELSE AIL(I)=AI1 AIL(I)=AI1 HBX(I)=H1 HBX(I)=H1 END IF END IFDO 180,I=2,JDO 180,I=2,J HBY(I)=HBX(I-1)+LBC(I-1)HBY(I)=HBX(I-1)+LBC(I-1)*AIL(I-AIL(I-1)1)将每个坡段赋予两个高程（将每个坡段赋予两个高程（HBXHBX前点，前点，HBYHBY后点）便于每个坡度段的修正后点）便于每个坡度段的修正180 CONTINUE180 CONTINUE DO 130,I=2,J-1DO 130,I=2,J-1 IF(ABS(HBY(I)-HBX(
38、I).LE.5.0)THEN IF(ABS(HBY(I)-HBX(I).LE.5.0)THEN HBX(I)=(HBY(I)+HBX(I)/2.0 HBX(I)=(HBY(I)+HBX(I)/2.0 HBY(I)=HBX(I)HBY(I)=HBX(I)AIL(I-1)=(HBY(I)-HBX(I-1)/LBC(I-1)AIL(I-1)=(HBY(I)-HBX(I-1)/LBC(I-1)AIL(I)=(HBY(I+1)-HBX(I)/LBC(I)AIL(I)=(HBY(I+1)-HBX(I)/LBC(I)END IF END IFIF(ABS(HBY(I)-HBX(I).GT.5.0.AND.A
39、IL(I-IF(ABS(HBY(I)-HBX(I).GT.5.0.AND.AIL(I-1).GT.AIL(I).AND.HBY(I).GT.HBX(I)THEN1).GT.AIL(I).AND.HBY(I).GT.HBX(I)THEN BPD(I)=BPD(I)-1 BPD(I)=BPD(I)-1 HBX(I)=HBY(I)-50 HBX(I)=HBY(I)-50*AIL(I-1)AIL(I-1)HBY(I)=HBX(I)HBY(I)=HBX(I)LBC(I-1)=LBC(I-1)-50 LBC(I-1)=LBC(I-1)-50 LBC(I)=LBC(I)+50 LBC(I)=LBC(I)+5
40、0 AIL(I)=(HBY(I+1)-HBX(I)/LBC(I)AIL(I)=(HBY(I+1)-HBX(I)/LBC(I)END IF END IF 50mHBY(i)HBX(i)新的BPDIF(ABS(HBY(I)-HBX(I).GT.5.0.AND.AIL(I-IF(ABS(HBY(I)-HBX(I).GT.5.0.AND.AIL(I-1).GT.AIL(I).AND.HBY(I).LT.HBX(I)THEN1).GT.AIL(I).AND.HBY(I).LT.HBX(I)THEN BPD(I)=BPD(I)+1 BPD(I)=BPD(I)+1 HBX(I)=HBX(I)-50 HBX(
41、I)=HBX(I)-50*AIL(I)AIL(I)HBY(I)=HBX(I)HBY(I)=HBX(I)LBC(I)=LBC(I)-50 LBC(I)=LBC(I)-50 LBC(I-1)=LBC(I-1)+50 LBC(I-1)=LBC(I-1)+50 AIL(I-1)=(HBY(I)-HBX(I-1)/LBC(I-1)AIL(I-1)=(HBY(I)-HBX(I-1)/LBC(I-1)END IF END IF50mHBY(i)HBX(i)新的BPD 另外还有几种情况，大家下去自己推导一另外还有几种情况，大家下去自己推导一下，不再一一陈述。下，不再一一陈述。七、重新计算坡度七、重新计算坡度
42、DO 140,I=1,J-1DO 140,I=1,J-1 AIL(I)=(HBY(I+1)-HBX(I)/LBC(I)AIL(I)=(HBY(I+1)-HBX(I)/LBC(I)140 CONTINUE140 CONTINUE八、超限坡调整八、超限坡调整 DO 150,I=2,JDO 150,I=2,J IF(ABS(AIL(I-1)IF(ABS(AIL(I-1)*1000)+0.5.GT.12.5)THEN1000)+0.5.GT.12.5)THEN HBX(I-1)=HBX(I-1)+LBC(I-1)HBX(I-1)=HBX(I-1)+LBC(I-1)*(12.5-(12.5-(ABS(A
43、IL(I-1)(ABS(AIL(I-1)*1000+0.5)/20001000+0.5)/2000 HBY(I)=HBY(I)-LBC(I-1)HBY(I)=HBY(I)-LBC(I-1)*(12.5-(12.5-(ABS(AIL(I-1)(ABS(AIL(I-1)*1000+0.5)/20001000+0.5)/2000 end if end ifContinueContinue 前点升高，后点降低！方法比较多！前点升高，后点降低！方法比较多！九、再次计算坡度九、再次计算坡度 DO 160,I=1,J-1DO 160,I=1,J-1 AIL(I)=(HBY(I+1)-AIL(I)=(HBY(I+1)-HBX(I)/LBC(I)HBX(I)/LBC(I)*10001000160 CONTINUE160 CONTINUE十、输出结果文件十、输出结果文件 DO 190,I=1,J-1DO 190,I=1,J-1 WRITE(2,300)LBC(I),AIL(I)WRITE(2,300)LBC(I),AIL(I)300 FORMAT(1X,I4,7X,F5.1)300 FORMAT(1X,I4,7X,F5.1)190 CONTINUE 190 CONTINUE
