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

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

基于响应面方法的桁架截面敏度分析和优化

把响应面方法引入桁架截面优化中,将应力和位移约束近似表达为桁架截面倒变量的线性函数。为拟合响应面,基于中心复合和单纯形试验设计方法开发了中心对称和拟单纯形试验设计两种方法,既可保证约束近似精度,又降低了结构分析计算量。对于桁架结构重量目标函数,直接推出倒变量的二阶形式,以桁架重量最小为目标的优化问题构造为标准的二次规划模型。将响应面方法计算的位移对设计变量的敏度与莫尔积分方法的近似显式进行了对比。以MSC.Patran为平台的数值算例表明:结合响应面方法,应用序列二次规划对问题进行寻优,其收敛精度及稳定性都可获得保障。     

基于响应面方法的结构碰撞优化

应用响应面方法可将结构响应表达为设计变量的显式形式．为拟合响应面，基于中心复合设计和单纯形设计简化出了中心对称和拟单纯形设计，其既可以保证约束近似精度，同时降低计算量．针对常见优化问题，使用响应面方法将其构造为线性或者非线性约束优化模型，并应用Matlab优化工具箱求解．方管碰撞优化数值算例表明文中提出的试验设计方法及优化模型有较好的求解效率和精度．     

应力和位移约束下桁架拓扑优化的有无复合体方法

按照独立连续拓扑变量的思想,将有无复合体模型由构造应力约束问题发展为应力与位移约束问题,用于桁架结构拓扑优化．基于该文提出的位移约束与拓扑变量关系的近似显式,建立了具有应力、位移约束下桁架结构拓扑优化的有无复合体模型,用序列二次规划算法求解．数值结果令人满意．这一工作表明独立连续拓扑变量的提出对于结构拓扑优化的研究是有价值的。     

频率约束板结构拓扑优化

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

板壳结构的应力约束拓扑优化设计

应用ICM 方法求解应力约束板壳结构拓扑优化问题,建立了寻求应力约束下结构重量极小化,每个设计变量控制多个单元的板壳结构拓扑优化近似显式的ICM 模型. 依据畸变能理论,将应力约束转化为畸变能约束,减少了约束数目. 采用精确对偶映射下的序列二次规划算法进行求解. 以MSC.Patran 及MSC.Nastran 软件作为二次开发平台,应用PCL 语言实现本文算法. 算例对于设计变量数等于单元数时的情况进行了计算,表明该方法有效可行.     

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

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

多约束膜结构截面优化及在MSC/Nastran上的程序实现

将准则法和数学规划法相结合，借助满应力准则将应力约束转化为动态尺寸约束，利用单位虚载荷法将位移约束转化为设计变量的显式表达式建立优化模型，然后用数学规划法求解；采用无量纲设计变量实现设计变量连接，对膜结构的厚度进行优化设计；根据对偶理论，应用对偶规划精确映射原问题，再按泰勒展式建立对偶问题的二阶近似。为了提高优化效率，采用射线步调整结构性态，运用粗选有效约束技术筛选约束，并采用主、被动变量循环确保收敛稳定。以MSC／Nastran软件作为结构分析的求解器，以MSC／Patran软件作为开发平台，完成了膜结构截面优化程序。对膜结构的单变位、多变位的结构优化问题进行了优化计算，并与MSC／Nastran优化模块的计算结果进行比较。算例结果表明程序的可靠性、高效性和稳定性以及理论算法的优越性。     

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

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

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

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

Topology optimization design of continuum structures under stress and displacement constraints

ＩｎｔｒｏｄｕｃｔｉｏｎＴｏｐｏｌｏｇｙｏｐｔｉｍｉｚａｔｉｏｎｏｆｃｏｎｔｉｎｕｕｍｓｔｒｕｃｔｕｒｅｓｄｉｄｎｏｔｄｅｖｅｌｏｐｒａｐｉｄｌｙｕｎｔｉｌｒｅｃｅｎｔｔｅｎｙｅａｒｓｏｗｉｎｇｔｏｔｈｅｓｐｅｃｉａｌｄｉｆｆｉｃｕｌｔｉｅｓｉｎｖｏ...     

互逆规划理论及其用于建立结构拓扑优化的合理模型

在数学规划的领域里定义了互逆规划——各自目标函数与约束条件位置相互交换的一对规划. 接着指出,尽管互逆规划与对 偶规划在表面上似乎类似,但是二者存在 5 点不同:(1) 是否为同一个问题的不同;(2) 存在``对偶间隙'与否的不同;(3) 设计变量数目的不同;(4) 是否单目标与多目标问题的不同;(5) 问题合理与否的不同. 然后,基于互逆规划的定义,用以审视结 构拓扑优化模型,给出如下结果:(1) 从这个角度洞悉,在结构优化中,确实有不合理的模型一直被沿用着;(2) 找到了修正不 合理模型使之合理化的方法;(3) 对于给定体积下的柔顺度最小化 (MCVC) 模型,指出了其不合理的原因;(4) MCVC 模型实际是互 逆规划的 m 方,由此建立起其对应的 s 方, 即给出了多个柔顺度约束的体积最小化 (MVCC) 模型;(5)给出了MVCC模型中的结构 柔顺度约束的物理解释和算法,论证了 ICM (independent continuous and mapping) 方法以往关于全局化应力约束的概念和方法;(6)数值算例表明了 MCVC 与 MVCC 模型作为互逆规划的差异,且印证了 MVCC 模型的合理性.MCVC 模型在不同体积约束及多工况下不同的权系数时,得到最优拓扑不同;但 MVCC 模型在多工况柔顺度约束下可得到唯一的最优拓扑.      

平面膜结构拓扑优化的有无复合体方法

将作者对桁架在应力约束下结构拓扑优化的有无复合体模型发展到平面膜结构在应力、位移约束下结构拓扑优化的建模与求解。同时提出了该模型的有效解法,获得了令人满意的数值结果。本文工作表明独立连续拓扑变量的提出对于结构拓扑优化的研究是有意义的。     

考虑性态约束时多工况桁架结构拓扑优化设计

本文提出了一种适用于桁架结构的拓扑优化设计方法。它以杆内力为设计变量,以结构重量为目标函数。该方法的主要特点是:第一,通过引入杆内力为设计变量,既克服了已有方法要求预定位移场这一主要困难,又为在拓扑优化过程中考虑应力、位移等性态约束创造了条件;第二,将多工况的拓扑优化问题描述为一个非光滑的数学规划问题,再通过一个变量代换将其转化为一般的规划为题,进而将原问题的求解又转化为几个线性规划问题的求解;第三,基于结构力学的三个基本方程,将位移与应力约束提成为线性不等式约束,这些约束同重量的目标函数一起构成了拓扑优化设计的线性规划模型。最后,将本方法应用于几个工程算例,得到了满意的数值计算结果。     

A UNIVERSAL APPROACH FOR CONTINUOUS OR DISCRETE NONLINEAR PROGRAMMINGS WITH MULTIPLE VARIABLES AND CONSTRAINTS

IntroductionExtensive research works have been published for solving nonlinear mathematicprogramming problems.Nonetheless,it is still difficult to find an effective and universalapproach for general programming problems with multiple design variables and …     

频率约束的三维连续体结构动力拓扑优化设计

针对频率约束和结构重量最小的动力拓扑优化问题, 基于(independent continuous mapping, ICM)独立、连续、映射方法,建立了频率约束下的三维连续体拓扑优化模型. 利用瑞利商和一阶泰勒展式对频率约束进行了显式化处理,并采用幂函数与复合指数函数作为过滤函数,将优化模型进行了标准化转换, 利用对偶理论及数学规划法进行了求解. 另外,利用质量矩阵和刚度矩阵过滤函数比值与动态约束处理了局部模态和模态交换等数值问题. 最后,通过应用两类不同过滤函数的数值算例表明了文中模型及方法在处理动力拓扑优化问题上的合理性与有效性.