(57)【要約】  広い観測条件によって収集されたデータに基づいて荷重の較正を行い、複数の風車ブレードの荷重の較正を迅速に行うことを目的とする。風車ブレード１０のピッチ角が可変とされている風車に適用される荷重計測装置１であって、風車ブレード１０の歪を求めるためのセンサと、風車ブレード１０の歪と風車ブレード１０の荷重との関係を表した関数を有し、関数にセンサの計測データに基づく歪を用いることで、風車ブレード１０の荷重を求める荷重算出部２０と、空力トルクの最大と最小とのばらつきが所定値以下となる風車ブレード１０のピッチ角範囲および回転数範囲で取得されたセンサの計測データに基づいて、関数を較正する較正部３０とを具備する荷重計測装置。
2010122658.tif 000002
【特許請求の範囲】
【請求項１】
風車ブレードのピッチ角が可変とされている風車に適用される荷重計測装置であって、
該風車ブレードの歪を求めるためのセンサと、
該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求める荷重算出手段と、
風速による空力トルクの最大と最小とのばらつきが所定値以下となる前記風車ブレードのピッチ角範囲および回転数範囲で取得された前記センサの計測データに基づいて、該関数を較正する較正手段と
を具備する荷重計測装置。
【請求項２】
前記較正手段は、前記空力トルクが所定値以下となる前記風車ブレードのピッチ角範囲および回転数範囲で取得された前記センサの計測データに基づいて、前記関数を較正する請求項１に記載の荷重計測装置。
【請求項３】
前記較正手段は、
無風時における前記風車ブレードの荷重及び前記風車ブレードのピッチ角並びにアジマス角が関連付けられたテーブルと、
前記センサによって計測データが取得されたときの前記風車ブレードのピッチ角及びアジマス角に対応する風車ブレードの荷重を前記テーブルから取得する荷重取得手段と、
前記センサの計測データから前記風車ブレードの歪を算出する歪算出手段と、
前記荷重取得手段によって取得された前記風車ブレードの荷重と前記歪算出手段によって算出された歪との関係に基づいて前記関数のパラメータを較正するパラメータ算出手段と
を具備する請求項１または請求項２に記載の荷重計測装置。
【請求項４】
前記較正手段は、
前記荷重取得手段によって取得された風車ブレードの荷重と、前記センサの計測データとに基づいて、無荷重時における前記センサの計測データを求め、この無荷重時における計測データを用いて前記センサの計測データをオフセット較正する請求項３に記載の荷重計測装置。
【請求項５】
前記センサは、
前記風車ブレードを挟んで対向する位置に設けられた一対の第１センサと、
該第１センサと異なる位置であって、該風車ブレードを挟んで対向する位置に設けられた一対の第２センサと
を具備する請求項１から請求項４のいずれかに記載の荷重計測装置。
【請求項６】
前記センサは、
前記風車ブレードを挟んで対向する位置であって、前記第１センサ及び前記第２センサと異なる位置であって、かつ、前記第１センサまたは前記第２センサのどちらかと平行する位置に設けられた一対の第３センサと
を具備する請求項１から請求項５のいずれかに記載の荷重計測装置。
【請求項７】
風車ブレードのピッチ角が可変とされている風車に適用される荷重計測装置であって、
該風車ブレードの歪を求めるためのセンサと、
該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求める荷重算出手段と、
風速が３メートル以下である場合に、第１のアジマス角及び該第１のアジマス角から１８０度回転させた第２のアジマス角からなる２地点において、ピッチ角を最小ピッチ角と最大ピッチ角とに設定したときの各前記センサの計測データに基づいて、該関数を較正する較正手段と
を具備する荷重計測装置。
【請求項８】
風車ブレードのピッチ角が可変とされている風車に適用される荷重計測方法であって、
該風車ブレードの歪を求め、
該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求め、
風速による空力トルクの最大と最小とのばらつきが所定値以下となる前記風車ブレードのピッチ角範囲および回転数範囲で取得された前記センサの計測データに基づいて、該関数を較正する
荷重計測方法。
【請求項９】
風車ブレードのピッチ角が可変とされている風車に適用される荷重計測プログラムであって、
該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求める第１処理と、
風速による空力トルクの最大と最小とのばらつきが所定値以下となる前記風車ブレードのピッチ角範囲および回転数範囲で取得された前記センサの計測データに基づいて、該関数を較正する第２処理と
をコンピュータに実行させるための荷重計測プログラム。
【請求項１０】
風車ブレードのピッチ角が可変とされている風車に適用される荷重計測方法であって、
該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求め、
風速が３メートル以下である場合に、第１のアジマス角及び該第１のアジマス角から１８０度回転させた第２のアジマス角からなる２地点において、ピッチ角を最小ピッチ角と最大ピッチ角とに設定したときの各前記センサの計測データに基づいて、該関数を較正する
荷重計測方法。
【請求項１１】
風車ブレードのピッチ角が可変とされている風車に適用される荷重計測プログラムであって、
該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求める第１処理と、
風速が３メートル以下である場合に、第１のアジマス角及び該第１のアジマス角から１８０度回転させた第２のアジマス角からなる２地点において、ピッチ角を最小ピッチ角と最大ピッチ角とに設定したときの各前記センサの計測データに基づいて、該関数を較正する第２処理と
をコンピュータに実行させる荷重計測プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、荷重測定装置及び方法並びにプログラムに関するものである。
【背景技術】
【０００２】
一般に、風力発電装置には、風車ブレードにかかる荷重を計測するセンサが翼根部等に取り付けられており、これらセンサによって計測されたデータを処理することで荷重を算出している。しかしながら、翼の製造時やセンサの取り付け時に生じる個体差に起因して、各風車ブレードに加わる荷重と歪との関係は一定でないこと等から、風車ブレード毎に荷重を測定し、荷重値の較正を行うことが提案されている。
【特許文献１】米国特許第６，９４０，１８６号公報明細書
【発明の開示】
【０００３】
従来、上記荷重の較正は人手によって（ロックピン等を用いて）風車ロータが回転しないように固定した状態で各風車ブレードの荷重を測定し、行っていた。
【０００４】
しかしながら、上記較正作業は風車毎に行う必要があり、全ての風車に対して較正を行うためには、人手によってロータを固定しなければならず、数百台もの風力発電装置が設置された大規模なウィンドファームの場合には、膨大な作業時間が必要であった。また、風車ブレードに加わる荷重と歪との関係は風車毎、風車ブレード毎に異なるため、それぞれの風車ブレードに対し、繰り返し較正作業を行う必要があった。また、ターニングモータを用いて所定の位置（角度）まで翼を移動させ、固定するまでに時間がかかってしまい、作業をスムーズに行えず、作業効率が悪いという問題があった。
【０００５】
本発明は、上記問題を解決するためになされたもので、観測条件によらずに、風車ブレードの荷重の較正を効率的に行うことができる荷重測定装置及び方法並びにプログラムを提供することを目的とする。
【０００６】
本発明の第１の態様は、風車ブレードのピッチ角が可変とされている風車に適用される荷重計測装置であって、該風車ブレードの歪を求めるためのセンサと、該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求める荷重算出手段と、風速による空力トルクの最大と最小とのばらつきが所定値以下となる前記風車ブレードのピッチ角範囲および回転数範囲で取得された前記センサの計測データに基づいて、該関数を較正する較正手段とを具備する荷重計測装置である。
【０００７】
このような構成によれば、荷重算出手段が保有している関数を較正する較正手段が設けられており、この較正手段が、風速による空力トルクの最大と最小とのばらつきが所定値以下となる風車ブレードのピッチ角範囲および回転数範囲で取得されたセンサの計測データに基づいて関数を較正するので、風速の条件を広範にすることができる。
【０００８】
上記荷重計測装置において、前記較正手段は、前記空力トルクが所定値以下となる前記風車ブレードのピッチ角範囲および回転数範囲で取得された前記センサの計測データに基づいて、前記関数を較正することとしてもよい。
【０００９】
これにより、空力トルクが所定値以下となる風車ブレードのピッチ角範囲で取得されたセンサの計測データを用いるので、空力トルクの影響を無視することができる。
【００１０】
上記荷重計測装置において、前記較正手段は、無風時における前記風車ブレードの荷重及び前記風車ブレードのピッチ角並びにアジマス角が関連付けられたテーブルと、前記センサによって計測データが取得されたときの前記風車ブレードのピッチ角及びアジマス角に対応する風車ブレードの荷重を前記テーブルから取得する荷重取得手段と、前記センサの計測データから前記風車ブレードの歪を算出する歪算出手段と、前記荷重取得手段によって取得された前記風車ブレードの荷重と前記歪算出手段によって算出された歪との関係に基づいて前記関数のパラメータを較正するパラメータ算出手段とを具備することとしてもよい。
【００１１】
このように、較正手段においてテーブルには無風時における風車ブレードの荷重及び風車ブレードのピッチ角並びにアジマス角が関連付けられ、計測データがセンサによって取得されたときの風車ブレードのピッチ角及びアジマス角に対応する風車ブレードの荷重が、荷重取得手段によってテーブルから取得され、歪算出手段によってセンサの計測データから風車ブレードの歪が算出され、荷重取得手段によって取得された風車ブレードの荷重と歪算出手段によって算出された歪との関係に基づいてパラメータ算出手段によって関数のパラメータが較正される。
【００１２】
これにより、較正手段が有するテーブルは、アジマス角とピッチ角とに荷重が対応付けられているので、計測データが取得されたときのアジマス角とピッチ角がわかれば、その時の風車ブレードの荷重を簡便に把握することができる。また、計測データに基づいて算出された風車ブレードの歪と、計測データに基づいて決定される荷重との関係に基づいて関数のパラメータが較正されるので、計測データの歪を精度高く較正することができる。
【００１３】
上記荷重計測装置において、前記較正手段は、前記荷重取得手段によって取得された風車ブレードの荷重と、前記センサの計測データとに基づいて、無荷重時における前記センサの計測データを求め、この無荷重時における計測データを用いて前記センサの計測データをオフセット較正することとしてもよい。
【００１４】
これにより、センサの計測データに含まれる無荷重時における計測データが求められ、オフセットするので、計測データの精度を向上させることができる。
【００１５】
上記荷重計測装置において、前記センサは、前記風車ブレードを挟んで対向する位置に設けられた一対の第１センサと、該第１センサと異なる位置であって、該風車ブレードを挟んで対向する位置に設けられた一対の第２センサとを具備することとしてもよい。
【００１６】
これにより、１つの風車ブレードにおける異なる向きの荷重を計測できる。例えば、第１センサを風車ブレードの腹と背とに設け、第２センサを風車ブレードのエッジ側に設ければ、これらのセンサによって風車ブレードがフェザー側となる向きにかかる荷重と、ファイン側となる向きにかかる荷重とを計測できる。
【００１７】
上記荷重計測装置において、前記センサは、前記風車ブレードを挟んで対向する位置であって、前記第１センサ及び前記第２センサと異なる位置であって、かつ、前記第１センサまたは前記第２センサのどちらかと平行する位置に設けられた一対の第３センサとを具備することとしてもよい。
これにより、第３のセンサによって、荷重以外の情報の計測に使用することができる。
【００１８】
本発明の第２の態様は、風車ブレードのピッチ角が可変とされている風車に適用される荷重計測装置であって、該風車ブレードの歪を求めるためのセンサと、該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求める荷重算出手段と、風速が３メートル以下である場合に、第１のアジマス角及び該第１のアジマス角から１８０度回転させた第２のアジマス角からなる２地点において、ピッチ角を最小ピッチ角と最大ピッチ角とに設定したときの各前記センサの計測データに基づいて、該関数を較正する較正手段とを具備する荷重計測装置である。
【００１９】
このような構成によれば、荷重算出手段が保有している関数を較正する較正手段が設けられており、この較正手段は、風速が３メートル以下である場合に、第１のアジマス角及び第１のアジマス角から１８０度回転させた第２のアジマス角からなる２地点において、ピッチ角を最小ピッチ角と最大ピッチ角とに設定したときの各センサの計測データに基づいて、関数を較正するので、少ない計測データに基づいて関数を較正できる。
【００２０】
本発明の第３の態様は、風車ブレードのピッチ角が可変とされている風車に適用される荷重計測方法であって、該風車ブレードの歪を求め、該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求め、風速による空力トルクの最大と最小とのばらつきが所定値以下となる前記風車ブレードのピッチ角範囲および回転数範囲で取得された前記センサの計測データに基づいて、該関数を較正する荷重計測方法である。
【００２１】
本発明の第４の態様は、風車ブレードのピッチ角が可変とされている風車に適用される荷重計測プログラムであって、該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求める第１処理と、風速による空力トルクの最大と最小とのばらつきが所定値以下となる前記風車ブレードのピッチ角範囲および回転数範囲で取得された前記センサの計測データに基づいて、該関数を較正する第２処理とをコンピュータに実行させるための荷重計測プログラムである。
【００２２】
本発明の第５の態様は、風車ブレードのピッチ角が可変とされている風車に適用される荷重計測方法であって、該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求め、風速が３メートル以下である場合に、第１のアジマス角及び該第１のアジマス角から１８０度回転させた第２のアジマス角からなる２地点において、ピッチ角を最小ピッチ角と最大ピッチ角とに設定したときの各前記センサの計測データに基づいて、該関数を較正する荷重計測方法である。
【００２３】
本発明の第６の態様は、風車ブレードのピッチ角が可変とされている風車に適用される荷重計測プログラムであって、該風車ブレードの歪と該風車ブレードの荷重との関係を表した関数を有し、該関数に該センサの計測データに基づく歪を用いることで、前記風車ブレードの荷重を求める第１処理と、風速が３メートル以下である場合に、第１のアジマス角及び該第１のアジマス角から１８０度回転させた第２のアジマス角からなる２地点において、ピッチ角を最小ピッチ角と最大ピッチ角とに設定したときの各前記センサの計測データに基づいて、該関数を較正する第２処理とをコンピュータに実行させる荷重計測プログラムである。
【００２４】
本発明によれば、観測条件によらずに、風車ブレードの荷重の較正を効率的に行うことができるという効果を奏する。
【図面の簡単な説明】
【００２５】
【図１】本発明の第１の実施形態に係る風力発電装置の全体の概略構成を示す図である。
【図２】翼根部を説明するための図である。
【図３】風車ブレードの翼根部から１．８メートルの位置における断面図の一例である。
【図４】翼根部から見たセンサ位置の配置を説明するための図である。
【図５】本発明の第１の実施形態に係る荷重測定装置の概略構成を示すブロック図である。
【図６】風速毎の空力トルク、発電機回転数、及びピッチ角の関係を示す図である。
【図７】較正部が有するテーブルの一例を示した図である。
【図８】アジマス角を説明するための図である。
【図９】計測データに基づいて歪と荷重との関係を示した一例の図である。
【図１０】翼根部における荷重と歪波長との関係を示した一例の図である。
【図１１Ａ】風速８メートル毎秒のときの領域Ａの範囲におけるピッチ角の計測データを含んで較正を行った場合の荷重（フラップ方向）と無風時の荷重とを比較の一例を示した図である。
【図１１Ｂ】風速８メートル毎秒のときの領域Ａの範囲におけるピッチ角の計測データを含んで較正を行った場合の荷重（エッジ方向）と無風時の荷重との比較の一例を示した図である。
【図１２Ａ】風速８メートル毎秒のときの領域Ａの範囲におけるピッチ角の計測データを含まずに較正を行った場合の荷重（フラップ方向）と無風時の荷重との比較の一例を示した図である。
【図１２Ｂ】風速８メートル毎秒のときの領域Ａの範囲におけるピッチ角の計測データを含まずに較正を行った場合の荷重（エッジ方向）と無風時の荷重との比較の一例を示した図である。
【符号の説明】
【００２６】
１  風力発電装置
７  センサ
２０  荷重算出部
３０  較正部
３１  テーブル
３２  荷重取得部
３３  歪算出部
３４  パラメータ算出部  
１００  荷重測定装置
【発明を実施するための最良の形態】
【００２７】
以下に、本発明に係る荷重測定装置及び方法並びにプログラムの一実施形態について、図面を参照して説明する。【００２８】
〔第１の実施形態〕
図１は、本実施形態に係る荷重測定装置１００を適用した風力発電装置の概略構成を示した図である。本実施形態に係る風力発電装置１は、風車ブレード１０のピッチ角が可変とされている風車である。
風力発電装置１は、図１に示されるように、支柱２と、支柱２の上端に設置されるナセル３と、略水平な軸線周りに回転可能にしてナセル３に設けられるロータヘッド（ハブ）４とを有している。ロータヘッド４には、その回転軸線周りに放射状に３枚の風車ブレード１０が取り付けられている。これにより、ロータヘッド４の回転軸線方向から風車ブレード１０に当たった風の力が、ロータヘッド４を回転軸線周りに回転させる動力に変換され、この動力が発電機によって電気エネルギーに変換されるようになっている。【００２９】
各風車ブレード１０には、風車ブレード１０の歪を求めるための複数のセンサ（センシング部）７が設けられている。センサ７は、例えば、ＦＢＧ（Fiber Bragg Grating sensor）である。ＦＢＧは、ブラッグ格子の刻まれた光ファイバセンサで、反射光の波長変化に基づいて、歪、熱膨張による格子間隔変化を検出する。
また、ロータヘッド４には、上記センサ７（センシング部）における計測結果を受信する信号処理部（図示略）を備えている。【００３０】
具体的には、各風車ブレード１０には、第１センサと第２センサと第３センサとが設けられている。第１センサ、第２センサ、及び第３センサは、それぞれ風車ブレード１０を挟んで対向する位置に設けられた一対のセンサを有する。好ましくは、第１センサ及び第２センサは、各々を構成する２つのセンサを結ぶ直線が略垂直に交わるように設けられているとよい。第３センサは、例えば、温度補償用に用いられるセンサであり、第１センサまたは第２センサの周辺に設けられている。【００３１】
図２は、風車ブレードに取り付けられたセンサ７（センシング部）の位置を説明するための図である。図２に示されるように、本実施形態において、センサ７は、例えば風車ブレード１０の付け根から１．８メートル離れた位置に設けられている。付け根とは、図２に示されるように風車ブレード１０とロータヘッド４との境界である。本実施形態においては、この付け根を「翼根部」という。【００３２】
図３は、センサ７が取り付けられた風車ブレード１０の翼根部から１．８メートルの位置における断面図を示した図である。図３において、風車ブレード１０の背側２１にはセンサＡ３、及び腹側２２にはセンサＡ１が設けられ、第１センサが構成されている。また、Ａ３と同じ位置にはセンサＡ５、Ａ１と同じ位置にはセンサＡ６が設けられ、第３センサが構成されている。また、風車ブレード１０の前縁２３の方向にはセンサＡ２、及び後縁２４の方向にはセンサＡ４が設けられ、第２センサを構成する。【００３３】
図４は風車ブレード１０の翼根部から見た場合の風車ブレード１０に取り付けられたセンサ７の配置を模式的に示した図である。図４に示すように、本実施形態では、センサＡ１が設けられている位置をＨＰ、センサＡ３が設けられている位置をＬＰ、センサＡ２が設けられている位置をＬＥ、センサＡ４が設けられている位置をＴＥと定義する。また、図４において、チルト角はタワーの鉛直軸線に対する風車ブレード１０の回転面の傾斜角度であり、運転時において風力によって風車ブレード１０が変形したとしても、風車ブレード１０とタワーとが接触するのを防ぐためにこのような傾斜角度が設けられている。このチルト角については、後述の計算においては無視してもよいし、考慮に入れても良い。【００３４】
次に、本実施形態に係る荷重測定装置１００の構成について詳細に説明する。
図５は、荷重測定装置１００が備える機能を展開して示した機能ブロック図である。
図５に示されるように、本実施形態に係る荷重測定装置１００は、荷重算出部（荷重算出手段）２０と、較正部（較正手段）３０とを備えている。【００３５】
荷重算出部２０は、風車ブレードの歪と風車ブレード１０の荷重との関係を表した関数を有し、この関数にセンサＡ１〜Ａ６の計測データに基づく歪を用いることで、風車ブレード１０の荷重を求める。【００３６】
較正部３０は、風速による空力トルクの最大と最小とのばらつきが所定値以下となる風車ブレード１０のピッチ角範囲および回転数範囲で取得されたセンサの計測データに基づいて、上述の関数を較正する。なお、空力トルクが所定値以下となる風車ブレード１０のピッチ角範囲の条件、および回転数範囲の条件を満たす期間に取得されたセンサの計測データを用いると、なおよい。【００３７】
具体的に、較正部３０によって関数の較正に用いられるセンサの計測データについて説明する。
図６は、風車ブレード１０をファイン側（ピッチ角２１度）から、フェザー側（ピッチ角１０９度）に変化させ、風車ブレード１０が停止されるまでの空力トルクの変化を風速毎（風速４ｍ毎秒から１２ｍ毎秒）に示した図である。ピッチ角が２１度や１０９度とは、風車ブレード１０を風車ロータ３に取り付けた時に決定される翼基準線の位置を０度とした場合の風車ブレード１０の角度である。ここで、ピッチ角０度とは、翼根断面上で定義される翼基準線上の角度であり、この線とロータ平面がなす角度とがピッチ角である。【００３８】
図６を得る際には、ピッチ角２１度から４５度の範囲では２．５度毎秒の速さ、ピッチ角４５度から１０９度の範囲では１．０度毎秒の速さで３つの風車ブレード１０を同時に変化させ、このときの３つの風車ブレード１０における計測データを取得する。また、風車ブレード１０は、ピッチ角を変化させることによって回転され、遊転の状態となっている。遊転とは、風力発電装置１が発電しない範囲で、風車ブレード１０を回転させる状態（例えば、低速で風車ブレード１０が回転する状態）である。
なお、３つの風車ブレード１０の全てに対して同様の処理を行うため、以下の説明では、１つの風車ブレード１０について説明することとする。【００３９】
図６に示されるように、風力発電装置が発電している期間と、停止等のためにピッチ角をフェザー側に開くことによりロータの回転を止めようとする力である空力ブレーキが強くかかる期間とを含めて領域Ａとする。図６に示されるように、領域Ａにおいては、風速に応じて異なる空力トルクがかかっている。その後、領域Ａより右側（換言すると、ピッチ角は６０度よりも大きく、かつ、発電機の回転数は０から３００ｒｐｍ（周波数６０Ｈｚの場合））である場合には、いずれの風速であってもトルクの値はおよそ−３００キロニュートンメートル以上となり、空力トルクの影響を無視できるような小さなトルクとなっている。【００４０】
なお、図６において、領域Ａより右側における回転数範囲の条件は、０から３００ｒｐｍとしていたが、これに限定されず、周波数に応じて設定するとよい。例えば、５０Ｈｚの場合には、回転数範囲の条件を０から２５０ｒｐｍとすることとしてもよい。【００４１】
このように、較正部３０は、領域Ａでない領域である部分、即ち、風速によらない空力トルクとなる範囲で取得された計測データを用いることとする。換言すると、空力トルクの最大と最小とのばらつきが所定値以下となる範囲のピッチ角である、ピッチ角６０度から１０９度の範囲で取得された計測データを用いることとする。なお、領域Ａでない領域のデータを用いる根拠については後述することとする。【００４２】
より具体的には、較正部３０は、テーブル３１、荷重取得部（荷重取得手段）３２、歪算出部（歪算出手段）３３、及びパラメータ算出部３４を備えている。【００４３】
較正部３０は、荷重取得部３２によって取得された風車ブレード１０の荷重と、センサ７の計測データとに基づいて無荷重時におけるセンサ７の計測データを求め、この無荷重時における計測データを用いてセンサ７の計測データをオフセット較正する。これにより、センサ７自身に含まれる測定誤差を考慮するので、較正の精度を向上することができる。【００４４】
テーブル３１は、無風時（較正を行う上で理想的な環境条件）における風車ブレード１０の翼根部における荷重及び風車ブレード１０のピッチ角並びにアジマス角が関連付けられたものであり、例えば、図７に示されるような一覧表（テーブル）として与えられ、風車ブレード１０のアジマス角とピッチ角との組み合わせに対応付けられた翼根部における荷重の値α１、α２、α３、α４・・・が格納されている。【００４５】
上記アジマス角とは、図８に示されるように、風車ブレード１０の回転面において、所定の基準と風車ブレード１０の軸線とのなす角をいい、本実施形態では、風車ブレード１０が最上部に位置したときを基準としている。従って、風車ブレード１０が風車の最上部に位置したときのアジマス角度は０度、最下部に位置したときのアジマス角度は１８０度である。【００４６】
図７に示されるような一覧表（テーブル）において、翼根部におけるモーメントは、以下に示す（１）式により、各センサＡ１〜Ａ４の位置における自重によるモーメントをそれぞれ計算し、これらのモーメントを座標変換することで、得ることができる。
Ｍ＝９．８×Ｗ×ｌｇ×ｓｉｎθ・ｃｏｓβ  ［Ｎｍ］    （１）【００４７】
上記（１）式において、Ｗは、風車ブレード１０の重量、ｌｇは、風車ブレード１０の翼根部から測定した重心位置（製造の段階で既知の値である）、θは、アジマス角とチルト角との関数、βは、ピッチ角とチルト角との関数である。【００４８】
荷重取得部３２は、センサ７によって計測データが取得されたときの風車ブレード１０のピッチ角及びアジマス角に対応する風車ブレード１０の翼根部における荷重をテーブル３１から取得する。【００４９】
歪算出部３３は、センサ７の計測データから風車ブレード１０の歪を算出する。具体的には、センサ７の計測データから歪波長を抽出し、歪波長を所定の関数に基づいて歪に変換する。より具体的には、センサ７の計測データは、ロータヘッド４に設けられる図示しない信号処理部において歪波長データが数値に変換され、数値として得られた歪波長を、歪εに変換する。また、歪εは、以下の（２）式で求められる。
ε＝Ｐｅ｛（λ−λｉ）−α（λＴ−λＴｉ）｝    （２）【００５０】
上記（２）式において、λは第１センサ（第２センサ）による計測データ、λｉは第１センサ（第２センサ）による無荷重時の計測データ、λＴは第３センサによる計測データ、λＴｉは第３センサによる無荷重時の算出データ、ｐｅは歪光学定数（８０９με／ｎｍ）、αは温度補償係数（２．２）である。なお、λｉは、計測データの平均値であり、以下の（３）式で求められる。
λｉ＝（λｍａｘ＋λｍｉｎ）／２      （３）【００５１】
上記（３）式において、λｍａｘは、データの最大ピーク値を示し、λｍｉｎは、最小ピーク値を示す。【００５２】
このように、第１センサを構成する各センサＡ１，Ａ３、第２センサを構成する各センサＡ２，Ａ４においてそれぞれ歪を算出するので、４つの歪が算出される。さらに、第１センサを構成する各センサＡ１，Ａ３から得られた歪の差分を算出することにより風車ブレード１０のフラップ方向（図４におけるＨＰ−ＬＰの方向）に対する歪εＦ、及び第２センサを構成する各センサＡ２，Ａ４における歪の差分を算出することにより風車ブレード１０のエッジ方向（図４におけるＬＥ−ＴＥの方向）に対する歪εＥを算出する。【００５３】
パラメータ算出部３４は、荷重取得部３２によって取得された風車ブレード１０の荷重と、歪算出部３３によって算出された歪との関係に基づいて関数のパラメータを較正する。具体的には、歪εＦ、εＥを算出した元データとなった計測データを、取得したタイミングにおけるアジマス角とピッチ角とに関連付けられた風車ブレード１０の翼根部における荷重と、歪εＦ、εＥとの関係に基づいて新たな関数を構成し、新たな関数の係数を用いて、荷重算出部２０に有する関数の係数を較正する。このとき、新たな関数は、フラップ方向について１つ、エッジ方向において１つ作成される。【００５４】
例えば、上記歪εをセンサの位置におけるモーメントに変換するには、以下の（４）式が用いられる。なお、ｄはセンサ７の設置位置（翼根部から１．８メートル）における風車ブレード１０の内径、Ｌは風車ブレード１０の番号（Ｌ＝１、２、３）、Ｅは翼材（ＦＲＰ）のヤング率、Ｉはセンサの設置位置における断面二次モーメント、Ｍｓｅｎｓｏｒはセンサ設置位置における曲げモーメント（荷重）、ε２Ｌ−１とε２Ｌとは一対のセンサ（第１センサまたは第２センサ）の計測データに基づく歪、ε２Ｌ−１，０とε２Ｌ，０とは第１センサまたは第２センサの歪の初期値である。
【数１】
2010122658.tif 000003
【００５５】
ここで、このモーメントを風車ブレード１０の翼根部におけるモーメントＭｒｏｏｔとセンサ７の設置位置（例えば、風車ブレード１０の翼根部から１．８ｍの場所）におけるモーメントＭｓｅｎｓｏｒとの比をβ（＞１）とすると、次式が得られる。
ただし、【００５６】
【数２】
2010122658.tif 000004
【００５７】
上記（５）式に示されるように、翼根部におけるモーメントＭｒｏｏｔは、係数ａ、ｂをパラメータとし、一対のセンサ（第１センサまたは第２センサ）による歪を変数とする一次関数として表わされる。
このことから、横軸に歪εＦまたはεＥとした場合に、縦軸にそのアジマス角及びピッチ角の場合のモーメントＭｒｏｏｔとするグラフを作成した場合に、この一次関数に基づいて得られる傾きａと切片ｂとを算出することにより、係数ａとｂとをパラメータとして算出することができる。【００５８】
次に、上記一次関数をグラフ化する方法について説明する。
横軸は各風車ブレード１０におけるフラップ方向の歪εＦとエッジ方向の歪εＥとし、縦軸はテーブル３１から取得した歪εＦと歪εＥとに対応した荷重（モーメント）Ｍｒｏｏｔとしてグラフ化し、このグラフから傾きａと切片ｂとを抽出する。より具体的には、図９に示されるようなグラフが作成され、このようなグラフに基づいてそれぞれの風車ブレード１０のフラップ方向の場合の係数ａ及びｂと、エッジ方向の場合の係数ａ及びｂを算出する。例えば、１番目の風車ブレード１０＃１におけるフラップ方向のパラメータとしてａ＝２．０１４×１０9、ｂ＝−０．４４８×１０３とする。同様にして、風車ブレード１０＃２、風車ブレード１０＃３に対してそれぞれフラップ方向とエッジ方向との係数ａ、およびｂを算出する。【００５９】
このように、パラメータ算出部３４は、係数ａ及びｂを算出すると、これを荷重算出部２０に出力する。これにより、荷重算出部２０の関数のパラメータは較正されることとなり、荷重算出部２０の関数にセンサから取得された計測データを用いることによって、得られる翼根部におけるモーメントが較正されることとなる。【００６０】
次に、本実施形態に係る荷重測定装置の較正部３０の作用について説明する。なお、各風車ブレード１０に対して行う処理は同じであるので、以下の説明においては、１つの風車ブレード１０に対する処理を例に挙げて説明する。
まず、本実施形態においては、ピッチ角を６０度から１０９度まで変化させ、このときの各センサＡ１〜Ａ６の計測データを取得する。計測データは、較正部３０の荷重取得部３２に与えられる。【００６１】
荷重取得部３２はテーブル３１を参照し、各センサＡ１〜Ａ６から計測された計測データのアジマス角とピッチ角との情報に関連付けられた翼根部における荷重を読み出す。続いて、荷重取得部３２は、センサ毎に、縦軸にセンサの計測データ、横軸に翼根部における荷重を示したグラフを作成し（図１０を参照）、これら各グラフから横軸の荷重の値が「０」であるときの歪波長をそれぞれ読み出す。この値は、無荷重時におけるセンサの計測データ、即ち、各センサのオフセット値に相当する。荷重取得部３２は、各センサのオフセット値とともに、テーブル３１から読み出した翼根部における荷重の情報及び該荷重が得られたときの各センサのデータを歪算出部３３に出力する。このように荷重取得部３２によって各センサのオフセット値が求められることにより、以降の処理において、センサ自身に含まれる測定誤差が較正され、荷重の測定精度を向上させることができる。【００６２】
歪算出部３３は、センサＡ１〜Ａ６の計測データから歪波長をそれぞれ抽出し、各センサから算出された計測データと、無荷重時の計測データとに基づいて、上記の（２）式を用いて、各センサ位置における歪（自重モーメント）εを算出する。例えば、センサＡ１における歪εＡ１は、以下の（２）´から求められる。
εＡ１＝Ｐｅ｛（λＨＰ−λＨＰｉ）−α（λＨＰＴ−λＨＰＴｉ）｝    （２）´【００６３】
上記（２）´において、λＨＰはセンサＡ１の歪波長データ、λＨＰＴは、センサＡ１の周辺に設置された温度補償用のセンサＡ５の歪波長データ、λＨＰｉはセンサＡ１のオフセット値（無荷重時の計測データ）、λＨＰＴｉはセンサＡ５のオフセット値（無荷重時の計測データ）である。【００６４】
歪算出部３３は、同様の算出処理により、センサＡ３，Ａ２，及びＡ４のそれぞれに対して歪εを算出する。これにより、歪εはεＡ１、εＡ２、εＡ３、εＡ４の合計４個算出される。歪算出部３３は、各センサＡ１〜Ａ４に対する歪εＡ１〜Ａ４を算出すると、これらの値と荷重取得部３２から入力された翼根部における荷重の情報とをパラメータ算出部３４に出力する。【００６５】
パラメータ算出部３４は、各センサに対する歪εＡ１〜Ａ４と、荷重取得部３２から入力された翼根部における荷重の情報との関係に基づいて、風車ブレード１０の歪と風車ブレード１０の翼根部における荷重との関係が示された関数のパラメータを較正する。
具体的には、風車ブレード１０の翼根部における荷重の情報と風車ブレード１０の歪との関係式は、上述した（５）式のように表される。具体的には、この関係式は、第１センサに対して一つ、第２センサに対して一つ求められ、一枚の風車ブレード１０に対して２つの関係式が作成される。例えば、第１センサの場合には、上記関係式は、以下の（５）´式で表わされる。
ＭＨＰ−ＬＰ＝ａ（εＡ１−εＡ３）×１０−６＋ｂ    ［Ｎｍ］    （５）´【００６６】
（５）式において、εＡ１は（２）´に基づいて算出されたセンサＡ１の位置における歪、εＡ３は同様の算出式に基づいて算出されたセンサＡ３の位置における歪である。【００６７】
ここで、上記（５）´式のフラップ方向（図４におけるＨＰ−ＬＰの方向）の歪の差分εＡ１−εＡ３を横軸とし、ＭＨＰ−ＬＰを縦軸とすると、図９に示されるようなグラフが得られる。このように、歪εＡ１〜Ａ４と翼根部における荷重の情報との関係をグラフ化すると、このグラフの切片と傾きを導くことができ、上記（５）´式の係数ａ（傾き）と係数ｂ（切片）とを算出できる。１つの風車ブレード１０において係数ａとｂの組み合わせがフラップ方向とエッジ方向との２つずつ算出される。【００６８】
そして、上記センサの計測から係数の算出までの処理を繰り返し行うことで、複数の係数を得、これらの係数の平均値を用いることにより、較正後の関数式を得る。例えば、最大値と最小値をそれぞれ省くことで、ノイズなどの突発したデータを取り除くことができる。なお、ここでノイズでないデータを取り除いても平均値の算出には影響を与えない。【００６９】
そして、このようにして、信頼度の高い関係式、換言すると、較正後の関係式が得られた後は、荷重計算部３２がこの関係式を用いて各センサの計測データから風車ブレード１０の翼根部における荷重を算出する。これにより、信頼度の高い荷重を算出することが可能となる。【００７０】
なお、本実施形態において、各センサＡ１〜Ａ６によって計測される計測データのうち、図６の領域Ａでない領域のデータを用いることとしていたが、これは較正の精度を向上させるためである。具体的には、上述した係数ａ及びｂを求める場合に、無風時に得られるモーメント（較正を行う場合の理想的なモーメント）に近似するためである。より具体的に、図１１Ａ、Ｂと、図１２Ａ、Ｂとを用いて説明する。【００７１】
図１１Ａは横軸に無風時のフラップ方向のモーメント、縦軸に風速８メートルの場合のフラップ方向のモーメントを示している。図１１Ｂは同様にエッジ方向のモーメントを比較している。図１１Ａ及び図１１Ｂは、図６の領域Ａを含む計測データによって係数ａ、ｂを求めた場合に得られる関数であり、理想的な無風時のデータと比較した場合である。理想的な無風時の関数の関数をｙ＝ｘであるのに対し、領域Ａを含む計測データを用いて係数ａ及びｂを算出すると、ｙ＝０．９７０１ｘ−３１．８８となる。また、無風時のモーメントと領域Ａを含む計測データから算出されるモーメントとを比較すると、ｙ＝ｘからずれた位置に点が多く存在している。これにより、領域Ａの計測データを含んで算出された係数ａ、ｂを用いて得られる荷重は、無風時の荷重と比較して誤差が大きくなることを示している。【００７２】
一方、図１２Ａ及び図１２Ｂは、図１１Ａ及び図１１Ｂと同様に横軸に無風時のモーメント、縦軸に風速８メートルの場合のモーメントを示しているが、図６の領域Ａにおいて得られた計測データを除外した計測データを用いて係数ａ及びｂを算出した場合の図である。例えば、図１２Ａに示されるように、領域Ａを含まない計測データを用いて係数ａ、ｂを算出した場合には、ｙ＝０．９９６８ｘ−３．３２２となる。また、これにより得られる荷重の結果は、グラフを見てもわかるように、無風時のモーメントと風速８メートルの場合のモーメントとがほぼ一致している。これにより領域Ａの計測データを含まないで得られた係数ａ、ｂによって関数が較正されることにより、無風時の荷重により近似させることができることとなる。【００７３】
なお、センサによる計測データを取得する場合は、風車ブレード１０を少なくともアジマス角０度から１８０度に１回、回転させればよいこととし、１８０度の回転で１つの較正用のデータファイルを生成することとするが、アジマス角の回転回数は特に限定されない。なお、本実施形態では、３６０度の回転で１つの較正用のデータファイルを生成するように構成されている。【００７４】
より具体的には、ピッチ角を１０９度から６０度、または６０度から１０９度に動かすことで、ロータは遊転をはじめる。ここでは、ピッチ角を１０９度から６０度に動かすことによって少なくともロータを１回転させ、１データファイルを取得する。同様に、ピッチ角を６０度から１０９度に動かすことによって少なくともロータを１回転させ１データファイルを取得する。【００７５】
即ち、本実施形態では、空力の影響を除く範囲において、ピッチ角をファイン側へ動作させる際に、少なくともロータが１回転する条件で一の較正用データファイルを取得する。同様に、ピッチ角をフェザー側へ動作させる際に一の較正用データファイルを取得する。【００７６】
本実施形態では、上述したようなピッチ角動作を行うことによって、１０の較正用データファイルを取得する。
さらに、このような１０回分の較正用のデータ（つまり、係数ａ及びｂ）が算出された場合には、これらの平均値を算出することにより、係数ａ及びｂの信頼性を検証するとよい。
また、Ｘは較正用のデータ（無荷重時の計測データ１８点、較正用データ（各風車ブレード１０のエッジ方向とフラップ方向のデータ）１２点）、Ｎはピッチ角を１０９度から６０度から１０９°度（１周期）と動作させた回数を示す。また、ｍを平均値とする。なお、ここでは、２Ｎ個の較正用データファイルから最大及び最小のデータを省き、２（Ｎ−１）個の平均値を求める。そして、キャリブレーション値ａ、ｂの平均値が下記の範囲条件を満たすか否かを検証する。基準値を満たす場合、較正用データの平均値を現地パラメータとする。【００７７】
【数３】
2010122658.tif 000005
〔平均値〕【００７８】
〔基準値の検証〕
1.７×１０９＜ａ＜２．７×１０９      （９）
−１００ｋＮｍ＜ｂ＜１００ｋＮｍ    （１０）【００７９】
以上説明してきたように、風車ブレード１０に設けられたセンサによって計測データを取得し、取得したデータに基づいて風車ブレード１０の歪と荷重とが算出される。荷重算出部２０に備えられた関数は、計測データに基づいて算出された風車ブレード１０の歪と荷重との関係から得られる新たな関数の係数によって較正される。これにより、計測データから簡便に新たな関数が算出できるので、簡便に較正のための係数を決定することができる。【００８０】
なお、このとき使用される計測データは、空力トルクの最大と最小とのばらつきが所定値以下となる風車ブレード１０のピッチ角範囲で取得されたセンサの計測データが用いられ、空力トルクの影響を無視できる計測データを使用するので、較正の精度を向上させることができる。
また、センサ自身の持つ歪を算出し、オフセットすることにより、より較正の精度を向上させることができる。さらに、計測データを取得するには、空力トルクの最大と最小のばらつきが所定値以下となる風車ブレード１０のピッチ角範囲（例えば、ピッチ角６０度から１０９度の範囲）で取得された計測データであればよく、特にアジマス角に限定されることがないので、広範な計測データを、較正のために用いることができる。【００８１】
また、このような計測データの取得、荷重及び歪の算出、及びその較正データの信頼性の検証を荷重測定装置１００によって行うので、較正作業にかかる時間を短縮できるとともに、ユーザにかかる負担を低減することができる。【００８２】
なお、上述した実施形態では、荷重測定装置としてハードウェアによる処理を前提としていたが、このような構成に限定される必要はない。例えば、各センサからの出力信号に基づいて別途ソフトウェアにて処理する構成も可能である。この場合、荷重測定装置は、ＣＰＵ、ＲＡＭ等の主記憶装置、上記処理の全て或いは一部を実現させるためのプログラムが記録されたコンピュータ読み取り可能な記録媒体を備えている。そして、ＣＰＵが上記記憶媒体に記録されているプログラムを読み出して、情報の加工・演算処理を実行することにより、上述の荷重測定装置と同様の処理を実現させる。
ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。【００８３】
〔変形例〕
また、本実施形態に係る荷重測定装置１００において、空力トルクの最大と最小とのばらつきが所定値以下となる風車ブレード１０のピッチ角範囲で取得されたセンサの計測データに基づいて、上述の関数を較正することとしていたが、これに限定されない。例えば、ピッチ角範囲に代えて風車ブレード１０の回転数の範囲とし、風速による空力トルクの最大と最小とのばらつきが所定値以下となる風車ブレードの回転数の範囲で取得されたセンサのデータに基づいて、関数を較正することとしてもよい。【００８４】
なお、本実施形態に係る風力発電装置１において、複数の風車ブレード１０の枚数は３枚としていたが、風車ブレード１０の枚数は特に限定されない。【００８５】
また、本実施形態に係る荷重測定装置１００において、１つの風車ブレード１０に対し、取り付けられるセンサは６個であったが、センサの個数は特に限定されない。【００８６】
また、本実施形態に係る荷重測定装置１００において、テーブル３１は、アジマス角とピッチ角とから算出されることとしていたが、これに限定されない。例えば、較正部３０に予めテーブルが与えられていることとしてもよい。【００８７】
〔第２の実施形態〕
次に、本発明の第２の実施形態について説明する。
本実施形態の荷重測定装置が第１の実施形態と異なる点は、アジマス角及びピッチ角の角度データを所定の値に限定し、風速を負の空力トルクが小さい範囲に限定してデータを取得する点である。以下、本実施形態の荷重測定装置について、第１の実施形態と共通する点については説明を省略し、異なる点について主に説明する。【００８８】
センサ７は、風速が３メートル以下である場合に、第１のアジマス角及び第１のアジマス角から１８０度回転させた第２のアジマス角からなる２地点において、ピッチ角を最小ピッチ角と最大ピッチ角とに設定したときの計測データを計測する。
より具体的には、センサは、風速が３メートル以下で、風車ブレード１０のアジマス角は９０度と２７０度の位置にある場合であり、かつ、それぞれのアジマス角においてピッチ角は２１度と１０９度に設定された状態である場合に計測データを取得する。【００８９】
１つの風車ブレード１０のデータを計測している場合には、他の２つの風車ブレード１０のピッチ角は、例えば、８５度等に設定し、遊転状態とする。【００９０】
このようにアジマス角を１８０度回転させた位置にある２地点から計測されたデータを用いることにより、パラメータ算出部３４から算出される関数に用いられる計測データは横軸に対して広範囲なデータを取得できることとなる。これにより、少ない計測データで係数ａとｂとを算出する場合の精度を向上させることができる。また、少ない計測データでパラメータが算出できるので、較正にかかる時間を短縮することができる。【００９１】
以上、本発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、本発明の要旨を逸脱しない範囲の設計変更等も含まれる。
【図１】
2010122658.tif 000006
【図２】
2010122658.tif 000007
【図３】
2010122658.tif 000008
【図４】
2010122658.tif 000009
【図５】
2010122658.tif 000010
【図６】
2010122658.tif 000011
【図７】
2010122658.tif 000012
【図８】
2010122658.tif 000013
【図９】
2010122658.tif 000014
【図１０】
2010122658.tif 000015
【図１１Ａ】
2010122658.tif 000016
【図１１Ｂ】
2010122658.tif 000017
【図１２Ａ】
2010122658.tif 000018
【図１２Ｂ】
2010122658.tif 000019
【手続補正書】
【提出日】平成21年12月3日(2009.12.3)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００５４
【補正方法】変更
【補正の内容】
【００５４】
例えば、上記歪εをセンサの位置におけるモーメントに変換するには、以下の（４）式が用いられる。なお、ｄはセンサ７の設置位置（翼根部から１．８メートル）における風車ブレード１０の内径、Ｌは風車ブレード１０の番号（Ｌ＝１、２、３）、Ｅは翼材（ＦＲＰ）のヤング率、Ｉはセンサの設置位置における断面二次モーメント、Ｍｓｅｎｓｏｒはセンサ設置位置における曲げモーメント（荷重）、ε２Ｌ−１とε２Ｌとは一対のセンサ（第１センサまたは第２センサ）の計測データに基づく歪、ε２Ｌ−１，０とε２Ｌ，０とは第１センサまたは第２センサの歪の初期値である。
【数１】
2010122658.tif 000020
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００５６
【補正方法】変更
【補正の内容】
【００５６】
【数２】
2010122658.tif 000021
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００６６
【補正方法】変更
【補正の内容】
【００６６】
（５）´式において、εＡ１は（２）´に基づいて算出されたセンサＡ１の位置における歪、εＡ３は同様の算出式に基づいて算出されたセンサＡ３の位置における歪である。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００７１
【補正方法】変更
【補正の内容】
【００７１】
図１１Ａは横軸に無風時のフラップ方向のモーメント、縦軸に風速８メートルの場合のフラップ方向のモーメントを示している。図１１Ｂは同様にエッジ方向のモーメントを比較している。図１１Ａ及び図１１Ｂは、図６の領域Ａを含む計測データによって係数ａ、ｂを求めた場合に得られる関数であり、理想的な無風時のデータと比較した場合である。理想的な無風時の関数をｙ＝ｘであるのに対し、領域Ａを含む計測データを用いて係数ａ及びｂを算出すると、ｙ＝０．９７０１ｘ−３１．８８となる。また、無風時のモーメントと領域Ａを含む計測データから算出されるモーメントとを比較すると、ｙ＝ｘからずれた位置に点が多く存在している。これにより、領域Ａの計測データを含んで算出された係数ａ、ｂを用いて得られる荷重は、無風時の荷重と比較して誤差が大きくなることを示している。
【国際調査報告】
2010122658.tif 000022
2010122658.tif 000023