【ATC优秀论文6】基于HyperWorks的六边形蜂窝板铺层等效建模方法研究

  本文摘要（由AI生成）：

本文提出了一种基于HyperWorks的六边形蜂窝板铺层等效建模方法，可快速准确地分析蜂窝板材的力学特性。该方法不仅与原蜂窝板力学性能有较好的一致性，而且可以极大地提高计算效率。同时，基于该铺层等效有限元模型，将复杂参数与目标量化，建立蜂窝胞元尺寸长度与前五阶模态频率的函数关系，可为结构的优化设计提供依据。

《基于HyperWorks的六边形蜂窝板铺层等效建模方法研究》

徐洋  王皓辉  盛晓伟

东华大学机械工程学院

蜂窝板特点是强度高、刚度高、质量轻、隔热性能好，且在受到高强度冲击载荷时其胞壁结构发生弹性和塑性变形的同时可吸收冲击能量，是一种具有良好力学特性的轻量化结构，因此被广泛地应用于航空航天结构等各个领域。在工程领域中，对于蜂窝板这样的复杂结构，通常是借助计算机利用有限元方法对结构进行分析，其优势是能够精确地处理复杂结构问题，但缺点是需要过多的时间成本。因此，如何快速高效的建立蜂窝板模型是分析结构特性的前提条件。

在大多数工程应用中，蜂窝板是由上下蒙皮和中间的蜂窝芯子构成，蒙皮与芯子采用粘结剂进行连接。通常六边形单胞构成的周期性板是用来作为蜂窝板芯层，具体结构如图1所示。

图1 蜂窝板结构示意图

目前，蜂窝板的有限元建模方法大致可分为两种[10-11]：全实体建模方法和等效板建模方法。

（1）蜂窝板的全实体建模方法：其本质是将蜂窝板按照实体模型建模，蒙皮、芯子等全部构建，如图2所示。

图2 蜂窝板全实体建模

该方法对于单块的蜂窝板件研究可以达到不错的效果与精度，但是在实际工程应用中，这样的建模方法所耗费的计算资源与时间成本非常巨大。

（2）等效板的建模方法：是目前使用最为广泛的方法，其特点是将蜂窝板等效为一块连续的各向同性/异性的板件，具体见图3。

图3 等效板结构示意图

等效板建模虽然可以节省计算成本与时间，但是在计算精度方面随着阶数的增加误差会越来越明显。另外，使用等效板方法建模只能得到宏观的力学响应结果，对于板件内部的微观性能是无法得到的（比如蒙皮与芯子的相互作用），通常需要设计师的经验才能进行判断。

本文在前两种建模方法的基础上提出了铺层等效有限元模型建模方法。其原理是：通过材料铺层将蜂窝板等效为一种铺层性质的材料。不仅能实现层与层的共节点，而且还可以输出层与层之间的响应（如图4），同时在分析过程中将胞壁纵向受载时的伸缩变形考虑在内，使得结果更为精确。

图4 铺层等效示意图

建立如图5所示的理论模型，考虑蜂窝材料在横向受载时的材料性质，假设等效后的匀质材料为均匀的单项拉伸状态（如图6），可得：

b是壁板宽度，通过弯曲理论，AB的挠度是：

其中

为惯性矩，负载P引起壁板AB的伸长量为：

等效应变为：

等效泊松比μ1为：

通过（2）和（5）可得横向杨氏模量为：

图5 蜂窝胞元

图6 蜂窝胞元的横向变形

对于纵向受载时的力学特性分析（如图7），同样假设等效后的匀质材料为均匀的单项受力状态，与横向受载不同的是将BC段的伸缩变形考虑在内，因此可得：

由弯曲理论AB的挠度是：

负载W引起壁板AB的伸长量为：

BC的伸长量为：

      等效应变为

等效泊松比μ2为:

则纵向杨氏模量为：

图7 蜂窝胞元的纵向变形（考虑伸缩变形）

因此蜂窝材料参数为：

对于一个远小于1的量a，可以有近似的公式：

考虑蜂窝胞元壁板（BC段）的伸缩变形对于高精度工程应用领域（如精编纺织机械、航天卫星部分结构等）效果比较明显，因为这类设备通常前几阶模态以伸缩扭转为主，而对于其他大型工程应用领域（如工程机械、海洋平台等）计算结果区别并不明显。

将等效参数输入铺层等效模型中，获得的计算结果相比于用实体网格建模的全实体建模方法，该方法计算规模得到大幅缩减。而相比于同样采用Shell单元建模的等效板建模方法，该方法可使层与层之间的响应输出，准确性得到极大提高。蜂窝材料的铺层等效有限元模型构建方法能够反应板件的宏观、微观力学特性，但由于蜂窝芯层的等效层与原六边形芯层有差异，因而不适用于胞元的微观变形分析。

以某卫星的一块长方形承重板为例进行计算。文中为证明铺层等效有限元法效果，依据“全实体建模的蜂窝板与实际蜂窝板有极高的一致性”结论[1]，特将蜂窝板全实体有限元模型计算结果作为进行比较的基准[12-14]。

该板具体示意图如图8所示，长为250mm，宽为200mm，蜂窝板总厚度为15mm，上下蒙皮厚度均为0.3mm，正六边形铝制蜂窝胞元厚度为0.04mm，边长为4mm。材料弹性模量72Gpa，密度2.8×103kg/m3，泊松比为0.3。

约束条件分别如下：

在计算模态时，AB与A1B1固支；

计算静力分析时，A、B、C、D与A1、B1、C1、D1固定；

计算冲击时A、B、C、D与A1、B1、C1、D1为载荷加载点。

图8 算例示意图

经计算，可得蜂窝板固有频率计算结果如表一所示，模型前五阶振型图如图9所示。

表一 前五阶固有频率

分析图9与表一可知，三种方法获得的前五阶振型基本一致。但铺层等效法的第1、2、4、5阶误差明显要小于等效板（5%以内），第三阶模态误差有所增加，近5%。其误差来自于两个方面，一个是忽略胶粘层所带来的影响增大，另一个是由来自于壳单元与上下蒙皮夹紧系统引起差异的边界条件。此外，随着频率段的增大，误差会逐渐明显。为了更进一步体现误差，在0-5000Hz范围内将计算结果进行比较，具体如图10所示。结果显示，铺层等效模型的计算结果是更加可靠。

 a.全实体有限元模型     b.等效板有限元模型     c.铺层等效板有限元模型

图9 前五阶振型图

图10 结果对比图

选择板的中心点作为位移参考点，（载荷为Fy=7.5g、Fz=1.2g）具体结果见表二。

表二 位移计算结果

分析表二的数据可知铺层等效模型的计算结果相较于等效板精度有明显提高，最大误差为-3.1%，而等效板误差最大达到6.25%。

选择板的中心作为应力参考点，具体结果见表三。

表三 应力计算结果

分析表三的数据可知铺层等效模型的计算结果相较于等效板精度有明显提高，但是由于等效过程中边界条件的差异，绝对误差仍达到了5.3%。

冲击响应分析属于峰值类的分析，一般只关注响应的峰值是多少，所以设置的工况大多为极限工况，主要关注变化趋势，其结果相当保守。本算例冲击加速度响应谱为100Hz-1500Hz是80g，1500Hz-5000Hz是2000g，加载方向是y方向。分析比较板的最大位移轮廓以及应力分布云图，如图11、12所示。 

a.全实体有限元模型       b.等效板有限元模型     c.铺层等效板有限元模型

图11 冲击响应位移云图

a.全实体有限元模型       b.等效板有限元模型     c.铺层等效板有限元模型

图12 冲击响应应力云图

从图中分析可知面板的变形特性稍有差别，靠近夹层结构边缘处的位移偏差有所增加，原因是结构边缘主要反应为蜂窝材料胞元的局部结构特性对面板的影响。但是面板变形的总体一致性较好。应力分布基本能够相互对应，在极限载荷工况下的冲击响应，应力的大小不关注。另外，冲击响应是运用模态叠加法进行计算的，在高阶频段下模态频率的误差会越来越明显，这样的影响在冲击响应中会进一步体现[15]。

对于蜂窝板这类的多孔板材，其均匀程度对板件的力学性能有着重要的影响，需要保证足够多的胞元数量以满足均匀性要求。而胞元尺寸对胞元数量有着决定性影响[16]。目前大部分有关于蜂窝板优化方面的文献很少有将胞元尺寸作为参数变量。其原因是胞元尺寸作为参数有其特殊性，与一般的形状变量不同，胞元尺寸的每一次改变都会牵涉到蜂窝芯子的重新布局，相当于每一次迭代都需要重新建模[17-18]。现有的有限元分析软件很难做到这一点[19-20]。

基于铺层等效有限元模型，将铺层等效有限元模型的材料参数作为变量，由[9]可知

等效参数公式是线性可逆的，胞元尺寸的每一次变化都会带动铺层板材料参数的变化，由此可以间接的知道胞元尺寸与固有频率之间的关系，为今后的蜂窝板优化提供了一定的参考意义。

频率与多参数方程可以简单地表示为：

式中：h—蜂窝芯子高度；l—胞元尺寸；t—蒙皮厚度；tcell—胞元厚度；

对（21）式进行全微分得到如下方程：

（2）式可以近似地表达为：

（23）式表示蜂窝板固有频率的变化是由蜂窝板内部其中一个或多个参数变化的结果。本次仅针对胞元尺寸与固有频率的关系，因此，保持其他参数不变，而胞元尺寸从=1mm到=10.5mm每隔0.5mm取点进行计算，总共选取20个点，将计算结果导入matlab画出关系图并拟合出函数，如表四与图13。

图13 胞元厚度与频率关系图

      参数常数项的确定依赖于一定区域范围内参数与频率的关系式，对于同一范围内的数据点可能不止一种函数关系表达式，但是这些函数表达式在这一区域所表达的参数意义是相同的，因此在确定参数范围所导出的数据点以拟合的方程中只需要选择形式简便连续性好的方程即多项式形式。

表四 函数关系

根据图13的计算结果可知在胞元尺寸不断增大的情况下，第1、2阶固有频率变化并不明显，从第3阶开始，固有频率的变化逐渐增大，第5阶固有频率初始变化最明显，但是随着尺寸的增大逐渐趋于平稳。总体而言，胞元尺寸对第3、4、5阶的固有频率影响要大于对1、2阶固有频率的影响，主要原因有两方面，一方面对于有限元模型，网格划分对计算频率结果较敏感，因此不同的网格划分方法 会直接影响到计算频率。另一方面，1、2阶模态主要以弯曲变形为主，而3、4、5阶模态则以伸缩、扭转为主，由前文分析考虑了蜂窝伸缩变形使得结果更精确，所以对于以伸缩、扭转为主的3、4、5阶模态影响更大。

（1）通过铺层等效法获得的模态分析、静力分析和冲击谱分析计算结果表明，与常用有限元建模方法相比，该方法能够在保证计算效率的前提下最大限度地兼顾计算精度。

（2）蜂窝胞元尺寸对各阶模态固有频率的影响程度主要取决于网格的划分方法以及各阶模态的振型，将胞壁伸缩变形考虑在内使得结果更精确。

（3）基于该铺层等效有限元模型所建立蜂窝胞元尺寸与前5阶模态频率的函数关系，可直观地反映出胞元尺寸对前五阶固有频率的影响程度，为结构优化设计提供依据。