JP5531350B2

JP5531350B2 - 張力測定装置

Info

Publication number: JP5531350B2
Application number: JP2011277763A
Authority: JP
Inventors: 広敏在原; 徹岡田
Original assignee: Kobe Steel Ltd
Current assignee: Kobe Steel Ltd
Priority date: 2011-12-20
Filing date: 2011-12-20
Publication date: 2014-06-25
Anticipated expiration: 2031-12-20
Also published as: JP2013130393A

Description

本発明は、帯状体の幅方向の張力分布を測定する張力測定装置に関する。

薄鋼板等の帯状体を通板して、圧延、矯正、焼鈍、表面処理、脱脂、スリッティング等の各種処理を連続的に行うプロセスラインでは、帯状体の走行時の蛇行を防止するためや、各種処理を良好に行うために、帯状体が走行する長手方向に張力を付与することが多い。また、これらのプロセスラインに通板される帯状体には、耳波、中伸び、片伸びや、これらが複合した複合伸び等の幅方向での不均一歪みが存在することがある。これらの幅方向の不均一歪みは、板製品の平坦度不良となるのみでなく、スムーズな通板を阻害することもあるので、ロールベンディング装置やクーラントのゾーンコントロール装置等によって、不均一歪みを除去または低減するように平坦度制御が行われている。

上述した帯状体の長手方向に張力が付与されるプロセスラインでは、幅方向の不均一歪みが潜在化して、見かけ上平坦になることがある。このため、平坦度制御を行うために、帯状体がロール等によって支持される長手方向の２箇所の支持部位間で、帯状体に振動荷重や静荷重を負荷して、そのときの帯状体の幅方向の変位分布を計測することにより、幅方向の張力分布を測定し、不均一歪みを間接的に検出する技術が開発されている（例えば、特許文献１、２参照）。

特許文献１、２に記載されたものは、いずれも幅方向の不均一歪みが全て潜在化する場合を想定しており、特許文献１では、振動荷重による変位分布として、一次共振モードと二次共振モードを板幅方向の複数点で計測し、計測した変位分布の形状を、予め区分した張力分布パターンのいずれかに決定し、決定した張力分布パターンに従って、別途に計測した全張力を板幅方向の各位置に割り振ることで、張力分布を測定するようにしている。また、特許文献２では、幅方向の複数点で振動モードを計測し、各振動モードの周波数と振幅の正規化値から、幅方向の張力分布を測定するようにしている。

特許文献１に記載されたものは、張力分布を予め設定したいずれかの張力分布パターンに当てはめるので、実際に生じる複雑な張力分布を精度よく測定できない問題があった。また、特許文献２に記載されたものは、加振方法や加振位置によって振動モード形状がばらつくので、測定精度が悪化する問題があった。

本発明者らは、これらの問題に対して、どのような張力分布であっても精度よく測定できるように、荷重が負荷された帯状体の幅方向の複数の測定点における変位を計測する変位計測手段と、帯状体について、複数の測定点に対応する節点と、該節点に接合された節点に作用する張力を模擬したばねとを有する２次元の多質点系モデルにモデル化するモデル化手段と、計測された測定点の変位量分布と、モデルの節点の変位量分布とが一致するようなばねのばね定数を算出するばね定数算出手段と、算出されたばね定数に基づいて張力を算出する張力算出手段とを有する張力測定装置を先に提案している（特願２０１０−２７４３４５）。この張力測定装置も、帯状体の幅方向の不均一歪みが全て潜在化する場合を想定している。

特開平６−４３０５１号公報特開平７−２１８３５８号公報

上述した薄鋼板等のプロセスラインには、帯状体に付与できる張力が低い値に制約されるものや、例えば、矯正ラインの上流側のように、帯状体に大きな幅方向の不均一歪みが存在するものがあり、不均一歪みが全て潜在化せずに一部が顕在することがある。このように幅方向の不均一歪みが一部顕在する場合は、特許文献１、２に記載された張力分布測定方法や特願２０１０−２７４３４５で提案した張力測定装置は、いずれも幅方向の不均一歪みが全て潜在化することを前提としているので、張力分布を精度よく測定できない問題がある。

そこで、本発明の課題は、帯状体に幅方向の不均一歪みが一部顕在しても、張力分布を精度よく測定できるようにすることである。

上記の課題を解決するために、本発明は、長手方向に張力を付与された帯状体の幅方向の張力分布を、長手方向の２箇所の部位で支持された支持部位間で測定する張力測定装置において、前記２箇所の支持部位間で前記帯状体に振動荷重を負荷する振動荷重負荷手段と、この振動荷重負荷手段によって生じる振動変位を帯状体の幅方向の複数の測定点で計測する変位計測手段とを備え、前記帯状体について、前記幅方向の複数の各測定点に対応する各節点に作用する張力を模擬する直線ばねと、該各節点における幅方向の曲げ剛性を模擬する回転ばねとを含む２次元の多質点系モデルにモデル化するモデル化手段と、前記多質点系モデルの固有値解析から得られる前記各節点での固有振動数および振動モードが、前記計測された各測定点での振動変位から得られる固有振動数および振動モードと一致するような前記直線ばねおよび回転ばねの各ばね定数を算出するばね定数算出手段と、算出された前記直線ばねおよび回転ばねの各ばね定数を、それぞれ前記各測定点における張力値および曲げ剛性値に換算するばね定数換算手段とを有し、換算された前記各測定点における張力値から前記帯状体の幅方向の張力分布を求める構成を採用した。

すなわち、２箇所の支持部位間で帯状体に振動荷重を負荷する振動荷重負荷手段と、振動荷重負荷手段によって生じる振動変位を帯状体の幅方向の複数の測定点で計測する変位計測手段とを備え、帯状体について、幅方向の複数の各測定点に対応する各節点に作用する張力を模擬する直線ばねと、該各節点における幅方向の曲げ剛性を模擬する回転ばねとを含む２次元の多質点系モデルにモデル化するモデル化手段と、多質点系モデルの固有値解析から得られる各節点での固有振動数および振動モードが、計測された各測定点での振動変位から得られる固有振動数および振動モードと一致するような直線ばねおよび回転ばねの各ばね定数を算出するばね定数算出手段と、算出された直線ばねおよび回転ばねの各ばね定数を、それぞれ各測定点における張力値および曲げ剛性値に換算するばね定数換算手段とを有し、換算された各測定点における張力値から帯状体の幅方向の張力分布を求めることにより、不均一歪みの一部顕在によって幅方向の曲げ剛性が変化する帯状体であっても、張力を模擬する直線ばねのばね定数と、幅方向の曲げ剛性を模擬する回転ばねのばね定数とを未知数として算出し、張力分布を精度よく測定できるようにした。また、この張力測定装置は、幅方向の曲げ剛性分布も併せて求めることができるので、不均一歪みの顕在形態の推定にも寄与することができる。

前記ばね定数算出手段は、前記固有値解析から得られる前記各節点での固有振動数および振動モードと、前記各測定点での振動変位から得られる固有振動数および振動モードとの各差の二乗和が最小となるような前記直線ばねおよび回転ばねの各ばね定数を、繰り返し計算を用いた最小二乗法によって算出するものとすることができる。

前記ばね定数算出手段が、前記帯状体の振動モードに対する前記直線ばねおよび回転ばねの各ばね定数の係数行列の擬似逆行列を計算し、前記固有値解析から得られる前記各節点での固有振動数および振動モードと、前記各測定点での振動変位から得られる固有振動数および振動モードとが一致するような前記直線ばねおよび回転ばねの各ばね定数を、前記擬似逆行列を用いた最小二乗法によって算出するものとすることもできる。

前記計算された擬似逆行列から特異値分解により誤差成分を除去し、この誤差成分を除去した擬似逆行列を前記最小二乗法による各ばね定数の算出に用いることにより、直線ばねおよび回転ばねの各ばね定数を求める演算を簡単にすることができるとともに、より精度よく張力分布を測定することができる。

本発明に係る張力測定装置は、２箇所の支持部位間で帯状体に振動荷重を負荷する振動荷重負荷手段と、振動荷重負荷手段によって生じる振動変位を帯状体の幅方向の複数の測定点で計測する変位計測手段とを備え、帯状体について、幅方向の複数の各測定点に対応する各節点に作用する張力を模擬する直線ばねと、該各節点における幅方向の曲げ剛性を模擬する回転ばねとを含む２次元の多質点系モデルにモデル化するモデル化手段と、多質点系モデルの固有値解析から得られる各節点での固有振動数および振動モードが、計測された各測定点での振動変位から得られる固有振動数および振動モードと一致するような直線ばねおよび回転ばねの各ばね定数を算出するばね定数算出手段と、算出された直線ばねおよび回転ばねの各ばね定数を、それぞれ各測定点における張力値および曲げ剛性値に換算するばね定数換算手段とを有し、換算された各測定点における張力値から帯状体の幅方向の張力分布を求めるようにしたので、不均一歪みの一部顕在によって幅方向の曲げ剛性が変化する帯状体であっても、張力分布を精度よく測定することができる。

帯状体の張力分布測定装置の実施形態を示す構成図図１のモデル化部で帯状体をモデル化する２次元多質点系モデルの概念図図１の張力分布測定装置で張力分布を測定する第１の実施形態を示すフローチャート図１の張力分布測定装置で張力分布を測定する第２の実施形態を示すフローチャート図４での特異値分解結果を示すグラフ図４での特異値分解の張力分布測定精度に対する効果を示すグラフ（ａ）、（ｂ）、（ｃ）は、それぞれ実施例での薄鋼板の張力分布測定結果をＦＥＭ解析結果と対比して示すグラフ（（ａ）、（ｂ）、（ｃ）は、それぞれ実施例でのアルミニウム薄板の張力分布測定結果をＦＥＭ解析結果と対比して示すグラフ（ａ）、（ｂ）、（ｃ）は、それぞれ実施例での薄鋼板の断面二次モーメント分布の算出結果を示すグラフ（ａ）、（ｂ）、（ｃ）は、それぞれ実施例でのアルミニウム薄板の断面二次モーメント分布の算出結果を示すグラフ

以下、図面に基づき、本発明の実施形態を説明する。この張力測定装置は、図１に示すように、走行する長手方向に張力を付与された帯状体１の張力分布を、長手方向の２箇所の部位で支持ロール２ａ、２ｂによって支持された支持部位間で測定するものであり、支持部位間の帯状体１に振動荷重を負荷する振動荷重負荷装置３と、振動荷重負荷装置３によって生じる振動変位を、帯状体１の幅方向の複数の測定点１ａで計測する非接触式の変位計４と、各変位計４の出力に基づいて、帯状体１の張力分布を演算する演算装置５とからなる。帯状体１の各測定点１ａは、支持部位間の長手方向の中間位置に等間隔で設定されている。

前記振動荷重負荷装置３は、帯状体１に空気を間歇的に噴射して振動荷重を負荷するものであるが、振動荷重負荷装置３は、水や油等の液体を間歇的に噴射するものや、磁力、静電力、電磁誘導による渦電流、音波等によって、振動荷重を負荷するものとすることもできる。また、帯状体１の１点を打撃する装置や、支持ロール２ａ、２ｂのいずれかを加振する装置とすることもできる。

前記変位計４は光反射式のレーザ変位計とされている。帯状体１が導電性を有するものである場合は、帯状体１に生じさせた渦電流の大きさを検出する渦電流式変位計や、帯状体１とセンサヘッド間の静電容量を検出する静電容量式変位計等とすることもできる。また、図１では、便宜上、変位計４を幅方向に等間隔で５台配置するように図示しているが、変位計４の配置台数（測定点１ａの数）は任意に設定することができ、幅方向での配置間隔も、例えば、幅端部を密に、幅中央部を粗くするように不等間隔で配置してもよい。さらに、一部の変位計４を幅方向に移動可能としてもよい。

前記演算装置５は、各変位計４で計測された各測定点１ａの振動変位から固有振動数と振動モードを算出する振動特性算出部５ａと、帯状体１を後述する２次元多質点系モデルにモデル化するモデル化部５ｂと、２次元多質点系モデルの固有値を解析する固有値解析部５ｃと、振動特性算出部５ａで算出された固有振動数および振動モードから２次元多質点系モデルの後述する直線ばねおよび回転ばねの各ばね定数を算出するばね定数算出部５ｄと、算出された直線ばねおよび回転ばねの各ばね定数を各測定点１ａにおける張力値および曲げ剛性値に変換するばね定数変換部５ｅと、変換された各測定点１ａの張力値から帯状体１の幅方向の張力分布を算出する張力分布算出部５ｆとで構成されている。

図２は、前記モデル化部５ｂで帯状体１をモデル化する２次元多質点系モデルを示す。この２次元多質点系モデルは、支持部位間の帯状体１について、振動変位の各測定点１ａに対応する各節点１１を、離間する固定面１２に垂直に直線ばね１３で連結し、隣接する節点１１をリンク１４で結合して、隣接するリンク１４同士を回転ばね１５で連結したものであり、各リンク１４はその幅方向部位での質量ｍと慣性モーメントＪを有する。

このように、前記２次元多質点系モデルは、帯状体１の張力の大小と固有振動数の大小との間に相関があることに着目し、振動荷重負荷装置３により加振された帯状体１の各計測点１ａにおける振動変位と、各直線ばね１３の各節点１１における振動変位とが等しいものとして、帯状体１の幅方向の張力分布を直線ばね１３のばね定数の変化として把握するとともに、各計測点１ａにおける幅方向の曲げ剛性を、各節点１１における回転ばね１５のばね定数の変化として把握するようにモデル化したものである。

図３は、上述した張力分布測定装置を用いて張力分布を測定する第１の実施形態の手順を示す。まず、帯状体１に振動荷重負荷装置３で振動荷重を負荷し（ステップ１）、変位計４によって帯状体１の振動変位を計測して（ステップ２）、計測された振動変位から、振動特性算出部５ａで帯状体１の固有振動数と振動モードを算出し（ステップ３）、モデル化部５ｂで帯状体１をモデル化する（ステップ４）。こののち、モデルの直線ばね１３と回転ばね１５の各ばね定数の初期値を設定し（ステップ５）、固有値解析部５ｃでモデルの固有振動数と振動モードを算出して（ステップ６）、後述する（１７）式の評価関数を用いて、計測から求められた帯状体１の固有振動数および振動モードとの各差の二乗和を計算する（ステップ７）。つぎに、この二乗和の評価値が最小となるように収束したか否かを判定し（ステップ８）、収束していない場合は、各ばね定数の値を修正してステップ６に戻り、モデルの固有振動数と振動モードを繰り返し計算する。この繰り返し計算のループはばね定数算出部５ｄで行われる。収束した場合は、ばね定数変換部５ｅで、同定された直線ばね１３と回転ばね１５の各ばね定数をそれぞれ張力値と曲げ剛性値に換算し（ステップ９）、張力分布算出部５ｆで張力分布を算出して（ステップ１０）、測定を終了する。

以下に、上述した第１の実施形態で帯状体１の張力分布を測定する方法を、具体的に説明する。ここでは、測定点１ａおよびモデルの節点１１の数を一般化してｎとする。

前記変位計４で計測された帯状体１の振動変位ｖは、（１）式のようにモードベクトルｖ_ｉで表される。
ｖ_ｉ＝｛ｖ_i,1 ｖ_i,2 ・・・ｖ_i,n｝^Ｔ（ｉ＝１〜ｍ）（１）
ここに、ｉは幅方向の振動モードのモード次数である。

一方、前記２次元多質点系モデルの運動方程式は、帯状体１の質量マトリクスをＭ、各節点１１の変位ベクトルをｚ、直線ばね１３のばね定数ｋに相当する張力剛性マトリクスをＫ_Ｔ、回転ばね１５のばね定数τに相当する曲げ剛性マトリクスをＫ_τとすると、（２）式で表される。
ｚ＝｛ｚ₁ ｚ₂ ・・・ｚ_n｝^Ｔ（３）

質量マトリクスＭは、帯状体１の寸法および物性から算出される既知行列であり、（４）式で表される。
質量マトリクスＭの各要素に含まれる質量ｍと慣性モーメントＪは、長手方向の振動モードを正弦波で近似して縮小した場合の等価質量と等価剛性から算出される。支持ロール２ａ、２ｂ間の距離をＬ、帯状体１の板厚をｔ、密度をρ、リンク１４の長さをｌとすると、支持ロール２ａ、２ｂ間の中心における等価質量ｍは、長手方向の振動モードを１次とすると、（５）式で求められる。
ｍ＝ρｔｌＬ／２（５）
また、慣性モーメントＪは、（５）式で求めた等価質量ｍから（６）式で求められる。
Ｊ＝ｍ（ｔ^２＋ｌ^２）／１２（６）

張力剛性マトリクスＫ_Ｔと曲げ剛性マトリクスＫ_τは未知行列であり、張力剛性マトリクスＫ_Ｔは直線ばね１３のばね定数ｋ_ｊ（ｊ＝１〜ｎ）を用いて（７）式で表される。
また、曲げ剛性マトリクスＫ_τは回転ばね１５のばね定数τ_ｊ（ｊ＝１〜ｎ−２）を用いて（８）式で表される。
曲げ剛性マトリクスＫ_τの各要素に含まれる幅方向の曲げ剛性に相当するばね定数τ_ｉが算出されたときの断面二次モーメントＩ_ｊへの換算式は（９）式で表される。
Ｉ_ｊ＝τ_ｊｌ／Ｅ（９）

つぎに、（１）式の帯状体１の振動変位ｖ_ｉから得られる固有振動数ωおよび幅方向の振動モードから、幅方向の張力分布と曲げ剛性分布を求める方法について述べる。剛性行列Ｋ_Ｔ、Ｋ_τと質量行列Ｍは、一般的に固有振動数と振動モードとの間に（１０）式の関係がある。
ＭΦΩ^２＝（Ｋ_Ｔ＋Ｋ_τ）Φ （１０）
Ω^２は幅方向振動モードの固有振動数ω_ｉ ^２（ｉ＝１〜ｎ）を対角要素とする行列で、（１１）式で表される。
Λはｉ次の固有値λ_ｉ（ｉ＝１〜ｎ）を対角要素とする行列である。また、Φは対応するモードベクトルを列ベクトルφ_ｉ（ｉ＝１〜ｎ）とするモード行列であり、（１２）式で表される。
Φ＝〔φ_１ φ_２・・・ φ_ｎ〕（１２）
φ_ｉ＝｛φ_i,1 φ_i,2 ・・・ φ_i,n｝^Ｔ（１３）
φ_ｉは固有値解析により得られるｉ次の振動モードベクトルで、（１）式で示した計測値のモードベクトルｖ_ｉに相当するが、計測値と解析値を区別するために、便宜上別の記号を用いる。

つぎに、未知の直線ばね１３のばね定数ｋ_ｊ（ｉ＝１〜ｎ）と回転ばね１５のばね定数τ_ｊ（ｉ＝１〜ｎ）を繰り返し計算によって求める方法を説明する。各ばね定数ｋ_ｊ、τ_ｊは未知数であるため、仮の初期値を設定して（１０）〜（１３）式から固有値λ_ｉと振動モードベクトルφ_ｉを計算し、繰り返し計算により計測値と解析値の誤差が最小となる各ばね定数ｋ_ｊ、τ_ｊを算出する。具体的には、（１４）式に示す評価関数を用いて、評価値Ｓが最小となる各ばね定数ｋ_ｊ、τ_ｊを求める。
この評価関数は、計測されたモードベクトルｖ_ｉと振動モードベクトルφ_ｉの各成分の差と、計測から得られた固有振動数ω_ｉの二乗と固有値λ_ｉの差を固有振動数ω_ｉの二乗で除算した値とを計算し、これらを二乗和するものであり、振動モードベクトルφ_ｉと固有値λ_ｉとが直線ばね１３と回転ばね１５の各ばね定数ｋ_ｊ、τ_ｊから計算されることから、パラメータとして各ばね定数ｋ_ｊ、τ_ｊを有している。多質点系モデルの節点１１での各ばね定数ｋ_ｊ、τ_ｊは物理的に正の値をとるので、ｋ_ｊ＞０、τ_ｊ＞０であることを拘束条件として、評価関数の評価値Ｓが最小となるように最適化問題を解く。最適化の方法としては、最急降下法、準ニュートン法等を用いることができ、これにより、各ばね定数ｋ_ｊ、τ_ｊが算出される。

つぎに、算出した直線ばね１３のばね定数ｋ_ｊを各節点１１すなわち各測定点１ａでの張力Ｔ_ｊに換算する。距離Ｌの支持ロール２ａ、２ｂ間の帯状体１を幅方向にリンク１４の長さｌで分割した１要素を弦と見なして、各測定点１ａでの張力Ｔ_ｊが（１５）式で与えられる。
また、帯状体１の総張力Ｔ_totalは、各測定点１ａでの張力Ｔ_ｊの合計として（１６）式で与えられる。

一方、算出した回転ばね１５のばね定数τ_ｊは、（９）式によって断面二次モーメントＩ_ｊに換算され、幅方向の曲げ剛性分布が算出される。

図４は、上述した張力分布測定装置を用いて張力分布を測定する第２の実施形態の手順を示す。帯状体１をモデル化するステップ４までは、第１の実施形態のものと同じである。第２の実施形態では、こののち、固有値解析部５ｃで、モデルの直線ばね１３と回転ばね１５の各ばね定数に対する係数行列の擬似逆行列を計算し（ステップ５）、特異値分解によって擬似逆行列の誤差成分を分離する（ステップ６）。つぎに、ばね定数算出部５ｄで、誤差成分を除去した擬似逆行列を用いて、最小二乗法によって直線ばね１３と回転ばね１５の各ばね定数を算出する（ステップ７）。各ばね定数を算出した後は第１の実施形態のものと同じであり、ばね定数変換部５ｅで、同定された直線ばね１３と回転ばね１５の各ばね定数をそれぞれ張力値と曲げ剛性値に換算して（ステップ８）、張力分布算出部５ｆで張力分布を算出し（ステップ９）、測定を終了する。

以下に、上述した第２の実施形態で帯状体１の張力分布を測定する方法を、具体的に説明する。ステップ１〜４およびステップ８〜９は第１の実施形態のものと同じであるので、ここではステップ５〜７についてのみ説明する。また、測定点１ａおよびモデルの節点１１の数は一般化してｎとする。

この実施形態では、（１０）式の関係を利用して、計測および固有値解析の各々から得られる固有振動数および振動モードが一致する直線ばね１３と回転ばね１５の各ばね定数を計算する。使用する固有振動数と振動モードの次数は、各ばね定数の同定計算を行う上ではいくつでもよく、ここでは、（１２）式におけるモード行列Φに、（１７）式に示すように、１〜ｍ次モードの列ベクトルφ_ｉ（ｉ＝１〜ｍ、ｍ≦ｎ）を使用する場合について述べる。
Φ＝〔φ_１ φ_２・・・ φ_ｍ〕（１７）
（１２）式に（１７）式を代入して書き直すと、モード毎に独立した自由度ｎのベクトル方程式が（１８）式に示すように得られる。
ω_i ^２Ｍφ_ｉ＝Ｋ_Ｔφ_ｉ＋Ｋ_τφ_ｉ（１８）
（１８）式の右辺第１項は、（１９）式に示すように、既知の係数行列Φ_T,iと、未知の直線ばね１３のばね定数ｋ_ｉのベクトルｋとの積にまとめることができる。
（１８）式の右辺第２項も同様に、（２０）式に示すように、既知の係数行列Φ_τ,iと、未知の回転ばね１５のばね定数τ_ｉのベクトルτとの積にまとめることができる。

（１９）式と（２０）式を（１８）式に代入してモード毎のベクトル方程式を、未知ベクトルｋとτについてまとめると、（２１）式に示すように、ｎ×ｍ個の独立した方程式が得られる。
（２１）式の左辺の既知ベクトルをｋ_Φ、右辺のベクトルｋの係数行列をΦ_Ｔ、ベクトルτの係数行列をΦ_τとおいて、（２２）式で表す。

つぎに、振動モードと固有値の関係、および振動モードの直交性を利用した定式化を行う。（１０）式に左側から転置行列Φ^Ｔを乗じると（２３）式が得られる。
Φ^Ｔ（Ｋ_Ｔ＋Ｋ_τ）Φ＝Φ^ＴＭΦΩ^２（２３）
（２３）式の関係を利用して、計測および固有値解析の各々から得られる固有振動数および振動モードが一致する直線ばね１３と回転ばね１５の各ばね定数ｋ_ｉ、τ_ｉを計算する。（２３）式の両辺は対称行列であるから、上三角成分または下三角成分についてのｍ（ｍ＋１）／２個の独立した等式が（２４）式のように得られる。
また、マトリクスＫ_Ｔは（７）式に示したように対角行列であるから、（２４）式の左辺の第１項について展開し、未知のばね定数ｋ_ｉ（ｉ＝１〜ｎ）についてまとめると、（２５）式のように既知の係数行列Ψ_Ｔと未知のベクトルｋとの積で表される。
同様に、（２４）式の左辺の第２項は、（２６）式のように既知の係数行列Ψ_τと未知のベクトルτとの積で表される。

（２５）式と（２６）式を（２４）式に代入し、右辺の既知ベクトルをｋ_Ψとすると（２７）式が得られる。
最終的に（２２）式と（２７）式を（２８）式のようにまとめることにより、未知数２ｎ−２個に対してｍ（ｍ＋１）／２＋ｎｍ個の等式が得られる。

（２８）式を最小二乗法で解けば直線ばね１３と回転ばね１５の各ばね定数ｋ_ｉ、τ_ｉを求めることができる。ここで、（２８）式の左辺における係数行列をＡ、未知ベクトルをｘ、右辺の既知ベクトルをｂとおいて、（２９）式のように表す。
Ａｘ＝ｂ（２９）
係数行列Ａの擬似逆行列Ａ^†を用いて（２９）式を書き直すと、（３０）式と（３１）式が得られる。
ｘ＝Ａ^†ｂ（３０）
Ａ^†＝（Ａ^ＴＡ）^−１Ａ^Ｔ（３１）
（３１）式で計算されるＡ^†を用いて（３０）式の未知ベクトルｘ、すなわち、各ばね定数ｋ_ｉ、τ_ｉを算出することができる。

ところで、計算に使用するモード数が少ない場合や、帯状体１の幅方向分割数が大きい場合は、行列Ａのランクが不足し、（３０）式の解から求められる張力が実際の張力分布と一致しなくなることがある。そこで、特異値分解を利用して望ましい解を求める。係数行列Ａは特異ベクトルを列ベクトルとする直行行列Ｐ、Ｑと特異値を対角成分に持つ対角行列Ｄを用いて（３２）式のように分解できる。
Ａ＝ＰＤＱ^Ｔ（３２）
対角行列Ｄを特異値の大きい部分行列Ｄ_ｓと特異値の小さい部分行列Ｄ_ｎに分けると、（３１）式の擬似逆行列Ａ^†は（３３）式のように表される。
Ａ^†＝Ｑ_ｓＤ_ｓ ^−１Ｐ_ｓ ^Ｔ（３３）
この（３３）式で計算されるＡ^†を（３０）式に用いることにより、より良好な結果を得るように、未知ベクトルｘ、すなわち、各ばね定数ｋ_ｉ、τ_ｉを算出することができる。

図５は、片伸びの不均一歪みがある帯状体１に対して、２次モードまで用いた（２８）式の係数行列Φ_Ｔを特異値分解した結果の例を示す。この特異値分解結果より、Ｄ_ｓ＝diag(σ₁，・・・σ₁₈)、Ｄ_ｎ＝diag(σ₁₉，σ₂₀，・・・)＝０とした（３３）式の擬似逆行列Ａ^†を用いて算出した張力分布を図６に示す。図６には、特異値分解をせずに（３１）式の擬似逆行列Ａ^†を用いて算出した張力分布と、３次元ＦＥＭ解析によって求めた張力分布も併せて示す。特異値分解なしの張力分布が、片伸びの存在する幅端部で大きく変動しているのに対して、特異値分解ありの張力分布は、正解値とするＦＥＭ解析結果とよく一致している。このことにより、特異値分解によって誤差成分を除去することにより、一般的な擬似逆行列を用いる最小二乗法よりも精度よく張力分布を測定できることが分かる。

実施例として、上述した第２の実施形態の張力測定装置を用いて、帯状体１としての板厚が１．０ｍｍの薄鋼板と板厚が０．５ｍｍのアルミニウム薄板とについて、不均一歪みが存在しないフラットな歪みパターンと、それぞれ中伸びと耳波の不均一歪みが一部顕在する歪みパターンの張力分布を測定した。いずれの場合も板幅は１０００ｍｍ、支持ロール２ａ、２ｂ間の距離Ｌは４５００ｍｍとし、長手方向に付与した総張力は、薄鋼板の場合は２９４００Ｎ、アルミニウム薄板の場合は４９００Ｎとした。また、比較例として、不均一歪みが全て潜在化すると仮定し、幅方向の曲げ剛性を既知とした特願２０１０−２７４３４５で提案した張力測定装置を用いた張力分布の測定も行った。これらの実施例と比較例の張力分布の測定結果を、３次元ＦＥＭ解析で張力分布を求めた解析結果と対比した。

図７（ａ）、（ｂ）、（ｃ）は、薄鋼板に対する張力分布の測定結果を示す。いずれの歪みパターンでも、不均一歪みの一部顕在による幅方向の曲げ剛性の変化を考慮した実施例の各張力分布が正解値とするＦＥＭ解析結果とよく一致しているのに対して、不均一歪みが全て潜在化すると仮定した比較例の各張力分布は、ＦＥＭ解析結果と大きくずれるように張力が変動している。この結果より、本発明に係る張力測定装置は、不均一歪みが一部顕在しても帯状体１の張力分布を精度よく測定できることが分かる。

図８（ａ）、（ｂ）、（ｃ）は、アルミニウム薄板に対する張力分布の測定結果を示す。この場合も、実施例の各張力分布は、比較例の各張力分布よりも正解値のＦＥＭ解析結果とよく一致している。ただし、比較例の各張力分布のＦＥＭ解析結果とのずれは、薄鋼板の場合よりも小さくなっている。これは、アルミニウムのヤング率が鋼よりも小さいことと、アルミニウム薄板の板厚が薄鋼板よりも薄いことのために、不均一歪みの一部顕在による幅方向曲げ剛性の変化が少ないためと思われる。

図９（ａ）、（ｂ）、（ｃ）および図１０（ａ）、（ｂ）、（ｃ）は、それぞれ薄鋼板とアルミニウム薄板の場合について、板幅方向の曲げ剛性分布に対応する板幅方向の断面二次モーメント分布を、本発明に係る張力測定装置で算出した結果を示す。いずれの場合も、中延びの歪みパターンでは、不均一歪みが顕在する幅中央部で断面二次モーメント、すなわち曲げ剛性が増大し、耳波の歪みパターンでは、不均一歪みが顕在する幅両端部で曲げ剛性が増大している。このことより、幅方向曲げ剛性を既知の一定値とする従来の張力測定装置では、不均一歪みが顕在する帯状体の張力分布の測定精度に限界があることが分かる。なお、アルミニウム薄板の場合は、各歪みパターンでの断面二次モーメントの増大量が薄鋼板の場合よりも１オーダー小さくなっており、図８において、比較例の張力分布とＦＥＭ解析結果とのずれが比較的小さかったことを裏付けている。

１帯状体
１ａ測定点
２ａ、２ｂ支持ロール
３振動荷重負荷装置
４変位計
５演算装置
５ａ振動特性算出部
５ｂモデル化部
５ｃ固有値解析部
５ｄばね定数算出部
５ｅばね定数換算部
５ｆ張力分布算出部
１１節点
１２固定面
１３直線ばね
１４リンク
１５回転ばね

Claims

長手方向に張力を付与された帯状体の幅方向の張力分布を、長手方向の２箇所の部位で支持された支持部位間で測定する張力測定装置において、前記２箇所の支持部位間で前記帯状体に振動荷重を負荷する振動荷重負荷手段と、この振動荷重負荷手段によって生じる振動変位を帯状体の幅方向の複数の測定点で計測する変位計測手段とを備え、前記帯状体について、前記幅方向の複数の各測定点に対応する各節点に作用する張力を模擬する直線ばねと、該各節点における幅方向の曲げ剛性を模擬する回転ばねとを含む２次元の多質点系モデルにモデル化するモデル化手段と、前記多質点系モデルの固有値解析から得られる前記各節点での固有振動数および振動モードが、前記計測された各測定点での振動変位から得られる固有振動数および振動モードと一致するような前記直線ばねおよび回転ばねの各ばね定数を算出するばね定数算出手段と、算出された前記直線ばねおよび回転ばねの各ばね定数を、それぞれ前記各測定点における張力値および曲げ剛性値に換算するばね定数換算手段とを有し、換算された前記各測定点における張力値から前記帯状体の幅方向の張力分布を求めることを特徴とする張力測定装置。
前記ばね定数算出手段が、前記固有値解析から得られる前記各節点での固有振動数および振動モードと、前記各測定点での振動変位から得られる固有振動数および振動モードとの各差の二乗和が最小となるような前記直線ばねおよび回転ばねの各ばね定数を、繰り返し計算を用いた最小二乗法によって算出するものである請求項１に記載の張力測定装置。
前記ばね定数算出手段が、前記帯状体の振動モードに対する前記直線ばねおよび回転ばねの各ばね定数の係数行列の擬似逆行列を計算し、前記固有値解析から得られる前記各節点での固有振動数および振動モードと、前記各測定点での振動変位から得られる固有振動数および振動モードとが一致するような前記直線ばねおよび回転ばねの各ばね定数を、前記擬似逆行列を用いた最小二乗法によって算出するものである請求項１に記載の張力測定装置。
前記計算された擬似逆行列から特異値分解により誤差成分を除去し、この誤差成分を除去した擬似逆行列を、前記最小二乗法による各ばね定数の算出に用いるようにした請求項３に記載の張力測定装置。