CN116451580A - 一种滚动轴承剩余使用寿命区间预测方法 - Google Patents

一种滚动轴承剩余使用寿命区间预测方法 Download PDF

Info

Publication number
CN116451580A
CN116451580A CN202310407813.1A CN202310407813A CN116451580A CN 116451580 A CN116451580 A CN 116451580A CN 202310407813 A CN202310407813 A CN 202310407813A CN 116451580 A CN116451580 A CN 116451580A
Authority
CN
China
Prior art keywords
prediction
rolling bearing
interval
predicting
algorithm
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310407813.1A
Other languages
English (en)
Other versions
CN116451580B (zh
Inventor
周萌
蔡畅
王晶
王昶
史运涛
董哲
翟维枫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianchang Chajian Building Installation Decoration Co ltd
Original Assignee
North China University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China University of Technology filed Critical North China University of Technology
Priority to CN202310407813.1A priority Critical patent/CN116451580B/zh
Priority claimed from CN202310407813.1A external-priority patent/CN116451580B/zh
Publication of CN116451580A publication Critical patent/CN116451580A/zh
Application granted granted Critical
Publication of CN116451580B publication Critical patent/CN116451580B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Computer Hardware Design (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Tourism & Hospitality (AREA)
  • Game Theory and Decision Science (AREA)
  • Geometry (AREA)
  • Development Economics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Artificial Intelligence (AREA)
  • General Business, Economics & Management (AREA)
  • Medical Informatics (AREA)

Abstract

本发明公开了一种滚动轴承剩余使用寿命区间预测方法。本发明中,以滚动轴承作为研究对象,以滚动轴承的剩余寿命区间预测为研究目标。在实际应用中,由于各种因素,如噪声传感器数据、可变的操作条件、未知的故障模式等,预测结果存在各种不确定性,点预测得到的预测结果虽然直观,但只能在某一预测时间提供某一预测值,因此,需要进行区间预测。首先,基于改进型相关向量机进行剩余使用寿命点预测,采用改进秃鹰搜索算法进行对相关向量机的参数进行寻优,提高预测精确度。进一步,利用粒子滤波对改进秃鹰搜索算法优化相关向量机预测模型得到的预测结果进行过滤优化,以提高预测精确度。最后,基于点预测结果,利用核密度估计进行区间预测。

Description

一种滚动轴承剩余使用寿命区间预测方法
技术领域
本发明涉及数据驱动的故障预测与健康管理领域,具体是一种基于相关向量机和核密度估计的滚动轴承剩余使用寿命区间预测方法。
背景技术
滚动轴承是旋转机械中最重要的部件之一。复杂的工作环境和长时间的运转可能会在一定程度上造成轴承的损坏,这将进一步影响机械的运行。滚动轴承RUL区间预测可以为预测性维修决策提供理论依据,提前对故障维修做好规划,从而提高整个机械设备的运行周期,避免事故发生。因此,RUL区间预测一直是热门的研究课题之一,吸引了众多学者的关注。
本发明以滚动轴承作为研究对象,以滚动轴承的剩余寿命区间预测为研究目标。在实际应用中,预测结果存在各种不确定性,由于各种因素,如噪声传感器数据,可变的操作条件,未知的故障模式。点预测得到的预测结果虽然直观,但只能在某一预测时间提供某一预测值。由于滚动轴承工作条件多变,点预测难以表达预测结果的可靠性,仅依靠点预测结果难以为后续健康管理提供可靠依据。因此,需要对寿命预测进行不确定性研究,进行相应的区间预测。
发明内容
本发明的目的在于解决滚动轴承剩余使用寿命区间预测问题,提出一种基于相关向量机与核密度估计的滚动轴承剩余使用寿命区间预测方法。提取滚动轴承振动信号的时域特征RMS作为健康指标进行后续的寿命预测。其次,由于相关向量机的核函数参数对相关向量机的预测性能有很大影响,所以采用改进秃鹰搜索算法进行参数寻优,提高预测精确度。进一步,利用粒子滤波对改进秃鹰搜索算法优化相关向量机预测模型得到的预测结果进行过滤优化,提高预测精确度。由于实际应用不确定性的存在,点预测结果与实际结果存在一定程度的偏差。因此,在点预测结果的基础上,采用核密度估计方法进行剩余使用寿命的区间估计。
本发明提出了一种数据驱动的基于相关向量机与核密度估计的滚动轴承剩余使用寿命区间预测方法,包括:
步骤一:针对秃鹰搜索算法容易陷入局部最优的缺点,采用变异策略对秃鹰搜索算法进行优化形成MBES。为了提高预测精度,采用帐篷混沌策略对MBES算法进行优化。采用改进的秃鹰搜索(IMBES)优化算法优化相关向量机的核函数参数。
步骤二:利用IMBES对相关向量机的核函数参数进行寻优,为了进一步提高点预测的预测精度,利用粒子滤波算法对改进秃鹰搜索算法优化相关向量机预测模型的预测结果进行过滤优化,将IMBES-RVM预测模型代替粒子滤波的状态转换方程进行滤波处理,也进一步解决了粒子滤波状态转换方程依靠机理模型的缺陷,也提高了预测精确度。改进秃鹰搜索算法优化相关向量机结合粒子滤波流程如图2所示。
步骤三:传统的点预测不能处理系统运行中的不确定性。为了得到数值估计及其可靠性,实际应用中需要对预测区间进行计算。区间预测是指在一定置信区间内预测值的估计区间。因此,本发明在点预测结束后,利用核密度估计进行区间预测。将预测结果利用核函数拟合分布,确定预测上下限。综上所述,由于采用了上述技术方案,本发明的有益效果是:本发明中,以滚动轴承作为研究对象,以滚动轴承的剩余寿命区间预测为研究目标。在实际应用中,预测结果存在各种不确定性,由于各种因素,如噪声传感器数据,可变的操作条件,未知的故障模式,点预测得到的预测结果虽然直观,但只能在某一预测时间提供某一预测值。采用改进秃鹰搜索算法进行参数寻优,提高预测精确度。进一步,利用粒子滤波对改进秃鹰搜索算法优化相关向量机预测模型得到的预测结果进行过滤优化,提高点预测精确度。在点预测结果的基础上,利用核密度估计进行区间预测。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1为本发明中整体系统的框架图;
图2为改进秃鹰搜索算法优化相关向量机结合粒子滤波流程图;
具体实施方式
以下结合附图和具体实施例对本发明作详细描述。
在步骤一中,滚动轴承数据集利用加速度传感器测量,即在滚动轴承的横向和纵向安装两个高灵敏度加速度传感器,分别测量滚动轴承水平方向和垂直方向的振动加速度数值。利用同一型号的滚动轴承水平振动信号进行剩余使用寿命预测。
对于相关向量机,建立模型如下:
设定输入为xi,输出为yi,M是数据的个数,退化过程描述如下:
其中,是y(xii)的回归值,x是输入向量,ξi~N(0,σ2)是独立同分布噪声,ω是权重向量,Φ是设计矩阵,K(x,xi)是核函数。
其中,η是需要被优化的核参数,采用SSA来优化核参数。在RVM中ω和ξi能被直接估计,但是,如果直接估计ω和ξi会造成过拟合,因此,引入超参数αi,假设ωi服从均值为0,方差为的高斯分布,
其中,α=[α01,…,αM]。
权重ω的方差V=[V0,V1,…,VM],均值μ=[μ01,…,μM]的计算公式为:
V=(σ-2ΦTΦ+A)-1
μ=σ-2TyT
其中A=diag(α01,…,αM),α和σ2是需要被更新的参数,更新过程如下:
其中γi=1-αiVii,Vii是V的对角元素。
训练相关向量机模型后,需要计算预测值和真实值的RMSE,作为IMBES的适应度函数。通过计算适应度值,得到核函数参数的个体最优解和全局最优解。对于用IMBES优化核函数参数η,使用相关向量机模型的均方根误差作为适应度函数。
其中yi表示真实值,表示预测值。
在步骤二中,利用步骤一估算完ω和ξi后,用IMBES算法优化核参数η。IMBES算法的过程如下:
秃鹰搜索优化算法是由Alsattar等在2020年提出的。主要模拟秃鹰在寻找猎物时的狩猎策略和智能社会行为。该算法可分为三个部分,即选择空间、搜索空间和俯冲。
(1)帐篷混乱策略。采用帐篷混沌策略对猎物初始位置进行优化,采用线性递减法对秃鹰迭代更新位置的控制参数进行改进。从而找到最优的模型参数,提高拟合质量。帐篷混沌映射函数描述为:
其中λ∈(0,1)。
(2)选择空间。秃鹰随机选择搜索区域,并通过判断猎物的数量来确定最佳位置。在这个阶段,秃鹰的位置Pi,new更新是通过随机搜索先验信息来决定。位置更新的数学模型描述为:
Pi,new=Pbest+C1·rand·(Pmean-Pi);
其中C1是控制位置变化的参数,rand是(0,1)之间的随机数,Pbest是当前最优位置。Pmean是经过前一次搜索后秃鹰位置的平均分布。Pi是第i只秃鹰的位置。
(3)搜索空间。秃鹰在选定的搜索空间内以螺旋模式搜索猎物,加快搜索过程,以找到最佳猎物。螺旋飞行的数学模型使用极性方程更新位置,如下:
θ(i)=π·ω·rand;
r(i)=θ(i)+C2·rand;
xr(i)=r(i)·sin[θ(i)];
yr(i)=r(i)·cos[θ(i)];
其中θ(i)和r(i)分别是螺旋方程的极角和直径。ω和C2是控制螺旋轨迹的参数。x(i)和y(i)分别为极坐标下秃鹰位置,取值为(-1,1)。秃鹰位置更新如下所示:
Pi,new=Pi+y(i)·(Pi-Pi+1)+x(i)·(Pi-Pmean);
其中Pi+1是第i只秃鹰的下一次更新位置。
(4)俯冲。秃鹰迅速从搜索空间的最佳位置俯冲向目标猎物,其余种群同时移动到最佳位置和攻击猎物。运动状态仍由极坐标方程描述如下:
θ(i)=π·ω·rand;
r(i)=θ(i);
xr(i)=r(i)·sinh[θ(i)];
yr(i)=r(i)·cosh[θ(i)];
秃鹰俯冲时位置更新公式为:
Pi,new=rand·Pbest+x1(i)·(Pi-C3·Pmean)+y1(i)·(Pi-C4·Pbest);
其中C3和C4为秃鹰向最佳中心位置移动的强度,其取值范围是(1,2)。
(5)突变策略。为了解决BES算法容易陷入局部优化的问题,在BES算法中引入了变异操作。当秃鹰距离锁定目标猎物的历史最优位置较近时,对满足突变条件的猎物位置进行突变,以增加种群的多样性,增强粒子的全局优化能力。变异操作如下:
Pi,mutation=Pi,new+pc·rand;
pc=pm·(Pi,max-Pi,min);
Pi,mutation变异为变异后的秃鹰位置,pc为变异因子其取值范围是所有猎物位置中最小的定义域。pm代表变化率。Pi,max表示最大集合搜索范围,Pi,min表示最小集合搜索范围,Pi,max-Pi,min表示猎物位置域。
猎物的位置受变化率的影响。在迭代早期,BES算法主要发挥自己的特点,变化率较小。随着迭代次数的增加,算法的多样性越来越差。计算变化率的公式为:
其中,σ为初始变化率集,i为当前迭代次数,N为最大迭代次数。
通过对相关参数进行初始化,通过计算秃鹰的个体适应度,以确定当前最佳和最差的个体。通过更新秃鹰位置,得到新的秃鹰位置和个体适应度值后,将其与之前的值进行比较,分析最佳适应度值。如果本轮较好,则更新全局最优适应度值。当达到最大迭代次数时,优化结束。
利用IMBES优化相关向量机模型得到退化趋势方程,并将其作为PF的状态转换方程。粒子滤波算法分为以下5个步骤:
(1)粒子滤波的状态转换方程和观测方程分别为:
xt=f(xt-1,ut);
zt=h(xt,vt);
其中,xt是系统状态,zt是观测值,ut是过程噪声,vt是观测噪声。
(2)生成粒子。相关向量机可以在退化轨迹上执行回归。在每个时间步,状态的回归值服从高斯分布:
由于相关向量机回归精度高,所以接近于后验分布,对于粒子滤波是一个很好的抽样分布,c表示粒子,c=1,2,...,n。t为状态轨迹的时间步长,t=1,2,...,T。那么粒子抽样的方法如下:
状态转换方程可以通过相关向量机回归得到每个时间步的
(3)更新权重。当获得新的观测值时,权重更新如下:
(4)重采样。按照公式(2-17)进行粒子重采样。
(5)预测。xt的后验分布可计算如下:
其中,如果n足够大,则可以有效地近似真值。
在步骤三中,利用核密度估计进行区间预测。核密度估计是一种非参数估计方法,它不预先假设数据服从某种概率分布,而是从数据本身的特征中获得概率分布。因此,与参数估计相比,例如假设数据服从高斯分布,它能更好地反映数据的真实分布,从而使概率密度函数能达到更好的拟合效果。设E1,E2,…,Eq,…,En为n个预测误差样本,核密度估计形式为:
其中:表示概率分布函数,E表示任意预测误差样本,h表示带宽,K(·)表示满足以下约束的核函数:
ρτ(u)=τuI(u…0)+(τ-1)uI(u<0);
其中I为指标函数,给定一定样本后,概率密度函数主要受核函数K(·)和带宽h的影响。
采用积分法计算累积分布函数,选取累积分布函数的5%分位数,计算95%置信度下预测误差的波动区间。将95%置信误差波动区间与IMBES-RVM-PF模型的预测结果叠加得到轴承寿命的预测区间。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (6)

1.一种滚动轴承剩余使用寿命区间预测方法,其特征在于:所述滚动轴承剩余使用寿命区间预测方法包括以下步骤:
S1:针对秃鹰搜索算法容易陷入局部最优的缺点,采用变异策略对秃鹰搜索算法进行优化形成MBES;为了提高预测精度,采用帐篷混沌策略对MBES算法进行优化;
S2:利用IMBES对相关向量机的核函数参数进行寻优,为了进一步提高点预测的预测精度,利用粒子滤波算法对改进秃鹰搜索算法优化相关向量机预测模型的预测结果进行过滤优化;
S3:为了提高预测可靠性,实际应用中需要进行区间预测;区间预测是指在一定置信区间内预测值范围的估计区间;在点预测结束后,利用核密度估计进行区间预测,即将预测结果利用核函数拟合分布,确定预测上下限;
S4:之后根据步骤S4中得到的上下限可以得到滚动轴承剩余使用寿命区间。
2.如权利要求1所述的一种滚动轴承剩余使用寿命区间预测方法,其特征在于:所述步骤S1中,滚动轴承数据集利用加速度传感器测量,即在滚动轴承的横向和纵向安装两个高灵敏度加速度传感器,分别测量滚动轴承水平方向和垂直方向的振动加速度数值,利用同一型号的滚动轴承水平振动信号进行剩余使用寿命预测。
3.如权利要求1所述的一种滚动轴承剩余使用寿命区间预测方法,其特征在于:所述步骤S2中,利用步骤一估算完ω和ξi后,用IMBES算法优化核参数η。
4.如权利要求1所述的一种滚动轴承剩余使用寿命区间预测方法,其特征在于:所述步骤S3中,利用核密度估计进行区间预测,核密度估计是一种非参数估计方法,它不预先假设数据服从某种概率分布,而是从数据本身的特征中获得概率分布,因此,与参数估计相比,例如假设数据服从高斯分布,它能更好地反映数据的真实分布,从而使概率密度函数能达到更好的拟合效果。
5.如权利要求1所述的一种滚动轴承剩余使用寿命区间预测方法,其特征在于:所述步骤S2中,将IMBES-RVM预测模型代替粒子滤波的状态转换方程进行滤波处理,也进一步解决了粒子滤波状态转换方程依靠机理模型的缺陷,也提高了预测精确度。
6.如权利要求1所述的一种滚动轴承剩余使用寿命区间预测方法,其特征在于:所述步骤S1中,采用改进的秃鹰搜索(IMBES)算法优化相关向量机的核函数参数。
CN202310407813.1A 2023-04-17 一种滚动轴承剩余使用寿命区间预测方法 Active CN116451580B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310407813.1A CN116451580B (zh) 2023-04-17 一种滚动轴承剩余使用寿命区间预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310407813.1A CN116451580B (zh) 2023-04-17 一种滚动轴承剩余使用寿命区间预测方法

Publications (2)

Publication Number Publication Date
CN116451580A true CN116451580A (zh) 2023-07-18
CN116451580B CN116451580B (zh) 2024-10-18

Family

ID=

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115933671A (zh) * 2022-12-17 2023-04-07 智联云大数据科技南京有限公司 基于mbesp优化的agv路径规划与自主避障方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115933671A (zh) * 2022-12-17 2023-04-07 智联云大数据科技南京有限公司 基于mbesp优化的agv路径规划与自主避障方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
石默涵等: "基于混沌优化和自适应反向学习的秃鹰搜索算法", 数学的实践与认识, vol. 52, no. 6, pages 149 - 159 *
贾鹤鸣等: "基于改进秃鹰搜索算法的同步优化特征选择", 控制与决策, vol. 37, no. 2, pages 445 - 454 *
黄维等: "基于改进秃鹰搜索算法的工业机器人几何参数标定方法", 机械设计与研究, vol. 38, no. 6, pages 21 - 33 *

Similar Documents

Publication Publication Date Title
CN107092582B (zh) 一种基于残差后验的异常值在线检测及置信度评估方法
CN110363349B (zh) 一种基于ascs的lstm神经网络水文预测方法及系统
CN111767517B (zh) 一种应用于洪水预测的BiGRU多步预测方法、系统及存储介质
CN109766745B (zh) 强化学习三态组合长短时记忆神经网络系统及训练和预测方法
Yang et al. A new hybrid optimization prediction model for PM2. 5 concentration considering other air pollutants and meteorological conditions
CN109325616A (zh) 一种基于高斯过程回归和萤火虫算法相结合的细颗粒物预测与溯源方法
CN109460608B (zh) 一种基于模糊时间序列的高陡边坡形变预测的方法
CN108832663A (zh) 微电网光伏发电系统的发电功率的预测方法及设备
CN108961460B (zh) 基于稀疏esgp与多目标优化的故障预测方法及装置
Dou et al. Top-down indoor localization with Wi-fi fingerprints using deep Q-network
CN114723188A (zh) 水质预测方法、装置、计算机设备和存储介质
CN113255963A (zh) 基于路元拆分和深度学习模型lstm的路面使用性能预测方法
CN114444561A (zh) 基于CNNs-GRU融合深度学习模型的PM2.5预测方法
CN115629160A (zh) 一种基于时空图的空气污染物浓度预测方法及系统
CN117349797A (zh) 基于人工智能的飞机故障检测方法及系统
CN115938104B (zh) 一种动态短时路网交通状态预测模型及预测方法
CN118261056A (zh) Rbf神经网络驱动的液态铅铋合金的管道冲蚀寿命预测方法及系统
Pourabdollah et al. A new dynamic approach for non-singleton fuzzification in noisy time-series prediction
CN116451580B (zh) 一种滚动轴承剩余使用寿命区间预测方法
CN116960962A (zh) 一种跨区域数据融合的中长期区域负荷预测方法
CN113191526A (zh) 一种基于随机敏感度的短期风速区间多目标优化预测方法及系统
CN116451580A (zh) 一种滚动轴承剩余使用寿命区间预测方法
CN111832818B (zh) 一种基于相关性分析的lstm网络发电多步预测方法
CN113987904A (zh) 一种输电项目修理成本测算方法、装置、设备及存储介质
CN115879569B (zh) 一种IoT观测数据的在线学习方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20240423

Address after: 518000 1002, Building A, Zhiyun Industrial Park, No. 13, Huaxing Road, Henglang Community, Longhua District, Shenzhen, Guangdong Province

Applicant after: Shenzhen Wanzhida Technology Co.,Ltd.

Country or region after: China

Address before: 100144 Beijing City, Shijingshan District Jin Yuan Zhuang Road No. 5

Applicant before: NORTH CHINA University OF TECHNOLOGY

Country or region before: China

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20240909

Address after: 239300 No.34 Tongxin Road, Tianchang City, Chuzhou City, Anhui Province

Applicant after: TIANCHANG CHAJIAN BUILDING INSTALLATION DECORATION CO.,LTD.

Country or region after: China

Address before: 518000 1002, Building A, Zhiyun Industrial Park, No. 13, Huaxing Road, Henglang Community, Longhua District, Shenzhen, Guangdong Province

Applicant before: Shenzhen Wanzhida Technology Co.,Ltd.

Country or region before: China