等效偶应力模量计算的代表体元模型.docVIP

基于偶应力理论的格栅材料等效介质模型 苏文政，刘书田 大连理工大学工程力学系工业装备结构分析国家重点实验室，116024 刘书田: stliu@；苏文政：wzhsu@126.com 摘要： 格栅材料的均质化等效介质模型通常是基于经典连续介质（Cauchy介质）理论建立的，而经典连续介质理论的前提是材料微结构尺寸相对宏观结构尺寸很小。实际工程结构中常使用的薄壁结构（板、壳、梁）某些方向的尺度与微结构尺度同量级或接近，此时Cauchy介质模型可能具有局限性。探讨格栅材料的高阶介质等效模型具有重要意义。本文基于偶应力理论，对格栅材料的弹性性能进行等效描述，建议了基于偶应力理论的格栅材料等效介质模型以及确定等效材料性质参数的代表体元模型，导出了计算各分量相应的位移边界条件。根据实际结构中广泛存在的格栅梁结构特点，导出了偶应力介质梁的抗弯刚度表达式，定义了偶应力影响因子δ以表征梁的偶应力效应；在此基础上具体计算了几种典型的格栅材料的等效偶应力模量，以及相应的梁结构的抗弯刚度，计算了格栅梁在一定边界条件和载荷作用下的挠曲线，并采用有限元离散模型分析了相应格栅材料结构在相同边界条件和载荷下的响应，将二者进行对比，结果表明，本文的等效结果具有较高精度。通过对不同厚度尺寸和不同几何构型的格栅梁的偶应力影响因子的分析，表明当结构的宏观尺寸和微结构尺寸相差不大时，梁表现出强烈的偶应力效应，而当宏观结构尺寸远大于微结构尺寸时，偶应力效应可以忽略，偶应力介质退化为经典连介质。 关键词：格栅材料，格栅梁，偶应力理论，等效模量，代表体元法 引言 轻质多孔固体以其低的相对密度、高的比刚度和比强度、相对低的热导率以及高压缩比等特点而作为轻质结构件、隔热体和吸能材料广泛应用于航空航天以及工业生产的其它领域[1]。其中的周期性二维格栅材料又称为蜂窝材料，由于其制造工艺相对成熟，结构规则价格低廉而得到更加广泛的应用。为充分发挥这些优异性质，需要研究格栅材料的性能以及结构的分析和设计技术。 格栅材料结构是由尺度远比结构尺度小的基本结构（单胞）周期性重复构成的，结构包含大量的基本结构单元。将格栅材料结构作为格栅结构，利用通常的有限元技术分析结构响应，将会因变量太多、计算工作量太大而不可行。因此，将格栅材料连续化为一种连续介质，通过分析连续介质的响应以确定格栅材料的响应，是有效的处理方法。等效介质的响应分析需要获得格栅材料的等效材料性质。关于格栅材料等效弹性性能表征方面的研究已开展了大量的工作。目前，大部分工作是将格栅材料等效连续化为经典连续介质（Cauchy介质），主要研究等效介质的杨氏模量和泊松比等[1][2]。然而，经典连续介质假定材料任意一点没有微结构，其变形可由3个位移变量完全描述，任何一个微元表面只传递力而不传递力矩。换言之，只有当宏观结构的最小尺寸远大于材料微结构尺寸时，非均质材料才可以等效为均匀的经典连续介质[3]。实际工程结构中的格栅材料结构件往往并不能满足这一假定，例如格栅梁结构，虽然其长度方向的尺寸远大于材料微结构尺寸，但是其厚度方向的尺寸与微结构尺寸相差并不大，对于这类材料，基于经典连续介质力学假设的均匀化方法以及代表体元法不能够准确的描述材料的等效性能。均质化的材料必须考虑微观结构的尺度特点，即包含微结构的高阶连续介质（high order continuum）[4]。 Mindlin[5]提出的偶应力理论考虑了材料微结构的影响，认为微元表面既能够传递力也能够传递力偶，成功的应用于圆孔应力集中问题。Ostoja-Starzewski等[6~8]计算了颗粒增强复合材料的等效偶应力模量，但是作者没有对等效结果进行验证。文献[9]简单计算了方形格栅材料的等效偶应力性能，认为偶应力就是截面上力偶的简单平均，显然会产生较大误差。本文基于偶应力理论研究了格栅材料的均质化问题，建立了确定本构关系各模量的代表体元边界条件，根据等效模量对格栅梁结构进行载荷—响应分析，并与实际的格栅结构的有限元离散解进行了比较。在此基础上定量分析了格栅梁结构问题中偶应力效应的影响。 偶应力理论简介 以微极理论为代表的极性理论认为材料微元的运动既包含微元体质心的平动又包含微元体内任意点绕质心的独立的刚性转动。如果假设材料点的微转动不是独立的，而是等同于传统弹性力学中一点邻域的刚性转动，即得到偶应力理论又称约束转动的微极理论[10]。平面偶应力理论中的微元体的应力状态可由四个应力分量σx、σy、τxy、τyx以及两个偶应力分量mxz、myz完全确定，如下图所示。 Figure 图 1 应力和偶应力分量指向 其变形方程为： 值得注意的是偶应力理论中，转角φ并不是独立的，即： 不计体力与体力偶，其平衡方程为： 由式可知，在偶应力理论中，剪应力互等关系并不成立（即τxy≠