一种考虑破损-安全的连续体频率约束拓扑优化设计方法技术

本发明专利技术公开一种考虑破损‑安全的连续体频率约束拓扑优化设计方法，首先基于连续体基结构建立有限元模型；确定破损工况数，并输入初始破损区域中心坐标值和破损区域参数；输入频率约束优化参数，形成考虑破损‑安全条件下的频率约束优化模型；依次对包含破损工况的破损‑安全结构进行模态分析，并提取各结构破损状况的模态分析结果；建立考虑破损‑安全条件下频率约束的拓扑优化模型的近似连续数学优化列式；采用数学规划算法，对优化模型进行求解；对最优拓扑构型进行反演处理，获得最优拓扑构型以及最优构型的基频。本发明专利技术可以明显提高结构的冗余度，为动力学领域考虑破损‑安全的拓扑优化问题提供一个参考。

A topology optimization design method of continuum with frequency constraint considering damage safety

【技术实现步骤摘要】
一种考虑破损-安全的连续体频率约束拓扑优化设计方法
本专利技术属于工程结构设计
，尤其涉及考虑破损-安全条件下的连续体频率约束拓扑优化设计方法。
技术介绍
拓扑优化是在给定的外荷载和边界条件下，以材料分布为优化对象，寻求结构的最佳传力路径，该优化方法是伴随着有限元理论和计算机软件技术的高速发展逐渐发展起来的。与尺寸优化和形状优化相比，结构拓扑优化能确定的设计参数更多，能更显著的节省材料，取得更大的经济效益，同时能在结构设计的初始阶段为工程设计人员提供一个综合性的概念性设计，把最优方案应用在结构的布局上，对工程设计人员具有很大的吸引力。经过拓扑优化所得的传统拓扑结构通常是在外力作用下的最佳传力路径，材料的性能得到充分发挥，但是该最优构型倾向于静定结构，所带来的弊端就是结构缺少冗余，在这种情况下结构的局部破损很容易导致整体的破坏。在安全性要求较高如航天航空、船舶、高速列车、桥梁和超高建筑等领域，拓扑结构对局部破损的高度敏感性会使其经济效果有余而安全保障不足，在安全性较高的工程领域常被认为经济但缺乏安全保障，因此考虑破损-安全的扑优化受到重视。工程结构普遍需要考虑振动问题，结构在动力载荷下的响应同结构的固有频率及外载频率有关，极端情况下，当外载频率与结构的某阶固有频率非常接近时，结构会发生共振现象。因此在进行结构动力设计时，需要对结构的固有频率进行限制，例如火箭、导弹这类飞行器的设计对结构的固有频率就有严格的要求，所以对结构频率进行拓扑优化设计具有重要意义。本专利技术针对传统的连续体结构频率约束拓扑优化设计问题，考虑了破损-安全的理念设计，提出了考虑破损-安全条件下的连续体频率约束拓扑优化设计方法。传统频率拓扑优化的结果类似于骨架结构，这种拓扑构型虽然在考虑频率的条件下使材料的利用率最大化，但在面临局部破损(碰撞、爆炸、疲劳失效引起的)，传统拓扑结构会发生完全破坏，造成灾难性的后果。引入破损-安全的是设计理念后，这种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法，结合结构模态有限元分析技术，在频率约束拓扑优化中考虑结构内部存在破损的条件，能够有效解决传统频率优化结构部件过于“高效”的问题。这种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法不仅可以实现传统频率拓扑优化的设计目的，更重要的是能明显提高最优拓扑构型的冗余度，在最优结构任意位置发生破损后，结构的频率会在约束频率附近，结构仍能够正常工作。本专利技术有助于降低最优结构频率对局部破损的敏感性，提高最优结构的冗余度，这对于航空、航天、汽车、桥梁、土木等许多领域的结构动力优化方面具有重要的理论意义和广泛的工程应用价值。
技术实现思路
本专利技术针对传统的连续体结构频率约束拓扑优化设计问题，考虑了破损-安全的设计理念，采用了层级铺设的策略来近似随机、无限的破损工况，首次将破损-安全的设计理念与频率约束拓扑优化结合起来。本专利技术以频率这种结构动力特性为约束对象，这种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法可以明显提高最优结构的冗余度，显著降低最优结构频率对局部破损的敏感性，在结构动力拓扑优化领域具有极强的理论意义和应用价值。为实现上述目的，本专利技术采用如下技术方案：一种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法包括以下步骤：第一步，确定连续体基结构，并建立有限元模型；第二步，确定破损工况数，并在有限元模型中输入初始破损区域中心坐标值和破损区域参数；第三步，输入频率约束优化参数，形成考虑破损-安全条件下的频率约束优化模型；第四步，依次对包含破损工况的破损-安全结构进行模态分析，并提取各个破损工况的模态分析结果；第五步，建立考虑破损-安全条件下频率约束的拓扑优化模型的近似连续数学优化列式；第六步，采用数学规划算法，对优化模型进行求解；第七步，对最优拓扑构型进行反演处理，获得最优拓扑构型以及最优构型的基频。本专利技术相比现有技术的优点在于：相比于传统的结构频率拓扑优化设计，充分考虑了结构破损-安全的设计理念，在优化设计中添加破损工况，可以明显提高最优结构的冗余度，显著降低最优结构对局部破损的敏感性，在结构动力拓扑优化领域具有极强的理论意义和应用价值。附图说明图1是考虑破损-安全条件下的连续体频率约束拓扑优化设计流程图。图2是考虑破损-安全条件的基结构及最优拓扑图。图3是考虑破损-安全条件的拓扑优化过程图。图4为迭代曲线图，(a)为结构总体积迭代曲线，(b)为各结构破损状况的基频迭代曲线。具体实施方式如图1所示，本专利技术中提供一种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法，具体解决方案如下：第一步，确定连续体基结构，并建立有限元模型：首先基于设计需求，确定所要优化连续体结构的最大设计边界，即确定基结构，并依据某些性能需求，确定破损设计域及非破损设计域；其次，基于MSC.Patran软件平台，建立基结构几何模型，并划分有限元网格；最后定义材料，赋予单元属性，施加边界条件；最后，设定模态分析输出选项。第二步，确定破损工况数，并输入初始破损区域中心坐标值和破损区域参数：在基于MSC.Patran软件内置语言PCL(PatranCommandLanguage)编写的主程序内，首先设置破损工况数，破损工况数可以根据实际问题确定；然后输入初始破损区域中心二维坐标值；然后输入破损区域参数：损伤区对角线一半的平方、损伤区边长一半的平方，这样正方形的初始破损区域就可以确定下来。第三步，输入频率约束优化参数，形成考虑破损-安全条件下的频率约束优化模型：依据结构固有频率性能设计需求，输入约束参数，包括约束基频大小、收敛精度以及过滤半径；这样破损区域大小、形状以及频率约束就确定下来，可以建立考虑破损-安全的以基频为约束、结构重量最小为优化目标的拓扑优化模型：式中：t表示拓扑变量向量，ti为第i个单元的拓扑变量，表征单元的有或无；t∈EN表示拓扑变量t是属于n维欧式空间的向量；W表示结构中重量，为第i个单元的初始重量；fw(ti)、fk(ti)、fm(ti)分别是ICM(Independent,Continuous,Mapping)方法中提出的重量过滤函数、刚度矩阵过滤函数以及质量矩阵过滤函数；λl表示第l个破损工况的基频对应的特征值，表示第l个破损工况的基频对应的特征值约束值；L表示的是破损工况的个数；N表示的是单元总数；第四步，依次对包含破损工况的破损-安全结构进行模态分析，并提取各个破损工况的模态分析结果：通常称含有一个结构局部破损区域的基结构为一个结构破损状况，可以在PCL主程序中实现从初始破损区域开始，根据破损工况数依次对结构进行破损，并对每一个结构破损状况进行模态分析。依次提取每一个结构破损状况单元的模态分析结果，分别是Uil和Vil，表示第l个破损状况下单元i的基频模态应变能，表示第l个破损状况下单元i的基频模态动能。这样，通过提取模本文档来自技高网...

【技术保护点】
1.一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，包括以下步骤：/n第一步，确定连续体结构的破损设计域和非破损设计域，并建立相应的有限元模型，设定边界条件，材料属性参数；/n第二步，输入破损工况数，并向有限元模型中输入初始破损区域中心坐标值，同时输入破损区域参数，这样就确定了初始破损区域的位置和形状；/n第三步，输入频率约束优化参数，包括约束基频的大小、收敛精度、过滤半径，形成考虑破损-安全条件下的频率约束优化模型；/n第四步，含有一个结构局部破损区域的基结构称为一个结构破损状况；从初始破损区域开始，根据破损工况数依次对结构进行破损，并对每一个结构破损状况进行模态分析，并提取每一个结构破损状况单元的模态分析结果，为建立显式优化方程提供结构单元的力学性能参数；/n第五步，形成考虑破损-安全的以频率为约束、结构重量最小为目标的近似连续数学优化列式；/n第六步，采用数学规划算法，对连续数学优化列式进行求解；/n第七步，对拓扑变量进行反演，获得最优拓扑构型及对应的基频。/n

【技术特征摘要】
1.一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，包括以下步骤：
第一步，确定连续体结构的破损设计域和非破损设计域，并建立相应的有限元模型，设定边界条件，材料属性参数；
第二步，输入破损工况数，并向有限元模型中输入初始破损区域中心坐标值，同时输入破损区域参数，这样就确定了初始破损区域的位置和形状；
第三步，输入频率约束优化参数，包括约束基频的大小、收敛精度、过滤半径，形成考虑破损-安全条件下的频率约束优化模型；
第四步，含有一个结构局部破损区域的基结构称为一个结构破损状况；从初始破损区域开始，根据破损工况数依次对结构进行破损，并对每一个结构破损状况进行模态分析，并提取每一个结构破损状况单元的模态分析结果，为建立显式优化方程提供结构单元的力学性能参数；
第五步，形成考虑破损-安全的以频率为约束、结构重量最小为目标的近似连续数学优化列式；
第六步，采用数学规划算法，对连续数学优化列式进行求解；
第七步，对拓扑变量进行反演，获得最优拓扑构型及对应的基频。


2.根据权利要求1所述的一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，第一步具体为：
首先基于设计需求，确定所要优化连续体结构的最大设计边界，即确定基结构，并依据性能需求，确定破损设计域及非破损设计域；其次，基于MSC.Patran软件平台，建立基结构几何模型，并划分有限元网格；最后定义材料，赋予单元属性，施加边界条件；最后，设定模态分析输出选项。


3.根据权利要求2所述的一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，第二步具体为：
在基于MSC.Patran软件内置语言PCL编写的主程序内，首先设置破损工况数，破损工况数根据实际问题确定；然后输入初始破损区域中心二维坐标值；然后输入破损区域参数：损伤区对角线一半的平方、损伤区边长一半的平方，确定正方形的初始破损区域。


4.根据权利要求3所述的一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，第三步具体为：
依据结构固有频率性能设计需求，输入约束参数，包括约束基频大小、收敛...

【专利技术属性】
技术研发人员：杜家政，孟凡伟，汤威，张颖，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11