一种基于刚度等效和模型降阶的减振优化快速计算方法技术

一种基于刚度等效和模型降阶的减振优化快速计算方法。首先，针对复杂加筋薄壁结构划分有限元周期性加筋单胞结构，建立原薄壁加筋结构同尺寸不含筋条的光筒薄壁结构。其次，基于有限元单胞得到等效刚度系数、等效密度、模型整体结构质量及结构承载力。再次，将等效刚度系数和等效密度参数赋予光筒薄壁结构建立等效模型，再通过建立降阶分析模型得到特定工作频段能量值。最后，建立有限元优化列式，以加筋薄壁结构的壳体厚度和筋条参数为设计变量，以承载力和能量值作为约束，质量最小化为目标进行轻量化设计。本发明专利技术对于多种类型的薄壁加筋结构、多种状况动力响应分析有较高的适用性，且具有满足精度和高效计算的优势，能够协同考虑多个约束条件。

【技术实现步骤摘要】
一种基于刚度等效和模型降阶的减振优化快速计算方法
本专利技术属于航空航天中主承力薄壁构件减振优化设计领域，提出一种基于刚度等效和模型降阶的减振优化快速计算方法。
技术介绍
薄壳结构广泛用于航空航天承力部件，如运载火箭燃料贮箱、飞机机身、航空发动机加力筒体等。典型的薄壳结构包括柱壳、锥壳、球壳、组合壳等。典型的服役工况包括轴压、外压、扭转、轴内压、轴外压、轴扭转等。同时薄壁结构除了承受静态载荷以外，还承受不同类型的动力载荷，在其使用期间，不可避免地受到外界环境的振动干扰，动力损伤是其主要的破坏形式之一。在部分航天部件的运行过程中，还会要求其固有频率远离发动机振动频率，并保证一定的抗击振动能力。对于大型复杂结构，其动载荷下的力学分析计算量比静力分析大几个数量级，求解过程极其耗费时间，构造合适的降阶模型能够在大幅度降低计算时间的同时保证计算动力分析的精度，并解决大型复杂结构的多约束多工况优化迭代中大量动力学计算非常耗时的问题。同时对于复杂薄壁结构的分析和优化，其快速承载力预测的需求也变得日益迫切。为了改善结构或者系统的动态特性，达到控制振动的目的，确保其动力环境下的工作安全可靠，开展动力优化设计方面的研究和应用变得日益迫切，专利技术一种在研发阶段就考虑动力性能并进行减振优化，并综合结构的承载力，多变量协同化设计，实现轻量化目标的快速优化计算方法尤为重要。
技术实现思路
本专利技术主要解决复杂薄壁加筋结构协同承载力和动力需求实现轻量化设计耗时巨大的问题，提出一种基于刚度等效和模型降阶的减振优化快速计算方法。针对优化区间范围内的复杂薄壁加筋结构，首先划分周期性单胞结构，并建立原薄壁加筋结构同尺寸不含筋条的光筒薄壁结构。然后，对单胞结构求解得到等效刚度系数，针对加筋薄壁结构整体屈曲载荷预测得到承载力，对单胞提取质量并经过乘积计算得到模型整体质量，并采用面积等效方法得到等效密度。进而，将等效密度和等效刚度赋予到光筒模型，使用模型降阶方法结合全局减缩基技术，建立降阶数值分析模型，对降阶后的模型进行动力学分析，计算得到特定工作频段能量值。最后，以加筋薄壁结构的壳体厚度和筋条参数为设计变量，结构承载力和工作频段能量值作为约束，质量最轻为目标进行优化设计。为了达到上述目的，本专利技术的技术方案为：一种基于刚度等效和模型降阶的减振优化快速计算方法，包括以下步骤：第一步，针对复杂加筋薄壁结构划分有限元周期性加筋单胞结构，并建立原薄壁加筋结构同尺寸不含筋条的光筒薄壁结构，具体如下：首先，建立复杂加筋薄壁结构，从加筋薄壁结构中划分出代表性周期加筋单胞结构，得到单胞有限元模型，其中单胞结构主要包括蒙皮和筋条两个部分。其次，根据原加筋薄壁结构，建立同样尺寸去除筋条的光筒薄壁结构。最后，对单胞结构和光筒薄壁结构有限元单元网格划分，并对单胞结构赋予材料属性。第二步，基于有限元单胞计算得到等效刚度系数、等效密度、模型整体结构质量以及结构承载力，具体如下：根据单胞有限元模型，基于渐进均匀化理论或者等效刚度法求解加筋单胞的等效刚度系数，将等效刚度系数带入瑞利-里兹公式，求解整体结构线性屈曲载荷值，即为结构承载力。对单胞有限元模型查询质量，用单胞质量除以单胞蒙皮面积得到等效密度，用单胞质量乘以单胞数量得到模型整体结构质量。第三步，将第二步求得的等效刚度系数和等效密度参数赋予划分网格后的光筒薄壁结构，建立等效模型，基于模型降阶方法结合全局减缩基技术建立降阶分析模型，进行动力学谐响应分析，得到特定工作频段能量值，具体如下：首先，将第二步求得的等效刚度和等效密度赋给第一步建立的光筒薄壁结构，获得复杂加筋薄壁结构的等效模型。然后，通过减缩基构造方法构造结构初始减缩基矩阵，再使用模型降阶方法结合全局减缩基技术建立全局减缩基矩阵，通过降低有限元数值计算模型的维度建立降阶分析模型，对有限元数值模型进行自由度减缩，使降阶模型在满足分析精度的同时减少计算时间。最后，对降阶后的模型进行动力学谐响应分析，根据谐响应分析得到的的动刚度和位移响应输出计算频点的能量，对特定工作频段内的所有能量值进行累加得到工作频段能量值。进一步的，所述减缩基矩阵构造方法包括里兹向量法、模态叠加法。进一步的，所述的模型降阶方法包括：特征正价本征正交分解(ProperOrthogonalDecomposition，POD)方法、动力学模态分解(DynamicModeDecomposition，DMD)方法。第四步，建立承载力与动响应协同优化列式并对优化问题进行求解。以加筋薄壁结构的壳体厚度和筋条参数为设计变量，第二步和第三步中得到的结构承载力和工作频段能量值作为约束，质量最小化为目标进行轻量化设计。建立的轻量化优化列式如下：Find：xMinimize：MassSubjectto：xilb≤xi≤xiub,i＝1,...IPcr(x)＞PiniTe(x)＜Teini式中，x为设计变量，包括加筋薄壁结构的壳体厚度和筋条参数；xilb和xiub分别是变量x的下限值和上限值。Mass为待求结构整体质量，Pcr(x)和Te(x)分别是待求结构整体承载力和工作频段内的总能量，Pini和Teini分别是原始构型的整体承载力和工作频段能量值。四步中在优化计算方法上，首先，采用聚类方法对优化设计空间减缩，以便于寻优。其次，根据缩小范围之后的参数变量区间进行局部寻优处理。最后，采用局部优化方法求解最优解，经过优化得到满足结构承载力许可，得到实现减振功能性需求和轻量化设计的最优构型。进一步的，所述的聚类方法包括:K均值聚类方法、模糊C均值聚类方法、SOM神经网络聚类方法。本专利技术的有益效果为：本专利技术针对复杂薄壁加筋结构协同承载力和动力需求实现轻量化设计耗时巨大的问题，通过对复杂加筋薄壁结构划分周期性加筋单胞结构，求解等效刚度系数和等效密度，基于整体屈曲载荷预测方法和快速质量计算方法快速获得结构承载力和质量信息，基于等效模型与模型降阶方法结合全局模型降阶技术实现快速动力学谐响应分析，并在优化求解过程中使用聚类方法对设计变量空间减缩，加速了优化进程，快速求解，实现综合考虑结构承载力、动力学响应协同优化设计的目的。本专利技术提出的方法能实现复杂薄壁加筋结构减振轻量化快速计算，实现结构承载力和动力学响应协同设计，在满足结构功能需求的同时，减轻结构质量，缩短设计时间周期，降低材料成本，提高经济效益。附图说明图1为本专利技术的实现流程图；图2为本专利技术加筋薄壁模型实例提供的初始设计构型：图2(a)为正置正交加筋圆柱壳模型，图2(b)为不含筋条的光筒薄壁模型，图2(c)为周期性加筋单胞结构模型。具体实施方式为使本专利技术解决的方法问题、采用的方法方案和达到的方法效果更加清楚，下面结合附图和实施例对本专利技术作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本专利技术，而非对本专利技术的限定。另外还需要说明的是，为了便于描述，附图中仅示出本文档来自技高网...

【技术保护点】
1.一种基于刚度等效和模型降阶的减振优化快速计算方法，其特征在于，包括以下步骤：/n第一步，针对复杂加筋薄壁结构划分有限元周期性加筋单胞结构，并建立原薄壁加筋结构同尺寸不含筋条的光筒薄壁结构；/n第二步，基于有限元单胞计算得到等效刚度系数、等效密度、模型整体结构质量以及结构承载力；/n第三步，将第二步求得的等效刚度系数和等效密度参数赋予划分网格后的光筒薄壁结构，建立等效模型，基于模型降阶方法结合全局减缩基技术建立降阶分析模型，进行动力学谐响应分析，得到特定工作频段能量值，具体如下：/n首先，将第二步求得的等效刚度和等效密度赋给第一步建立的光筒薄壁结构，获得复杂加筋薄壁结构的等效模型；然后，通过减缩基矩阵构造方法构造结构初始减缩基矩阵，再使用模型降阶方法结合全局减缩基技术建立全局减缩基矩阵，通过降低有限元数值计算模型的维度建立降阶分析模型，对有限元数值模型进行自由度减缩，使降阶模型在满足分析精度的同时减少计算时间；最后，对降阶后的模型进行动力学谐响应分析，根据谐响应分析得到的的动刚度和位移响应输出计算频点的能量，对特定工作频段内的所有能量值进行累加得到工作频段能量值；/n第四步，建立承载力与动响应协同优化列式并对优化问题进行求解；以加筋薄壁结构的壳体厚度和筋条参数为设计变量，第二步和第三步中得到的结构承载力和工作频段能量值作为约束，质量最小化为目标进行轻量化设计；/n建立的轻量化优化列式如下：/nFind：x/nMinimize：Mass/nSubject to：x...

【技术特征摘要】
1.一种基于刚度等效和模型降阶的减振优化快速计算方法，其特征在于，包括以下步骤：
第一步，针对复杂加筋薄壁结构划分有限元周期性加筋单胞结构，并建立原薄壁加筋结构同尺寸不含筋条的光筒薄壁结构；
第二步，基于有限元单胞计算得到等效刚度系数、等效密度、模型整体结构质量以及结构承载力；
第三步，将第二步求得的等效刚度系数和等效密度参数赋予划分网格后的光筒薄壁结构，建立等效模型，基于模型降阶方法结合全局减缩基技术建立降阶分析模型，进行动力学谐响应分析，得到特定工作频段能量值，具体如下：
首先，将第二步求得的等效刚度和等效密度赋给第一步建立的光筒薄壁结构，获得复杂加筋薄壁结构的等效模型；然后，通过减缩基矩阵构造方法构造结构初始减缩基矩阵，再使用模型降阶方法结合全局减缩基技术建立全局减缩基矩阵，通过降低有限元数值计算模型的维度建立降阶分析模型，对有限元数值模型进行自由度减缩，使降阶模型在满足分析精度的同时减少计算时间；最后，对降阶后的模型进行动力学谐响应分析，根据谐响应分析得到的的动刚度和位移响应输出计算频点的能量，对特定工作频段内的所有能量值进行累加得到工作频段能量值；
第四步，建立承载力与动响应协同优化列式并对优化问题进行求解；以加筋薄壁结构的壳体厚度和筋条参数为设计变量，第二步和第三步中得到的结构承载力和工作频段能量值作为约束，质量最小化为目标进行轻量化设计；
建立的轻量化优化列式如下：
Find：x
Minimize：Mass
Subjectto：xilb≤xi≤xiub,i＝1,...I
Pcr(x)＞Pini
Te(x)＜Teini
式中，x为设计变量，包括加筋薄壁结构的壳体厚度和筋条参数；xilb和xiub分别是变量x的下限值和上限值；Mass为待求结构整体质量，Pcr(x)和Te(x)分别是待求结构整体承载力和工作频段内的总能量，Pini和Te...

【专利技术属性】
技术研发人员：可钊，田阔，李增聪，石云峰，周才华，杜凯繁，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21