发明名称 基于自适应遗传算法的机床切削用量能耗优化方法
摘要 本发明公开了一种基于自适应遗传算法的机床切削用量能耗优化方法,包括以下步骤:1)确定模型优化变量步骤;2)确定优化目标函数步骤;3)确定模型中约束条件步骤;4)使用自适应遗传算法来确定切削用量。本发明的优点是:本发明由于采用了自适应的遗传算法方案,所以在切削量选取方面更加的合理,有效的提高了机床的利用效率,减少了能源消耗。
申请公布号 CN105844356A 申请公布日期 2016.08.10
申请号 CN201610170152.5 申请日期 2016.03.22
申请人 江南大学 发明人 王艳;彭竹清;纪志成
分类号 G06Q10/04(2012.01)I;G06N3/12(2006.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 无锡市大为专利商标事务所(普通合伙) 32104 代理人 曹祖良;韩凤
主权项 基于自适应遗传算法的机床切削用量能耗优化方法,其特征是,包括以下步骤:(1)确定模型优化变量:在切削优化数学模型中,切削加工的三个要素:切削速度v<sub>c</sub>、进给量f、背吃刀量a<sub>p</sub>是加工过程中影响碳排放、加工能耗的3个最活跃并且相互独立的变量,在数控编程中a<sub>p</sub>是由用户根据加工余量确定,v<sub>c</sub>由主轴转速n和切削直径来确定,因此,优化变量选取主轴转速n和进给量f两个参数;(2)确定优化目标函数:在离散制造系统加工中,一个零件的加工能耗分为工步层能源消耗、工序层能源消耗、零件层能源消耗、产品层能源消耗,(2.1)工步层的能量消耗函数模型为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>E</mi><mn>1</mn></msub><mo>=</mo><msub><mi>E</mi><mrow><mi>s</mi><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>E</mi><mrow><mi>s</mi><mo>-</mo><mi>s</mi></mrow></msub><mo>+</mo><msub><mi>E</mi><mrow><mi>i</mi><mi>e</mi></mrow></msub><mo>+</mo><msub><mi>E</mi><mi>c</mi></msub><mo>=</mo><mrow><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>t</mi><mi>w</mi></msub></msubsup><mrow><msub><mi>P</mi><mi>w</mi></msub><mi>d</mi><mi>t</mi></mrow></mrow><mo>+</mo><mrow><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>t</mi><mrow><mi>t</mi><mi>e</mi></mrow></msub></msubsup><mrow><msub><mi>P</mi><mrow><mi>i</mi><mi>e</mi></mrow></msub><mi>d</mi><mi>t</mi></mrow></mrow><mo>+</mo><msubsup><mo>&Integral;</mo><mn>0</mn><msub><mi>t</mi><mi>c</mi></msub></msubsup><msub><mi>P</mi><mi>c</mi></msub><mi>d</mi><mi>t</mi></mrow>]]></math><img file="FDA0000947135320000011.GIF" wi="927" he="86" /></maths>其中E<sub>1</sub>表示工步能耗,E<sub>st</sub>表示机床的启动能耗,E<sub>s‑s</sub>表示工步的待机能耗,E<sub>ie</sub>表示工步的空载能耗,E<sub>c</sub>表示切削能耗;t<sub>w</sub>为待机时间,P<sub>w</sub>为设备的总功率,t<sub>ie</sub>表示的是空载时间,P<sub>ie</sub>表示的是空载功率,P<sub>c</sub>表示切削功率,t<sub>c</sub>为加工时间,车削加工过程中,P<sub>c</sub>表示为:<maths num="0002"><math><![CDATA[<mrow><msub><mi>P</mi><mi>c</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>6</mn><mo>&times;</mo><msup><mn>10</mn><mn>4</mn></msup></mrow></mfrac><msub><mi>C</mi><mrow><mi>F</mi><mi>C</mi></mrow></msub><msubsup><mi>a</mi><mi>p</mi><msub><mi>x</mi><mrow><mi>F</mi><mi>C</mi></mrow></msub></msubsup><msup><mi>f</mi><msub><mi>y</mi><mrow><mi>F</mi><mi>C</mi></mrow></msub></msup><msubsup><mi>v</mi><mi>C</mi><msub><mi>n</mi><mrow><mi>F</mi><mi>C</mi></mrow></msub></msubsup><msub><mi>K</mi><mrow><mi>F</mi><mi>C</mi></mrow></msub><msub><mi>v</mi><mi>C</mi></msub></mrow>]]></math><img file="FDA0000947135320000012.GIF" wi="836" he="99" /></maths>其中x<sub>fc</sub>、y<sub>fc</sub>、n<sub>fc</sub>分别是背吃刀量a<sub>p</sub>、进给量f、切削速度v<sub>c</sub>的指数,K<sub>FC</sub>表示各种因素对切削力的修正系数,C<sub>FC</sub>表示决定于被加工金属和金属条件的系数;(2.2)工序层的能量消耗:在工序层,每一个工序的能耗包括其各个工步的能耗加上用于零件运输的能耗,<maths num="0003"><math><![CDATA[<mrow><msub><mi>E</mi><mn>2</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>1</mn></msub></munderover><msub><mi>E</mi><msub><mn>1</mn><mi>i</mi></msub></msub><mo>+</mo><msub><mi>E</mi><mrow><mi>t</mi><mi>s</mi></mrow></msub></mrow>]]></math><img file="FDA0000947135320000013.GIF" wi="414" he="103" /></maths>式中:E<sub>2</sub>为工序能耗,N<sub>1</sub>表示工序中工步个数,<img file="FDA0000947135320000017.GIF" wi="55" he="47" />表示每个工步的能耗,E<sub>ts</sub>表示运输能耗,运输能耗E<sub>ts</sub>由下式来计算得到:<maths num="0004"><math><![CDATA[<mrow><msub><mi>E</mi><mrow><mi>t</mi><mi>s</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>P</mi><mrow><mi>t</mi><mi>s</mi></mrow></msub><mo>&times;</mo><msub><mi>t</mi><mrow><mi>t</mi><mi>s</mi></mrow></msub></mrow><mi>Q</mi></mfrac></mrow>]]></math><img file="FDA0000947135320000014.GIF" wi="198" he="96" /></maths>式中:p<sub>ts</sub>为运输设备的额定功率,t<sub>ts</sub>为运输设备的运行时间,Q为一次运输装载的零件数量;(2.3)零件层的能量消耗:在这一层,当零件供应不足时,会出现机床被迫等待的情况,同时,考虑到零件清洗和上漆的耗能,生产出一个零件的能耗表示为<maths num="0005"><math><![CDATA[<mrow><msub><mi>E</mi><mn>3</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>2</mn></msub></munderover><msub><mi>E</mi><msub><mn>2</mn><mi>i</mi></msub></msub><mo>+</mo><msub><mi>E</mi><mrow><mi>p</mi><mi>t</mi></mrow></msub><mo>+</mo><msub><mi>E</mi><mrow><mi>c</mi><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>E</mi><mrow><mi>s</mi><mo>-</mo><mi>p</mi></mrow></msub></mrow>]]></math><img file="FDA0000947135320000015.GIF" wi="492" he="115" /></maths>式中:E<sub>3</sub>为零件能耗,N<sub>2</sub>为加工零件的工序数,E<sub>2i</sub>表示每个工序的能耗,E<sub>pt</sub>为上漆能耗,E<sub>cn</sub>为清洗能耗,E<sub>s‑p</sub>表示机器等待能耗;(2.4)产品层的能量消耗:<maths num="0006"><math><![CDATA[<mrow><msub><mi>E</mi><mn>4</mn></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mn>3</mn></msub></munderover><msub><mi>E</mi><msub><mn>3</mn><mi>i</mi></msub></msub><mo>+</mo><msub><mi>E</mi><mrow><mi>a</mi><mi>e</mi></mrow></msub><mo>+</mo><mfrac><msub><mi>E</mi><mrow><mi>f</mi><mi>y</mi></mrow></msub><mrow><mi>K</mi><mi>M</mi></mrow></mfrac></mrow>]]></math><img file="FDA0000947135320000016.GIF" wi="468" he="127" /></maths>式中:E<sub>4</sub>表示产品的能耗,N<sub>3</sub>表示加工产品零件的个数,E<sub>3i</sub>表示每个零件的能耗,E<sub>ae</sub>表示装配产品的耗能,E<sub>fy</sub>辅助设备的耗能,K和M分别为一定时间内生产产品的种类数和特定种类产品的个数;(3)确定模型中约束条件,包括所选设备的切削速度、最大进给量、表面粗糙度、机床功率、切削力中的一种或多种;(4)使用自适应遗传算法来确定切削用量;(4.1)设置自适应遗传算法参数,包括种群的个数、最大运行的代数、个体染色体的长度、遗传操作的交叉概率p<sub>c</sub>和变异概率p<sub>m</sub>;(4.2)物种的编码:在切削量的优化参数设置好后,用字长为m的二进制串s<sub>m</sub>来表示进给量f,进给量f的取值范围为[f<sub>min</sub>,f<sub>max</sub>],用m位二进制数x<sub>1</sub>表示,其关系为:<maths num="0007"><math><![CDATA[<mrow><mi>f</mi><mo>=</mo><msub><mi>f</mi><mi>min</mi></msub><mo>+</mo><mfrac><msub><mi>x</mi><mn>1</mn></msub><mrow><msup><mn>2</mn><mi>m</mi></msup><mo>-</mo><mn>1</mn></mrow></mfrac><mrow><mo>(</mo><msub><mi>f</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>f</mi><mi>min</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000947135320000021.GIF" wi="597" he="109" /></maths>主轴转速n的取值范围为[n<sub>min</sub>,n<sub>max</sub>],用字长为K的二进制串s<sub>k</sub>来表示,它的二进制数为x<sub>2</sub>,其关系为:<maths num="0008"><math><![CDATA[<mrow><mi>n</mi><mo>=</mo><msub><mi>n</mi><mi>min</mi></msub><mo>+</mo><mfrac><msub><mi>x</mi><mn>2</mn></msub><mrow><msup><mn>2</mn><mi>k</mi></msup><mo>-</mo><mn>1</mn></mrow></mfrac><mrow><mo>(</mo><msub><mi>n</mi><mi>max</mi></msub><mo>-</mo><msub><mi>n</mi><mi>min</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000947135320000022.GIF" wi="553" he="103" /></maths>根据上述2个公式可对主轴转速和进给量进行编码,以进给量和主轴转速为优化变量的个体染色体表示为s<sub>m</sub>s<sub>k</sub>;(4.3)适应值的计算:适应值函数如下:<maths num="0009"><math><![CDATA[<mrow><msub><mi>F</mi><mrow><mi>f</mi><mi>i</mi><mi>t</mi><mi>n</mi><mi>e</mi><mi>s</mi><mi>s</mi></mrow></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>&lambda;</mi><mrow><mo>(</mo><msub><mi>c</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><mi>f</mi><mo>(</mo><mrow><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>,</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mn>1</mn><mo>,</mo></mrow></msub><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>c</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>c</mi><mi>max</mi></msub></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000947135320000023.GIF" wi="982" he="143" /></maths>式中:λ为线性变换系数,f(x<sub>1</sub>,x<sub>2</sub>)为切削参数目标函数,c<sub>max</sub>为适应值阈值;(4.4)切削量的自适应遗传操作,包括选择、交叉和变异。
地址 214122 江苏省无锡市蠡湖大道1800号