本发明涉及蜂窝板建模技术领域，特别是涉及一种蜂窝材料的铺层等效有限元模型构建方法。
背景技术：
蜂窝板的特点是强度高、刚度高、质量轻、隔热性能好，且在受到高强度冲击载荷时其胞壁结构发生弹性和塑性变形的同时可吸收冲击能量，是一种具有良好力学特性的轻量化结构，因此被广泛地应用于航空航天结构等各个领域。在工程领域中，对于蜂窝板这样的复杂结构，通常是借助计算机利用有限元方法对结构进行分析，其优势是能够精确地处理复杂结构问题，但缺点是需要过多的时间成本。因此，如何快速高效的建立蜂窝板模型是分析结构特性的前提条件。
在大多数工程应用中，蜂窝板是由上下蒙皮和中间的蜂窝芯子构成，蒙皮与芯子采用粘结剂进行连接，通常六边形单胞构成的周期性板是用来作为蜂窝板芯层。
目前，蜂窝板的有限元建模方法大致可分为两种：全实体建模方法和等效板建模方法。
(1)蜂窝板的全实体建模方法：其本质是将蜂窝板按照实体模型建模，蒙皮、芯子等全部构建。
该方法对于单块的蜂窝板件研究可以达到不错的效果与精度，但是在实际工程应用中，这样的建模方法所耗费的计算资源与时间成本非常巨大。
(2)等效板的建模方法：是目前使用最为广泛的方法，其特点是将蜂窝板等效为一块连续的各向同性/异性的板件。
等效板建模虽然可以节省计算成本与时间，但是在计算精度方面随着阶数的增加误差会越来越明显。另外，使用等效板方法建模只能得到宏观的力学响应结果，对于板件内部的微观性能是无法得到的(比如蒙皮与芯子的相互作用)，通常需要设计师的经验才能进行判断。
技术实现要素：
本发明所要解决的技术问题是提供一种蜂窝材料的铺层等效有限元模型构建方法，不仅能实现层与层的共节点，而且还可以输出层与层之间的响应，同时在分析过程中将胞壁纵向受载时的伸缩变形考虑在内，使得结果更为精确。
本发明解决其技术问题所采用的技术方案是：提供一种蜂窝材料的铺层等效有限元模型构建方法，通过材料辅层将蜂窝材料等效为一种辅层性质的材料，其中，所述材料辅层的等效参数包括：对于蜂窝材料在横向受载时的力学特性分析后得到的横向等效泊松比和横向杨氏模量，以及对于蜂窝材料在纵向受载时的力学特性分析后得到的纵向等效泊松比和纵向杨氏模量；将材料辅层的等效参数输入有限元模型中，逐层辅层，构建蜂窝材料的铺层等效有限元模型。
所述横向等效泊松比
所述蜂窝材料的蜂窝胞元为正六边形的蜂窝胞元，对于远小于1的量，假设其为量a，则有等效公式
所述胞元斜边ab与所述胞元竖边bc为蜂窝胞元的一组相邻壁板。
所述胞元竖边bc垂直于蜂窝材料在横向的负载。
所述胞元竖边bc平行于蜂窝材料在纵向的负载。
所述的蜂窝材料的铺层等效有限元模型构建方法还包括将材料辅层的等效参数作为变量，得到蜂窝材料的胞元尺寸与固有频率的函数关系的步骤。
所述得到蜂窝材料的胞元尺寸与固有频率的函数关系具体为：将蜂窝材料的固有频率与多个参数的方程简化表示为：fi＝f(h,l,m,tcell)，其中，h为胞元高度，l为胞元边长，m为蒙皮厚度，tcell为胞元厚度，对方程进行全微分得到
所述的蜂窝材料的铺层等效有限元模型构建方法还包括将得到的蜂窝材料的铺层等效有限元模型，与蜂窝材料的全实体有限元模型进行对比分析，判断建模精度的步骤。
所述对比分析至少包括：固有频率比对、静力分析和冲击谱的计算。
有益效果
由于采用了上述的技术方案，本发明与现有技术相比，具有以下的优点和积极效果：本发明可快速准确分析蜂窝板材力学特性。基于该模型所计算的模态、静力和冲击谱结果与目前常用有限元建模方法计算结果进行对比，可以得出结论：铺层等效有限元模型不仅与原蜂窝板力学性能有较好的一致性，而且可以极大地提高计算效率。同时，基于该铺层等效有限元模型，将复杂参数与目标量化，建立蜂窝胞元尺寸长度与前五阶模态频率的函数关系，可为结构的优化设计提供依据。
附图说明
图1为本发明方法的具体实施方式的辅层等效示意图。
图2为本发明方法的具体实施方式的蜂窝胞元的示意图。
图3和图4为本发明方法的具体实施方式的蜂窝胞元的横向变形的受力分析图。
图5～图7为本发明方法的具体实施方式的蜂窝胞元的纵向变形(考虑伸缩变形)的受力分析图。
图8为本发明方法的具体实施方式中与蜂窝板全实体有限元模型的对比算例的建模示意图。
图9为本发明方法的具体实施方式中与蜂窝板全实体有限元模型的对比算例的固有频率对比图。
图10为本发明方法的具体实施方式的胞元尺寸与固有频率的关系图。
具体实施方式
下面结合具体实施例，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解，在阅读了本发明讲授的内容之后，本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。
迄今为止，已经有很多学者研究了蜂窝板等效理论，如：a.boudjemai、r.amri等人利用了各向同性的整块板等效方法，方法虽可行实用，但却只能从宏观上描述其特性，无法从微观上进行有效分析——即无法反应蜂窝板蒙皮与蜂窝芯子之间的影响。董彦鹏、吕振华等人提出了一种结构相似有限元模型等效方法，等效性得到了验证，但是建模过程仍然非常复杂，并没有得到有效简化。蔡园武、徐亮等人提出了一种周期性板件刚度等效方法，但是其计算结果误差相对较大，只能在相对狭小的特定范围内才能使用。allen等人建立的蜂窝板模型则是忽略了面内刚度和弯曲刚度，虽然蜂窝芯层很软，但是对于蒙皮而言具有较大的厚度，所以忽略面内刚度和弯曲刚度必然导致明显的误差。ha与崔光育等人通过将蜂窝芯层视为服从剪切变形理论的正交异性层建立了分析模型，在部分领域内效果很好。gisbon等人通过理论方法给出了蜂窝材料的等效参数公式，虽然其解析形式便于应用，但是也有缺点，在推导过程中仅仅只考虑了蜂窝板的弯曲变形，而没有考虑蜂窝胞元壁板的伸缩变形，因此其公式带有一定的误差。夏利娟等人将蜂窝板各种等效方法进行计算并将计算结果进行了比较，但是却没有具体介绍等效方法。
为了快速准确分析蜂窝板材力学特性，本发明提出了一种基于hypermesh的铺层等效有限元建模方法，其计算结果与目前常用有限元建模方法比较可知，该方法具有更高的效率且准确可靠。
本发明的蜂窝材料的铺层等效有限元模型建模方法，其主要原理是：通过材料铺层将蜂窝板等效为一种铺层性质的材料。不仅能实现层与层的共节点，而且还可以输出层与层之间的响应，如图1，同时在分析过程中将胞壁纵向受载时的伸缩变形考虑在内，使得结果更为精确。其中，在图1中，11为上蒙皮层，12为上胶粘层、13为蜂窝芯层，14为下胶粘层，15为下蒙皮层，箭头方向a为胞元辅层方向。
建立如图2所示的理论模型，考虑蜂窝材料在横向受载时的材料性质，假设等效后的匀质材料为均匀的单项拉伸状态，如图3和图4所示，可得
p＝σ1(h+lsinθ)b(2)
其中，b是蜂窝材料的蜂窝胞元的壁板宽度，p是蜂窝材料在横向的负载，m是蜂窝胞元的胞元斜边对载荷点的弯矩，l是胞元斜边的长度，θ是胞元斜边与水平方向夹角，h是胞元竖边的长度，σ1是胞元在均匀的单项拉伸状态下受到的应力(方向为横向)；
通过弯曲理论，壁板ab的挠度是：
其中，i为惯性矩，es是蜂窝材料的剪切模量，t是蜂窝胞元的壁板厚度。
负载p引起壁板ab的伸长量为：
等效应变为：
其中，δ1是横向等效应变，δ2是纵向等效应变。
横向等效泊松比μ1为：
通过式(2)和式(5)可得横向杨氏模量为：
对于纵向受载时的力学特性分析，如图5～7所示，同样假设等效后的匀质材料为均匀的单项受力状态，与横向受载不同的是，将壁板bc段的伸缩变形考虑在内，因此可得：
w＝lσ2sinθ(10)
其中，w是蜂窝材料在纵向的负载，σ2是胞元在均匀的单项拉伸状态下受到的应力(方向为纵向)
根据弯曲理论，壁板ab的扰度是：
负载w引起壁板ab的伸长量为：
壁板bc的伸长量为：
等效应变为：
纵向等效泊松比μ2为：
则纵向杨氏模量为：
因此，蜂窝材料参数为：
进一步地：
对于一个远小于1的量a，可以有近似的公式
在公式中可以注意到t2/l2是很小的，忽略其高阶无穷小量。对于正六边形的蜂窝胞元h＝l、θ＝π/6，则式(18)简化为：
考虑蜂窝胞元壁板(bc段)的伸缩变形对于高精度工程应用领域(如精编纺织机械、航天卫星部分结构等)效果比较明显，因为这类设备通常前几阶模态以伸缩扭转为主，而对于其他大型工程应用领域(如工程机械、海洋平台等)计算结果区别并不明显。
将等效参数输入铺层等效模型中，获得的计算结果相比于用实体网格建模的全实体建模方法，该方法计算规模得到大幅缩减。而相比于同样采用shell单元建模的等效板建模方法，该方法可使层与层之间的响应输出，准确性得到极大提高。蜂窝材料的铺层等效有限元模型构建方法能够反应板件的宏观、微观力学特性，但由于蜂窝芯层的等效层与原六边形芯层有差异，因而不适用于胞元的微观变形分析。
下面以某卫星的一块长方形承重板为例进行计算。为证明铺层等效有限元法效果，依据“全实体建模的蜂窝板与实际蜂窝板有极高的一致性”结论，特将蜂窝板全实体有限元模型计算结果作为进行比较的基准。
该板具体示意图如图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为载荷加载点。
一、固有频率计算结果
经计算，可得蜂窝板固有频率计算结果如表1所示。
表1.前五阶固有频率
分析表1和前五阶段对应的阵型图可知，三种方法获得的前五阶振型基本一致。但铺层等效法的第1、2、4、5阶误差明显要小于等效板(5％以内)，第三阶模态误差有所增加，近5％。其误差来自于两个方面，一个是忽略胶粘层所带来的影响增大，另一个是由来自于壳单元与上下蒙皮夹紧系统引起差异的边界条件。此外，随着频率段的增大，误差会逐渐明显。为了更进一步体现误差，在0-5000hz范围内将计算结果进行比较，具体如图9所示。结果显示，铺层等效模型的计算结果是更加可靠。
二、静力计算结果
(1)位移计算
选择板的中心点作为位移参考点，(载荷为fy＝7.5g、fz＝1.2g)具体结果见表2。
表2.位移计算结果
分析表二的数据可知铺层等效模型的计算结果相较于等效板精度有明显提高，最大误差为-3.1％，而等效板误差最大达到6.25％。
(2)应力计算
选择板的中心作为应力参考点，具体结果见表3。
表3.应力计算结果
分析表3的数据可知铺层等效模型的计算结果相较于等效板精度有明显提高，但是由于等效过程中边界条件的差异，绝对误差仍达到了5.3％。
三、冲击谱计算结果
冲击响应分析属于峰值类的分析，一般只关注响应的峰值是多少，所以设置的工况大多为极限工况，主要关注变化趋势，其结果相当保守。本算例冲击加速度响应谱为100hz-1500hz是80g，1500hz-5000hz是2000g，加载方向是y方向，分析比较板的最大位移轮廓以及应力分布云图。
分析可知面板的变形特性稍有差别，靠近夹层结构边缘处的位移偏差有所增加，原因是结构边缘主要反应为蜂窝材料胞元的局部结构特性对面板的影响。但是面板变形的总体一致性较好。应力分布基本能够相互对应，在极限载荷工况下的冲击响应，应力的大小不关注。另外，冲击响应是运用模态叠加法进行计算的，在高阶频段下模态频率的误差会越来越明显，这样的影响在冲击响应中会进一步体现。
基于本发明的蜂窝材料的铺层等效有限元模型构建方法还可以得到胞元尺寸与固有频率的函数关系。
对于蜂窝板这类的多孔板材，其均匀程度对板件的力学性能有着重要的影响，需要保证足够多的胞元数量以满足均匀性要求。而胞元尺寸对胞元数量有着决定性影响。目前大部分有关于蜂窝板优化方面的文献很少有将胞元尺寸作为参数变量。其原因是胞元尺寸作为参数有其特殊性，与一般的形状变量不同，胞元尺寸的每一次改变都会牵涉到蜂窝芯子的重新布局，相当于每一次迭代都需要重新建模。现有的有限元分析软件很难做到这一点。
基于铺层等效有限元模型，将铺层等效有限元模型的材料参数作为变量，而等效参数公式是线性可逆的，胞元尺寸的每一次变化都会带动铺层板材料参数的变化，由此可以间接的知道胞元尺寸与固有频率之间的关系，为今后的蜂窝板优化提供了一定的参考意义。
具体地：
蜂窝材料的固有频率与多参数方程可以简单地表示为：
fi＝f(h,l,m,tcell)，i＝1,2,3,4,5(21)
其中，h为胞元高度，l为胞元边长，m为蒙皮厚度，tcell为胞元厚度，对式(21)进行全微分得到：
式(22)可以近似地表达为：
式(23)表示蜂窝板固有频率的变化是由蜂窝板内部其中一个或多个参数变化的结果。本次仅针对胞元尺寸与固有频率的关系，因此，保持其他参数不变，胞元尺寸从l＝1mm到l＝10.5mm间隔0.5mm取点进行计算，总共选取20个点，将计算结果导入matlab软件，绘制关系图并拟合出胞元尺寸与固有频率的函数关系，如表4和图10所示。
表4.函数关系表
其中，参数常数项的确定依赖于一定区域范围内参数与频率的关系式，对于同一范围内的数据点可能不止一种函数关系表达式，但是这些函数表达式在这一区域所表达的参数意义是相同的，因此在确定参数范围所导出的数据点以拟合的方程中只需要选择形式简便连续性好的方程即多项式形式。
根据图10的计算结果可知在胞元尺寸增大的情况下，第1、2阶固有频率变化并不明显，从第3阶开始，固有频率的变化逐渐增大，第5阶固有频率初始变化最明显，但是随着尺寸的增大逐渐趋于平稳。总体而言，胞元尺寸对第3、4、5阶的固有频率影响要大于对1、2阶固有频率的影响，主要原因有两方面，一方面对于有限元模型，网格划分对计算频率结果较敏感，因此不同的网格划分方法会直接影响到计算频率。另一方面，1、2阶模态主要以弯曲变形为主，而3、4、5阶模态则以伸缩、扭转为主，由前文分析考虑了蜂窝伸缩变形使得结果更精确，所以对于以伸缩、扭转为主的3、4、5阶模态影响更大。
不难发现，本发明提出的一种蜂窝材料的铺层等效有限元模型构建方法，可快速准确分析蜂窝板材力学特性，其计算结果与目前常用有限元建模方法比较可知，该方法具有更高的效率且准确可靠，具体体现在：
(1)通过铺层等效法获得的模态分析、静力分析和冲击谱分析计算结果表明，与常用有限元建模方法相比，该方法能够在保证计算效率的前提下最大限度地兼顾计算精度；
(2)蜂窝胞元尺寸对各阶模态固有频率的影响程度主要取决于网格的划分方法以及各阶模态的振型，将胞壁伸缩变形考虑在内使得结果更精确；
(3)基于该铺层等效有限元模型所建立蜂窝胞元尺寸与前5阶模态频率的函数关系，可直观地反映出胞元尺寸对前五阶固有频率的影响程度，为结构优化设计提供依据。