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

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

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

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

连续体结构拓扑优化应力约束凝聚化的ICM方法

为克服应力约束下拓扑优化问题约束数目多、应力敏度计算量大的困难,提出 了应力约束化凝聚化的ICM方法. 在利用Mises强度理论将应力约束转换成应变能约束后, 提出了应力约束凝聚化的两条途径：其一为应力全局化的方法,其二为应力约束集成化的方 法. 由此建立了多工况下以重量为目标、以凝聚化应变能为约束的连续体结构优化模型,并 利用对偶理论对优化模型进行了求解. 4个数值算例表明：该方法具有较高的计算效率,得 到的拓扑结构比较合理,不仅适用于二维连续体结构,也适用于三维连续体结构.     

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

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

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

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

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

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

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

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

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

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

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

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

频率约束板结构拓扑优化

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

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

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

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

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

The ICM method with objective function transformed by variable discrete condition for continuum structure

ICM (Independent Continuous Mapping) method can solve topological optimization problems with the minimized weight as the objective and subjected to displacement constraints. To get a clearer topological configuration, by introducing the discrete condition of topological variables and integrating with the original objective, an optimal model with multi-objectives is formulated to make the topological variables approach 0 or 1 as near as possible, and the model reduces the effect of deleting rate on the result. The image-filtering method is employed to eliminate the checkerboard patterns and mesh dependence that occurred in the topology optimization of a continuum structure. The computational efficiency is enhanced through selecting quasi-active displacement constraints and a design region. Numerical examples indicate that this algorithm is robust and practicable, though the number of iterations is slightly increased with respect to the original algorithm. The project supported by the National Natural Science Foundation of China (10472003), Beijing Natural Science (3002002) and Beijing Educational Committee Foundations (KM200410005019) and Suspensoried by American MSC Company. The English text was polished by Keren Wang.     

响应面方法在二维连续体形状优化中的应用

简单介绍了响应面方法的实质,并将其应用于二维连续体结构形状优化中,结合有限元方法,近似拟合出了状态约束(应力和位移)对设计变量的显函数,代替了复杂的解析灵敏度分析,以MSC/Patran为平台采用PCL语言开发了二维连续体形状优化的程序．文中算例通过与差分法的比较表明了这一方法的有效性和可靠性．     

基于概率的平面连续体结构拓扑优化

构建了随机载荷作用下具有随机参数的平面连续体结构拓扑优化设计模型,以结构总质量均值极小化为目标函数,以满足单元应力可靠性为约束条件;对应力可靠性约束进行了等价显示化处理;导出了随机参数结构的应力反应数字特征的计算表达式;提出了相应的求解策略和方法。算例的优化结果表明文中模型的合理性和方法的有效性。     

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

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