一种基于水平集拓扑优化的局部模态识别方法技术

本发明专利技术公开了一种基于水平集拓扑优化的局部模态识别方法，其包括以下步骤：采用Delaunay三角剖分算法划分非结构化背景网格；通过水平集函数定义的设计边界来修改所述非结构化背景网格，以得到实际结构内部的非结构化三角网格；利用实际结构内部的非结构化三角网格进行有限元模态分析，取出前n阶模态振型与振动频率，其中n为正整数；对每阶模态振型定义未振动区域位移；将三角网格单元中心位移符合未振动区域位移定义的单元面积叠加求和以得到面积Vtiny；计算面积Vtiny与三角网格单元总面积Ventire的比值，并将所述比值与设定阀值进行比较，进而判断局部模态。

【技术实现步骤摘要】

本专利技术属于结构优化相关
，更具体地，涉及一种基于水平集拓扑优化的局部模态识别方法。
技术介绍
结构拓扑优化设计方法是新型的数字化结构设计方式，通过建立包含目标函数与约束方程的数学模型，对结构进行有限元数值分析，在设计域中按照优化准则或者数学规划方法迭代出满足目标要求的材料分布。结构拓扑优化设计方法广泛运用于车身设计、航空航天、通信电子等结构设计。在动力学结构拓扑优化领域，有效的避免局部模态出现的优化方法是研究的重要课题，是实现和模态相关目标函数优化的关键，能够保证优化目标的一致性与正确性，避免目标函数的改变带来的优化失败。在动力学结构拓扑优化中，如常见的对一阶频率的优化或则关于临界屈曲载荷的优化中，一个关键的步骤是进行结构的模态分析，得到的某阶或者某几阶模态振型用于相关的目标或者约束的敏度计算。局部模态是指模态振型在结构的大部分自由度上的幅值接近零，仅在少部分自由度上存在明显的幅值。在有限元模态分析计算中，对于没有局部模态的结构，可以通过模态阶数表达目标模态。然而，拓扑优化中涉及结构的演化不可避免会产生局部模态，此时有限元计算后的模态阶数不可直接用于特定模态的表达，需要剔除其中的局部模态，否则会产生涉及目标的不一致而致使优化失去意义。目前，常用的拓扑优化方法是SIMP方法(固体各向同性微结构惩罚法)，通过材料相对密度表示结构拓扑。设计域中材料相对密度为1表示为实际材料，相对密度为0表示为孔洞，SIMP优化过程中材料相对密度介于0～1之间，优化中施加惩罚因子使相对密度趋近0～1区间的整数。材料性质与相对密度满足幂函数关系，在有限元计算中由于弱材料(低密度材料)刚度矩阵的惩罚因子大于质量矩阵的惩罚因子，导致伴随低频率值的局部模态的出现。这些局部模态只反映弱材料区域的振动，并非是需要的实际结构振动，且占据有限元模态分析的前几阶次。通常的处理方式是调节惩罚因子的比值或者采用不同的插值方式使局部模态出现在高阶振动，这种办法并没有完全去除局部模态，并且调节惩罚因子会增加目标函数收敛的难度。
技术实现思路
针对现有技术的以上缺陷或改进需求，本专利技术提供了一种基于水平集拓扑优化的局部模态识别方法，其基于水平集拓扑优化及有限元分析，针对局部模态识别方法进行了设计。所述局部模态识别方法在实际结构内部使用与设计边界相适应的非结构化三角网格进行有限元分析，避免了人工弱材料带来的局部模态现象；同时，通过计算每阶模态的未振动区域的体积占总体积的百分比，有效筛选出由于结构边界变化带来的局部模态，确保了优化的准确性。为实现上述目的，按照本专利技术的一个方面，提供了一种基于水平集拓扑优化的局部模态识别方法，其包括以下步骤：(1)采用Delaunay三角剖分算法划分非结构化背景网格；(2)通过水平集函数定义的设计边界来修改所述非结构化背景网格，以得到实际结构内部的非结构化三角网格；(3)利用实际结构内部的非结构化三角网格进行有限元模态分析，取出前n阶模态振型与振动频率，其中n为正整数；(4)对每阶模态振型定义未振动区域位移；(5)将三角网格单元中心位移符合未振动区域位移定义的单元面积叠加求和以得到面积Vtiny；(6)计算面积Vtiny与三角网格单元总面积Ventire的比值，并将所述比值与设定阀值进行比较，进而判断局部模态。进一步的，所述设定阀值为0.2。进一步的，所述比值小于等于所述设定阀值的模态为主模态；所述比值大于所述设定阀值的模态为局部模态。进一步的，对每阶模态振型定义未振动区域位移包括以下步骤：(4-1)对于每阶振动模态，利用有限元计算得到振型，并找出每个结点的所有自由度，求出每个结点的位移并遍历以得到最大结点位移；(4-2)最小结点位移为所述最大结点位移乘以设定系数，所述未振动区域位移定义为小于等于所述最小结点位移的位移。进一步的，所述设定系数为0.001。进一步的，对所述非结构化背景网格的修改包括删除不符合设定条件的三角网格单元及对应的结点的步骤。总体而言，通过本专利技术所构思的以上技术方案与现有技术相比，本专利技术提供的基于水平集拓扑优化的局部模态识别方法，其在实际结构内部使用与设计边界相适应的非结构化三角网格进行有限元分析，避免了人工弱材料带来的局部模态现象；同时，通过计算每阶模态的未振动区域的体积占总体积的百分比，有效筛选出由于结构边界变化带来的局部模态，确保了优化的准确性。附图说明图1是本专利技术较佳实施方式提供的基于水平集拓扑优化的局部模态识别方法的流程图。图2是图1中的基于水平集拓扑优化的局部模态识别方法涉及的集中质量悬臂梁一阶频率优化的示意图。图3是图1中的基于水平集拓扑优化的局部模态识别方法涉及的背景网格的示意图。图4是图1中的基于水平集拓扑优化的局部模态识别方法涉及的用于有限元分析的结构内部与设计边界相适应的网格示意图。图5是图1中的基于水平集拓扑优化的局部模态识别方法涉及的修改背景网格得到结构内部网格方法的示意图。图6是图1中的基于水平集拓扑优化的局部模态识别方法涉及的有限元分析得到的第5阶模态振型示意图。在所有附图中，相同的附图标记用来表示相同的元件或结构，其中：1-集中质量块，2-零水平集，3-背景网格，4-未振动区域，5-最大结点位移。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。此外，下面所描述的本专利技术各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。请参阅图1至图4，本专利技术较佳实施方式提供的基于水平集拓扑优化的局部模态识别方法，其利用实际结构内部的非结构化网格进行模态分析，同时对每阶候选模态进行计算，计算未振动区域的体积占总体积的百分比，避免了弱材料引起的局部模态及结构边界变化引发的局部模态。本实施方式中，以带有集中质量悬臂梁一阶频率最大化的优化问题为例来解释本专利技术。在给定的矩形设计域内给定材料体积约束，区域左侧边界施加位移约束，区域右下角施加集中质量块1。对所述悬臂梁进行的优化迭代进行到15步时出现了局部模态，此时有限元计算的一阶频率并非是优化目标需要的主结构一阶频率，而是局部模态相应的频率，需要剔除。请一并参阅图5及图6，本实施方式中，所述基于水平集拓扑优化的局部模态识别方法，其包括以下步骤：步骤一，采用Delaunay三角剖分算法划分非结构化背景网格3。具体地，利用Delaunay三角剖分算法在所述矩形设计域内划分所述非结构化背景网格3。步骤二，通过水平集函数定义的设计边界来修改所述非结构化背景网格3，以得到实际结构内部的非结构化三角网格。具体地，首先，找出所述非结构化背景网格3中被设计边界(即零水平集)2截断的单元边，并收集单元边界上的交叉点集P0；找出所述非结构化背景网格3在实际结构Ω内的结点集P1，将所述交叉点集P0与所述结点集P1进行合并，同时剔除重复点以得到点集P；随后，利用Delaunay三角剖分算法剖分所述点集P以得到与所述设计边界2相适应的三角网格单元；随后，利用所述三角网格单元的三个结点坐标与所述水平集函数插值以得到每个所述三角网格单元中心处的水平集函数值，当所述三角网格单元中心处的水平集函数值小于设定值时，删除此本文档来自技高网...

【技术保护点】
一种基于水平集拓扑优化的局部模态识别方法，其包括以下步骤：(1)采用Delaunay三角剖分算法划分非结构化背景网格；(2)通过水平集函数定义的设计边界来修改所述非结构化背景网格，以得到实际结构内部的非结构化三角网格；(3)利用实际结构内部的非结构化三角网格进行有限元模态分析，取出前n阶模态振型与振动频率，其中n为正整数；(4)对每阶模态振型定义未振动区域位移；(5)将三角网格单元中心位移符合未振动区域位移定义的单元面积叠加求和以得到面积Vtiny；(6)计算面积Vtiny与三角网格单元总面积Ventire的比值，并将所述比值与设定阀值进行比较，进而判断局部模态。

【技术特征摘要】
1.一种基于水平集拓扑优化的局部模态识别方法，其包括以下步骤：(1)采用Delaunay三角剖分算法划分非结构化背景网格；(2)通过水平集函数定义的设计边界来修改所述非结构化背景网格，以得到实际结构内部的非结构化三角网格；(3)利用实际结构内部的非结构化三角网格进行有限元模态分析，取出前n阶模态振型与振动频率，其中n为正整数；(4)对每阶模态振型定义未振动区域位移；(5)将三角网格单元中心位移符合未振动区域位移定义的单元面积叠加求和以得到面积Vtiny；(6)计算面积Vtiny与三角网格单元总面积Ventire的比值，并将所述比值与设定阀值进行比较，进而判断局部模态。2.如权利要求1所述的基于水平集拓扑优化的局部模态识别方法，其特征在于：所述设定阀值为0.2。3.如权利要求2所述的基于水平集拓扑优化的局部模...

【专利技术属性】
技术研发人员：夏奇，李振华，史铁林，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：湖北;42