《【数学与应用数学】论文——大塘路段开山建路优化方案》由会员分享，可在线阅读，更多相关《【数学与应用数学】论文——大塘路段开山建路优化方案（7页珍藏版）》请在金锄头文库上搜索。
1、54 大塘路段开山建路优化方案大塘路段开山建路优化方案摘要摘要:本文以大塘路的修建为背景, 对局部路段所实施的开山建路工程为主要问题, 并对该路段进行深入的实地研究, 提出了以工程的总造价最少为主要判断依据的最优方案. 在搜寻最优方案的过程中, 本文从边坡的整体稳定安全系数着手来分析边坡的稳定性, 对在给定的了边坡的坡角和圆弧滑动面的曲线方程的情况下，并采取了一定的简化，运用条分法及刚体平衡法来进行边坡局部的稳定性分析，再进行整体分析，求出最危险可能圆弧滑动面得出了使得总造价最少的最优坡角的数学规划模型，运用计算机强大的运算能力进行搜索，应用数学软件 Maple 进行求解，得出了在各假定的坡角
2、下相应的总造价，综合比较，得出总造价最少的最优坡角（=）下相应的总造价，最终求出整个开山建路工程的总费用（=176，4417 元）063.关键词关键词:条分法；边坡稳定性；下滑力矢量；等值面法1问题的提出问题的提出 大塘路连接着大塘与韶关市区，由于这条路路面较差，特别是在下雨时，更是影响到 这条路上来往的行人，而且路间有三间学校（铁中、南方技校、韶关学院） ，为了方便人们 及学生的来往，现政府决定对这条路进行修建，由于路边较多山丘，开路时，要对它进行 评估，对这一条路段的费用进行决策，并提出一种优化的方案，使得费用最少. 大塘路开山建路中的某一路段（大概离火车站 6、7 公里处）长是 L 米，
3、开挖剖面简 化如下图，在山体开挖一底宽（路宽）为 S 米，两侧边坡高为 H（不定）米，其中山体至高点高为 H，坡角为 米（待求）的公路.max山体一挖后，由于重力和地应力释放作用，边坡可能沿某一圆弧滑动.为了保证开挖后 道路两侧边坡的安全稳定性，必须进行加固，以免发生两侧滑坡，对道路上行人生命安全造成危险.对加固有一个要求，即是要保证边坡稳定的安全系数不低于某一值,在此情况sof下确定最佳开挖坡角(当此路段建路费最少时).HHSS以上是山体立体图和某一山体剖面图:2 问题的分析问题的分析以上问题主要是要确定这一个路段两侧的边坡的坡角,在符合某一安全条件下,使得费用 最少的问题.在一般土坡稳定分
4、析的工程实践中，一般采用条分法.对此问题，我们亦可以运 用条分法对这一路段进行分析, 条分法虽然有事先假定滑裂面,不考虑土体实际的应力-应变 状态等不合理之处,但它毕竟抓住了问题的主要方面,而且经过数十年的研究,应用和修正,已 积累了很多经验.对于本文的问题,我们首先需要解决的问题是对某一山体剖面进行研究时,第一期（2002 年 10 月） 韶关学院学生数学建模论文集 No.155 如何确定潜在最危险滑动面,而后就是给出早就这一剖面费用研究最优化模型,然后把它延伸 到 n 个剖面的情况,把这些剖面叠加在一起,即是在山体上开山建路后的立体图,最终给出一 个路段开山建路的总费用最少的优化模型.3
5、模型的假设及符号的约定模型的假设及符号的约定首先, 我们假设山体是两边对称的, 而后对山体所受的力进行分析, 之后得出以下假设:山体地应力沿高度分布,且与山体深度成正比,即山体地应力 ,其中,侧压力系数为 H0,山体的容重为,山体的磨擦系数为 K ,山体的粘聚力为 C.我们必须先对山体以垂直于道路的方向切割山体,得出它的剖面图,之后就对这一剖面 进行单独的研究, 假设滑动面是圆弧, 并假设以 1m 为单位,用条分法把边坡分成 n 等份,忽 略各土条之间的相互作用力的情况下,取任意第 i 个土条 j 为隔离体, 如图(一), 运用刚体平 衡法对它进行受力分析, 假设水平应力(地应力)与垂直应力(
6、自重应力)都随着岩体的深度增 大而线性增大,它们的比值即为侧压力系数 ,并约定以下符号: : 路宽; S: 边坡高; H ijN: 边坡的坡角; ijC ijF:圆弧滑动面圆心坐标,r 00, yx iF0为半径;:土条 j 底部中间的横坐标、纵 yx, ivF坐标的值; 图一: 滑动范围内的土条数;n: 土条 j 圆弧底部的切向方向与水平方向的夹角;: 土条 j 圆弧底部的长度,为简化计算,我们在这里取;jlcos/1:土条 j 所受的垂直应力大小,取值为; ivFiilhl:土条 j 所受的水平应力大小,取值为; iF0jjlhl0:土条 j 在法向方向上所受怕支持力； ijN：土条所受的
7、粘聚力，这里即；jjilC ：土条所受的摩擦力；这里为 ijFjjNk ：土条所受的阻滑力; ijfj: 土条所受的下滑动力; ijsj56 : 土条的稳定安全系数; isfj,:分别为一个单位上的加固、开挖费用;1R2R:分别为加固、开挖及它们的总费用;MMM,214 条分法静力分析离散模型的建立及求解条分法静力分析离散模型的建立及求解传统上沿用的评判边坡稳定性的量化指标是稳定安全系数,它最初是以刚体滑块模型 上的阻滑力与滑动力的比值来定义的, 从可能滑体的抗滑稳定性看, 采用阻滑力与滑动力的比值定义稳定安全系数,并以,是否大于某一个常数(=1.82)为抗滑稳定判据是合理的.该sfsf准则被
8、现行规范所采用,但实施的前提是事先要给出滑动面.在条分法静力分析下给出下列讨 论: I.刚体平衡分析:对第 条土条的圆弧底部的法向方向上达到力的平衡,即ijcossin0iiivjFFN则我们可以求出土条所受的支持力、粘聚力及摩擦力,它们分别为:支持力:cossin0 iiivjFFN粘聚力: jjlCC i摩擦力: iijjNKF根据这一土条水平与坚直方向上力的分析得出:j阻滑力为: iiijjjFCf水平阻滑力为:cos ijf垂直阻滑力为:sin ijf下滑力为:sincos0 iiiFFsvj水平下滑力:cos ijs垂直下滑力:sin ijsII. 具体实现过程: 在最危险滑弧的搜索
9、过程中,滑弧的形状是不断变化 的,在这里是指它的滑裂面的圆心坐标不确定,对于不同的第一期（2002 年 10 月） 韶关学院学生数学建模论文集 No.157 滑弧, 滑弧内边坡所受的下滑力与阻滑力都是不同的,因此,对于不同的圆弧滑动面,相应的安 全系数也是不同的. 1) 首先,我们假定坡角是某一个固定的值,再确定圆弧滑动面的方程,假设圆弧滑动面的坐标(),则滑弧的方程为.00, yx)( ,)()( :2 02 022 02 01yxrryyxxL则对于,第 个土条来说,它的横坐标即为 ，代入上滑弧方程即可得纵坐标值, 它ijxi1Ly的高度就可以求出, ;h cot,cot,)tan( Hx
10、yHHxyxh它的圆弧底部切向角.)arctan(00 yyxx 2) 最危险滑向的确定: 在滑体(边坡)沿圆弧滑动面滑移时,滑面上每一点下滑力、阻滑力 (即安全系数)与滑移方向是都不相同的,但它总可以合成为整体的下滑力矢量,其方向称为整体的合滑向. 使为最小植的整体的合滑向必然是最危险的滑向.sf滑体沿滑动面滑移的水平、垂直合阻滑力各为: nijjiff1cos 1 nijjiff1sin 2则总的合阻滑力为:. 21jjjfff滑体沿滑动面滑移的水平、垂直合下滑力各为: nijjiss1cos 1 nijjiss1sin 2则,总的合下滑力为:21jjjsss3) 稳定安全系数的正确定义:
11、根据摩擦理论, 滑面上每点的阻滑力与该点相对位移方向相 逆,而相对位移方向就为该点的滑动方向. 有些文献认为阻滑力与相应点合剪应力方向相逆, 我们认为这个观点是不妥的, 因每一点的相对位移方向是不同, 故每点的阻滑力方向也是不 同的. 由于稳定安全系数是个标量, 因此在求阻滑力与滑动力的比值时, 应取其在同一个方向分量的比值. 所以按最危险合滑向来求解才能获得最小的稳定安全系数. 这里的稳定sf58 安全系数其实是边坡的整体抗滑安全系数. 4）计算安全系数: 基于摩擦理论和矢量几何概念, 可以导出稳定安全系数公式:.212121212222)sin()cos()sin()cos()()()()
12、(2121 nijnijnijnijjjjjsiiiissffssfff这样,我们就可以求出在固定坡角下, 得出最小的稳定安全系数, 也就意味着唯一确定 了该山体剖面两侧的边坡的最危险可能圆弧滑动面的方程.事先给定一个加圆心坐标的搜索范围(),().rlxxx0rlyyy0其中求解最小稳定安全系数的数学模型规划模型(I):minsf. .tsrlrlyyyxxx00,5) 然后, 我们对在该坡角上的开挖和加固工程量的总造价进行研究, 计算出总造价最少的最优坡角. 当我们算出的最小稳定安全系数小于时, 我们就要对该工程进行加固,并把 0sf该费用计算入总造价里.对山体加固需要增加的力是:;jis
13、fsf 0加固费用:;11)( 0RfsfMjjs开挖费用:;22)cot(RHHSM总费用:.221)cot()(22 0RHHSRfsfMMMjjs最终,我们由前面的一系列推理、论证及演算,可以得出一个求解某一山体剖面的总施工费用最少的数学规划模型():min M.ts206) 至此,我们求解出了具体对于某一山体剖面的最少费用及相应的坡角大小，和需要提供 的加固力.接下来,我们就要对整个山体的开挖、建路及加固的总费用的大小进行计算,得出 一个衡量它的模型出来. 我们知道,对山体进行垂直切割后的切片(切片厚度为 1m)的计算,上面已经给出了详细 的分析及求解,因为各个山体剖面是离散的,如果把
14、它们用积分的方法,似乎不太现实,而且 真正岩土工程的研究上并不会用这种方法计算,而只需把它们全部累加起来,即可以. 但是,我们还必须知道每一个山体剖面的高度,对山体的受力进行分析知道，山体的各 个剖面的两侧边坡的高度并不相同,故必须找出描述山体形状的数学表达式来,经过我们的 实地考察,可以粗略的认为山体是一个上表面凸起的矩形带状馒头形模型,对问题中给出的 山体立体图仔细分析,可看出修好后道路后任一侧边坡的顶部是一个椭圆弧曲线,方程可近第一期（2002 年 10 月） 韶关学院学生数学建模论文集 No.159 似地视为:.)0,2, 0( ; 1 )2()()2()(00 max2 022 0 2yLxyHyy LxxL在这里,我们可以根据边坡顶部曲线方程求出各个山体剖面(共 L 个) 的施工费用,最后 对它们进行累加,就是我们所要求的最少总费用大小,这也确定了最优坡角,则大塘路段开山 建路的最优方案是:在开挖山体修建路面时,要求路两侧的边坡坡角等于我们求出的最优坡角,并要对两侧边坡进行加固,加固量及费用为最优坡角下相应的值的.1M4 4 数据仿真及求解数据仿真及求解我们对大塘路段某一开山建路进行实地考察，并进行研究，得出以下描述路段中此山体的数据
