基于子结构法的周期直梁有限元单元等效表征与低频动力学验证

《Finite Elements in Analysis and Design》：Substructuring-based characterization of periodic straight beams from finite-element unit cells

　　

研究人员基于有限元（FE, Finite Element）轴向单元对周期性直梁进行表征。从FE组装矩阵中提取等效梁的有效刚度与惯性参数。在统一框架内对实体（Solid）与壳（Shell）单元模型进行降阶处理。准确复现了周期梁的文献基准结果。等效梁模型通过与详细FE仿真的对比验证了其预测能力。

工程中存在大量沿跨度方向呈周期性异质分布的梁结构，如夹芯构件、波纹板及带重复内部细节的轻质航空结构。尽管此类体系本质上属于三维连续体，其全局响应往往由一维梁式运动学主导。因此，构建能够保留周期微架构效应并允许高效分析的等效梁（equivalent beam）模型具有重要意义。然而，对于周期性直梁（periodic straight beams），其代表性实体不再是孤立截面，而是沿梁轴变化的轴向单元（axial cell）。由于几何、厚度及内部加劲随轴线变化，经典基于截面的表征方法不再充分，降阶本构描述必须考虑三维周期单元的本质。现有渐近均匀化（asymptotic homogenization）与计算型有限元（FE, Finite Element）单元方法虽已建立理论框架，但在复杂工程单元和标准FE工作流中的操作实现仍显繁琐。近期发展的三维轴向单元公式虽能识别剪切可变形梁本构，但多聚焦于刚度均匀化，缺乏对惯性特性的联合表征及低频模态验证。因此，研究人员亟需一种直接操作组装FE刚度与质量矩阵、统一处理实体与壳离散化、并将本构与惯性特性转移至等长梁元素的降阶工作流。

研究人员基于先前发展的基于矩阵的梁结构表征策略，将其专门化并应用于轴周期设置，提出了一种基于子结构（substructuring）的周期直梁有限元单元表征框架。该框架通过Craig-Bampton固定界面成分模态合成（fixed-interface component mode synthesis）对内部自由度进行凝聚，将两个端界面抽象为刚性手柄（handle），通过六种平衡的广义力状态识别6×6等效本构矩阵，并从约化边界质量矩阵中提取实际惯性参数。研究在三类文献基准族（带周期通孔的方孔与圆孔夹芯梁、带与不带横向加劲肋的薄壁箱形梁、非对称波纹条芯夹芯梁）以及一个sp-like周期梁和矩形穿孔实体梁上进行了系统验证。结果表明，该方法能够准确复现文献中的有效刚度系数，成功统一处理实体（Solid）与壳（Shell）单元离散化，并揭示等效梁模型在静力与低频模态预测中的能力与局限。该研究为周期微结构梁的一维高效分析提供了可直接从标准FE矩阵出发的实用化途径，对航空航天等领域的轻量化周期结构设计具有重要意义。

研究人员主要采用基于子结构的矩阵降阶与刚性手柄等效识别技术。将周期单元离散为FE模型后，通过Craig-Bampton方法凝聚内部自由度，保留两端界面；将界面运动分解为刚性手柄广义位移与局部翘曲（warping）模式，凝聚补余自由度；施加六种平衡广义载荷状态，由界面位移响应识别等效柔度矩阵并求逆得到本构矩阵；同时从约化质量矩阵投影得到单位长度质量、重心位置及单位长度惯性矩阵。该方法仅依赖于FE组装刚度与质量矩阵，统一适用于实体与壳离散化。验证基于三类文献基准（带孔夹芯梁、带/不带肋箱形梁、非对称波纹条芯梁）、一个sp-like周期壳梁及一个矩形穿孔实体梁，所有样本数据均来自公开文献基准配置。

引言与文献综述
研究人员指出，周期性直梁的经典截面表征无法处理沿轴变化的周期单元。近期基于平衡的三维轴向单元公式虽能识别耦合剪切变形梁本构，但多集中于静态刚度恢复，缺乏对惯性参数的系统提取及低频动力验证。

周期直梁的有限元表征框架
研究人员建立了直接操作组装FE刚度与质量矩阵的表征流程。将周期单元在两端界面处划分为刚体手柄变量与补余变形模式，通过静态凝聚消去内部自由度与补余界面变量，基于六种平衡广义载荷状态识别等效柔度矩阵，进而求得其逆矩阵作为等效梁本构矩阵。惯性表征则在相同约化框架下，通过刚体运动投影得到单位长度质量、截面重心位置及单位长度惯性矩阵。

基准配置与数值设置
研究人员选取三类文献基准族进行验证：带周期通孔的方孔与圆孔夹芯梁、带与不带横向加劲肋的薄壁箱形梁，以及非对称波纹条芯夹芯梁。此外，进一步考察了一个具有局部凸起特征的sp-like周期C型薄壁梁和一个带矩形通孔的实体周期梁。所有分析均基于单个周期单元开展。

文献基准验证结果
对于方孔与圆孔夹芯梁，通过单周期三维实体单元表征，所得有效抗弯刚度与文献结果高度吻合。对于薄壁箱形梁，无肋与带肋配置的轴向、扭转及弯曲刚度均与参考值一致；带肋配置正确再现了因横向加劲引起的刚度增量。壳单元与实体单元的系统对比表明，Kirchhoff-Love与Mindlin-Reissner壳公式均能复现实体单元的主导刚度项，且薄壁构型下壳模型在保持主导刚度的同时显著降低计算规模。对于非对称波纹条芯夹芯梁，所有六个对角系数及非零扭转-弯曲耦合系数D₄₆与参考结果的差异均小于2%，且网格加密后主刚度项变化不超过0.665%。

sp-like周期梁的表征
针对具有薄壁开截面和局部凸起细节的sp-like周期梁，研究人员系统研究了周期切割位置的影响。采用六个等长但切割位置不同的周期单元（包含反射对称的Cut A、Cut B及四种中间平移构型），基于壳模型建立。结果表明，轴向刚度、弯曲刚度、牵拉-弯曲耦合及剪切中心（shear center）参考的扭转刚度对切割位置高度不敏感；横向剪切块及对称性破缺产生的牵拉-剪切耦合对界面位置表现出可测量的依赖性，其中横向剪切系数D₂₂的残余敏感性最强，无法通过参考点变换消除。质量特性方面，单元总质量严格守恒，横向惯性张量分量随局部特征在有限区间内的分配而改变。

等效特性在一维分析中的应用
在静力悬臂梁对比中，研究人员将识别的本构矩阵赋予一维梁单元，通过剪切中心位置施加补偿扭矩，4周期与8周期sp-like梁的尖端位移与三维详细FE模型偏差很小。在自由-自由（free-free）模态分析中，八周期sp-like梁的前几阶整体弯扭模态与三维壳模型差异在9%以内，模态置信准则（MAC, Modal Assurance Criterion）值与展开MAC值均较高，但部分含局部壳变形的高阶模态无法找到对应的一维模态。对于十二周期矩形穿孔实体梁，前七阶整体模态中，前六阶频率差异大多在2%以内，第七阶约为3.25%，约化与展开MAC值均接近1.000，且多周期一致性检验表明将1至4个连续周期合并表征不改变本构系数，质量及轴向惯性按线性比例缩放，横向惯性包含预期的平行轴（parallel-axis）贡献。

讨论与结论
讨论指出，该框架提供了从标准FE矩阵到可直接用于降阶一维模型的等效梁属性的有效途径，其实质为一次性离线计算，识别后的属性可在多种静力、低频动力、参数研究或优化循环中重复利用。然而，等效梁模型仅适用于全局静力响应与整体梁式变形主导的低频动力行为，不能替代完整三维或壳结构谱。研究结论强调，该方法准确复现了文献基准刚度，统一了实体与壳离散化处理，成功识别了惯性参数，并在一维静力与模态分析中验证了其预测能力；未来工作将针对更完整的结构装配以及周期惯性特性在动力降阶模型中的作用展开。