[发明专利]一种多材料结构拓扑优化方法

摘要

本发明属于结构优化设计领域，并具体公开了一种多材料结构拓扑优化方法。该方法包括：对待处理对象的设计域进行有限元网格划分获得预设数量的网格单元并建立与实体材料一一对应的水平集函数；计算弹性矩阵、刚度矩阵和位移向量；计算速度场和敏度数，并根据敏度数对网格单元中的实体材料进行转换或删除，获得当前迭代下第一次优化后的水平集函数；利用速度场，更新获得当前迭代下第二次优化后的水平集函数；根据柔度和体积误差判别第二次优化的结果是否收敛。本发明克服了经典水平集拓扑优化方法对初始设计的依赖性，以及双向渐进优化方法无法获得光滑边界的缺陷，并在一定程度上提高了结构刚度，提供的方法稳定并且有效。

主权项

1.一种多材料结构拓扑优化方法，其特征在于，该方法包括如下步骤：(a)对待处理对象的设计域进行有限元网格划分获得预设数量的网格单元，利用N种实体材料和一种弱材料占据每个所述网格单元，其中被所述实体材料占据的区域定义为实体区域，被所述弱材料占据的区域定义为孔洞区域，然后根据所述设计域中实体区域和孔洞区域的分布，建立与所述实体材料一一对应的水平集函数；(b)根据待处理对象的材料属性，分别设定所述实体材料和弱材料的弹性模量和泊松比，以此计算所述实体材料和弱材料对应的弹性矩阵，并获得所述实体区域中各网格单元对应的刚度矩阵，然后计算所述设计域的位移向量；(c)利用所述实体材料的弹性矩阵和设计域的位移向量计算速度场，并利用所述网格单元的刚度矩阵和位移向量计算所述实体区域内各网格单元的敏度数，根据所述实体区域中各网格单元的敏度数对该网格单元中的实体材料进行转换或删除，从而更新所述实体区域和孔洞区域的分布，由此获得当前迭代下第一次优化后的水平集函数；(d)利用步骤(c)中获得的所述速度场，更新当前迭代下所述第一次优化后的水平集函数，获得当前迭代下第二次优化后的水平集函数；(e)计算当前迭代下所述第二次优化后的实体区域的柔度和体积误差，并分别与其预设误差阈值进行比较，判断所述柔度和体积误差是否均小于其所述预设误差阈值，若是，则将当前迭代下所述第二次优化后的水平集函数作为最优结构输出，若否，则返回步骤(b)。