位移约束集成化处理的连续体结构拓扑优化

为解决多工况下多位移约束的连续体结构拓扑优化问题,引入了K-S函数对位移约束进行集成化处理.在建立优化模型时,基于莫尔定理按ICM方法导出约束点位移与设计变量之间的近似显函数关系,然后采用Lagrange乘子法进行求解.给出三个算例对该方法进行验证,并与ESO法、BESO法、MD法以及均匀化方法的结果进行比较.结果表明:该方法计算效率较高,并且能够计算出更合理的结构拓扑.     

结构拓扑优化局部性能约束下轻量化问题的互逆规划解法

瞄准应力和疲劳两类局部性能约束的结构拓扑优化问题,概括为分部、化整和集成3种解法和交融的3种解法.类比应力约束推导了疲劳寿命情况,一为满疲劳公式,二为疲劳寿命约束全局化的结构寿命概念和相应解法.在倒寿命概念下,实现了疲劳寿命约束与应力约束的规格统一.补充和完整了已有的局部性能约束解法,属于单目标模型,有分部、化整、集成三种.基于互逆规划理论的定理2,提出了交融优化解法,是单目标与多目标模型的交替迭代,有分部-集成、化整-集成和集成-集成三种.上述6种解法皆基于ICM方法进行建模.算例表明,新提出的交融优化方法提高了求解效率.     

连续体结构屈曲约束的ICM方法拓扑优化

基于ICM(独立、连续、映射)方法解决具有屈曲约束的连续体拓扑优化问题。建立以结构重量为目标,以屈曲临界力为约束的拓扑优化模型;采用独立的连续拓扑变量,借助泰勒展式将目标函数作二阶近似展开;借助瑞利商、泰勒展式、过滤函数将约束化为近似显函数,避免了灵敏度的计算;将优化模型转化为对偶规划,并利用序列二次规划求解,减少了设计变量的数目,缩小了模型的求解规模。给出三个算例,结果表明:该方法可有效地解决屈曲约束的连续体拓扑优化问题,能够得到合理的拓扑结构,并有较高的计算效率。     

化可分离凸规划为对偶规划显式模型的普适解法研究比较

文章旨在提升对偶规划显式模型(dual programming-explicit model, DP-EM)的建模和求解的境界. DPEM模型从一类变量可分离凸规划的特点出发,突破了对偶目标二阶采用近似的定势,推导得出显式的对偶目标函数;应用于ICM方法求解连续体结构拓扑优化问题时,其求解效率比对偶序列二次规划方法 (DSQP)和可移动渐近线方法 (MMA)求解效率更高.文章进一步把常见的一类显式模型抽象为普适的可分离凸规划列式,在需要满足的一些条件下,转换为DP-EM模型,并且提出4种处理方法:(1)对偶变量迭代逼近法;(2)指数函数形式的解法;(3)幂函数形式的解法;(4)基于变换的精确解法.为了进行数值验证,做了广泛的计算,限于篇幅,文章列出了5个具有代表性的算例,除了算例1属于纯数学问题,其余4个算例皆基于ICM方法,分别对于位移、应力、疲劳等约束和破损-安全的连续体结构拓扑优化问题,基于所提出的方法进行建模和求解,都显示了所提出方法的普适性及更高的求解效率.工作的意义在于:(1)深度方面,加深了结构优化对偶解法的研究;(2)广度方面,对数学规划对偶理论的发展做出了新的贡献.     

频率约束板结构拓扑优化

本文利用ICM（独立、连续、映射）方法建立了频率约束下平板结构重量最轻的拓扑优化模型。采用指数函数作为单元重量、单元质量及单元刚度的过滤函数。通过瑞利商对刚度过滤函数倒变量的泰勒一阶展式,将频率约束近似显式化。利用对偶理论将含有大量设计变量的约束优化模型转化为易于求解的少设计变量拟无约束优化模型,通过序列二次规划将转化模型进行求解,提高了求解的效率。本文选择MSC.Patran&Nastran软件及PCL二次开发语言构架了平板结构频率约束拓扑优化问题的软件。数值算例表明：本文的方法具有迭代稳定性和收敛高效性。     

互逆规划的拓宽和深化及其在结构拓扑优化的应用

本文的工作涉及数学与力学两方面,数学方面:(1) 将数学规划论中新提出的互逆规划,从 s-m 型 (或称为 m-s 型) 发展出 s-s 型和 m-m 型互逆规划 (其中 s 意为单目标,m 意为多目标),从而使互逆规划的定义完备成为 3 种;(2) 从 KKT 条件审视互逆规划的两方面,得到了互逆规划双方求解涉及拟同构和拟同解的 3 个定理,并且予以证明,提供了在求解中对于互逆规划双方在求解中相互借鉴的理论基础;(3) 对一对互逆规划双方皆合理的情况和某一方不合理的情况,皆给出了求解策略和具体解法. 力学方面:(1) 给出结构优化设计模型合理与否的诠释;(2) 在互逆规划对结构拓扑优化的应用中,提出了不合理结构拓扑优化模型的求解策略;(3) 给出了借助 MVCC 模型 (多个柔顺度约束下的体积最小化) 解决 MCVC 模型 (对于给定体积下的多个柔顺度的最小化) 的途径,其中的建模基于 ICM (独立连续映射) 方法. 用 Matlab 编程给出了数值算例. 其中两个数学问题图示了互逆规划的双方关系;其中,结构拓扑优化问题是众多结构拓扑优化中的两个个案;数值结果均支持了本文提出的互逆规划理论.      

比较基于化整交融应力拓扑优化诸解法的效果

由于单元应力属于局部性能约束,导致相应的结构拓扑优化存在难以承受的大量约束条件;尽管化整方法极大地减少了约束数量,但是优化结果中有少数应力超限现象.为此,本文在应力约束的结构拓扑优化中,瞄准克服应力超限和提高求解效率两个目标,进行了探索.提出了乘子法及序列二次规划(SQP)法两种解法,首先在化整交融(即化整-集成)解法...     

多约束作用下连续体结构的拓扑优化

基于ICM(独立、连续、映射)方法建立了以结构重量最小为目标,以屈曲临界力、位移及应力三种约束同时作用的连续体拓扑优化模型:采用独立的连续拓扑变量,借助泰勒展式、过滤函数将目标函数作二阶近似展开。借助瑞利商、泰勒展式、过滤函数将屈曲约束化为近似显函数,借助于过滤函数,将位移约束用莫尔定理显式化;将应力这种局部性约束采用全局化策略进行处理,即借助第四强度理论、过滤函数将应力局部性约束转化为应变能约束,大大减少了灵敏度分析的计算量;将优化模型转化为对偶规划,减少了设计变量的数目,并利用序列二次规划求解,缩小了模型的求解规模。数值算例表明,ICM方法在解决屈曲、位移及应力三种约束共同作用的连续体拓扑优化问题上有优势。     

同时满足刚度和强度约束的框架拓扑优化

基于ICM（Independent Continuous Mapping,即独立、连续、映射）方法．对单元重量、单元许用应力和单元刚度分别引入不同的过滤函数,把0-1型离散拓扑变量转化为[0．1]区间上的连续变量．建立了拓扑变量连续的优化模型。借助满应力准则将应力约束转化为拓扑变量的动态下限．用单位虚载荷法将位移约束显式化．得到拓扑优化的近似显式模型。为了提高模型的求解效率,根据对偶理论求解原问题的对偶模型,通过在对偶空间迭代求解对偶模型得到原模型的解。引入结构非奇异、结构响应不被违背和结构重量不改变三个准则判断迭代收敛。并根据这三个准则自适应的调整折减系数来搜索最佳阁值。然后根据闻值将连续拓扑变量回归为0-1型离散拓扑变量。利用MSC／Nastran的开放性,借助MSC／Patran提供的PCL（Patran Command Language）开发环境．完成了满足刚度和强度的多变量的框架拓扑优化程序。算例结果表明．用ICM方法解决多变量框架拓扑优化问题是快速、有效的。     

基于改进的双向渐进结构优化法的应力约束拓扑优化

工程结构设计时经常需要限制最大名义应力,以避免发生断裂或疲劳破坏,一个有效的策略是采用拓扑优化方法. 常规的双向渐进结构优化法(bi-evolutionary structural optimization, BESO)不能有效求解应力约束拓扑优化问题,为此本文提出一种改进的双向渐进结构优化方法,处理体积和应力约束下的最小柔顺性问题. 引入基于K-S函数的全局应力度量,以减小大量局部应力约束引起的计算代价. 采用拉格朗日乘子法将应力约束函数引入到目标函数,然后由二分法确定合适的拉格朗日乘子的值使得应力约束得到满足. 而且,详细推导了基于BESO方法的应力约束拓扑优化模型及其灵敏度列式,最后通过三个典型拓扑优化算例验证改进方法的有效性. 为展示考虑应力约束的优点,将应力约束设计与传统的基于刚度的设计进行了比较. 结果表明, 改进的BESO方法优化迭代过程稳健,获得了边界灰度单元很少的清晰的拓扑构型,并实现了有效降低应力集中效应的设计.      

基于多相材料的稳态热传导结构轻量化设计

在多相材料的结构拓扑优化问题中,通常给定各相材料体积约束或材料总重量约束作为材料的控制用量.在结构轻量化设计的实际工程背景下,以结构总重量最小化为目标的优化模型具有明确的工程意义.针对含多相材料的稳态传热结构拓扑优化问题,提出了以结构总重量最小化为目标和给定热柔顺度为约束的多工况连续体结构拓扑优化建模方法.遵循独立连续映射建模方式,采用两类独立拓扑变量分别表征单元热传导矩阵和单元重量状态.推导了热柔顺度和总重量对设计变量的敏度,基于一阶和二阶泰勒展开得到各自的近似表达式.通过求解偏微分方程,实现了约束函数一次项过滤,消除了棋盘格现象和网格依赖性问题,并保证了约束方程在过滤后严格成立.建立的近似优化模型具有二次函数形式的目标函数和一次函数形式的约束函数.基于对偶序列二次规划方法对优化模型进行求解直至收敛.通过四个三维结构数值算例分析对比了热柔顺度约束限值、不同材料混合及多工况、多约束条件对优化结果的影响.数值算例结果表明,本文提出的优化方法在基于多相材料的多工况稳态热传导结构轻量化设计中具有可行性和有效性.     

近似评价函数确定运动极限及其在形状优化中的应用

针对二维连续体形状优化过程中解析灵敏度求解困难的情况,利用响应面方法将目标函数和约束函数近似显式化,建立序列二次规划模型。为了保证优化模型的可靠性,结合试验设计方法,以较少的结构分析代价构造约束响应面,从而能够快而稳定地收敛,本文重点研究并建立了二次近似评价函数用以计算自适应运动极限的策略。算例说明这种策略是有效而稳定的。     

斜拉桥索力优化的强次可行序列二次规划法

提出了一种斜拉桥索力优化的实用方法-强次可行序列二次规划法. 该方法通过建立斜拉 桥索力优化的非线性规划模型，以主梁和索塔的弯曲应变能为目标函数，斜拉索的索力为设 计变量，结构应力和索力为约束条件，计入大跨度斜拉桥各种几何非线性因素的影响，采用 强次可行序列二次规划算法进行优化求解，确定斜拉桥成桥合理状态的索力. 运用该方法对 某斜拉桥进行索力优化，结果表明该方法简单、有效.     

应力和位移约束下的板壳结构截面优化

将准则法和数学规划法相结合，根据满应力准则将应力约束化为动态下限，借助单位虚载荷法将位移约束转化为设计变量的近似显函数，建立了满足应力和位移约束的优化模型. 为了解决多变量的大型优化问题，根据对偶理论将上千设计变量的优化模型转化为几个变量的对偶模型，并通过二次规划求解. 以MSC/Nastran软件为结构分析的求解器，借助MSC/Patran软件为开发平台，完成了板壳结构截面优化程序. 程序完全和Patran及Nastran融为一体，在Patran中建立模型，利用Nastran分析计算，根据优化结果对设计变量调整，再用Nastran进行结构重分析，反复迭代直到结构重量收敛. 算例表明程序的合理性和有效性，能够满足工程设计的需要.     

面向应力约束的独立连续映射方法

大多数已有的拓扑优化研究为系统刚度最大化设计,尤其以体积比约束下的静态柔顺度最小化问题为典型.从工程角度出发,结构强度设计至关重要.以往的应力研究表明,应力约束拓扑优化存在着奇异性、约束数目庞大、高度非线性特性等诸多数值困难.为了实现应力约束下的拓扑优化设计,采用归一化p范数应力指标以减少单元应力约束数目.遵循独立连续映射建模方式,引入密度变量的倒变量函数作为设计变量.推导了应力约束函数和体积目标函数对设计变量的敏度,并基于一阶和二阶泰勒近似得到各自的显式表达式.通过构造的系列二次规划子问题,原拓扑优化问题采用序列二次规划算法高效求解.二维数值算例考察了结构刚度和强度设计结果的异同,以及不同应力约束上限值对应力约束拓扑优化结果的影响.通过提出方法与传统变密度法结果的比较,说明提出的独立连续映射方法在应力约束下具有可行性和有效性.优化结果也表明了考虑应力约束的连续体拓扑优化具有必要性.      

Plate/shell structure topology optimization of orthotropic material for buckling problem based on independent continuous topological variables

The purpose of the present work is to study the buckling problem with plate/shell topology optimiza-tion of orthotropic material.A model of buckling topology optimization is established based on the independent,con-tinuous, and mapping method, which considers structural mass as objective and buckling critical loads as constraints. Firstly, composite exponential function (CEF) and power function(PF)as filter functions are introduced to recognize the element mass,the element stiffness matrix,and the ele-ment geometric stiffness matrix.The filter functions of the orthotropic material stiffness are deduced. Then these fil-ter functions are put into buckling topology optimization of a differential equation to analyze the design sensitiv-ity.Furthermore,the buckling constraints are approximately expressed as explicit functions with respect to the design vari-ables based on the first-order Taylor expansion.The objective function is standardized based on the second-order Taylor expansion. Therefore,the optimization model is translated into a quadratic program.Finally,the dual sequence quadratic programming(DSQP)algorithm and the global convergence method of moving asymptotes algorithm with two different filter functions(CEF and PF)are applied to solve the opti-mal model.Three numerical results show that DSQP&CEF has the best performance in the view of structural mass and discretion.