本发明属于工程结构设计技术领域，尤其涉及考虑破损-安全条件下的连续体频率约束拓扑优化设计方法。
背景技术：
拓扑优化是在给定的外荷载和边界条件下，以材料分布为优化对象，寻求结构的最佳传力路径，该优化方法是伴随着有限元理论和计算机软件技术的高速发展逐渐发展起来的。与尺寸优化和形状优化相比，结构拓扑优化能确定的设计参数更多，能更显著的节省材料，取得更大的经济效益，同时能在结构设计的初始阶段为工程设计人员提供一个综合性的概念性设计，把最优方案应用在结构的布局上，对工程设计人员具有很大的吸引力。
经过拓扑优化所得的传统拓扑结构通常是在外力作用下的最佳传力路径，材料的性能得到充分发挥，但是该最优构型倾向于静定结构，所带来的弊端就是结构缺少冗余，在这种情况下结构的局部破损很容易导致整体的破坏。在安全性要求较高如航天航空、船舶、高速列车、桥梁和超高建筑等领域，拓扑结构对局部破损的高度敏感性会使其经济效果有余而安全保障不足，在安全性较高的工程领域常被认为经济但缺乏安全保障，因此考虑破损-安全的扑优化受到重视。
工程结构普遍需要考虑振动问题，结构在动力载荷下的响应同结构的固有频率及外载频率有关，极端情况下，当外载频率与结构的某阶固有频率非常接近时，结构会发生共振现象。因此在进行结构动力设计时，需要对结构的固有频率进行限制，例如火箭、导弹这类飞行器的设计对结构的固有频率就有严格的要求，所以对结构频率进行拓扑优化设计具有重要意义。
本发明针对传统的连续体结构频率约束拓扑优化设计问题，考虑了破损-安全的理念设计，提出了考虑破损-安全条件下的连续体频率约束拓扑优化设计方法。传统频率拓扑优化的结果类似于骨架结构，这种拓扑构型虽然在考虑频率的条件下使材料的利用率最大化，但在面临局部破损(碰撞、爆炸、疲劳失效引起的)，传统拓扑结构会发生完全破坏，造成灾难性的后果。引入破损-安全的是设计理念后，这种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法，结合结构模态有限元分析技术，在频率约束拓扑优化中考虑结构内部存在破损的条件，能够有效解决传统频率优化结构部件过于“高效”的问题。这种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法不仅可以实现传统频率拓扑优化的设计目的，更重要的是能明显提高最优拓扑构型的冗余度，在最优结构任意位置发生破损后，结构的频率会在约束频率附近，结构仍能够正常工作。本发明有助于降低最优结构频率对局部破损的敏感性，提高最优结构的冗余度，这对于航空、航天、汽车、桥梁、土木等许多领域的结构动力优化方面具有重要的理论意义和广泛的工程应用价值。
技术实现要素：
本发明针对传统的连续体结构频率约束拓扑优化设计问题，考虑了破损-安全的设计理念，采用了层级铺设的策略来近似随机、无限的破损工况，首次将破损-安全的设计理念与频率约束拓扑优化结合起来。本发明以频率这种结构动力特性为约束对象，这种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法可以明显提高最优结构的冗余度，显著降低最优结构频率对局部破损的敏感性，在结构动力拓扑优化领域具有极强的理论意义和应用价值。为实现上述目的，本发明采用如下技术方案：
一种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法包括以下步骤：
第一步，确定连续体基结构，并建立有限元模型；
第二步，确定破损工况数，并在有限元模型中输入初始破损区域中心坐标值和破损区域参数；
第三步，输入频率约束优化参数，形成考虑破损-安全条件下的频率约束优化模型；
第四步，依次对包含破损工况的破损-安全结构进行模态分析，并提取各个破损工况的模态分析结果；
第五步，建立考虑破损-安全条件下频率约束的拓扑优化模型的近似连续数学优化列式；
第六步，采用数学规划算法，对优化模型进行求解；
第七步，对最优拓扑构型进行反演处理，获得最优拓扑构型以及最优构型的基频。
本发明相比现有技术的优点在于：
相比于传统的结构频率拓扑优化设计，充分考虑了结构破损-安全的设计理念，在优化设计中添加破损工况，可以明显提高最优结构的冗余度，显著降低最优结构对局部破损的敏感性，在结构动力拓扑优化领域具有极强的理论意义和应用价值。
附图说明
图1是考虑破损-安全条件下的连续体频率约束拓扑优化设计流程图。
图2是考虑破损-安全条件的基结构及最优拓扑图。
图3是考虑破损-安全条件的拓扑优化过程图。
图4为迭代曲线图，(a)为结构总体积迭代曲线，(b)为各结构破损状况的基频迭代曲线。
具体实施方式
如图1所示，本发明中提供一种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法，具体解决方案如下：
第一步，确定连续体基结构，并建立有限元模型：
首先基于设计需求，确定所要优化连续体结构的最大设计边界，即确定基结构，并依据某些性能需求，确定破损设计域及非破损设计域；其次，基于msc.patran软件平台，建立基结构几何模型，并划分有限元网格；最后定义材料，赋予单元属性，施加边界条件；最后，设定模态分析输出选项。
第二步，确定破损工况数，并输入初始破损区域中心坐标值和破损区域参数：
在基于msc.patran软件内置语言pcl(patrancommandlanguage)编写的主程序内，首先设置破损工况数，破损工况数可以根据实际问题确定；然后输入初始破损区域中心二维坐标值；然后输入破损区域参数：损伤区对角线一半的平方、损伤区边长一半的平方，这样正方形的初始破损区域就可以确定下来。
第三步，输入频率约束优化参数，形成考虑破损-安全条件下的频率约束优化模型：
依据结构固有频率性能设计需求，输入约束参数，包括约束基频大小、收敛精度以及过滤半径；这样破损区域大小、形状以及频率约束就确定下来，可以建立考虑破损-安全的以基频为约束、结构重量最小为优化目标的拓扑优化模型：
式中：t表示拓扑变量向量，ti为第i个单元的拓扑变量，表征单元的有或无；t∈en表示拓扑变量t是属于n维欧式空间的向量；w表示结构中重量，
第四步，依次对包含破损工况的破损-安全结构进行模态分析，并提取各个破损工况的模态分析结果：
通常称含有一个结构局部破损区域的基结构为一个结构破损状况，可以在pcl主程序中实现从初始破损区域开始，根据破损工况数依次对结构进行破损，并对每一个结构破损状况进行模态分析。依次提取每一个结构破损状况单元的模态分析结果，分别是uil和vil，
第五步，建立考虑破损-安全条件下频率约束的拓扑优化模型的近似连续数学优化列式：
将设计变量取为
其中α＝γw/γk，γw＝2是重量过滤函数的幂指数，γk＝5是单元刚度过滤函数的幂指数；dil和el是显式约束不等式的常系数；ν是迭代次数；
为使拓扑变量t尽量取0或1，加入离散性条件并与原目标按线性加权组成单目标，对目标二次近似后略去常数项得：
这样更新近似连续数学优化列式为：
第六步，采用数学规划算法，对优化模型进行求解：
考虑到优化模型(4)中的设计变量数目远大于约束数目，因此，基于库恩-塔克条件将其转化成对偶模型，之后采用序列二次规划算法求解。如果得到的解，满足收敛精度条件，则输出最优解，进行第七步；如果不满足收敛精度条件，则修改拓扑变量，返回第四步，进行下一轮的破损状况有限元分析及优化模型处理，直到满足收敛精度条件。
第七步，对最优拓扑构型进行反演处理，获得最优拓扑构型以及最优构型的基频：
根据拓扑优化结果的拓扑变量云图，选择合适的阈值对拓扑变量进行反演，实现拓扑优化变量从连续到离散0/1的映射，获得最优拓扑构型以及对应的基频。
本发明公开一种考虑破损-安全条件下的连续体频率约束拓扑优化设计方法，包括以下步骤：(1)确定连续体基结构，并建立有限元模型；(2)确定破损工况数，并输入初始破损区域中心坐标值和破损区域参数；(3)输入频率约束优化参数，形成考虑破损-安全条件下的频率约束优化模型；(4)依次对包含破损工况的破损-安全结构进行模态分析，并提取各个破损工况的模态分析结果；(5)建立考虑破损-安全条件下频率约束的拓扑优化模型的近似连续数学优化列式；(6)采用数学规划算法，对优化模型进行求解；(7)对最优拓扑构型进行反演处理，获得最优拓扑构型以及最优构型的基频。本发明可以有效解决考虑破损-安全条件下的连续体频率约束的拓扑优化设计问题，相比较于传统的频率拓扑优化，可以明显提高结构的冗余度，为动力学领域考虑破损-安全的拓扑优化问题提供一个参考。
实施算例：
下面结合一个考虑破损-安全条件的频率拓扑优化算例对本发明的具体实施步骤进行详细说明。
第一步，首先基于msc.patran建立尺寸如图2的(a)所示的基结构，基结构为200mm×100mm×9mm的矩形板，一端固定，另一端中间位置布置一个集中质量块，质量块大小是0.001g；材料弹性模量为68890mpa，泊松比为0.3，密度为1.0×10-9mg/mm3；划分为20×40个矩形单元，左边界全部采用固定约束，同时约束所有单元x、y方向的转动和z方向的平动；基结构的第一阶频率是1958hz，频率约束为第一阶频率大于900hz；有限元模型如图2的(b)所示。考虑8个破损工况，其破损区域为边长50mm×50mm的正方形均匀的分布在基结构中，如图2的(c)所示。最后在模态分析设置选项中选择输出单元应变能；
第二步，在基于pcl语言编写的主程序内，设置破损工况参数。首先设置破损工况数，令破损工况数等于8；然后输入初始破损区域中心二维坐标值(25,25)；然后输入破损区域参数：损伤区对角线一半的平方是225、损伤区边长一半的平方是225，这样正方形的初始破损区域就可以确定下来；
第三步，在优化主程序中输入频率约束优化参数。频率优化参数：频率约束方向≥，约束值是900hz；收敛精度是0.001，过滤半径是7.5。这样破损区域大小、形状以及频率约束就确定下来；
第四步，采用msc.nastran对初始结构破损状况进行模态分析，并通过for循环依次对结构破损状况进行模态分析，并通过程序实现单元模态动能和模态应变能的自动提取，为建立显式优化方程提供结构单元的力学性能参数；
第五步，基于灵敏度分析及泰勒展式等方法，建立考虑破损-安全条件下频率约束的拓扑优化模型的近似连续数学优化列式；
第六步，采用数学规划算法，对优化模型进行求解。如果满足给定的收敛条件，则输出优化结果，如果不满足收敛条件，则修改拓扑变量，更新有限元模型，返回到第四步，进行下一轮的破损状况有限元分析及优化模型处理，直到满足收敛条件；
第七步，对拓扑变量进行反演，获得最优拓扑构型。
考虑破损-安全条件的频率拓扑优化算例最优拓扑图如图2的(e)所示，阈值选取0.95反演后的最优结构如图2的(f)所示，最优拓扑结构的具体值如表1所示。相比于图2的(d)的传统频率优化最优拓扑结构，可以明显提高结构的冗余度，显著降低最优结构对局部破损的敏感性。从图3的拓扑优化迭代历史可以看出，承受频率响应的路径被保留下来；同时从图3目标体积迭代曲线可以看到，考虑破损条件下的结构总体积经过31次迭代后最终稳定收敛到78180mm3。从图3的各结构破损状况基频的迭代曲线可以看到，各结构破损状况的基频经过迭代后都趋于收敛。因此，本发明提出的方法可以实现考虑破损-安全条件的连续体频率约束拓扑优化设计。
表1考虑破损-安全条件的频率拓扑优化算例最优拓扑结果
1.一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，包括以下步骤：
第一步，确定连续体结构的破损设计域和非破损设计域，并建立相应的有限元模型，设定边界条件，材料属性参数；
第二步，输入破损工况数，并向有限元模型中输入初始破损区域中心坐标值，同时输入破损区域参数，这样就确定了初始破损区域的位置和形状；
第三步，输入频率约束优化参数，包括约束基频的大小、收敛精度、过滤半径，形成考虑破损-安全条件下的频率约束优化模型；
第四步，含有一个结构局部破损区域的基结构称为一个结构破损状况；从初始破损区域开始，根据破损工况数依次对结构进行破损，并对每一个结构破损状况进行模态分析，并提取每一个结构破损状况单元的模态分析结果，为建立显式优化方程提供结构单元的力学性能参数；
第五步，形成考虑破损-安全的以频率为约束、结构重量最小为目标的近似连续数学优化列式；
第六步，采用数学规划算法，对连续数学优化列式进行求解；
第七步，对拓扑变量进行反演，获得最优拓扑构型及对应的基频。
2.根据权利要求1所述的一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，第一步具体为：
首先基于设计需求，确定所要优化连续体结构的最大设计边界，即确定基结构，并依据性能需求，确定破损设计域及非破损设计域；其次，基于msc.patran软件平台，建立基结构几何模型，并划分有限元网格；最后定义材料，赋予单元属性，施加边界条件；最后，设定模态分析输出选项。
3.根据权利要求2所述的一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，第二步具体为：
在基于msc.patran软件内置语言pcl编写的主程序内，首先设置破损工况数，破损工况数根据实际问题确定；然后输入初始破损区域中心二维坐标值；然后输入破损区域参数：损伤区对角线一半的平方、损伤区边长一半的平方，确定正方形的初始破损区域。
4.根据权利要求3所述的一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，第三步具体为：
依据结构固有频率性能设计需求，输入约束参数，包括约束基频大小、收敛精度以及过滤半径；这样破损区域大小、形状以及频率约束就确定下来，建立考虑破损-安全的以基频为约束、结构重量最小为优化目标的拓扑优化模型：
式中：t表示拓扑变量向量，ti为第i个单元的拓扑变量，表征单元的有或无；t∈en表示拓扑变量t是属于n维欧式空间的向量；w表示结构中重量，
5.根据权利要求4所述的一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，第四步具体为：
在pcl主程序中实现从初始破损区域开始，根据破损工况数依次对结构进行破损，并对每一个结构破损状况进行模态分析；依次提取每一个结构破损状况单元的模态分析结果，分别是uil和vil，
6.根据权利要求5所述的一种考虑破损-安全的连续体频率约束拓扑优化设计方法，其特征在于，第五步具体为：
将设计变量取为
其中α＝γw/γk，γw＝2是重量过滤函数的幂指数，γk＝5是单元刚度过滤函数的幂指数；dil和el是显式约束不等式的常系数；ν是迭代次数；
为使拓扑变量t尽量取0或1，加入离散性条件并与原目标按线性加权组成单目标，对目标二次近似后略去常数项得：
这样更新近似连续数学优化列式为：
