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

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

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

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

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

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

多工况应力和位移约束下连续体结构拓扑优化

将文「１」所提出的对拓扑变量的独立连续映射（ＩＧＭ）的拓扑估化方法应用于连续体结构,从而建立了统一的以重量为目标,考虑应力和位移约束的连续体结构拓扑优化模型。通过对位移－应力拓扑解和各工况下应力拓扑解的综合协调,进而对于协调拓扑解按照阈值完成从离散到连续的反演,并且采用分层与加权策略克服了“荷载病态”困难。给出的经典的二维平面问题和三维连续体结构拓扑优化算例表明,这种统一的模型由骨架结构发展到连续     

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

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

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

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

基于可靠性的桁架结构拓扑优化设计

建立了以杆截面为设计变量、结构重量极小化为目标、具有位移、应力等性态可靠性约束的桁架结构拓扑优化设计数学模型．通过引入可靠性安全系数，并利用结构力学的三个基本方程，将结构的位移和杆件应力可靠性约束等价显示化为设计变量的线性函数，使原基于可靠性的优化模型转化为常规的序列线性规划问题，利用修正的单纯形法求解．算例表明文中提出的方法既简单又有效．     

连续体结构考虑离散性目标的ICM方法

为了使求解位移和应力约束下重量最轻的连续体结构拓扑优化问题的ICM(独立连续映射)方法得到更清晰的拓扑形式,引入拓扑变量离散性条件作为目标之一,与原目标组成多目标规划模型,使拓扑变量向0或1两端靠近,减小了删除率对结果的影响;用图形处理方法消除了棋盘格现象及网格依赖性;通过初选准有效约束及选择设计区域等,提高了求解效率。算例表明:改进后的算法虽然迭代次数比原算法略多,但更稳键,更实用。     

基于相场描述的拉压不对称强度结构拓扑优化

运用了基于相场描述的拓扑优化方法，来寻找在拉伸和压缩中表现出不对称强度行为的连续体结构的最优布局。依据Drucker-Prager屈服准则和幂率插值方案，优化问题可以描述为在局部应力约束下的最小化结构的体积。用qp放松法来解决应力约束的奇异性，并采用基于P-norm函数的聚合方法对应力约束进行凝聚，该方法实现了约束个数的降低，同时引入了稳定转化法来处理大量的局部应力约束和高度非线性的应力行为，以修正应力，提高优化收敛的稳定性。在优化问题求解时，使用拉格朗日乘子法对目标函数和应力约束进行处理。利用伴随变量法进行灵敏度分析，并通过求解Allen-Cahn方程更新相场函数设计变量。数值算例证明了该优化模型和相应数值技术的有效性，相关算例还揭示了考虑拉压不同强度和考虑同拉压强度约束时得到的结构优化拓扑构型具有显著的差异。     

A Method for Topological Optimization of Structures with Discrete Variables under Dynamic Stress and Displacement Constraints

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂａｓｉｃｍｅｔｈｏｄｆｏｒｔｈｅｓｔｒｕｃｔｕｒｅｔｏｐｏｌｏｇｙｏｐｔｉｍｉｚａｔｉｏｎｉｓｂａｓｅｄｏｎｔｈｅｇｒｏｕｎｄｓｔｒｕｃｔｕｒｅｍｅｔｈｏｄ .Ａｓｅｔｏｆｎｏｄｅｓｉｓｃｏｎｓｔｒｕｃｔｅｄａｃｃｏｒｄｉｎｇｔｏｇｉｖｅｎｓｕｐｐｏｒｔｉｎｇｃｏｎｄｉｔｉｏｎｓ,ｌｏａｄｃａｓｅｓａｎｄｏｔｈｅｒｃｏｎｄｉｔｉｏｎｓ.Ｔｈｅｎｏｄｅｓａｒｅｌｉｎｋｅｄｂｙｍｅｍｂｅｒｓｔｏｆｏｒｍｔｈｅｉｎｉｔｉａｌｇｒｏｕｎｄｓｔｒｕｃｔｕｒｅ .Ａｌｌｔｏｐｏｌｏ…     

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

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

基于ICM方法三维连续体结构拓扑优化

基于ICM方法,建立了在应力和位移约束下以重量为目标的多工况下的三维连续体结构拓扑优化模型.利用von Mises强度理论,提出了应力全局化的方法,从而将局部的应力约束转化为全局的应变能约束问题,减少了约束数目,避免了敏度分析的困难;利用单位虚载荷法,将位移约束表示为设计变量的显式关系.为减小由于不同物理量在数量级上相差太大引起数值计算的误差,将应变能约束和位移约束进行无量纲化,由此建立了包含两类约束的无量纲化的优化模型.同时处理了多工况下的最佳传力路径的问题.利用对偶规划理论对模型进行了求解.另外,利用PCL语言在MSC/Patran的开发平台上实现了该文算法.数值算例表明了该方法的可行性和有效性.     

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

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

统一骨架与连续体的结构拓扑优化的ICM理论与方法

技术了基于ＩＣＭ方法的结构拓扑优化新模型并应用于骨架与连续体结构。ＩＣＭ方法意指独立、连续变量与映射及其反演。新模型将两种结构统一建立了具有重量目标函数和多工况下应力与位移约束下的优化问题,提出的过滤函数是ＩＣＭ方法的关键技术之一。说明了优化策略与算法。     

强迫谐振动下连续体结构拓扑优化

应用结构拓扑优化ICM(独立连续映射)方法,对强迫谐振动下结构拓扑优化问题建立了以重量极小为目标,位移幅值为约束的优化模型.位移幅值采用一阶泰勒展式近似,由于拓扑优化中设计变量数目通常很多,对强迫谐振动位移幅值的敏度分析推导了伴随法公式,使得一次敏度分析可以计算出对所有设计变量的偏导数,克服了采用直接法敏度分析中一次只能计算出对一个设计变量的偏导数的不足.算例表明用伴随法分析敏度在结构拓扑优化中可以大幅提高计算效率,ICM方法采用独立于截面及形状参数的拓扑优化设计变量更清晰地反映了拓扑优化的本质.     

METHOD BASED ON DUAL-QUADRATIC PROGRAMMING FOR FRAME STRUCTURAL OPTIMIZATION WITH LARGE SCALE

The optimality criteria (OC) method and mathematical programming (MP) were combined to found the sectional optimization model of frame structures. Different methods were adopted to deal with the different constraints. The stress constraints as local constraints were approached by zero-order approximation and transformed into movable sectional lower limits with the full stress criterion. The displacement constraints as global constraints were transformed into explicit expressions with the unit virtual load method. Thus an approximate explicit model for the sectional optimization of frame structures was built with stress and displacement constraints. To improve the resolution efficiency, the dual-quadratic programming was adopted to transform the original optimization model into a dual problem according to the dual theory and solved iteratively in its dual space. A method called approximate scaling step was adopted to reduce computations and smooth the iterative process. Negative constraints were deleted to reduce the size of the optimization model. With MSC/Nastran software as structural solver and MSC/Patran software as developing platform, the sectional optimization software of frame structures was accomplished, considering stress and displacement constraints. The examples show that the efficiency and accuracy are improved.     

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

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

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

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

基于ICM方法的多工况位移约束下板壳结构拓扑优化设计

基于ICM方法,建立了多种载荷工况下受位移约束极小化结构重量的拓扑优化近似显式模型,采用精确对偶映射下的序列二次规划算法进行求解,得到了结构的最优拓扑.数值算例表明:ICM方法也可以很好地应用于多工况下板壳结构拓扑优化设计,且收敛快,稳定性好；边界条件和位移约束值都会影响结构的最优拓扑；多工况最优传力路径不是各个单工况最优传力路径的简单叠加.     

SOME RESULTS ON OPTIMAL TOPOLOGY DESIGN OF TRUSS WITH UNILATERAL CONTACT CONSTRAINTS

The optimal topology design of truss structures concerning stress and frictionless unilateral contact displacement constraints is investigated. The existence of ununique optimal solution under contact gaps is found. This shows that the contact conditions have an effect on structural topology, and different initial contact gaps may lead to different structural topologies. To avoid the singular optima in structural topology optimization in multiple loading cases, an ε-relaxed method is adopted to establish the relaxing topology optimization formulations. The problem is solved by means of a two-level optimization method. In the first sublevel, the solution of the frictionless unilateral contact problem is obtained by solving an equivalent quadratic programming. In the second sublevel, topology optimization of truss is carried out by an ε-relaxed method. The validity of the method proposed is verified by computational results. Supported by the Postdoctoral Science and Research Foundation of Shanghai, China.