CN106407669A - 一种切削加工表面粗糙度的预测方法 - Google Patents
一种切削加工表面粗糙度的预测方法 Download PDFInfo
- Publication number
- CN106407669A CN106407669A CN201610807001.6A CN201610807001A CN106407669A CN 106407669 A CN106407669 A CN 106407669A CN 201610807001 A CN201610807001 A CN 201610807001A CN 106407669 A CN106407669 A CN 106407669A
- Authority
- CN
- China
- Prior art keywords
- surface roughness
- variable
- function
- distribution function
- parameter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Abstract
本发明公开了一种切削加工表面粗糙度的预测方法,属于切削加工技术领域,S1、选取一定数量的变量参数和对应的表面粗糙度作为样本数据;S2、对样本数据进行预处理,确定出变量参数和表面粗糙度的变量概率分布函数;S3、对待选Copula函数的模型参数进行估计分析,选择出最优的Copula函数;S4、把最优Copula函数和变量概率分布函数合成,确定出变量参数和表面粗糙度的联合分布函数,并推导出基于变量参数的表面粗糙度条件概率分布函数;S5、通过条件概率分布函数,对变量参数与表面粗糙度进行局部相关性分析,校正预测模型;S6、通过校正后的预测模型来计算表面粗糙度的预测值。本方法模型简单,所需样本较少,自适应与自组织能力较强,泛化能力较好,精度较高。
Description
技术领域
本发明属于切削加工技术领域,尤其是涉及到一种切削加工表面粗糙度的预测方法。
背景技术
在机械加工中,表面粗糙度不仅是进行零件设计的重要要求之一,同时也是衡量产品质量的重要指标。但是,在切削过程中,由于工况的不平稳性,很难准确地控制零部件表面粗糙度,而产品的耐腐蚀性、润湿性、抗疲劳性以及配合性质等等都与表面粗糙度有关,因此,如何有效地控制工件的表面粗糙度成为了机械制造领域的一个关键问题。
目前,表面粗糙度的控制方法主要还是依据经验来拟定加工参数,然后通过离线测量手段在不同时间段内完成相应的表面粗糙度测量工作。但是,由于参数设定值较为保守以及离线测量后产生的重新安装误差,导致生产效率低,同时还无法准确地保证产品的表面质量。于是,为了保证零部件的加工表面质量,减少不可预知的偶然性错误或者偏差,国内外许多学者对表面粗糙度预测技术进行了深入研究。目前关于表面粗糙度预测研究方法主要分为三种:基于切削理论模型、基于切削参数的经验模型以及基于人工智能的预测方法。这些方法基本上都是通过切削用量来预测表面粗糙度,并且一般都没有对特征变量与表面粗糙度输出变量进行合理的相关性分析。
发明内容
针对智能制造中通过在线监测切削加工过程来判定系统的状态及预测,提出了基于Copula函数的表面粗糙的预测新方法。
本发明是通过如下技术方案得以实现的:
一种切削加工表面粗糙度的预测方法,以监测信号作为变量参数,
S1、选取一定数量的变量参数和对应的表面粗糙度作为样本数据;
S2、对样本数据进行预处理,确定出变量参数和表面粗糙度的变量概率分布函数;
S3、对待选Copula函数的模型参数进行估计分析,选择出最优的Copula函数;
S4、把最优Copula函数和变量概率分布函数合成,确定出变量参数和表面粗糙度的联合分布函数,并推导出基于变量参数的表面粗糙度条件概率分布函数;
S5、通过条件概率分布函数,对变量参数与表面粗糙度进行局部相关性分析,校正预测模型;
S6、通过校正后的预测模型来计算表面粗糙度的预测值。
进一步的,所述步骤S1中,所述变量参数为切削力。
进一步的,所述步骤S1中,样本数据的数量不少于30个。
进一步的,所述步骤S1中,对样本数据采用取自然对数的方法进行预处理。
进一步的,所述步骤S2中,变量参数和表面粗糙度的变量概率分布函数为正态分布函数。
进一步的,所述步骤S3中,利用极大似然法计算对待选的Copula函数中的建模参数进行估计分析,通过计算不同Copula函数与经验分布函数的欧式距离d2作为检验其拟合优度的评价标准,选用d2值最小的为最优的Copula函数。
进一步的,步骤S6中,把基于变量参数的条件概率分布的条件期望值作为表面粗糙度的预测值。
进一步的,所述步骤S5中,在进行相关性分析之前,选用相同水平α范围内的分位数来对变量参数和表面粗糙度的大小进行分类,根据分类结果,应用条件概率计算模型就可以得到不同取值范围内的变量参数与表面粗糙度的相关关系。
更进一步的,所述步骤S5中,根据样本数据局部相关性分析结果,确定出监测切削力的水平范围,校正预测模型的积分范围,然后通过校正后的预测模型来计算表面粗糙度的预测值。
有益效果:
1.本方法模型简单,所需样本较少,自适应与自组织能力较强,泛化能力较好,精度较高,为自动化加工和机械加工过程在线监测提供了一个新的视角。
2.通过在线监测运行状态下的工艺系统,能够实现表面加工质量的预测预报,对于实现智能制造和工业4.0具有重要的科学意义和较大的实用价值。
3.根据样本数据局部相关性分析结果,确定出监测切削力的水平范围,从而校正预测模型的积分范围。
附图说明
图1为本发明的表面粗糙度预测方法的流程图;
图2为本发明的F与Ra两变量观测数据经验累积概率与理论累积概率的散点图。
具体实施方式
下面结合附图1和2以及具体实施例对本发明作进一步的说明,但本发明的保护范围并不限于此。
本实例中选用切削力为监测变量参数。
1)样本数据的获取
为了较好地描述切削力与表面粗糙度两变量的相关性,选取不同切削条件下所对应的切削力和表面粗糙度值作为样本数据,如表1所示。
表1 车削试验结果
2)数据预处理
切削力和表面粗糙度数值大小及变化范围,相差较大。为了保证样本数据的平稳性,得到满意的分析结果,需要对切削力F和表面粗糙度Ra数据进行预处理,取自然对数处理,之后所有分析过程都利用预处理过后的数据。
3)变量边缘分布函数的确定
切削加工过程一般都服从正态分布,因此假设切削力与表面粗糙度变量都服从正态分布。为了验证假设的合理性,首先利用极大似然法估计两变量概率分布函数的未知参数值,然后运用Kolmogorov-Smirnov(简称K-S检验法)检验法来进行相应的分析。其参数值及检验统计量见表2。
表2 样本的概率分布参数与K-S检验统计量
在显著水平α为0.05时,则统计量D(64,0.05)=0.17,而切削力和表面粗糙度变量所对应的统计量D都小于0.17,且相伴概率P值都大于0.05,则不能拒绝原假设,从而说明切削力与表面粗糙度变量都服从正态分布,各变量的参数值见表3。
表3 变量边缘分布参数
4)Copula函数的参数估计及选择
选用Gaussian Copula、t-Copula、Gumbel Copula、Clayton Copula和FrankCopula五种常用的Copula函数来拟合切削力与表面粗糙度两变量的相关结构分布形式,由于备选的五种Copula函数中都含有未知参数,因此,首先利用极大似然法计算其中的参数。由于切削力与表面粗糙度变量都服从正态分布,所以首先将各变量样本值转换为[0,1]区间内的概率分布值(用U,V表示),然后计算得到Copula函数中所包含的未知参数估计值,如表4所示。
表4 切削力变量与表面粗糙度变量的不同Copula函数建模参数
然后通过计算不同Copula函数与经验分布函数的欧式距离d2作为检验其拟合优度的评价标准,选用d2值最小的为最优的Copula函数。将经概率积分变换后的样本值代入到各个待选Copula函数以及样本的经验分布函数中,就可以获得样本数据对应的理论概率值和经验概率值,然后通过计算两者的欧式距离,结算结果如表5所示。由表5计算结果可知,Frank Copula函数可以较好地表征切削力与表面粗糙度的相关结构。
表5 F与Ra两变量各建模模型的拟合评价指标
5)联合分布函数模型的建立
结合Sklar定理及上述分析结果可知,切削力与表面粗糙度两变量的连接函数为Frank Copula函数,此外,它们的边缘分布都服从正态分布,则切削合力F变量(用x表示)与表面粗糙度Ra变量(用y表示)的联合分布函数分别如下式所示:
式中,α=4.1185
为了验证各切削力变量与表面粗糙度变量所建立的联合分布函数模型的合理性,可以通过比较各个观测样本数据(xi,yi)所对应的经验累积概率值与理论累积概率值是否近似相同。将切削力与表面粗糙度两相关变量各个观测样本数据(xi,yi)的理论累计概率Fthe值作为纵坐标,经验累计概率Femp值作为横坐标,分别画出各变量对应的散点图,如图2所示。各相关变量所对应的理论值与经验值基本在一条直线上,拟合效果好,说明建立的联合分布函数模型是合理的。
6)切削力与表面粗糙度的相关性分析
由于切削力与表面粗糙度所对应的分布函数并不相同,相同的数值对应的概率并不一致,则其在各自样本中所表征的大小程度也不一致,这样可能会导致错误的结论。因此,在进行相关性分析之前,将选用相同水平α范围内的分位数来对切削力和表面粗糙度的大小进行分类,其分类结果见表6。根据分类结果,应用条件概率计算公式(2)就可以得到不同取值范围内的切削力与表面粗糙度的相关关系,如表7所示。
表6 切削力和表面粗糙度样本数据的大小分类结果
表7 切削合力F与表面粗糙度Ra的条件概率
注:表7中qα表示水平值α所对应的分位数大小。
7)预测模型的建立
根据概率学理论推导出基于切削力X的表面粗糙度Y条件概率分布函数FY|X(y),其表达形式如式(3)所示。然后,通过计算Y的条件概率分布FY|X(y)的条件期望值作为表面粗糙度的预测值Rapre,其表达形式如式(4)所示。
式中fY|X(y)、c(u,v)及fY(y)分别是FY|X(y)、Copula函数和变量Y的概率密度分布函数。
由于切削力与表面粗糙度变量的相关函数是Frank Copula函数,其密度函数表达如式(5)所示
表面粗糙度预测模型的最终表达形式如式(6)所示。
式中α为常数,为v=FY(y)的反函数。
8)预测模型的校正
由于不同大小的切削力与表面粗糙度的相关程度不一致,因此在切削力F输入之前,需要根据局部相关性分析结果对所建立的预测模型进行校正,然后根据判定结果,通过调整式(6)中的积分上下限,从而校准了预测模型,实现精确的预测,其校正结果如表8所示。然后在利用校正过后的模型计算得到最终的表面粗糙度预测值Rapre,其测试结果如表9所示。
表8 模型校正结果
表9 表面粗糙度预测模型的分析结果
从表9中可以看出,该预测模型的平均误差率是5.6536%,且模型最高误差率为9.3589%(低于10%),低误差率说明了表面粗糙度预测值Rapre与测量值Ra的拟合情况较好。
所述实施例为本发明的优选的实施方式,但本发明并不限于上述实施方式,在不背离本发明的实质内容的情况下,本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。
Claims (9)
1.一种切削加工表面粗糙度的预测方法,其特征在于,以监测信号作为变量参数,
S1、选取一定数量的变量参数和对应的表面粗糙度作为样本数据;
S2、对样本数据进行预处理,确定出变量参数和表面粗糙度的变量概率分布函数;
S3、对待选Copula函数的模型参数进行估计分析,选择出最优的Copula函数;
S4、把最优Copula函数和变量概率分布函数合成,确定出变量参数和表面粗糙度的联合分布函数,并推导出基于变量参数的表面粗糙度条件概率分布函数;
S5、通过条件概率分布函数,对变量参数与表面粗糙度进行局部相关性分析,校正预测模型;
S6、通过校正后的预测模型来计算表面粗糙度的预测值。
2.根据权利要求1所述的表面粗糙度的预测方法,其特征在于,所述步骤S1中,所述变量参数为切削力。
3.根据权利要求1所述的表面粗糙度的预测方法,其特征在于,所述步骤S1中,样本数据的数量不少于30个。
4.根据权利要求1所述的表面粗糙度预测方法,其特征在于,所述步骤S1中,对样本数据采用取自然对数的方法进行预处理。
5.根据权利要求1所述的表面粗糙度的预测方法,其特征在于,所述步骤S2中,变量参数和表面粗糙度的变量概率分布函数为正态分布函数。
6.根据权利要求1所述的表面粗糙度的预测方法,其特征在于,所述步骤S3中,利用极大似然法对待选的Copula函数中的建模参数进行估计分析,通过计算不同Copula函数与经验分布函数的欧式距离d2作为检验其拟合优度的评价标准,选用d2值最小的为最优的Copula函数。
7.根据权利要求1所述的表面粗糙度的预测方法,其特征在于,步骤S6中,把基于变量参数的条件概率分布的条件期望值作为表面粗糙度的预测值。
8.根据权利要求1所述的表面粗糙度的预测方法,其特征在于,所述步骤S5中,在进行相关性分析之前,选用相同水平α范围内的分位数来对变量参数和表面粗糙度的大小进行分类,根据分类结果,应用条件概率计算模型就可以得到不同取值范围内的变量参数与表面粗糙度的相关关系。
9.根据权利要求1至8任一项所述的表面粗糙度的预测方法,其特征在于,所述步骤S5中,根据样本数据局部相关性分析结果,确定出监测切削力的水平范围,校正预测模型的积分范围,然后通过校正后的预测模型来计算表面粗糙度的预测值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610807001.6A CN106407669A (zh) | 2016-09-07 | 2016-09-07 | 一种切削加工表面粗糙度的预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610807001.6A CN106407669A (zh) | 2016-09-07 | 2016-09-07 | 一种切削加工表面粗糙度的预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106407669A true CN106407669A (zh) | 2017-02-15 |
Family
ID=57999931
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610807001.6A Pending CN106407669A (zh) | 2016-09-07 | 2016-09-07 | 一种切削加工表面粗糙度的预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106407669A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108805142A (zh) * | 2018-05-31 | 2018-11-13 | 中国华戎科技集团有限公司 | 一种犯罪高危人员研判方法及系统 |
CN108920812A (zh) * | 2018-06-28 | 2018-11-30 | 江苏大学 | 一种机械加工表面粗糙度预测方法 |
CN109856337A (zh) * | 2019-02-28 | 2019-06-07 | 北京航空航天大学 | 一种碳纤维复合材料加工表面质量评定方法及装置 |
CN111781824A (zh) * | 2020-05-26 | 2020-10-16 | 华东理工大学 | 一种基于vine copula分位数回归的自适应软测量方法及系统 |
CN116578832A (zh) * | 2023-07-13 | 2023-08-11 | 兰州理工大学 | 预测机械加工零件表面粗糙度的装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101537431A (zh) * | 2008-03-21 | 2009-09-23 | 宝山钢铁股份有限公司 | 冷轧薄带钢表面粗糙度在线预测和控制方法 |
CN102172990A (zh) * | 2011-01-11 | 2011-09-07 | 哈尔滨工业大学 | 一种单点金刚石车削加工中表面粗糙度的预测方法 |
CN102880771A (zh) * | 2012-10-31 | 2013-01-16 | 贵州大学 | 高速切削加工中工件表面粗糙度的预测方法 |
CN103761429A (zh) * | 2014-01-10 | 2014-04-30 | 大连理工大学 | 铣削加工工件表面粗糙度的预测方法 |
CN104318022A (zh) * | 2014-10-28 | 2015-01-28 | 湘潭大学 | 一种预测工件表面粗糙度并提高切削效率的方法 |
-
2016
- 2016-09-07 CN CN201610807001.6A patent/CN106407669A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101537431A (zh) * | 2008-03-21 | 2009-09-23 | 宝山钢铁股份有限公司 | 冷轧薄带钢表面粗糙度在线预测和控制方法 |
CN102172990A (zh) * | 2011-01-11 | 2011-09-07 | 哈尔滨工业大学 | 一种单点金刚石车削加工中表面粗糙度的预测方法 |
CN102880771A (zh) * | 2012-10-31 | 2013-01-16 | 贵州大学 | 高速切削加工中工件表面粗糙度的预测方法 |
CN103761429A (zh) * | 2014-01-10 | 2014-04-30 | 大连理工大学 | 铣削加工工件表面粗糙度的预测方法 |
CN104318022A (zh) * | 2014-10-28 | 2015-01-28 | 湘潭大学 | 一种预测工件表面粗糙度并提高切削效率的方法 |
Non-Patent Citations (1)
Title |
---|
陈少峰: ""基于Copula函数的精密车削表面粗糙度预测研究"", 《中国学位论文全文数据库》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108805142A (zh) * | 2018-05-31 | 2018-11-13 | 中国华戎科技集团有限公司 | 一种犯罪高危人员研判方法及系统 |
CN108920812A (zh) * | 2018-06-28 | 2018-11-30 | 江苏大学 | 一种机械加工表面粗糙度预测方法 |
CN109856337A (zh) * | 2019-02-28 | 2019-06-07 | 北京航空航天大学 | 一种碳纤维复合材料加工表面质量评定方法及装置 |
CN111781824A (zh) * | 2020-05-26 | 2020-10-16 | 华东理工大学 | 一种基于vine copula分位数回归的自适应软测量方法及系统 |
CN111781824B (zh) * | 2020-05-26 | 2022-08-09 | 华东理工大学 | 一种基于vine copula分位数回归的自适应软测量方法及系统 |
CN116578832A (zh) * | 2023-07-13 | 2023-08-11 | 兰州理工大学 | 预测机械加工零件表面粗糙度的装置 |
CN116578832B (zh) * | 2023-07-13 | 2023-09-12 | 兰州理工大学 | 预测机械加工零件表面粗糙度的装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106407669A (zh) | 一种切削加工表面粗糙度的预测方法 | |
CN110303380B (zh) | 一种数控机床刀具剩余寿命预测方法 | |
CN101863088B (zh) | 一种橡胶混炼过程中门尼粘度的预报方法 | |
CN106681183B (zh) | 监测制造装置的方法、装置、系统与计算机可读储存媒介 | |
CN111508216A (zh) | 一种大坝安全监测数据智能预警方法 | |
TW201642147A (zh) | 具取樣率決定機制的量測抽樣方法 與其電腦程式產品 | |
CN104932488A (zh) | 一种模型预测控制性能评估与诊断方法 | |
CN107627152B (zh) | 一种基于bp神经网络的数控加工切屑控制方法 | |
CN102880809A (zh) | 基于关联向量回归模型的聚丙烯熔融指数在线检测方法 | |
CN105867341B (zh) | 一种烟草加工设备的在线设备健康状态自检方法及系统 | |
Kim et al. | Wafer edge yield prediction using a combined long short-term memory and feed-forward neural network model for semiconductor manufacturing | |
Ak et al. | Data analytics and uncertainty quantification for energy prediction in manufacturing | |
CN105955198B (zh) | 基于最小二乘迭代算法的机床工步能耗监测方法 | |
CN110298765B (zh) | 一种基于客观关联因素的配电网用电功率异常检测方法 | |
CN113539382B (zh) | 一种亚磷酸二甲酯关键工艺参数的预警定位方法及系统 | |
CN103605323B (zh) | 化工生产的离散控制方法及装置 | |
EP3906447B1 (en) | Method, system and computer program product for evaluation of energy consumption in industrial environments | |
KR102110319B1 (ko) | 학습 데이터 생성 시스템 | |
CN116540523A (zh) | 一种基于关键指标的pid回路参数智能整定方法及系统 | |
Stebel et al. | Performance improvement for quasi periodical disturbances in PH control | |
CN113312587A (zh) | 基于arima预测和回归预测的传感器采集数据缺失值处理方法 | |
CN113536512A (zh) | 一种弹体质量智能预测分析方法 | |
Bhinge et al. | An Intelligent Machine Monitoring System Using Gaussian Process Regression for Energy Prediction | |
CN114971599B (zh) | 食品制造加工工艺参数设置方法、装置、设备及介质 | |
Karandikar et al. | Remaining useful tool life predictions using Bayesian inference |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170215 |
|
RJ01 | Rejection of invention patent application after publication |