企业官网建设流程全解析

做铝合金板材的无损检测或者结构健康监测，想在正式投样做实验之前，先在上游把"一发一收"的信号链路摸清楚，这是很多项目启动时最实际的需求。压电圆片贴在铝板上，一片负责激励、一片负责接收，中间传播的板波就是Lamb波。这个场景用COMSOL做全波仿真，这些年几乎是业内标准动作：一方面要回答"能不能收到信号"这种可行性问题，另一方面要为后面的传感器布局、频率选择、缺陷定位算法提供数据支撑。

这篇博文不打算把某个现成的mph文件扔出来就完事，而是把从建模思路到参数选择的完整链条摊开讲。标题里的几个关键词——160kHz圆片、一发一收、铝板Lamb波——每一个都对应具体的物理问题和仿真决策。我会把激励频率为什么取160kHz、模型尺寸怎么定、网格和时间步怎么协同、接收信号怎么从压电方程里提出来，以及跑完之后怎么验证结果这些环节都串起来讲一遍。适合正在做压电超声仿真、结构健康监测方向，或者刚接触COMSOL瞬态压电仿真的朋友参考。

1. 项目要做什么：160kHz一发一收Lamb波检测的仿真逻辑

1.1 一发一收到底在测什么

先把这个实验布局说清楚。所谓一发一收，就是在同一块铝板的表面布置两个压电圆片，中间隔着一段距离。一个片接激励电压，另一个片只做接收。激励片在厚度方向加上交变电压后，由于逆压电效应产生厚度方向的伸缩振动，振动通过压电片的底面传入铝板，在板内激发出Lamb波。波沿着板面传播到接收片下方的时候，板表面的应变施加到接收片上，正压电效应把机械应变转换成电荷，电荷在电极上建立起电压。这一整条链路，从电压到机械波再到电压，中间经过了两次压电转换和一段板内传播。

为什么要用一发一收而不是单探头自发自收？模拟几次实验就能体会。自发自收配置里，激励脉冲还没结束，反射回波就到了，近场区域里发射信号和接收信号混在一起，很难切开。一发一收的几何布局把激励和接收在空间上分开，激励脉冲的电磁串扰和直达波在时间上可以隔开，信号干净得多。而且一发一收对应的是结构健康监测里最常见的"主动传感"方式——一条声路径一次测量，后续做损伤定位时需要多对这样的路径交叉覆盖。所以这个仿真做出来，不止是验证一个探头能不能工作，而是给后续路径分析和阵列设计打个底。

这个实验布局里有一个隐蔽但决定成败的因素：接收片测到的电压是多少，取决于接收片下表面这片铝板的应变状态。换句话说，仿真里建的"接收电路"和你设置的探针位置，会直接决定你看到的A扫信号是几毫伏还是几十毫伏。很多初学者建完模型后遇到"接收电压为零"或"信号特别小"，问题往往出在电极和边界条件上，而不是波没传过去。

1.2 为什么是160kHz和铝板

Lamb波和体波不一样，它在板内传播时存在频散现象，而且随着频厚积（频率×板厚）的变化，传播模式的数量和特性都在变。对3mm厚的铝板来说，160kHz对应的频厚积是480kHz·mm，在这个区间里，最低阶的对称模式S0和反对称模式A0都可以被激发出来，二阶以上的高阶模式还没有出现，信号相对干净，适合用来做模式识别和声速标定。

频率选低了或者选高了都有问题。频率再低，波长变大，对缺陷的灵敏度下降；频率再高，高阶模式大量出现，接收信号变成一堆模式的叠加，想分离出单一声路径信息就很麻烦。160kHz这个量级在铝板检测里属于"中低频"区间，波长大概在十几毫米到三十几毫米这个范围，既不会在3mm板里迅速衰减，又能保持对厘米级缺陷的响应能力。铝材的声阻抗比钢材低，波在板内传播衰减相对小，加上铝板在航空、轨道交通、船舶领域应用极广，这个材料选择本身就是贴近工程实际的。

另外一个隐含的问题是：标题里说的"160kHz圆片传感器"，这个频率通常不是随便定的。市面上很多压电圆片的径向振动谐振频率，就和直径、厚度直接相关。直径20mm左右、厚度1~2mm的PZT圆片，径向谐振频率往往落在100kHz~250kHz这个区间。模拟的时候我们不必真的去复现谐振腔的阻抗特性，但模型里压电片的直径和厚度应该参考这种实际规格来建，这样仿真结果才和实验有可比性。我自己的习惯是激励频率确定之后，再回头查一下对应尺寸压电片的谐振频率范围，确保激励频带没有离谱地偏离。

1.3 COMSOL为什么适合这个仿真

做Lamb波仿真的工具其实不少，纯有限元软件都能跑结构力学，但压电片场景最核心的需求是压电材料的多物理场耦合：电场、机械场、压电本构三个关系要在同一个单元里解出来。COMSOL里的压电设备接口（Piezoelectric Devices）把固体力学和静电耦合在一起，材料矩阵里直接给弹性常数、压电应力常数和介电常数，求解时自动形成双向耦合。激励片用的是逆压电耦合，接收片用的是正压电耦合，两个物理过程在同一个求解器里完成，不需要手工折算等效载荷。

我用过一些版本，早期版本里压电耦合的边界条件设置比较繁琐，物理场接口分离，需要手动加多物理场耦合节点。到了6.x版本，压电设备接口默认把正压电和逆压电效应都打开，电极设置也有专门的终端和浮电位边界，做一发一收模型省了不少事。还有一个很实际的好处：COMSOL后处理里可以直接把接收片的浮电位按时间画出来，得到的就是和实验示波器上形态接近的电压信号，方便和实测数据对照。这是个很顺手的验证闭环——仿真信号和实验信号对得上，模型参数的可信度就上来了。

2. 模型搭建前的关键决策：几何、材料与频率参数

2.1 几何尺寸的取舍

几何是这类仿真里第一个坑。很多教程喜欢拿一块大板子举例，800mm×800mm的铝板，乍一看很干净，但算起来足够让工作站冒烟。实际上，Lamb波仿真完全不需要把板建得无穷大，关键是保证在你要观察的时间窗口内，边界反射不会混进直达波信号里，同时自身对外传播的衰减路径要够。

我推荐的第一版几何是400mm长、300mm宽、3mm厚的铝板。激励圆片中心放在距短边100mm的位置，接收圆片中心放在距离150mm的位置，也就是x方向上距另一边也剩150mm。这样激励片到接收片之间有150mm的声路径，而波从激励片出发到达最近边界，反射回来再到达接收片，路径差足够让直达波先被完整接收。以S0模式大约5.3km/s的群速度算，直达波约在28μs到达，而从激励片出发碰到x=0边界再反射到接收片，路径大约是100+150+100+150=500mm，到达时间约94μs，两者错开了66μs，足够分离。

压电片的尺寸我建议取直径20mm、厚度1mm。这个规格有两点考量：一是径向谐振频率范围接近160kHz激励的频带，二是压电片直径相对板厚比较大，耦合面积大，接收信号幅值可观。如果压电片直径太小，比如6mm的小圆片，160kHz的波长相对于片径很大，片内应变梯度小，接收灵敏度会很弱，仿真出来的电压幅值小到让人怀疑模型跑错了。尺寸选20mm是经验范围内的稳妥选择。

宽度方向取300mm也有讲究。Lamb波在板内是二维扩散传播的，不是纯平面波沿一条线传播，波前会在板宽方向展开。宽度太窄的时候，侧向边界反射会以曲面波前形式混进接收信号；300mm宽度配合200mm以内的传播距离，主观察窗口内侧向反射基本进不来。如果后续想进一步减小计算量，可以沿传播中线用对称面把模型切半，这样板宽其实只建150mm就够，自由度直接少掉将近一半。

2.2 材料参数：铝和PZT怎么给

材料参数这块看起来简单，其实坑都埋在后面。铝板用6061-T6的参数就够，各向同性弹性：弹性模量E=68.9GPa，泊松比ν=0.33，密度ρ=2700kg/m³。需要注意的是Lamb波传播速度对剪切模量非常敏感，剪切模量和杨氏模量、泊松比直接挂钩，所以哪怕铝的牌号换一下，E和ν稍微变一点，S0和A0的到达时间就会差好几个微秒。仿真要对比实验的时候，最好以实际板材的超声波纵波和横波声速为准反推弹性常数，而不是照抄手册。

压电材料有一个重要的认知要建立：不要自己在材料节点里手敲那几十个弹性常数和压电常数，直接用COMSOL材料库里的PZT-5H或PZT-5A。材料库里的数据格式是完整的，弹性矩阵、耦合矩阵、介电矩阵按压电本构关系组织好了，选完就能用。但我强烈建议建完模型之后，从材料节点里把压电矩阵调出来看一眼，确认极化方向。COMSOL材料库默认极化方向通常是z轴，如果你的压电片是贴在铝板顶面上、厚度方向沿z轴，那刚好对得上；如果几何摆放时让厚度沿y轴或x轴，那就必须在材料坐标系里做旋转。

PZT-5H和PZT-5A的区别主要在压电常数和介电常数上。PZT-5H灵敏度高、介电常数大，适合做接收；PZT-5A温度稳定性好一些，发射效率也不差。仿真第一版我习惯选PZT-5H，因为接收电压幅值更好看，后处理时信号特征明显。到了要和实验做定量对比的阶段，再换成实际探头的型号参数。压电片密度通常取7500kg/m³左右，材料库会给，不用额外设置。

2.3 电极、接地与初始条件

电极在COMSOL里不需要建立真实的金属层几何。压电片的上表面和下表面，通过边界条件直接设置成电极即可。激励片的顶面电极设置为终端边界，类型选电压，然后在这个终端上施加你定义的时变电压；底面电极设置为接地，即电位恒为0V。注意终端和接地不能同时放在同一个面上，一上一下刚好构成厚度方向的电场回路。

接收片的电学边界是另外一个逻辑。接收片顶面设置浮电位，意味着这个电极上没有外部电路约束，电荷在整个面上自由分布，电压由压电效应自发产生；底面接地。这样接收片就相当于一个开路传感器，输出的是开路电压。为什么要用浮电位而不是直接测试某个点的电势？因为压电片电极面上的电势并不是均匀的，实验中探头电极是整片金属箔，测到的电压是电极面上的平均效果。浮电位边界正好模拟了这个"电荷自由移动、表面等电势"的条件，提取出来的电压才有实验对比价值。

初始条件方面，所有位移场和电势场初始值保持默认的0即可。整个系统从静止状态开始，激励电压从第一个时间步开始作用，不涉及预应力或初始应变，不需要额外设置。需要注意的一点是模型要避免刚体位移：铝板和压电片在空间里如果没有任何约束，瞬态求解会出现整体的刚体漂移，波形直接飞到自由空间里。所以铝板的某个角点或一条边要做低约束，比如在远离接收探头的一个角点固定一个小区域，约束三个方向位移，让板不能整体平移或转动。但约束点离传播路径足够远，避免影响波场。

3. 仿真核心细节：激励信号、网格与时间步的三角关系

3.1 激励信号：为什么用汉宁窗调制

激励信号不是简单给一个160kHz的正弦波持续振荡就行。你需要的是带宽可控的窄带脉冲，让接收信号里能分辨出不同模式的到达时刻。如果信号是一个连续的稳态正弦，波会不停地从激励片向外辐射，反射波和直达波搅在一起，A扫信号完全没法看。

标准做法是用汉宁窗调制的正弦脉冲。以5周期为例，时间长度是5×T=5×6.25μs=31.25μs，幅值包络从零平滑上升再平滑下降到零。公式可以写成：

V_in(t) = A × sin(2π f t) × 0.5[1 − cos(2π f t / 5)]

其中 f=160kHz，t的范围从0到5/f=31.25μs。A可以取10V或100V。仿真里电压幅值不影响波传播速度，只影响信号幅值，所以第一版用10V就可以，计算得到的接收电压按线性比例放大到实验值再对比。

为什么非要用汉宁窗？从频谱角度看，不加窗的矩形截断正弦，频谱旁瓣很高，能量泄漏到远离中心频率的频带。而Lamb波的频散特性意味着不同频率分量传播速度不同，旁瓣能量传播后会在时间轴上被拉开，形成拖尾噪声，污染S0和A0的波包形态。汉宁窗把频谱主瓣集中到中心频率附近，同时压低了旁瓣。调制周期越多，频带越窄，信号越纯净；但周期太多，脉冲时间变长，S0和A0两个波包在时间轴上靠得近的话会重叠。5个周期是160kHz下比较折中的选择，频带半宽大约在±16kHz左右，不足以激发高阶模式。

需要特别提醒的是，COMSOL里定义这种分段函数不要用if语句硬切。建议用解析函数加平滑阶跃，或者直接利用内置的窗函数类型。如果用手写的分段函数，在t=0和t=31.25μs两端，函数值虽然连续但导数不连续，瞬态求解器在这里会产生局部的高频振荡，相当于额外注入噪声。实际做法是把函数定义成0到T_pulse的有效区间，在边界处让函数值归零，同时让求解器的时间步长足够捕捉脉冲起始沿。

3.2 网格尺寸：按最短波长定

网格是Lamb波仿真里最容易让人抓狂的部分。网格太粗，波传播速度和衰减都算不准，波形畸变；网格太细，自由度爆炸，一夜都跑不完。这里有一个明确的下限逻辑：要分辨一个空间波长为λ的波，至少要保证每个波长内有6到10个单元节点。少于6个点，波会出现明显的数值频散，即不同频率分量传播速度被求解器扭曲，波形在传播过程中逐渐"散架"。

那么160kHz激励在3mm铝板上的最短波长是多少？A0模式在这个频厚积下相速度大约在2.2~2.6km/s之间，波长约为14~16mm；S0模式相速度约5.0~5.3km/s，波长约31~33mm。最严格的情况应该按A0的波长来定网格。按每波长8个节点，需要单元尺寸不超过约2mm。但这个结论只考虑了面内的波长分辨，厚度方向还有另一重约束：铝板厚度3mm，至少保证4层单元，否则厚度方向上的应力应变梯度无法正确表示，尤其是A0这种反对称模式，上下表面的应变相位是相反的，厚度方向单元不够，模式形态会失真。所以厚度方向建议以0.75mm为层厚，整个板厚划分4层。

综合起来，铝板区域我用面内1.5mm、厚度0.75mm的扫掠六面体网格，压电片区域因为尺寸小、且要解析厚度模式，网格加密到0.5mm。整个模型规模大概是：板面内267×200个单元，乘以4层，大约21万个单元，再加上两块压电片区域的加密网格，总量在25万单元左右。自由度数量大约是单元的4到5倍，也就是百万量级。这个规模在COMSOL里用瞬态求解器跑是完全可控的，64GB内存的工作站一晚上能出结果。

新手常犯的错误是全局统一用极细网格，比如0.3mm，结果自由度冲到几千万，还没看到波形就把内存耗尽了。正确的顺序是先按波长估算出单元尺寸，然后用"映射+扫掠"的方式控制厚度分层，而不是让自由四面体网格在厚度方向乱分。

3.3 时间步长与求解时长

时间步长和网格尺寸是绑定的。显式求解要求时间步满足CFL条件，COMSOL的瞬态默认使用隐式BDF求解器，理论上可以放宽步长，但波传播问题如果时间步长太大，时间方向上的数值色散会让波形变钝、相位偏移。经验法则：时间步长等于最小单元尺寸除以最大波速再乘一个安全系数。

这里最大波速取铝中纵波速度约6320m/s，最小单元尺寸取厚度方向的0.75mm。0.75mm除以6320m/s约等于0.119μs，再乘0.8的安全系数，得到时间步长约0.095μs。实际设置可以直接取0.1μs，既满足分辨率要求，又不至于让总步数太多。

求解时长怎么定？主要看你要观察多少信号。接收片距离激励片150mm，S0约28μs到达，A0约60多μs到达，如果要看到完整的直达S0和A0波包，再留出一些余量观察可能的边界反射，取120μs比较合适。120μs除以0.1μs等于1200步，每一步都是一个完整的瞬态求解过程，这个步数在25万单元的规模下，普通工作站几个小时到一晚上能跑完。如果板再大、网格再细，建议先只算200μs以内的第一波信号，不要贪多把反射全收进来。

求解器设置里有一个细节很容易被忽略：COMSOL的瞬态BDF求解器默认允许内部时间步进器自动伸长步长，以保证收敛和效率。但对波传播问题，这个自适应拉伸会破坏信号的简单谐波结构，让波形出现莫名其妙的畸变。建议在时间步进设置里把"允许步长变大"或者等效的选项关掉，强迫求解器按指定的最大步长时间网格输出，代价是计算时间略增，但波形质量明显提升。

4. 实操走一遍：COMSOL里的模型构建与求解流程

4.1 物理场和几何操作

建模从新建一个三维组件开始。第一步用块体命令建立铝板，尺寸X=400mm、Y=300mm、Z=3mm；再建立两个圆柱，直径20mm、高度1mm，分别放置在激励位置和接收位置。圆柱底面和铝板顶面完全重合，它们之间不需要额外设置接触，因为压电设备接口里的固体力学方程默认所有连续实体通过共享边界传递位移和应力。只要几何布尔操作里保证压电片底面和铝板顶面是同一个面，网格划分后会形成连续网格，机械波可以无缝穿过界面。

物理场添加"压电设备"接口，在结构力学模块下。COMSOL 6.x里这个接口默认包含固体力学和静电两个子接口，并且自动生成了压电效应的多物理场耦合。在物理场设置里，域选择要把铝板和两块压电片全部包含进去，但注意铝板不需要压电本构，所以后续要在材料分配上做区分：铝板用线性弹性材料，压电片用压电材料。

几何和物理场就位后，先做一次稳态小测试验证耦合是否正常。可以在激励片的终端上加一个恒定的直流电压，然后跑一个一步静态求解，检查接收片浮电位是否出现一个非零的微小电压。如果这一步接收电压是零，说明压电耦合方向或者材料极化方向有问题，趁早排查，不要等到瞬态跑完再后悔。

4.2 激励加载、接收探针与求解器

激励电压函数先在组件级别的定义节点里建立解析函数。函数类型选"分段"，在0到31.25μs区间定义V=A×sin(2π×160kHz×t)×0.5×(1−cos(2π×160kHz×t/5))，区间外函数值均为0。如果函数起始端有小幅跳变，可以在解析函数设置窗口的平滑选项里打开连续一阶导平滑，给脉冲前沿加一个几纳秒的过渡。

物理场里对激励片上表面电极用"终端"边界条件，终端类型选"电压"，终端电压关联到刚才定义的解析函数。底面电极用"接地"边界。接收片上表面用"浮电位"边界，底面同样接地。注意浮电位边界的求解变量是一个未知常数，在结果数据集里可以提取这个常数的值，它就是接收电极上的开路电压。

研究步骤选择"瞬态"，时间范围用"range(0,0.1e-6,120e-6)"，代表从0到120μs、每0.1μs输出一个结果点。求解器配置里，线性求解器选直接求解器PARDISO。三维压电多物理场瞬态问题，刚度矩阵是非对称的，迭代求解器在波动问题上收敛性不好，PARDISO最稳妥。如果模型规模超过200万自由度且内存吃紧，可以改用MUMPS，它分了块内存管理，在某些机器上反而比PARDISO省内存。求解容差用默认的0.01就可以，若发现波形尾部有毛刺，再收紧到0.001，但计算时间会有明显增加。

并行方面，COMSOL默认会用满所有物理核心，但对波动问题不建议盲目追求核心数量。PARDISO在超过16核后加速比提升有限，尤其是小规模模型，线程切换开销反而拖慢速度。我自己习惯控制在物理核心数的75%左右，让其他程序也有余量。

4.3 后处理：A扫信号和波场

结果分析的第一步是做A扫信号图。在派生值里添加"全局计算"，表达式选接收片浮电位边界的电压变量，或者直接对接收片顶面做面积分平均电位。1D绘图组里把这条电压随时间变化的曲线画出来，你应当看到两个明显的波包：第一个约在28μs附近，幅值相对小，这是S0模式；第二个约在60μs附近，幅值明显更大，波包更宽，这是A0模式。两个波包之间有一段相对平静的区间，说明模型里信号是干净的。

如果想进一步观察波场，在结果里创建二维截面试图，截取铝板中面或顶面，画位移绝对值或者面外位移分量。从云图动画里能直观看到波前从激励片出发向外扩散，S0波前面外位移分量小、面内位移分量大，A0波前则相反，面外位移非常明显。这种模式形态的差异可以作为判断模态类型的经验依据。

提一个后处理的小技巧：如果接收信号里S0波包和A0波包靠得太近，分辨吃力，可以换一个提取量。比如把接收片顶面的平均面外位移也和浮电位曲线一起画出来，在面外位移里A0占主导，在面内位移里S0占主导，对照两张图可以更清楚地确定各波包到达时间。

5. 波形验证：怎么确认仿出来的Lamb波是对的

5.1 用群速度对照频散曲线

仿真跑完不是终点，关键一步是验证结果的物理正确性。最简单有效的方法是测量波包到达时间，换算成群速度，再和理论频散曲线对比。

具体操作是：在A扫曲线上找到S0波包的峰值时刻t1，和A0波包的峰值时刻t2。已知传播距离L=150mm，则群速度vg=L/t。我常用频散曲线计算程序事先算好3mm铝板在160kHz附近的S0群速度约5.3km/s、A0群速度约2.5km/s左右，然后用仿真测得的到达时间去反推实测群速度。如果误差在5%以内，基本可以认为模型的弹性参数、几何尺寸、材料密度都对得上。误差偏大时，优先检查铝板的杨氏模量和泊松比是否有偏差，其次是网格和时间步长的分辨率是否足够。

这里特别提醒一个时间参考点的问题：波包峰值到达时间记的是"从激励脉冲起始瞬间到接收波包峰值"的绝对时间差，而激励脉冲本身有一个包络上升过程，所以实测到达时间比理论群速度到达时间会偏大半个脉冲宽度左右。更严谨的做法是用两个不同传播距离的模型做差分，或者直接对比整段波形而不是只看峰值。如果只做单距离仿真，就把误差容忍度放宽到10%以内，重点看两个模式到达时间之间的相对间隔。

5.2 波场云图与动画里的模态特征

单凭A扫信号还不够全面。打开波场云图动画，你会看到激励后大约10μs时，波前刚刚离开压电片向四周扩散。此时注意观察波前形状和位移分量云图的分布：S0传播速度快，波前面内位移显著，波前较平直；A0速度慢，面外位移显著，波前面伴随明显的弯曲形态，波长相对更长。到了60μs以后，A0波前在云图里形成清晰的大波浪形，沿板面向四周扩散，这非常符合低频A0模式的特征。

再看接收片附近的局部应力分布。激励片正下方铝板表面有强烈的应力集中，应力方向沿厚度和径向都有分量；接收片下方的应变场则是入射波导致的均匀面应变。如果这里的应力分布出现锯齿状条纹或者棋盘格状图案，说明网格分辨率不够，波被数值方式散射了，要加密网格重新算。

动画里另一个值得观察的现象是边界反射。如果你把求解时间拉长到200μs以上，会看到波前碰到板边缘后产生反射波向板内传播。在云图里，边界反射的波前弧度和直达波明显不同，可以直观判断边界反射什么时候到达接收区域，从而确定你关心的信号是否被污染。

5.3 尺寸缩比与对称性建模技巧

仿真验证通过后，你可能还想变换参数做系列分析，比如不同板厚、不同探头间距、不同频率。在这些场景下，不一定每次都建全尺寸三维模型。板结构Lamb波沿传播路径的物理特性，很大程度上可以用二维平面应变模型来复现。二维模型把压电片简化成一条线状的机电载荷区域，虽然无法模拟圆片径向波前扩散和接收面的电荷平均效应，但可以快速扫参，把大致趋势摸清楚。

三维模型里也可以利用对称性。沿激励片和接收片中心的连线建立对称平面，只建一半模型，对称面上施加对称边界条件。这样做需要几何布局也对称，所以板宽方向最好取偶数尺寸。半模型自由度少一半，内存和计算时间差不多能省掉一半，这在做参数扫描时会让你明显感觉轻松。

还有一种常见的方案是在模型外围设置吸收层。COMSOL里可以做完美匹配层（PML），但PML在固体波里要设置得足够厚才能有效吸收低频弯曲波，而且参数设置稍有不慎会产生虚假反射。如果只是做单次验证，用大板子加短时间窗口更省心；只有当你需要长时间观察多次反射或者模拟无限板条件时，才值得花时间调PML参数。

6. 常见问题与避坑清单

6.1 计算慢、内存爆的排查

最典型的问题是模型跑不动，或者跑起来内存占用一路飙到顶。遇到这种情况不要急着换机器，先按这个顺序排查：建模区域是否缩到最小、对称条件是否用上、网格是否严格按波长控制、求解时间是否设置过长、时间步长是否过小。我见过不少情况是把求解时间拉到500μs，步长设到0.01μs，这样步数高达5万步，神仙机器也扛不住。记住目标只是"看到直达S0和A0波包"，找到满足这个目标的最短时间窗口。

6.2 信号形态异常的一般排查顺序

接收信号如果出现三种典型异常，排查优先级不同。第一种是完全没有信号，先查压电片的材料域选择是否正确，再查压电本构是否被分配给了铝板，最后检查压电材料极化方向。第二种是信号极其微弱，比正常情况小几个数量级，重点查浮电位边界是否设置成功，以及压电片底面接地和激励片底面的电位参考是否正确。第三种是波形有剧烈振铃或者拖尾振荡，重点查激励函数的平滑性、时间步长是否太大、网格在压电片区域是否太粗。

6.3 极化方向这个隐蔽的坑

压电片极化方向错了，仿真不会报错，但信号特征会和理论完全对不上。COMSOL材料库里的压电材料数据，默认情况下极化方向沿坐标系的z轴。你组建几何的时候，如果压电片厚度方向沿z轴，那没问题；但如果旋转了几何或者自己定义了不同朝向，材料本构矩阵不会自动跟着旋转。解决方法是给压电片单独分配一个坐标系或者旋转后的材料方向，让极化方向始终沿片的厚度方向。我建议建完模型后做一个静态验证：激励片加直流电压，看接收片浮电位是不是同号（如果激励片和接收片极化方向一致，直流下接收电压应该是正或负的对应状态）。一旦静态验证通过了，极化方向的正确性基本可以放心。

6.4 边界反射污染怎么缓解

即使板子尺寸选得合理，接收窗口尾部还是会混入部分反射波。特别是A0模式低频成分传播慢，板边反射回来之后和直达波在时间上可能重叠。缓解方法有三层：第一层，保证你分析信号的时间窗口在反射到达之前截止；第二层，把模型加大或加吸收层，但这会增加计算量；第三层，在接收信号后处理时做时间窗截断，只提取直达波窗口内的数据。对于一发一收路径分析来说，只要首到波包的时间足够清晰，反射污染并不影响核心参数的提取。

我做这一类仿真下来，最深的体会是：压电Lamb波仿真不像静力学结构分析，跑完看个应力云图就完事，它特别强调"信号处理思维"。建模的时候就要想着这段波形后面要用来做什么——是要看到达时间、要看模式幅值、还是要做缺陷散射分析。带着这个目标去定模型尺寸、网格密度、时间窗口和探针位置，你建出来的模型才不是花架子，后面和实验对照或者做参数研究时，能少走非常多的弯路。