• 综合性科技类中文核心期刊
    • 中国科技论文统计源期刊
    • 中国科学引文数据库来源期刊
    • 中国学术期刊文摘数据库(核心版)来源期刊
    • 中国学术期刊综合评价数据库来源期刊

基于独立连续变量和复合指数函数的位移约束平面连续体结构拓扑优化

叶红玲, 尹芳放, 王伟伟, 隋允康

叶红玲, 尹芳放, 王伟伟, 隋允康. 基于独立连续变量和复合指数函数的位移约束平面连续体结构拓扑优化[J]. 北京工业大学学报, 2016, 42(12): 1810-1817. DOI: 10.11936/bjutxb2016080005
引用本文: 叶红玲, 尹芳放, 王伟伟, 隋允康. 基于独立连续变量和复合指数函数的位移约束平面连续体结构拓扑优化[J]. 北京工业大学学报, 2016, 42(12): 1810-1817. DOI: 10.11936/bjutxb2016080005
YE Hongling, YIN Fangfang, WANG Weiwei, SUI Yunkang. Planar Continuum Structure Topology Optimization With Displacement Constraint Based on Independent Continuous Variables and Composite Exponential Function[J]. Journal of Beijing University of Technology, 2016, 42(12): 1810-1817. DOI: 10.11936/bjutxb2016080005
Citation: YE Hongling, YIN Fangfang, WANG Weiwei, SUI Yunkang. Planar Continuum Structure Topology Optimization With Displacement Constraint Based on Independent Continuous Variables and Composite Exponential Function[J]. Journal of Beijing University of Technology, 2016, 42(12): 1810-1817. DOI: 10.11936/bjutxb2016080005

基于独立连续变量和复合指数函数的位移约束平面连续体结构拓扑优化

基金项目: 国家自然科学基金资助项目(11072009);北京市教育委员会资助项目(KM201610005001);北京工业大学基础研究基金(001000514313003)
详细信息
    作者简介:

    叶红玲(1972—), 女, 副教授, 主要从事结构优化与多学科结构分析方面的研究, E-mail: yehongl@bjut.edu.cn

  • 中图分类号: TP31;O34

Planar Continuum Structure Topology Optimization With Displacement Constraint Based on Independent Continuous Variables and Composite Exponential Function

  • 摘要:

    为了进一步研究连续体结构拓扑优化模型的合理性和可行性,基于独立、连续、映射(independent continuous mapping,ICM)方法,在满足结构位移约束的条件下,通过引入复合指数形式过滤函数对位移约束下质量最小化(minimum weight with a displacement constraint,MWDC)模型进行了改进,建立了基于独立连续变量和复合指数函数的位移约束平面连续体结构拓扑优化模型,并进行了优化求解. 同时,利用M语言,基于Matlab软件平台,开发了相应的拓扑优化计算程序,并针对4种典型平面连续体结构进行了数值验证,分别比较分析了体积约束下的柔顺度最小化(minimum compliance with a volume constraint,MCVC)模型、MWDC模型以及改进的MWDC模型所得到的最优拓扑结构. 数值结果表明:采用复合指数形式过滤函数改进的MWDC优化模型迭代次数更少,优化求解计算效率更高.

    Abstract:

    To study the rationality and feasibility of the continuum structure topology optimization model, an improved minimum weight with a displacement constraint (MWDC) model by using exponential function was studied, which was based on the independent continuous mapping (ICM) method. A new topology optimization model for the problem of planar continuum structure with independent continuous variables and displacement constraints was established and solved. At the same time, a calculator program was developed and compiled based on the MATLAB in accordance with the new method. In addition, four typical numerical examples were adopted to verify the presented method. The topological results by taking advantage of MCVC model, MWDC model and improved MWDC model were compared with the view of structural mass and iterative numbers.Numerical results show that there is obvious advantage to solve the problem of planar continuum structure topology optimization with the improved MWDC optimization model in terms of calculation efficiency.

  • 根据结构设计变量范围的不同,结构优化问题可以分为尺寸优化、形状优化和拓扑优化3个层次. 结构拓扑优化相对于尺寸优化与形状优化具有更多的设计自由度,从而获得了更大的设计空间. 此外,拓扑优化在节省材料方面比尺寸优化和形状优化更为显著,因而可以取得更大的经济效益. 且拓扑优化待确定的参数较多,求解难度较大,因此更具有发展前景和挑战性.

    连续体结构拓扑优化设计是指在设计区域中搜索出一个优化子集,在满足预设约束条件(如应力约束、应变约束或频率约束等)的基础上,寻求使目标函数取得极小值时材料在平面或空间内的最优分布及结构的最佳传力路径的一类优化问题. 目前,连续体结构拓扑优化中已经建立了多种较为完善的拓扑表达形式与材料插值模型方法 [ 1] ,其中经典方法包括:均匀化方法 [ 2] (homogenization method)、相对密度法 [ 3] (relative denisity method)、水平集方法 [ 4] (level set method)、渐进结构优化法 [ 5] (evolutionary structural optimization method)、双向渐进优化法 [ 6] (bidirectional evolutionary structural optimization method)、相场法 [ 7] (phase field method)、节点变量法 [ 8] 、独立连续映射法 [ 9] (independent continuous mapping method,ICM)等.

    根据连续体结构拓扑优化问题中所包含的“三要素”即目标函数、设计变量和约束条件中的目标函数的属性,可以将其划分为不同的连续体结构拓扑优化模型. Yi等 [ 10] 按照结构的物理量将连续体拓扑优化问题划分为结构经济指标和结构性能指标两大类. 其中,结构经济指标是指结构的体积、质量或造价等;结构性能指标是指结构的力学性能,如结构柔顺度、位移、应力、固有频率、振幅等. 由此结构拓扑优化模型可以分为体积约束下的柔顺度最小化(minimum compliance with a volume constraint,MCVC)的结构性能指标模型和位移约束下质量最小化(minimum weight with a displacement constraint,MWDC)的结构经济指标模型. 模型不同,将导致其适用性和合理性也不尽相同,彭细荣等 [ 11] 针对包括MCVC模型和MWDC模型的连续体结构拓扑优化合理模型进行了一步探讨.

    本文基于ICM拓扑优化方法,针对MCVC模型和MWDC模型进一步深入研究,并引入新的过滤函数,即复合指数过滤函数对MWDC拓扑优化模型进行了改进. 采用单位虚载荷法将结构目标点位移显式化处理,并利用拉格朗日乘子法,推导出了采用新过滤函数的连续体结构MWDC模型的求解算法. 同时,基于MATLAB软件平台开发并编译了136行平面连续体结构拓扑优化程序,并利用该程序,通过4个数值算例对改进的模型进行了结果验证与比较.

    ICM拓扑优化方法 [ 12- 13] ,是以一种以独立于单元具体物理参数的变量来表征单元“有”与“无”的拓扑优化模型. 在优化过程中,利用过滤函数将本质上属于0/1的离散拓扑变量映射为[0,1]的连续变量,使得单元上的物理量同拓扑变量之间的关系由不确定、不连续、不可导变成确定、连续、可导,从而将离散拓扑优化问题转化成了光滑的数学模型. 优化迭代结束后,再通过逆映射将连续的变量反演为离散变量. ICM结构拓扑优化方法通常以结构整体质量最小(或结构整体体积最小)为优化目标,以结构响应作为约束条件,是一种典型的MWDC结构经济指标模型. 其以位移为约束的结构整体质量最小的拓扑优化模型可以表示为

    min W ( t ) = i = 1 N w i 0 s.t. s il s ¯ i 0 0 < t i 1 ( i = 1,2 , , N ; l = 1,2 , , L ) (1)

    式中: t= ( t 1 , t 2 , , t N ) T 为设计区域内单元拓扑设计变量; w i 0 为单元的固有质量; s il 为编号 i单元在工况 l条件下的结构最大位移值; s ¯ i 0 为编号 i单元的位移约束上限值; N为结构单元拓扑设计变量值; L为工况数量值 .

    典型的MCVC优化模型是基于固体各向同性材材惩罚(solid isotropic microstructured with penalization,SIMP)方法 [14 - 15] 的拓扑优化模型,其表达式为

    min c ( x ) = U T KU = e = 1 N ( x e ) p u T e k 0 u e s.t. V ( x ) V 0 f KU = F 0 < x min x e x max 1 (2)

    式中: c( x)为结构的柔顺度; x为设计变量; U为位移列向量; K为整体刚度矩阵; x e为单元设计变量; u e为单元位移列向量; k 0为单元刚度矩阵; V( x)为设计变量对应状态下的结构体积; V 0为初始结构体积; f为体积比; F为载荷列向量; x minx max分别为单元内设计变量的下限值与上限值; p为惩罚因子 .

    由于工程中最常用的结构物理量指标为结构质量(或体积)、强度、刚度及稳定性等,因此相比于MCVC模型,MWDC模型更符合工程问题的惯常提法. 在MWDC模型的发展中,为了便于与由Sigmund [ 15] 所提出的典型的MCVC优化模型所得最优拓扑结构进行分析比较,Yi等 [ 10] 针对幂函数形式过滤函数的MWDC结构经济指标模型进行了研究.

    在ICM优化方法中,过滤函数发挥着极其重要的作用. 过滤函数不仅可以实现对于拓扑变量的过滤和筛选,完成对拓扑变量由连续模型向离散模型的回归,而且过滤函数在建模中还起到了识别几何或物理量的作用,比如

    w i=f w( t i ) w i 0 (3)

    k i=f k( t i ) k i 0 (4)

    所示的识别单元质量、单元刚度矩阵等 . 式中: w i k i 分别为 i单元全程变化的质量和刚度矩阵; w i 0 k i 0 分别为 i单元的固有质量和刚度矩阵 . f w( t i )、 f k( t i )分别为结构质量过滤函数和结构刚度过滤函数 .

    综上所述,可以得到引入过滤函数的位移约束下结构质量最小的连续体结构拓扑优化MWDC模型表达式为

    min W = i = 1 N f w ( t i ) w i 0 s.t. f k ( t i ) s il s ¯ i 0 0 < t i 1 ( i = 1,2 , , N ; l = 1,2 , , L ) (5)

    式中: f w( t i )和 f k( t i )为过滤函数 .

    在连续体结构拓扑优化问题中,选取不同形式的过滤函数,直接会导致优化求解计算效率与优化结果的差异 . ICM方法中所采用的过滤函数主要包括幂函数形式的过滤函数和复合指数函数形式的过滤函数等 [ 13] ,其数学表达式分别为

    f( t i α ) = t i α , α≥1

    f( t i ) = e t i / γ - 1 e 1 / γ - 1 , γ > 0

    式中 αγ为常数,可以通过数值试验的方法获得 .

    由文献[16]可知,对于幂函数形式的过滤函数和复合指数形式的过滤函数2种常见过滤函数而言,在设计变量[0,1]区间中,采用复合指数形式的过滤函数相比于幂函数形式的过滤函数更易使对应设计变量迅速趋近于1 . 本文将研究采用

    f w( t i ) = e t i / γ w - 1 e 1 / γ w - 1 , γ w > 0(6)

    f k( t i ) = e t i / γ k - 1 e 1 / γ k - 1 , γ k > 0(7)

    所示的复合指数函数形式过滤函数对式(1)中单元质量与单元刚度矩阵进行过滤,从而进一步探讨过滤函数对MWDC模型拓扑优化结果的影响 .

    由莫尔定理可知,结构任意节点在某一方向上的广义位移可表示为

    s= i = 1 N s i= i = 1 N σ T ig ε i vd v(8)

    式中: s i 为第 i个单元对于结构目标点位移的影响值; σ i g为设计区域结构中实载荷下的单元应力向量; ε i v为设计区域结构中虚载荷下的单元应变向量 .

    根据虚功原理“外力在虚位移上做的功等于内力在虚位移导致的虚变形上所做的虚功”,可以得到

    s= ( P R i ) T s i v (9)

    s= i = 1 N ( P R i ) T s i v (10)

    式中: P R i 为实载荷下第 i个单元节点力向量; s v i 为虚载荷下的第 i个单元节点位移向量 .

    根据有限元原理,设计域结构的单元刚度

    K i s v i = P v i (11)

    式中: K i 为第 i个单元的刚度矩阵; P v i 为虚载荷下结构单元节点力列向量 .

    引入复合指数型过滤函数,则单元刚度矩阵可以表示为

    K i = k i 0 E i= K i 0 e t i / γ k - 1 e 1 / γ k - 1 (12)

    将设计域结构的单元刚度代入虚功原理公式(10),可以得到设计区域中目标点位移为

    s= i = 1 N ( P R i ) T ( K i ) - 1 P v i (13)

    将得到的单元刚度矩阵代入设计区域中目标点位移公式,可以得到目标点位移的显式化表达式为

    s= i = 1 N ( P R i ) T ( K i 0 ) - 1 P v i e 1 / γ k - 1 e t i / γ k - 1 = i = 1 N D i 0 e 1 / γ k - 1 e t i / γ k - 1 (14)

    其中,由于静定结构的位移不依赖于节点力,因此,单元对目标位移点的作用分量常系数 D i 0 不依赖于实载荷下第 i个单元节点力向量 P i 和虚载荷下结构单元节点力列向量 P v i .

    通过上述转化,借助刚度矩阵过滤函数,位移约束由设计变量的隐函数转化为了显函数 .

    将上述得到的目标点位移的显式化表达式(14)代入如式(1)所示的 MWDC结构性能指标位移约束下结构整体质量最小的优化模型,可以得到改进的 MWDC结构性能指标拓扑优化模型为

    min W ( t ) = i = 1 N f w ( t i ) w i 0 s.t. i = 1 N D i 0 / f k ( t i ) s ¯ 0 < t min t i 1 ( i = 1,2 , , N ) (15)

    为了获得上述优化模型的设计变量最优解,利用拉格朗日乘子法进行计算分析. 引入拉格朗日乘子 λ,得到拉格朗日方程为

    L( t i , λ) = - (16)

    拉格朗日函数的一阶导数可以表示为

    L t i = ( 1 - e 1 / γ w ) e t i / γ w γ w ( e t i / γ w - 1 ) 2 w i 0 - λ D i 0 ( 1 - e 1 / γ k ) e t i / γ k γ k ( e t i / γ k - 1 ) 2 =0(17)

    L λ = i = 1 N D i 0 e 1 / γ k - 1 e t i / γ k - 1 + s ¯ =0(18)

    消去拉格朗日乘子 λ可得最终解为

    t i = i = 1 N D i 0 e 1 / γ k - 1 e t i / γ k - 1 [ 1 / ( γ k + 1 ) ] ( w i 0 ) γ w / ( γ w + 1 ) + e ( 1 / γ w - 1 ) s ¯ e 1 / γ k ( D i 0 ) ( 1 / γ k - 1 ) (19)

    至此,完成了该改进的MWDC结构性能指标拓扑优化模型设计变量的分析求解 .

    为了保持拓扑优化模型程序的可比性,采用求解MCVC模型相同的优化收敛准则,即第 n+1次迭代得到的最优解 t i ( n + 1 ) 与第 n次迭代所得到的最优解 t i ( n ) 满足

    max(max( | t i ( n + 1 ) - t i ( n ) |))≤0 .01(20)

    基于以上位移约束下改进的MWDC优化模型,将整个求解过程利用M语言在Matlab软件平台进行了编译与实现,其结构拓扑优化程序流程如 图1所示.

    图  1  改进的MWDC优化模型拓扑优化流程
    Figure  1.  Improved MWDC optimization model flow chart

    为了进一步验证该改进的MWDC优化模型的有效性与可行性,本文给出4个经典数值算例,分别对采用幂函数形式过滤函数的MCVC模型、传统的MWDC模型与采用复合指数形式过滤函数的改进的MWDC模型所得到的最优拓扑结构进行分析比较.

    算例1左端固支悬臂梁,其结构如 图2所示. 基结构为32mm×20mm×1mm,载荷为 F=1N,弹性模量 E 0=1.0MPa,泊松比为 μ=0.3.

    图  2  左端固支悬臂梁结构示意
    Figure  2.  Diagram of cantilever beam structure

    采用MCVC模型、传统的MWDC模型以及本文改进的MWDC模型的3种优化模型所得到的平面连续体最优拓扑结构及迭代次数分别如 表1所示,与该结构相对应的拓扑优化结构质量迭代曲线以及结构位移迭代过程曲线如 图3、4所示.

    表  1  左端固支悬臂梁平面连续体结构拓扑优化结果
    Table  1.  Plane continuum structure topology optimization results for cantilever beam
    优化模型 MCVC模型 [ 15] MWDC模型 [ 10] 本文改进的MWDC模型
    输入参数 Top99(32,20,0.4,3.0,1.2) Top120(32,20,57.35,3.0,1.2) Top120(32,20,57.35,3.0,1.2)
    优化结果
    位移约束/mm 57.35 57.35
    结构质量/kg 250.89 255.81 238.07
    迭代次数 71 70 35
    下载: 导出CSV 
    | 显示表格
    图  3  左端固支悬臂梁平面连续体结构质量迭代曲线
    Figure  3.  Massiterative curve for cantilever beam structure
    图  4  左端固支悬臂梁平面连续体结构位移迭代曲线
    Figure  4.  Displacement iterative curve for cantilever beam structure

    根据以上优化结果,采用MCVC模型、传统的MWDC模型以及改进的MWDC模型3种优化模型所得到的最优拓扑结构基本相同. 其中,采用MCVC模型与传统的MWDC模型所得到的优化后结构总质量差别较小,而采用改进的MWDC模型所得到的优化后结构总质量较MCVC模型与传统的MWDC模型分别降低5.11%和6.93%. 此外,采用3种模型进行拓扑优化求解的迭代次数分别为71、70、35次. 可见,采用改进后的MWDC模型较MCVC模型与传统的MWDC模型在优化求解效率上分别实现了50.70%和50.00%的显著提升.

    算例2左端固支带孔悬臂梁,其结构如 图5所示. 结构为45mm×30mm×1mm,载荷为 F=1N,弹性模量 E 0=1.0MPa,泊松比为 μ=0.3,孔的位置为从左往右横向长度的1/3和纵向长度的1/2,半径为纵向长度的1/3.

    图  5  左端固支带孔悬臂梁结构示意
    Figure  5.  Diagram of cantilever beam structure with a hole

    采用MCVC模型、传统的MWDC模型以及本文改进的MWDC模型的3种优化模型所得到的平面连续体最优拓扑结构及迭代次数分别如 表2所示,与该结构相对应的拓扑优化结构质量迭代曲线以及结构位移迭代过程曲线如 图6、7所示.

    表  2  左端固支带孔悬臂梁平面连续体结构拓扑优化结果
    Table  2.  Plane continuum structure topology optimization results for cantilever beam with a hole
    优化模型 MCVC模型 [ 15] MWDC模型 [ 10] 本文改进的MWDC模型
    输入参数 Top99(45,30,0.5,3.0,1.5) Top120(45,30,60,3.0,1.5) Top120(45,30,60,3.0,1.5)
    优化结果
    位移约束/mm 60.00 60.00
    结构质量/kg 633.25 673.29 620.09
    迭代次数 34 51 11
    下载: 导出CSV 
    | 显示表格
    图  6  左端固支带孔悬臂梁平面连续体结构质量迭代曲线
    Figure  6.  Mass iterative curve for cantilever beam structure with a hole
    图  7  左端固支带孔悬臂梁平面连续体结构位移迭代曲线
    Figure  7.  Displacement iterative curve for cantilever beam structure with a hole

    根据以上优化结果可知,采用MCVC模型、传统的MWDC模型以及改进的MWDC模型3种优化模型所得到的最优拓扑结构基本相同. 其中,采用改进的MWDC模型所得到的优化后结构总质量较MCVC模型与传统的MWDC模型分别降低了2.08%和7.90%. 此外,采用3种优化模型进行拓扑优化求解的迭代次数分别为34、51、11次. 可见,采用改进的MWDC模型较MCVC模型与传统的MWDC模型在优化求解效率方面具有显著提升.

    算例3左端固支右端中部加载悬臂梁,其结构如 图8所示. 基结构为80mm×50mm×1mm,载荷为 F=9kN,弹性模量 E 0=1.0MPa,泊松比为 μ=0.3.

    图  8  左端固支右端中部加载悬臂梁示意
    Figure  8.  Diagram of cantilever beam structure with a load at the middle of right eadge

    采用MCVC模型、传统的MWDC模型以及本文改进的MWDC模型的3种优化模型所得到的平面连续体最优拓扑结构及迭代次数分别如 表3所示,与该结构相对应的拓扑优化结构质量迭代曲线以及结构位移迭代过程曲线如 图9、10所示. 由 表3可见,采用改进的MWDC模型迭代次数最少,结构质量最轻,由 表3可见.

    表  3  左端固支右端中部加载悬臂梁平面连续体结构拓扑优化结果
    Table  3.  Plane continuum structure topology optimization results for cantilever beam with a load at the middle of right eadge
    优化模型 MCVC模型 [ 15] MWDC模型 [ 10] 本文改进的MWDC模型
    输入参数 Top99(80,50,0.5,3.0,1.5) Top120(80,50,0.35,3.0,1.5) Top120(80,50,0.35,1.9,1.5)
    优化结果
    位移约束/mm 0.35 0.35
    结构质量/kg 1966.02 1990.17 1926.08
    迭代次数 78 68 12
    下载: 导出CSV 
    | 显示表格
    图  9  左端固支右端中部加载悬臂梁结构质量迭代曲线
    Figure  9.  Mass iterative curve for cantilever beam structure with a load at the middle of right eadge
    图  10  左端固支右端中部加载悬臂梁结构位移迭代曲线
    Figure  10.  Displacement iterative curve for cantilever beam structure with a load at the middle of right eadge

    算例4 MBB简支梁的一半,其结构如 图11所示. 基结构为60mm×20mm×1mm,载荷为 F=1N,弹性模量 E 0=1.0MPa,泊松比为 μ=0.3.

    图  11  MBB简支梁结构示意图
    Figure  11.  Diagram of MBB beam structure

    采用MCVC模型、MWDC模型以及改进的MWDC模型的3种优化模型所得到的平面连续体最优拓扑结构及迭代次数分别如 表4所示,与该结构相对应的拓扑优化结构质量迭代曲线如 图12所示. 可以看到采用改进的MWDC模型较MCVC模型与传统的MWDC模型在优化求解效率方面具有显著提升,迭代次数最少.

    表  4  MBB平面梁连续体结构拓扑优化及迭代次数
    Table  4.  Plane continuum structure topology optimization for MBB beam and the number of iterations
    优化模型 MCVC模型 MWDC模型 改进的MWDC模型
    输入参数 Top99(60,20,0.5,3.0,1.5) Top120(60,20,203.30,3.0,1.5) Top120(60,20,203.30,3.0,1.5)
    优化结果
    位移约束/mm 203.30mm 203.30mm
    结构质量/kg 582.57 599.15 606.23
    迭代次数 94 121 18
    下载: 导出CSV 
    | 显示表格
    图  12  MBB平面梁结构质量迭代曲线
    Figure  12.  Mass iterative curve for MBB plane beam structure

    1) 采用改进的MWDC结构经济指标模型所得到的最优拓扑结构与采用MCVC结构性能指标模型和传统的MWDC结构经济指标模型所得到的最优拓扑结构基本相同. 因此,改进的MWDC结构经济指标模型能够满足解决以目标点位移作为约束条件的拓扑优化问题.

    2) 采用改进的MWDC结构经济指标模型所得到的优化后结构的总质量较采用MCVC结构性能指标模型和传统的MWDC结构经济指标模型均有所降低.

    3) 在拓扑优化求解收敛速率方面,采用改进的MWDC结构经济指标模型较采用MCVC结构性能指标模型和传统的MWDC结构经济指标模型均实现了显著提升,大幅降低了结构拓扑优化问题的求解成本.

    The authors have declared that no competing interests exist.
  • 图  1   改进的MWDC优化模型拓扑优化流程

    Figure  1.   Improved MWDC optimization model flow chart

    图  2   左端固支悬臂梁结构示意

    Figure  2.   Diagram of cantilever beam structure

    图  3   左端固支悬臂梁平面连续体结构质量迭代曲线

    Figure  3.   Massiterative curve for cantilever beam structure

    图  4   左端固支悬臂梁平面连续体结构位移迭代曲线

    Figure  4.   Displacement iterative curve for cantilever beam structure

    图  5   左端固支带孔悬臂梁结构示意

    Figure  5.   Diagram of cantilever beam structure with a hole

    图  6   左端固支带孔悬臂梁平面连续体结构质量迭代曲线

    Figure  6.   Mass iterative curve for cantilever beam structure with a hole

    图  7   左端固支带孔悬臂梁平面连续体结构位移迭代曲线

    Figure  7.   Displacement iterative curve for cantilever beam structure with a hole

    图  8   左端固支右端中部加载悬臂梁示意

    Figure  8.   Diagram of cantilever beam structure with a load at the middle of right eadge

    图  9   左端固支右端中部加载悬臂梁结构质量迭代曲线

    Figure  9.   Mass iterative curve for cantilever beam structure with a load at the middle of right eadge

    图  10   左端固支右端中部加载悬臂梁结构位移迭代曲线

    Figure  10.   Displacement iterative curve for cantilever beam structure with a load at the middle of right eadge

    图  11   MBB简支梁结构示意图

    Figure  11.   Diagram of MBB beam structure

    图  12   MBB平面梁结构质量迭代曲线

    Figure  12.   Mass iterative curve for MBB plane beam structure

    表  1   左端固支悬臂梁平面连续体结构拓扑优化结果

    Table  1   Plane continuum structure topology optimization results for cantilever beam

    优化模型 MCVC模型 [ 15] MWDC模型 [ 10] 本文改进的MWDC模型
    输入参数 Top99(32,20,0.4,3.0,1.2) Top120(32,20,57.35,3.0,1.2) Top120(32,20,57.35,3.0,1.2)
    优化结果
    位移约束/mm 57.35 57.35
    结构质量/kg 250.89 255.81 238.07
    迭代次数 71 70 35
    下载: 导出CSV

    表  2   左端固支带孔悬臂梁平面连续体结构拓扑优化结果

    Table  2   Plane continuum structure topology optimization results for cantilever beam with a hole

    优化模型 MCVC模型 [ 15] MWDC模型 [ 10] 本文改进的MWDC模型
    输入参数 Top99(45,30,0.5,3.0,1.5) Top120(45,30,60,3.0,1.5) Top120(45,30,60,3.0,1.5)
    优化结果
    位移约束/mm 60.00 60.00
    结构质量/kg 633.25 673.29 620.09
    迭代次数 34 51 11
    下载: 导出CSV

    表  3   左端固支右端中部加载悬臂梁平面连续体结构拓扑优化结果

    Table  3   Plane continuum structure topology optimization results for cantilever beam with a load at the middle of right eadge

    优化模型 MCVC模型 [ 15] MWDC模型 [ 10] 本文改进的MWDC模型
    输入参数 Top99(80,50,0.5,3.0,1.5) Top120(80,50,0.35,3.0,1.5) Top120(80,50,0.35,1.9,1.5)
    优化结果
    位移约束/mm 0.35 0.35
    结构质量/kg 1966.02 1990.17 1926.08
    迭代次数 78 68 12
    下载: 导出CSV

    表  4   MBB平面梁连续体结构拓扑优化及迭代次数

    Table  4   Plane continuum structure topology optimization for MBB beam and the number of iterations

    优化模型 MCVC模型 MWDC模型 改进的MWDC模型
    输入参数 Top99(60,20,0.5,3.0,1.5) Top120(60,20,203.30,3.0,1.5) Top120(60,20,203.30,3.0,1.5)
    优化结果
    位移约束/mm 203.30mm 203.30mm
    结构质量/kg 582.57 599.15 606.23
    迭代次数 94 121 18
    下载: 导出CSV
  • [1] 罗震,陈立平,黄玉盈,等.连续体结构的拓扑优化设计[J].力学进展,2004,34(4):463-476.

    LUOZ,CHEN LP,HUANG YY,et al.Topological optimization design for continuum structures[J].Advances in Mechanics,2004,34(4):463-476. (in Chinese)

    [2] BENDSØE MP,KIKUCHIN.Generating optimal topologies in structural design using a homogenization method[J].Computer Methods in Applied Mechanics & Engineering,1988,71(2):197-224.
    [3] MLEJNEK HP,SCHIRRMACHERR.An engineer's approach to optimal material distribution and shape finding[J].Computer Methods in Applied Mechanics & Engineering,1993,106(1/2):1-26.
    [4] MAKHIJAD,MAUTEK.Numerical instabilities in level set topology optimization with the extended finite element method[J].Structural & Multidisciplinary Optimization,2014,49(2):185-197.
    [5] WANG BQ,WANG BL,HUANG ZY.Topology optimization for constrained layer damping plates using evolutionary structural optimization method[J].Advanced Materials Research,2014,894:158-162.
    [6] 匡兵,李应弟,刘夫云,等.基于单元密度进化步长控制的双向渐进结构优化方法[J]. 计算力学学报,2016(1):15-21.

    KUANGB,LI YD,LIU FY,et al.Bi-directional evolutionary structural optimization method based on control for evolutionary step length of element density[J]. Chinese Journal of Computational Mechanics,2016(1):15-21. (in Chinese)

    [7] TAKEZAWAA,NISHIWAKIS,KITAMURAM.Shape and topology optimization based on the phase field method and sensitivity analysis[J].Journal of Computational Physics,2010,229(7):2697-2718.
    [8] 占金青,杨康,黄志超.基于节点变量法的连续体结构拓扑优化设计[J].农业机械学报,2014,45(9):329-332.

    ZHAN JQ,YANGK,HUANG ZC.Topology optimization of continuum structures using node variable method[J].Transactions of the Chinese Society for Agricultural Machinery,2014,45(9):329-332. (in Chinese)

    [9] 隋允康.建模变换优化——结构综合方法新进展[M].大连:大连理工大学出版社,1996:2-65.
    [10] YI GL,SUI YK.An adaptive approach to adjust constraint bounds and its application in structural topology optimization[J].Journal of Optimization Theory & Applications,2016,169(2):656-670.
    [11] 彭细荣,隋允康.对连续体结构拓扑优化合理模型的再探讨[J]. 固体力学学报,2016(2):181-191.

    PENG XR,SUI YK.A further discussion on rational topology optimization models for continuum structures[J]. Chinese Journal of Solid Mechanics,2016(2):181-191. (in Chinese)

    [12] SUIY,YANGD.New method for structural topological optimization based on the concept of independent continuous variables and smooth model[J].Acta Mechanica Sinica,1998,14(2):179-185.
    [13] 隋允康,叶红玲.连续体结构拓扑优化的ICM方法[M].北京:科学出版社,2014:27-59.
    [14] BENDSØE MP,SIGMUNDO.Material interpolation schemes in topology optimization[J].Archive of Applied Mechanics,1999,69(9/10):635-654.
    [15] SIGMUNDO.A 99 line topology optimization code written in Matlab[J].Structural and Multidisciplinary Optimization,2001,21(2):120-127.
    [16] 叶红玲,沈静娴,隋允康.频率约束的三维连续体结构动力拓扑优化设计[J].力学学报,2012,44(6):1037-1045.

    YE HL,SHEN JX,SUI YK.Dynamic topological optimal design of three-dimensional continuum structures with frequencies constraints[J].Chinese Journal of Theoretical and Applied Mechanics,2012,44(6):1037-1045. (in Chinese)

图(12)  /  表(4)
计量
  • 文章访问数:  31
  • HTML全文浏览量:  3
  • PDF下载量:  13
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-07-31
  • 网络出版日期:  2023-05-28
  • 刊出日期:  2016-12-09

目录

/

返回文章
返回