非线性规划法优化雨水管道设计.doc

?　314? 2009年　第2期第35卷　总第148期 非线性规划法优化雨水管道设计 王 珂,赵 鑫 2 (11中国建筑西南设计研究院有限公司,四川成都　610081;21大港油田公司滩海开发公司,天津　300280) 　　摘　要:尽管给排水管网系统设计计算中存在着关系 错综复杂的约束条件,但是只要对其中的某些条件适当取舍,合理地应用数学工具,就可以把它简化、抽象为容易解决的数学模型,通过计算得出最优解。非线性规划法基于求导原则,即目标函数的导数为零的点就是所求的最优解。本文拟采用一种新的精确罚函数,使得含约束的非线性规划问题能采用无约束优化方法中许多有效的解析方法,对雨水管道进行优化设计。　　关键词:雨水管道;优化;;　　中图分类号:TU99文章编号:02-02 的解x3满足S(x3)=0,且对任意满足S(x)=0的点x都有f(x)≥f(x3)。 对于问题(112),: (113)(x,μ)(x-S(x) 易知x,μ)。 : (114)min()x∈Rn 33 μ)x,有0≤F(xμ,则对任意的x、μ,μ F(x,μ),且 333 )=0的充要条件是f(xF(x且S(xμ,μμ)≤μμ)= 0。 定理1　如果x3是问题(112)的最优解且μ=f(x3),则x3是罚问题(114)的最优解且F(x3,μ)=0。 3 定理2　设x3是问题(112)的最优解,对某个μ,xμ是 3 罚问题(114)的最优解和问题(112)的可行解,并且F(xμ,μ)≠0。如果μ≤f(x3),那么xμ3是问题(112)的最优解。 1　在一定约束条件下求目标函数的最优值的问题,通常 称为规划问题,通常写为: 　　　　mf(x),x∈Rn g(x)≤0,i=1,2,…,m 　　　　s,i　　(111) hj(x)=0,j=1,2,…,l 当目标函数和约束函数中至少有一个是x的非线性函数时,称式(111)为非线性规划。 由于雨水管道的约束条件是不等式约束,因此,本文只介绍不等式约束问题(112)的最优性条件。 (112)　　　　minf(x),x∈Rn 　　　　s1t1gi(x)≤0,i=1,2,…,m 记D1={x∣gi(x)≤0,i=1,2,…,m},称为问题(112)的可行域。 定义　设x3∈D1,若对某个i有gi(x3)=0,称gi(x)≤0为点x3处的起作用约束,也称有效约束、积极约束。记33 I(x)={i∣gi(x)=0}为起作用约束下标集,简记为I。若 33 gi(x)<0,则称gi(x)≤0为在点x处的不起作用约束。 罚函数法又称为序列无约束极小化方法,是解决非线性约束最优化问题的一个重要的途径,其基本思想是:利用问题中约束函数作出适当的惩罚函数,由此构造出带参数的辅助函数,将约束问题转化为求解一系列无约束的最优化问题。但是,罚函数法的主要缺点是目标函数的病态性质,这种病态性质是由罚因子的无限增大或减小而引起的,从而使得辅助函数的无约束非线性规划的求解出现数值困难。为了改造罚函数法提出的各种方法中,精确罚函数法是一种重要的算法,它吸引了许多人的注意。这种罚函数的构造都是按约束给出的。本文拟采用一种新的精确罚函数,使得含约束的非线性规划问题能采用无约束优化方法中许多有效的解析方法。这种新的精确罚函数不同于已经研究的罚函数形式,在一定条件下同时具有精确性和光滑性。 对于问题(112),其中可行解的集合是D1={x∣gi(x)≤0,i=1,2,…,m}。 令Φ(t)=[max(0,t)]2,于是对于t有连续的一阶导 (t)=2t;当t<0时,Φ′(t)=0。数。当t≥0时,Φ′ 再令s(x)= Φ[g(x)],易知S(x)∑ i i=1m 2　实例分析 某飞机维修机库,屋面面积为5000m2,屋面为拱形, 由中线向两侧倾斜。设计为压力流(虹吸式)屋面雨水排水系统,管材为塑料管。 该机库位于南方某市,取设计重现期=5a,查设计手册q=5138L/s?100m2。机库雨水总量Q=F×q=50×5138=269L/s。选用DN75压力流(虹吸式)雨水斗,取单斗排水量q=12L/s,则需用雨水斗的个数N=Q/q=269/12=2214只。 设计使用雨水斗24只,当取雨水斗的排水量12L/s时,设计重现期略大于5a。每侧12只分成两个系统,每个系统6只雨水斗。斗前水深0107m。雨水系统平面布置见图1。 在已布置管线下,各管段的长度是确定的,如图1所示。根据雨水管道的总投资最小为优化目标建立的目标函数可知,在已定线的情况下,各管段的管径较小时,造价较低。但是,各管段还受一些约束条件的限制。以管段13～14为例,优化计算管道的管径,寻求造价最低: 图1　压力流(虹吸式)屋面雨水排水系统平面布置 ≥0,且S(x)=0 的充要条件是gi(