不同压电分流电路对声子晶体梁带隙影响的研究报告
作者及发表信息
本研究由国防科技大学三院机电工程研究所的陈圣兵、韩小云、郁殿龙、温激鸿共同完成。研究成果以题为《不同压电分流电路对声子晶体梁带隙的影响》的论文形式，发表于《物理学报》（*Acta Physica Sinica*）第59卷第1期，出版于2010年1月。该研究得到了国家自然科学基金（批准号：50875255）的资助。
研究的学术背景
本项研究属于声子晶体（Phononic Crystal）这一前沿交叉学科领域。声子晶体是一种具有弹性波带隙特性的周期性复合材料，其核心特性在于，特定频率（即带隙频率）范围内的弹性波或振动无法在其中传播。这一特性为减振降噪、声波滤波器、新型换能器等工程应用提供了全新的理论基础和设计思路。声子晶体的带隙形成机理主要分为两类：布拉格（Bragg）散射型和局域共振（Local Resonance）型。其中，局域共振型声子晶体由刘正猷等人于2000年提出，其优势在于能够在波长远大于晶格常数的低频范围产生带隙，突破了布拉格散射机理的限制。
随着智能材料的发展，将压电材料等智能结构引入声子晶体，通过电路等主动控制手段实现带隙的实时、可调控制，已成为该领域的一个重要研究方向。基于此背景，本研究旨在探索一种新颖的声子晶体结构设计：在匀质的基体梁上周期性粘贴压电陶瓷片（PZT），并将每个压电片连接外部分流电路。其核心科学问题在于，不同的分流电路（如纯电阻、LC振荡电路、LCR电路）如何通过电-机耦合效应影响压电片的等效力学参数，进而调控整个“声子晶体梁”的弹性波带隙特性（包括带隙的产生机制、位置、宽度和衰减强度）。研究的目标是为设计具有可调、主动控制能力的声子晶体减振结构提供理论指导和方法论依据。
详细的研究流程
本研究主要采用理论建模与数值计算相结合的方法，流程清晰，逻辑严谨，具体步骤如下：
1. 物理模型建立与理论推导： * 研究对象： 研究构建了一个一维声子晶体梁的简化物理模型。模型以环氧树脂材料制成的匀质基体梁为主体，在其上周期性地对贴压电陶瓷片（PZT-5H）。每个压电片的上下电极分别连接一个外部分流电路，构成一个周期性单元。单元由贴有压电片的梁段（a段）和未贴压电片的纯基体梁段（b段）交替组成，晶格常数为两段长度之和。 * 关键方程： 研究从压电方程出发，这是连接机械场（应力、应变）与电场（电场强度、电位移）的核心物理方程。对于文中特定的压电片构型（极化方向沿厚度方向，主要考虑一维弯曲振动），采用了简化的压电方程。通过将分流电路的复阻抗 ( Z ) 作为边界条件引入电学方程，并结合力学方程，推导出在连接任意分流电路 ( Z ) 的情况下，压电片的等效弹性模量 ( E_p ) 的解析表达式（见原文公式(3)）。这个表达式是连接电路参数与声子晶体梁力学性能的桥梁，其中 ( s ) 为拉普拉斯算子，( C_p ) 为压电片本身的等效电容。 * 分流电路代入： 将三种具体分流电路的阻抗表达式分别代入上述通用公式(3)。 * 电阻电路（R）： 阻抗 ( Z = R )，代入后得到等效弹性模量 ( E_p ) 与电阻 ( R ) 和振动频率 ( \omega ) 的关系式。 * LC电路： 阻抗 ( Z = sL + 1/(sC) )，代入后得到 ( E_p ) 与电感 ( L )、电容 ( C ) 和频率 ( \omega ) 的关系式。 * LCR电路： 阻抗 ( Z = R + sL + 1/(sC) )，代入后得到 ( E_p ) 与 ( R, L, C, \omega ) 的关系式。
2. 带隙特性计算方法： * 计算方法： 对于所建立的一维周期性梁结构，研究采用传递矩阵法（Transfer Matrix Method, TM） 来计算其弯曲振动的带隙结构。该方法基于梁的弯曲振动方程，结合周期性结构的布洛赫（Bloch）定理，通过计算波在单个周期单元内的传递矩阵及其本征值，来确定弹性波能否在结构中传播（通带）或被禁止（带隙）。 * 参数设置： 研究采用了固定的材料和几何参数进行数值计算，以确保不同电路情况下结果的可比性。具体参数详见原文表1（基体环氧树脂梁参数）和表2（PZT-5H压电片参数）。电路参数（R, L, C）作为变量在计算中变化。
3. 数值计算与结果分析流程： * 将步骤1中得到的、对应于不同电路的等效弹性模量 ( E_p(\omega) ) 表达式，代入到步骤2的传递矩阵法计算框架中。 * 首先，分析等效弹性模量本身 随电路参数和频率的变化。通过绘制 ( E_p ) 的幅值和相位图，直观展示分流电路如何改变压电材料的“有效刚度”。 * 其次，计算有限周期结构的振动传输特性。为了更直观地展示带隙效果（即振动衰减），研究计算了一个包含6个周期单元的有限长梁的振动传输特性曲线。方法是在梁的一端施加单位简谐位移激励，计算另一端在频域上的位移响应，响应幅值（通常取对数）的谷值即对应带隙频率范围。 * 最后，系统性地改变电路参数（如电阻R从0到无穷大变化，电感L变化，固定电容C等），重复上述计算，对比分析不同参数对 ( E_p ) 曲线和传输特性曲线（即带隙）的影响。
研究的主要结果
本研究通过详尽的数值计算，清晰揭示了三种分流电路对声子晶体梁带隙特性的不同影响机制：
1. 电阻分流电路（R电路）的影响： * 对等效弹性模量的影响： 计算结果显示，压电片的等效弹性模量 ( E_p ) 的幅值随分流电阻 ( R ) 的增大，从“短路弹性模量” ( E^E ) 单调增加至“开路弹性模量” ( E^D )。相位变化表明，电阻主要引起能量耗散（对应于复模量的虚部）。但其对模量实部的调节范围有限，仅在两极值之间变化。 * 对梁带隙/振动特性的影响： 传输特性曲线表明，纯电阻电路不能产生新的局域共振带隙。其作用主要体现在两个方面：一是对梁在固有模态频率处的共振峰有一定的抑制（阻尼）作用，这是由于电阻耗散了部分振动能量；二是由于改变了a段（贴片段）与b段（纯梁段）之间的阻抗匹配情况，可以对原本存在的布拉格散射型带隙的位置和宽度产生微调。随着电阻增大，带隙有向高频移动的趋势。但总体而言，调节能力有限。
2. LC分流电路的影响： * 对等效弹性模量的影响： 这是本研究的关键发现之一。当分流电路为纯LC电路（R=0）时，等效弹性模量 ( E_p ) 的频率响应曲线出现奇异点（理论上趋于无穷大）。更重要的是，在两个奇异点之间的频率范围内，( E_p ) 的相位变为180°，这意味着其值变为负实数，即压电片表现出负弹性模量。这是产生局域共振的典型特征。 * 对梁带隙特性的影响： 传输特性曲线清晰地显示，在 ( E_p ) 为负值的频率区间内，梁产生了显著的振动衰减带隙，即局域共振带隙。其物理机理是：梁的振动通过压电效应在LC电路中激励起电磁振荡，而该振荡通过反压电效应反作用于梁，相当于在系统中引入了一个“电磁振子”，与梁的机械振动发生强烈的局部共振耦合，从而抑制了弹性波的传播。带隙的中心频率由LC电路的谐振频率 ( f = 1/(2\pi\sqrt{LC}) ) 决定，因此可以通过改变电感L或电容C来主动、精确地调控带隙的位置。
3. LCR分流电路（引入电阻的LC电路）的影响： * 对等效弹性模量的影响： 在LC电路中引入电阻R后，( E_p ) 频率响应曲线的奇异点消失，幅值的峰值降低，相位变化也变得平缓。电阻R起到了阻尼的作用，削弱了LC电路的振荡强度。 * 对梁带隙特性的影响： 传输特性曲线表明，引入电阻后，局域共振带隙的衰减幅值（即带隙深度）明显减小，因为能量耗散增加，共振效应减弱。但与此同时，带隙的频率范围（带宽）得到了一定的展宽。当电阻增大到一定程度（如文中示例R=100Ω），局域共振效应被过度抑制，带隙几乎消失。这揭示了损耗（电阻）与共振强度（带隙深度）及带宽之间存在权衡关系，为工程设计中根据需求优化带隙特性提供了依据。
研究的结论与价值
本研究系统地论证了通过压电分流电路调控一维声子晶体梁带隙的可行性，并得出了明确结论：
本研究的科学价值在于，它提出并验证了一种非机械式的局域共振声子晶体实现方案。与传统的通过在基体中嵌入硬芯软包层散射体等机械振子方案相比，这种“电磁振子”方案极大地减小了附加质量，更易于在轻质结构上实现。其应用价值尤为突出：电路参数（L、C、R）易于在工程中实现实时、大范围的调节，这使得基于该原理的声子晶体结构能够动态适应环境振动频率的变化，实现自适应减振，在航空航天精密设备隔振、车辆主动悬架、建筑结构减震等领域具有广阔的应用前景。
研究的亮点
其他有价值的内容
研究中对“短路弹性模量” ( E^E ) 和“开路弹性模量” ( E^D ) 的区分和讨论，体现了压电材料在电路边界条件不同时的本征特性差异，这是分析所有压电分流问题的理论基础。文中使用的传递矩阵法（TM）是计算一维周期性结构波传播特性的经典且高效的方法，保证了计算结果的可靠性。