频率响应位移幅值敏度分析的伴随法.pdfVIP

维普资讯 第 25卷 第 2期 应 用 力 学 学 报 Vo1．25 No．2 2008年 6月 CHINESEJoURNALOFAPPLIEDMECHANICS Jun．2008 文章编号：1000—4939(2008)02～0247-06 频率响应位移幅值敏度分析的伴随法 彭细荣 隋允康 (清华大学深圳研究生院 518055 深圳) (北京工业大学 100022 北京) 摘要：频率响应位移幅值的敏度分析通常采用直接法，一次敏度分析只能计算出对一个设计变量 的偏导数，这在设计变量很多的拓扑优化中因敏度分析计算量太大而显得不适用，本文推导了频率 响应位移幅值敏度分析的伴随法，一次敏度分析可计算出对所有设计变量的偏导数，算例表明伴随 法计算结果与直接法及差分法结果符合得很好，用伴随法分析敏度在结构拓扑优化中可以大幅提 高计算效率。 关键词：频率响应；敏度分析；伴随法；拓扑优化 中图分类号：0343．1 文献标识码： A 式，对其求导数归结为对特征值及特征向量求导数， 1 引 言 一 般需要计算前几对特征值及特征向量的导数，计 算的模态数越多其计算精度也越高，但计算量更大。 随着对结构动态性能要求的提高，结构动力优 双模态叠加法对运动方程直接求导数，形成关于敏 化问题越来越受到重视。现有的动力优化问题主要 度的一个形式与运动方程相同的方程，右端项设定 考虑结构固有频率的优化 【1]，使其固有频率远离激 为虚载荷，对响应及其敏度都用模态叠加法求解，一 振频率，以免发生共振。但工程中有时需要对结构 般要分三步计算，先用模态叠加法计算响应，再用响 某些位置的位移响应幅值进行限制。结构在简谐扰 应及运动方程中矩阵的导数计算虚载荷，最后用模 力作用下，控制其位移幅值就是很重要的结构动力 态叠加法计算响应敏度。双模态加速法是基于模态 优化领域，但 目前已有的工作主要集中在截面及形 加速法的敏度计算方法，其做法与双模态叠加法类 状优化层次，拓扑优化方面的研究相对较少 2【]。 似，不同之处在于其运动方程的求解是用模态加速 在过去近30年内，人们对结构的特征值及特征 法进行的。以上各种敏度计算方法都是一次敏度分 向量的敏度计算进行了详细的研究并取得了非常丰 析只能计算出状态变量对一个设计变量的导数，这 富的成果3【]，然而对频率响应敏度分析进行的研究 些方法只适用于设计变量数不多的情况，对结构拓 相对较少。对结构系统进行频率响应分析通常有两 扑优化中设计变量数 目成千上万的情况，用上述各 种方法：直接法和模态叠加法，与频率响应分析对应 种方法计算敏度所产生的计算量都是无法接受的。 的敏度分析方法也主要有直接法 4【]、模态叠加法 5【】、 本文对无阻尼系统及阻尼系统频率响应的敏度分析 双模态叠加法[62及双模态加速法[73等。直接法是基 推导了伴随法公式，采用伴随法，单位虚载荷工况可 于频率响应分析的直接法，其不足之处也在于计算 以与原载荷工况同时分析，一次敏度分析可计算出 量太大。模态叠加法将频率响应表达为模态叠加形 一 个状态变量对所有设计变量的导数，在此方法基 * 基金项目：北京市自然科学基金(3042002)；北京市教委基金(KM200410005019) 来稿日期：2006—10-20 修回日期：2oO7-o59 第一作者简介：彭细荣，男，1972年生，清华大学深圳研究生院，博士后；研究方向——