结构优化设计基本方法
作者:小易发布日期:2020-02-17

教学规划法


  数学规划法的命题是:求n个变量xi(i=l,2,…,n),满足m个约束条件Gj(xi)≤0 (j=l,2,…,m),且使目标函数W(xi)为最小(或最大)。如果约束条件和目标函数都是xi的线性函数,这便是线性规划问题,已有成熟的解法;如果在这些函数中有一个是非线性函数,便成为非线性规划问题。随着非线性函数的性质和形式的不同,非线性规划问题有很多类型,特殊的解法很多,在应用上各有局限性,没有普遍适用的最好解法。


  用数学规划法来作结构优化设计,变量xi便代表可以变化的各种结构参数,如元件截面积或厚度、节点位置、材料性质等;约束条件Gj(xi)≤0代表设计必须满足的各种限制,例如结构各部位的静应力,动应力或变位不得超过规定的容许值,元件的截面或厚度尺寸不得超出给定的范围,结构的频率不应落在某个禁区,结构的失稳临界力或飞行器的颤振速度不得小于某一下限,等等;而目标函数则代表结构优化所追求的指标,例如,结构重量最小和成本最低等可以定量的指标;也可将重量、造价作为约束条件,而把某种结构性能,例如刚度作为目标函数。


  数学规划法的基本目的是,在以设计变量为坐标的多维空间里搜索最优点。如果有n个设计变量,则相应的n维设计变量空间中的每个点都代表一个设计方案。在无限多的点中要尽快地搜索出既满足所有的约束条件,又能使目标函数尽量接近最小值(或最大值)的点,就是数学规划设计法的任务,这种搜索的过程称为“优化过程”。


  附图表示一个二维设计空间,图中的一簇曲线是目标函数W(x1,x2)为常数的等值线。约束函数Gj(x1,x2)为零的曲线所围成的区域是可行域。A、B、C点各代表一个可行的方案.围线以外的点(如D)不满足约束条件,所以是不可行方案。显然,满足约束条件并使目标函数W最小的最优方案点是M。数学规划就是要以最迅速的方式找到点M。


  这好比在山坡上—个用栅栏围起来的区域里找最低点,如果这个山坡不是凹的,则可以断定最低点必在栅栏所在的边界上。数学规划提供了很多搜索的办法,基本原则都是在选好一个出发点后,经过分析判断,找出一个迈步的有利方向,沿这个方向跨出有利的步长以到达新的一点。再从此点出发,重复上述过程,一步一步走下去,直到再也找不到可走的有利方向,就是达到了最低点。从第n点到第(n+1)点这一步可表达为:


  式中  为有利方向,为有利步长系数,它们依靠在点进行的分析所提供的信息来确定。例如,从可行点A出发,沿着等高线的梯度负向,即最陡下降方向逐步走到边界点1,然后再沿着边界逐步走到最低点M,这个方法叫作梯度投影法。实际上还有很多其他的方法。可以看出,如果初始出发点选的是B,用同样的走法也可以走到最低点M;但如果初始点选的是C,那就会走到另一个局部最低点N。M点代表全局最优解,因为它是全部可行域中的最低点。


  N点只是在它附近的可行域中的最低点,所以是局部最优解。现在还没有一个可靠的实用方法能保证搜索到的解一定是全局最优解。一般是在可能的情况下取若干不同的出发点作几次搜索,以期找到全局最优解。


  如果是线性规划问题,搜索过程就简便得多。所以有时把非线性问题转化成一系列线性问题来逼近。为此,在某一设计点附近将目标函数和约束函数都线性化,也就是在该点将函数作泰勒展开,并只保留它们的线性项。然后作有一定步长限制的线性规划,得到新的一点。如此重复下去,直到收敛于最优点为止。


  由于不带约束的规划问题比较容易作,所以有时也把有约束问题转化成一个序列的无约束问题。为此,可以把约束表示成一个罚函数加到目标函数上去,构成一个新的目标函数,即式中即为罚函数,r是个相当小的正数,它在序列无约束问题中,逐次减小。因为r值很小,当代表某一设计方案的点在离开边界较远的可行域内部行动时,;但是当接近可行域的边界.某约束函数Gj(xi)将由负值趋近于零,于是罚函数急剧增大,因此,的最小点不可能越过可行域边界。


  r越小,无约束问题的W最小点越接近于有约束问题的W最小点。但是如果一开始就取很小的r,无约束问题将遇到收敛上的困难,所以有必要将有约束问题化成一个序列的无约束问题,让系数r在这个序列中逐渐减小到适当的程度。


  此外,还有一些非线性规划的特殊方法,如几何规划和动态规划,各有其适应的范围,在结构优化设计中也得到应用。


  优化准则法


  以满足某种准则来代替目标函数在约束条件下取极值的方法,叫作优化准则法。最简单的一个优化准则法,便是前面提到的满应力设计方法。只有对于内力分布不随设计变量改变而变化的静定结构,而且容许应力与设计变量无关的情况下,才能通过一次结构分析和修改设计得出满应力结构。对于其他情况,为使各元件趋向于满应力,必须进行下列的选代:


  式中  和  为第n次迭代的第i元件的截面积和最大应力,为第i元件的容许应力。公式给出经过修正的第i元件的截面积。迭代收敛时,就达到的满应力准则。满应力准则和结构最小重量之间没有必然的联系,但是一般的满应力设计可能相当接近于甚至就等于最轻设计。当然,这个方法只适用于受应力约束的最轻设计问题。