COMSOL—弹性波 模型介绍： 激励信号为高斯子波和雷克子波，雷克子波通过时间分布，高斯子波通过x,y俩个变元在源空间进行定义。 本模型采用弹性波，时域显示物理场，中心频率为170kHz，通过体载荷来添加激励信号。

弹性波仿真里激励信号的设定总能让刚入坑的同学犯迷糊。今天咱们来拆解一个典型的COMSOL弹性波模型，重点看看高斯子波和雷克子波这两种激励源怎么玩转。模型参数挺有意思，中心频率170kHz，相当于超声波检测常用的频段，搞无损检测的朋友应该很眼熟。

先说说雷克子波的实现。这玩意儿在时域里通过体载荷加载，核心在于时间函数的表达式要写对。下面这段代码看着简单，但参数搞错分分钟出鬼：

ricker_wave = (1 - 2*(t - t0)**2 * fc**2) * exp(-(t - t0)**2 * fc**2)

这里t0=3/(2*fc)是个关键参数，相当于时间偏移量。新手常犯的错是直接照搬网上代码不调这个偏移量，结果波形直接怼在时间零点，导致仿真初始阶段出现异常震荡。记得把脉冲突起位置挪到时间轴中间区域，就像煎饼果子得摊匀了才好吃。

转到高斯子波这边，玩法就换套路了。这货得在空间维度上做文章，x和y两个方向都要加料。看这段空间分布的定义：

% 高斯空间分布
gauss_xy = exp(-((x-x0)/sigma_x)^2 - ((y-y0)/sigma_y)^2);

这里藏着三个玄机：中心点坐标(x0,y0)、横向扩散系数sigmax和纵向sigmay。搞地震波模拟的兄弟应该知道，sigma_x调太小会导致声源像个尖锥，能量过于集中容易数值发散；调太大又像摊大饼，分辨率扑街。建议初始值取波长的一半，比如170kHz对应波长约8.9mm（假设材料是铝），sigma设4-5mm比较稳妥。

COMSOL—弹性波 模型介绍： 激励信号为高斯子波和雷克子波，雷克子波通过时间分布，高斯子波通过x,y俩个变元在源空间进行定义。 本模型采用弹性波，时域显示物理场，中心频率为170kHz，通过体载荷来添加激励信号。

实际在COMSOL里组装体载荷时，得把时间和空间分量耦合起来。有个骚操作是用条件判断语句实现动态加载：

// 体载荷表达式
(body_load, (t < pulse_width) ? ricker_wave * gauss_xy : 0)

这个三目运算符相当于给激励源加了个时间闸门，避免持续加载导致边界反射干扰。注意脉冲宽度pulse_width别设太短，至少包含主瓣和两个旁瓣，否则就跟用菜刀切牛排似的——使不上劲。

调试时建议先跑个简化模型：把材料设成各向同性弹性体，网格划得粗一点，时间步长放大到1e-7秒量级。重点看波前传播是否呈同心圆状（二维情况）——要是出现棱角分明的多边形，八成是空间离散参数没调好或者激励信号频带太宽。这时候该掏出频域分析工具，看看激励信号的傅里叶变换是不是在170kHz附近有尖峰，别整成宽带噪声了。

最后唠叨句参数联动：当改变中心频率fc时，记得同步调整雷克子波的t0和空间高斯分布的sigma。这就像吃重庆火锅，毛肚涮的时间和辣度得匹配，参数之间不协调的话，轻则波形畸变，重则数值爆炸。下次遇到仿真结果像抽象画的时候，先别砸键盘，检查下这几个参数的勾稽关系准没错。