基于COMSOL的复能带与凋落波研究：大、小单元嵌套声学黑洞结构PDE建模与文献复现

comsol实能带、复能带(PDE)建模 文献复现-“周期嵌套声学黑洞结构的复能带和凋落波研究”-“Complex band structure and evanescent Bloch wave propagation of periodic nested acoustic black hole phononic structure” 包括大单元、小单元、嵌套单元PDE模型文件与英文期刊PDF声学黑洞结构在减振降噪领域的潜力这几年越来越受关注最近在复现某篇关于嵌套声学黑洞复能带的论文时发现用COMSOL搞周期性嵌套结构的复能带分析比想象中更有意思。尤其是处理嵌套单元时如何在同一个模型里协调大单元和小单元的边界条件这个坑值得拿出来聊聊。先扔个模型框架上来在COMSOL的PDE模块里手动敲波动方程。比如对于二维声压场核心方程得写成这种形式% 弱形式PDE设置 damp_coeff 1e3; % 黑洞区域阻尼系数 test(p)*p_tt test(p_x)*c^2*p_x test(p_y)*c^2*p_y damp_coeff*test(p)*p_t 0重点在于黑洞区域的阻尼系数设置——这里用了随半径变化的渐变函数。论文里嵌套结构的关键在于单元内部又套娃式地排列了更小的黑洞这时候网格划分必须上尺寸函数不然小单元边缘的曲率突变直接让求解器崩掉。comsol实能带、复能带(PDE)建模 文献复现-“周期嵌套声学黑洞结构的复能带和凋落波研究”-“Complex band structure and evanescent Bloch wave propagation of periodic nested acoustic black hole phononic structure” 包括大单元、小单元、嵌套单元PDE模型文件与英文期刊PDF做复能带分析时传统实能带的Floquet周期边界已经不够用了。得在频域求解器中启用复数特征值搜索把参数扫描改成复平面上的路径积分。实际操作中发现用特征值求解器手动指定搜索路径比默认的自动扫描更稳solver.study(std1).feature(eig).set(shift, 1e4-2000i); solver.study(std1).feature(eig).set(search_radius, 5000);跑出来的复能带图里虚部超过实部的那些模式对应着论文里的强衰减区。有意思的是嵌套结构在低频段比如论文中的500Hz以下会出现能带折叠——原本单一条带突然分裂成多个枝杈这时候得用介电常数类比法判断哪些是物理解。有个取巧的方法是观察场分布真正的凋落波会在相邻单元间呈现相位反转而伪解往往伴随不合理的能量聚集。最后提一嘴模型验证的野路子把PDE模块的结果和结构力学模块做交叉对比。当黑洞壁厚趋近于零时两者在200-800Hz频段的误差居然小于3%这侧面说明手动写波动方程的可靠性。不过遇到网格敏感问题时在材料属性里加个虚刚度项比如弹性模量乘以10.01i能显著改善收敛性虽然物理上有点作弊嫌疑...