发明名称 一种基于任务可靠度的自组织配作加工方法
摘要 本发明提供了一种基于任务可靠度的自组织配作加工方法,它利用控制图原理建立配作加工过程的卡方控制区域,并给出对应配作过程的失效阈值判定标准。通过历史数据构建匹配特征值到失效阈值的距离与任务可靠度的关系模型,进而利用该模型评估任意配作加工过程在给定置信水平下的任务可靠度,同时可估计待匹配样本匹配成功的概率。最后,对比生产实际的任务可靠度要求,确定配作加工方案;它有五大步骤。本发明能精确评估样本和总体的可靠度,并给出一套指导配作加工过程的工具表格,并能根据实际配作加工数据不断更新修正模型,为改进配作加工过程提供了新的依据和方法。
申请公布号 CN104090494A 申请公布日期 2014.10.08
申请号 CN201410312351.6 申请日期 2014.07.02
申请人 北京航空航天大学 发明人 戴伟;匡芬;赵宇
分类号 G05B17/00(2006.01)I 主分类号 G05B17/00(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种基于任务可靠度的自组织配作加工方法,其特征在于:它包括如下具体步骤:步骤一:对零件i对应零件j的第k个质量特性的n个历史数据<img file="FDA0000531835040000011.GIF" wi="826" he="86" />进行正态性检验,若该特性为正态随机变量,即Minitab正态性检验结果显示p值大于0.05则转步骤二;步骤二:根据多变量卡方控制图理论,建立配作过程的卡方控制超椭球区域,如下式所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>&chi;</mi><mn>0</mn><mn>2</mn></msubsup><mo>=</mo><mi>n</mi><msup><mrow><mo>(</mo><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>&prime;</mo></msup><msup><mi>&Sigma;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000531835040000012.GIF" wi="1307" he="82" /></maths>上式中,<img file="FDA0000531835040000013.GIF" wi="259" he="79" />p为关联质量特性个数;特别地,当p=2时,有<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>&chi;</mi><mn>0</mn><mn>2</mn></msubsup><mo>=</mo><mfrac><mi>n</mi><mrow><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>&sigma;</mi><mn>12</mn><mn>2</mn></msubsup></mrow></mfrac><mo>[</mo><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup><msup><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><msup><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mn>2</mn><msub><mi>&sigma;</mi><mn>12</mn></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000531835040000014.GIF" wi="1746" he="131" /></maths>此时,<img file="FDA0000531835040000015.GIF" wi="246" he="84" />由于配作过程的关联质量特性不独立,σ<sub>12</sub>≠0,将上式转化为标准型,令<img file="FDA0000531835040000016.GIF" wi="524" he="78" />椭圆旋转角度是θ,则有:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>x</mi><mn>1</mn><mo>&prime;</mo></msubsup><mo>=</mo><msubsup><mi>x</mi><mn>1</mn><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mi>cos</mi><mi>&theta;</mi><mo>-</mo><msubsup><mi>x</mi><mn>2</mn><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mi>sin</mi><mi>&theta;</mi></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mn>2</mn><mo>&prime;</mo></msubsup><mo>=</mo><msubsup><mi>x</mi><mn>1</mn><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mi>sin</mi><mi>&theta;</mi><mo>+</mo><msubsup><mi>x</mi><mn>2</mn><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mi>cos</mi><mi>&theta;</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000531835040000017.GIF" wi="1260" he="171" /></maths>代入得<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfrac><msubsup><mi>x</mi><mn>1</mn><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mfrac><mrow><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>&sigma;</mi><mn>12</mn><mn>2</mn></msubsup><mo>)</mo></mrow><msubsup><mi>&chi;</mi><mrow><mi>&alpha;</mi><mo>,</mo><mn>2</mn></mrow><mn>2</mn></msubsup></mrow><mrow><mi>n</mi><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup><msup><mi>cos</mi><mn>2</mn></msup><mi>&theta;</mi><mo>+</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><msup><mi>sin</mi><mn>2</mn></msup><mo>-</mo><msub><mi>&sigma;</mi><mn>12</mn></msub><mi>sin</mi><mn>2</mn><mi>&theta;</mi><mo>)</mo></mrow></mrow></mfrac></mfrac><mo>+</mo><mfrac><msubsup><mi>x</mi><mn>2</mn><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mfrac><mrow><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>&sigma;</mi><mn>12</mn><mn>2</mn></msubsup><mo>)</mo></mrow><msubsup><mi>&chi;</mi><mrow><mi>&alpha;</mi><mo>,</mo><mn>2</mn></mrow><mn>2</mn></msubsup></mrow><mrow><mi>n</mi><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup><msup><mi>sin</mi><mn>2</mn></msup><mi>&theta;</mi><mo>+</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><msup><mi>cos</mi><mn>2</mn></msup><mi>&theta;</mi><mo>+</mo><msub><mi>&sigma;</mi><mn>12</mn></msub><mi>sin</mi><mn>2</mn><mi>&theta;</mi><mo>)</mo></mrow></mrow></mfrac></mfrac><mo>=</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000531835040000018.GIF" wi="1801" he="221" /></maths>在化简过程中,令<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><mi>sin</mi><mn>2</mn><mi>&theta;</mi><mo>-</mo><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup><mi>sin</mi><mn>2</mn><mi>&theta;</mi><mo>-</mo><mn>2</mn><msub><mi>&sigma;</mi><mn>12</mn></msub><mrow><mo>(</mo><msup><mi>cos</mi><mn>2</mn></msup><mi>&theta;</mi><mo>-</mo><msup><mi>sin</mi><mn>2</mn></msup><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>,</mo></mrow>]]></math><img file="FDA0000531835040000019.GIF" wi="975" he="86" /></maths>得到旋转角的计算关系式,如式所示:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>tan</mi><mn>2</mn><mi>&theta;</mi><mo>=</mo><mfrac><mrow><mn>2</mn><msub><mi>&sigma;</mi><mn>12</mn></msub></mrow><mrow><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup></mrow></mfrac><mo>,</mo><mi>tan</mi><mn>2</mn><mi>&theta;</mi><mo>&NotEqual;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00005318350400000110.GIF" wi="1345" he="141" /></maths>故建立控制区域即椭圆区域;步骤三:构建距离可靠度模型,通过统计拟合,确定控制区域到约束匹配公差区域的距离与配作过程可靠度之间的解析关系,给出解析表达式;其中约束公差区域匹配前后对各特性公差要求的交集,为六边形区域AEFCGH,在匹配前待匹配零部件的关联质量特性应落在矩形区域ABCD范围内,匹配目标为匹配后关联特性的观测值落在两条平行线EF、GH之间;根据工程实践经验,当输入样本特性分别满足各自的公差要求时,配作过程的任务可靠度随着控制区域到约束匹配公差区域距离的减小而减小,因此用递增函数来拟合可靠度R与距离d之间的定量关系为:R=1‑λe<sup>‑d/η</sup>   (6)其中,λ,η均为待估参数,根据工程经验或统计拟合得到;步骤四:对于某匹配样本,确定匹配前的质量特性观测值,计算观测值落在哪个超椭球范围内确定该样本的配作过程可靠度和置信度;步骤五:给出具体的配作加工方案;给定某一零件的待匹配特性值,过该点作垂直于坐标轴的直线与控制域相交,则匹配零件特性落在相交区域范围内则配作过程可靠;当配作过程只涉及一个关联特性时,若特性1取值为x,则选取特性2的值在x~y<sub>U</sub>的零件来匹配x;同理若特性2取值为y,则选取特性1的值在x<sub>L</sub>~y的零件来匹配y。
地址 100191 北京市海淀区学院路37号