Application of Active Magnetic Bearing in Waste Heat Generator
-
摘要:
针对传统余热发电设备存在发电效率低、成本高、体积大等问题，设计了一款装备5自由度（5-DOF）主动磁悬浮轴承的余热发电机，用于余热发电有机朗肯循环系统. 首先，基于电机转子轴径限制和最大承载力要求，确定径向磁悬浮轴承和轴向磁悬浮轴承的结构形式，采用一维磁路模型和二维有限元分析，计算并校核磁悬浮轴承的尺寸和性能参数；其次，为保证磁悬浮轴承的稳定裕度，采用三电平PWM功率放大器降低输出电流纹波，利用不完全微分PID控制器和不平衡补偿算法实现5自由度稳定悬浮；最后，将该磁悬浮余热发电机应用到实际客户现场，从稳定裕度、承载力和轴振峰峰值3个维度验证磁悬浮余热发电机的可靠性. 现场试验结果表明：磁悬浮余热发电机系统运行稳定可靠，实现了满功率发电和长期运行考核；在全工况运行范围，磁悬浮轴承系统的灵敏度函数小于12 dB，余热发电机转子振动峰峰值小于53 μm，满足ISO14839-3规定的长期稳定运行要求；轴向承载力最大达到3 600 N，满足实际工况需求.
-
关键词:
- 有机朗肯余热发电机 /
- 主动磁悬浮轴承 /
- 三电平PWM功率放大器 /
- 转子不平衡补偿
Abstract:Aiming at the problems of low efficiency, high cost and large volume of the traditional waste heat power generation equipment, a waste heat generator equipped with the 5-DOF (5-degree of freedom) active magnetic bearing for Organic Rankine waste heat generation system was designed. Firstly, based on the rotor diameter limit and maximum bearing capacity requirements, the structural forms of radial and axial magnetic bearings were determined. Furthermore, the dimensions and performance parameters of magnetic bearings were calculated and checked by one-dimensional magnetic circuit model and two-dimensional finite element analysis. Secondly, in order to ensure the stability margin of the magnetic bearing, a three-level PWM (pulse-width modulation) power amplifier was used to reduce the output current ripple, and an incomplete differential PID controller incorporated with unbalance compensation were used to realize the 5-DOF stable suspension. Finally, the waste heat generator equipped with magnetic bearing was applied to the actual customer site. The reliability of the waste heat generator was verified from three dimensions of stability margin, bearing capacity and shaft peak-to-peak value. The field test results show that the waste heat generator system has operated stably and reliably in all operation conditions, and can achieve full-power generation and long-term operation. The sensitivity function of the magnetic bearing system is less than 12 dB, and the peak-to-peak vibration value of the rotor of the waste heat generator is less than 53 μm, which meet the long-term stable operation requirements of ISO14839. The maximum axial bearing capacity reaches 3 600 N, which meets the requirements of actual working conditions.
-
随着三维扫描仪等硬件设备的成本降低和市场增长，其应用越来越广泛. 在三维扫描应用场景中，由于每个摄像头或扫描仪在其自身的空间中生成点云，而并非在对象空间中，所以2个点云模型之间没有对应关系，使点云的数据融合[1]、交叉纹理映射[2]、三维人体扫描[3]、即时定位与高精度地图构建[4]等后处理过程面临巨大挑战. 因此，构建点云模型之间的对应关系成为一个亟待解决的问题，特别是构建非刚性模型间的对应关系更具挑战性. 非刚性模型发生的局部形变会导致模型间特征差异过大. 此外，离散的点云模型缺少几何先验知识，使得无监督学习在计算对应关系时很难实现，在实际应用中，往往需要大量的带标注数据. 因此，如何采用无监督学习准确构建非刚性三维点云模型间密集的对应关系是一个极具挑战性的问题.
近年来，主流的三维模型对应关系计算方法主要有2类：基于模型的方法[5-8]和数据驱动的方法[9-10]. 基于模型的方法使用人工制作的特征描述符作为输入，优化预定义的函数映射；而基于数据驱动的方法使用大规模数据集训练神经网络来预测对应关系. 近期的研究表明，数据驱动方法计算精度明显高于基于模型的方法[9]. 但数据驱动方法需要大量的标注数据和使用模型中顶点之间的连接信息（例如三角网格模型）. 研究者通过计算模型在拉普拉斯-贝尔特拉米算子特征基[9-11]上的投影，来实现三角网格模型间对应关系映射，该方法已被广泛使用，且证明是高效的.
研究者致力于非刚性模型的密集对应关系研究，其目的是计算发生非刚性形变的三维模型间逐点对应关系. 其中，针对于三角网格模型对应关系计算的一个重要方法是函数映射法. 该方法使用简化后的基函数将模型的特征描述符表示为函数映射矩阵，将模型间对应关系计算问题转化为求解相应的函数映射矩阵问题[8,12]. 文献[13]结合局部函数映射和局部流形谐波算子计算三维模型特征描述符，并构建模型间对应关系. 文献[14]提出同时计算完整模型与多个发生非刚性形变残缺模型的对应关系计算方法，并且对带有噪声的模型得到更准确的匹配结果. 该类方法的优点是使用函数映射进行结构化对应预测. 然而，频谱基的计算需要较高的算力和预处理时间，并且需要点之间的连通性信息. 文献[15]提出一种点云频谱匹配方法，利用神经网络学习模型基函数，从而避免对拉普拉斯算子进行特征分解的计算. 文献[16]提出基于上采样细化和全频谱特征值对齐的方法，提高了非刚性残缺模型对应关系结果的平滑度和准确性，对带有拓扑噪声的三维模型也具有较好的鲁棒性，但是该方法需要大量特征值来定位正则化项. 文献[17]提出三维模型之间的基矩阵校准方法，通过尖端点遍历匹配即可获取较优的对应关系结果. 这些方法在谱域计算模型拉普拉斯-贝尔特拉米算子特征基的预处理过程中非常耗时，占用过多计算资源. 此外，拉普拉斯-贝尔特拉米算子的计算需要模型顶点之间连通性信息，但当数据来源于三维扫描设备时，点云模型采样点之间是离散的.
为解决上述问题，研究者提出了多种三维点云模型对应关系计算方法[18-20]，这些方法通常采用编码器-解码器结构. 文献[18]使用给定的模板模型以重建不同的输入点云，解码器返回每个点的坐标，并根据重建模型与模板的最近邻点确定对应关系；但是该方法严重依赖人工设计的参数化模板质量. 文献[20]使用交叉重建的方式摒弃了参数化模板，但是非刚性模型中对应点特征差别较大，其使用最近邻搜索寻找匹配特征的方法泛化能力较差. 文献[21]提出利用特征序列匹配的思想来预测2个点云模型的刚性变换，但该方法无法计算非刚性模型的对应关系.
综上所述，目前针对三维点云模型对应关系计算的研究中，由于离散的点云模型缺少几何先验知识，难以实现无监督学习；此外，当模型发生非刚性形变时，现有的特征提取方法难以获得唯一且具有丰富的局部与全局信息的三维点云模型本征描述符，导致对应关系计算准确率不高. 为此，本文提出一种基于特征序列注意力机制和结构化构造的无监督三维点云模型对应关系计算方法.
1. 序列到序列的无监督点云模型对应关系计算
设输入的三维点云源模型和目标模型分别为${\mathcal X} $,$ {{\mathcal Y}} \in {\mathbb{R}}^{n \times 3} $，其中，n为点的个数. 本文的目标是在源模型和目标模型之间寻找一个映射，即在目标模型$ {\mathcal{Y}} $中寻找与点$ {x_i} \in {\mathcal X} $的对应点$ {y_j} \in {\mathcal{Y}} $，$ 1 \leqslant i,j \leqslant n $. 为此，构建如图1所示的神经网络架构，主要由特征提取模块、Transformer模块、对应关系预测模块和点云重构模块构成. 图中：$ {{\boldsymbol{F}}_{\mathcal X}} $、$ {{\boldsymbol{F}}_{\mathcal{Y}}} $为点云特征，$ {\hat{\boldsymbol{ F}}_{\mathcal X} } $、$ {\hat{\boldsymbol{ F}}_{\mathcal Y}} $为施加了注意力的点云特征，$ {{\boldsymbol{S}}_{ \mathcal{XY}}} $、$ {{\boldsymbol{S}}_{\mathcal{YX}}} $为软映射矩阵，$ {\hat {\mathcal X}_{\mathcal Y}} $、$ {\hat {\mathcal{Y}}_ {\mathcal X}} $为重建点云. 首先，使用特征提取模块提取三维点云模型特征，通过Transformer模块对源模型和目标模型的特征序列施加注意力机制；然后，通过对应关系预测模块得到软映射矩阵；最后，通过软映射矩阵确定目标点及其邻域，并对目标点邻域加权融合，以进行点云重构.
1.1 特征提取模块
本文的核心任务是寻找2个模型间相匹配的特征，所以获取点云模型准确的局部特征尤为重要. 为获得较好的点云模型局部特征，使用具有多个边缘卷积层[22]的动态图卷积网络，将源模型$ \mathcal X $和目标模型$ {\mathcal{Y}} $嵌入到高维特征空间中. 具体来说，对于点云中第i个点，首先计算该点特征到其最近邻特征集合$ \varOmega _{i,l} $中各特征的欧氏距离，其中，l为边缘卷积层的序号. 然后，把该点特征$ {\boldsymbol{f}}_{i,l} $和该点特征到最近邻特征的距离$ \{ {\boldsymbol{f}}_{j,l} - {\boldsymbol{f}}_{i,l}|{\boldsymbol{f}}_{j,l} \in \varOmega _{i,l}\} $相拼接（$ {\boldsymbol{f}}_{i,l} $,$ {\boldsymbol{f}}_{j,l} \in {\mathbb{R}^{1 \times d}，} $d为特征维度，$ i \ne j $），通过多层感知机和最大池化获得新的特征表示：
fi,l+1=Γ(Ml(fi,l,fj,l−fi,l)),fj,l∈Ωi,l, (1) 式中：$ \varGamma $(•)为最大池化操作，Ml (•)为第l个边缘卷积中的多层感知机.
边缘卷积中使用的最近邻域特征聚合操作可以使得$ {\boldsymbol {f}}_{i,l + 1} $更好地捕捉点云的局部几何信息.
1.2 Transformer模块
由特征提取模块得到的逐点特征虽然在一定程度上可以较好地表征其局部特征，但是在非刚性模型中，局部特征差异较大，难以获得较好的匹配结果. 本文使用Transformer模块中的注意力机制[23]完成特征序列到特征序列的学习，从而优化非刚性模型间的特征匹配. 经特征提取模块获得的模型$ {\mathcal X} $、$ {\mathcal {Y}} $的特征序列$ {{\boldsymbol{F}}_ {\mathcal X}} $,$ {{\boldsymbol{F}}_{\mathcal Y}} \in {\mathbb{R}^{n \times d}} $. 对应关系可以通过神经网络在目标模型$ {\mathcal {Y}} $中预测一个与源模型$ {\mathcal X} $中一个点对应的点来构建，此时可以利用注意力机制使$ {{\boldsymbol{F}}_{\mathcal Y}} $对$ {{\boldsymbol{F}}_ {\mathcal X}} $施加注意力，得到兼顾模型$ {\mathcal {Y}} $的模型$ {\mathcal X} $的特征信息；同理，也可得到兼顾模型$ {\mathcal X} $的模型$ {\mathcal {Y}} $的特征信息.
Transformer模块的任务是学习一个函数$ g:{\mathbb{R}^{n \times P}} \times {\mathbb{R}^{n \times P}} \to {\mathbb{R}^{n \times P}} $，P为特征嵌入的维度，并由函数g(·)得到点云新的特征嵌入，具体过程可表示为
{ˆFX=FX+g(FX,FY),ˆFY=FY+g(FY,FX). (2) 如图2所示，Transformer模块为编码器-解码器结构. 前后传播模块由1个线性层和1个ReLu激活函数层组成. 编码器的功能为：输入特征序列$ {{\boldsymbol{F}}_ {\mathcal X}} $，经过自注意力机制和共享多层感知机将其编码并嵌入到特征空间中. 解码器分为两部分：第一部分输入另一个特征序列$ {{\boldsymbol{F}}_{\mathcal Y}} $，并以与编码器相同的方式对其进行编码；第二部分使用交叉注意力机制关联2个嵌入的特征序列. 具体而言，注意力机制将输入的线性编码表示分别建模为查询向量$ {\boldsymbol{Q}} \in {\mathbb{R}^{n \times P}} $、关键向量$ {\boldsymbol{K}} \in {\mathbb{R}^{n \times P}} $和值向量$ {\boldsymbol{V}} \in {\mathbb{R}^{n \times P}} $，从而可以计算出注意力得分矩阵W=Softmax(QKTd−0.5)，最终得到输出特征$ \hat{\boldsymbol{ F}} = {\boldsymbol{WV}} $. 其中，Softmax(·)可以将输出值转换为0～1且和为1的概率分布. 图2中，编码器和解码器的输入分别为$ {{\boldsymbol{F}}_ {\mathcal X}} $和$ {{\boldsymbol{F}}_{\mathcal Y}} $，可得到点云模型$ {\mathcal X} $的新特征$ {\hat{\boldsymbol{ F}}_ {\mathcal X}} $；同理，若将$ {{\boldsymbol{F}}_ {\mathcal X}} $和$ {{\boldsymbol{F}}_{\mathcal Y}} $调换输入位置（即编码器的输入为$ {{\boldsymbol{F}}_{\mathcal Y}} $，解码器的输入为$ {{\boldsymbol{F}}_ {\mathcal X}} $），可得到点云模型$ {\mathcal {Y}} $的新特征$ {\hat{\boldsymbol{ F}}_{\mathcal Y}} $. 因此，输出的特征序列$ {\hat{\boldsymbol{ F}}_ {\mathcal X}} $和$ {\hat{\boldsymbol{ F}}_{\mathcal Y}} $具有来自于2个特征序列（$ {{\boldsymbol{F}}_ {\mathcal X}} $和$ {{\boldsymbol{F}}_{\mathcal Y}} $）的上下文信息.
1.3 对应关系预测模块
在大多数特征匹配的任务中，直接使用特征最近邻搜索即可得到对应关系，但是在网络训练之初，预测结果远未达到较好效果的时候，这种方法容易使神经网络模型陷入局部最优. 为解决这个问题，本文使用概率来定位对应点，即生成从一组特征序列到另一组特征序列的概率映射矩阵，称为软映射矩阵. 软映射定义为源模型的特征嵌入和目标模型的特征嵌入之间的相似度，所以，软映射矩阵明确了源模型上任一点与目标模型上任一点之间的相似度. 软映射矩阵由式（3）求出.
{SXY=Softmax(ˆFXˆFTY),SYX = Softmax(ˆFYˆFTX), (3) 式中：$ {{\boldsymbol{S}}_{ {\mathcal X}{\mathcal Y}}} $,$ {{\boldsymbol{S}}_{{\mathcal Y} {\mathcal X}}} \in {\mathbb{R}^{n \times n}} $，矩阵中第i行第j列的元素${S_{ij}} \in \left[ {0,1} \right]$，表示$ {\mathcal X} $中点i和$ {\mathcal {Y}} $中点j的对应概率.
本文最终的对应关系由软映射矩阵中最大对应概率确定. 但由于使用的是无连接信息的点云模型，所以采用重建的方式构造损失函数来监督训练，从而优化软映射矩阵.
1.4 点云重构模块
在使用谱方法计算三角网格模型的对应关系时，可以根据三角网格模型自身的几何拓扑结构应用无监督学习方法，但是本文算法以点云模型作为输入，采样点之间不具有连接信息，很难应用自身几何性质完成无监督训练. 为了在没有真实对应关系的情况下学习点云模型的特征，设计一个点云重构模块，通过重构输入点云，并利用重构模型与输入模型的相似度来间接监督训练.
早期针对点云模型的对应关系计算方法是通过模板和解码器完成三维重建，利用最近邻搜索进行无监督训练，但是需要大量训练数据，并且难以在跨数据集评估中推广. 本文使用特征序列相似性和输入坐标来重构点云并确定对应关系，以取代使用坐标回归的解码器模块.
如图3所示，首先，k近邻特征搜索模块接收2组三维点云模型和软映射矩阵作为输入，寻找点云模型$ {\mathcal {Y}} $上与模型$ {\mathcal X} $上第i个点对应的概率最大k个点，可以表示为$ {\{ {y_i}\} _{i \in [1,k]}} $. 其次，根据软映射矩阵计算出$ {\mathcal X} $上第i个点与点集{yi}中各点的对应概率Sij，并使用Softmax进行归一化处理再求和，得到权重集合{wij}，如式（4）所示.
wij=eSij/∑h∈N(xi)eSih, (4) 式中：$ N({x_i}) $为软映射矩阵中点云模型$ {\mathcal {Y}} $上与模型$ {\mathcal X} $上第i个点对应的概率最大k个点的序号.
最后，根据得到的权重集合{wij}和点集{yi}计算重构后的点$ {\hat y_{{x_i}}} $如式（5）所示.
ˆyxi=∑h∈N(xi)wihyh. (5) 同理，根据对应概率将模型$ {\mathcal X} $上与模型$ {\mathcal {Y}} $上第i个点对应的k个点$ {\{ {x_i}\} _{i \in [1,k]}} $和$ {\mathcal {Y}} $上第i个点与{xi}的对应概率集合{Sij}使用式（4）得到权重集合{wij}，根据式（5）得到重构后的点$ {\hat x_{{y_i}}} $. 此时，式（4）、（5）中的N(xi)需修改为软映射矩阵中模型$ {\mathcal X} $上与模型$ {\mathcal {Y}} $上第i个点对应的概率最大k个点的序号.
点云模型$ {\mathcal X} $和$ {\mathcal {Y}} $上所有点经过重构模块后得到重构点云$ {\hat {{\mathcal Y}}_ {\mathcal X}} $,$ {\hat { {\mathcal X}}_{\mathcal Y}} \in {\mathbb{R}^{M \times 3}} $，且$ {\hat y_{{x_i}}} \in {\hat {{\mathcal Y}}_ {\mathcal X}} $，$ {\hat x_{{y_i}}} \in {\hat { {\mathcal X}}_{\mathcal Y}} $.
1.5 损失函数
在对应关系计算任务中，有监督训练依赖于人工标注的点对点真实对应关系来监督特征匹配结果. 而无监督训练则在未知真实点对点对应关系的情况下进行，通过利用模型数据自身的几何性质进行监督，从而节约人工标注成本.
定义损失函数为
L=λ1LCD+λ2LSC+λ3LNC, (6) 式中：$ {\lambda _1} $、$ {\lambda _2} $和$ {\lambda _3} $为系数，依次取为1、10和1，LCD、LSC和LNC分别为倒角距离损失函数、促进特征平滑的正则化项和相邻点对应的正则化项.
输入点云模型经过重构点云步骤，$ {\mathcal X} $和$ {\hat { {\mathcal X}}_{\mathcal Y}} $、$ {\mathcal {Y}} $和$ {\hat {{\mathcal Y}}_ {\mathcal X}} $2组点云应分别保持近似一致，所以使用$ {\mathcal X} $与$ {\hat { {\mathcal X}}_{\mathcal Y}} $、$ {\mathcal {Y}} $与$ {\hat {{\mathcal Y}}_ {\mathcal X}} $之间的倒角距离作为损失函数，如式（7）所示.
LCD=C(X,ˆXY)+C(Y,ˆYX), (7) C(X,Y)=1|X|∑x∈Xminy∈Y‖x−y‖22+1|Y|∑y∈Yminx∈X‖y−x‖22, (8) 式中：C(·)为倒角距离，||·||2为两点间欧式距离.
为使特征提取模块和Transformer模块得到的特征更平滑，即使得相邻点的特征更加相似，本文利用点云重构模块将$ {\mathcal X} $和$ {\mathcal {Y}} $分别根据其自身进行重构，得到重构模型$ {\hat { {\mathcal X}}_ {\mathcal X}} $和$ {\hat {{\mathcal Y}}_{\mathcal Y}} $，并由此计算正则化项LSC（式（9））. 如果重建结果与自身相似度较大，则模型上的点到其邻域点的对应概率较大. 因为采用了特征匹配的思想，对应概率大意味着特征相似度高.
LSC=C(Y,ˆYY)+C(X,ˆXX). (9) 此处的重构模型是通过点云重构模块，根据模型自身及其自身到自身的软映射矩阵计算得到的. 其自身到自身的软映射矩阵为
{SXX=Softmax(ˆFXˆFTX),SYY=Softmax(ˆFYˆFTY). (10) 将点云$ {\mathcal X} $和对应的软映射矩阵$ {{\boldsymbol{S}}_{ {\mathcal X} {\mathcal X}}} $输入重构模块后，由式（11）即可求得模型$ {\mathcal X} $的重构模型$ {\hat { {\mathcal X}}_ {\mathcal X}} $.
ˆxi=∑h∈N(xi)wc,ihxh, (11) 式中：$ {\hat x_i} \in {\hat { {\mathcal X}}_ {\mathcal X}} $，为模型$ {\mathcal X} $中第i个点的重构点；wc,ih为第i个点到N(xi)中各点对应分配的权重，由式（4）计算得到.
同理，$ {\hat {{\mathcal Y}}_{\mathcal Y}} $也可根据点云$ {\mathcal {Y}} $和软映射矩阵$ {{\boldsymbol{S}}_{{\mathcal Y}{\mathcal Y}}} \in {\mathbb{R}^{M \times M}} $求出.
为了使源点的邻点与源点对应的目标点的邻点相对应，本文设置了正则化项LNC，如式（12）所示. 该正则化项通过源模型中一点到其邻域点的距离对预测对应点到其邻域点的距离进行加权，从而促使源模型中相邻点在目标模型中也保持相邻.
LNC=1nk∑i∑h∈NX(xi)vih‖ˆyxi−ˆyxh‖22, (12) 式中：$ {N_ {\mathcal X}}({x_i}) $为点云模型$ {\mathcal X} $中与xi相邻点的集合； $ {v_{ih}} = {{\mathrm{e}}^{{{ - \left\| {{x_i} - {x_h}} \right\|_2^2} / \alpha }}} $，为根据邻域内各点到中心点的距离所计算的权重，$\alpha $=8.
2. 实验结果与分析
2.1 实验数据集与评估标准
2.1.1 数据集
在人体和非人体数据集上对本文方法进行测试. 对于人体模型，训练时采用从参数化人体模型SMPL[24]生成的Surreal数据集，共包含
230000 个训练模型. 将数据集划分为115000 对，按对输入网络模型进行训练；测试时，使用FAUST数据集[3]和SHREC’19数据集[25]. FAUST数据集包含10个真人扫描模型，每个人体模型包含10个不同姿态. SHREC’19数据集包含44个真人扫描模型，共有430个带标签的测试对，模型间具有较大的非刚性形变. 因此，对应关系可分为同一个人体不同姿态的模型间对应关系（类内）和不同人体模型之间的对应关系（类间）.对于非人体模型，采用SMAL数据集[26]进行训练和评估. SMAL包含猫、狗、牛、马和河马的动物模型. 本文使用其提供的模板为每种动物类型创建
2000 个不同姿态的模型，获得包含共10000 个模型的数据集，每个类别又分别生成40个测试模型. 在训练时，从同一动物类别中随机选择模型对作为神经网络的输入；在测试时，将测试数据中同类的两两模型组成一对，形成一个100对的测试集.实验从每个模型中随机抽取
1024 个点，随机排列点的顺序，创建用于训练和测试的点云模型.2.1.2 定量评估标准
模型对应的通用评估指标是平均测地误差[27]，但是计算测地误差时需要已知模型中点的连接信息. 由于本文使用的是点云数据，采样点之间缺少连接信息，因此，需要使用欧氏距离来代替测地距离，构建对应误差进行评估.
给定源模型$ {\mathcal X} $和目标模型$ {\mathcal {Y}} $，对应误差为
e=1n∑x∈X,ygt∈Y‖φ(x)−ygt‖2, (13) 式中：$ {y_{{\mathrm{gt}}}} \in {{\mathcal Y}} $，为点$ x $的真实对应点；$ \varphi (x) $为点$ x $根据计算对应关系映射到模型$ {\mathcal {Y}} $上的一个点，即$ \varphi (x) \in {{\mathcal Y}} $.
为进一步评估对应关系的计算准确率，使用式（14）计算在不同的对应误差容忍度量下对应关系精度.
a=1n∑x∈X,ygt∈Yξ(‖φ(x)−ygt‖2⩽ε), (14) 式中：$ \xi $(·)为指示器函数，$ \varepsilon $为对应误差容忍度量.
2.2 实验结果与分析
将本文方法与最新的针对点云模型的对应关系计算方法进行定性和定量对比，包括3D-CODED[18]、Elementary Structures[19]和CorrNet3D[20]，其中，3D-CODED和Elementary Structures为不使用点之间连接信息的有监督方法.
2.2.1 定性实验
为测试本文方法的跨数据集泛化能力，在FAUST数据集上测试，可视化对比结果如表1所示. 表中，源模型点的颜色是根据点的空间坐标值设定的，此外，根据真实对应关系和各算法计算得到的对应关系，将源模型点的颜色映射到Ground-truth和其余4种方法的模型上. 其中，Ground-truth为对应关系的真值. 可以看出，由于3D-CODED的重建模型在躯干伸展的方向和源模型之间存在角度误差，所以在手部和腿部产生错误的对应关系；Elementary Structures在手部、脚部和头部存在较多错误对应，这是因为在特征提取阶段没有强制施加特征平滑措施；CorrNet3D依赖于解码器网络来回归重建模型坐标，所以模型的手部、腹部和腿部均存在较多错误的对应关系. 与上述方法相比，本文方法既不需要使用模板，也不需要使用解码器重建，而是借助于Transformer完成点云的序列到序列对应，并实现了结构化重建，获得更好的对应关系结果.
表 1 FAUST数据集上各方法的实验结果对比Table 1. Comparison of experimental results of methods onFAUST dataset示例 源模型 Groud-truth 3D-CODED Elementary Structures CorrNet3D 本文算法 1 
2 
为更好地测试本文方法的跨数据集泛化能力和在非刚性模型上对应关系的构建能力，在SHREC’19数据集上进行了实验，可视化对比结果如表2所示. 3D-CODED和Elementary Structures都依赖于模板分析源点云模型和目标点云模型的变形来构建他们之间的对应关系，但是当模型发生较大的非刚性形变时，由解码器得到的重建模型效果较差，使得计算的对应关系准确率下降，如图中人体模型的腿部和手部的错误对应；CorrNet3D依赖于解码器网络回归重建模型坐标，没有对特征进行平滑处理，所以泛化能力弱，且在计算软映射矩阵时容易产生错误对应. 以上3种方法均需要使用模型的全局特征进行重建，这导致其跨数据集的泛化能力较差，而本文方法既不需要使用模型的全局特征，也不需要使用解码器重建，而是借助于Transformer模块完成点云的序列到序列对应和结构化重建，提高了非刚性模型间对应关系计算的准确度和跨数据集的泛化能力.
表 2 SHREC’19数据集上各方法的实验结果对比Table 2. Comparison of experimental results of methods on SHREC’19 dataset示例 源模型 Groud-truth 3D-CODED Elementary Structures CorrNet3D 本文算法 3 
4 
表1和表2的实验都是在人体模型数据集上进行的，为进一步评估本文方法的灵活性和准确性，在SMAL动物数据集上进行训练和测试，可视化对比结果见表3. 3D-CODED和Elementary Structures算法结果存在较多错误的对应点，因为他们都依赖于模板对源模型和目标模型进行变形来预测对应关系. 相比之下，CorrNet3D在这方面有一定改进，但由于依赖训练解码器网络来重建模型，在多种尺度混合数据集中表现不佳. 而本文方法采用对应的最近邻点加权融合的重建方法，增强了算法的健壮性和对多尺度数据集的适应能力.
表 3 SMAL数据集上各方法的实验结果对比Table 3. Comparison of experimental results of methods on SMAL dataset示例 源模型 Groud-truth 3D-CODED Elementary Structures CorrNet3D 本文算法 5 
6 
7 
2.2.2 定量实验
为进一步评估本文方法的有效性，对不同方法的对应误差变化曲线进行计算，结果如图4所示. 由图可知：1） 在SHREC’19数据集上，本文方法有14.1%的对应关系没有对应误差，而3D-CODED、Elementary Structures和CorrNet3D在没有对应误差时的对应关系准确率明显低于本文方法；且在对应误差增大时，本文方法的对应关系准确率始终高于这3种方法. 2） 在FAUST数据集上，本文方法有24.6%的对应关系没有对应误差，高于其他3种方法；且在对应误差为0.04时，CorrNet3D和本文方法的对应关系准确率接近100%，结果优于3D-CODED和Elementary Structures. 3） 在SMAL数据集上，本文方法有20.2%的对应关系没有对应误差，高于其他3种方法的结果；且在对应误差为0.03时，本文方法的对应关系准确率接近100%，3D-CODED、Elementary Structures和CorrNet3D均未能达到100%. 这主要归功于基于注意力的特征序列对应和提取到的平滑特征，使本文方法取得了更好的结果.
实验进一步计算测试集中的平均对应误差，对比结果如表4所示. 可以看出，本文方法在数据集SGREC’19、FAUST 和SMAL上均取得了更低的平均对应误差，尤其是在非刚性变形程度较大的SHREC’19数据集中平均对应误差下降幅度更大.
表 4 不同方法平均对应误差Table 4. Average correspondence errors of different methods数据集 3D-
CODEDElementary
StructuresCorrNet
3D本文
算法FAUST 7.2 6.6 5.3 5.1 SHREC’19 8.1 7.6 6.9 5.8 SMAL 6.8 6.4 5.7 5.4 2.2.3 消融实验
通过设置有Transformer模块和无Transformer模块两组实验，验证本文方法中Transformer模块的作用，在数据集FAUST、SHREC’19和SMAL上的对比结果如表5所示. 可以看出，有Transformer模块的网络平均对应误差低于无Transformer模块的网络，进一步验证了本文提出的特征序列匹配思想可以显著提升对应关系准确率.
表 5 Transformer模块消融实验误差Table 5. Errors of ablation experiment on Transformer module数据集 无 Transformer 模块 有 Transformer 模块 FAUST 5.6 5.1 SHREC’19 6.1 5.8 SMAL 7.6 5.4 此外，设置2两组实验验证对应关系预测模块的有效性，结果如表6所示，其中，无对应关系预测模块实验中使用最近邻搜索策略进行特征匹配. 可以看出，有对应关系预测模块的网络平均对应误差均低于无对应关系预测模块的网络. 所以，基于概率的方法可以显著提升对应关系准确率，且在当模型发生较大尺度非刚性变形时表现更好.
表 6 对应关系预测模块消融实验误差Table 6. Errors of ablation experiment on correspondence prediction module数据集 无对应关系
预测模块有对应关系
预测模块FAUST 6.1 5.1 SHREC’19 7.3 5.8 SMAL 8.0 5.4 3. 结束语
本文结合注意力机制，使用特征序列注意力匹配的方法生成软映射矩阵，有效提高了无监督深度学习方法计算三维点云模型间稠密对应关系的准确率，在多个数据集上表现出较好的泛化能力. 然而，本文方法也存在需要改进的地方，如使用Transformer模块导致了网络参数增加；与针对三角网格模型的对应关系计算结果相比，仍有一定的提升空间；离散点云模型目前必须借助于重建才能实现无监督训练等，这些都是今后需要继续研究的方向.
致谢： 兰州交通大学天佑创新团队项目（TY202002）的资助.
-
图 1 磁悬浮轴承余热发电机的结构示意
Figure 1. Schematic diagram of magnetic bearing waste heat generator
图 2 径向磁轴承E型电磁铁结构及尺寸定义
Figure 2. Structure and dimension definition of E-type electromagnet for radial magnetic bearing
图 4 上线圈电流最大、下线圈电流为0的磁场分布
Figure 4. Magnetic field distribution when the top coil current is maximum and the bottom coil is zero
图 5 电磁力与控制电流的关系曲线
Figure 5. Relationship between the electromagnetic force and the control current
图 6 电磁力与转子位移的关系曲线
Figure 6. Relationship between the electromagnetic force and the rotor displacement
图 7 轴向磁轴承电磁铁结构及尺寸定义
Figure 7. Structure and dimension definition of the electromagnet for the axial magnetic bearing
图 9 上端线圈电流为0、下线圈电流最大的磁场分布
Figure 9. Magnetic field distribution when the top coil current is zero and the bottom coil current is maximum
图 10 电磁力随控制电流变化曲线
Figure 10. Relationship between the electromagnetic force and the control current
图 11 电磁力与转子位移的关系曲线
Figure 11. Relationship between the electromagnetic force and the rotor displacement
图 12 磁悬浮轴承余热发电机现场安装
Figure 12. Site installation drawing of waste heat generator with magnetic bearing
图 16 额定转速时空载与负载条件下灵敏度函数对比
Figure 16. Sensitivity function comparison between no-load and full load at rated speed
图 17 实际测量的轴向磁轴承电流增益曲线
Figure 17. Measured curves of the current-force coefficient of axial magnetic bearing
图 18 整个转速范围内的转子轴振峰峰值曲线
Figure 18. Peak-to-peak curves of rotor shaft vibration in the whole speed range
表 1 径向磁轴承设计输入条件
Table 1. Input conditions of radial magnetic bearing
变量 符号 数值 备注 额定气隙/mm x0 0.3 最大磁感应强度/T Bmax 1.20 最大承载力/N Fmax 2 250 径向转子外径/mm Dr2 170.0 接口尺寸 定子内径/mm Ds1 170.6 
下载: 导出CSV
表 2 轴向磁轴承设计输入条件
Table 2. Input conditions of axial magnetic bearing
变量 符号 数值 额定气隙/mm x0 0.4 最大磁感应强度/T Bmax 1.25 最大承载力/N Fmax 2 600 电机内部安装外径/mm dlim 394
下载: 导出CSV
表 3 额定转速条件下，空载与负载运行的敏感性函数峰值
Table 3. Sensitivity function peaks of no-load and load operations at rated speed
类别 工况 峰值/dB 转速/krpm 非推力侧径向磁轴承 X 方向 无工质 9.34 4.44 有工质 8.47 4.50 非推力侧径向磁轴承 Y 方向 无工质 11.65 4.62 有工质 10.07 4.92 推力侧径向磁轴承 X 方向 无工质 8.97 3.90 有工质 7.30 4.86 推力侧径向磁轴承 Y 方向 无工质 9.54 3.66 有工质 7.58 3.24 轴向磁轴承 Z 方向 无工质 10.19 2.58 有工质 10.00 2.76
下载: 导出CSV
-
[1] 孟祥睿, 马新灵. 有机朗肯循环低品位热能发电技术[M]. 郑州: 河南科技出版社, 2015. [2] 纪历. 低温余热磁悬浮发电机转子轴系振动机理研究[J]. 电机与控制学报,2019,23(11): 67-75.JI Li. Vibration mechanism analysis of magnetic levitation rotor system for low temperature waste heat power generation[J]. Electric Machines and Control, 2019, 23(11): 67-75. [3] 孙立佳,任小坤,高元景,等. 高速透平发电机的特点及相关技术研究[J]. 低温与超导,2015,43(8): 23-26.SUN Lijia, REN Xiaokun, GAO Yuanjing, ea al. Study on feature and related technology of high speed turbine generator[J]. Cryogenics & Superconductivity, 2015, 43(8): 23-26. [4] 梁欣. 余热发电高速永磁电机设计与分析[D]. 武汉: 华中科技大学, 2019. [5] 戴其城. 涡轮直驱余热发电机组中锥形永磁同步电机的研究[D]. 长沙: 湖南大学, 2018. [6] 胡业发, 周祖德, 江征风. 磁力轴承的基础理论与应用[M]. 北京: 机械工业出版社, 2006. [7] 王大彪,段 捷,胡哺松,等. 有机朗肯循环发电技术发展现状[J]. 节能技术,2015,33(191): 235-242.WANG Dabiao, DUAN Jie, HU Bu song, et al. Status of organic Rankine cycle power generation technology[J]. Energy Conservation Technology, 2015, 33(191): 235-242. [8] 占智军,祝长生. E型径向电磁轴承的参数设计及特性分析[J]. 机电工程,2013,30(3): 267-272.ZHAN Zhijun, ZHU Changsheng. Parameters design and characteristics analysis of E-core radial magnetic bearings[J]. Journal of Mechanical & Electrical Engineering, 2013, 30(3): 267-272. [9] 杨贝,张宁,张海波,等. 径向磁悬浮轴承的结构设计和数值模拟[J]. 低温与超导,2015,43(4): 31-35,42.YANG Bei, ZHANG Ning, ZHANG Haibo, et al. Design and simulation of radial active magnetic bearings[J]. Cryogenics & Superconductivity, 2015, 43(4): 31-35,42. [10] 代燕杰. 磁悬浮推力轴承的电磁性能分析及结构设计[D]. 济南: 山东大学, 2008. [11] 吴宝贵. 轴向磁悬浮轴承的研究[D]. 南京: 南京航空航天大学, 2008. [12] 邹望蠡. 轴向磁轴承的拓扑结构及参数优化设计[D]. 南京: 南京航空航天大学, 2014. [13] 刘程子. 适用于高速电机的混合型磁悬浮轴承设计及控制策略的研究[D]. 南京: 南京航空航天大学, 2014. [14] 魏坚. 磁悬浮飞轮电池支承控制系统的硬件设计与算法研究[D]. 武汉: 武汉理工大学, 2010. [15] 刘建明. 电磁轴承结构设计与转子动力学分析[D]. 合肥: 合肥工业大学, 2007. [16] 李冰. 电磁轴承系统集成化技术的研究[D]. 南京: 南京航空航天大学, 2003. [17] 郑坚强. 电磁轴承有限元分析结构设计及控制[D]. 杭州: 浙江大学, 2004. [18] 周传月,JOACHIM S. 基于Madyn_2000的磁悬浮轴承控制器的设计[J]. 风机技术,2017,59(4): 51-59,78.ZHOU Chuanyue, JOACHIM S. Design of magnetic bearing controller using Madyn 2000[J]. Chinese Journal of Turbomachinery, 2017, 59(4): 51-59,78. [19] 杨先,李洪斌. 有机朗肯循环发电技术的原理及应用[J]. 大众科技,2020,22(9): 61-63. doi: 10.3969/j.issn.1008-1151.2020.09.019YANG Xian, LI Hongbin. Principle and application organic Rankine cycle power generation technology[J]. Popular Science & Technology, 2020, 22(9): 61-63. doi: 10.3969/j.issn.1008-1151.2020.09.019 [20] HERZOG R, BUHLER P, GAHLER C, et al. Unbalance compensation using generalized notch filters in the multivariable feedback of magnetic bearings[J]. IEEE Transactions on Control Systems Technology, 1996, 4(5): 580-586. doi: 10.1109/87.531924 [21] 宋腾,韩邦成,郑世强,等. 基于最小位移的磁悬浮转子变极性 LMS 反馈不平衡补偿[J]. 振动与冲击,2015,34(7): 24-32.SONG Teng, HAN Bangcheng, ZHENG Shiqiang, et al. Variable polarity LMS feedback based on displacement nulling to compensate unbalance of magnetic bearing[J]. Journal of Vibration and Shock, 2015, 34(7): 24-32. 期刊类型引用(5)
1. 杨杰，杨帅，周发助，李智鑫，宋成林. 多点永磁电磁混合悬浮实验平台设计. 实验技术与管理. 2025(03): 153-160 . 
百度学术2. 李至勇，杨杰，周发助，胡海林. 电磁悬浮系统的改进滑模控制方法. 电子测量技术. 2024(04): 87-94 . 百度学术3. 张明，谢延松，李洪涛，孙凤，徐方超，张磊. 一种高静-低动刚度磁弹簧的建模与特性分析. 西南交通大学学报. 2023(04): 933-939+946 . 
本站查看4. 张晓友，尚利冬，赵川，孙凤，金俊杰，徐方超. 混合磁悬浮平台的零功率特性分析. 机械工程与自动化. 2023(06): 40-42+45 . 百度学术5. 魏发南，刘英. 面向微细加工的二自由度磁悬浮平台. 西南交通大学学报. 2023(06): 1318-1327 . 本站查看其他类型引用(6)
-
下载:
下载:
百度学术
点击查看大图
计量
- 文章访问数: 392
- HTML全文浏览量: 274
- PDF下载量: 25
- 被引次数: 11