压电执行器安装偏心怎么测？倾斜、横向力与Python误差分析

压电执行器安装偏心不能只靠目测判断。可在执行器前后两个测点布置位移传感器，在低速小行程条件下记录两点位移差，再用位移差除以测点间距估算倾斜角；同时用力传感器或导向机构反力检查横向载荷。该方法能区分轴向位移、安装倾斜和偏载，但不能单独证明内部陶瓷缺陷或整机定位精度。

本文给出测量基准、装夹流程、误差预算和 Python 离线分析示例。示例数据为教学输入，不代表任何产品实测指标；代码只读取 CSV，不连接或控制真实压电驱动器。

从驱动器、执行器、导向夹具到双点位移和力传感器，说明测量基准与数据流。

1. 为什么要测安装偏心

压电执行器通常沿轴向输出位移。若推杆端面与负载基准不平行，或夹具中心线与执行器轴线不重合，位移会被分解为轴向分量和转动分量。转动会使导向件承受横向力，表现为摩擦变化、回程差增大、位移传感器读数不一致，甚至出现卡滞。

在纳米定位装配中，安装误差往往与压电陶瓷本身的迟滞、蠕变和温漂叠加。先把机械误差单独测出来，才能避免把装配问题误判为控制算法问题。

2. 术语定义和直接结论

用执行器轴线、负载基准、两测点和横向力箭头区分三个概念。

3. 测量原理与适用边界

设前后两个位移测点的读数分别为 d1 和 d2，测点间距为 L。小角度条件下，倾斜角约等于 d1 与 d2 的差除以 L。若位移单位为微米、间距单位为毫米，结果的数值单位可直接表示为毫弧度。中心位移可取两点读数的平均值，用于与单点传感器结果交叉核对。

该近似要求倾斜角较小、两传感器量程和方向已校准、测点刚性连接且读数同步。若存在明显转动、柔性结构变形或大角度摆动，应使用刚体位姿拟合或多点几何标定，不能继续套用小角度公式。

横向力的测量需要明确结构模型。若用导向件位移乘以横向刚度估算，刚度必须通过独立标定得到；未知摩擦、间隙和接触预紧会使该估算失效。

4. 测试架构和装夹要点

推荐把双点位移、力信号和驱动命令统一时间戳，便于定位机械误差来源。

• 先定义机械基准面和执行器作用线，再安装传感器，避免以夹具外形代替测量基准。
• 前后测点应位于同一刚性部件，传感器测头预载和方向保持一致。
• 使用低速、小幅值往返动作检查卡滞，确认没有撞限位和侧向刮擦。
• 力传感器若只测轴向力，不得把轴向读数直接当作横向力。

5. 测试流程

1. 记录样品编号、夹具版本、传感器型号、校准日期和环境温度。
2. 断电安装，确认高压驱动器处于放电完成状态，并设置机械限位。
3. 用量块或标准位移检查两个传感器的零点、比例和方向。
4. 在空载或低负载下做三次往返，确认位移差没有明显突变。
5. 逐步增加轴向负载，每个工况等待读数稳定，再记录双点位移和力。
6. 计算倾斜角、中心位移、峰峰值和重复性；改变装夹方向做交叉验证。
7. 若需要估算横向力，先完成导向刚度和摩擦的独立标定，再报告模型假设。

断电装夹、传感器校准、低速预检、负载扫描和交叉复测构成闭环流程。

6. 软件、硬件、编译及测试条件

7. Python离线分析代码

代码读取 CSV 并输出每个工况的位移差、倾斜角、中心位移和轴向力。横向力不在代码中虚构计算，因为没有提供经过标定的横向刚度模型。代码只做离线分析，不控制驱动器。

python alignment_error.py --self-test
python alignment_error.py alignment_data.csv

from __future__ import annotations
import argparse
import csv
import json
import math
import statistics
from pathlib import Path

def number(value: str, field: str, line: int) -> float:
    try:
        result = float(value)
    except (TypeError, ValueError) as exc:
        raise ValueError(f"invalid {field} at line {line}") from exc
    if not math.isfinite(result):
        raise ValueError(f"non-finite {field} at line {line}")
    return result

def load(path: Path):
    with path.open("r", encoding="utf-8-sig", newline="") as handle:
        reader = csv.DictReader(handle)
        required = {"case_id", "front_displacement_um", "rear_displacement_um", "sensor_spacing_mm", "axial_force_n"}
        missing = required - set(reader.fieldnames or [])
        if missing:
            raise ValueError(f"missing columns: {sorted(missing)}")
        rows = []
        for line, raw in enumerate(reader, start=2):
            row = {"case_id": raw["case_id"].strip()}
            for field in required - {"case_id"}:
                row[field] = number(raw[field], field, line)
            if row["sensor_spacing_mm"] <= 0:
                raise ValueError(f"sensor_spacing_mm must be positive at line {line}")
            rows.append(row)
    if not rows:
        raise ValueError("CSV contains no data rows")
    return rows

def analyze(rows):
    result = []
    for row in rows:
        difference_um = row["front_displacement_um"] - row["rear_displacement_um"]
        tilt_mrad = difference_um / row["sensor_spacing_mm"]
        center_um = (row["front_displacement_um"] + row["rear_displacement_um"]) / 2.0
        result.append({"case_id": row["case_id"], "difference_um": difference_um, "tilt_mrad": tilt_mrad, "center_um": center_um, "axial_force_n": row["axial_force_n"]})
    tilts = [item["tilt_mrad"] for item in result]
    return {"samples": len(result), "results": result, "tilt_mean_mrad": statistics.fmean(tilts), "tilt_peak_to_peak_mrad": max(tilts) - min(tilts)}

def self_test():
    rows = [{"case_id": "teach-01", "front_displacement_um": 10.20, "rear_displacement_um": 10.00, "sensor_spacing_mm": 20.00, "axial_force_n": 5.0}]
    output = analyze(rows)
    assert output["samples"] == 1
    assert abs(output["results"][0]["tilt_mrad"] - 0.01) < 1e-12
    assert abs(output["results"][0]["center_um"] - 10.10) < 1e-12
    print("SELF_TEST_OK")
    print(json.dumps(output, ensure_ascii=False, indent=2))

def main():
    parser = argparse.ArgumentParser(description="Analyze piezo actuator installation alignment")
    parser.add_argument("csv_path", nargs="?", type=Path)
    parser.add_argument("--self-test", action="store_true")
    args = parser.parse_args()
    if args.self_test:
        self_test()
    elif args.csv_path is None:
        parser.error("csv_path is required")
    else:
        print(json.dumps(analyze(load(args.csv_path)), ensure_ascii=False, indent=2))

if __name__ == "__main__":
    main()

示例计算中，前点 10.20 微米、后点 10.00 微米、间距 20.00 毫米，位移差为 0.20 微米，倾斜角为 0.010 毫弧度，中心位移为 10.10 微米。这些数字只用于程序自检。

8. 教学结果、误差预算与判读

主要不确定度来源包括位移传感器分辨率和线性度、测点间距测量误差、同步偏差、夹具热膨胀、接触预载、导向摩擦和数据滤波。报告中应逐项给出量级或校准证据，不能只保留一个倾斜角数字。

建议同时绘制倾斜角、中心位移和轴向力，标注教学数据与未测横向力边界。

9. 失效场景和排查顺序

1. 两个测点差值恒定但换向突变，优先检查传感器方向、采样同步和机械间隙。
2. 增加轴向负载后倾斜角增大，检查端面平行度、预紧力和导向刚度。
3. 位移正常但力信号漂移，检查力传感器零点、温漂和安装预载。
4. 重复装夹结果差异大，记录夹具定位面、螺钉顺序和拧紧力矩。
5. 出现卡滞或异常噪声时立即停机，卸载并确认驱动器放电后再拆线。

10. 常见问题

10.1 单个位移传感器能测安装偏心吗

只能观察轴向位移，难以区分整体移动和倾斜。至少需要双点位移，或配合独立的几何标定。

10.2 位移差除以间距得到的是什么

在小角度、刚体近似下得到倾斜角的近似值。它不是位置准确度，也不是横向力。

10.3 为什么不能用倾斜角直接算横向力

横向力取决于结构刚度、摩擦、间隙和接触状态。没有独立刚度标定时，直接相乘会产生无依据的数值。

10.4 测点间距越大越好吗

间距增大有利于提高角度分辨能力，但会受到结构弯曲、空间和传感器量程限制，应以刚体假设成立为前提。

10.5 安装偏心会影响纳米定位吗

可能影响。偏心引起的横向载荷和摩擦变化会叠加到迟滞、蠕变和温漂上，最终表现为重复性或稳定时间变化。

10.6 如何确认是夹具问题而不是执行器问题

可做空夹具、换向、重复装夹和标准件对比。若误差随夹具或装夹方向变化，应先排查机械链路。

10.7 文章中的 0.010 毫弧度是产品指标吗

不是。它由一组教学输入计算得到，不含具体型号、仪器、样本数和校准记录。

10.8 什么时候需要多点位姿拟合

当结构存在明显弯曲、大角度转动、多个自由度或测点不同步时，应使用三点以上测量和刚体位姿拟合。

11. 结论

压电执行器安装误差测试的关键，是先定义作用线和机械基准，再用双点位移测量倾斜，用独立力测量或已标定模型评估横向力。倾斜角、轴向位移、横向力和定位精度不能混为一个指标。

本文代码和数据均为可复现教学示例。正式发布某型号的安装误差或纳米定位性能时，必须补充传感器型号、校准状态、负载、采样率、样本数量、环境、统计方法和原始记录。

参考资料

1. ISO 230-2:2014，Test code for machine tools — Part 2: Determination of accuracy and repeatability of positioning of numerically controlled axes。ISO标准页面
2. ISO 5725-1:1994，Accuracy (trueness and precision) of measurement methods and results。ISO标准页面
3. JCGM 100:2008，Evaluation of measurement data — Guide to the expression of uncertainty in measurement。BIPM官方页面
4. ANSI/IEEE Std 176-1987，IEEE Standard on Piezoelectricity。IEEE官方来源
5. ISO/IEC 17025:2017，General requirements for the competence of testing and calibration laboratories。ISO标准页面
6. Python Software Foundation，csv、statistics、math模块文档，Python 3.12。Python官方文档