[发明专利]一种可重入制造系统瓶颈设备预测方法

摘要

本发明公开一种基于ANFIS(自适应模糊推理系统)的可重入制造系统瓶颈设备预测方法。首先确定影响系统瓶颈设备的参数及ANFIS输入与输出；其次采用固定工件移动量间隔的方法采集原型系统数据并采用数据比对方法对数据进行预处理得到ANFIS的输入输出数据集；最后用训练数据确定ANFIS参数，用测试数据进行预测得到下一时刻瓶颈设备编号，设备利用率及队列长度并对预测结果进行分析得出预测精度，并且每次进行瓶颈预测之前，采用模式匹配方法对ANFIS模型参数在线调整。利用本发明能够构建可重入制造系统瓶颈设备预测方案，准确地预测出下一时刻系统的瓶颈设备及相关参数，并且对后续基于瓶颈设备的调度方法有指导意义。

主权项

1.一种可重入制造系统瓶颈设备预测方法，其特征在于，包含以下步骤：(1)确定影响系统中瓶颈设备的关键因素为工件类型、工件投料方式、设备加工时间、设备故障间隔时间、设备平均维修时间；以及设备利用率、缓存区队列长度、生产线在制品水平WIP；(2)获取以下实际生产线数据：(2.1.1)每台设备的设备加工时间、设备故障间隔时间、设备平均维修时间；(2.1.2)产品工艺路径：总共O道工序；(2.1.3)工件类型、工件投料方式；根据上述数据信息构建模型，模型由M台设备及工件投料装置、工件回收装置构成，投料装置按照某种投料方式进行投料，投放W种类型工件进入生产线中，工件按照工艺路径在设备中根据加工时间进行加工，最后由回收装置收集成品；对该模型进行仿真，设工件在任意一台设备上完成一道工序的加工为一步，每完成一步进行一次数据采集，得到R₀组数据样本：{QL_s(1)，QL_s(2)，L，QL_s(t)，L，QL_s(R₀)}，{η_s(1)，η_s(2)，Lη_s(t)L，η_s(R₀)}，{W(1)，W(2)，L，W(t)，L，W(R₀)}，{Q(1)，Q(2)，L，Q(t)，L，Q(R₀)}，{T(1)，T(2)，L，T(t)，L，T(R₀)}，t＝1，2，L，R₀，s为设备编号，s＝1，2，L，M；其中，QL_s(t)、η_s(t)分别为第t步数据采集所得到的设备s的缓冲区队列长度与设备利用率，W(t)为第t步数据采集得到的生产线WIP，Q(t)、T(t)分别为第t步数据采集得到的生产线的投料方式及产品类型；除({Q(1)，Q(2)，L，Q(t)，L，Q(R₀)}，{T(1)，T(2)，L，T(t)，L，T(R₀)}外，对每一个数据样本进行如下数据修正：设采集得到的数据样本为{g(t)}，t＝1，2，L，R₀，且实际生产线中对应数据的样本均值为ξ，方差为σ，定义数据偏差为D(t)＝g(t)-ξ，若|D(t)|≤2σ，t＝1，2，L，R₀，则认为g(t)能够反映实际生产线状况，否则剔除第t步所采集的所有数据，得到修正后的数据；(3)对上述修正过的数据按下式计算每个设备的瓶颈度，并进行数据再修正：BVs(n)=θ1ηs(n)+θ2QLs(n)T---(1)]]>式中QL_s(n)，η_s(n)，分别为第n步采集得到的设备s的缓冲区队列长度及设备利用率，BV_s(n)第n步计算得到的设备s瓶颈度，s为设备编号，n＝1，2，L，R₀，但其中不包括步骤(2)剔除的数据；T为缓冲数预估上限，θ₁与θ₂为经验参数，θ₁，θ₂∈[0，1]，θ₁+θ₂＝1；设瓶颈度阈值B∈[0.8，0.95]，若第n步每台设备的瓶颈度BV_s(n)均小于B，则剔除第n步所采集数据，得到再修正的数据样本；若存在BV_s(n)≥B，取瓶颈度最大值所对应的设备作为系统中瓶颈设备，得到第n步瓶颈设备编号及相应的设备利用率、缓冲区队列长度数据；(4)确定每台设备的队列长度、设备利用率、加工时间、故障时间间隔、平均维修时间，以及生产线中WIP，投料方式，工件类型为三个ANFIS的输入，下一时刻瓶颈设备编号、缓冲区队列长度、设备利用率分别为三个ANFIS的输出；(5)采样数据经过数据预处理后，得到H+1组数据，按如下方式构建三个数据集[x1j,x2j,x3jL,x5×M+1j,patternj,typej,y1j+1],]]>[x1j,x2j,x3jL,x5×M+1j,patternj,typej,y2j+1],]]>[x1j,x2j,x3jL,x5×M+1j,patternj,typej,y3j+1],]]>及确定需要三个ANFIS；式中j代表第j组数据，j＝1，2，L，H；分别为第j组数据中设备1的队列长度、设备利用率、加工时间、故障时间间隔、平均维修时间；分别为第j组数据中设备2的队列长度、设备利用率、加工时间、故障时间间隔、平均维修时间；以此类推；第j组数据中生产线WIP，M为设备数量；第j+1组数据中瓶颈设备编号；第j+1组数据中瓶颈设备缓冲区队列长度；第j+1组数据中瓶颈设备的设备利用率；pattern^j：第j组数据中的投料方式；type^j：第j组数据中的工件类型；以上3个数据集均由H组数据组成，取每个数据集V组为训练数据，V＝[0.6H，0.7H]，剩余H-V组为测试数据，按如下步骤对数值型数据进行聚类，确定3个ANFIS结构，并对类别型数据进行编码，ANFIS结构确定包括每个输入变量对应隶属度函数及规则的个数；(5.1)对每个数据集中的类别型数据按如下方式进行编码，形成一个行向量：(5.1.1)设投料方式有A种，对第j组数据，令若第j组数据中生产线采用某种投料方式，则对应的行向量中的元素为1，其余均为0；(5.1.2)设工件类型有C种，对第j组数据，令第j组数据中生产线存在几种类型工件，则对应的行向量中的元素为1，其余均为0；(5.1.3)将pattern^J，type^j  行向量构成向量s^j，即sj=[patternj,typej]=[p1j,p2j,L,pAj,q1j,q2j,L,qCj],]]>则完成对类别数据的编码；(5.2)对每个数据集的训练数据中数值型数据采用二叉树及模糊C均值进行聚类，并根据聚类后子空间形成初始模糊推理系统，步骤如下：(5.2.1)令二叉树的根节点为LN(1)，LN(1)为训练数据中的数值型数据，LN(1)＝[X，Y1]，X＝[x¹L x^jL x^V]^T，X为由x^j构成的矩阵，x^j为由构成的行向量，为数据集中的第j组数据，分别代表每台设备的队列长度、设备利用率、加工时间、故障时间间隔、平均维修时间以及WIP，Y₁为由构成的列向量，j＝1，2，L，V，为第j+1组数据中的瓶颈设备编号；初始化叶节点集合为Ψ＝{LN(M)}，整棵树的节点M＝1；(5.2.2)设定误差限R，R∈[0.001，0.05]，对叶节点集合Ψ中所有节点采用二叉树法进行判定，集合中每个节点是否需要被继续划分，若节点M的线性拟合误差e(M)大于R，则按照步骤(5.2.3)继续划分，若所有节点的线性拟合误差均小于R，则转到步骤(5.2.4)，；(5.2.3)采用模糊C均值对e(M)大于误差限R的节点LN(M)进行聚类，详细聚类步骤如下：(5.2.3.1)给定聚类中心数为c＝2，设定容许误差E_max及加权指数m＝2，E_max∈[0.001，0.05]，随机初始化聚类中心；(5.2.3.2)若模糊C均值的聚类准则函数值J_m≤E_max，则结束聚类；(5.2.3.3)节点LN(M)被划分成两个叶节点LN(M+1)，LN(M+2)；(5.2.3.4)更新节点集合：将Ψ与{LN(M+1)，LN(M+2)}合并，更新Ψ；(5.2.3.5)返回第(5.2.2)步；(5.2.4)得到叶节点数据集Ψ＝{LN(1)，LN(2)，L，LN(K₁)}，LN(1)，LN(2)，L，LN(K₁)为聚类后形成的数据子空间，K₁为通过上述聚类过程得到的聚类中心的个数；(5.2.5)根据聚类中心个数确定模糊系统的模糊规则数为K₁；将(5.2.1)中的Y1=[y12Ly1j+1Ly1V+1]T,]]>分别替换成Y2=[y22Ly2j+1Ly2V+1]T,]]>Y3=[y32Ly3j+1Ly3V+1]T,]]>Y₁为由构成的列向量，Y₃为由构成的列向量，j＝1，2，L，V，为第j+1组数据中的瓶颈设备利用率，为第j+1组数据中的瓶颈设备缓冲区队列长度，重复上述聚类过程，得到每个ANFIS的中模糊规则数K₁，K₂，K₃；(6)按照如下步骤确定三个ANFIS中的参数：(6.1)初始化参数集其中称为前件参数，i＝1，2，L，5×M+1，k＝1，2，L K₁，K₁为模糊规则数量，v_i，k，σ_i，k均为(0，1)中的随机数；为(A+C)×K₁维矩阵，矩阵中每个元素均为(0，1)中的随机数；称为后件参数，k＝1，2，L K₁，均为(0，1)中的随机数；(6.2)第一个数据集中训练数据输入到ANFIS中，确定第一个ANFIS中参数：(6.2.1)对数值型数据从ANFIS第一层输入，第一层有5×M+1个输入量，输入第j组数据为j＝1，2，L，V，为数据集中的第j组数据，分别代表每台设备的队列长度、设备利用率、加工时间、故障时间间隔、平均维修时间以及WIP，第一层对输入变量模糊化，输出为对应模糊集的隶属度，这里隶属度函数采用高斯型函数：O1,hj=μA,ik(xij)=exp(-12(xij-vi,kσi,k)2)---(2)]]>其中i＝1，2，L，5×M+1；k＝1，2，L K₁；h＝(5×M+1)×K₁，K₁为模糊规则数量，为隶属度函数，代表第一层输出，其中为前件参数；(6.2.2)第二层，采用乘法规则计算每条模糊规则的激活强度w_k，第二层各节点的输出为k＝1，2L，K₁；O2,kj=wk=Πi=126μA,ik(xij)---(3)]]>(6.2.3)第三层，归一化得到各条模糊规则的激活强度输出为k＝1，2L，K₁；O3,k′j=wk‾=wkΣk=1K1wk---(4)]]>对编码后的类别输入经过一个激活强度转换矩阵T₁，且输入到第三层F^j＝s^jT₁F^j＝s^jT₁                 (5)其中j＝1，2，L，V，sj=[patternj,typej]=[p1j,p2j,L,pAj,q1j,q2j,L,qCj]]]>为1×(A+C)行向量，Fj=[f1j,f2j,L,fK1j]]]>为K₁维行向量，T₁为(A+C)×K₁维矩阵；则第三层的总输出为；k＝1，2L，K₁，j＝1，2，L，V；(6.2.4)第四层，每个节点的传递函数为线性函数，计算得到每条规则的输出O_4，kO4,kj=O3,kj(d0k+d1kx1j+d2kx2j+L+d5×M+1kx5×M+1j)=O3,kj·gk(xj)---(6)]]>式中gk(xj)=(d0k+d1kx1j+d2kx2j+L+d5×M+1kx5×M+1j),]]>P1c={d0k,d1k,d2k,L,d5×M+1k}]]>为后件参数，xj=[x1j,L,x5×M+1j],]]>，k＝1，2，L，K₁，j＝1，2，L，V；(6.2.5)计算所有规则输出之和O5j=Σk=1K1O4,kj---(7)]]>(6.3)设定最大训练次数epoc及误差限errot，epoc∈[100，500]，error∈[0，0.005]，训练数据的总输出误差为：E0=Σj=1Jet=Σj=1J12[O5j-y1j+1]2---(8)]]>式中为Y₁中第j+1组数据中的瓶颈设备编号，J＝V，Y₁为第一个ANFIS中训练数据中的输出数据，为第j组训练数据输入到ANFIS中得到的输出；若E₀≤error，或者达到最大训练次数epoc，结束训练并转到(6.5)，否则进行(6.4)(6.4)训练参数集(6.4.1)固定参数集采用误差反向传播算法对进行修正；(6.4.2)固定参数集采用最小二乘法对进行修正；(6.4.3)固定参数集采用最小二乘法对进行修正；(6.4.4)得到修正后的参数集返回步骤(6.2)；(6.5)得到确定参数集用第二个数据集的训练数据对第二个ANFIS进行训练，重复上述训练过程，保持ANFIS的输入数据不变，输出数据将替换Y₁为Y₂，J＝V，为Y₂中第j+1组数据中的瓶颈设备利用率，K₁将替换为K₂，确定第二个ANFIS的模糊规则数为K₂，训练ANFIS，确定第二个ANFIS的参数集用第三个数据集的训练数据对第三个ANFIS进行训练，重复上述训练过程，保持ANFIS的输入数据不变，输出数据将替换Y₁为Y₂，J＝V，为Y₃中第j+1组数据中的瓶颈设备缓冲区队列长度，K₁将替换为K₃，确定第二个ANFIS的模糊规则数为K₃，训练ANFIS，确定第三个ANFIS的参数集按照以上步骤用训练数据完成对三个ANFIS的训练，训练完成后；将后H-V组测试数据输入到ANFIS中，得到预测输出。