发明名称 一种考虑机器多维特征属性的作业车间瓶颈识别方法
摘要 本发明提出了一种考虑机器多维特征属性的作业车间瓶颈识别方法。首先采用优化算法进行作业调度优化求解,其次改变传统依据机器单方面因素识别瓶颈的做法,综合考虑瓶颈的多方面特征属性,提出多属性瓶颈识别方法,综合全面评价机器性能进行瓶颈识别。本发明将瓶颈利用和瓶颈识别问题放在统一框架下进行集成求解,解决了传统瓶颈识别与优化方案相割离而导致瓶颈识别不准、优化方案不优等不足;得到调度优化方案后,综合考虑机器的不同特征属性,采用TOPSIS多属性瓶颈识别方法进行瓶颈识别,克服了传统瓶颈识别指标片面、识别结果偏颇的缺陷。
申请公布号 CN102768737B 申请公布日期 2015.10.28
申请号 CN201210232168.6 申请日期 2012.07.06
申请人 西北工业大学 发明人 王军强;陈剑;康永;张仲田;崔福东;张映锋;孙树栋
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 西北工业大学专利中心 61204 代理人 陈星
主权项 一种考虑机器多维特征属性的作业车间瓶颈识别方法,其特征在于:包括以下步骤:步骤1:采用优化算法确定作业车间的最优调度方案Ω:Ω={{B<sub>11</sub>,…,B<sub>v1</sub>,…,B<sub>e1</sub>;C<sub>11</sub>,…,C<sub>v1</sub>,…,C<sub>e1</sub>;Z<sub>11</sub>,…,Z<sub>v1</sub>,…,Z<sub>e1</sub>},…,{B<sub>1i</sub>,…,B<sub>vi</sub>,…,B<sub>ei</sub>;C<sub>li</sub>,…,C<sub>vi</sub>,…,C<sub>ei</sub>;Z<sub>li</sub>,…,Z<sub>vi</sub>,…,Z<sub>ei</sub>},…,{B<sub>1m</sub>,…,B<sub>vm</sub>,…,B<sub>em</sub>;C<sub>lm</sub>,…,C<sub>vm</sub>,…,C<sub>em</sub>;Z<sub>lm</sub>,…,Z<sub>vm</sub>,…,Z<sub>em</sub>}}其中,B<sub>vi</sub>表示第v∈E={1,2,…,e}个工件在第i∈M={1,2,...,m}台机器上的开始加工时间,C<sub>vi</sub>表示第v个工件在第i台机器上的加工完成时间,Z<sub>vi</sub>表示第v个工件在第i台机器上的工装准备时间;步骤2:建立机器集为A={A<sub>1</sub>,A<sub>2</sub>,...,A<sub>m</sub>},评价属性集为X={X<sub>1</sub>,X<sub>2</sub>,...,X<sub>n</sub>};根据步骤1得到的最优调度方案Ω,计算每个机器各评价属性值,得到用于瓶颈识别的决策矩阵D=(x<sub>ij</sub>)<sub>m×n</sub>,x<sub>ij</sub>是第i个的机器的第j个评价属性值;步骤3:采用向量标准化方法将步骤2的决策矩阵D=(x<sub>ij</sub>)<sub>m×n</sub>转换成标准决策矩阵R=(r<sub>ij</sub>)<sub>m×n</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>r</mi><mi>ij</mi></msub><mo>=</mo><mmultiscripts><mrow><msub><mi>x</mi><mi>ij</mi></msub><mo>/</mo><msqrt><mmultiscripts><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>x</mi><mi>ij</mi><mn>2</mn></msubsup></mrow></mmultiscripts></msqrt></mrow></mmultiscripts></mrow>]]></math><img file="FDA00001857425600011.GIF" wi="319" he="140" /></maths>步骤4:采用熵权法确定评价属性权重,评价属性权重向量W=(ω<sub>1</sub>,ω<sub>2</sub>,...,ω<sub>j</sub>,...,ω<sub>n</sub>):步骤4.1:将各评价属性值标准化,其中第j个属性下第i台机器的无量纲化属性值p<sub>ij</sub>为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>p</mi><mi>ij</mi></msub><mo>=</mo><msub><mi>x</mi><mi>ij</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>x</mi><mi>ij</mi></msub></mrow>]]></math><img file="FDA00001857425600012.GIF" wi="286" he="124" /></maths>步骤4.2:计算第j个评价属性的熵值E<sub>j</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>E</mi><mi>j</mi></msub><mo>=</mo><mo>-</mo><mi>&delta;</mi><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>p</mi><mi>ij</mi></msub><mi>ln</mi><msub><mi>p</mi><mi>ij</mi></msub></mrow>]]></math><img file="FDA00001857425600013.GIF" wi="376" he="113" /></maths>其中δ=1/lnm,保证了0≤E<sub>j</sub>≤1;步骤4.3:计算第j个属性的信息偏差度d<sub>j</sub>:d<sub>j</sub>=1‑E<sub>j</sub>步骤4.4:对信息偏差度进行归一化计算权重:当决策者没有评价属性间的偏好时:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>j</mi></msub><mo>=</mo><msub><mi>d</mi><mi>j</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>d</mi><mi>j</mi></msub></mrow>]]></math><img file="FDA00001857425600021.GIF" wi="286" he="131" /></maths>当决策者有评价属性间的偏好时:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>j</mi></msub><mo>=</mo><msub><mi>&lambda;</mi><mi>j</mi></msub><msub><mi>d</mi><mi>j</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&lambda;</mi><mi>j</mi></msub><msub><mi>d</mi><mi>j</mi></msub></mrow>]]></math><img file="FDA00001857425600022.GIF" wi="374" he="131" /></maths>其中λ<sub>j</sub>为决策者的主观权重;步骤5:将步骤4得到的评价属性权重纳入标准决策矩阵R=(r<sub>ij</sub>)<sub>m×n</sub>中,得到加权标准化决策矩阵V=(v<sub>ij</sub>)<sub>m×n</sub>=(ω<sub>j</sub>r<sub>ij</sub>)<sub>m×n</sub>;步骤6:确定评价属性的正、负理想解A<sup>+</sup>和A<sup>‑</sup>:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msup><mi>A</mi><mo>+</mo></msup><mo>=</mo><mo>{</mo><mrow><mo>(</mo><munder><mi>max</mi><mi>i</mi></munder><msub><mi>v</mi><mi>ij</mi></msub><mo>|</mo><mi>j</mi><mo>&Element;</mo><mi>J</mi><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><munder><mi>min</mi><mi>i</mi></munder><msub><mi>v</mi><mi>ij</mi></msub><mo>|</mo><mi>j</mi><mo>&Element;</mo><msup><mi>J</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>|</mo><mi>i</mi><mo>=</mo><mi>M</mi><mo>}</mo><mo>=</mo><mo>{</mo><msubsup><mi>v</mi><mn>1</mn><mo>+</mo></msubsup><mo>,</mo><msubsup><mi>v</mi><mn>2</mn><mo>+</mo></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>v</mi><mi>j</mi><mo>+</mo></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>v</mi><mi>n</mi><mo>+</mo></msubsup><mo>}</mo></mrow>]]></math><img file="FDA00001857425600023.GIF" wi="1343" he="69" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msup><mi>A</mi><mo>-</mo></msup><mo>=</mo><mo>{</mo><mrow><mo>(</mo><munder><mi>min</mi><mi>i</mi></munder><msub><mi>v</mi><mi>ij</mi></msub><mo>|</mo><mi>j</mi><mo>&Element;</mo><mi>J</mi><mo>)</mo></mrow><mo>,</mo><mrow><mo>(</mo><munder><mi>max</mi><mi>i</mi></munder><msub><mi>v</mi><mi>ij</mi></msub><mo>|</mo><mi>j</mi><mo>&Element;</mo><msup><mi>J</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>|</mo><mi>i</mi><mo>&Element;</mo><mi>M</mi><mo>}</mo><mo>=</mo><mo>{</mo><msubsup><mi>v</mi><mn>1</mn><mo>-</mo></msubsup><mo>,</mo><msubsup><mi>v</mi><mn>2</mn><mo>-</mo></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>v</mi><mi>j</mi><mo>-</mo></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>v</mi><mi>n</mi><mo>-</mo></msubsup><mo>}</mo></mrow>]]></math><img file="FDA00001857425600024.GIF" wi="1341" he="67" /></maths>其中J为效益型评价属性集合,J′为成本型评价属性集合;步骤7:通过n维Euclid距离计算各机器与正理想解以及负理想解之间的距离:每台机器到正理想解的距离为:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>S</mi><msup><mi>i</mi><mo>+</mo></msup></msub><mo>=</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>v</mi><mi>ij</mi></msub><mo>-</mo><msubsup><mi>v</mi><mi>j</mi><mo>+</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow>]]></math><img file="FDA00001857425600025.GIF" wi="394" he="146" /></maths>每台机器到负理想解的距离为:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>S</mi><msup><mi>i</mi><mo>-</mo></msup></msub><mo>=</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>v</mi><mi>ij</mi></msub><mo>-</mo><msubsup><mi>v</mi><mi>j</mi><mo>-</mo></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mrow>]]></math><img file="FDA00001857425600026.GIF" wi="393" he="146" /></maths>步骤8:计算每台机器与正理想解A<sup>+</sup>的贴近度C<sub>i</sub>:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>C</mi><mi>i</mi></msub><mo>=</mo><msub><mi>S</mi><msup><mi>i</mi><mo>-</mo></msup></msub><mo>/</mo><mrow><mo>(</mo><msub><mi>S</mi><msup><mi>i</mi><mo>+</mo></msup></msub><mo>+</mo><msub><mi>S</mi><msup><mi>i</mi><mo>-</mo></msup></msub><mo>)</mo></mrow><mo>,</mo><mn>0</mn><mo>&lt;</mo><msub><mi>C</mi><mi>i</mi></msub><mo>&lt;</mo><mn>1</mn></mrow>]]></math><img file="FDA00001857425600027.GIF" wi="625" he="61" /></maths>得到贴近度C<sub>i</sub>最大的机器为作业车间的瓶颈。
地址 710072 陕西省西安市友谊西路127号