摘要

压电陶瓷谐振频率通常通过阻抗分析仪或具备扫频功能的LCR仪测量。先对夹具进行开路、短路和必要的负载校准，再在规定的机械边界、预紧、温度和小信号激励条件下扫频，同时保存频率、阻抗实部和阻抗虚部。单一模态附近的阻抗幅值局部最小点可作为串联谐振候选频率，局部最大点可作为并联谐振候选频率，但最终判定还要结合相位、导纳、等效电路拟合和所采用的标准定义。

压电执行器装入机构以后，测到的是压电堆、预紧结构、负载、夹具和传感链路共同形成的系统模态，不能直接当成裸陶瓷材料常数。本文给出一套只依赖Python标准库的离线分析代码，用于从真实CSV数据中寻找候选点和电抗过零点。代码自带Butterworth Van Dyke等效电路的教学自检数据，教学数据不代表任何产品实测性能。

1. 压电陶瓷谐振频率测试在解决什么问题

多层压电陶瓷和压电执行器也称压电促动器，本质上是电气域与机械域耦合的器件。阻抗随频率变化时，机械振动会通过压电耦合反映到电端，因此阻抗扫频可用于识别谐振区域、检查装配前后模态变化，并为驱动频率和控制带宽留出工程裕量。

这项测试不能单独回答位移、推力、定位准确度或闭环带宽问题。它回答的是：在给定安装、负载、预紧、温度、激励电平和测量链路下，电端阻抗在哪些频率附近出现与机电模态相关的明显特征。

1.1 材料谐振、器件谐振和结构谐振不是同一个量

• 材料或标准试样测试强调规定几何形状、极化方向和边界条件。
• 压电堆或压电执行器测试包含电极、粘接层、外壳和预紧结构。
• 装机测试还包含负载质量、连接件、柔性铰链和安装刚度。

因此，脱离对象和边界条件只报告一个“谐振频率”没有可比性。

1.2 谐振频率不等于允许驱动频率

允许驱动频率还受电容量、驱动电流、介质损耗、温升、机械应力、控制稳定性和寿命要求限制。靠近结构谐振可能放大位移，也可能带来振动、冲击、发热和失稳风险。实际驱动上限应由器件规格、驱动器能力和系统验证共同确定。

2. 串联谐振和并联谐振怎样理解

在单一模态附近，常用Butterworth Van Dyke等效电路描述压电器件。静态电容C0与运动支路并联，运动支路由Rm、Lm和Cm串联组成。它是局部频段的等效模型，不是器件内部真实机械结构图。

2.1 串联谐振候选点

运动支路中的感抗和容抗在某一频率附近相互抵消时，运动支路接近电阻性，整个器件常出现阻抗幅值局部最小值。实际器件存在损耗、寄生参数和多个模态，所以本文把该点称为串联谐振候选点，而不是未经复核的最终串联谐振频率。

2.2 并联谐振候选点

运动支路与静态电容支路共同作用时，阻抗幅值可能在更高频率附近出现局部最大值。该点可作为并联谐振候选点。候选点是否对应目标振型，需要结合相位、电抗过零、导纳曲线、结构振型或等效电路拟合判断。

2.3 公式应放在明确模型边界内

单模态BVD模型的基本关系如下。公式只用于解释和教学自检，不应直接替代标准规定的材料参数测试。

Z_C0 = 1 / (j * 2 * pi * f * C0)
Z_motion = Rm + j * 2 * pi * f * Lm + 1 / (j * 2 * pi * f * Cm)
Z_total = 1 / (1 / Z_C0 + 1 / Z_motion)

3. 压电阻抗扫频测试系统怎样搭建

3.1 必要硬件

• 可输出频率扫描并测量复数阻抗的阻抗分析仪、网络分析仪配阻抗夹具，或具备足够频段的LCR仪。
• 适合器件电容量和目标频段的四端对、四端子或厂家规定夹具。
• 可记录安装状态、预紧力、负载质量和环境温度的辅助测量工具。
• 具备放电路径和触电防护的安全连接。压电器件断电后仍可能保留电荷。

3.2 校准位置

开路、短路和负载校准的目的，是把线缆和夹具的残余阻抗从测量结果中尽量移除。校准面距离器件端子越远，线缆电感、接触电阻和杂散电容对高频结果的影响越明显。具体校准方式应遵循仪器和夹具手册。

3.3 机械边界

若目标是材料或裸器件参数，应采用对应标准规定的试样和边界。若目标是装机风险，应按实际预紧、安装刚度和负载测试。两类结果用途不同，不能混在同一数据表中比较。

4. 扫频条件应该怎样设置

项目	必须记录的内容	常见风险
被测对象	器件类型、样品编号、极化方向、安装与预紧状态	不同装配状态的数据被混为同一规格
频率范围	起止频率、线性或对数扫频、点数和步进	范围太宽时错把其他模态当目标模态
激励条件	交流电平、直流偏置、积分时间和平均次数	激励过大引起非线性、发热或应力变化
电气连接	夹具型号、线缆、校准方式和校准面	寄生参数主导高频曲线
机械条件	自由、夹持或装机状态，预紧力和负载质量	结构模态被误认为材料谐振
环境条件	温度、稳定时间和测试顺序	温漂和自热造成重复性下降
原始数据	频率、阻抗实部、阻抗虚部，必要时保存相位和导纳	只留截图，无法复算和追溯

4.1 先宽范围粗扫，再对单一模态细扫

第一次测试可以用较宽频率范围识别主要特征，再把一个目标模态放入更窄窗口，以更密的点间隔重新扫描。代码中的局部最小值和局部最大值算法假设分析窗口主要包含一个目标模态。窗口内存在多个强模态时，应分段分析。

4.2 扫频分辨率没有统一固定值

分辨率应根据峰宽、目标不确定度和仪器稳定性确定。若相邻频点间隔大于特征峰宽，候选频率只能落在稀疏采样点上。增加点数不能补偿错误的机械边界、校准或过低的信噪比。

5. 软件、编译与测试条件

项目	本文代码条件
操作系统	Windows 11 64位，其他支持Python 3.12的平台可运行
Python	Python 3.12，代码仅使用标准库
输入文件	UTF-8 CSV，列名为frequency_hz、resistance_ohm、reactance_ohm
输出文件	impedance_summary.json和impedance_analyzed.csv
硬件连接	代码不连接仪器，只分析仪器导出的复数阻抗数据
验证状态	自检模式已做离线运行验证；真实硬件数据仍需按仪器和夹具验证

6. Python阻抗扫频分析代码

下面代码先限制分析频率窗口，再寻找阻抗幅值的局部最小值、局部最大值和电抗过零点。它输出的是候选点，不负责判定振型，也不计算材料机电耦合系数。

from __future__ import annotations

import argparse
import csv
import json
import math
from pathlib import Path

REQUIRED_COLUMNS = {
    "frequency_hz",
    "resistance_ohm",
    "reactance_ohm",
}

def read_impedance_csv(path: Path) -> list[dict[str, float]]:
    rows: list[dict[str, float]] = []
    with path.open("r", encoding="utf-8-sig", newline="") as handle:
        reader = csv.DictReader(handle)
        if reader.fieldnames is None:
            raise ValueError("CSV has no header")
        missing = REQUIRED_COLUMNS.difference(reader.fieldnames)
        if missing:
            raise ValueError(f"missing columns: {sorted(missing)}")

        previous_frequency = -math.inf
        for line_number, raw in enumerate(reader, start=2):
            try:
                frequency = float(raw["frequency_hz"])
                resistance = float(raw["resistance_ohm"])
                reactance = float(raw["reactance_ohm"])
            except (TypeError, ValueError) as exc:
                raise ValueError(f"invalid number at line {line_number}") from exc

            if not all(math.isfinite(v) for v in (frequency, resistance, reactance)):
                raise ValueError(f"non-finite value at line {line_number}")
            if frequency <= 0 or frequency <= previous_frequency:
                raise ValueError("frequency_hz must be positive and strictly increasing")
            previous_frequency = frequency

            magnitude = math.hypot(resistance, reactance)
            phase_degree = math.degrees(math.atan2(reactance, resistance))
            rows.append({
                "frequency_hz": frequency,
                "resistance_ohm": resistance,
                "reactance_ohm": reactance,
                "magnitude_ohm": magnitude,
                "phase_degree": phase_degree,
            })

    if len(rows) < 7:
        raise ValueError("at least 7 frequency points are required")
    return rows

def bvd_impedance(frequency: float) -> complex:
    # Teaching values for code self-test only, not product parameters.
    r_m = 20.0
    l_m = 0.032
    c_m = 5e-9
    c_0 = 50e-9
    angular_frequency = 2.0 * math.pi * frequency
    z_c0 = 1.0 / (1j * angular_frequency * c_0)
    z_motion = (
        r_m
        + 1j * angular_frequency * l_m
        + 1.0 / (1j * angular_frequency * c_m)
    )
    return 1.0 / (1.0 / z_c0 + 1.0 / z_motion)

def create_self_test_csv(path: Path) -> None:
    path.parent.mkdir(parents=True, exist_ok=True)
    with path.open("w", encoding="utf-8", newline="") as handle:
        writer = csv.writer(handle)
        writer.writerow([
            "frequency_hz",
            "resistance_ohm",
            "reactance_ohm",
        ])
        for frequency in range(10000, 16001, 10):
            impedance = bvd_impedance(float(frequency))
            writer.writerow([
                frequency,
                f"{impedance.real:.9f}",
                f"{impedance.imag:.9f}",
            ])

def local_extreme_indices(values: list[float], mode: str) -> list[int]:
    indices: list[int] = []
    for index in range(1, len(values) - 1):
        left, current, right = values[index - 1:index + 2]
        if mode == "minimum" and current <= left and current <= right:
            indices.append(index)
        if mode == "maximum" and current >= left and current >= right:
            indices.append(index)
    return indices

def reactance_zero_crossings(rows: list[dict[str, float]]) -> list[float]:
    crossings: list[float] = []
    for first, second in zip(rows, rows[1:]):
        x1 = first["reactance_ohm"]
        x2 = second["reactance_ohm"]
        if x1 == 0.0:
            crossings.append(first["frequency_hz"])
        elif x1 * x2 < 0.0:
            f1 = first["frequency_hz"]
            f2 = second["frequency_hz"]
            fraction = -x1 / (x2 - x1)
            crossings.append(f1 + fraction * (f2 - f1))
    return crossings

def analyze(
    all_rows: list[dict[str, float]],
    frequency_min: float | None,
    frequency_max: float | None,
) -> tuple[list[dict[str, float]], dict[str, object]]:
    rows = [
        row for row in all_rows
        if (frequency_min is None or row["frequency_hz"] >= frequency_min)
        and (frequency_max is None or row["frequency_hz"] <= frequency_max)
    ]
    if len(rows) < 7:
        raise ValueError("analysis window must contain at least 7 points")

    magnitudes = [row["magnitude_ohm"] for row in rows]
    minimum_indices = local_extreme_indices(magnitudes, "minimum")
    maximum_indices = local_extreme_indices(magnitudes, "maximum")
    if not minimum_indices or not maximum_indices:
        raise ValueError(
            "no complete local minimum and maximum found; refine the sweep window"
        )

    minimum_index = min(minimum_indices, key=lambda i: magnitudes[i])
    maximum_index = max(maximum_indices, key=lambda i: magnitudes[i])
    minimum_row = rows[minimum_index]
    maximum_row = rows[maximum_index]
    crossings = reactance_zero_crossings(rows)

    warnings: list[str] = []
    if len(minimum_indices) > 1 or len(maximum_indices) > 1:
        warnings.append("multiple local extrema found; the window may contain several modes")
    if maximum_row["frequency_hz"] < minimum_row["frequency_hz"]:
        warnings.append("maximum-impedance candidate is below minimum candidate; verify mode selection")

    summary: dict[str, object] = {
        "data_status": "candidate frequencies from supplied data",
        "point_count": len(rows),
        "frequency_window_hz": [
            rows[0]["frequency_hz"],
            rows[-1]["frequency_hz"],
        ],
        "minimum_impedance_candidate": minimum_row,
        "maximum_impedance_candidate": maximum_row,
        "reactance_zero_crossings_hz": crossings,
        "warnings": warnings,
        "interpretation_boundary": (
            "Confirm resonance definitions with phase, admittance, mode shape, "
            "equivalent-circuit fitting, and the applicable standard."
        ),
    }
    return rows, summary

def write_analyzed_csv(path: Path, rows: list[dict[str, float]]) -> None:
    columns = [
        "frequency_hz",
        "resistance_ohm",
        "reactance_ohm",
        "magnitude_ohm",
        "phase_degree",
    ]
    with path.open("w", encoding="utf-8", newline="") as handle:
        writer = csv.DictWriter(handle, fieldnames=columns)
        writer.writeheader()
        writer.writerows(rows)

def main() -> None:
    parser = argparse.ArgumentParser(
        description="Find candidate extrema in piezoelectric impedance sweep data."
    )
    source = parser.add_mutually_exclusive_group(required=True)
    source.add_argument("--input", type=Path, help="input impedance CSV")
    source.add_argument(
        "--self-test",
        action="store_true",
        help="create and analyze illustrative BVD data",
    )
    parser.add_argument("--f-min", type=float, default=None)
    parser.add_argument("--f-max", type=float, default=None)
    parser.add_argument("--output-dir", type=Path, default=Path("."))
    args = parser.parse_args()

    args.output_dir.mkdir(parents=True, exist_ok=True)
    if args.self_test:
        input_path = args.output_dir / "bvd_self_test.csv"
        create_self_test_csv(input_path)
        frequency_min = 10000.0
        frequency_max = 16000.0
    else:
        input_path = args.input
        frequency_min = args.f_min
        frequency_max = args.f_max

    rows = read_impedance_csv(input_path)
    analyzed_rows, summary = analyze(rows, frequency_min, frequency_max)
    write_analyzed_csv(args.output_dir / "impedance_analyzed.csv", analyzed_rows)
    with (args.output_dir / "impedance_summary.json").open(
        "w", encoding="utf-8"
    ) as handle:
        json.dump(summary, handle, ensure_ascii=False, indent=2)

    print(json.dumps(summary, ensure_ascii=False, indent=2))

if __name__ == "__main__":
    main()

6.1 自检命令

python impedance_sweep_analyzer.py --self-test --output-dir self_test_output

自检模式根据代码中明确给出的BVD教学参数生成CSV，再执行同一套分析流程。它只验证程序流程，不验证测量仪器，也不代表压电陶瓷或压电执行器的实测性能。

6.2 真实CSV输入格式

frequency_hz,resistance_ohm,reactance_ohm
12000,101.2,-850.5
12010,98.7,-812.4
12020,95.1,-770.8

真实数据分析时应先通过粗扫确定目标模态窗口，然后运行：

python impedance_sweep_analyzer.py --input actual_impedance.csv --f-min 10000 --f-max 16000 --output-dir result

上面的频率窗口只是命令格式示例，不能直接作为任何器件的测试范围。

7. 可验证数据与测试方法

7.1 本文可验证的只有程序自检

项目	条件	数据状态	用途
输入模型	代码内BVD教学参数，10 kHz到16 kHz，步进10 Hz	理论生成数据	验证CSV、极值和过零分析流程
运行环境	Python 3.12.13，601个频率点	已验证	确认代码可编译并完成自检流程
阻抗局部最小候选	12580 Hz，阻抗幅值约19.8868 ohm	理论自检结果	检查局部最小值提取
阻抗局部最大候选	13200 Hz，阻抗幅值约2927.2190 ohm	理论自检结果	检查局部最大值提取
电抗过零候选	约12586.226 Hz和13192.234 Hz	线性插值自检结果	检查符号变化和过零插值
产品性能	没有连接真实阻抗仪和压电器件	未测试	不得形成产品规格或对外性能结论

7.2 真实测试至少需要三类重复

1. 同一连接不拆装连续重复，用于观察仪器短期重复性。
2. 重新连接或重新装夹重复，用于观察接触和夹具影响。
3. 不同样品重复，用于评估样品差异，样本数和统计口径必须报告。

结果建议同时报告候选频率、阻抗幅值、相位、电抗过零点、扫频分辨率、重复次数及离散程度。没有原始文件、仪器编号和校准记录的单次截图，不应作为可验证性能数据。

7.3 曲线判断顺序

1. 先检查频率是否严格递增，单位是否一致。
2. 确认目标窗口内是否只有一个主要模态。
3. 查看阻抗幅值局部最小点和局部最大点。
4. 查看相位和电抗过零点是否与特征相互支持。
5. 比较重复扫频、重新装夹和温度变化后的漂移。
6. 需要材料参数时，按对应标准和振型公式计算，而不是套用通用公式。

8. 常见误判与排查方法

8.1 把全频段最小值直接当串联谐振

电容量背景、线缆寄生和其他模态可能让全频段极值落在错误位置。应先选定单一模态窗口，再结合相位和导纳复核。

8.2 忽略夹具和预紧

压电执行器的预紧和安装刚度会改变系统模态。若装配前后测试目的不同，应分别建档，不要用裸件结果替代装机结果。

8.3 只保存阻抗幅值

仅有阻抗幅值会丢失复数方向信息。优先保存阻抗实部和虚部，幅值与相位可由原始复数阻抗复算。

8.4 用大信号结果替代小信号阻抗

激励电平变化可能引起材料非线性、自热和机械状态变化。小信号阻抗与高压、大位移动态响应属于不同测试，不能直接替代。

8.5 从两个候选频率直接计算耦合系数

不同振型、几何形状和标准可能采用不同公式。没有确认试样模式和适用标准时，本文代码不计算机电耦合系数。

9. 常见问题

9.1 用普通万用表能测压电陶瓷谐振频率吗

通常不能。万用表无法提供所需的频率扫描和复数阻抗测量。应使用覆盖目标频段并适配器件阻抗范围的阻抗分析仪、LCR仪或合适的网络分析方案。

9.2 测一个固定频率的电容量能确定谐振频率吗

不能。固定频率电容量只反映该测试条件下的等效电容信息，不能代替跨频率的阻抗幅值和相位扫描。

9.3 阻抗最小点一定是串联谐振频率吗

只能先作为候选点。损耗、寄生参数、多模态和扫频窗口都会影响极值。正式判定应结合相位、导纳、等效电路和所采用标准。

9.4 压电陶瓷加预紧后谐振频率会变化吗

可能变化。预紧会改变接触状态、等效刚度和边界条件。报告必须写明预紧方法、力值或可追溯设置，以及是否重新装夹。

9.5 扫频激励电压越高越好吗

不是。激励应在仪器、夹具和器件允许范围内，并与测试目的一致。过高激励可能引起非线性、自热、机械冲击或仪器过载。

9.6 压电执行器谐振频率等于闭环带宽吗

不等于。闭环带宽还受传感器、驱动器、采样、延迟、控制器和其他结构模态限制，通常需要单独的频率响应和稳定性验证。

9.7 能用串联和并联谐振频率直接计算材料参数吗

只有在试样几何、振动模式、极化方向、边界条件和公式都符合对应标准时才可以。装机压电执行器的系统模态通常不能直接反推出材料常数。

9.8 本文Python代码能直接控制阻抗分析仪吗

不能。代码只读取CSV并离线分析，不包含VISA通信、仪器触发、校准控制或高压联锁。需要自动采集时，应按具体仪器编程手册另行实现和验证。

10. 结论

压电陶瓷谐振频率测试的关键，不是从曲线上挑出一个最小值，而是把被测对象、机械边界、校准面、扫频条件和原始复数阻抗同时保存。对于单一模态，阻抗局部最小值和局部最大值可以分别作为串联、并联谐振候选点，但必须结合相位、导纳、等效电路和标准定义复核。

压电执行器装机后测得的是系统模态。谐振频率也不等于允许驱动频率或闭环带宽。本文Python代码适合做CSV数据质量检查和候选点提取，不负责产品性能判定，更不能替代实际仪器校准与机械振型验证。

参考资料

1. IEEE Std 176-1987，IEEE Standard on Piezoelectricity，Institute of Electrical and Electronics Engineers，1988。DOI：https://doi.org/10.1109/IEEESTD.1988.79638
2. IEC 60483:1976，Guide to dynamic measurements of piezoelectric ceramics with high electromechanical coupling，International Electrotechnical Commission。标准页面：https://webstore.iec.ch/en/publication/2185
3. Keysight Technologies，Impedance Measurement Handbook，Application Note 5950-3000。资料页面：https://www.keysight.com/us/en/assets/7018-06840/application-notes/5950-3000.pdf
4. Python Software Foundation，csv - CSV File Reading and Writing，Python 3.12文档。https://docs.python.org/3.12/library/csv.html
5. Python Software Foundation，math - Mathematical functions，Python 3.12文档。https://docs.python.org/3.12/library/math.html

说明：标准可能需要通过正式渠道购买或查阅，并可能存在更新、替代或撤销状态。正式试验应核对企业当前采用版本、客户规范以及具体仪器和夹具手册。