当前位置：首页 > news >正文

基于传递矩阵法计算多层结构声表面波声速 - 教程

news 2026/5/12 19:06:46

一、理论基础与模型建立

1.1 传递矩阵法基本原理

传递矩阵法通过将多层介质分解为独立单元，建立各层界面处的声压和法向速度连续性条件，形成系统级传递关系：

[ P_{m+1
} ] = [ T_m ] [ P_m ]
[ V_{m+1
} ] [ ] [ V_m ]

其中，T_m为第m层介质的传递矩阵，P为声压向量，V为法向速度向量。

1.2 声表面波特性的数学描述

对于压电多层介质，声表面波(SAW)的传播满足耦合波动方程：

∇×(C:∇u) + ρ∂²u/∂t² = 0
∇×E = -∂B/∂t

通过引入机电耦合系数k² = (C₁₁/C₃₃)·(ε₁₁/ε₃₃)，可推导出SAW相速度频散方程。

二、多层结构建模步骤

2.1 单层介质传递矩阵推导

对于厚度为h、密度为ρ、弹性刚度矩阵为C的介质层，其传递矩阵T可表示为：

function T = layer_transfer_matrix(h, rho, C)
% 计算波数矩阵
k = sqrt(C(1,1)/rho);
% 纵波波数
% 构建传递矩阵（以ST切石英为例）
T = [exp(-1j*k*h), 0;
(k/h)*exp(-1j*k*h), exp(1j*k*h)];
end

2.2 多层结构总传递矩阵

通过级联各层传递矩阵获得系统总传递矩阵：

function T_total = multi_layer_matrix(layers)
n = length(layers);
T_total = eye(2);
for i = 1:n
T_total = T_total * layer_transfer_matrix(layers(i).h, layers(i).rho, layers(i).C);
end
end

三、声速计算与频散分析

3.1 相速度计算公式

通过总传递矩阵的本征值分析，得到SAW相速度：

v_p = ω / Re(k_eff)

其中，k_eff为总传递矩阵的特征波数。

3.2 机电耦合系数计算

机电耦合系数K²反映能量转换效率：

K2 = (abs(det(T_total)) - 1) / (2*abs(det(T_total)));

四、MATLAB仿真实现

4.1 参数定义

layers = {
struct('h', 0.5e-6, 'rho', 2650, 'C', [390 120 120;
120 390 120;
120 120 390]), % ST-cut石英
struct('h', 1e-6, 'rho', 7500, 'C', [120 120 120;
120 120 120;
120 120 120]) % ZnO
};

4.2 传递矩阵计算

T_total = multi_layer_matrix(layers);

4.3 相速度频散曲线

f = linspace(100e6, 5e9, 1000);
% 频率范围
omega = 2*pi*f;
k_eff = zeros(size(omega));
for i = 1:length(omega)
[~, eigenvectors] = eig(T_total);
k_eff(i) = angle(eigenvectors(1,1)) / (h_total);
end
v_p = omega ./ real(k_eff);
plot(f, v_p*1e6);
% 转换为μm/s
xlabel('Frequency (Hz)');
ylabel('Phase Velocity (μm/s)');