COMSOL三维弹性波能带计算实战技巧

1. 弹性波三维能带计算实战指南

在周期性结构波动控制领域，三维弹性波能带计算一直是研究者们又爱又恨的课题。最近我在复现Smart Materials and Structures期刊2015年那篇经典论文《Three dimensional periodic foundations for base seismic isolation》时，深刻体会到了这一点。经过两周的反复调试，终于用COMSOL 5.4成功构建了钢/橡胶复合周期性基础的完整能带计算模型。下面就把我的实战经验完整分享给大家，特别是那些容易踩坑的关键设置。

2. 模型构建基础

2.1 物理场选择与初始化

新手最容易犯的第一个错误就是物理场选择。在COMSOL 5.4中新建模型时，一定要选择"声学模块"下的"压力声学，频域"接口，而不是直觉上更接近的"固体力学"模块。这是因为：

1. 声学模块内置了Bloch周期边界条件的完整支持
2. 固体力学模块默认不考虑相位延迟边界条件
3. 弹性波方程在声学模块中可以通过压力-位移耦合形式完美表达

正确的初始化命令应该是：

matlab复制model = ModelUtil.create('Model');
model.modelNode.create('comp1');
model.physics.create('waeq1', 'WaveEquation', 'geom1');

2.2 材料参数设置

论文中的钢/橡胶复合材料需要特别注意橡胶的材料模型。普通弹性模型无法准确描述橡胶的阻尼特性，必须使用复弹性模量：

matlab复制% 钢材参数
model.material.create('steel');
model.material('steel').propertyGroup('def').set('youngs_modulus', '200e9');
model.material('steel').propertyGroup('def').set('poissons_ratio', '0.3');
model.material('steel').propertyGroup('def').set('density', '7850');

% 橡胶参数（含3%损耗因子）
model.materi