非线性结构周期解共振峰值的优化方法.docxVIP

非线性结构周期解共振峰值的优化方法 非线性模态性质的求解方法 工程中非线性结构的循环解算方法大致可分为两类：时间间隔方法（如射击法）和频率范围方法（如波形平衡法及其变体）。 除了周期解的求解方法研究外, 周期解的稳定性分析是很重要的 当非线性结构某个参数连续变化时, 连续延拓方法通常用于跟踪周期解。最近, 文献[4]采用打靶法和伪弧长连续方法研究非线性系统的非线性模态性质。在渐近方法 (Asymptotic Numerical Meth-od) 框架内, 文献[5]提出了一种结合谐波平衡法和HILL法的连续方法。在文献[6]中, 谐波平衡法和伪弧长方法被用于分析几何非线性叶盘结构的自由和强迫振动特性。 研究确定非线性结构共振峰值的方法是很有必要的。例如, Petrov应用谐波平衡法计算包括摩擦阻尼影响的失谐叶盘结构最坏振动情形 为了克服重复求根计算和处理参数不确定问题, 本文提出了一种非线性结构共振峰值求解方法。下面首先介绍确定非线性结构共振峰值的新方法, 然后通过典型Duffing振子算例验证本文方法并通过几何非线性叶盘结构数值算例演示本文方法的优点, 最后给出相关结论。 1 非线性系统周期说 本节提出了确定非线性结构共振极值的方法, 将确定非线性结构共振峰值问题转换为非线性约束优化问题, 首次采用非线性代数方程组等式约束和稳定性不等式约束计算非线性系统的周期解。 下面首先研究基于时域打靶法的非线性等式约束, 其次分析基于状态转移矩阵稳定性分析方法的非线性不等式约束条件, 然后综合非线性等式约束和不等式约束限制条件, 给出基于时域打靶法和状态转移矩阵稳定性分析方法的共振峰值求解方法, 最后采用OQNLP多重启全局优化算法求解该非线性约束优化问题。 1.1 刚度运动方程 采用打靶法求解非线性系统的周期解, 考虑具有n个自由度机械系统的运动方程 式中M , C和K分别表示质量、阻尼和刚度矩阵; 引入状态向量 打靶法的实质是求解边界值问题, 该边界值问题通过如下打靶函数定义 式中 打靶法的关键是寻找满足式 (3) 所示打靶函数的初始条件z 1.2 基于状态转移矩阵特征乘子求解 采用状态转移矩阵法来判定周期解的稳定性。定义D为导算子, 则状态转移矩阵 式中 通过对式 (4) 在一个周期内数值积分便得到状态转移矩阵。最终, 便可计算得到状态转移矩阵在周期T处的N=2n个Floquet特征乘子:ρ 其中ρ=[ρ 采用Floquet理论, 通过解式 (4) 所示的初始值问题, 得到状态转移矩阵的特征乘子。最终, 式 (5) 表示的周期解稳定性条件便构成了非线性约束优化问题的非线性不等式约束条件。 1.3 使用oqnlp多时间节点计算 本文目标是求解非线性结构的共振峰值, 所以式 (3) 表示的非线性代数方程组和式 (5) 表示的周期解稳定性条件必须联立。因而, 寻求非线性结构中具有最大振动幅值的周期解可转化为下述非线性约束优化问题: 式中 ‖u‖ 非线性约束优化问题式 (6) 的求解和有效计算是很重要的问题。在本文研究中, 采用文献[8]的OQNLP多重启算法求解式 (6) 。 OQNLP多重启算法优化过程分为两个主要阶段。在第一阶段计算所有试点的罚函数, 从OptQuest数据库选择具有最好罚函数值的试点作为SQP算法的起始点。在主要迭代循环第二阶段, 选取满足位移过滤和绩效过滤条件的试点启动局部搜索算法进行优化求解。OQNLP多重启算法的具体原理论述可参考文献[8]。 2 计算值的示例 为验证本文方法, 并演示其能力, 本节给出2个数值算例。 2.1 duffing振子 采用Duffing振子作为算例。为验证本文方法并分析其精度, 本文方法结果将与文献[6]中HBM-ANM-HILL方法得到的结果进行比较。Duffing振子的运动方程为 式中 μ, β表示阻尼系数和非线性刚度系数;f表示力幅值。 给定f=1.25, μ=0.1 和β=1, 使用HBM-ANM-HILL方法得到的Duffing振子的频率响应曲线如图1所示, 其中H 从图1中可观察到典型的骨架曲线。图1中在激励频率ω=2.44时, 系统达到共振峰值2.624, 并存在2个分叉点S1和S2。在激励频率ω=2存在多解C, D和E, 其中D是不稳定的周期解。下面研究三种情形以验证本文方法: a找到共振峰值p 为搜寻共振峰值P, 式 (6) 中优化目标设置为具有稳定周期解的Duffing振子振动幅值最大化。需要确定的未知优化变量为初始条件z 稳定性正则方程 对于求解分叉点S1和S2, 根据Floquet理论, 分叉点的所有Floquet乘子的模的最大值等于1, 考虑到数值精度, 式 (6) 稳定性不等式约束条件更改为|max (|ρ|) -1|10 仿真结果及分析 在计算多解集