Independent Cascade Control Method for Permanent Magnetic Levitation Platform
-
摘要:
针对各悬浮单元磁力特性差异导致的倾斜问题，设计一种用于无接触传送的永磁悬浮平台，提出了具有气隙偏差积分反馈的分散式串级控制方法. 首先，分析了悬浮平台的动力学模型与平衡条件，获得各悬浮单元气隙与平台三自由度之间的变换关系；其次，根据悬浮单元系统特点设计了双闭环串级控制系统，外环是以气隙为被控对象的主调节回路，内环是以转角为被控对象的随动回路，并引入气隙偏差进行积分反馈；最后，通过悬浮实验进行验证. 结果表明：引入积分反馈后，起浮后各磁悬浮单元气隙一致，向不同位置施加0.1 kg重物，各悬浮单元气隙同步增大0.12 mm；悬浮平台能够通过调节内环转角弥补磁力特性差异并实现偏载下水平悬浮，但系统的调节时间比无积分反馈的分散串级控制系统增加约1.4倍.
Abstract:Aiming at tilt caused by the difference of magnetic characteristics of each MLU (magnetic levitation unit), a permanent magnetic levitation platform is designed for contactless transportation; and an independent cascade control method is presented, which allows integral feedback of air gap deviation. Firstly, the dynamic model and equilibrium conditions of the levitation are analyzed to establish the transformation relationship between the air gap of each suspension unit and three degrees of freedom of the platform. Secondly, given each MLU features, the cascade control system with a double closed-loop is designed, in which the outer loop plays a major regulating role for the air gap, the inner loop works as a follow-up control for the rotation angle, and air gap error among each MLU runs as the integral feedback. Finally, the levitation experiment is conducted to verify the proposed controller. The results show that, when introducing the integral feedback, the air gap of each MLU is the same after levitating, and the air gap increases 0.12 mm synchronously when 0.1 kg of weight is applied to different positions. The platform adjusts the rotational angle to control the difference of magnetic characteristics and fulfill the horizontal levitation under eccentric load. However, the adjusting time of the system is about 1.4 times that of the independent cascade control system without integrated feedback.
-
Key words:
- magnetic levitation /
- PD control /
- cascade control systems /
- integrated feedback
-
随着三维扫描仪等硬件设备的成本降低和市场增长，其应用越来越广泛. 在三维扫描应用场景中，由于每个摄像头或扫描仪在其自身的空间中生成点云，而并非在对象空间中，所以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. Three-dimensional model of permanent magnetic levitation platform
图 2 永磁悬浮平台的运动分析简图
Figure 2. Schematic diagram of motion analysis for permanent magnetic levitation platform
图 3 具有偏差积分反馈的分散串级控制系统框图
Figure 3. Block diagram of independent cascade control method with integral feedback of air gap deviation
图 4 永磁悬浮平台的实验系统
Figure 4. Experimental system of permanent magnetic levitation platform
图 5 无积分反馈时系统起浮实验结果
Figure 5. Experimental results during system floating without integral feedback
图 6 有积分反馈时系统起浮实验结果
Figure 6. Experimental results during system floating with integral feedback
图 8 无积分反馈时在Ⅰ加载的实验结果
Figure 8. Experimental results of loading at position Ⅰ without integral feedback
图 9 有积分反馈时在Ⅰ加载的实验结果
Figure 9. Experimental results of loading at position Ⅰ with integral feedback
图 10 悬浮平台的气隙阶跃实验结果
Figure 10. Experimental results of air gap step for levitated platform
表 1 施加偏载后悬浮系统的稳态气隙
Table 1. Steady-state air gap length of levitation system under eccentric load
mm 加载位置编号 无气隙偏差补偿 有气隙偏差补偿 1 2 3 4 1 2 3 4 0 3.94 3.92 4.11 4.17 4.00 4.00 4.00 4.00 Ⅰ 4.15 4.02 4.08 4.24 4.12 4.12 4.12 4.12 Ⅱ 4.05 4.16 4.21 4.11 4.12 4.12 4.12 4.12 Ⅲ 3.90 4.0 4.33 4.26 4.12 4.12 4.12 4.12 Ⅳ 4.00 3.87 4.20 4.39 4.12 4.12 4.12 4.12 
下载: 导出CSV
表 2 施加偏载后永磁体的转角
Table 2. Rotational angle of permanent magnet under eccentric load
（°） 加载位置编号 无气隙偏差补偿 有气隙偏差补偿 1 2 3 4 1 2 3 4 0 59.5 52.9 55.2 59.1 58.4 56.3 53.1 59.6 Ⅰ 64.3 54.7 54.6 60.9 61.4 57.8 53.6 61.1 Ⅱ 61.9 57.4 57.0 57.7 60.4 58.6 54.5 59.7 Ⅲ 58.6 54.3 59.3 61.2 58.8 57.7 55.3 60.8 Ⅳ 60.9 51.9 57.0 64.2 60.2 56.2 54.4 61.9
下载: 导出CSV
-
[1] 王念先,胡业发,吴华春,等. 大气隙混合磁悬浮轴承承载力模型的研究[J]. 机械工程学报,2015,51(1): 153-160. doi: 10.3901/JME.2015.01.153WANG Nianxian, HU Yefa, WU Huachun, et al. Research on bearing capacity model of large-air-gap hybrid magnetic bearings[J]. Journal of Mechanical Engineering, 2015, 51(1): 153-160. doi: 10.3901/JME.2015.01.153 [2] 邓自刚,王家素,王素玉,等. 高温超导磁悬浮轴承研发现状[J]. 电工技术学报,2009,24(9): 1-8. doi: 10.3321/j.issn:1000-6753.2009.09.001DENG Zigang, WANG Jiasu, WANG Suyu, et al. Research and development status of high temperature superconducting magnetic bearings[J]. Transactions of China Electrotechnical Society, 2009, 24(9): 1-8. doi: 10.3321/j.issn:1000-6753.2009.09.001 [3] LI H, DENG Z, HUANG H, et al. Experiments and simulations of the secondary suspension system to improve the dynamic characteristics of HTS maglev[J]. IEEE Transactions on Applied Superconductivity, 2021, 31(6): 1-8. [4] ZHU H Y, JOO T T, KHIANG P C. Magnetically levitated parallel actuated dual-stage (maglev-PAD) system for six-axis precision positioning[J]. IEEE/ASME Transactions on Mechatronics, 2019, 24(4): 1829-1838. doi: 10.1109/TMECH.2019.2928978 [5] 王莉,张昆仑. 基于零功率控制策略的混合磁悬浮系统[J]. 西南交通大学学报,2005,40(5): 667-672. doi: 10.3969/j.issn.0258-2724.2005.05.022WANG Li, ZHANG Kunlun. Hybrid magnetic suspension system based on zero power control strategy[J]. Journal of Southwest Jiaotong University, 2005, 40(5): 667-672. doi: 10.3969/j.issn.0258-2724.2005.05.022 [6] 张钢,孟庆涛,钟永彦,等. 五自由度全永磁轴承系统的稳定悬浮特性分析[J]. 机械工程学报,2015,51(5): 56-63. doi: 10.3901/JME.2015.05.056ZHANG Gang, MENG Qingtao, ZHONG Yongyan, et al. Stable levitation performance analysis of five degrees of freedom all permanent magnetic bearing system[J]. Journal of Mechanical Engineering, 2015, 51(5): 56-63. doi: 10.3901/JME.2015.05.056 [7] 胡坤,王爽,李德永,等. 一种永磁磁垫式带式输送机的气隙特性研究[J]. 煤炭学报,2015,40(3): 701-706.HU Kun, WANG Shuang, LI Deyong, et al. Research on air gap characteristics of permanent magnetic cushion belt conveyor[J]. Journal of China Coal Society, 2015, 40(3): 701-706. [8] OKA K, YAMAMOTO K, HARADA A. Magnetic suspension mechanism using rotary permanent magnets[J]. International Journal of Applied Electromagnetics and Mechanics, 2020, 64(1/2/3/4): 977-983. [9] AKABAYSHI A, MIZUNO T, TAKASAKI M, et al. Development of flux-path control magnetic suspension system using flux concentration[J]. Transactions of the Japan Society of Mechanical Engineers, 2013, 79(801): 1483-1494. doi: 10.1299/kikaic.79.1483 [10] 孙凤,韦伟,金俊杰,等. 可控磁路式并联型永磁悬浮系统[J]. 仪器仪表学报,2017,38(7): 1714-1722. doi: 10.3969/j.issn.0254-3087.2017.07.019SUN Feng, WEI Wei, JIN Junjie, et al. Parallel permanent magnetic suspension system with flux path control[J]. Chinese Journal of Scientific Instrument, 2017, 38(7): 1714-1722. doi: 10.3969/j.issn.0254-3087.2017.07.019 [11] 李强,唐敬虎,孙凤,等. 可变磁路式永磁悬浮系统的防跌落防吸附控制[J]. 仪器仪表学报,2019,40(3): 246-254.LI Qiang, TANG Jinghu, SUN Feng, et al. Anti-fall and anti-adsorption control of permanent magnetism levitation system with flux path control[J]. Chinese Journal of Scientific Instrument, 2019, 40(3): 246-254. [12] 孙凤,唐敬虎,李强,等. 可变磁路式永磁悬浮系统刚度特性分析及变刚度控制[J]. 振动与冲击,2020,39(7): 132-139.SUN Feng, TANG Jinghu, LI Qiang, et al. Stiffness characteristics analysis and variable stiffness control for a permanent magnetic levitation system with variable magnetic circuit[J]. Journal of Vibration and Shock, 2020, 39(7): 132-139. [13] ZHAO C, SUN F, JIN J, et al. Analysis of quasi-zero power characteristic for a permanent magnetic levitation system with a variable flux path control mechanism[J]. IEEE/ASME Transactions on Mechatronics, 2021, 26(1): 437-447. [14] 寇宝泉,张赫. 多自由度短行程超精密平面电机技术发展综述[J]. 电工技术学报,2013,28(7): 1-11. doi: 10.3969/j.issn.1000-6753.2013.07.001KOU Baoquan, ZHANG He. Development of multi-DOF micro-motion ultra-precision planar motors[J]. Transactions of China Electrotechnical Society, 2013, 28(7): 1-11. doi: 10.3969/j.issn.1000-6753.2013.07.001 [15] 李广,张鸣,朱煜,等. 六自由度磁浮微动台电磁场建模与计算[J]. 清华大学学报(自然科学版),2010,50(12): 1978-1982.LI Guang, ZHANG Ming, ZHU Yu, et al. Six degree-of-freedom maglev microstage electromagnetic field modeling and calculation[J]. Journal of Tsinghua University (Science and Technology), 2010, 50(12): 1978-1982. [16] 李黎川,丁玉成,卢秉恒. 超精密磁悬浮工作台及其解耦控制[J]. 机械工程学报,2004,40(9): 84-88,94. doi: 10.3321/j.issn:0577-6686.2004.09.018LI Lichuan, DING Yucheng, LU Bingheng. High-precision magnetically suspended table and its decoupling control[J]. Chinese Journal of Mechanical Engineering, 2004, 40(9): 84-88,94. doi: 10.3321/j.issn:0577-6686.2004.09.018 [17] 陈启会,李群明,胥晓. 平面型磁悬浮平台的结构与分散控制[J]. 机械设计与研究,2013,29(2): 21-25. doi: 10.3969/j.issn.1006-2343.2013.02.006CHEN Qihui, LI Qunming, XU Xiao. Structure and decentralized control of planar magnetic levitation platform[J]. Machine Design & Research, 2013, 29(2): 21-25. doi: 10.3969/j.issn.1006-2343.2013.02.006 [18] 孙兴伟,路英园,孙凤,等. 可变磁路式永磁悬浮系统的动力学特性研究[J]. 中国机械工程,2014,25(20): 2782-2787. doi: 10.3969/j.issn.1004-132X.2014.20.015SUN Xingwei, LU Yingyuan, SUN Feng, et al. Research on dynamics performance of permanent magnet suspension system with variable magnetic circuit[J]. China Mechanical Engineering, 2014, 25(20): 2782-2787. doi: 10.3969/j.issn.1004-132X.2014.20.015 期刊类型引用(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)
-
下载:
下载:
百度学术
点击查看大图
计量
- 文章访问数: 234
- HTML全文浏览量: 215
- PDF下载量: 13
- 被引次数: 11