CN102608302A - 一种橡胶硬度的测量方法 - Google Patents
一种橡胶硬度的测量方法 Download PDFInfo
- Publication number
- CN102608302A CN102608302A CN2012100743611A CN201210074361A CN102608302A CN 102608302 A CN102608302 A CN 102608302A CN 2012100743611 A CN2012100743611 A CN 2012100743611A CN 201210074361 A CN201210074361 A CN 201210074361A CN 102608302 A CN102608302 A CN 102608302A
- Authority
- CN
- China
- Prior art keywords
- hardness
- new
- theta
- rubber
- value
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种橡胶硬度的测量方法,当用于模型更新的新样本Xnew的数目达到阈值T时,硬度在线预测模型f:x →y更新一次;当没有达到阈值T时,不进行更新,直到达到阈值T。获取更新后的硬度在线预测模型,通过更新后的硬度在线预测模型预测下一车次的新样本Xnew对应的橡胶硬度值,获取橡胶硬度值并与实际测得的橡胶硬度值的差值,获得误差;本方法通过不断的对硬度在线预测模型进行更新,从而得到精确的橡胶硬度值,提高了橡胶的质量,并且通过实际的试验验证了采用本方法可以得到较高的精度,减少了安全隐患。
Description
技术领域
本发明涉及橡胶轮胎制造领域,特别涉及一种橡胶硬度的测量方法。
背景技术
近年来,随着生产和生活中对橡胶需求的不断增加,橡胶行业有着突飞猛进的发展。加上中国汽车行业的不断发展,对橡胶行业更起到了巨大的推动作用。不论是在生产还是生活领域,越来越多的使用到橡胶制品。新产品的不断出现,应用领域的不断扩大,均对橡胶产品的生产提出了更新更高的要求。
橡胶混炼是橡胶生产最主要的工序之一,对橡胶混炼过程的有效控制直接影响到橡胶产品的质量。橡胶混炼是典型的间歇过程,具有明显的时变特性和非线性,因此对混炼胶质量检测的及时、准确,即对混炼过程的实时控制,是保证橡胶产品质量的关键。但是橡胶混炼过程很短,一般只有2~3分钟左右的时间,而且混练过程复杂,所以对混炼过程直接进行控制具有很大的挑战,因此也很难保障橡胶的质量。硬度作为衡量胶料粘度特性、可加工特性和压延特性等多方面性能的综合物性指标,是橡胶质量检测最关键的指标之一。
目前,橡胶行业的测量过程主要分为以下几个步骤:
(1)橡胶通过密炼、挤出、压延、开炼后,经过隔离剂,再进行风冷,风冷后进行垫片;(2)停放一定时间(一般在四小时以上)以保证胶料物性趋于稳定后,工艺人员手动取样并记录其车次信息;(3)将取得的终炼胶样品拿到快检室进行手动冲样;(4)用橡胶硬度计测量橡胶样品的硬度。
经过上述一系列的复杂工序后,才能得到该车次橡胶的硬度。
近年来,基于软模型的质量预测方法受到广泛的应用,软测量的核心问题是建立测量模型,也就是建立待估计变量与其它直接测量变量间的关联模型。其中部分最小二乘法(Partial Least Square PLS)不依赖于过程机理模型,训练时不需要故障样本,近年来在化工过程的质量控制和在线监测等方面得到广泛应用和研究。现有PLS技术在生产过程中,采集包含大量相关信息的数据,根据PLS相关理论只建立一次软测量模型,对整个生产过程的产品进行质量预测。由于橡胶生产过程是一个典型的间歇性过程,通常只进行一次模型的建立不能完成对整个生产过程的质量实时监控,造成工厂橡胶产品质量受到严重影响,最终影响产品的质量和寿命。而且橡胶产品可能会有安全性隐患,给企业和产品使用者造成经济损失。
发明内容
本发明要解决的技术问题在于提供一种橡胶硬度的测量方法,该方法可以实现实时的测量硬度,测量结果准确和生产成本低等优点,详见下文描述:
一种橡胶硬度的测量方法,所述方法包括以下步骤:
(1)在不同置信水平α下设定Q统计量的三个控制限Qas、Qal和Qah,其中Qas<Qal<Qah;
(2)将当前车次混炼过程结束后质量监测系统采集到的终炼胶质量指标中的流变参数与硬度在线预测模型f:x →y自动连接,并且设定流变参数的新样本Xn ew;
(3)计算所述新样本Xnew对应的Q统计量的值Qnew;
(4)判断是否满足Qnew>Qah,如果是,执行步骤(5);如果否,执行步骤(6);
(5)与所述硬度在线预测模型f:x →y相连接的警报器发生警报,操作工人查看生产线的工况,然后进行相应的工艺调整,使当前工况恢复到正常状况,流程结束;
(6)判断是否满足Qas≤Qnew≤Qal,如果是,执行步骤(7);如果否,执行步骤(9);
(7)则所述新样本Xnew用于所述硬度在线预测模型f:x →y的更新,当用于模型更新的新样本Xnew的数目达到阈值T时,所述硬度在线预测模型f:x →y更新一次;当没有达到阈值T时,不进行更新,直到达到阈值T,获取更新后的硬度在线预测模型,执行步骤(8);
(8)通过所述更新后的硬度在线预测模型预测下一车次的新样本Xnew对应的橡胶硬度值,获取所述橡胶硬度值与实际测得的橡胶硬度值的差值,获得误差;
(9)流程结束。
式中: (i=1,2,3...)
λj为X的协方差矩阵的特征值,Cα为正态分布在检验水平为α下的临界值,k是主元模型中保留的主元个数,n是全部主元个数。
本发明提供的技术方案的有益效果是:本发明实施例给出了一种橡胶硬度的在线测量方法,本方法通过不断的对硬度在线预测模型进行更新,从而得到精确的橡胶硬度值,提高了橡胶的质量,并且通过实际的试验验证了采用本方法可以得到较高的精度,减少了安全隐患。
附图说明
图1为本发明提供的一种橡胶硬度的测量方法的流程图;
图2为本发明提供的橡胶硬度模型预测试验示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
为了提高生产效率和产品的质量,实现在线自动测量,本发明实施例提供了一种橡胶硬度的测量方法,参见图1,详见下文描述:
橡胶的质量指标流变参数(主要包括最低扭矩ML,最高扭矩MH和硫化时间参数TC30、TC60以及TC90)是每车橡胶在混炼结束后必须要测量的参数,每次测量只需要2~5分钟,检测快速简便。考虑到邵氏硬度计和流变仪工作原理的相似性,以及流变参数与硬度本身具有较强的相关性,那么如何通过流变参数拟合出硬度成为解决硬度测量瓶颈问题的一个可行方法。
目前为止,多变量统计回归方法在数据驱动问题中已经广泛应用,而且对提高预测性能起到了明显的作用,解决了很多生产生活中的实际问题。橡胶工业上橡胶硬度的检测多为多变量统计回归过程。用于回归的方法很多种,如传统的统计分析方法以及神经网络方法等。这些方法虽然在实际应用中占据主导地位,但人们也发现它们还存在着许多不足之处。比如,传统的统计方法一般需要事先知道样本的先验分布,并要求有足够多的样本数据,而这些要求在实际应用中往往难以达到,这就使其在实际应用中效果往往并不理想。神经网络方法虽然很好地解决了非线性问题,但由于其自身存在着结构不易确定、易陷入局部极小等固有的缺陷,从而限制了其实际应用。另外,神经网络的学习算法仅仅试图使经验风险最小化,并没有使期望风险最小化,与传统的最小二乘法相比,在原理上缺乏实质性的突破,这也是神经网络过拟合现象产生的原因,从而导致了其推广能力的下降。
101:在不同置信水平α下设定Q统计量的三个控制限Qas、Qal和Qah,其中Qas<Qal<Qah;
在多变量控制过程中,需要建立一个反映正常运行的主元模型,将反映正确运行的历史数据收集起来,对这些数据进行主元分析,建立主元模型。由于主元分析的结构受数据本身的影响,因此在进行主元分析时,需要先将数据进行标准化,即将每个变量的均值减掉然后除以它的标准差。假设X∈Rn×m为收集到的过程正常运行数据,先将X做如下标准化,
式中M=[m1,m2...mm]为变量X的均值,s=[s1,s2...,sm]为变量的标准差。对进行主元分析可以得到
如果可以用前k(k<m)个主元来代表数据中的主要变化,那么可以得到下面的主元模型
其中
这个主元在i时刻的平方预测误差(SPE,Square Prediction Error)可以写为
为了利用主元模型对生产过程进行监控,需要工业过程正常运行的数据来确定工业过程运行的SPE控制限,当主元模型的SPE超过它们的控制限时,就认为工业过程中出现了不正常情况。
主元模型的SPE也称作Q统计量,对第i个采样点来说
其中ei是式E的第i行,Pk=[p1,p2,Lpk],I是n×n的单位矩阵。SPE代表了数据没有被主元模型解释的变化。当SPE过大时,说明过程中出现了不正常情况,从而过程正常运行时的模型不再试用。控制线的计算是建立在一定的假设基础上的,当置信水平为α时,SPE控制限可以按下式计算:
式中: (i=1,2,3...)
λj为X的协方差矩阵的特征值,Cα为正态分布在置信水平为α下的临界值,k是主元模型中保留的主元个数,n是全部主元个数。其中,三个控制限Qas、Qal和Qah对应于置信α取三个不同值时得到的函数表达式。
Q统计量度量了数据点偏离主元模型的程度,当Q统计量超过控制限时,可以判断操作出现了不正常情况。
102:将当前车次混炼过程结束后采集到的终炼胶质量指标中的流变参数与硬度在线预测模型f:x →y自动连接,并且设定流变参数的新样本Xnew;
其中,硬度在线预测模型f:x→y具体为:
1)采集终炼胶质量指标数据样本,建立数据库,数据样本包括流变参数和硬度值,且规定流变参数对应自变量X={x1,x2,...,xn},硬度值对应因变量Y={y1,y2,...,yn},n为建模所需样本数,xi和yi分别对应流变参数和硬度值,i=1,2,...,n;X与Y之间存在线性关系f:x→y
Y=XC+V
式中:V和C分别代表相应维数的矩阵和系数矩阵。
2)运用PLS回归算法通过对X和Y进行线性分解建立如下回归模型:
PLS回归系数矩阵为:
C=(XTX)+XTY=W(P′W)-1Q′
式中:算子(XTX)+表示广义逆运算(伪逆)。
3)将输入带入函数f:x →y中,得到相应的模型预测值。
103:计算新样本Xnew对应的Q统计量的值Qnew;
104:判断是否满足Qnew>Qah,如果是,执行步骤105;如果否,执行步骤106;
105:与在线预测模型f:x →y相连接的警报器发生警报,操作工人查看生产线的工况,然后进行相应的工艺调整,使当前工况恢复到正常状况,执行步骤109;
其中,相应的工艺调整包括:原料反应温度和混料顺序进行调节。
106:判断是否满足Qas≤Qnew≤Qal,如果是,执行步骤107;如果否,执行步骤109;
107:则新样本Xnew用于硬度在线预测模型f:x →y的更新,当用于模型更新的新样本Xnew的数目达到阈值T时,硬度在线预测模型f:x →y更新一次;当没有达到阈值T时,不进行更新,直到达到阈值T,获取更新后的硬度在线预测模型,执行步骤108;
其中,阈值A取值根据实际生产状况确定,具体实现时,本发明实施例对此不做限制。例如:当阈值T的取值为5时,当前车次用于模型更新的新样本Xnew的数目为2,第二车次的模型更新的新样本Xnew的数目为3,则用于更新的新样本Xnew的数目为5,达到了阈值T,对硬度在线预测模型进行更新,获取更新后的硬度在线预测模型。
其中,对硬度在线预测模型f:x →y更新一次具体为:将用于更新的新样本输入到硬度在线预测模型数据库中,获取更新后的硬度在线预测模型。
其中,该步骤中的统计量的值Qnew满足:Qnew<Qa或Qal<Qnew≤Qah。
108:通过更新后的硬度在线预测模型预测下一车次的新样本Xnew对应的橡胶硬度值,获取橡胶硬度值与实际测得的橡胶硬度值的差值,获得误差;
109:流程结束。
下面以一个具体的试验来验证本发明实施例提供的一种橡胶硬度的测量方法的可行性,详见下文描述:
随机抽取8个样本做实验分析,编号1到8,分别用现有技术和本专利方法计算橡胶的硬度值,并分别得到预测值与真实值误差的平方(Square Error,SE),结果比较见下图2。
从图2中对比分析可以得出,本方法得到的橡胶硬度预测值与真实值的误差平方可以更小,即橡胶硬度测量值更接近真实值。本方法在一定程度上解决了橡胶硬度预测值偏差过大的问题,因而减小了因为硬度测量不准确带来的损失。
综上所述,本发明实施例给出了一种橡胶硬度的测量方法,本方法通过不断的对硬度在线预测模型进行更新,从而得到精确的橡胶硬度值,提高了橡胶的质量,并且通过实际的试验验证了采用本方法可以得到较高的精度,减少了安全隐患。
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种橡胶硬度的测量方法,其特征在于,所述方法包括以下步骤:
(1)在不同置信水平α下设定Q统计量的三个控制限Qas、Qal和Qah,其中满足关系式Qas<Qal<Qah;
(2)将当前车次混炼过程结束后,采集到的终炼胶质量指标中的流变参数与硬度在线预测模型f:x →y自动连接,并且设定流变参数的新样本Xnew;
(3)计算所述新样本Xnew对应的Q统计量的值Qnew;
(4)判断是否满足Qnew>Qah,如果是,执行步骤(5);如果否,执行步骤(6);
(5)与所述硬度在线预测模型f:x →y相连接的警报器发生警报,操作工人查看生产线的工况,然后进行相应的工艺调整,使当前工况恢复到正常生产状况,流程结束;
(6)判断是否满足Qas≤Qnew≤Qal,如果是,执行步骤(7);如果否,执行步骤(9);
(7)则所述新样本Xnew用于所述硬度在线预测模型f:x →y的更新,当用于模型更新的新样本Xnew的数目达到阈值T时,所述硬度在线预测模型f:x →y更新一次;当没有达到阈值T时,不进行更新,直到达到阈值T,获取更新后的硬度在线预测模型,执行步骤(8);
(8)通过所述更新后的硬度在线预测模型预测下一车次的新样本Xnew对应的橡胶硬度值,获取所述橡胶硬度值与实际测得的橡胶硬度值的差值,获得误差;
(9)流程结束。
2.根据权利要求1所述的一种橡胶硬度的测量方法,其特征在于,
式中: (i=1,2,3...)
λj为X的协方差矩阵的特征值,Cα为正态水平在检验水平为α下的临界值,k是主元模型中保留的主元个数,n是全部主元个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100743611A CN102608302A (zh) | 2012-03-20 | 2012-03-20 | 一种橡胶硬度的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012100743611A CN102608302A (zh) | 2012-03-20 | 2012-03-20 | 一种橡胶硬度的测量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102608302A true CN102608302A (zh) | 2012-07-25 |
Family
ID=46525842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012100743611A Pending CN102608302A (zh) | 2012-03-20 | 2012-03-20 | 一种橡胶硬度的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102608302A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102621018A (zh) * | 2012-03-20 | 2012-08-01 | 天津大学 | 一种橡胶硬度的在线自动测量方法 |
CN108491683A (zh) * | 2018-03-26 | 2018-09-04 | 河北工业大学 | 一种青霉素发酵过程故障检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006349495A (ja) * | 2005-06-15 | 2006-12-28 | Central Res Inst Of Electric Power Ind | 高分子材料の劣化解析方法および劣化解析プログラム |
CN101872163A (zh) * | 2010-05-21 | 2010-10-27 | 杭州电子科技大学 | 一种基于递推非线性部分最小二乘的过程监控方法 |
CN101893884A (zh) * | 2010-06-30 | 2010-11-24 | 浙江大学 | 密炼机的橡胶混炼过程中质量指标数据的软测量方法 |
CN102357934A (zh) * | 2011-08-29 | 2012-02-22 | 天津大学 | 一种基于橡胶混炼过程的质量监控软测量方法 |
-
2012
- 2012-03-20 CN CN2012100743611A patent/CN102608302A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006349495A (ja) * | 2005-06-15 | 2006-12-28 | Central Res Inst Of Electric Power Ind | 高分子材料の劣化解析方法および劣化解析プログラム |
CN101872163A (zh) * | 2010-05-21 | 2010-10-27 | 杭州电子科技大学 | 一种基于递推非线性部分最小二乘的过程监控方法 |
CN101893884A (zh) * | 2010-06-30 | 2010-11-24 | 浙江大学 | 密炼机的橡胶混炼过程中质量指标数据的软测量方法 |
CN102357934A (zh) * | 2011-08-29 | 2012-02-22 | 天津大学 | 一种基于橡胶混炼过程的质量监控软测量方法 |
Non-Patent Citations (2)
Title |
---|
宋凯等: "基于递推部分最小二乘自适用质量监控策略及其在橡胶混炼过程中的应用", 《化工学报》 * |
赵刚等: "软测量在橡胶混炼过程建模中的应用", 《橡胶工业》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102621018A (zh) * | 2012-03-20 | 2012-08-01 | 天津大学 | 一种橡胶硬度的在线自动测量方法 |
CN108491683A (zh) * | 2018-03-26 | 2018-09-04 | 河北工业大学 | 一种青霉素发酵过程故障检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109459993B (zh) | 一种流程工业过程在线自适应故障监测与诊断方法 | |
CN102601881B (zh) | 一种橡胶硬度的在线质量监控和预测模型更新的方法 | |
CN102357934B (zh) | 一种基于橡胶混炼过程的质量监控软测量方法 | |
CN101863088B (zh) | 一种橡胶混炼过程中门尼粘度的预报方法 | |
CN108549908B (zh) | 基于多采样概率核主成分模型的化工过程故障检测方法 | |
CN102042848A (zh) | 基于多元混合时序分析的多性能参数加速退化试验产品寿命预测方法 | |
CN107862324B (zh) | 一种基于mwspca的cbr预测模型智能化预警方法 | |
CN103838202A (zh) | 参数控制方法和参数控制系统 | |
CN102357933B (zh) | 一种基于橡胶混炼过程的质量实时监控方法 | |
CN114282398A (zh) | 一种基于大数据的桥梁健康监测系统及方法 | |
CN105095990A (zh) | 一种预测维修方法和装置 | |
CN105334185A (zh) | 基于光谱投影判别的近红外模型维护方法 | |
CN102621953B (zh) | 一种橡胶硬度的在线自动质量监控和预测模型更新的方法 | |
CN104880217A (zh) | 一种基于测量值关联度的故障传感器信息重构方法 | |
CN102608303B (zh) | 一种橡胶硬度的在线测量方法 | |
CN102390096A (zh) | 一种橡胶门尼粘度的在线自动测量方法 | |
CN104615123B (zh) | 基于k近邻的传感器故障隔离方法 | |
CN109978059B (zh) | 一种原油蒸馏过程初馏塔冲塔故障预警方法 | |
CN103389360B (zh) | 基于概率主元回归模型的脱丁烷塔丁烷含量软测量方法 | |
CN102621018A (zh) | 一种橡胶硬度的在线自动测量方法 | |
CN102608302A (zh) | 一种橡胶硬度的测量方法 | |
CN102521698B (zh) | 大型建筑工程质量的关键工序识别与监控方法 | |
CN101872163B (zh) | 一种基于递推非线性部分最小二乘的过程监控方法 | |
CN103279030A (zh) | 基于贝叶斯框架的动态软测量建模方法及装置 | |
CN102567386A (zh) | 在线检测数据库更新方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20120725 |