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

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

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

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

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

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

Independent continuous mapping for topological optimization of frame structures

Based on the Independent Continuous Mapping method (ICM), a topological optimization model with continuous topological variables is built by introducing three filter functions for element weight, element allowable stress and element stiffness, which transform the 0-1 type discrete topological variables into continuous topological variables between 0 and 1. Two methods for the filter functions are adopted to avoid the structural singularity and recover falsely deleted elements: the weak material element method and the tiny section element method. Three criteria (no structural singularity, no violated constraints and no change of structural weight) are introduced to judge iteration convergence. These criteria allow finding an appropriate threshold by adjusting a discount factor in the iteration procedure. To improve the efficiency, the original optimization model is transformed into a dual problem according to the dual theory and solved in its dual space. By using MSC/Nastran as the structural solver and MSC/Patran as the developing platform, a topological optimization software of frame structures is accomplished. Numerical examples show that the ICM method is very efficient for the topological optimization of frame structures.The project supported by the National Natural Science Foundation of China (10472003) and Beijing Natural Science Foundation (3042002) The English text was polished by Yunming Chen.     

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

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

Topology optimization design of continuum structures under stress and displacement constraints

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

膜结构截面的离散变量优化方法

通过将截面离散变量映射为连续变量的途径,建立了力学模型,在连续变量截面优化最优解附近构造两节模型,并采用两节无限小单元的无穷组合的方法和变量无量纲化技术对膜结构截面离散变量进行了优化。建立了多工况下具有尺寸和应力、位移约束的优化模型。同时,在MSC.Patran和MSC.Nastran上进行了二次开发。数值算例结果表明该理论的高效性和二次开发的可行性。     

对偶法求解应力与位移约束不同显式的优化模型

基于响应面法,通过内力将应力约束显式化为正变量的一阶近似,位移约束采用倒变量一阶近似,建立以桁架结构重量最小为目标的截面优化模型,采用对偶规划理论转化为拟无约束模型,采用序列二次规划算法求解。此方法同应力约束零阶近似和一阶近似(使用响应面方法显式化)等方法进行了算例比较,表明本文的优化模型有较好的求解效率。     

频率约束板结构拓扑优化

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

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

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

非光滑优化的双层规划模型及内点算法

对于包含接触约束的非光滑结构优化问题,其非光滑性体现在状态函数并不是处处可微的,针对含有应力约束及接触约束的非光滑结构优化问题,建立了一种双层规划模型,避免了求解时非光滑性所带来的问题,同时提出了一种迭代算法,用对偶内点二次规划进行分析,线性规划进行优化,算例表明这种方法十分有效。     

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

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

应力约束下三维连续体结构拓扑优化分析

本文基于ICM(独立、连续、映射)方法,建立了以连续体结构重量为目标的应力约束下多工况的三维连续体结构拓扑优化模型。利用von Mises强度理论,提出了应力约束全局化方法,从而将局部的应力约束转化为全局的应变能约束,由此减少了约束数目,降低了问题的规模,避免了敏度分析。另外,利用对偶理论,将建立的优化模型转化为对偶模型,用序列二次规划法进行了求解,从而减少了设计变量的数目,提高了求解效率。借助MSC．Patran和MSC．Nastran提供的软件平台,利用PCL语言实现了本文的优化算法。数值算例验证了方法的可行性与有效性。     

应力约束全局化策略下的连续体结构拓扑优化

利用Mises强度理论，提出了应力约束全局化策略，将局 部的应力约束问题转化为结构整体的应变能约束问题. 基于ICM(独立、连续、映射)方法， 引入了独立、连续的拓扑变量，对单元重量、单元刚度和单元许用应力的过滤函数进行了选 择，建立了以重量为目标，以结构应变能代替应力约束的多工况下连续体结构拓扑优化模型， 寻找到了多工况下的最佳传力路径. 运用对偶二次规划方法对上述优化模型进行了求解. 另 外，利用PCL语言，在MSC/PATRAN的开发平台上，实现了应用应力约束全局化策略进行连 续体结构拓扑优化的模块化处理. 数值算例表明了该方法的可行性和有效性.     

基于ICM方法的刚架拓扑优化

将ICM(独立、连续、映射)优化方法应用于刚架结构，建立了以重量为目标函数、考虑多工况下应力和位移约束的刚架结构拓扑优化模型。借助于磨光函数和过滤函数，实现了离散变量和连续变量之间的相互转换。考虑多工况应力拓扑解之间、应力和位移约束下拓扑解之间的协调，得到了合理的结果。文中给出的刚架结构拓扑优化算例表明，刚架结构的拓扑结构不同于桁架或连续体。     

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

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

基于层次分解方法的桁架结构形状优化

对于桁架结构形状优化,应用层次分解优化方法,将设计变量分成杆件截面积和节点位置两类变量。求解时分为两层,第一层在给定节点位置下对杆件截面进行优化,同时考虑了应力、局部稳定约束和位移约束的重量最轻;第二层假定截面层的有效位移约束作用不变,求解一个使桁架刚度增强的二次规划问题,获得既不违反约束,又使目标函数不上升的新的节点位置,再返回第一层。两层交替进行直至收敛。     

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

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

APPROACH FOR LAYOUT OPTIMIZATION OF TRUSS STRUCTURES WITH DISCRETE VARIABLES UNDER DYNAMIC STRESS, DISPLACEMENT AND STABILITY CONSTRAINTS

A mathematical model was developed for layout optimization of truss structures with discrete variables subjected to dynamic stress, dynamic displacement and dynamic stability constraints. By using the quasi-static method, the mathematical model of structure optimization under dynamic stress, dynamic displacement and dynamic stability constraints were transformed into one subjected to static stress, displacement and stability constraints. The optimization procedures include two levels, i.e., the topology optimization and the shape optimization. In each level, the comprehensive algorithm was used and the relative difference quotients of two kinds of variables were used to search the optimum solution. A comparison between the optimum results of model with stability constraints and the optimum results of model without stability constraint was given. And that shows the stability constraints have a great effect on the optimum solutions.     

Plate/shell topological optimization subjected to linear buckling constraints by adopting composite exponential filtering function

In this paper, a model of topology optimization with linear buckling constraints is established based on an independent and continuous mapping method to minimize the plate/shell structure weight. A composite exponential function (CEF) is selected as filtering functions for element weight, the element stiffness matrix and the element geomet-ric stiffness matrix, which recognize the design variables, and to implement the changing process of design variables from“discrete”to“continuous”and back to“discrete”. The buck-ling constraints are approximated as explicit formulations based on the Taylor expansion and the filtering function. The optimization model is transformed to dual programming and solved by the dual sequence quadratic programming algo-rithm. Finally, three numerical examples with power function and CEF as filter function are analyzed and discussed to demonstrate the feasibility and efficiency of the proposed method.