基于无网格局部PetrovGalerkin法板结构拓扑优化.docVIP

基于无网格局部PetrovGalerkin法板结构拓扑优化 　　摘 要：为将无网格法的优势集成到结构拓扑优化中，基于无网格局部Petrov-Galerkin（Meshless Local Petrov-Galerkin，MLPG）法进行板结构的拓扑优化.基于带惩罚的各向同性固体微结构(Solid Isotropic Microstructure with Penalization，SIMP)的拓扑优化模型和优化准则法建立设计变量的优化修正方案.位移场和相对密度场均采用自然邻接点插值形函数进行离散插值.几种典型的拓扑优化算例证明该数值算法的正确性和有效性. 　　关键词：板结构； 拓扑优化； 无网格法； 自然邻接点插值； 带惩罚的各向同性固体微结构模型 　　中图分类号：??TB115.2?? 文献标志码：??A?? 　　Topology optimization of plate structure based on 　　meshless local Petrov-Galerkin method 　　LI Shunli, LONG Shuyao, LI Guangyao, DING Canhui 　　(State Key Laboratory of Advanced Design and Manufacture for Vehicle Body,Hunan University, Changsha 410082, China) 　　Abstract： To apply the advantage of meshless method into the structural topology optimization, the topology optimization of plate structure is performed by using the Meshless Local Petrov-Galerkin(MLPG) method. The optimization correction program of the design variables is built by the topology optimization model based on Solid Isotropic Microstructures with Penalization(SIMP) and the optimality criteria method. The natural neighbour interpolation shape function is used to discretize both displacement field and relative density field. Several typical topology optimization examples are used to verify the validity and effectiveness of the numerical method. 　　Key words： plate structure; topology optimization; meshless method; natural neighbour interpolation; solid isotropic microstructure with penalization model 　　 　　0 引 言 　　一种典型的拓扑优化问题目标为用给定的有限材料，在指定的空间内找到能支撑给定载荷的刚性最强的结构布局.自从BENDS??E等［1］提出均匀化法以来，许多拓扑优化方法被相继提出，其中具有代表性的有带惩罚的各向同性固体微结构(Solid Isotropic Microstructure with Penalization, SIMP)模型［2-4］、进化法［5］和水平集法［6］等. 　　目前，拓扑优化中主要的分析方法是有限元法.但有限元法也有其缺点，在处理如大变形和移动边界问题时需不断重新划分和重构网格以解决网格畸变和网格移动等问题.在拓扑优化中，精确的结构响应分析至关重要.近年来，为克服对网格的依赖性，许多学者致力于无网格方法的研究.目前，无网格方法已有数十种之多，其中具有代表性的有光滑粒子动力学法［7］、无单元伽辽金法［8］、无网格局部Petrov-Galerkin（Meshless Local Petrov-Galerkin，MLPG）法［9］和自然元法［10-11］等.近年来，蔡永昌等［12］和WANG等［13］将基于自然邻接点插值的MLPG法成功应用于求解弹性力学的平面问题. 　　为将无网格数值方法的优势集成到结构拓扑优化的应用中，本