基于偶应力模型的连续体结构拓扑优化设计.pdfVIP

第30卷第3期 固体力学学报 V01．30No．3 CHINESE SoLIDMECHANICS JoURNALOF 2009 2009年6月 June 基于偶应力模型的连续体结构拓扑优化设计。 刘书田” 苏文政 (大连理工大学工业装备结构分析国家重点实验室／T_程力学系，大连，116024) 摘 要经典连续介质理论不包含材料尺度参数，因而基于经典理论的结构拓扑优化无法显现尺度效应．本 文在偶应力理论的框架下，构造了四节点四边形离散偶应力单元，将传统的SIMP方法推广至偶应力介质．结果表 明，在以结构的最大刚度为目标的设计中，偶应力介质的最优结果取决于宏观结构尺寸与材料微结构尺寸(或者特 征长度)的比值，最优结果具有明显的尺度效应，具体为，二者比值较大将产生与传统理论相似的结构，而二者比值 相当则产生独特的偶应力主导的结构． 关键词偶应力，高阶介质，拓扑优化，有限元，尺度效应 的尺度效应和变形局部化现象[9’10。．微极介质理论， 0 引言 20世纪80年代逐渐兴起的拓扑优化技术与有成为微极理论．该理论假设物体质点既具有平动自 限元数值方法相结合在工程中得到了广泛的应 由度，又具有独立的转动自由度，质点之间既能够传 用[1。4]．基于均匀化理论[53以及基于变密度法[61的拓递力又能够传递力偶．需要指出的是，如果进一步假 扑优化技术是目前拓扑优化的主流方法，而后者由 设转动自由度依赖于平动自由度，等同于该点的局 于格式简单更是得到了广泛的应用．其主要的思想 部刚体转动，称为偶应力理论[1¨，或约束转动的微 为，假设单元的弹性模量和密度之间满足幂指数关 极理论．研究表明，生物的骨组织，颗粒夹杂复合材 系，通过改变单元的密度获得最优的结构拓扑，该方 料以及多孑L介质等在一定条件下需要采用高阶介质 MaterialwithPe- 法也称为SIMP(Solid Isotropie nalization)方法[1酒]．Bendsoe等[7]证明了该假设存高阶介质的材料常数． 在明确的物理意义． 高阶理论较之传统理论更为复杂，有限元法仍 传统的拓扑优化技术均基于经典连续介质(又 然是有效的求解方法．微极介质有限元可由现有的 称Cauchy介质)理论，该理论认为物体是具有三个有限元格式进行简单的推广得到[17|．而偶应力理论 自由度的质点的集合，质点之间只传递力而不传递 由于宏观转角和微观转角相等，由此导致总势能泛 力偶，是一种一阶的近似．对于大多数问题，宏观结 函中将出现位移的二阶导数，从而对位移提出了导 构的尺寸远大于材料微结构的尺寸，上述假定能够 数连续的要求，这就使偶应力有限元要求单元之间 准确描述物理问题；然而，材料的微观结构因为夹 满足C，连续．如何克服C，连续性要求是偶应力有 杂、缺陷以及晶格的存在总是不均匀的，材料的不均 限元构造的核心问题，国内外学者做了大量的工作， 匀性产生的尺度效应以及变形局部化现象是经典连 主要可以归结为两大类．第一类是直接构造满足C· 续介质力学无法解释的．在多孔材料／结构的一体化 连续的单元，例如Soh和Chenc”3借用薄板有限元 设计中，张卫红等L81给出的最优结果也具有明显的 单元格式进行组合叠加得到了一种18节点平面三 尺度关联现象．以微极理论、偶应力理论为代表的广 角形单元；第二类则是