首页/文章/ 详情

UM-GearEccDataset | 澳门大学齿轮偏心故障数据集

3天前浏览321

齿轮箱是使用范围很广的部件,而且工况复杂多变,目前大部分公开数据集是集中在断齿、缺齿、齿根裂纹和表面磨损等齿面的故障,对齿轮偏心故障研究较少。本期给大家推荐一个包含不同转速不同传动系统声音电流电压振动数据的齿轮偏心数据集,大大提高了诊断难度。该数据集由澳门大学杨志新教授SMART团队在2024年公开的,写论文的同学们赶紧用起来吧!

论文基本信息

论文题目A comprehensive gear eccentricity dataset with multiple fault severity levels: Description, characteristics analysis, and fault diagnosis applications

论文期刊Mechanical Systems and Signal Processing
论文日期2024
论文链接https://doi.org/10.1016/j.ymssp.2024.112068
数据链接https://github.com/LeeJMJM/UM-GearEccDataset
作者Jiaming Li (a), Hao Chen (a), Xian-Bo Wang (b), Zhi-Xin Yang (a),∗
 
机构

a State Key Laboratory of Internet of Things for Smart City (UM), Department of Electromechanical Engineering and Centre of Artificial Intelligence and Robotics, University of Macau, 999078, Macao Special Administrative Region of China.  

b The Hainan Institute of Zhejiang University, Sanya, 572025, China.

团队介绍澳门大学杨志新教授SMART团队,长期从事机电设备的安全监控工作。主要研究方向有数据驱动的故障诊断和寿命预测、安全监控导向的机器人控制和计算机视觉等。

注:上述链接若打不开,文末有百度网盘链接。

 摘要

我们开发了一个名为UM-GearEccDataset的齿轮偏心多故障程度数据集,以促进故障机理的研究和数据驱动的故障诊断。其他现有数据集没有充分考虑齿轮偏心诊断的故障严重程度FSL, fault severity level),为弥补这一差距,我们提出了一种新颖的齿轮偏心模拟结构,实现了FSL的连续调整。该综合全面的数据集包含多种类故障信号,通过记录五种传感器收集的11个通道信号考虑到了各实验变量传动系结构转速FSL同时故障多模态信号。这个丰富的数据集可成为各种研究应用的宝贵资源。

在获得了该数据集之后,我们还对UM-GearEccDataset进行了细致的检查,以尽量消除其可能含有的其他现有数据集中广泛存在的可靠性问题:

  • 1)首先,检查数据本身。通过分析信号的频谱,计算特征频率与FSL之间的相关系数,并研究不同变量的影响,研究信号特征是否合理。
  • 2)然后,通过应用卷积神经网络(CNN)和梯度加权类激活映射加加(GradCAM++)等深度学习技术来验证数据集的可靠性。FSL的分类任务由CNN模型完成,以分析诊断准确率随试验变量的变化。GradCAM++实现了显著性分析,以找出输入频谱的哪些区域贡献更大。

结果表明,该数据集具有明显的齿轮偏心故障特征。不同信号的特征和所有试验变量的影响也是合理的。因此,所提出的数据集以其精心的设计和可靠性,可以加强各种新兴的智能故障诊断研究,为该领域的进一步研究提供坚实的基础。  

关键词故障诊断,多传感器信号,直齿轮,齿轮偏心,数据集,频谱分析,显著性分析

目录

1 引言
2 数据集建立
    2.1 偏心齿轮结构
    2.2 试验台
    2.3 实验计划和过程
3 信号本身的分析
    3.1 信号总览
    3.2 高敏感信号
    3.3 低敏感信号
    3.4 信号检查与分析结论
4 该数据集的故障诊断示范
    4.1 基准CNN模型的诊断准确率
    4.2 GradCAM++的显著性分析
5 结论    

1 引言

齿轮广泛应用于各行各业,其故障诊断的重要性无需多言。但关于齿轮的故障诊断往往集中于破齿、缺齿、齿根裂纹和表面磨损等齿面的故障,关于偏心故障这种更系统性的故障研究较少齿轮偏心指的是齿轮实际和理论的旋转中心不重合的一种现象。机加工公差、来料质量、热处理过程等对此均有贡献。而且,机器装配后长轴的弯曲和轴承的劣化更是使得齿轮偏心难以避免。偏心后的齿轮在振动噪声等各方面都会得到恶化,因此是一个相当需要引起重视的问题。现有的理论研究多从大幅简化的动力学模型和有限元模型出发来分析故障特征,其中一个结论便是:齿轮偏心会在振动频谱上造成故障齿轮的啮合频率及其倍频周围出现大量间距为转频的边带成分。然而,目前却缺乏相关的坚实实验研究,而这正是本文将要填补的一个窗口。

另一研究窗口出现在对故障程度的诊断。表1列举了一些学界常用的数据集,可以看到他们的故障严重程度(FSL, fault severity level)数量不超过5。事实上,这些数据集都往往被用于故障诊断的分类任务。直觉上说,用这种很有限FSL的数据集去做FSL的故障诊断研究,是不够好的。现有的数据集FSL很有限,大概是因为实验不好做,因为一般一种故障程度就得对应一个故障件,相当麻烦。而那些加速寿命实验,往往又很难去控制最后的故障件坏成什么样子。因此,一个包含了大量故障严重程度的数据集是非常有利于相关研究的。

表1 一些数据集FSL数量比较

除去上述两个主要的问题之外,现有用于故障诊断的数据集还有以下不足:

  • 1)信号源不足。对于近年来很火的多信号融合诊断,很多数据集提供的信号信道不够多,信号种类也不够丰富。
  • 2)没考虑实际的噪声干扰。在进行诊断模型的抗噪性研究时,往往通过给信号添加白噪声和粉红噪声等噪声来实现,没考虑到现实中那些分布更加复杂多变的噪声和干扰。
  • 3)未对数据本身进行足够的检查来避免隐藏变量。隐藏变量指的是机器学习中意外引入的变量,应尽量避免。比如,有文献指出CWRU数据集中存在非一致的电磁噪声、机械部件的松动和滚动体间卡滞的材料,难说这样训练出来的诊断模型学到的是真的故障特征还是这些隐藏变量。所以数据集本身的可靠性需要得到验证。

为了解决上述问题,我们提出了UM-GearEccDataset,一个考虑了大量实验变量和历经严格检查的全面数据集。本文贡献如下:

  • 1)我们提出了一种新颖的齿轮偏心故障模拟装置,并基于此开发了一个全面的数据集。该装置可以帮助验证现有的齿轮偏心理论研究;该数据集则可以用于数据驱动的诊断方法研究。
  • 2)我们对数据集本身的信号进行了深入的分析以保证数据集的可靠性。这些检查包括反映了偏心故障的故障特征和设定各实验变量的影响模式。结果表明隐藏变量基本不存在,研究者们可以放心地使用该数据集而无需担心可靠性的问题。
  • 3)另外我们还用两个深度学习的案例来侧面验证了数据集的可靠性并展示了数据集的价值。从诊断准确率上来看,准确率是和数据检查部分的推测重合的。显著性分析则表明故障相关频带对诊断结果的影响很大。

2 数据集建立

为了收集不同FSL的齿轮偏心实验数据,我们提出了一种新颖的偏心齿轮结构,可以在设计的最大故障程度内调节任意FSL。基于该结构,搭建了实验平台并布置了11个传感器。通过深思熟虑的实验方案,考虑到了以下实验变量齿轮偏心的FSL转速传动系结构信号来源偏心齿轮是否同时还有破齿故障。包括实验数据,一些其他补充材料也一并发布于https://github.com/LeeJMJM/UM-GearEccDataset

(为了更好的阅读体验,以下文章中的部分细节请读者参考原文,这里仅作简介)

2.1 偏心齿轮结构

所提出的结构主要原理是改变齿轮旋转中心和节圆中心的距离。为实现该效果,齿轮被设计成具有一定的偏心量,并且在齿轮和轴之间设计了一个偏心轴套,如图1所示。齿轮和轴套间的轴承仅为后续调节方便,与上述原理无关。轴和齿轮的锁紧螺钉可在偏心调节完毕后拧紧以保持所调整的偏心量,即FSL。  

图1 所提出的偏心齿轮结构

具体来说,FSL的调整通过改变如图1 (b)中的所示的角度θ来实现。θ=0°时理论上不存在偏心,θ=180°时理论上FSL最大。其中,在偏心轴套的一侧设计有刻度盘,在齿轮的同侧设计有刻度,如图2 (a)所示。

(具体细节请读者参考原文)

图2 偏心齿轮的照片

2.2 试验台

试验台由一疲劳寿命试验台改装而来,如图3所示。负载电机的联轴器已拆除,来模拟一种齿轮负载很小的情况。测试齿轮箱中的原95T齿轮被替换为了上一节所述的偏心FSL可调节的齿轮,如图3中标红的齿轮所示。

图3 试验台的布置。偏心齿轮已标红

整个传动系包含一个驱动电机和三个齿轮箱。测试齿轮箱为减速箱负载齿轮箱1和2为增速箱
为了采集多源信号,共布置了11个传感器。他们的位置和参数等信息请参考图4和表2。  

图4 传感器位置

表2 数采和传感器的信息

2.3 实验计划和过程

该数据集所包含的实验变量于表3中列出。包括理论无偏心的情况,共有11种FSL(表示为E00到E20)。有4种电机转速(表示为S600到S1500)。有4种传动系结构(表示为T1到T4,如图3中彩框圈出的部分):
  • T1:移除测试齿轮箱的36T齿轮,传动系终止于轴2

  • T2:移除联轴器1。传动系终止于测试齿轮箱的输出轴;

  • T3:移除联轴器2。传动系终止于负载齿轮箱1的输出轴;

  • T4:传动系终止于负载齿轮箱2的输出轴。

此外,还准备了两个偏心齿轮,一个是没有破齿故障的(表示为G1)。另一个则是有的(表示为G2),如图2 (b)所示。上述四种实验变量也称为“工况”(operating condition)。实验中,每个工况运行一次试验台,记录一个文件。最后一个变量——信号源,在文件中由11个通道的信号同时记录

表3 该数据集所考虑的实验变量

实验总共记录了 4 × 4 × 11 × 2 = 352 个文件,每个文件时长600 s。以文件N139_G1_E16_T3_S1200.mat为例。其序号为139。工况为:无破齿齿轮,理论偏心量为0.16 mm传动系结构3,理论驱动电机转速1200 rpm

3 信号本身的分析

3.1 信号总览

数据集中所存储的所有信号,其原始数值都已转为了真实的物理单位。各传感器的灵敏度见表2。

这里取两个工况来总览信号:N352_G2_E20_T4_S1500 和 N1_G1_E00_T1_S600,如图5所示。总体上来说这些信号的数值都在合理的范围内。虽然这两个工况的时域波形相差很大,但频域上都有个共同点,即8000 Hz及其倍频的电机开关频率明显可见。除了传感器⑨到⑪因为数采的采样频率不足频谱到不了8000 Hz。

图5 信号总览。请注意那些合理的特征,如合理的数值范围、波形上同步的波纹、电机开关频率

另一个明显的特征是③到⑤的蓝色时域波形波纹明显同步。同样地,⑨到⑪也明显同步。但是遗憾的是受限于实验硬件条件(具体细节请读者参考原文),这两组信号之间是没有同步的
一个常见的展示两个信号相关性的方法是将两个信号分别作为横纵坐标数值直接绘出散点图。若这两信号有完美的线性关系,则会出现沿直线分布的散点,比如一个信号和它自己。前面说到该数据集的信号并没有专门做同步处理,故将频谱如此绘图而非原始信号直接绘图,得到各传感器信号频域之间的线性关系如图6所示。从图中可以看出,这些信号大体可以分成两种。
  • 红色区域的信号彼此比较相关,大部分散点都位于对角线周围区域;
  • 绿色区域的信号彼此十分相关,散点基本上都位于对角线上;
  • 黄色的区域信号相关性低,散点基本都分布在横纵坐标轴处。
如此,可将这些信号分成两种,且称之为高敏感信号和低敏感信号吧。

 

图6 各传感器信号频谱之间的散点图,其中传感器信号标于各横纵坐标轴。各个被标红、绿、黄的部分,代表高敏感信号之间、低敏感信号之间、高低敏感信号之间的散点图

为了更有逻辑地刻画这些信号的特点,本文接下来从高敏感信号和低敏感信号两个方面来介绍。总体来说,高敏感信号有更强的故障特征,并且故障诊断中的准确率也会更高,这些后续将会一一印证。

3.2 高敏感信号

高敏感信号指的是5个加速度传感器的信号(③到⑤,⑨和⑩),1个麦克风的信号(⑥)和一个转矩传感器的信号(⑪)。

(以下的内容比较繁复,具体细节请读者参考原文)

3.2.1 FSL的影响

现有文献中的理论分析表明,齿轮偏心时会产生一些故障特征频率,并且偏心越严重这些频率的幅值会越大。以工况G1_T3_S1500的振动信号为例,如图7所示,这种趋势大致是存在的。

图7 偏心齿轮啮合频率周围的频谱,明显可见故障特征频率

为了更好研究上述趋势,我们将一些故障特征的幅值,称为“故障频率特征”,作为纵坐标,FSL的数值作为横坐标,作图,如图8所示。并且,常见的一些时域特征,如均方根值、峭度、均值,也一并绘出以作比较。其中阴影代表±一个std的范围。  

图8 两种特征对FSL的变化趋势。其中一种展现出了更强的单调性

然后,我们接着用了Spearman’s rank correlation coefficient (SRCC)来定量分析这种单调性。用大量表格列出高敏感性信号在不同工况下的SRCC,结果说明那三个时域特征单调性都不如故障频率特征。且偏心齿轮转频的一倍边带这种故障频率特征单调性最高,也符合现有文献中的结论。(这里略去了表格内容和具体的数值对比)
因此,可以总结,从故障频率特征的角度来看,该数据集是合理的,符合其他做理论分析的文献的描述。且没有发现明显的隐藏变量。

3.2.2 转速的影响

转速最明显的就是特征频率的变化,如图9所示。  

图9 不同转速下频谱的变化。注意转速增加后的边带范围会散开

而且RMS会显著受到转速的影响,如图10所示。  

图10 不同信号的RMS和mean的散点图。注意:1) 更高的转速产生更大的RMS。2) 更长的传动系结构产生更高的转矩信号mean

此外,声压级的计算结果表明,声音信号的声压级在75到82dB的范围。长期无防护地暴露在这样的实验条件下确实对人体有害。
另一个细节是,通过观察特征频率来推测实际的驱动电机转速,发现推测出的转速并不严格吻合设定转速,会大个若干rpm。(具体细节请读者参考原文。这也算是该数据集的不足,未用编码器、转速计等设备记录真实的转速

3.2.3 传动系结构的影响

引入传动系结构这一变量主要是想添加更贴近实际的噪声干扰。一个实际的场景就是,厂家有大量的某设备的故障诊断数据,但使用者会在该设备上加装一些结构,比如添加一些变速装置。这些添加的结构会带来新的频率成分且会改变负载状态。多齿轮箱结构在现有文献中很常见,但并未见到变传动系长度的研究。以该数据集为例,不同传动系结构的特征频率如表4所示。

表4 不同传动系结构下的齿轮啮合频率和轴旋转频率

接下来用一个振动信号和一个声音信号,在齿轮啮合频率和轴旋转频率的频谱中说明了:
  • 有些特征频率振动信号中很明显,声音信号中不明显;
  • 有些特征频率则反过来在声音中明显,振动中不明显。
说明不同信号的特征某种程度上说是互补的,这也是一个该数据集要记录多信号源数据的原因。(这里略去了若干频谱的图片和一些具体内容的解释)
不同传动系结构下的RMS和mean,如图11所示。  

图11 不同信号的RMS和mean的散点图。注意:1) 转矩信号有一定的零漂。2) T4传动系结构下转矩的mean格外高

3.2.4 缺齿故障的影响

引入缺齿故障这一变量是想模拟多故障同时发生的场景。另一原因是想降低数据的可诊断性,因为现有文献指出缺齿故障的故障特征频率和偏心相同,两者同时发生时势必会彼此影响,使得诊断更困难。
类似前面对G1齿轮的故障频率特征趋势分析,我们同样对G2列出了SRCC的表格。结果表明,引入缺齿后,SRCC也出现了下降。可以理解为可诊断性也出现了下降。(这里略去了表格和一些具体数值的比较)
此外,缺齿势必也会导致振动信号RMS的增加。这里列出了具体的RMS数值,如表5所示。结果确实也是这样的。

表5 G1和G2的RMS数值的比较。两者中较大的数值已加黑

3.3 低敏感信号

低敏感信号包括两个电流钳(①和②)采集的信号和两个驻极体传感器(⑦和⑧)采集到的信号。这两种传感器在旋转设备的状态监测里相比于高敏感信号确实更少用到。经过这一节的分析,它们的故障特征确实没那么明显,但各实验变量对他们的影响依然是合理的,且没有发现明显的隐藏变量。

电流钳测量的是驱动电机U/V两相的电流。虽然电流在电机的故障诊断中很常见,但一般做的都是电机本身的故障,如绕组断路短路、导条断裂、电池轴承损坏等。但当做的是电机之外的部件,比如整个传动系中其他轴承的故障,准确率往往低很多,比如发布PU数据集的会议论文中的描述。因此,可以认为,用UM-GearEccDataset中的电流信号来做故障诊断应该是很有挑战性的任务。
驻极体传感器本质上相当于电容。少有几篇文献用这种传感器来做设备的故障诊断,因为这种传感器采集到的信号毕竟与振动相比相差太多。但这种传感器有其优势:价格低、轻薄、可以自供电。相信在物联网场景下这种传感器有其不可替代的独特优势。(该传感器的细节请参考原文)

3.3.1 电流传感器和驻极体传感器信号间的相似性

总的来说,这两种传感器是很像的。这里绘出放大的频谱,如图12所示。  

图12 电流传感器和驻极体传感器信号的放大频谱。可以看到在    及其倍频的相似性

然后,我们用欧式距离和余弦距离计算了各信号频谱之间的距离,发现电流传感器驻极体传感器信号这两种信号的距离确实很近。这两种信号的距离比两个振动传感器之间的距离还近。(计算过程和相关数据请参考原文)

因此,简洁起见,接下来我们仅对传感器①进行分析

3.3.2 实验变量的影响

总体来说,在高敏感信号里找到的那些故障特征频率在低敏感信号中并不明显,甚至找不到。比如在图13和图14中,间隔为转频的边带可以找到但齿轮啮合频率很难找到。  

图13 T4工况下两种速度的电流频谱低频部分。    的边带特征没有振动频谱的啮合频率的边带特征稳定  

图14 S600工况下两种传动系结构的电流频谱低频部分。    的边带特征没有振动频谱的啮合频率的边带特征稳定

转速的影响是显然的,就像图9一样。但其他明显且有物理意义的特征就不太好找了。
这里利用电流频谱中的频率信息,估计了电机的滑差率,其数值基本上可以认为正比于电机的输出扭矩。结果发现在T4工况下,滑差率明显大于其他三个工况。这与3.2.3小节中对转矩信号的分析是吻合的。(计算过程和相关数据请参考原文)

3.4 信号检查与分析的结论

这里列了一张表格,用了展示上述对信号检查和分析的结论,如表6所示。该表总结了支持数据集可靠性的结论:数据有合理的故障特征且无明显的隐藏变量,使用者可以放心使用而无需担心可靠性的问题。

表6 数据检查的简短结论

4 该数据集的故障诊断示范

前一章针对数据本身进行了检查,这一章利用两个深度学习的案例来分析数据集。

4.1 基准CNN模型的诊断准确率

除了FSL外,该数据集还考虑了其他4中实验变量:转速、传动系结构、信号源、同时存在破齿故障。在这一节中,这4个变量的每一种组合都用相应数据训练一个CNN诊断模型来开展11个FSL的分类任务。然后,分类准确率和混淆矩阵用来分析这些变量的影响、数据的可诊断性、数据集的可靠性。
首先,数据进行了分段和信号处理,用频谱作为输入。然后,一部分数据用来训练模型,另一部分则用来测试。(具体细节请参考原文)
图15展示了G1和G2两种工况下不同传感器的测试准确率,可以得到以下结论:
  • 1)G2情况下的准确率确实低一些,对应着2.2.4小节中所述的引入破齿会降低数据可诊断性的推测。

  • 2)高敏感信号的准确率基本都大于90%,与低敏感信号显著不同。

  • 3)(一些其他描述性的结论这里略去)

 

图15 G1和G2两种工况下不同传感器的测试准确率

该数据的另一个特点是相近FSL的数据分布差异比较远FSL的数据小。这一点与常见的不同故障类型的数据集很不一样。这一特点在混淆矩阵中可以看出,如图16所示。从这一点上,也可以侧面反映出该数据集的可靠性。

图16 所有转速和传动系结构条件下的混淆矩阵。注意到误分类的样本大多按对角线分布

该数据集的一个优势在于其使用了很多数量和类型的传感器。这样一来,该数据集就具有较大的可诊断性范围。以测试准确率为例,与其他数据集的比较如图17所示。研究者们可以按照各自的需求自行选择不同的信号来验证所提出的诊断方法。  

图17 与其他数据集的可诊断性比较。我们提出的数据集具有更大的可诊断性范围

转速和传动系结构对测试准确率的影响并没有明显的规律。但这里要强调传动系结构这一变量的意义。传动系结构对信号的影响远比简单给信号添加白噪声这么简单。如今研究中对故障诊断任务的迁移任务研究火热,常集中于同一设备跨工况的迁移任务和跨设备的迁移任务。大家的共识是前者过于简单且实际意义不够大,后者具有很强的实际意义但难度很高。这里,我们提供了一种折中的任务,即跨传动系的迁移。我们用T1的数据训练了CNN,再用T1到T4的数据直接做测试,称为Task 1到Task 4,准确率结果如图18所示。可以看到,准确率几乎是随着传动系的变长而单调下降的。数据集的该特点,不但可以为那些研究抗噪性的诊断方法提供一个除了单纯添加白噪声、粉红噪声等理想噪声外的验证途径,也在跨转速和跨设备迁移任务之间填补了一种空白。  

图18 跨传动系结构任务的准确率结果。更长的传动系结构准确率更低

4.2 GradCAM++的显著性分析

类激活映射(CAM)是一种显著性分析方法,在分类任务中用来指出模型输入的那些部分对指定分类结果的重要程度。CAM相关技术近年来在故障诊断中运用非常广泛,在原始信号输入、频谱输入、时频特征输入等模型中都可以用来指出这些输入的高激活区域。这里以转矩传感器⑪的信号为例,说明信号故障特征。我们对11个FSL对应标签的数据分类结果分别使用CAM分析对应标签的CAM结果,即重要性权重,如图18所示。
由于转矩信号的特点(在3.2.3节中分析过),T4条件下    极为突出。但显著性分析的结果表明,诊断模型其实更关心齿轮啮合频率及其倍频的区域,并非是    所在的区域。这说明:
  • 1)输入频谱数值高的地方未必一定会产生高的重要性权重。

  • 2)按照CAM技术对模型的诠释,其实故障特征频段对分类结果的重要性更大。

这些都从侧面说明了数据具有合理的故障特征。

图18 不同标签频谱的重要性权重。      虽然在频谱中很突出,但并非在重要性权重中总是明显。相对地,齿轮啮合频率及其倍频的区域却总是很明显

此外,我们还展示了更多工况和传感器的结果,结论是类似的。(详情请参考原文)

5 结论

本文提供了一个大量故障严重程度的齿轮偏心故障数据集,提出了一种新颖的齿轮偏心模拟结构。实验结果可以支持大量的齿轮偏心理论研究,数据集则对数据驱动的故障诊断方法研究大有益处。

此外,我们还对数据从两方面进行了细致的检查:

  • 1)对数据本身,用故障特征相关系数实验变量的具体影响,说明了数据集是可靠的。

  • 2)从深度学习技术的角度出发,分类结果和显著性分析的结果也从侧面说明了数据集的可靠性。

总之,该数据集全面而可靠,研究者们请放心使用。
使用该数据集请引用:J. Li, H. Chen, X.-B. Wang, and Z.-X. Yang, “A comprehensive gear eccentricity dataset with multiple fault severity levels: Description, characteristics analysis, and fault diagnosis applications,” Mechanical Systems and Signal Processing, vol. 224, p. 112068, 2025.

编辑:李正平

校核:陈凯歌、赵栓栓、曹希铭、赵学功、白亮、任超

该文资料搜集自网络,仅用作学术分享,不做商业用途,若侵权,后台联系小编进行删除


来源:故障诊断与python学习
ACTMechanicalSystem振动疲劳断裂旋转机械航空UGpython声学裂纹理论电机材料传动机器人数字孪生
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-11-13
最近编辑:3天前
故障诊断与python学习
硕士 签名征集中
获赞 64粉丝 66文章 140课程 0
点赞
收藏
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈