(57)【特許請求の範囲】
【請求項１】
フレキシブル材料ローラツーローラ加工製造システムモジュールのうちの巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラのリアルタイムの張力に関する方程式を作成するステップであって、前記巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラのリアルタイムの張力に関する方程式はそれぞれ
0006755320.tif 000035
であるステップと、
前記方程式より張力変動の影響要因を抽出し、前記影響要因を連続的な公差解析モデルの入力変数として、連続的な公差解析モデルを作成するステップであって、前記張力変動の影響要因は、前記巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラのそれぞれのローラ半径Ｒ（ｔ）、ローラ回転角速度ω（ｔ）であるステップと、
前記巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラのそれぞれのローラ半径Ｒ（ｔ）、ローラ回転角速度ω（ｔ）を用いて、張力伝達システム方程式を作成するステップであって、前記巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラの張力伝達システム方程式はそれぞれ
0006755320.tif 000036
であるステップと、
逐次回帰解析法、Ｔｏｂｉｔアルゴリズムによって張力伝達システム方程式の構造とパラメータを識別し、Ｊｉ、ａｉ、ｂｉが確定された張力伝達システム状態方程式を取得するステップと、
張力伝達システム状態方程式によって前記巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラのそれぞれの張力値を算出するステップと、を含み、
式の左側のＴｕ（ｔ）は巻き戻しローラと次のローラとの間の張力、Ｍｆｕ（ｔ）は巻き戻しローラの摩擦抵抗モーメント、Ｍｕ（ｔ）は巻き戻しローラに印加された電磁トルク、
0006755320.tif 000037
は巻き戻しローラの角速度変化率、
0006755320.tif 000038
は巻き戻しローラの角速度、Ｊ０は巻き戻しローラの慣性モーメントを示し、固定値であり、ρはフレキシブル薄膜材料の密度、ｂはフレキシブル薄膜材料の幅、Ｒ０は材料を巻いていない時の巻き戻しローラの半径、Ｒｕ（ｔ）は巻き戻しローラに巻いた材料のリアルタイム半径、ｈはフレキシブル薄膜の厚さ、Ｔｗ（ｔ）は巻取りローラと前のローラとの間の張力、
0006755320.tif 000039
は巻取りローラのリアルタイム角速度、
0006755320.tif 000040
は巻取りローラの角速度変化率、Ｒｗ（ｔ）は巻取りローラに巻いた材料のリアルタイム半径、Ｍｗ（ｔ）は巻取りローラに印加された電磁トルク、Ｍｆｗ（ｔ）は巻取りローラに印加された摩擦抵抗モーメント、Ｔｉはローラの張力値、Ｔｉ−１は前記ローラの前のローラの張力値、Ｊｉ（ｔ）はローラの等価慣性モーメント、
0006755320.tif 000041
はローラのリアルタイム角速度、
0006755320.tif 000042
はローラの角速度変化率、ａｉはモータから駆動ローラに伝達されるトルクの百分率、Ｍｍｏｔｏｒ＿ｉはモータの電磁トルク、ｂｉは摩擦係数、Ｒｉは現ローラに巻いた材料のリアルタイム半径を示し、
式の左側Ｔｉ＊では、ｉが位置番号を示し、Ｒｉ（ｔ）はｉ番目のローラ半径、Ｒｉはｉ番目のローラが材料を巻いていない時のローラ半径、Ｊｉはｉ番目のローラの慣性モーメント、Ｗｉはランダムバリエーションを示し、ｉ番目のローラの電磁トルクは
0006755320.tif 000043
、ｉ番目のローラの摩擦抵抗モーメントは
0006755320.tif 000044
であり、ａｉはモータから駆動ローラに伝達されるトルクの百分率、ｂｉは摩擦係数であり、Ｔ０（Ｒ１（ｔ））は巻き戻しローラの初期張力であり、ローラ半径が変化すると、初期張力も変化し、
0006755320.tif 000045
はｉ番目のローラの角速度変化率、
0006755320.tif 000046
はｉ番目のローラの角速度、Ｍｍｏｔｏｒ＿１は１番目のモータの電磁トルク、ａ１は１番目のモータから駆動ローラに伝達されるトルクの百分率、ｂ１は１番目のローラの摩擦係数、
0006755320.tif 000047
は１番目のローラの角速度、Ｍｍｏｔｏｒ＿ｉはｉ番目のモータの電磁トルク、Ｔｉ−１はいｉ−１番目のローラの張力値、ρはフレキシブル薄膜材料の密度、ｂはフレキシブル薄膜材料の幅、Ｒ０は材料を巻いていない時の巻き戻しローラの半径、ｈはフレキシブル薄膜の厚さを示す、ことを特徴とするフレキシブル材料のローラツーローラ加工過程の張力予測用の連続的な公差解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本願は２０１７年０９月１４日に中国特許庁に提出し、出願番号が２０１７１０８２８５８０．７、発明名称が「フレキシブル材料のローラツーローラ加工過程の張力予測用のＳｏＶ方法」の中国特許出願の優先権を主張し、その全内容は引用により本願に組み込まれる。
【０００２】
本願の実施例は加工システムの各位置のローラの張力予測の技術分野に関し、特にフレキシブル材料のローラツーローラ加工過程の張力予測用の連続的な公差解析（英語：Ｓｔｒｅａｍ  ｏｆ  Ｖａｒｉａｔｉｏｎ、略称：ＳｏＶ）方法に関する。
【背景技術】
【０００３】
近年、フレキシブル薄膜をベース材料とするウェアラブルセンサー、有機発光ダイオード（英語：Ｏｒｇａｎｉｃ  Ｌｉｇｈｔ−Ｅｍｉｔｔｉｎｇ  Ｄｉｏｄｅ、略称：ＯＬＥＤ）、薄膜太陽電池についての技術研究は飛躍的な進歩を遂げ、その産業化、量産化が課題として挙げられている。しかし、フレキシブル薄膜材料１５は異方性材料であり、変形には多様性や不確定性を有し、しわ、層間スリップ、破損などの品質問題が生じやすい。ローラツーローラ連続製造システムは自動化程度が高く、生産効率が高く、人的要因による生産過程への影響が小さいなどの利点を有するため、フレキシブル材料の加工製造に最適である。
【０００４】
フレキシブル材料のローラツーローラ（Ｒ２Ｒ）製造は通常、数個又は数十個のローラによって連続生産システムを構成してフレキシブル薄膜にプレス、印刷などの加工を行い、フレキシブル薄膜材料の加工過程では、ローラ間の張力変動によって材料の変形を引き起こし、さらに加工品質問題を招いてしまい、Ｒ２Ｒ加工過程の張力変化は複数の要因が存在し、公差が累積するという特徴を有する。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記技術的問題を解決するために、本発明は、Ｒ２Ｒ加工システムの各位置のローラの張力予測に用いられ、張力予測によって加工品質の評価を間接的に実現するフレキシブル材料のローラツーローラ加工過程の張力予測用の連続的な公差解析（Ｓｔｒｅａｍ  ｏｆ  Ｖａｒｉａｔｉｏｎ）方法を提供することを目的する。
【課題を解決するための手段】
【０００６】
本発明の目的は、以下の技術案によって実現される。
【０００７】
フレキシブル材料のローラツーローラ加工過程の張力予測用の連続的な公差解析方法であって、
Ａフレキシブル材料のローラツーローラ加工製造システムモジュールを動的解析するステップと、
Ｂシステムモジュールのうちの巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラの動的方程式を作成するステップと、
Ｃ巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラの動的方程式の解析によって張力変動影響要因を抽出し、前記影響要因を連続的な公差解析モデルの入力変数とするステップと、
Ｄそれぞれ巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラの張力伝達システム状態方程式を作成するステップと、
Ｅ逐次回帰解析法、Ｔｏｂｉｔアルゴリズムによって張力伝達状態方程式の構造とパラメータを識別し、完全な張力伝達システム状態方程式を取得するステップと、
Ｆ張力伝達システム状態方程式によって各位置のローラの張力値を算出し、前記張力値に基づきローラツーローラ製造システムの加工過程品質評価を実現するステップと、を含む。
【０００８】
従来技術に比べて、本発明の１つ又は複数の実施例は以下の利点を有する。
【０００９】
１、最も少ないセンサによって最も効果的な監視能力を実現し、それにより製造加工コストを最大限に削減させる。
２、ローラツーローラ製造システムの加工過程オンライン品質評価の難題を克服する。
３、張力予測の平均相対誤差はわずか０．３９％であり、フレキシブル材料のローラツーローラ高精度予測製造に対して重要な意味がある。
【図面の簡単な説明】
【００１０】
本願の実施例又は従来技術の技術案をより明確に説明するために、以下、本願の実施例又は従来技術の説明に必要な図面を簡単に説明し、明らかなように、後述する図面は本願のいくつかの例示的な実施例に過ぎず、当業者は創造的な努力をせずにこれらの図面に基づきほかの図面を想到し得る。
【００１１】
【図１】図１はフレキシブル材料のローラツーローラ加工過程の張力予測用の連続的な公差解析方法のブロック図である。
【図２】図２は逐次回帰解析方法のフローチャートである。
【図３】図３はＴｏｂｉｔ最尤推定アルゴリズムのフローチャートである。
【発明を実施するための形態】
【００１２】
本発明の目的、技術案及び利点をより明確にするために、以下、実施例及び図面を参照して本発明をさらに詳細に説明する。
【００１３】
図１に示すように、フレキシブル材料のローラツーローラ加工過程の張力予測用の連続的な公差解析方法であり、
フレキシブル材料ローラツーローラ加工製造システムモジュールを動的解析するステップと、
システムモジュールのうちの巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラの動的方程式を作成するステップと、
巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラの動的方程式の解析によって張力変動影響要因を抽出し、前記影響要因を連続的な公差解析モデルの入力変数とするステップと、
それぞれ巻き戻しローラ、巻取りローラ、駆動ローラ及びガイドローラの張力伝達システム状態方程式を作成するステップと、
逐次回帰解析法、Ｔｏｂｉｔアルゴリズムによって張力伝達状態方程式の構造とパラメータを識別し、完全な張力伝達システム状態方程式を取得するステップと、
張力伝達システム状態方程式によって各位置のローラの張力値を算出し、前記張力値に基づきローラツーローラ製造システムの加工過程品質評価を実現するステップと、を含む。
【００１４】
上記巻き戻しローラの動的方程式、巻取りローラの動的方程式、駆動ローラの動的方程式、ガイドローラの動的方程式はそれぞれ、式（１）、式（２）、式（３）、式（１）である。
0006755320.tif 000002
ここで、式の左側のＴｕ（ｔ）は巻き戻しローラと次のローラとの間の張力、Ｍｆｕ（ｔ）は巻き戻しローラの摩擦抵抗モーメント、Ｍｕ（ｔ）は巻き戻しローラに印加された電磁トルク、
0006755320.tif 000003
は巻き戻しローラの角速度変化率、
0006755320.tif 000004
は巻き戻しローラの角速度、Ｊ０は巻き戻しローラの慣性モーメントを示し、固定値であり、ρはフレキシブル薄膜材料の密度、ｂはフレキシブル薄膜材料の幅、Ｒ０は材料を巻いていない時の空き軸の半径、Ｒｕ（ｔ）は巻き戻しローラに巻いた材料のリアルタイム半径、ｈはフレキシブル薄膜の厚さ、Ｔｗ（ｔ）は巻取りローラと前のローラとの間の張力、
0006755320.tif 000005
は巻取りローラのリアルタイム角速度、
0006755320.tif 000006
は巻取りローラの角速度変化率、Ｒｗ（ｔ）は巻取りローラに巻いた材料のリアルタイム半径、Ｍｗ（ｔ）は巻取りローラに印加された電磁トルク、Ｍｆｗ（ｔ）は巻取りローラに印加された摩擦抵抗モーメント、Ｔｉは現ローラの張力値、Ｔｉ−１は現ローラの前のローラの張力値、Ｊｉ（ｔ）は現ローラの等価慣性モーメント、
0006755320.tif 000007
は現ローラのリアルタイム角速度、
0006755320.tif 000008
は現ローラの角速度変化率、ａｉはモータから駆動ローラに伝達されるトルクの百分率、Ｍｍｏｔｏｒ＿ｉはモータの電磁トルク、ｂｉは摩擦係数、Ｒｉは現ローラに巻いた材料のリアルタイム半径を示す。
【００１５】
ローラ半径Ｒ（ｔ）、ローラ回転角速度ω（ｔ）などを含む、張力変動を影響する主要要因をモデルの入力変数とし、動的方程式に基づき、連続的な公差解析（ＳｏＶ）理論を導入して巻き戻しローラ、巻取りローラ、駆動ローラ、ガイドローラの初期張力伝達システム状態方程式を作成し、
巻き戻しローラの初期張力伝達システム状態方程式は、
0006755320.tif 000009
巻取りローラの初期張力伝達システム状態方程式は、
0006755320.tif 000010
、
駆動ローラの初期張力伝達システム状態方程式は、
0006755320.tif 000011
ガイドローラの初期張力伝達システム状態方程式は、
0006755320.tif 000012
であり、
ここで、式の左側Ｔｉ＊はｉで位置番号を示し、Ｒｉ（ｔ）はｉ番目のローラ半径、Ｒｉはｉ番目のローラの半径、Ｊｉはｉ番目のローラの慣性モーメント、Ｗｉはランダムバリエーションを示し、ｉ番目のローラの電磁トルクは、
0006755320.tif 000013
ｉ番目のローラの摩擦抵抗モーメントは
0006755320.tif 000014
であり、ａｉはモータから駆動ローラに伝達されるトルクの百分率、ｂｉは摩擦係数である。Ｔ０（Ｒ１（ｔ））は巻き戻しローラの初期張力であり、ローラ半径が変化すると、初期張力も変化する。
0006755320.tif 000015
はｉ番目の角速度変化率、
0006755320.tif 000016
はｉ番目のローラの角速度、Ｍｍｏｔｏｒ＿１は１番目のモータの電磁トルク、ａ１は１番目のローラのモータから駆動ローラに伝達されるトルクの百分率、ｂ１は１番目のローラの摩擦係数、
0006755320.tif 000017
は１番目のローラの角速度、Ｍｍｏｔｏｒ＿ｉはｉ番目のローラのモータの電磁トルク、Ｔｉ−１は現ローラの前のローラの張力値、ρはフレキシブル薄膜材料の密度、ｂはフレキシブル薄膜材料の幅、Ｒ０は材料を巻いていない時の空き軸の半径、ｈはフレキシブル薄膜の厚さを示す。
【００１６】
初期張力伝達システムの状態方程式を作成した後、逐次回帰解析法によってモデル構造を画定する。図２に示すように、逐次回帰解析法による変数選択過程は、回帰方程式から不顕著な変数を除去するステップと、新たな変数を回帰方程式に導入するステップと、を含む。
ｙを従属変数、ｘ１、ｘ２、Ｌ、ｘｍをすべての独立変数、ｙｉ、ｘｉ１、ｘｉ２、Ｌ、ｘｉｍ（ｉ＝１，２，Ｌ，ｎ）を独立に抽出したｎ組のサンプルとする。独立変数がモデルに選択される顕著性レベルをα１、モデルから除去される顕著性レベルをα２（０＜α１≦α２＜１）とする。逐次スクリーニング法の一般なステップは以下の通りである。
【００１７】
ｌ回のスクリーニングを行い、且つｐ個の独立変数を選択し、対応した残差平方和を
0006755320.tif 000018
、共分散行列を
0006755320.tif 000019
とすると、
第ｌ＋１回のスクリーニング過程は、以下のステップを含む。
【００１８】
ａ）回帰方程式から変数を除去できるか否かを考慮し、独立変数の貢献を算出する。
0006755320.tif 000020
【００１９】
ｂ）選択された独立変数が顕著であるか否かを検証する。
0006755320.tif 000021
【００２０】
ｃ）新たな変数をさらに導入できるか否かを考慮し、回帰方程式にない変数の偏回帰平方和を算出し、
0006755320.tif 000022
【００２１】
0006755320.tif 000023
【００２２】
その後、モデル訓練過程では、逐次回帰解析方法によって行う。
【００２３】
モデル構造を確認した後、その未知パラメータを求める必要がある。加工過程で張力は多くの要因に影響され、打切りデータに属する。従って、Ｔｏｂｉｔ回帰モデル最尤推定方法によってモデルの未知パラメータを求め、実現方法は以下の通りである。
【００２４】
Ｙは制限付き従属変数であり、下限０を有するとする。
0006755320.tif 000024
ここで、Ｘ１〜Ｘｍは従属変数、εは打ち切ったランダム変数、β０〜βｍは係数ベクトルを示す。
φ（ｘ）は標準密度関数、Φ（ｘ）は標準正規分布関数を示す、対応した確率分布関数は、
0006755320.tif 000025
であり、
該モデルに対して、ＭＬＥ推定尤度関数を用いて、ニュートン反復法によって尤度関数が最大値であるオイラー方程式を求め、プロセスは以下の通りである。
【００２５】
１）（ａ０、ａ１、...、ａｍ、ａ）を
0006755320.tif 000026
の推定値とする。
【００２６】
２）サンプルの尤度関数を求める。
0006755320.tif 000027
0006755320.tif 000028
0006755320.tif 000029
【００２７】
３）尤度関数に対して対数を取る。
0006755320.tif 000030
【００２８】
４）一次条件式を求める。
0006755320.tif 000031
0006755320.tif 000032
【００２９】
一次条件式に基づき、初期値を代入し、ニュートン反復法で計算し、Δａの値の変化が非常に小さくなるまで反復し、取得する推定値はＴｏｂｉｔモデルの推定値
0006755320.tif 000033
であり、アルゴリズムのフローチャートは図３に示される。
Ｔｏｂｉｔに基づく連続的な公差解析モデル構造を決定した後、上記最尤推定（ＭＬＥ）を用いて、未知パラメータＪｉ、ａｉ及びｂｉを推定する。ここまで、フレキシブル材料ローラツーローラ製造システムの張力予測連続的な公差解析モデルの作成は完了する。
【００３０】
フレキシブル材料ローラツーローラシステムの各位置のローラの張力値を算出するために、巻き戻し過程のデータをランダムに収集して入力として上記モデル訓練を行い、モデリングを実現する。２秒ごとに収集し、主に初期張力値と最後位置の張力値、巻き戻し速度、ガイドローラ速度、駆動ローラ速度、ローラ半径、初期張力、巻き戻しモータトルク、駆動ローラモータトルクを含む。すなわち、
0006755320.tif 000034
を必要とせずにパラメータＪ１、ａｑ、ｂ１、Ｊ２、ｂ２、Ｗ１、Ｗ２の推定を行い、限られたセンサで各位置のローラの張力を推定でき、それを根拠にローラツーローラ製造システムの加工過程の品質評価を実現する。
【００３１】
以上、本発明の実施形態を説明したが、本発明を理解するためのものに過ぎず、本発明を限定するものではない。当業者は、本発明の開示する精神及び範囲を逸脱せずに、実施形態及び詳細に様々な変更や変化を行うことができるが、本発明の保護範囲は添付特許請求の範囲に定められる。
【図１】
0006755320.tif 000048
【図２】
0006755320.tif 000049
【図３】
0006755320.tif 000050