[论文解读] Wave propagation and homogenization in 2D and 3D lattices: a semi-analytical approach

Anton A. Kutsenko, Adam J. Nagy+3 more|arXiv (Cornell University)|Mar 10, 2017

Acoustic Wave Phenomena Research参考文献 16被引用 8

一句话总结

本论文提出了一种半解析方法，用于通过支持纵向和弯曲波的梁单元，对二维和三维周期性晶格中的波传播与均质化进行建模。通过施加Bloch-Floquet条件，该方法利用动态刚度矩阵将问题简化为广义厄米特特征值问题，从而实现准静态极限下色散关系和Christoffel矩阵的解析推导。该方法在六方和四面体晶格中准确预测了能带结构和单波行为，且与FEM模拟结果高度一致。

ABSTRACT

Wave motion in two- and three-dimensional periodic lattices of beam members supporting longitudinal and flexural waves is considered. An analytic method for solving the Bloch wave spectrum is developed, characterized by a generalized eigenvalue equation obtained by enforcing the Floquet condition. The dynamic stiffness matrix is shown to be explicitly Hermitian and to admit positive eigenvalues. Lattices with hexagonal, rectangular, tetrahedral and cubic unit cells are analyzed. The semi-analytical method can be asymptotically expanded for low frequency yielding explicit forms for the Christoffel matrix describing wave motion in the quasistatic limit.

研究动机与目标

• 开发一种高精度与计算效率的半解析框架，用于周期性二维和三维梁晶格中的波传播建模。
• 通过基于动态刚度矩阵的广义特征值问题求解，推导出矩形、六方和立方晶格的显式色散关系。
• 通过闭式Christoffel矩阵提取准静态极限行为，实现有效材料属性的精确均质化。
• 针对复杂晶格几何结构（包括六方和四面体晶胞）的FEM模拟结果，对模型进行验证。
• 展示对称晶格中单波（声学类）行为的出现，及其在声子晶体与五模材料设计中的潜力。

提出的方法

• 采用Timoshenko梁理论，对晶格中梁单元的纵向和弯曲波模态进行建模。
• 显式推导出动态刚度矩阵，以关联端部力与位移，确保该矩阵为厄米特矩阵并具有正特征值。
• 在晶胞上施加Bloch-Floquet周期性条件，将问题简化为关于波频率与波数的广义特征值问题。
• 在低频渐近极限下，通过解析推导准静态极限下的Christoffel矩阵，实现均质化。
• 通过在离散化的k空间网格上寻找系统矩阵的最小特征值，获得数值解，收敛性由容差阈值控制。
• 利用COMSOL对FEM模拟结果进行验证，直接对比多种晶格类型的色散曲线。

实验结果

研究问题

• RQ1 是否能够通过半解析方法在不依赖完整FEM模拟的前提下，准确预测二维和三维梁晶格的完整色散谱？
• RQ2 弯曲波贡献如何影响晶格结构的准静态均质化，特别是对有效质量的影响？
• RQ3 晶格对称性在实现单波行为与准声波传播中起何种作用？
• RQ4 动态刚度矩阵形式化在三维晶格（如四面体与立方结构）中的适用范围在多大程度上可被扩展？
• RQ5 能否在低频极限下，为不同晶格对称性推导出Christoffel矩阵的显式解析表达式？

主要发现

• 半解析模型准确再现了立方晶格的前八阶色散分支，与FEM模拟结果高度一致，尤其在传播模态与共振模态方面表现优异。
• 在立方晶格中，前两个共振频率分别位于51.949 kHz与92.354 kHz，其在色散图中表现为平坦分支，模型对其预测效果良好。
• 对于四面体晶格，在5至20 kHz范围内存在一个宽广的单波区域，由于高对称性，波几乎非色散且各向同性，与有效标量波行为一致。
• 在准静态极限下，矩形、六方和立方晶格的Christoffel矩阵均以闭式表达式推导得出，实现了对有效刚度与质量属性的精确均质化。
• 研究表明，忽略弯曲波效应将导致静态极限下有效质量的低估，证实了在准确均质化中必须包含弯曲模态。
• 该方法成功拓展至三维晶格，首次基于此梁基半解析框架，对六方与四面体晶胞晶格实现了解析处理。