基于阶跃函数的连续体拓扑优化问题.docxVIP

基于阶跃函数的连续体拓扑优化问题 1 基于边界类条件的拓扑优化 作为结构优化的高水平领域，在国内外科学家的共同努力下，做出了迅速的发展。目前在连续体拓扑优化中占据主导地位的是均匀化法和SIMP (Solid Isotropic Material with Penalization) 法。因为它们都是主要研究单元的有和无, 所以把它们统称为单元类方法。这类方法从本质上讲是一个0-1整数规划问题。但是由于这一问题是不能保证存在最优解的病态问题, 所以不能直接求解。为了求解, 通常采用松弛设计变量和增加约束两类办法。均匀化方法是把拓扑优化的基结构分成许多带有孔洞的胞元, 以每一个胞元孔洞的几何参数为设计变量, 从而把拓扑优化问题转化为尺寸优化问题来求解。SIMP法则是把每一个单元内的材料属性作为常数, 单元密度作为设计变量, 通过对单元密度惩罚的方式来得到趋近0-1的拓扑结构。此外, Beckers等人通过引进周长约束直接求解整数0-1规划, 隋允康等人把拓扑变量从尺寸优化变量中提取出来的独立映射反演 (ICM) 方法也都属于单元类方法。这一类方法在拓扑优化的研究过程中起到了重大的作用, 现在已经被广泛用于求解多种约束, 多种物理场耦合等很多拓扑优化新兴领域之中。但是由于它们都是以结构的单元为研究对象, 所以不容易处理拓扑构型演化过程中的边界问题。 而以边界移动来确定结构构型的边界类算法在边界的处理上则具有非常明显的优势。但是标准的边界类算法只是通过边界形状改变的灵敏度来确定构型边界的移动速度, 不能起到改变结构拓扑的作用, 所以大多数情况下它只能作为一种形状优化算法来应用。为了能够更好地应用边界类算法的优点来处理拓扑优化中的边界问题, J.A. Sethian和Andreas Wiegmann把用于描述材料边界运动的水平集概念引入到结构优化之中。水平集算法的主要过程是由优化的目标和约束来确定边界的移动速度, 然后由Hamilton-Jacobi方程求解水平集函数, 从而确定结构的拓扑构型。由于水平集函数可以使边界在移动的同时产生拓扑改变, 所以水平集方法在拓扑优化中也取得了很好的数值结果。但是因为水平集方法不能使远离边界的结构内部的拓扑发生改变, 所以当初始构型取的不好时不易收敛到最优解。尽管水平集方法在拓扑优化的应用中还存在一些问题, 但是它可以使拓扑优化构型的边界处理变得简单。 最近, Belytschko和郭旭分别提出了基于水平集函数的求解方法。他们把有限元离散后的水平集函数节点值作为拓扑设计变量, 通过阶跃函数把水平集函数转化为0-1形式建立优化模型, 然后通过有限元离散和优化算法进行求解, 从而确定结构的拓扑构型。其实, 这种算法从本质上讲是求解整数0-1规划的一种连续化方法。数值算例证明了这种方法的有效性。但是在这两篇文献中也可以发现一对矛盾。在这两篇文献中都引进了三角函数来光滑阶跃函数。如果其中的光滑区域取的比较小, 问题的求解过程就接近于水平集方法。但是通过对拓扑导数的概念进行分析后, 可以看出这种接近于水平集的求解方法很难改变离边界较远的结构内部的拓扑。如果光滑区域取的比较大, 拓扑导数就可以在离边界比较远的结构内部发挥作用, 结构的拓扑改变就变得更加容易, 优化初始阶段数值收敛速度也就越快, 但是这样做却会影响收敛到最终拓扑构型的速度。 为了解决这个矛盾, 在文的基础上, 采用sigmoid函数来光滑阶跃函数, 并通过对sigmoid函数中陡度参数的调整来更加方便地控制拓扑优化的进程。这样既可以在拓扑优化的初始阶段采用较小的陡度参数, 使结构内部的拓扑改变起来更加容易, 加快收敛速度;又可以在拓扑优化的后期采用比较大陡度参数来确保以较快的速度收敛到最优拓扑构型。 2 等效积分拟合 拓扑优化模型 在连续体拓扑优化模型中, 需要引进一个0-1变量的特征函数χ(x) 来确定基结构内的点是否属于最优拓扑构型 χ(x)={0x∈D\Ω1x∈Ω(1) 式中:D代表基结构;Ω代表拓扑构型所占区域。 对于二维线弹性体, 不考虑体积力, 以整体柔顺性为目标, 材料体积为约束的连续体拓扑优化模型P原为 Μinimizeu,?L(u,χ)=∫?ΩitudS=∫DEijklεij(u)εkl(u)χ(x)dΩ subject to: a (u,v,χ) =L(v,χ) ,u|?Dn=u0, ?v∈U V(χ) =∫Dχ(x) dΩ≤Vmax(2) 其中a(u,v,χ) =∫DEijklεij(u)εkl(v)χ(x) dΩ,L(v,χ) =∫?ΩitvdS,L(u,χ) =∫DEijklεij(u)εkl(u)χ(x) dΩ为由最小位能原理推出的柔顺性在整个基结构上的表达,a(u,v,χ) =L(v,χ) 为平衡方程的等效积分弱形式