冲击载荷下结构拓扑优化设计与动态响应分析.docxVIP

? ? 冲击载荷下结构拓扑优化设计与动态响应分析 ? ? 史峰源，李世强，刘志芳 （1. 太原理工大学 机械与运载工程学院，应用力学研究所, 山西, 太原 030024；2. 太原理工大学 山西省材料强度与结构冲击重点实验室, 山西, 太原 030024） 拓扑优化的目的是运用结构拓扑信息寻求结构材料在给定设计域内的最佳材料分布，从而满足一个或多个设计约束的同时最大化或最小化某些设计目标[1]. 拓扑优化是研究结构最优拓扑形式的优化问题，在优化的过程中去除多余或者影响小的材料从而构建满足约束条件的最优结构. 目前，比较成熟的连续体拓扑优化方法主要有均匀化方法[2?3]（homogenization method）、变密度法[4?5]（solid isotropic material penalization，SIMP）、水平集方法[6?7]（level set method）以及渐进结构优化方法[7?8]（evolutionary structural optimization，ESO）等. 渐进结构优化方法作为拓扑优化的一种重要方法，最早是由XIE 等[8]基于结构通过逐渐删除低应力材料而使结构进化到最优的概念的基础上提出的. 渐进结构优化方法不会产生中间密度，优化过程中会得到清晰的拓扑构型，但渐进结构优化方法只能不断删除低应力单元，QUERIN等[9?10]提出了双向渐进结构拓扑优化方法（BESO），不但可以逐渐删除低应力单元，同时也可以逐渐添加已经被删除的单元. 这样既能删减单元同时又可以添加单元，解决了最初只能删减单元的问题.HUANG 等[11]提出了解的收敛性与网格无关的一种BESO 方法，这种方法结合了灵敏度滤波和利用历史信息稳定优化方案，使优化过程更加高效平稳，在结构设计问题中表现出了良好的性能. 目前，双向渐进结构拓扑优化方法已经显示出高效和稳健的优化性能，成为一种无论是在学术研究还是工程应用中都广泛使用的设计方法[12]. 从SIGMUND[13]开始将杨氏模量、导电率或导热率不同的两种材料同时引入设计域，在某些情况下用以显著提升结构特性，从此多相材料拓扑优化设计得到广大学者的关注[14]，复合材料由于多种材料特性之间的相互补充，从而在性能方面要比传统材料更加优良，具有比强度高、比模量大、可设计性等优点，并且能满足结构的一些特殊功能，如传热性[15]、导电性等. 目前拓扑优化方法多用于结构静态载荷优化，而工程实际中结构常受到动态载荷[16]，直接对结构进行非线性动力响应优化分析虽然在理论上可行，但是计算过程异常复杂、计算效率低，而且计算难以收敛，无法得到较为理想的拓扑构型[17]. 为了可以高效地进行结构非线性动态优化设计，PARK 等[18]提出了等效静力载荷优化方法，通过构造等效静力载荷的方法将动力优化设计问题转化成为静力优化设计问题[19]，并且验证了方法的可行性，提高了非线性动力优化的设计效率，这种方法也因此得到了广泛应用[20?21]. 本文针对结构非线性动态响应拓扑优化优化流程复杂，计算效率低等问题，将双向渐进结构拓扑优化方法和等效静力载荷优化方法相结合，对非线性动态响应拓扑优化流程进行改进，采用内外双层循环的方式，外层循环实现对结构的非线性动态响应分析，依据结构的非线性动态响应构造等效静态载荷，内层循环依据构造的等效静态载荷，并将权值法引入其中，采用改进的双向渐进结构优化方法进行优化，简化优化流程，优化过程更加平稳，整体优化时间大幅减小，计算效率大幅提高. 采用双向渐进结构优化方法能够避免中间密度的产生，提高了工程可制备性，同时也将材料插值模型引入其中，实现冲击载荷下双相复合材料拓扑优化设计，在保证结构质量相同的情况下使结构具有更高的刚度. 最后通过两个算例来验证改进方法的可行性和高效性，在得到优化结果后对最终优化结构进行对比和动态响应分析，说明本文提出的方法更加适用于动载荷下结构拓扑优化设计. 1 等效静力载荷优化法 如果各个时刻下的结构响应通过积分计算的方法来获取，这会使得结构分析变得十分复杂并且计算效率低. 因此等效静力载荷优化方法提出了一种通过构造等效静力载荷可以高效的获得非线性结构响应的等效分析方法，这种方法要求线弹性静力分析在关心时刻下要与非线性动力分析具有相同的结构响应. 图1 为等效静力载荷的等效过程，在相对应的时间步中，等效静力载荷是能够与结构施加动载荷产生相同位移场的静载荷. 在等效时，将结构的非线性动力学分析按所需的时间步分开，每一个时间步的动态响应等效为静力载荷的一个工况，并且要求相对应时间步结构在动载荷下产生的位移和相对应等效静态载荷产生的位移相等[19]. 图1 等效静态载荷的等效过程Fig. 1 Equivalent process of equivale