CN107729287B - 一种用于柔性材料卷对卷加工过程张力预测的sov方法 - Google Patents
一种用于柔性材料卷对卷加工过程张力预测的sov方法 Download PDFInfo
- Publication number
- CN107729287B CN107729287B CN201710828580.7A CN201710828580A CN107729287B CN 107729287 B CN107729287 B CN 107729287B CN 201710828580 A CN201710828580 A CN 201710828580A CN 107729287 B CN107729287 B CN 107729287B
- Authority
- CN
- China
- Prior art keywords
- roll
- roller
- tension
- sov
- unwinding
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Controlling Rewinding, Feeding, Winding, Or Abnormalities Of Webs (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
本发明公开了一种用于柔性材料卷对卷加工过程张力预测的SOV方法,所述方法包括以下步骤:对柔性材料卷对卷加工制造***模块进行动力学分析;建立***模块中放卷辊、收卷辊、驱动辊和导向辊的动力学方程;通过放卷辊、收卷辊、驱动辊和导向辊的动力学方程分析提取出张力波动影响因素,并将所述影响因素作为SOV模型的输入变量;分别建立放卷辊、收卷辊、驱动辊和导向辊的张力传递***状态方程;通过逐步回归法、Tobit算法对张力传递状态方程的结构和参数进行识别,得到完整的张力传递***状态方程;通过张力传递***状态方程计算各工位的张力值,并根据所述张力值实现卷对卷制造***加工过程质量评估。
Description
技术领域
本发明涉及加工***各工位张力预测技术领域,尤其涉及一种用于柔性材料卷对卷加工过程张力预测的SOV方法。
背景技术
现有近年来,以柔性薄膜作为衬底材料的可穿戴传感器、OLED、薄膜太阳能电池技术研究取得突破,其产业化、规模化生产需求已经提上日程。但是柔性薄膜材料属于各向异性材料,变形具有多样性和不确定,极易出现褶皱、层间滑移、破损等质量问题。卷对卷连续制造***由于具有自动化程度高,生产效率高,生产过程受人为因素影响小等优点,从而成为柔性材料加工制造的最佳选择。
柔性材料R2R制造通常由几个或数十个辊轴组成连续生产***进行柔性薄膜的压制、印刷等加工,柔性薄膜材料加工过程由于卷辊之间张力波动而引起材料变形,进而造成加工质量问题,考虑到R2R加工过程张力变化具有多因素耦合、误差累积的特点。
发明内容
为解决上述技术问题,本发明的目的是提供一种用于柔性材料卷对卷加工过程张力预测的SOV(Stream Of Variation)方法,该方法用于R2R加工***各工位的张力预测,并通过张力的预测来间接实现加工质量的评估。
本发明的目的通过以下的技术方案来实现:
一种用于柔性材料卷对卷加工过程张力预测的SOV方法,包括:
A对柔性材料卷对卷加工制造***模块进行动力学分析;
B建立***模块中放卷辊、收卷辊、驱动辊和导向辊的动力学方程;
C通过放卷辊、收卷辊、驱动辊和导向辊的动力学方程分析提取出张力波动影响因素,并将所述影响因素作为SOV模型的输入变量;
D分别建立放卷辊、收卷辊、驱动辊和导向辊的张力传递***状态方程;
E通过逐步回归法、Tobit算法对张力传递状态方程的结构和参数进行识别,得到完整的张力传递***状态方程;
F通过张力传递***状态方程计算各工位的张力值,并根据所述张力值实现卷对卷制造***加工过程质量评估。
与现有技术相比,本发明的一个或多个实施例可以具有如下优点:
1、实现了以最少的传感器产生最有效的监测能力,从而最大限度地降低制造加工的成本。
2、克服了卷对卷制造***加工过程在线质量评估的困难。
3、张力预测平均相对误差为0.39%,这对实现柔性材料卷对卷高精度预测性制造具有重要意义。
附图说明
图1是用于柔性材料卷对卷加工过程张力预测的SOV方法框图;
图2是逐步回归分析方法流程图;
图3是Tobit最大似然估计算法流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合实施例及附图对本发明作进一步详细的描述。
如图1所示,为用于柔性材料卷对卷加工过程张力预测的SOV方法,包括以下步骤:
对柔性材料卷对卷加工制造***模块进行动力学分析;
建立***模块中放卷辊、收卷辊、驱动辊和导向辊的动力学方程;
通过放卷辊、收卷辊、驱动辊和导向辊的动力学方程分析提取出张力波动影响因素,并将所述影响因素作为SOV模型的输入变量;
分别建立放卷辊、收卷辊、驱动辊和导向辊的张力传递***状态方程;
通过逐步回归法、Tobit算法对张力传递状态方程的结构和参数进行识别,得到完整的张力传递***状态方程;
通过张力传递***状态方程计算各工位的张力值,并根据所述张力值实现卷对卷制造***加工过程质量评估。
上述放卷辊动力学方程、收卷辊动力学方程、驱动辊动力学方程、导向辊动力学方程分别为:
其中公式左侧T表示当前卷辊的张力值,Ti-1表示上一个卷辊的张力值;R(t)表示卷材的实时半径,M(t)表示摩擦阻力矩,M(t)表示施加在卷辊上的电磁转矩,ω(t)表示卷辊的实时角速度,J(t)表示等效转动惯量;J0为放卷辊轴的转动惯量,为固定值;用ρ表示柔性薄膜材料的密度,b表示宽度,R0表示未缠绕材料时的空轴半径。
将影响张力波动的主要因素作为模型的输入变量,包括卷材半径R(t)、卷辊转动角速度ω(t)等,根据动力学方程,引入误差流(SOV)理论建立了放卷辊、收卷辊、驱动辊、导向辊初始的张力传递***状态方程:
放卷辊初始张力传递***状态方程:
收卷辊初始张力传递***状态方程:
驱动辊初始张力传递***状态方程:
导向辊初始张力传递***状态方程:
其中公式左侧Ti*表示用i表示工位序号,Ri(t)表示料卷的半径,Ri表示第i个工位辊轴的半径,Ji表示第i个工位辊轴的转动惯量,Wi、Vi表示随机误差,电磁转矩Mi=aiMmotor_i,摩擦阻力矩Mfi=biωi(t),ai是从电机传给驱动辊的扭矩的百分比,bi为摩擦系数。T0(R1(t))是放卷辊的初始张力,当料卷半径变化时,初始张力也会变化.
初始的张力传递***状态方程建立完成后,通过逐步回归分析法来实现模型结构确定。其中,逐步回归法选择变量的过程包含两个基本步骤(见图2),一是从回归方程中剔除不显著的变量,二是引入新变量到回归方程中。
设y是因变量,x1,x2,…,xm是所有自变量,yi,xi1,xi2,…,xim(i=1,2,…,n)是独立抽取的n组样本。设自变量被选进模型的显著性水平为α1,被剔除模型的显著性水平为α2,且0<α1≤α2<1。逐步筛选法的一般步骤为:
则第l+1步的筛选过程为:
a)考虑能否从回归方程中剔除变量,计算自变量的贡献:
b)检验已选入的自变量是否显著:
检验变量xk能否被剔除:
检验变量xk能否被引入:
在确认T0(R1(t))、L0(R1(t))的模型结构的时候,根据测量得到的T0、R1(t)以及L0、R1(t),进行回归方程的选择和自变量的选择,然后对回归模型进行求解。首先假定T0(R1(t))、L0(R1(t))均为R1(t)的3次多项式:
然后在模型训练过程中通过逐步回归分析方法完成。
模型结构确认以后,需要对其未知参数进行求解。由于加工过程中张力受多因素的影响,且属于删失数据。因此,采用Tobit回归模型最大似然估计的方法来进行模型未知参数求解,实现方法如下所述:
假设Y是受限因变量,具有下限0。
用φ(x)表示标准密度函数,Φ(x)表示标准正态分布函数,相应的概率分布函数为:
对该模型用MLE估计似然函数,并用牛顿迭代法求解似然函数最大值时的欧拉方程,过程如下:
②求样本的似然函数:
其中,I′i=Y′iα=a0+a1X′1i+a2X′2i+…+amX′mi
Ii=Yiα=a0+a1X1i+a2X2i+…+amXmi
③对似然函数取对数:
④求一阶条件表达式:
k=0,1,2,…,m
基于Tobit的SOV模型结构确定以后,采用上述最大似然估计(MLE)来估计Ji、ai和bi这些未知参数。到此,柔性材料卷对卷制造***张力预测SOV模型构建完成。
为了计算柔性材料卷对卷***各工位张力值,需要随机采集放卷过程的数据作为输入来进行上述的模型训练来完成建模。每隔2秒采集一次,主要包括:初始张力值和最后工位张力值,放卷速度、导向辊速度、驱动辊速度、料卷半径、初始张力、放卷电机转矩、驱动辊电机转矩。即:在不需要T1 *的基础上进行参数J1、a1、b1、J2、b2、W1、W2的估计,实现在有限传感器情况下估计每个工位的张力,并以此为依据实现卷对卷制造***加工过程质量评估。
虽然本发明所揭露的实施方式如上,但所述的内容只是为了便于理解本发明而采用的实施方式,并非用以限定本发明。任何本发明所属技术领域内的技术人员,在不脱离本发明所揭露的精神和范围的前提下,可以在实施的形式上及细节上作任何的修改与变化,但本发明的专利保护范围,仍须以所附的权利要求书所界定的范围为准。
Claims (5)
1.一种用于柔性材料卷对卷加工过程张力预测的SOV方法,其特征在于,所述方法包括以下步骤:
A对柔性材料卷对卷加工制造***模块进行动力学分析;
B建立***模块中放卷辊、收卷辊、驱动辊和导向辊的动力学方程;
C通过放卷辊、收卷辊、驱动辊和导向辊的动力学方程分析提取出张力波动影响因素,并将所述影响因素作为SOV模型的输入变量;
D分别建立放卷辊、收卷辊、驱动辊和导向辊的张力传递***状态方程;
E通过逐步回归法、Tobit算法对张力传递***状态方程的结构和参数进行识别,得到完整的张力传递***状态方程;
F通过张力传递***状态方程计算各工位的张力值,并根据所述张力值实现卷对卷制造***加工过程质量评估;
所述放卷辊、收卷辊、驱动辊和导向辊的动力学方程分别为:
其中公式左侧T表示当前卷辊的张力值,Ti-1表示上一个卷辊的张力值;R(t)表示卷材的实时半径,Mfu(t)表示放卷辊的摩擦阻力矩,Mfw(t)表示收卷辊的摩擦阻力矩,Mu(t)表示施加在放卷辊上的电磁转矩,Mw(t)表示施加在收卷辊上的电磁转矩,ω(t)表示卷辊的实时角速度,J(t)表示等效转动惯量;J0为放卷辊轴的转动惯量,为固定值;用ρ表示柔性薄膜材料的密度,b表示宽度,R0表示未缠绕材料时的空轴半径;
步骤D具体包括:根据动力学方程,引入误差流理论分别建立放卷辊、收卷辊、驱动辊和导向辊初始的张力传递***状态方程,所述方程如下:
其中公式左侧Ti*表示用i表示工位序号,Ri(t)表示料卷的半径,Ri表示第i个工位辊轴的半径,Ji表示第i个工位辊轴的转动惯量,Wi、Vi表示随机误差,电磁转矩Mi=aiMmotor_i,摩擦阻力矩Mfi=biωi(t),ai是从电机传给驱动辊的扭矩的百分比,bi为摩擦系数;T0(R1(t))是放卷辊的初始张力,当料卷半径变化时,初始张力也会变化。
2.如权利要求1所述的用于柔性材料卷对卷加工过程张力预测的SOV方法,其特征在于,所述张力波动影响因素包括卷材半径R(t)、卷辊转动角速度ω(t)。
3.如权利要求1所述的用于柔性材料卷对卷加工过程张力预测的SOV方法,其特征在于,所述步骤E具体包括:通过逐步回归法来实现卷对卷加工制造***模型结构的确定,在模型结构的确定后,采用Tobit回归模型最大似然估计的方法来进行模型未知参数求解。
4.如权利要求3所述的用于柔性材料卷对卷加工过程张力预测的SOV方法,其特征在于,所述逐步回归法选择变量的过程包含:
从回归方程中剔除不显著的变量;
并引入新变量到回归方程中。
5.如权利要求1所述的用于柔性材料卷对卷加工过程张力预测的SOV方法,其特征在于,所述柔性材料卷对卷加工制造***各工位张力值是通过随机采集放卷过程的数据作为输入来进行模型训练完成建模;所述采集的数据包括初始张力值和最后工位张力值,即放卷速度、导向辊速度、驱动辊速度、料卷半径、初始张力、放卷电机转矩、驱动辊电机转矩。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710828580.7A CN107729287B (zh) | 2017-09-14 | 2017-09-14 | 一种用于柔性材料卷对卷加工过程张力预测的sov方法 |
PCT/CN2018/075508 WO2019052117A1 (zh) | 2017-09-14 | 2018-02-07 | 一种用于柔性材料卷对卷加工过程张力预测的误差流方法 |
JP2018536849A JP6755320B2 (ja) | 2017-09-14 | 2018-02-07 | フレキシブル材料のローラツーローラ加工過程の張力予測用の連続的な公差解析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710828580.7A CN107729287B (zh) | 2017-09-14 | 2017-09-14 | 一种用于柔性材料卷对卷加工过程张力预测的sov方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107729287A CN107729287A (zh) | 2018-02-23 |
CN107729287B true CN107729287B (zh) | 2021-05-18 |
Family
ID=61206232
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710828580.7A Active CN107729287B (zh) | 2017-09-14 | 2017-09-14 | 一种用于柔性材料卷对卷加工过程张力预测的sov方法 |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP6755320B2 (zh) |
CN (1) | CN107729287B (zh) |
WO (1) | WO2019052117A1 (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108622707B (zh) * | 2018-03-20 | 2020-12-01 | 武汉科技大学 | 柔性基材卷绕***的分布式跟随误差控制模型构建方法 |
CN110502043B (zh) * | 2018-05-18 | 2024-05-24 | 伊顿公司 | 幅材卷绕机张力控制***中的动态性能和主动阻尼方法 |
CN110704884B (zh) * | 2019-09-30 | 2023-05-12 | 中国人民解放军海军工程大学 | 一种对数正态型通用件备件需求量计算方法及装置 |
CN112561317B (zh) * | 2020-12-15 | 2023-07-28 | 广东工业大学 | 一种柔性材料卷对卷加工控制决策生成方法及装置 |
CN113151989B (zh) * | 2021-04-19 | 2022-10-18 | 山东大学 | 布料处理方法、***及缝制机器人 |
CN115146492B (zh) * | 2022-09-06 | 2022-11-18 | 楚能新能源股份有限公司 | 一种评估卷绕极片张力的方法 |
CN118306859A (zh) * | 2024-06-11 | 2024-07-09 | 成都飞机工业(集团)有限责任公司 | 送料盘电机转速控制方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101306773A (zh) * | 2008-06-30 | 2008-11-19 | 陕西北人印刷机械有限责任公司 | 一种非线性收卷张力控制方法 |
CN103440416A (zh) * | 2013-08-27 | 2013-12-11 | 西北工业大学 | 基于扩展误差流的叶片加工过程误差预测方法 |
CN105887286A (zh) * | 2016-05-31 | 2016-08-24 | 天津工业大学 | 一种碳纤维多层角联织机纱线张力控制方法 |
CN106391725A (zh) * | 2016-09-28 | 2017-02-15 | 燕山大学 | 一种适用于冷轧过程轧制压力随轧制速度变化预报方法 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100582986C (zh) * | 2008-07-11 | 2010-01-20 | 华中科技大学 | 一种卷到卷间歇式柔性基材传输的张力控制方法及装置 |
CN102556737B (zh) * | 2011-12-23 | 2014-11-26 | 深圳众为兴技术股份有限公司 | 一种张力控制方法 |
KR101247184B1 (ko) * | 2011-12-30 | 2013-03-25 | 건국대학교 산학협력단 | 롤투롤 라미네이션 공정에서 박리현상과 리와인딩 구간에서의 슬립현상을 방지하기 위한 운전 장력 설정 방법 |
CN102615923B (zh) * | 2012-04-16 | 2015-01-07 | 姜洪波 | 一种实现卷对卷定位复合的设备和方法 |
CN202758329U (zh) * | 2012-04-16 | 2013-02-27 | 姜洪波 | 一种制备用于柔性触摸屏的触控电极卷材的设备 |
US9557374B2 (en) * | 2014-10-08 | 2017-01-31 | Eastman Kodak Company | Vision-guided alignment system |
CN105775806A (zh) * | 2014-12-24 | 2016-07-20 | 重庆鑫仕达包装设备有限公司 | 一种圆盘式收放卷复合机 |
CN205114638U (zh) * | 2015-10-15 | 2016-03-30 | 苏州伍洲设计包装有限公司 | 一种纸卷收放卷机 |
CN106626532B (zh) * | 2016-11-23 | 2018-11-27 | 福建东南艺术纸品股份有限公司 | 一种采用水性油墨印刷的餐巾纸生产方法 |
-
2017
- 2017-09-14 CN CN201710828580.7A patent/CN107729287B/zh active Active
-
2018
- 2018-02-07 WO PCT/CN2018/075508 patent/WO2019052117A1/zh active Application Filing
- 2018-02-07 JP JP2018536849A patent/JP6755320B2/ja active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101306773A (zh) * | 2008-06-30 | 2008-11-19 | 陕西北人印刷机械有限责任公司 | 一种非线性收卷张力控制方法 |
CN103440416A (zh) * | 2013-08-27 | 2013-12-11 | 西北工业大学 | 基于扩展误差流的叶片加工过程误差预测方法 |
CN105887286A (zh) * | 2016-05-31 | 2016-08-24 | 天津工业大学 | 一种碳纤维多层角联织机纱线张力控制方法 |
CN106391725A (zh) * | 2016-09-28 | 2017-02-15 | 燕山大学 | 一种适用于冷轧过程轧制压力随轧制速度变化预报方法 |
Non-Patent Citations (3)
Title |
---|
Output Regulation of Nonlinear Systems With Application to Roll-to-Roll Manufacturing Systems;Pramod R.Raul 等;《IEEE/ASME Transactions on Mechatronics》;20141120;第20卷(第3期);1089-1098 * |
制造过程质量控制中误差流理论的研究;罗振璧 等;《机械工程学报》;19950715(第4期);62-69 * |
柔性材料加工变形影响因素提取层次分析方法;邓耀华 等;《机械工程学报》;20160630;第52卷(第11期);161-169 * |
Also Published As
Publication number | Publication date |
---|---|
CN107729287A (zh) | 2018-02-23 |
JP2019535609A (ja) | 2019-12-12 |
JP6755320B2 (ja) | 2020-09-16 |
WO2019052117A1 (zh) | 2019-03-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107729287B (zh) | 一种用于柔性材料卷对卷加工过程张力预测的sov方法 | |
JP2019535609A5 (zh) | ||
CN111829782B (zh) | 一种基于自适应流形嵌入动态分布对齐的故障诊断方法 | |
TWI267012B (en) | Quality prognostics system and method for manufacturing processes | |
CN111209708B (zh) | 基于机器学习的桩-土相互作用预测分析方法 | |
CN102674052B (zh) | 卷绕装置及卷绕控制方法 | |
CN107885928A (zh) | 考虑测量误差的步进应力加速性能退化可靠性分析方法 | |
CN107464023A (zh) | 生命周期预测模型的在线学习方法和设备 | |
CN102795497B (zh) | 卷筒物卷取装置及其控制方法 | |
JP6080297B2 (ja) | 塗工乾燥シミュレーション装置および塗工乾燥装置 | |
CN111709577B (zh) | 基于长程相关性gan-lstm的rul预测方法 | |
CN111579993A (zh) | 一种基于卷积神经网络的锂电池容量在线估计方法 | |
CN108492026B (zh) | 一种基于集成正交成分最优化回归分析的软测量方法 | |
Jia et al. | State of health prediction of lithium-ion batteries based on bidirectional gated recurrent unit and transformer | |
CN110728314B (zh) | 一种大规模免调度***活跃用户检测方法 | |
Zhou et al. | A wavelet neural network informed by time-domain signal preprocessing for bearing remaining useful life prediction | |
CN114770607A (zh) | 一种基于大数据的机器人健康监控方法及*** | |
Ma et al. | A missing manufacturing process data imputation framework for nonlinear dynamic soft sensor modeling and its application | |
Miao et al. | Fault diagnosis of wheeled robot based on prior knowledge and spatial-temporal difference graph convolutional network | |
Li et al. | Intelligent fault diagnosis of rotating machinery based on deep recurrent neural network | |
CN112527547A (zh) | 一种基于自动卷积神经网络的机械智能故障预测方法 | |
CN116861256A (zh) | 一种固废焚烧过程的炉温预测方法、***、设备及介质 | |
JP2005242803A (ja) | 機械の性能推定器、性能推定方法及び性能推定プログラム | |
CN116565851A (zh) | 基于聚类算法的风电场功率预测方法、装置、设备和介质 | |
CN116960962A (zh) | 一种跨区域数据融合的中长期区域负荷预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |