12 土石坝断面最优化设计土石坝是坝工中重要的坝型之一。

它具有就地取材、施工方便、工期短、造价低、节约大量水泥、钢材，能适应较差的地质条件，安全性能好等优点，在国内外被广泛采用。

目前，土石坝的设计大多采用传统的设计方法，其设计仍处在经验阶段，设计理论也在进一步成熟过程中。

探讨土石坝的优化设计对提高土石坝的设计效率、优化结构布局及料区分布，充分发挥坝料的作用，降低工程造价具有重要的实际意义。

开展土石坝工程的优化设计研究是提高土石坝工程设计水平的一个重要发展方向。

土石坝断面最优化设计是在结构类型、材料布局已定的条件下，根据设计者的某一目标(如用料最省或造价最低等)，并在满足规范要求的前提下，利用数学手段从众多可行方案中通过自动寻优技术，获得结构断面的最优尺寸。

本章主要介绍混凝土面板堆石坝与土质心墙堆石坝的断面优化设计数学模型，并给出几个工程算例。

12.1 岩基上混凝土面板堆石坝断面优化设计12.1.1 岩基上混凝土面板堆石坝断面优化设计数学模型12.1.1.1 设计变量图12.1给出了岩基上面板堆石坝典型断面及设计变量示意图。

对于岩基上面板堆石坝断面，坝高、坝顶宽度都是根据工程规划要求确定的不变参数，故一般选取描述土石坝断面形状的其他一些几何特征量为设计变量。

(1)取上下游坝坡角x1、x4为设计变量，以反映整个断面的大小。

(2)取堆石料分界坡角x2、x3为设计变量，使各种石料得到合理利用。

(3)取下游坝坡变坡角x5、x6以及相应高度x7、x8为设计变量，以反映下游坝坡变坡要求与设置马道的需要。

一般情况下，岩基上面板堆石坝典型断面的优化设计变量可表示为X=[x1 x2 (x)8]T自然，针对不同的具体工程，设计变量的选取亦有所不同。

12.1.1.2 目标函数岩基上面板堆石坝断面有不同材料分区组成，对于特定的工程而言，综合考虑开采、运输及施工等因素可以确定所用不同材料的方量单价。

这样对单位长断面，其目标函数可表示为F(X)=∑断面内某种坝料方量单价×断面上该料区的面积如果将各料区的单价进行比较，选取某材料的单价为1.O，则其他材料单价与它的比值，即为单价比(记为Ci)。

从而目标函数可进一步表示为F(X)=∑某种坝料方量单价比Ci ×断面上该料区的面积Si12.1.1.3 约束条件岩基上面板堆石坝断面优化设计约束条件有几何约束和性态约束。

(1)几何约束条件。

几何约束条件是设计变量应遵循的取值范围，亦即对设计变量几何尺寸限制。

土石坝优化设计中，一般可根据国内外已建工程的经验拟定。

面板坝上下游边坡 1:1.8≤tanxi≤1:1.3保证分区不会重叠 x4≤x3≤x2≤x1(2)性态约束。

性态约束是保证岩基上面板堆石坝在各种工况下正常工作，安全运行所要满足的稳定、应力、变形等的限制条件。

根据岩基上面板堆石坝的结构特点和工作条件，面板坝的水压力作用于上游坝坡，由坝体自重、面板上水重所产生的抗滑力远大于水的水平推力，一般不存在整体滑动问题，故一般不做整体滑动验算。

面板坝由于其防渗体--混凝土面板在上游表面，坝体堆石为自由排水体，故坝内不存在水的渗透压力及孔隙水压力问题，当然也不存在地震时产生附加孔隙压力问题。

对岩基上面板堆石坝而言，一般应考虑以下的条件。

1)上下游边坡稳定条件FS (X)≥[FS]各工况下边坡稳定安全系数不小于规范规定的允许稳定安全系数。

2)最大应力条件σ1(X)≤[σ1]堆石体、混凝土面板的应力分别不超过允许值。

3)最大变位条件δ(X)≤[δ]坝体最大沉降或面板最大挠度，或周边接缝张开及错动变位分别小于允许值。

4)防止塑性剪切破坏条件SL<1.0堆石体应力水平不超过1.0。

12.1.1.4 优化数学模型综合上述分析，针对具体工程，取定一些约束条件，构成了相应工程的断面优化设计数学模型。

12.1.2 算例12.1.2.1 工程概况某水电站枢纽主要任务是发电，兼顾灌溉供水。

正常蓄水位为2005.0m，校核库水位为2008.0m。

混凝土面板堆石坝方案最大坝高为129.0m。

总库容量为5.5亿m3，最大发电水头为107.4m，装机容量为50万kW，电站保证出力48.6万kW，年发电量为49.5亿kW·h，工程规模属一等大(1)型工程。

典型断面如图12.2所示。

(1)设计变量。

设计变量取上游坝坡角为x1，主堆石与任意料区分界坡角为x 2，任意料区下游堆石分界坡角为x3，下游坝坡角为x4。

这样，设计变量可表示为X=[ x1 x2x3x4]T(2)目标函数。

各料综合比价Ci取为：主堆石料1.00，垫层料、过渡料1.20，砂砾料O.90，强分化料O.70。

由断面的几何参数确定相应的面积si，即可得到目标函数表示式。

(3)几何约束条件2.485≤x1≤2.544(rad)1.500≤x2≤2.300(rad)O.785≤x3≤1.200(rad)O.500≤x4≤O.656(rad)(4)该工程主要考虑的性态约束。

1)上下游坝坡静力稳定安全系数FS≥1.30。

2)堆石体最大主应力σ1(X)≤[σ1]=2γH(按国外工程经验，H为坝高，m)。

3)堆石体最大坝体最大沉降δ(X)≤[δ]=H／150(m)(按国外工程经验，H为坝高，m)。

4)面板最大挠度δ(X)≤[δ]=H／200(m)(按国外工程经验，H为坝高，m)。

5)堆石体应力水平SL<I.O。

12.1.2.2 最优设计方案的求解计算工况，上游水位为2005.0m，下游枯水。

结构分析时土石料本构关系采用邓肯E—μ模型。

分级加荷模拟施工过程，面板与垫层之间均设置Goodman单元，以反映相互间的变形。

坝坡稳定分析用瑞典圆弧法。

计算参数见表12.1。

面板、趾板的弹性模量取25.0GPa，泊松比取O.17，容重取24.0kN／m3。

接触面参数K1、n、Rf、δ依次取4800、O.56、O.74、O.64(rad)。

初始设计方案X0 = [2.520 2.200 1.090 0.530]T (rad)初始设计的目标函数值为F(X0)=22644，经过罚函数优化过程的15次迭代得到优化解X* = [2.491 2.293 1.191 0.573]T (rad)= [142.706 131.374 68.228 32.821]T (°) 优化解的目标函数值为F*(X)=20712，比初始设计降低了8.53％，蓄水期下游坝坡的稳定安全系数为1.336，坝体最大沉降为O.777m，为坝高的O.6％，坝体最大主应力为2027.0kN／m2，坝体最大应力水平为0.656。

蓄水后面板最大挠度为30.5cm，优化方案与原设计方案的断面比较见图12.3。

由图12.3可见，优化后上下游坝坡稍增大，价格比较低的强风化料区在断面中的比重增大。

从优化设计坝体的应力位移及稳定安全系数看，所得方案的坝型是合理的，是在造价最低这个特定目标下的最佳方案。

12.2 覆盖层地基上混凝土面板堆石坝断面优化设计近年来，我国在建和拟建的在一定深度覆盖层地基上的面板堆石坝增多，在覆盖层地基上建混凝土面板堆石坝，一般采用混凝土防渗墙处理坝基渗流，将趾板建在防渗墙顶部，坝体直接建在覆盖层地基上。

覆盖层地基上面板堆石坝与岩基上面板堆石坝的优化设计有相同的地方，亦有其特殊之处。

表12.1 坝体堆石料参数坝 料γ(kN ·m -3) △φ (°) φ(°)k n R t G F D K ur n ur垫层料 22.0 12.0 55 1000 0.32 0.60 0.45 O.15 0.50 2100 0.32 堆石料 21.5 12.0 54 1000 O.32 O.66 0.46 O.16 5.20 2100 0.32 砂砾料 22.0 10.7 53 900 O.50 0.77 O.50 O.26 13.0 1890 O.50 强分化料21.O11.350450O.30 0.650.22 0.066.189500.30图12.2 某水电站混凝土面板堆石坝典型断面及设计图12.3 优化方案与原设计断面比较示意图12.4 土石坝广义模糊优化设计土石坝作为水利工程中的一种重要的水工建筑物，存在着大量随机性和模糊性不确定性因素，如设计的优劣标准、荷载、设计参数、计算模型、材料强度及物理量的允许区间范围，结构的刚度、稳定性、频率等因素，它们既具有随机性，又具有模糊性，而这些因素又直接影响着结构的安全性和经济性。

在进行土石坝结构设计中通常考虑的不确定性因素主要包括：大坝上下游水位的随机性；地震烈度及场地等级划分的不确定性；岩土地基或筑坝材料参数的随机性；土石坝稳定分析和应力应变分析的计算模型的不确定性；土石坝结构设计中的容许应力、允许位移、尺寸限制、频率禁区的模糊性；土石坝的坝料分区、刚度、稳定性等要求都存在模糊性。

本章前三节描述的土石坝断面优化设计，没有考虑土石坝结构设计中的众多物理量的边界具有的中间过渡性，称为普通优化设计或确定性优化设计。

如果仅考虑约束允许范围的模糊性，从而得到土石坝模糊优化设计模型，称之为土石坝的普通模糊约束优化设计。

本节中将要讨论的是土石坝的广义模糊约束优化设计，即同时考虑土石坝优化模型中的约束允许范围、约束函数和目标函数的模糊性。

12.4.1 土石坝广义模糊优化设计数学模型根据土石坝断面确定性优化模型，建立土石坝广义模糊约束优化模型。

⎪⎭⎪⎬⎫=⊂=M m G X g t s X F x x x m m Tn ,...,2,1)~(~..)(min ]....[X 21求 式中目标函数仍然采用工程造价最低为目标函数，同时考虑价格比系数及坝料分界线的模糊性，并考虑设计变量尺寸限制边界的模糊性。

约束函数g m (x)和约束允许区间G m 包括模糊性和非模糊性因素，如考虑土石坝稳定的模糊随机可靠指标为模糊约束函数，土石坝坝坡稳定的允许可靠指标的模糊性，土石坝坝体或坝基允许渗透坡降J 的模糊性以及土石坝坝料价格比系数的模糊性等。

当然，普通约束函数和普通约束允许区间可以看作模糊约束函数和模糊约束条件特殊型式。

在实际应用过程中，根据具体的坝型(心墙坝或面板坝)和不同的分析方法(静力分析或动力分析)选取不同的约束条件。

上式中不仅目标函数具有模糊性，而且约束函数的物理量及其取值允许范围具有模糊性，故称为广义模糊约束优化设计问题。

若完全不考虑各个量的模糊性，则广义模糊约束优化模型退化为常规的优化设计模型。

12.4.2 隶属函数的选取隶属度表示在模糊集合中每一个元素隶属于模糊集合的程度，可以通过隶属函数来表示。

在土石坝的广义模糊优化问题中，统一采用斜线型隶属函数进行讨论。

各模糊量的隶属函数的斜线表达式主要有几种形式：降半梯形分布、升半梯形分布和梯形分布。