技术摘要：
本发明公开了考虑频带约束的结构动力学拓扑优化方法、系统及应用，其技术方案为：包括建立频带约束表达式，并计算频带约束表达式的敏度；根据敏度信息将频带约束表达式集成到动力学拓扑优化模型中，获得设定频率范围下禁带结构。本发明基于改进的Heaviside函数，以建立  全部
背景技术：
机床在高速切削过程中，往复运动部件在改变切削方向时会引发惯性冲击、切削 冲击等外部载荷，造成机床结构振动。机床的振动会影响加工质量、减少刀具寿命并产生噪 声污染，因而拟制加工过程中的振动现象成为重要的研究方向。当外激励频率等于或接近 结构的固有频率时，会产生共振现象，因此合理的调整结构固有频率使其避开加工常用频 率，可有效的拟制机床结构的加工振动变形，以提升机床加工精度并减少加工噪音。 随着计算机技术的飞速发展，机床结构的设计除了依靠传统经验外，采用结构优 化设计的思想进行改进已成为高端机床研究的重点。其中，拓扑优化方法因无需初始构型， 可获得特定工况或性能需求下的最优概念构型，已成为机床结构概念设计的重要手段。 结构频率优化设计一直是拓扑优化领域的重点研究内容。发明人发现，现有技术 中提出了一种基频最大的结构拓扑优化方法，但此方法仅仅针对结构的第一阶固有频率进 行最大化设计，并未考虑结构的固有频率是否避开常用的加工频率。现有技术还提出了频 率间隙最大化的自由振动拓扑优化方法，虽然此方法考虑了相邻阶频率差的最大，但是所 考虑的结构模态阶数需预先确定，且无法给出确定的频带范围。
技术实现要素：
针对现有技术存在的不足，本发明的目的是提供一种考虑频带约束的结构动力学 拓扑优化方法、系统及应用，基于改进的Heaviside函数，以建立具有连续、可微等特点的频 带约束数学表达列式；通过推导该函数的敏度，将其集成到传统的动力学拓扑优化模型中， 实现了指定频率范围下结构禁带设计。 为了实现上述目的，本发明是通过如下的技术方案来实现： 第一方面，本发明的实施例提供了一种考虑频带约束的结构动力学拓扑优化方 法，包括： 建立频带约束表达式，并计算频带约束表达式的敏度； 根据敏度信息将频带约束表达式集成到动力学拓扑优化模型中，获得设定频率范 围下禁带结构。 作为进一步的实现方式，首先确定待设计结构并划分网格，定义拓扑优化的设计 域；添加边界条件，计算单元刚度矩阵和单元质量矩阵，并建立总体刚度阵和总体质量阵； 求解结构特征频率和特征向量，根据要施加的频率禁带构造频带约束函数。 作为进一步的实现方式，对频带约束函数进行光滑化处理，得到相应频率对应的 频带约束函数值；采用P范数求解所有频带约束函数值的最大值。 作为进一步的实现方式，将频带约束函数公式放入优化列式中，在对应的拓扑优 3 CN 111597724 A 说　明　书 2/7 页 化问题中定义目标函数为结构的基频最大，约束条件是满足频带约束和体积约束，建立拓 扑优化模型。 作为进一步的实现方式，求解所述最大值对设计变量的一阶导数，采用MMA算法对 拓扑优化模型中的设计变量进行迭代更新，完成对结构的优化设计。 作为进一步的实现方式，采用SIMP插值格式组装建立总体刚度阵K和总体质量阵 M，其中，SIMP插值中的刚度惩罚系数为pK、质量惩罚系数为pM。 第二发明，本发明实施例还提供了一种考虑频带约束的结构动力学拓扑优化方法 在高端机床结构设计中的应用。 第三方面，本发明实施例还提供了一种考虑频带约束的结构动力学拓扑优化系 统，包括： 频带约束表达式建立模块，用于建立频带约束表达式，并计算频带约束表达式的 敏度； 禁带结构获取模块，用于根据敏度信息将频带约束表达式集成到动力学拓扑优化 模型中，获得设定频率范围下禁带结构。 第四方面，本发明实施例还提供了一种电子设备，包括存储器、处理器及存储在存 储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的考虑频 带约束的结构动力学拓扑优化方法。 第五方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机 程序，该程序被处理器执行时实现所述的考虑频带约束的结构动力学拓扑优化方法。 上述本发明的实施例的有益效果如下： (1)本发明的一个或多个实施方式基于改进的Heaviside函数建立具有连续、可微 等特点的频带约束数学表达列式；通过推导该函数的敏度，将其集成到传统的动力学拓扑 优化模型中，实现了指定频率范围下结构禁带设计； (2)本发明的一个或多个实施方式获得的设计结构，其固有频率出现在外激励载 荷频段范围之外，保证了设计结构在实际工作运行过程中避免共振现象，提高加工精度并 减振降噪； (3)本发明的一个或多个实施方式可应用到机床结构的设计中，以提升机床结构 的动力学性能，实现高端机床装备的进一步升级。 附图说明 构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示 意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。 图1是本发明根据一个或多个实施方式的流程图； 图2是本发明根据一个或多个实施方式的设计域模型示意图； 图3(a)是本发明根据一个或多个实施方式结构无频带约束的拓扑优化结果图； 图3(b)是本发明根据一个或多个实施方式结构无频带约束的材料拓扑图； 图4是本发明根据一个或多个实施方式结构无频带约束的前九阶频率的迭代曲线 图； 图5(a)-图5(i)是本发明根据一个或多个实施方式结构无频带约束的前九阶模态 4 CN 111597724 A 说　明　书 3/7 页 图； 图6(a)是本发明根据一个或多个实施方式结构施加双频带约束的拓扑优化结果 图； 图6(b)是本发明根据一个或多个实施方式结构施加双频带约束的材料拓扑图； 图7是本发明根据一个或多个实施方式结构施加双频带约束的前九阶频率的迭代 曲线图； 图8(a)-图8(i)是本发明根据一个或多个实施方式结构施加双频带约束的前九阶 模态图。