最新新闻
我要投稿
联系电话:027-87592219/20/21转188
投稿邮箱:tb@e-works.net.cn
您所在的位置:首页 > 智库 > 智能装备

车削切削力不确定性的模糊

发布时间:2012-02-28 作者:数控  来源:360cnc
加工过程中出现的不确定性现象已成为一个突出的问题。切削力作为过程特征量之一,其不确定性来自许多过程因素(包括可量化、难量化因素以及可预测、难预测因素)。若忽略难量化和难预测因素的影响,可近似认为切削力的 ...

 

       加工过程中出现的不确定性现象已成为一个突出的问题。切削力作为过程特征量之一,其不确定性来自许多过程因素(包括可量化、难量化因素以及可预测、难预测因素)。若忽略难量化和难预测因素的影响,可近似认为切削力的不确定性主要源于两个方面:实际的力测量系统的影响(包括标定方法、漂移补偿方法、软件中的平均值和圆整方法、测力仪温度等的影响)和过程参数变化的影响(包括切削用量三要素的影响)。 常见的对系统不确定性的估计方法有概率分布法、模糊估计法和最小二乘法等,但至今所见关于加工过程不确定性研究的文章或报道为数不多。为了评价切削力的不确定性,考虑测量系统的标定误差和切削参数误差,根据最小二乘回归法和误差补偿理论,有研究者提出了TLSRM模型,此模型产生的不确定性预测结果是一个确定值而不是一个范围;考虑磨削参数变化和忽略测量系统的标定误差影响,基于模糊理论和线性规划理论,有关文献提出了一个模糊预测模型GFPM,而此模型产生的磨削力不确定性预测结果是一个范围而非一个确定值。 模糊预测模型GFPM使用区间形式来预测和评价磨削力的不确定性,从工程实践的角度看,要比预测一个确定点更有意义。但是模型GFPM不确定性预测的相对误差仍然较大。为了进一步研究这个问题,笔者注意到灰色理论擅长于估计任意系统的不确定性,因而试图建立一种模糊—灰色预测方法,并将其应用于车削力不确定性预测。为了验证所提出的方法,作者还建立了一个车削力检测系统。 2 车削力不确定性的模糊—灰色预测原理车削中,主切削力的经验数学模型可表达为

Fz=Cz(ap)x(f)y (1)

式中,ap为切削深度(mm),f 为进给量(mm/r),x和y为指数,Cz为系数。 对式(1)取对数,可得

log(Fz)=log(Cz)+xlog(ap)+ylog(f) (2)

令:Y=log(Fz),b0=log(Cz),b1=x,X1=log(ap),b2=y,X2=log(f),代入式(2)得线性模型分明形式为

Y=b0+b1X1+b2X2 (3)

定义式(3)对应的模糊形式为

Y(Xp)=A0+A1X1p+A2X2p (4)

这里选择对称三角形形式的模糊隶属度函数。设三角形底边的中心和宽度分别为C和W,则对称模糊数表达为Y(Xp)=[C(Xp),W(Xp)],A0=[C0,W0],A1=[C1,W1],A2=[C2,W2]。分明变量为X1p=log(ap)p,X2p=log(f)p(p表示试验代码,且p=1,2,3,…,n)。 为了使力的测量值落在式(4)所得预测区间内,并且使对称模糊数的宽度Wi之和最小,应当满足以下线性规划关系: 规划目标为

∑W(Xp)→min. (5)

约束条件为

log(Fz)p≤C(Xp)+(1-h)W(Xp) (6)

log(Fz)p≥C(Xp)-(1-h)W(Xp) (7)

式中,W(Xp)=W0+W1X1p+W2X2p,C(Xp)=C0+C1X1p+C2X2p,宽度Wi≥0,且中心Ci≥0(下标i=0,1,2)。模糊隶属度h 满足0Ys0={Ys0(1),Ys0(2),…,Ys0(j),…,Ys0(q)}其中j=1,2,…,qYs1={Ys1(1),…,Ys1(j),…,Ys1(q)}其中j=1,2,…,q

Yk1(j)=∑Yk0(i) (8)

其中j=1,2,…,q

Zu1={Zu1(2),…,Zu1(i),…,Zu1(q)} (i=2,…,q

Zs1={Zs1(2),…,Zs1(i),…,Zs1(q)}

其中,元素Zu1(i)和Zs1(i)可表达为 Zk1(i)=0.5Yk1(i)+0.5Yk1(i-1)(k=u 或s) 假设列向量Y=[Yk0(2),Yk0(3),…,Yk0(q)]T,且矩阵B定义为 B=[-Zk1(2),1;-Zk1(3),1;…;-Zk1(q),1](k=u 或s) 灰色微分方程dYk1/dt+aYk1(t)=b 中的参数a和b的估计值确定为:

[a,b]T=(BTB)-1 BTY (9)

根据式(8),有Yk1(0)=Yk0(1)。 灰色微分方程的解为

Yk1(1)=Yk0(1) (10a)

Yk1(i)=[Yk0(1)-b/a]exp[-a(i-1)]+b/a (10b)

式中,i=2,3,…,q。 式(10)可用于求Yk1的模拟序列Yk1。因此,Yk0的模拟序列Yk0可确定为

Yk0(1)=Yk0(1)(11a)

Yk0(i)=Yk1(i)-Yk1(i-1) (11b)

其中i=2,3,…,q;k=u 或s。i≤q 用于原始序列Yu0和Ys0的模拟;i>q 用于切削力不确定性预测。 3 切削力检测试验

表1 检测系统主要组成部分   组成部分说明

1 测量传感器电阻应变式

2 信号放大器 YD-4A电阻应变式电压放大

3 数字测量仪器 12位A/D 转换卡

4 检测软件用C++和MASM 语言开发

5 采集参数采样频率50Hz;样本500~1000

6 计算机 586微机

 

表2 车削力试验记录

No. 切削深度ap
(mm) 进给量f
(mm/r) 主切削力Fz
(N)

1 2 0.1 439

2 2 0.2 878

3 2 0.3 1129

4 2 0.4 1443

5 2 0.5 1756

 

  • 切削力检测试验条件 为验证模糊—灰色预测方法的有效性,建立了切削力检测系统,以获取建模和评价数据。检测系统的主要组成部分见表1。其它试验条件包括:工件材料45钢(正火处理,硬度HB187),工件直径81mm;CA6140车床,主轴转速380r/min;切削速度96m/min;YT15硬质合金刀具。
    • 切削力试验与检测 基于上述切削力检测试验条件和表2所列切削参数进行切削力试验,并将切削力检测结果记录在表2中。4 预测和试验评价
    • 切削力不确定性的模糊回归模型 根据式(4)~(7),用表2中全部数据进行切削力不确定性模糊建模(考虑到在建模数据范围以外模糊估计较差)。为便于编程计算,令m=2且X1p=log(10mf)p。表3 列出了切削力不确定性的模糊估计结果,对应的模糊回归模型为 Y(Xp)=[logFz]p=A0+A1 [log(ap)]p+A2 [log(100f)]p式中,A0=[1.272603,0.1092808],A1=[0.8845206,0],A2=[0.856115,0]。模糊隶属度h=0.5。相对误差∆u(或∆s)定义为 ∆=(估计值-测量值)/测量值 #p#分页标题#e#

      表3 车削力不确定性模糊估计结果

      No. 实测力模糊估计下限∆u(%) 上限∆s(%)

      1 439 [415,516] -5.5 +17.5

      2 878 [749,932] -14.7 +6.20

      3 1129 [1063,1324] -5.8 +17.3

      4 1443 [1363,1696] -5.5 +17.5

      5 1756 [1646,2047] -6.35 +16.6

       

      表4 车削力不确定性模糊—灰色预测
      (含3个建模数据的部分信息模型)

      No. 实测力模糊—灰色预测下限∆u(%) 上限∆s(%)

      1 439 [415.00,516.00] -5.5 +17.5

      2 878 [740.10,920.80] -15.7 +4.9

      3 1129 [1046.6,1303.4] -7.3 +15.4

      4 1443 [1480.1,1845.1] +2.6 +27.9

      5 1756 [2093.2,2611.8] +19.2 +48.7

       

      表5 车削力不确定性模糊—灰色预测
      (含4个建模数据的新信息模型)

      No. 实测力模糊—灰色预测下限∆u(%) 上限∆s(%)

      1 439 [415.00,516.00] -5.5 +17.5

      2 878 [765.60,953.30] -12.8 +8.5

      3 1129 [1020.6,1270.5] -9.6 +12.5

      4 1443 [1360.4,1693.2] -5.7 +17.3

      5 1756 [1813.5,2256.6] +3.3 +28.5

       

      表6 车削力不确定性模糊—灰色预测
      (含3个建模数据的新陈代谢模型)

      No. 实测力模糊—灰色预测下限∆u(%) 上限∆s(%)

      2 878 [749.00,932.00] -14.7 +6.2

      3 1129 [1056.9,1316.4] -6.4 +16.6

      4 1443 [1353.4,1684.2] -6.2 +16.7

      5 1756 [1733.2,2154.6] -1.3 +22.7

       

    • 车削力不确定性的模糊—灰色预测
      • 选择表3中No.1~3的数据用于模糊—灰色预测建模,则车削力不确定性的模糊—灰色预测结果如表4 所示,其中数据No.1~3对应于模糊—灰色建模数据的模拟值(建模数据来自模糊估计数据No.1~3),而数据No.4~5表示模糊—灰色建模数据范围之外的预测值。此外,灰色微分方程中的参数确定为 [as,bs]=[-0.3475,590.7376]
        [au,bu]=[-0.3466,475.3764] 相对误差定义为 ∆=(预测值-测量值)/测量值

         

        • 选择表3中的数据No.1~4用于模糊—灰色预测建模,将车削力不确定性模糊—灰色预测结果列于表5,其中数据No.1~4对应于模糊—灰色建模数据的模拟值(建模数据来自模糊估计数据No.1~4),而数据No.5 表示模糊—灰色建模数据范围之外的预测。灰色微分方程中的参数确定为 [as,bs]=[-0.2872,674.6824]
          [au,bu]=[-0.2874,541.5835]
          • 选择表3中的数据No.2~4用于模糊—灰色预测建模,则车削力不确定性模糊—灰色预测结果列于表6,其中数据No.2~4对应于模糊—灰色建模数据的模拟值(建模数据来自模糊估计数据No.2~4),而数据No.5 表示模糊—灰色建模数据范围之外的预测值。灰色微分方程中的参数确定为 [as,bs]=[-0.2464,931.306]
            [au,bu]=[-0.2473,746.3059]
          5 预测结果分析由表4数据可知:数据No.1~3是原始数据的较好模拟(原始数据来自表3中No.1~3),但预测区间稍微向左漂移;虽然数据No.4在建模数据范围之外,但它的预测较好且含较小的相对误差(+2.6%),即实测力落在预测区间下限的左侧邻近;数据No.5的预测较差,即实测力落在预测区间下限的左侧且含较大的相对误差(+19.2%)。 由表5数据可知:数据No.1~4是原始数据的较好模拟(原始数据来自表3中数据No.1~4);虽然数据No.5在建模数据范围之外,但它的预测较好且含较小的相对误差(+3.3%),即实测力落在预测区间下限的左侧邻近。 由表6数据可知:数据No.2~4是原始数据的较好模拟(原始数据来自表3中数据No.2~4);虽然数据No.5是在建模数据范围之外,但它的预测较好且含较小的相对误差(-1.3%);实测力落在预测区间下限的右侧邻近。 基于上述分析可知:模糊—灰色预测方法可以较好地模拟建模数据,同时可以在建模数据样本范围之外有效地预测切削力不确定性区间,不过真正的有效预测点仅局限于最后那个建模数据。此外还可发现:表6中的新陈代谢模型的预测结果优于表4中的部分信息模型,而表4的部分信息模型又优于表5中的新信息模型。与常规的最小二乘回归法、神经网络法和模糊回归法明显不同,模糊—灰色预测法仅需要少数试验数据样本来建模,并能有效地预测将来的不确定性区间。 6 结论在建立的车削力检测系统中验证了所提出的切削力不确定性模糊—灰色预测方法的有效性。模糊—灰色预测结果与实测切削力数据的比较表明:所提出的预测新方法不仅能较好地模拟建模数据样本,而且还能以较小的相对误差预测建模数据范围以外的其它数据点(即预测未来数据点)。此外,与传统建模方法明显不同的是:模糊—灰色预测方法在建模中仅需要少量试验数据样本,可大大减少切削试验的成本和工作量。

  • 相关资料推荐