结构优化和可靠度分析迭代算法的混沌动力学分析和控制.docVIP

结构优化和可靠度分析迭代算法的混沌动力学 分析和控制 杨迪雄 易平 程耿东 大连理工大学工程力学系,工业装备结构分析国家重点实验室,大连,116023 Email: yangdx@dlut.edu.cn; Tel: 0411摘要：绝大多数实际结构和机械系统都是非线性的，非线性体系的迭代求解是当前工程系统分析的一种普遍性方法。然而，对于一些非线性问题，结构优化的牛顿法、满应力法、序列线性规划、序列二次规划法和可靠度分析的一次可靠度方法以及基于可靠度的优化设计等迭代算法并不能得到期望的收敛解（即不动点解），而得到振荡甚至混沌解。本文结合工程力学迭代算法振荡和混沌控制的特点，介绍稳定转换法和松弛牛顿法为代表的两类（全局和局部）混沌反馈控制方法。然后，综合阐述我们在非线性函数最优化的牛顿法、一次可靠度分析方法和基于功能度量法的概率结构优化设计等迭代算法的混沌动力学机理分析和收敛控制方面的研究进展。 关键词：结构优化和可靠度分析，迭代算法，收敛失败，混沌动力学分析，混沌控制 1 引言 绝大多数实际结构和机械系统都是非线性的，非线性体系的迭代求解是当前工程系统分析的一种普遍性方法。由于迭代方法具有易于理解、编程方便的优点，而被大家广泛使用。然而，对于一些非线性问题，迭代方法并不能得到期望的收敛解（即不动点解），而得到周期振荡解甚至貌似随机、杂乱无章的混沌解。目前，计算力学和工程界对迭代方程数值解的振荡、混沌等复杂现象的产生机理缺乏相应的知识，更不用说如何控制解的收敛行为并获得需要的不动点解。幸运的是，混沌动力学理论为我们全面、深刻理解这些复杂现象并控制收敛失败提供了强有力的方法。 结构优化的牛顿法、满应力法、序列线性规划、序列二次规划法[1]和可靠度分析的一次可靠度方法[2]以及基于可靠度的优化设计[3]等迭代算法对于一些非线性问题都可能出现收敛振荡、混沌现象。过去针对这类不收敛问题的原因解释和处理方法是基于几何意义、直观和经验的，比如认为一次可靠度方法(FORM)不收敛的原因在于极限状态函数设计点的曲率较大，相应的解决方法着眼于如何减少和限制设计点的移动幅度，人为设定一些经验性的限制参数，但对某些问题求解仍不收敛[4]。然而，迭代算法实际上构成了一个非线性映射，即：，本质上是一个离散动力系统。因此，应该从混沌动力学理论出发才能全面准确地理解迭代解的周期振荡、混沌等收敛失败问题[4]。对于非线性函数最优化的牛顿法，杨迪雄等[5]从混沌动力学角度阐明了牛顿映射迭代计算可能的复杂动力学现象，并采用混沌控制的自定位法控制非线性函数最优化的牛顿法出现的周期振荡现象，算例结果表明它是有效的和高效率的。针对可靠度分析的不收敛问题，杨迪雄等[4]构造了典型极限状态函数的非线性映射显式表达式，基于混沌动力学探讨了一次可靠度方法在求解极限状态函数可靠指标时产生周期振荡和混沌解的根本原因。而且，利用稳定转换法使可靠指标的振荡解和混沌解稳定化，获得了期望的收敛解[5, 6]。进一步，作者将稳定转换法应用于基于功能度量法的概率结构优化设计的收敛控制，控制效果令人满意[7]。 本文结合工程力学迭代算法振荡和混沌控制的特点，介绍稳定转换法和松弛牛顿法（自定位法是其特例）为代表的两类（全局和局部）混沌反馈控制方法。然后，综合阐述我们在结构优化、可靠度分析和基于可靠度的优化设计迭代算法的混沌动力学机理分析和收敛控制方面的研究进展。 2 离散动力系统的两类混沌反馈控制方法 混沌反馈控制主要是研究如何将混沌吸引子中嵌入的所需要的不稳定不动点或周期轨道进行目标引导和定位，然后对动力系统施加控制使之运行到需要的目标轨道。同时，它也能镇定动力系统周期轨道中嵌入的不稳定不动点，实现系统的振荡和分岔控制[8-10]。一般地，n维非线性映射即离散动力系统的表达式为： ，, , (1) 式中，迭代函数向量f是可微的，X是n×1维状态向量，p是动力系统的控制参数向量。 对离散系统实施混沌反馈控制，系统加上控制项后变成： (2) 不同的反馈混沌控制方法具有不同的控制项表达式。根据表达式(2)是否含有Jacobi矩阵的形式，可将混沌反馈控制方法分成两类：全局混沌反馈控制法和局部混沌反馈控制法。全局混沌反馈控制法无需求系统Jacobi矩阵，收敛速度较慢，但具有全局收敛性。而局部混沌反馈控制法需要求系统Jacobi矩阵，收敛速度较快，只具有局部收敛性。 2.1 全局混沌反馈控制——稳定转换法 Schmelcher和Diakonos(1997)[8]引入一个适当的线性变换改变系统Jaco