最优化方法探索.docx

上海电力学院最优化方法报告院　　系：　能机学院　　　专业：　动力机械及工程　学生姓名：　张焱 任课教师：　　薛文娟　　　　2015年　6　月　21　日1动力机械及工程简介动力机械及工程是“动力工程及工程热物理”一级学科的重点组成部分，它以工程热物理为主要理论基础，与工程力学、机械工程学、自动控制、计算机、环境科学、微电子技术等学科互相交融，密切相关。动力系统与机械建模、仿真、优化，动力机械与设备的气动热力学，燃烧理论与技术，动力机械工作过程及排放净化，动力机械的控制理论与技术，热力机械的结构分析及设计方法，新型动力机械是该学科研究的范围。动力机械及工程学科其侧重点主要以燃气轮机、汽轮机、内燃机和正在发展中的其它新型动力机械及其系统为对象，研究如何把燃料的化学能和流体动能安全、高效、低污染地转换成动力的基本规律和过程，研究转换过程中的系统和设备，以及与此相关的控制技术，同时，对于风能、太阳能的利用也积极探索，如何能高效地获得这些能源。它涉及能源、交通、电力、航空、农业、环境等与国民经济、社会发展及国防工业密切相关的领域。2问题及解决方案2.1 非线性结构优化问题我所研究的方向是燃气轮机转子非线性动力学，燃气轮机转子在运转的情况下，其工作条件非常恶劣，承受高温、高压。另外燃气轮机转子的结构复杂也是一大难题，复杂的结构必然导致模型参数维数高，所以在对结构进行优化的过程非常困难，所以我们需要对其非线性特性进行优化，然而结构非线性优化的梯度计算代价太高，而且当考虑结构的动态效应时，处理时域内的函数也相当复杂。因此，利用传统的基于梯度的数值优化方法进行结构动态非线性优化是比较困难的。近年来，全局优化算法得到了快速发展，如粒子群优化算法(Particle swarm optimization，PSO)、遗传算法(Genetic algorithm，GA)、自适应模拟退火法(Adaptive simulated annealing，ASA)等。这些方法不需要进行梯度分析，具有较好的全局寻优性，一般支持并行计算。但是这些方法也有各自的不足：PSO 方法搜索精度不高，不适合处理离散优化问题；GA 方法效率较低，可能发生早熟收敛，且需要形成大量的种群数目和多次种群进化，计算资源消耗较大；ASA 方法多次退火计算时间较长，占用大量计算资源，不适合求解设计空间平滑或局部最优解较少的问题。全局优化算法与基于梯度的优化算法相比，消耗的计算资源较大，加之结构动态非线性分析本身的计算代价很高，因此不适合求解大规模结构动态非线性优化问题。等效静载荷法(Equivalent static loads method，ESLM)由KANG等于2001年针对动态线性优化问题而提出，从最初考虑单个时间步等效到考虑所有时间步等效，并从数学上证明了该方法解的最优性。等效静载荷法将非线性优化分为分析域和设计域：分析域进行动态非线性分析得到节点位移、加速度、应变等输出响应，设计域利用节点位移、等效静态载荷等作为输入条件进行静态线性优化，并返回其优化解作为分析域下一次分析计算的初始值。等效静载荷法利用结构线性静态优化技术求解结构动态非线性优化问题，可以看成一种基于半梯度的优化方法。与传统的基于梯度的数值优化算法和全局优化算法相比，该方法可以极大地提高优化效率，这是等效静载荷法的最大优势。当结构的非线性很显著或设计变量较多时，非线性分析与线性分析的差异将增大，相邻两次静态线性分析的设计变量相对变化值减小，优化效率降低。为提高优化效率，一个可行的方法是加大设计变量的迭代步长，但是这种方法可能导致优化不收敛。本文利用结构线性静态优化得到最优解处的梯度信息，结合最速下降法，针对大变形和多变量结构动态非线性优化问题提出基于梯度的等效静载荷法 (Equivalent static loads method based on gradient，ESLMG)。将静态线性结构优化视为对动态非线性结构优化进行的一次梯度分析，利用该梯度信息并结合最速下降法得到设计变量迭代准则，在梯度方向上加大设计变量迭代步长，同时将更新的设计变量作为下次静态线性优化的初始值。ESLMG与ESLM相比，可以显著提高优化效率，同时保证算法的稳健性与结果的可靠性。2.11 结构动态非线性优化数学模型（1）式中，分别为质量矩阵和刚度矩阵，均为设计变量与位移矢量的函数；为加速度矢量；为速度矢量；下标表示响应由非线性分析获得。式(1)是基于有限元法的结构动态非线性分析控制方程。常数表示总时间步；与分别表示约束函数与设计变量的个数为第个时间步的外载；和分别表示第个设计变量的下限与上限。在第次外层循环中，令为非线性分析的位移响应，为内层循环第一次线性分析的位移响应，根据等效静载荷的定义，有（2）定义以替代与，并