约束优化设计

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第四章 约束优化设计

● 概述

● 约束坐标轮换法 ● 随机方向法 ● 罚函数法

概述

结构优化设计的问题,大多属于约束优化设计问题,其数学模型为:

根据求解方式的不同,可分为直接解法和间接解法两类。

直接解法是在仅满足不等式约束的可行设计区域内直接求出问题的约束最优解。属于这类方法的有:随机实验法、随机方向搜索法、复合形法、可行方向法等。其基本思路:

在由m 个不等式约束条件g u (x )≤0所确定的可行域φ内,选择一个初始点0

X 然后确定一个可行搜索方向S ,且以适当的步长沿S 方向进行搜索,取得一个目标函数有所改善的可行的新点1

X 即完成了一次迭代。以新点为起始点重复上述搜索过程,每次均按如下的基本迭代格式进行计算:

k+1k k k =+S (k=0,1,2,..)X X α逐步趋向最优解,

直到满足终止准则才停止迭代。

直接解法的原理简单,方法实用,其特点是:

1) 由于整个过程在可行域内进行,因此,迭代计算不论何时终止,都可以获得比初始点好

的设计点。

2) 若目标函数为凸函数,可行域为凸集,则可获得全域最优解,否则,可能存在多个局部

最优解,当选择的初始点不同,而搜索到不同的局部最优解。 3) 要求可行域有界的非空集

1,2,...,1,2,...,u m v p n ==

f x n x R ∈..s t ()0

()0u v g x h x ≤=

间接解法

间接解法是将约束优化问题转化为一系列无约束优化问题来解的一种方法。由于间接解法可以选用已研究比较成熟的无约束优化方法,并且容易处理同时具有不等式约束和等式约束的问题。因而在机械优化设计得到广泛的应用。

间接解法中具有代表性的是惩罚函数法。将约束函数进行特殊的加权处理后,和目标函数结合起来,构成一个新的目标函数,即将原约束优化问题转化为一个或一系列的无约束优化问题。

然后对新目标函数进行无约束极小化计算。

间接法是结构优化设计中广泛使用的有效方法,其特点:

1) 由于无约束优化方法的研究日趋成熟,为间接法提供可靠基础。这类算法的计算效率和

数值计算的稳定性大有提高;

2) 可以有效处理具有等式约束的约束优化问题;

3) 目前存在的主要问题,选取加权因子较为困难,选取不当,不仅影响收敛速度和计算精

度,甚至导致计算失败。

a) 可行域是凸集;b)可行域是非凸集

()

()()()121211

,,m l j k j k X F X G g X H h X φμμμμ==⎡⎤=++⎡

⎤⎣⎦⎣⎦∑∑

新目标函数 加权因子

● 约束坐标轮换法

约束坐标轮换法是在无约束坐标轮换法的基础上,再加上由约束条件构成的可行性逻辑判断而构成的方法,这样可以使搜索点保持在可行域内,求得最优解。迭代步长不是采用最优步长,而是加速步长。其基本思路: 在可行域任取一点0X ,取一个初始步长0α,

按10

11X X e α=+,取得沿1x 坐标轴第一个迭

代点,检查该点是否满足可行性和适用性:

11X D ∈(可行性条件)

101()()F X F X <(适用性条件)

若两者均满足,步长加倍,迭代计算

1

0212X X e α=+,只要迭代点满足条件,加

倍增大步长,继续迭代获得新点;当迭代点不

满足条件,取前一个迭代点,转而沿2x 坐标轴方法搜索,不满足条件时,取负步长进行,如图所示,直到逼近最优点X *。

约束坐标轮换法虽然方法简单、算法明确,便于设计,但当维数较高时收敛速度慢,还会出现“死点”,导致出现伪最优点。

● 约束随机方法

随机方向法的基本思路:

在可行域内选择一个初始点,利用随机数的概率特性,产生若干个随机方向,并从中选择一个能使目标函数值下降最快的随机方向作为搜索方向S 。

从初始点0

X 出发,沿S 方向以一定步长进行搜索,得到新点X ,新点X 应满足约束条件且0

()()f X f X <,至此完成一次迭代。随机方向法程序设计简单,搜索速度快,是解决小型机械优化问题的十分有效的算法。如图所示。

1. 随机数的产生

下面介绍一种常用的产生随机数的数学模型

首先令353637

1232,2,2r r r ===取r=2657863,按一下步骤计算:

令 若 则 若 则 若 则 则 (0,1)之间的随机数 在任意(a,b)区间内的随机数

2. 初始点的选择

随机方向法的初始点0X 必须是一个可行点,既满足全部不等式约束条件。 初始点可以通过随机选择的方法产生。

1)输入设计变量的下限值和上限值,即 2)在区间(0,1)内产生n 个伪随机数i q

3)计算随机点x 的各分量 4)判别随机点x 是否可行,若随机点可行,用x 代替x0为初始点;若非可行点,转到步骤

2)重新产生随机点,只到可行为止。

3. 可行搜索方向的产生

产生可行随机方向的方法:从k 个随机方向中, 选取一个较好的方向。其计算步骤为: 1)在(-1,1)区间内产生伪随机数j

i r , 得随机单位向量j

e

()121

211...j j j n j j i n i r r e r r =⎡⎤⎢⎥⎢⎥=⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥

⎣⎦

⎣⎦∑ 2) 取一试验步长0α,按下式计算k 个随机点

00j j

X X a e =+

3)检验k 个随机点是否为可行点,除去非可行点,计算余下的可行点的目标函数值,比较其大小,选出目标函数最小的点L

X 。

4)比较L X 和0X 两点的目标函数值,若0

()()L

F X F X <,则取L X 和0

X 连线方向为可行

搜索方向;若0()()L F X F X >,则步长0α缩小,转步骤1)重新计算,直至0

()()

L F X F X <5r r ←3r r ≥3r r r ←-2r r ≥2r r r ←-1r r ≥1r r r ←-1

/q r r =()x a q b a =+-i i i a x b ≤≤()

i i i i i x a q b a =+-

相关文档
最新文档