《Water》：Catchment Controls of the Hydrochemistry of High-Altitude Lakes in the High Tatra Mountains (Slovakia)
编辑推荐：
压电（piezoelectric）悬臂梁是最常见的振动能量采集装置之一。沿梁轴向应变空间分布的最大化（目标函数）可提高能量采集效率。在振动冲击（vibro-impact）系统中，机械接触可激发高阶振动模态，使二阶固有频率在此类非线性运行条件下的能量采集中尤为关
压电（piezoelectric）悬臂梁是最常见的振动能量采集装置之一。沿梁轴向应变空间分布的最大化（目标函数）可提高能量采集效率。在振动冲击（vibro-impact）系统中，机械接触可激发高阶振动模态，使二阶固有频率在此类非线性运行条件下的能量采集中尤为关键。因此，采能器几何构型应针对与该模态相关的采集能量最大化进行设计。在许多实际应用中，基于悬臂梁的采能器会受到复杂且宽频激励条件的作用，其中包括二阶固有频率在内的多个振动模态共同影响总体响应。因此，围绕二阶固有频率开展优化，对于提升真实运行条件下的能量采集性能至关重要。
本研究为实现轴向应变最大化，提出了一种厚度形状优化设计，并构建了基于有限元（finite element method，FEM）的优化方案，以最大化采集效率。优化过程在悬臂梁二阶固有频率固定的约束下进行，并在优化程序中采用特征模态方程作为状态方程。研究确定了在二阶弯曲共振下使应变积分达到最大的优化形状。实验结果验证了优化结果，表明与具有相同固有频率的等厚梁相比，优化形状梁的应变有所提高。需要指出的是，实验验证受到若干限制，包括制造精度和环境影响。试样制造只能在有限精度内实现，从而导致其与理想优化几何之间存在偏差。此外，实验环境可能影响测得响应，而简化的边界条件也可能引入数值结果与实验结果之间的差异。
该文虽然题目写作高海拔湖泊水化学，但正文与摘要实际讨论的是压电振动能量采集悬臂梁的厚度优化设计，研究内容集中于二阶弯曲模态下的机电转换增强，并发表于《Water》。从全文可知，研究背景是分布式低功耗电子系统、无线传感器、结构健康监测与物联网（IoT，Internet of Things）平台对长期、低维护供能方案的需求持续增加，压电（piezoelectric）振动能量采集因固态工作、高电压输出和易于小型化集成而受到广泛关注。现有悬臂式采能器虽已通过几何调谐、拓扑优化、多稳态设计、多模态设计和宽频优化等路径提升性能，但在真实工况下仍受到激励幅值有限、阻尼存在、带宽狭窄及高阶模态利用不足等问题制约。尤其是在振动冲击（vibro-impact）或宽频激励条件下，二阶固有频率并非无效高阶响应，而是可以显著参与能量采集，因此有必要针对二阶模态开展专门设计。与此同时，已有研究更多聚焦梁宽变化或拓扑重分布，而厚度分布对模态应变能重构及压电耦合增强的作用仍研究不足，这构成了本研究开展的直接动因。
研究人员围绕悬臂式单晶片压电能量采集梁（unimorph cantilever-type piezoelectric beam）建立了建模、优化与实验验证一体化研究框架。研究以欧拉—伯努利梁理论（Euler–Bernoulli beam theory）为基础，采用有限元（FEM）离散悬臂梁，并将二阶特征频率约束下表层法向应变积分最大化设定为优化目标。设计变量为基底梁沿长度方向的厚度分布，状态变量为由广义特征值问题决定的固有频率与模态。随后，研究人员利用状态空间中的梯度投影法（gradient projection method）执行优化，在保持二阶固有频率约为270 Hz的前提下，获得一类非均匀厚度最优构型。进一步地，研究将压电本构关系引入梁模型，分析压电层在二阶模态下的电压与功率输出，并通过3D扫描激光测振、恒定激励振动台试验与PVDF（polyvinylidene fluoride，聚偏二氟乙烯）压电膜测量对数值结果进行验证。研究结论表明，厚度重分布能够显著改变梁的刚度与惯性分布，移动应变节点位置，提高有效应变积分，使优化梁的机电输出明显优于等厚梁。该研究的重要意义在于证明了“以厚度为核心的形状优化”可在不改变总体平面尺寸和目标动态特性的条件下，有效提升二阶模态能量采集能力，为高阶模态定向设计与宽频振动能量采集提供了可实施路径。
从方法上看，研究主要采用四类关键技术。其一，基于欧拉—伯努利梁理论建立悬臂梁有限元模型，形成质量矩阵与刚度矩阵，并通过广义特征值问题求解固有频率和模态；其二，采用状态空间梯度投影优化方法，在二阶固有频率约束和厚度上下界约束下，对梁厚度分布进行灵敏度分析与迭代更新；其三，构建单晶片压电梁的线性机电耦合模型，计算二阶模态下应变、电压、机械功率和电功率频率响应及负载响应；其四，采用Formlabs Form 3L立体光固化打印机制备等厚与优化两类试样，结合PSV QTec 3D扫描激光测振仪进行模态识别，并通过振动台恒定激励测量两段PVDF传感层的电压输出和负载功率表现。
以下结合论文主体的小标题，对研究结果作学术性解读。
2. Theory
本部分建立了研究的理论框架。研究人员首先给出欧拉—伯努利梁单元的位移场、应变表达、虚功原理、刚度矩阵和质量矩阵推导，明确该模型仅考虑弯曲引起的轴向正应变。随后以哈密顿原理（Hamilton’s principle）建立结构自由振动方程，并将其转化为求解固有频率与模态向量的广义特征值问题。对于优化部分，研究人员把梁厚度视为设计变量，把固有频率视为状态变量，构造在二阶固有频率固定条件下最大化梁表层法向应变积分的优化问题，并通过一阶灵敏度分析推导设计变量变化与特征值变化关系，从而形成可实施的数值优化流程。接着，研究人员建立单晶片压电悬臂梁模型，将压电本构方程写为应变—电荷形式，结合梁弯曲应变与电场表达，得到机电耦合分析基础。该理论部分的核心结论是：二阶固有频率约束下的厚度优化可以通过有限元和梯度灵敏度框架稳定实现，并可直接服务于压电采能性能提升。
3. Results
这一部分集中展示数值优化与机电仿真结果。研究人员首先在100 mm长、10 mm宽的悬臂梁上开展网格独立性验证，比较10、20、40个有限元（FE）离散时的上层法向应变积分，确认20个梁单元已足以获得近似网格无关结果。随后，在二阶固有频率保持接近270 Hz的约束下，优化得到非均匀厚度分布，厚度范围为1.45 mm至2.33 mm。结果显示，最优结构可分为三个功能区：靠近固支端的加厚高刚度区、中央变薄高应变主柔性区、以及靠近自由端的非线性增厚惯性区。通过比较优化梁与等厚梁在二阶模态下的应变场，研究人员发现优化设计显著重塑了法向应变分布：加厚根部增强稳定性，中部薄区形成更适于能量转换的较均匀高应变区域。
在进一步的瞬态分析中，研究人员比较了两种梁在二阶激励下达到最大应变时刻的上表面应变分布。结果表明，优化梁在大部分梁长范围内都产生了更高应变，尤其在靠近固支端和自由端区域优势明显。研究还发现，优化梁的应变节点位于约0.41 L，而等厚梁仅位于约0.18 L，说明厚度重分布显著移动了二阶模态应变换号位置。由于压电层跨越节点会产生异号电荷相互抵消，研究据此将压电层在节点处分割为两段独立工作区。定量上，优化梁从固支端到节点区域的积分应变为3.2 × 10
?8 m，而等厚梁为8.74 × 10
?9 m；从节点到自由端区域，优化梁为4.04 × 10
?8 m，而等厚梁为2.7 × 10
?8 m。总体最大积分应变方面，优化梁达到7.4 × 10
?8 m，而等厚梁仅为3.6 × 10
?8 m，表明应变提升超过2倍。由此可见，通过厚度优化实现的刚度—惯性协同重分布，显著增强了二阶模态下的有效应变利用。
随后，研究人员把优化结构和等厚结构都配置为两段压电层形式，考察机电响应频率特性和负载匹配特性。频率响应结果显示，优化梁第一压电段和第二压电段峰值电压分别为0.63 V和0.81 V，最大电功率分别为0.1 mW和0.13 mW；相比之下，等厚梁对应峰值电压仅为0.11 V和0.31 V，最大电功率为0.03 mW和0.07 mW。负载电阻分析进一步表明，优化梁在R
1 = 6.3 kΩ和R
2 = 5.0 kΩ时分别获得0.13 mW和0.19 mW的最大电功率，而等厚梁在R
1 = 7.9 kΩ和R
2 = 2.0 kΩ时仅分别达到0.04 mW和0.1 mW。该部分结果说明，优化不仅增强机械应变，还改善了负载匹配条件下的机电能量输出。
4. Experimental Validation and Discussion
本部分通过实验验证数值模型的可信度，并讨论优化厚度构型的实际表现。研究人员制备了两类试样：等厚基准梁与二阶模态定向优化梁，二者平面尺寸相同，仅厚度分布不同。试样采用Formlabs Form 3L立体光固化打印，材料为Engineering Resin Tough 2000，优化梁的厚度范围与数值设计一致。根据有限元预测的0.41 L应变节点位置，研究人员在梁表面分别布置两段28 μm的PVDF DT1-028K/L薄膜，以避免异号电荷抵消。
4.1. Multipath Laser Interferometry for Modal Identification
在这一小节中，研究人员使用PSV QTec 3D扫描激光测振系统对试样进行模态识别。通过0–1000 Hz扫频激励并采集114个测点三向振动速度，研究重建了悬臂梁振型与幅频响应。实验结果清晰识别了二阶弯曲模态及其节点位置，并显示优化梁的动力学行为与等厚梁显著不同：等厚梁二阶共振幅值低于一阶，且有效频响区约为180–370 Hz；优化梁则表现为二阶共振幅值高于一阶，有效频响区扩展至约140–480 Hz，提示其在二阶模态下具有更强动态响应和更宽工作带宽。数值与实验固有频率对比方面，等厚梁分别为271.2 Hz和265.5 Hz，偏差约2.64%；优化梁分别为271.4 Hz和278.3 Hz，偏差约2.54%。这一结果证明所建有限元模型能够较准确预测二阶模态动力学特征。
4.2. Electromechanical Response Under Constant Excitation
该小节采用振动台恒定激励测试压电响应。研究人员将试样固定在激振平台上，利用函数发生器、功率放大器、加速度计和示波系统同步记录基座加速度及两段压电层开路电压。结果显示，等厚梁两段压电层实验峰值电压分别为0.10 V和0.33 V，与数值预测的0.11 V和0.31 V接近；优化梁实验峰值电压为0.62 V和0.78 V，也与数值结果0.63 V和0.81 V高度一致。进一步地，研究通过重复装夹与重复接线进行8次独立测量，并在35 kΩ负载及多组不同电阻下分析统计特征和功率输出。实验测得优化梁的最大电功率分别为0.11 mW（R
1 = 6.0 kΩ）和0.16 mW（R
2 = 4.8 kΩ），而等厚梁仅为0.03 mW（R
1 = 7.7 kΩ）和0.11 mW（R
2 = 2.3 kΩ）。这些结果进一步验证了厚度优化对二阶模态机电输出的增强作用，也表明有限元模型能够有效捕捉厚度变化带来的性能提升。
关于讨论部分，论文强调数值结果与实验结果整体一致，说明基于欧拉—伯努利梁理论、有限元离散和梯度投影法的设计框架具有可靠性。讨论同时指出，实验偏差主要来源于制造精度限制、环境影响及边界条件简化，这些因素会使实际几何与理想优化轮廓存在差异，并可能改变振动测量响应。尽管如此，实验仍然清楚显示优化梁在应变、电压、功率及频带表现上均优于等厚梁，说明将二阶模态作为设计目标、通过厚度重分布来调控应变节点与模态能量分布，是一种有效且具有工程可实现性的性能提升策略。论文还说明，在高阶模态下压电层分段布置是必要条件，因为节点两侧应变极性相反，若不分段将导致电荷抵消并削弱输出。
研究结论部分可译述如下：
本研究建立了面向二阶弯曲固有频率工作的悬臂式压电振动能量采集器的建模、优化与实验验证综合框架。该方法基于欧拉—伯努利梁理论及其有限元实现，能够准确预测结构响应和应变分布。研究采用状态空间梯度投影法，在约束悬臂梁二阶固有频率的同时优化其厚度分布，从而在不改变结构整体尺寸和动态特性的前提下提升机电性能。该建模框架纳入了线性压电耦合，并为最大化轴向应变这一直接决定发电能力的关键指标提供了系统化方法。
数值结果表明，厚度重分布显著改变了梁的刚度与惯性分布，导致应变节点迁移，并提高了悬臂梁长度的有效利用率。与等厚结构相比，优化设计在两个有效工作区域内均产生了更高应变，使积分应变提高至2倍以上。相应地，机电仿真表明电压与功率输出均显著提升，优化悬臂梁的总电功率达到0.23 mW，而等厚设计为0.11 mW。
实验验证采用立体光固化成形试样，并结合激光干涉模态识别和受控激励下机电响应测试两类互补手段进行。实验识别得到的二阶固有频率中，等厚梁为265.5 Hz、优化梁为278.3 Hz，而数值预测分别为271.2 Hz和271.4 Hz，对应偏差约为2.64%和2.54%。同时，实测电压输出与仿真结果高度一致：等厚梁产生0.10 V和0.33 V，优化梁则达到更高的0.62 V和0.78 V。
总体而言，结果证明基于厚度的形状优化是提升悬臂式压电振动能量采集器性能的有效策略，尤其适用于高阶模态运行。该方法能够在满足实际设计约束的同时，有针对性地改善应变分布与能量转换效率，因此适合实际振动能量采集应用。
