一种含局部振动模态识别的结构优化方法和系统与流程

本发明属于空气过滤器技术领域，具体涉及一种波纹滤筒式空气过滤器，可用于暖通空调或其他需要对空气净化除尘产品中。

背景技术：

技术领域

本发明涉及一种结构优化方法，尤其涉及一种含局部振动模态识别的结构优化方法和系统。

背景技术

现代结构优化方法结合了数学规划理论和力学分析方法，以计算机为工具，可以自动地优化和改进结构设计，并满足各种条件限制，被广泛应用于航空航天、汽车和建筑等领域，其潜在经济效益十分显著，相关结构优化方法的研究已受到国内外的高度重视。

在工程结构优化过程中，往往需要对不同方向上局部振动的频率进行约束。但是，随着设计变量的迭代变化，局部振动模态容易发生模态置换现象，使原来的频率约束错误地附加到其他模态上。这种模态置换的反复出现，导致优化问题不能收敛；即使能够收敛，最后的优化结构也是不可信的。

现有的结构优化方法无法解决上述局部振动模态的模态置换的问题，因此，有必要设计一套新的结构优化方法和系统，以克服上述问题。

技术实现要素：

本发明的目的是提供一种新的含局部振动模态识别的结构优化方法和系统，以在发生局部模态置换现象时，自动识别目标模态，修改频率约束，保证迭代过程的顺利进行，避免模态置换。

为了实现上述目的，本发明提供的含局部振动模态识别的结构优化方法包括以下步骤：步骤S1，确定有限元模型和优化模型，其中，初始的有限元模型 和优化模型按照初始设计的结构参数建立，并在优化模型中添加局部模态频率约束；步骤S2，对有限元模型进行模态分析；步骤S3，根据模态分析的分析结果搜寻局部模态频率约束中指定的局部振动模态；步骤S4，若未搜寻到局部模态频率约束中指定的所有自由度的局部振动模态，则增加模态分析求解阶数，返回步骤S2重新进行模态分析；步骤S5，若搜寻到所有自由度的局部振动模态，采用所有自由度的局部振动模态更新优化模型中的局部模态频率约束；步骤S6，建立近似优化问题；步骤S7，采用更新局部模态频率约束后的优化模型进行寻优，其中，若近似优化问题收敛，则停止寻优，否则返回步骤S1，根据当前寻优结果更新优化模型的设计变量并修改有限元模型的对应结构参数，以重新确定有限元模型和优化模型。

进一步地，步骤S2中模态分析的分析结果为每个自由度的多个模态的单元模态应变能信息；步骤S3具体为：判断一个自由度的局部振动区域的多个模态的单元模态应变能与整体的单元模态应变能的比值是否大于或等于预定阈值，确定大于素数预定阈值的比值所对应的模态为局部振动模态。

进一步地，步骤S6具体为：对有限元模型进行结构分析和敏度分析；根据结构分析和敏度分析的分析结果建立近似优化问题。

进一步地，步骤S1中建立的优化模型的数学表达式为：

其中，X＝{x₁，x₂，...，x_n}^T，xi为独立的设计变量，f(X)和g(X)分别表示目标函数和约束函数，n为总的约束个数，和对应设计变量xi的上、下限。

进一步地，约束函数为：

其中，为预期达到的局部模态频率，ω_local为被约束的局部模态的频率。

进一步地，步骤S7中寻优时采用二级多点逼近优化算法。

进一步地，该方法还包括：在停止寻优之后，重新对有限元模型进行模态分析和后处理，以检验优化结果是否满足设计要求。

为了实现上述目的，本发明提供的含局部振动模态识别的结构优化包括： 有限元前后处理模块，用于确定有限元模型和优化模型，其中，初始的有限元模型和优化模型按照初始设计的结构参数建立，并在优化模型中添加局部模态频率约束；结构分析模块，用于对有限元模型进行模态分析；模态识别模块，用于根据模态分析的分析结果搜寻局部模态频率约束中指定的局部振动模态，并采用搜寻到的所有自由度的局部振动模态更新优化模型中的局部模态频率约束；寻优模块，用于建立近似优化问题，并对更新局部模态频率约束后的优化模型进行寻优，其中，若近似优化问题收敛，则停止寻优；其中，结构分析模块还用于在模态识别模块未搜寻到所有自由度的局部振动模态时，增加模态分析求解阶数，重新进行模态分析；有限元前后处理模块还用于在近似优化问题不收敛时，根据当前寻优结果更新优化模型的设计变量并修改有限元模型的对应结构参数，以重新确定有限元模型和优化模型。

进一步地，结构分析模块进行模态分析的分析结果为每个自由度的多个模态的单元模态应变能信息；模态识别模块具体用于判断一个自由度的局部振动区域的多个模态的单元模态应变能与整体的单元模态应变能的比值是否大于或等于预定阈值，确定大于素数预定阈值的比值所对应的模态为局部振动模态。

进一步地，结构分析模块还用于对有限元模型进行结构分析和敏度分析，以使寻优模块根据结构分析和敏度分析的分析结果建立近似优化问题。

本发明通过一种含局部振动模态识别的结构优化方法，与现有技术相比，该方法对建立的有限元模型进行模态分析和模态识别，以识别出局部振动模态，然后采用识别出的局部振动模态更新优化模型的局部模态频率约束，使优化模型的局部模态频率约束附加到正确的模态上，克服了优化过程中局部模态置换现象导致的约束错误问题，保证优化时迭代过程的顺利进行，可以应用于各种工程结构的优化设计。

附图说明

图1为本发明第一实施例所提供的含局部振动模态识别的结构优化方法的流程图；

图2为本发明第二实施例所提供的含局部振动模态识别的结构优化方法的流程图；

图3为本发明第三实施例所提供的含局部振动模态识别的结构优化系统的框图；

图4为本发明第三实施例所提供的含局部振动模态识别的结构优化系统的工作流程示意图。

具体实施方式

为使本发明的上述目的、特征和优点能够更加明显易懂，在下面的具体实施方式中，将对本发明作进一步详细的说明。

第一实施例

该实施例提供了一种含局部振动模态识别的结构优化方法，如图1所示，该方法包括如下的步骤。

步骤S1，确定有限元模型和优化模型。

该步骤中，在优化初始，先按照初始设计的结构参数建立有限元模型和优化模型，并在优化模型中添加局部模态频率约束；在优化过程中，如果在步骤S7中寻优时近似优化问题不能收敛，则根据当前寻优结果更新优化模型的设计变量并修改有限元模型的对应结构参数，重新确定有限元模型和优化模型。

步骤S2，对有限元模型进行模态分析。

例如，调用现有的有限元分析软件执行模态分析。

步骤S3，根据模态分析的分析结果搜寻局部模态频率约束中指定的局部振动模态。

一般来说，在一个自由度上，局部振动区域的单元模态应变能与整体的单元模态应变能的比值相对较大，基于这一特性，在搜寻局部振动模态时，可根据分析结果中的单元模态应变能信息来确定局部振动模态。

步骤S4，若未搜寻到局部模态频率约束中指定的所有自由度的局部振动模态，则增加模态分析求解阶数。

为了避免任意模态的模态置换问题，在该实施例中，需要搜寻到所有自由度的局部振动模态，如果未搜寻到，需要在增加模态分析求解结束之后，返回上述步骤S2重新进行模态分析，然后步骤S3根据新的模态分析结果进行搜寻。由此，可以识别出各个自由度上的局部振动模态分别位于哪一阶。

步骤S5，若搜寻到所有自由度的局部振动模态，采用所有自由度的局部振动模态更新优化模型中的局部模态频率约束。

在该步骤中，根据识别出来的结果更新局部模态频率约束，使优化模型的局部模态频率约束附加到正确的模态上。

步骤S6，建立近似优化问题。

以建立的近似优化问题作为寻优迭代过程的收敛条件。

步骤S7，采用更新局部模态频率约束后的优化模型进行寻优。其中，若近似优化问题收敛，则停止寻优，得到的寻优结果作为优化结果进一步检验结构设计是否满足要求。否则，也即近似优化问题未收敛时，返回步骤S1，根据当前寻优结果更新优化模型的设计变量并修改有限元模型的对应结构参数，以重新确定有限元模型和优化模型，后续步骤S2至步骤S7完成重新寻优的过程。

采用该实施例提供的方法，对建立的有限元模型进行模态分析和模态识别，以识别出局部振动模态，然后采用识别出的局部振动模态更新优化模型的局部模态频率约束，使优化模型的局部模态频率约束附加到正确的模态上，克服了优化过程中局部模态置换现象导致的约束错误问题，保证优化时迭代过程的顺利进行。

第二实施例

在上述第一实施例的基础上，该第二实施例提供了一种含局部振动模态识别的结构优化方法的优选实施例，如图2所示，该方法包括的步骤详细说明如下。

根据初始设计方案，建立结构的有限元模型，并建立优化模型，添加局部模态频率约束，其中，优化模型的数学表达式为：

其中，X＝{x₁，x₂，...，x_n}^T，xi为独立的设计变量，该设计变量指具体设计结构的参数，例如梁的横截面尺寸或者板壳的厚度等。f(X)和g(X)分别表示目标函数和约束函数，n为总的约束个数，和是对应设计变量xi的上、下限，min：f(X)表示f(X)取最小值，s.t.g(X)≤0表示受约束条件为g(X)≤0。

随着实际航天器总体设计复杂程度的提高，对结构局部振动的模态频率的要求更加严格。因此，优化模型中的约束函数为

其中为预期达到的局部模态频率，ω_local为被约束的局部模态的频率。

确定模态应变能比的阈值，用以搜寻局部振动模态。由于局部振动模态在局部区域表现出大振幅，应变被集中在这块局部区域的点上，亦即相应的应变能，相反地，总体模态的应变能很均匀地分布于整个系统。一般来说，局部振动区域的单元模态应变能与整体的单元模态应变能的比值相对较大，该实施例通过设定一个阈值以区分局部振动的模态。如果局部区域的单元模态应变能与整体的单元模态应变能的比值不小于该阈值，则认为该模态为局部振动模态。

调用有限元分析软件执行模态分析，得到当前设计点下模态应变能比的信息。分析结束后，根据分析结果搜寻约束中指定的局部振动模态。具体地，采用上述阈值进行搜寻。

判断是否找到所有需要约束的局部振动模态，也即局部模态频率约束中指定的所有自由度的局部振动模态。若是，则更新约束，执行下述步骤5)；否则，增加模态分析求解阶数，返回步骤3)再次进行模态分析；

进行结构分析和敏度分析，以建立近似优化问题，供寻优程序求解。

优化模型采用预定的优化算法进行优化，检查近似问题是否收敛。若是，则停止，进行后处理；否则，根据当前寻优结果更新设计变量，修改对应结构参数，更新有限元模型和优化模型，返回步骤3)重新进行模态分析。

近似问题收敛得到优化结果之后，重新对有限元模型进行模态分析，并进行后处理，以检验优化结果是否满足要求。

该实施例提供的含局部振动模态识别的结构优化方法，基于单元模态应变能信息进行模态识别，具有计算量少、易于使用的优点，可以应用到各种工程结构的设计过程中，以获得最优方案，缩短设计周期，节约成本，同时，应变能是模态振幅的平方，这种应变能比的方法精确性高，花费的计算精力小，实际工程意义更大。

第三实施例

如图3和4所示，该实施例提供了一种含局部振动模态识别的结构优化系统，包括有限元前后处理模块、结构分析模块、模态识别模块和寻优模块四个模块。各个模块之间通过PCL(Patran Connand Language)环境下编写的框架程序相互协调和对接。

有限元前后处理模块采用Patran软件，作为整个优化系统的工作平台。模型的前后处理也在该软件中完成。结构分析模块采用Nastran软件，鉴于其和Patran环境兼容性最佳，用于执行必要的模态分析、结构分析和敏度分析，得到相关的结果信息用于寻优。

模态识别模块用于自动提取模态分析结果文件中的单元模态应变能信息，由此可以识别出各个自由度上的局部振动模态位于哪一阶。并根据识别出来的结果，更新优化模型的局部模态频率约束，使其附加到正确的模态上。

寻优模块为基于二级多点逼近优化算法，采用Fotran语言开发。该模块的输入信息包括设计变量值、目标函数和约束函数的当前值、敏度信息以及寻优控制参数；寻优之后，会自动修改文件中储存的设计变量值，以供后续的迭代使用。

在优化时，先通过有限元前后处理模块建立有限元模型和优化模型，然后调用结构分析模块进行模态分析，模态识别模块读取分析结果，识别局部模态后更新约束条件，之后结构分析模块执行必要的结构分析和敏度分析，供寻优模块寻优求解，如果结果满足收敛条件，则停止，否则返回有限元前后处理模块更新设计变量和结构参数，进行下一次循环。

本申请中的结构优化方法，包括模态识别和约束更新的步骤，保证优化模型的局部模态频率约束附加到正确的模态上，因而能够解决优化过程中局部模态置换现象导致的约束错误问题，使得优化时迭代过程的顺利进行。

对于前述的各方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域的技术人员应该知悉，本发明并不受所描述的动作顺序的限制，因为根据本发明，某些步骤可以采用其他顺序或同时执行；其次，本领域技术人员也应该知悉，上述方法实施例均属于优选实施例，所涉及的动作和模块并不一定是本发明所必须的。

对于前述的各系统实施例，为了简单描述，故将其都表述为一系列的模块组合，但是本领域的技术人员应该知悉，本发明并不受所描述的模块组合的限制，因为根据本发明，某些模块可以采用其他模块执行；其次，本领域技术人员也应该知悉，上述系统实施例均属于优选实施例，所涉及的模块并不一定是本发明所必须的。

本说明书中的各个实施例均采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。 对于系统实施例而言，由于其与方法实施例基本相似，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

以上对本发明所提供的一种含局部振动模态识别的结构优化方法和系统进行了详细介绍，本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。