微结构固体中钟型与扭结孤立波的演化-力学学报.DOCVIP

微结构固体中钟型与扭结孤立波的演化1) 双山 那仁满都拉2) （内蒙古民族大学物理与电子信息学院，通辽，028043） 摘要：根据Mindlin理论和Murnaghan模型，首先建立了描述耗散、频散及非线性微结构固体中一维纵波传播的一种简单模型。然后利用有限差分方法，数值模拟了微结构效应对钟型与扭结孤立波演化的影响。结果表明，随着微结构效应的减弱，钟型孤立波的幅度衰减以及非对称特征变得越来越明显；随着微结构效应的增强，扭结孤立波顶部出现的“帽子”状变化以及由此产生的非对称特征变得越来越明显。 关键词：微结构固体，Mindlin理论，Murnaghan模型，孤立波 中图分类号： O347 文献标识码： A 文章编号： 引 言 随着工业器件的微型化，人们越来越关注微细结构材料的力学行为的研究及其性能的评价问题。近几十年来微细结构理论有了很好的发展[1-9]，并在某些方面得到了较好的应用。Engelbrecht[10]等根据Mindlin理论，建立了考虑微尺度强非线性效应的一种波模型，并分析了该模型的主要特性；Salupere[11]等用包含高阶频散项和高次非线性项的KdV类方程作为描述微结构固体中波传播模型，研究了非线性波在微结构固体中的传播特性；Janno[12]等基于Mindlin理论，建立了描述微结构固体中波传播模型，并用严格的数学理论证明了在一定条件下微结构固体中可以形成一种非对称孤立波；Porubov [13]等考虑固体材料的耗散（或能量注入）、频散及非线性效应，建立了一维应变纵波传播模型，并分析了弱耗散和弱频散效应对钟型和扭结孤立波的影响；文[14]中我们利用数值方法，详细分析了细观结构固体层中孤立波的传播及相互作用特性；文[15]中Berezovski等基于Mindlin模型建立了描述双尺度微结构固体的几种可能模型，并分析了它们的频散特性；文[16]中Casasso等向量微结构固体中波传播问题进行了研究。这些研究的重要意义在于固体材料性能的无损检测与评价。因为非线性波，特别是孤立波在微细结构固体材料中传播时，其形状、幅度以及速度中携带着反映固体材料内部结构特征的重要信息，所以对微细结构固体中孤立波形成与传播规律的研究对固体材料的无损检测与评价具有重要意义。目前已有些这方面的工作，如文[17]中提出了利用微结构固体中形成的非对称孤立波可以确定固体材料结构参数的论点，并在理论上给出了严格的证明。文[18]中我们利用一种简单模型验证了微结构频散固体中孤立波形成与传播的可能性，并给出了利用所形成的孤立波来确定固体材料结构参数的可能方案。 但目前提出的微结构模型通常包含很多参数而且这些参数的值都很难甚至不可能确定。到目前很少有文献报导关于微结构参数数据[19]。数据的缺少，使得这些微结构模型在实际应用中受到了很大的限制。因此对于微结构固体建立起简单实用的理论模型是迫切而重要的工作。本文首先根据Mindlin理论和Murnaghan模型，建立耗散、频散及非线性固体中一维纵波传播的一种简单模型。然后利用有限差分方法，模拟钟型与扭结孤立波的演化过程，给出反映孤立波演化特征的模拟图像并分析微结构效应对钟型与扭结孤立波演化的影响。 1 模型方程的建立 根据Mindlin理论[1]，在忽略体力情况下，微结构固体一维运动方程为 （1） （2） 这里为宏观密度，为微惯性，为宏观位移，为微形变，，和分别为Cauchy应力，相对应力和二重应力。在一维纵波传播情况下，，和的平衡部分分别由下式给出 ， ， （3） 这里为宏观应变，为相对形变，为微形变梯度，而为自由能，是，和的函数。这里我们考虑固体材料微结构效应足够弱的情况，因此根据Mindlin理论，相对形变线性化地取为，，微形变梯度线性化地取为，并可忽略微尺度非线性效应。另外，我们采用Murnaghan模型[20]，只考虑物理非线性而忽略几何非线性，故宏观应变线性化为。此时，利用文[1]中给出的自由能表达式（5.5）并考虑到宏观应变的三次项，则可有 （4） 这里,为Lame系数，，为Murnaghan三阶弹性常数，和为微结构有关的参数。此式类似于文[12]给出的自由能表达式。由式（4）和（3）计算得三种应力的平衡部分分别为 ， ， （5） 这里。由于固体材料的微结构效应足够弱，我们可忽略相对形变和微形变有关的耗散效应，只考虑宏观应变有关的耗散效应。此时宏观应力的考虑耗散影响的非平衡部分为（为常数）。因此，三种总应力分别为 ， ， （6） 把（6）式，代入方程（1