搞过超声仿真或者接触过压电换能器的人应该都听过“横波激励”这个词。刚入行那会儿，我对着COMSOL里那堆物理场接口和边界条件看了好几天，愣是没搞明白怎么让模型里产生一列干净的横波。后来踩了不少坑，翻了无数篇案例文档，才算是把横波激励这条路走通了。这篇就聊聊我在COMSOL里做横波激励的那些实操经验，从物理概念到建模设置，再到后处理看结果，一次性说清楚。
1. 横波激励是个啥：先弄清物理本质
1.1 纵波和横波的区别
先说物理基础。弹性介质里的波分为两大类：纵波（P波、压缩波）和横波（S波、剪切波）。纵波传播时质点振动方向和波的传播方向平行，像弹簧一样一伸一缩。横波则正好相反，质点振动方向垂直于波的传播方向，类比一下就是抖绳子的那个效果，手上下抖动，波沿着绳子水平传。
在COMSOL仿真里，这个区别直接决定了你要用哪个物理场接口、怎么设置边界条件。纵波的激励相对简单，直接在某条边上给一个法向位移或力载荷就行。横波麻烦得多，你需要在边界上施加一个切向的力或切向位移，让介质产生剪切变形。而且横波的传播速度比纵波慢，波长更短，网格尺寸的要求也更高，这个后面细说。
1.2 横波激励的三种典型路径
我在实际项目中遇到过的横波激励方式，大体可以归为三类：
路径一：压电换能器直接激励。这是最自然也最常见的一种。压电陶瓷在逆压电效应下产生厚度剪切振动模式（比如常用的d15模式），通过粘接层把剪切位移传给被测介质。这种方式的优势是频率可控、波形干净，无损检测领域用得非常多。
路径二：边界载荷直接施加切向力。适用于纯理论研究或验证算法。不考虑换能器的具体结构，直接在模型边界上设置横向位移或切向应力，观察介质内部的横波传播行为。
路径三：激光超声或热弹效应。激光照射材料表面产生热膨胀，热弹机制下会同时激发纵波和横波。这种方式在COMSOL里需要耦合固体传热和固体力学两个物理场，设置相对复杂。
从我个人的经验看，如果你刚开始接触横波激励，建议从路径二入手，物理场少，变量清晰，方便你把横波传播的基本特征摸透。等理解了波场特征，再上手压电换能器的完整建模。
1.3 为什么COMSOL特别适合干这个
COMSOL Multiphysics的优势在于多物理场耦合。横波激励这个问题，表面看只是固体力学的事，但只要涉及真实的换能器，就必然是压电效应（AC/DC模块的静电接口）和弹性波传播（固体力学接口）的联合仿真。再加上你要模拟有限尺寸试件里的波场，还存在边界反射问题，需要设置完美匹配层或者低反射边界。COMSOL的结构化建模流程——几何、材料、物理场、网格、求解器、后处理——对这类多物理场耦合问题的支持确实非常成熟。
2. COMSOL里横波激励的建模思路
2.1 物理场选型
横波激励的核心物理场是固体力学（Solid Mechanics）。但根据你的激励方式，需要搭配其他接口。
如果是压电换能器激励，需要用到压电器件（Piezoelectric Devices）接口。COMSOL 6.x版本的压电接口内置了d15模式的剪切耦合选项，比老版本方便不少。这个接口是把静电（Electrostatics）和固体力学耦合在一起，压电本构方程里的应力—电场耦合矩阵就是关键。
如果是纯理论验证，只选固体力学接口就够了。在边界上施加切向位移载荷或者切向应力，甚至在域内条件里设置初始剪切应力分布，都能完成横波激励。
还有一个容易被忽略的选择：弹性波接口（Elastic Waves interface）。这是COMSOL 5.4之后新加入的高频弹性波专用接口，主要面向时谐分析，计算效率高，PML（完美匹配层）设置也简化了很多。如果你的目标是计算频域响应、研究波的反射折射，弹性波接口会比固体力学接口好用很多。但如果是瞬态时域分析，想看波形传播过程，固体力学接口配合时域求解器更直观——因为弹性波接口的默认框架是频域，要做瞬态需要额外改求解器设置，反而多一道手续。
我的话，分情况：验证算法的基本概念用固体力学时域，分析换能器频率特性用压电接口频域。两者各干各的活。
2.2 压电换能器的核心参数设置
先建换能器模型，PZT压电陶瓷是最常用的材料。COMSOL材料库自带PZT-5A、PZT-5H、PZT-4等常见压电材料，材料属性比较全，直接用也行。不过我建议你关注一下极化方向这个参数。
在压电物理场中，极化方向（极化轴）决定了耦合矩阵的取向。对于厚度剪切模式的换能器，你要确保压电材料的极化方向与后续设置的激励电场呈特定关系。比如常见的d15剪切模式，需要把极化方向设定为x方向，然后在y方向施加电场，这样才能激发xz平面的剪切应变，进而向外辐射横波。这个方向和坐标轴的对应关系搞错了，模型算出来就是一团乱麻，介质里啥波都激发不出来。
具体的做法是，在压电材料定义里有一个“极化方向”选项，把坐标轴选对。同时要设置好压电材料坐标系的旋转。COMSOL材料坐标系默认与模型x、y、z对齐，如果你的换能器切型是旋转过的，需要增加坐标旋转操作。我吃过这个亏——一开始模型里看起来是剪切变形，但波的传播方向总是反的，后来发现是材料坐标系没旋转到位。
2.3 边界条件布置的几个坑
边界条件有两个方向容易踩坑。一个是激励边界本身，另一个是人工截断边界。
激励边界。如果你用边界载荷法，直接在边界上选“边界载荷”，类型选“力”或“位移”，方向需要设定为切向。COMSOL的边界载荷默认根据边界法向来确定方向，你需要在分量设置里把切向力写成类似{0; Fy; 0}的分量形式（具体根据边界朝向调整）。这里有个隐蔽的问题：如果边界法向不完全对齐某个坐标轴，切向分量的分解就会出错。所以我的习惯是，在画几何时就把激励面的法向对齐到坐标轴方向，省得后续设置时还要做矢量投影计算。
人工截断边界。实际物理模型中，波应该向无限远处传播，但仿真模型一定是有限尺寸的。边界上波会反射，干扰横波的波形。解决手段有两个：一个是低反射边界（Low-Reflecting Boundary），一个是完美匹配层（PML）。
低反射边界在COMSOL里就是边界条件里直接勾选的一个选项，适合大多数一般用途，能吸收大部分入射波，但做不到完全无反射，尤其对于斜入射的横波效果打折。PML则是在模型外圈额外画的一层“吸收域”，里面设置了与物理属性匹配的阻尼参数，能显著吸收全向入射波。代价是要多画几何域、多剖分网格。如果你的模型需要看到非常干净的主波、不影响波形分析，PML才是正解。
在COMSOL 6.x中，PML设置已经模板化了，在域设置里把材料属性“材料域类型”改成“匹配层”，然后选“弹性波”类型就可以了。但注意：PML的厚度至少需要覆盖一个波长，层数建议6到8层网格，这是收敛的关键。很多人的模型反射严重，一查就是PML层太薄了。
2.4 网格剖分：横波仿真最容易翻车的地方
横波仿真对网格的要求远远高于纵波。原因很简单：相同频率下，横波波速约为纵波的0.5到0.6倍，所以横波波长只有纵波的60%左右。要解析一个波，每个波长至少需要5到6个二阶单元，这是COMSOL默认规则。
举个例子，假设材料是钢，纵波波速5900m/s，横波波速3200m/s。设定激励频率1 MHz，横波波长就是3.2mm，纵波波长5.9mm。要算横波，网格尺寸得控制在0.5mm左右，才能保证每波长有6个以上单元。如果是把网格按纵波波长尺度画到1mm，那横波波形基本就已经被离散数值耗散吃干净了。
还有一种情况是使用一阶单元，那网格要求更苛刻，每波长至少要10个单元。我的经验是：直接在网格节点里设置“最大单元尺寸”为横波波长的六分之一，同时保证时间步长满足CFL条件，即时间步长不大于网格尺寸除以介质最大波速（通常取纵波波速）。两者配合，波形才不会畸变。
3. 实操：一个完整的横波激励仿真流程
为了让你更有体感，我把一个真实的仿真案例从头到尾过一遍：一个长方体钢制试件（长40mm，宽20mm，厚10mm），在上表面用一个压电陶瓷片（厚度剪切模式）产生横波，激励中心频率为1MHz。目标是观察试件内部横波从顶部向底部传播的过程。
3.1 第一步：几何建模与材料定义
几何分两部分：压电片和试件。压电片尺寸设为8mm×8mm×0.5mm，贴在试件上表面中央。用COMSOL的“阵列”功能把两个域组合一下，或者直接通过工作平面拉伸生成，都可以。
材料选取上，压电片选PZT-5H，试件选结构钢（Structural steel）。这两个材料都在材料库中内置，选完即用。注意把PZT-5H的密度、弹性矩阵、压电矩阵都检查一遍。材料库自带的数据可能在某些版本里压电矩阵符号约定不一样。COMSOL默认的压电本构关系用的是stress-charge形式，这个形式跟很多教科书的表达略有出入，计算前确认一下符号约定，否则计算收敛了结果也是错的。
3.2 第二步：物理场与激励信号设置
我推荐在“压电器件”接口里完成激励设置。在压电域上定义“静电”的“电势”边界条件：下表面接地（0V），上表面施加交流电压信号。
信号形式有两种主流的做法：一种是频域分析，也就是直接设置固定频率1 MHz正弦激励，求解稳态响应；另一种是时域分析，用脉冲信号激励，观察瞬态波传播。如果你用瞬态分析，信号可以用高斯调制正弦脉冲宽度5个周期，幅值50V。具体可以写为：
V(t) = 50*sin(2*pi*1[MHz]*t)*exp(-((t-5[t_us])/(1.5[t_us]))^2)在COMSOL里，直接把这个表达式写在电势边界条件的“电势”栏里50[V]*sin(2*pi*1[MHz]*t)*exp(-((t-5[us])/(1.5[us]))^2)就行。
3.3 第三步：求解器配置
瞬态求解我一般用BDF（向后差分法，COMSOL默认）或者广义alpha方法。广义alpha在高频瞬态波传播问题上表现更好，相位误差小，阻尼少。BDF二阶梯度的数值耗散在比较长的传播距离下会导致波形振幅和速度明显衰减，做瞬态大模型时推荐切换广义alpha。
时间步长设定用自适应模式就行，但为了保险，在求解器配置里把“最大时间步长”限制在0.02[us]（大约为CFL条件的二分之一）。如果你还想更图省事，也可以开启“基于网格尺寸限制时间步长”选项，COMSOL会自己算。
求解范围设定为0到20μs，足够横波从顶部传播到底部再反射回来。其实横波在钢中传播40mm只需要12.5μs左右，给点余量。
3.4 第四步：后处理与结果读取
最重要的一步在后处理。计算完成后，你先在默认的“表面”图里查看位移场分布。能不能看到横波？横波的特征是质点位移垂直于传播方向。所以你可以在结果里绘制水平位移分量（比如在x方向和z方向之间），检查波前形态是否为“剪切型”，并在波传播路径上放置一组探针来逐点显示位移时间历程。
我最常用的是“一维绘图组”里的“点图”，取试件中心线上几个点，用“全局”或“数据集”方式在点上同时显示位移分量。横波到达某点时应能清晰看到一个波包，它的幅度突变、起振时间对应横波到达时间。可以根据这个时间差来反推波速，和理论值对照。如果波包不明显，波形混乱，大概率是激励方向设置错了或者网格太粗。
还有一个经验是：横波到达前，受纵波前驱波影响，波形往往会先有个小波动。不要把这个小波动当成横波。横波信号的幅度通常比纵波大，出现的时刻晚一些，而且振动方向与传播方向垂直。分析数据时要把纵波的前驱效应区分开。
4. 常见问题与排查技巧实录
4.1 问题1：零能量输入，模型没反应
症状：算了半天，整个结构纹丝不动，位移场全是零。
排查思路：先检查激励有没有真的加载到模型上。在“研究”里计算之前，先跑一次“稳态”或“特征频率”看看能不能算；如果能算但没位移，多半是激励方向写错了。边界载荷里如果力的分量全写成了法向，横波激发不出来——法向力产生的是纵波，只有切向分量才能产生剪切波。另外，电势边界条件如果压电域上下面的电压都设置为0V，或者电极面选错，逆压电效应不成立，自然没有激励。
这个问题的另一个隐蔽原因是压电材料坐标系不对。我遇到过好多次，压电矩阵里d15系数被映射到坐标轴的错误组合，导致电场方向与实际剪切变形方向不匹配，结果就是模态看起来是“死的”。解决方法是输出压电材料的耦合矩阵，手动验证d15对应的分量方向和你的几何坐标是否对应。
4.2 问题2：网格泄能导致波速异常
症状：横波传播速度比理论值偏慢，波形展宽，幅度下降明显。
排查思路：网格太粗。前面提的每波长6个单元的规则必须严格执行。时间步长过大也会造成类似问题，需要把最大时间步长缩短到横波通过一个网格所需时间的二分之一以内。更好的方式是保证网格在频率上限对应波长内都足够细。多频激励下，中心频率1 MHz的脉冲，频谱扩展到了1.5 MHz甚至2 MHz，网格如果你只罩住了中心频率，高频成分没法被解析，波形自然失真。
4.3 问题3：PML怎么配才能做到无反射
症状：波传播到边界后反射回来，波形叠加非常严重，看不出主波包了。
排查思路：PML的域必须足够厚，至少一个横波波长。如果只为一个最小波长厚的PML，对掠入射波吸收差，反射就会来捣乱。第二要检查PML域的材料属性设置，确保已经把材料域类型切换为“匹配层”，而不是普通的弹性体。COMSOL里如果你忘了设，那PML只是一个普通材料层，波照常反射。
第三就是网格。PML内部要求网格边长从内到外逐渐变大，层数至少4-6层。直接使用“扫掠网格”配合“分布”功能，把PML层均匀分成6份，是最稳妥的实操路径。
4.4 问题4：模态混叠与结果解读
症状：时域结果里，前驱信号和后续波包混杂，主波包形状不再对称，看起来像有多个波。
排查思路：这其实是正常现象。纵波前驱波速度快、先到达，虽然幅度通常低，但在高灵敏度检测时会被看到。横波到达后，后续还有表面波、转换波和边界反射波。不要指望整个时间历程只有一列干净横波。如果只想看纯横波波型，可以考虑在激励频率上做窄带滤波，或者在边界上设置低反射边界条件来减少反射。
后处理上，我常用的一个技巧是：把速度场的旋度（curl）或者应变场的偏量部分绘制出来，因为纵波是无旋场，横波是无散场。旋度图里纵波成分被自动压制，横波轮廓就非常清晰。COMSOL里可以用衍生值功能，在体上算速度场的旋度分量并显示。
5. 从横波激励延伸到更复杂的场景
5.1 横波换能器的声场仿真
一旦完成了横波激励这块的基础建模，你可以往上叠很多实际工程需求。最常见的是换能器声场仿真。压电片辐射的横波通过楔块或直接接触进入被测介质，在介质中形成横波声场。
用COMSOL的“声-固耦合”能力，把压电模块产生的位移场作为声场激励源，计算介质内部的声压、位移势等。这里横波激励做得好不好直接决定了声场仿真的可信度。如果你在激励这层就错了，声场计算得再漂亮也没有意义。
5.2 横波激励在无损检测中的应用
无损检测（NDT）里横波应用的路子很宽：利用横波在界面上的反射、折射特性判断缺陷位置。COMSOL可以建模横波倾斜入射到裂纹缺陷后的散射场。设置一个带半椭圆形裂纹试件的模型，让横波沿着特定角度入射，之后在接收点测量回波信号。这类仿真需要小心设置入射波：通过边界位移施加倾斜方向位移分量，而不是简单勾选一个边界力。我做过一个类似的模型，横波斜入射到裂纹时，波场出现了明显的模式转换，纵波和横波在缺陷处互相耦合，单靠看位移云图里分析会很痛苦。这时候要在后处理里做时频分析或者指向性分析，才能定量判断缺陷尺寸。
5.3 与其他模块的耦合
横波激励的应用场景可不止无损检测。压电超声还能应用在医学超声、能量收集、微流控声学操控等领域。COMSOL的多物理场耦合灵活性在这里就体现出来了：压电—传热—固体力学—流体动压，甚至电化学都可以组合成一套复杂仿真链路。
以声流控为例：横波激励的压电换能器带动微流道内液体振动，产生声辐射力和声流现象。这时你需要压电器件接口、层流接口、动网格或者基于声辐射力的粒子追踪模块一起工作。横波激励在固体中的传播形态直接影响耦合到流体中的声压场分布，很多细节在仿真时都要注意匹配边界——声固耦合界面上位移连续、应力连续这两个条件必须同时满足。这也意味着在界面两侧的网格不能差太多，否则耦合矩阵病态，收敛困难。
我做这类耦合项目时，习惯分步走：第一步先单独验证横波激励，看换能器的频响是否跟实验对标；第二步再做声固耦合，看声压值是否合理；第三步才加入流体或粒子。每一层都锁定变量，排查起来轻松很多。
6. 最后再分享几个实战小技巧
几个容易忽视但特别能提升效率的点，收个尾。
一个是关于材料库参数。COMSOL材料库里PZT系列参数大多来自供应商数据表，不同版本可能会略有更新。每次做新模型前，建议跟供应商提供的最新材料手册核对一遍弹性矩阵和压电矩阵。压电材料的老化、温度漂移会影响仿真与实验的偏差，手里参数越准，后续分析越靠谱。
另一个是关于频域与瞬态的选择。用频域分析算横波传播时，虽然计算快、结果整洁，但看不出波包的瞬态特征和边界反射时间顺序。很多新手为了图省事直接跑稳态频域，结果把横波激励的物理过程完全给“稳态”掉了。项目前期建议先用瞬态把传播过程看清楚，之后再做频域精确计算，这样既有物理图景又有数值精度。
最后一个关于复制性。项目做完后，把关键波长、网格尺寸、时间步长、中心频率、材料坐标系这几个核心参数单独记录到一个配置文档里。下次要做不同尺寸、不同频率的仿真，直接按比例缩放修这几个核心变量就行，模型重建效率能翻一番。我自己就吃过好几次“新项目重新来一遍”的亏，数据记录做好了，整个团队都能受益。
