考虑频带约束的结构动力学拓扑优化方法、系统及应用技术方案

本发明专利技术公开了考虑频带约束的结构动力学拓扑优化方法、系统及应用，其技术方案为：包括建立频带约束表达式，并计算频带约束表达式的敏度；根据敏度信息将频带约束表达式集成到动力学拓扑优化模型中，获得设定频率范围下禁带结构。本发明专利技术基于改进的Heaviside函数，以建立具有连续、可微等特点的频带约束数学表达列式；通过推导该函数的敏度，将其集成到传统的动力学拓扑优化模型中，实现了指定频率范围下结构禁带设计。

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

【技术保护点】
1.考虑频带约束的结构动力学拓扑优化方法，其特征在于，包括：/n建立频带约束表达式，并计算频带约束表达式的敏度；/n根据敏度信息将频带约束表达式集成到动力学拓扑优化模型中，获得设定频率范围下禁带结构。/n

【技术特征摘要】
1.考虑频带约束的结构动力学拓扑优化方法，其特征在于，包括：
建立频带约束表达式，并计算频带约束表达式的敏度；
根据敏度信息将频带约束表达式集成到动力学拓扑优化模型中，获得设定频率范围下禁带结构。


2.根据权利要求1所述的考虑频带约束的结构动力学拓扑优化方法，其特征在于，首先确定待设计结构并划分网格，定义拓扑优化的设计域；添加边界条件，计算单元刚度矩阵和单元质量矩阵，并建立总体刚度阵和总体质量阵；
求解结构特征频率和特征向量，根据要施加的频率禁带构造频带约束函数。


3.根据权利要求2所述的考虑频带约束的结构动力学拓扑优化方法，其特征在于，对频带约束函数进行光滑化处理，得到相应频率对应的频带约束函数值；采用P范数求解所有频带约束函数值的最大值。


4.根据权利要求3所述的考虑频带约束的结构动力学拓扑优化方法，其特征在于，将频带约束函数公式放入优化列式中，在对应的拓扑优化问题中定义目标函数为结构的基频最大，约束条件是满足频带约束和体积约束，建立拓扑优化模型。


5.根据权利要求4所述的考虑频带约束的结构动力学拓扑优化方法，其特征在于，求解所述最大值对设计变量的一阶导数，采...

【专利技术属性】
技术研发人员：李取浩，吴强波，张松，李建勇，满佳，
申请(专利权)人：山东大学，
类型：发明
国别省市：山东;37