结构拓扑优化设计的一个模型更加合理的120行程序.docVIP

结构拓扑优化目标和约束函数选取的探讨 易桂莲, 隋允康 （北京工业大学，100124） 摘 要：本文通过数值算例证明了在MCVC模型中结构的最优拓扑构型依赖于体积比约束值的大小。由此得到启发，从优化模型的构建角度，本文讨论了结构拓扑优化中目标和约束函数的选取问题。将与结构优化相关的目标和约束函数分为经济指标与性能指标，以MCVC模型代表经济指标约束下性能指标最优化建模，以MWDC模型代表性能约束下经济指标最小化建模，通过数值算例比较这两种模型，发现两个模型得到相似的结果，但是MWDC模型在建立工程问题模型时更具可行性，满足安全与经济的需求。 关键词：结构拓扑优化；结构柔顺度最小化;结构重量最轻化; 99行拓扑优化程序;120行拓扑优化程序 1 引言 结构拓扑优化的目的是寻求结构构件或子域的某种布局，在满足指定约束条件的情形下，使某种性能指标达到最优，它作为结构优化研究的重要分支，已经出现了相当多的研究方法。对于结构拓扑优化的研究进展，Bends?e和Sigmund[1]作了详尽的综述，Rozvany[2]也对数值建模方法在结构拓扑优化中的应用和发展进行了阐述。 在以往的研究中，通常能取为目标函数的有柔顺度[3-4]、刚度[5]、应变能[6]、重量[7-8]和造价[9]等，能取为约束条件的有：结构体积约束[3-4]、重量约束[10]、位移约束[11]、应力约束[12]、屈曲约束[13]和频率约束[14]等。在实际工程问题中，作为工程设计追求的目标，其选取通常是非常明确的。但是，如果只有目标函数而没有约束条件的限定，无论目标函数是什么，优化的结果全部为0或者全部为1，那么结构拓扑优化就毫无意义；如果在选定的目标函数下，约束条件不合理，那么无论采用哪种数值方法进行计算，结构就会在不断地增加或删除单元的过程中朝着与预期相反的方向发展。可见，合理地选取目标和约束条件是实现结构拓扑优化设计的关键。 在众多的结构拓扑优化模型中，传统而经典的建模是体积约束下结构柔顺度最小化 (MCVC，Minimum compliance with a Volume Constraint)问题，也是结构拓扑优化算法研究中使用最多的模型。这主要是因为MCVC模型在灵敏度分析上具有得天独厚的优势，备受研究者的青睐。MCVC问题通常是基于变密度法当中的SIMP(Solid Isotropic Material with Penalization)模型[1,3]： (1) 来进行建模。其中，下标e为单元号，为单元相对密度，也为设计变量。Ee和E0分别是人工材料和实际材料的Young’s模量。p是惩罚因子，通常取p=3.0。基于式(1)，MCVC问题的数学表达式可表示如下： (2) 其中，N是设计变量总数，是相对密度向量，为了避免刚度阵奇异，设定最小相对密度= 0.001。K、U、F分别是结构总刚度矩阵、总位移向量和总外力向量；ue是e单元的位移向量；是e单元刚度矩阵，是与实际材料的Young’s模量E0对应的单元刚度矩阵。和V*分别是使用材料的总体积和基结构的总体积，f是体积比值。 对于式(2)的求解，通常采用准则方法。Sigmund编译了一个99行拓扑优化程序，该程序可免费下载使用：http://www.topopt.dtu.dk/?q=node/2。Rozvany对99行拓扑优化程序予以高度评价[2]，称其“played an important role in SIMP’s general acceptance”。短短99行程序代码包含了有限元分析和拓扑优化计算，用户可根据实际问题定义结构形状大小、修改边界载荷条件和计算问题，该程序对初学者有很大的启发和帮助。 我们在试用该程序的过程中，计算了一批算例，发现MCVC模型中预先设定的体积比对于拓扑优化的结果有决定性的影响：不同的体积比约束值会产生不同的拓扑构型。 2 MCVC模型中最优解对体积约束值的依赖性 我们取图1中的四个算例： 图1 算例模型示意图 Fig. 1 Ground structures and boundary conditions for the four test examples 例1、见图1(a). MBB简支梁的一半(基结构为60×20×1，载荷F=1，Young’s模量E0=1.0)； 例2、见图1(b). 左端固支悬臂梁(基结构为32×20×1，载荷F=1，Young’s模量E0=1.0)； 例3、见图1(c). 左端固支带孔悬臂梁(基结构为45×30×1，载荷F=1，Young’s模量E0=1.0，孔的圆心位置为从左往右横向长度的1/3和纵向长度的1/2的交点处，半径为纵向长