JP7346462B2 - ターゲット物体との関係において複数のトランスデューサの空間的構成を決定する方法 - Google Patents

Description

一態様において、本開示は、ターゲット物体との関係における複数のトランスデューサ、具体的には、加速度計、の空間的構成を決定する方法に関し、この場合に、複数のトランスデューサは、ターゲット物体に装着されている。

構造解析や振動計測などのような、多くの用途においては、通常、例えば、接着剤又は適切な取付装置を使用することにより、加速度計などの複数のトランスデューサがターゲット物体に装着されている。通常、有用な計測結果を促進するべく、妥当な精度により、ターゲット物体との関係におけるそれぞれのトランスデューサの位置及び空間的向きを確立する必要がある。具体的には、ターゲット物体の、例えば、ＣＡＤモデル、ＦＥモデル、又はこれらに類似したものなどの、デジタルモデルとの関係においてトランスデューサの空間的構成を確立することが望ましい場合がある。

本説明を目的として、ボディの、具体的には、トランスデューサの、「空間的構成」という用語は、適切な基準座標系との関係におけるボディの位置及び向きを意味している。

大きな計測セットアップの複数のトランスデューサの、空間的構成、即ち、位置及び向き、のチェック及び検証は、構造解析や振動計測などのような、多数のタイプの用途において、大きな課題として残されている。ターゲット物体の寸法が大きくなるのに伴って、ターゲット物体上の計測場所の数も増大する。これには、通常、計測セットアップ内に含まれるトランスデューサの数の対応する増大が伴っている。更には、多くのターゲット物体は、トランスデューサの空間的構成の決定を困難にする複雑な表面構造を含む、複雑な形状を有する。

この結果、誤ったトランスデューサの配置や誤った空間的構成などのような、計測セットアップにおける誤りや間違いの機会が相対的に迅速に増大しうる。最近、不揮発性の半導体メモリを有する、所謂「スマートセンサ」が進歩を遂げている。半導体メモリは、連番、較正値、及び場所などの、様々なタイプの有用なトランスデューサ情報をＩＥＥＥ－Ｐ１４５１．４規格において定義されている標準化フォーマットにおいて電子的に保存する能力を有する。そのフォーマットは、現在、ＴＥＤＳ（Ｔｒａｎｓｄｕｃｅｒ Ｅｌｅｃｔｒｏｎｉｃ Ｄａｔａ Ｓｈｅｅｔ）と表記されている。それぞれのＴＥＤＳ準拠トランスデューサは、標準化通信プトロコルを介して、その保存されているトランスデューサ情報をリモート計測システム又は機器に送信する能力を有する。リモート計測システムは、自動的に、トランスデューサ情報を計測装置のセットアップ記述又はファイルに直接的に読み込むことができる。従って、この特徴は、トランスデューサデータの計測システムへの手動入力と関連したヒューマンエラーを低減しうる。

但し、このような計測セットアップにおけるそれぞれのトランスデューサの位置及び／又は空間的向きのチェック及び検証は、大きな課題として残されている。トランスデューサの適切な空間的構成は、例えば、加速度などの意図された又は少なくとも既知である成分が、意図された又は少なくとも既知である位置において、トランスデューサによって計測されることを保証するべく、多数のタイプの振動計測において重要である。

特許文献１は、ハンドヘルド型の光学スキャニング装置により、トランスデューサの空間的向きを検出する方法を開示している。これを目的として、この従来技術の方法は、トランスデューサの外側ハウジング表面上において空間的向き特徴を使用している。

上述の従来技術の方法は、トランスデューサの空間的向きを決定するべく必要とされる時間及び努力を相当に低減している一方で、この方法は、依然として、それぞれのトランスデューサの手動的な決定を必要としている。更には、しばしば、位置及び向き座標の両方を取得することが望ましい。更には、例えば、空間的向き特徴をスキャニングすることが困難である、小さなオリフィス内などの、到達することが困難であるターゲット物体上の位置においてトランスデューサを位置決めすることが望ましい場合もある。

国際公開第２０１６／１３５１９８号

従って、時間及び費用効率に優れた方式により、且つ、十分に高い精度を伴って、トランスデューサが装着されているターゲット物体との関係におけるトランスデューサの相対的に大きな組の空間的構成、即ち、位置及び向き、の決定を許容する、改善された方法に対するニーズが依然として残されている。

第１の態様によれば、本明細書において開示されているのは、ターゲット物体に装着された複数のトランスデューサの、それぞれの位置及び向きを示すそれぞれの空間的構成を決定するコンピュータ実装された方法の実施形態であり、この方法は、
－複数のトランスデューサのそれぞれからセンサ信号を受け取るステップであって、センサ信号は、ターゲット物体の誘発された運動の際のそれぞれのトランスデューサのそれぞれの運動を示している、ステップと、
－複数のトランスデューサの空間的構成を決定するステップであって、トランスデューサの少なくともサブセットの、ターゲット物体との関係における位置及び向きを示す空間的構成が、受け取られたセンサ信号に基づいて決定される、ステップと、
を有する。

本発明者らは、装着されているターゲット物体との関係におけるその位置及び向きを含む、複数のトランスデューサの空間的構成が、複数のトランスデューサ、並びに、恐らくは、ターゲット物体に装着されたその他のトランスデューサ、からのセンサ信号に基づいて、十分高精度において効率的に決定されうることを認識した。その結果、ボティに装着されたすべてのトランスデューサの位置及び向きを手動的に計測するニーズが、除去される、或いは、少なくとも大幅に低減される。その結果、特に、トランスデューサのサブセットが多数のトランスデューサを含んでいる際に、すべてのトランスデューサの空間的構成を確立するべく必要とされる時間及び努力が大幅に低減される。

好適な実施形態において、方法は、複数のトランスデューサの空間的構成の一部分のみを表す情報を入力として受け取っている。従って、複数のトランスデューサのうちのすべてのトランスデューサの完全な空間的構成を確立するべく、プロセスが決定する必要があるのは、空間的構成のうちの残りの部分のみである。前記情報は、トランスデューサの第１サブセットの空間的構成、具体的には、第１サブセットのトランスデューサの位置及び向きの両方を定義する情報、を含みうる。これを目的として、情報は、トランスデューサの第１サブセットのそれぞれのトランスデューサごとの６つの自由度（３つの位置及び３つの回転自由度）に沿った座標を表す情報を含みうる。或いは、この代わりに、又はこれに加えて、情報は、少なくとも第１サブセットのトランスデューサの間の空間的関係を含むこともできる。空間的関係は、トランスデューサのペアの間における相対距離、例えば、トランスデューサの計測座標系の１つの軸の向きなどのトランスデューサの向きに関する部分的情報、或いは、１つ又は２つ又は３つの自由度に沿った座標などの６つ未満の自由度に沿った座標を表すその他の部分的な空間的構成であってよい。従って、複数のトランスデューサの空間的構成の一部分のみを表す、受け取られた情報は、トランスデューサの第１サブセットのみの完全な空間的構成及び／又はトランスデューサの少なくとも第１サブセットの部分的な空間的構成を表しうる。

従って、複数のトランスデューサは、トランスデューサの第１サブセットと第２サブセットに分割されてもよく、この場合に、第２サブセットは、第１サブセットとは異なっており、これらのサブセットは、具体的には、オーバーラップしていなくてもよい。従って、いくつかの実施形態によれば、複数のトランスデューサの空間的構成を決定するステップは、
－第１サブセットのトランスデューサの空間的構成を取得するステップと、
－受け取られたセンサ信号から、且つ、第１サブセットのトランスデューサの取得された空間的構成から、第２サブセットのトランスデューサの空間的構成を決定するステップと、
を有することができる。

従って、トランスデューサの第１サブセットの空間的構成が判明したら、それらのそれぞれのセンサ信号から、且つ、第１サブセットの既知の空間的構成から、トランスデューサの第２サブセットの空間的構成を自動的に決定することができる。その結果、すべてのトランスデューサの空間的構成を確立するべく必要とされる時間及び努力が大幅に低減され、その理由は、具体的には、トランスデューサの第１サブセットが、トランスデューサの第２サブセットよりも、格段に小さくなりうるからである。例えば、１００個以上のトランスデューサがターゲット物体に装着される場合に、例えば、３つのトランスデューサのみ、などの、わずかな数のトランスデューサを第１サブセットに含むことで十分となりうる。これらのトランスデューサの空間的構成が、例えば、従来の方式により、手動的に確立されたら、或いは、第１サブセットのトランスデューサの既知の空間的関係から、且つ、受け取られたセンサ信号から、演算されたら、すべての残りのトランスデューサの空間的構成を自動的に決定することができる。本説明を目的として、第１サブセットのトランスデューサは、基準トランスデューサとも呼称することとする。一般に、ターゲット物体に装着されたトランスデューサの空間的構成は、トランスデューサの計測座標系とターゲット物体のローカル座標系の間の変換として表現されうる。或いは、この代わりに、ターゲット物体に装着されたトランスデューサの空間的構成は、トランスデューサの計測座標系とグローバル座標系の間の変換として表現することもできる。トランスデューサの計測座標系は、加速度がそれに沿ってトランスデューサによって計測される軸を定義しうる。

それぞれのトランスデューサは、単一軸の、或いは、例えば、３軸などの、多軸の、加速度計であってよく、或いは、これを有しうる。トランスデューサからのセンサ信号は、前記トランスデューサによる１つ又は複数の一連の計測値を表しうる。例えば、それぞれの計測値は、例えば、３つの方向におけるものなどの、１つ又は複数の方向に沿った加速度を示しうる。３軸加速度計は、時系列の計測値を提供してもよく、この場合に、それぞれの計測値は、加速度計の計測座標系との関係における加速度ベクトルを表している。トランスデューサのその他の例は、歪ゲージ、近接性プローブ、及び傾斜計、或いは、運動、加速度、力、又はターゲット物体の運動が導出されうるその他の量を計測するべく動作可能なその他の計測装置を含む。

第１サブセットのトランスデューサの空間的構成は、例えば、手動的に決定及び入力される、などのように、任意の適切な方式により、取得することができる。例えば、第１サブセットのトランスデューサは、例えば、トランスデューサを配置するための「係留点」として使用されうるターゲット物体の形状のいくつかの別個の特徴において、などのように、ターゲット物体の表面上の容易に識別可能な位置においてターゲット物体に取り付けることができる。トランスデューサは、重力の方向との関係においてトランスデューサをアライメントするように、スピリットレベル及び／又はスイベルベースを使用することにより、前記位置において取り付けられてもよい。このようにして識別された位置及び向きは、例えば、ターゲット物体の、例えば、ＣＡＤモデルなどの、適切なデジタル３Ｄモデルを使用することにより、適切な計測システムに手動的に入力することができる。従って、いくつかの実施形態において、第１サブセットのトランスデューサの空間的構成を取得するステップは、第１サブセットのトランスデューサの空間的構成を入力として受け取るステップを有する。その他の実施形態においては、プロセスは、第１サブセットのトランスデューサの正確な空間的構成を受け取ってはおらず、第１サブセットのトランスデューサの間の空間的関係の組を示す情報のみを受け取っている。次いで、プロセスは、受け取られた情報から、且つ、受け取られたセンサ信号から、第１サブセットのトランスデューサの空間的構成を演算している。

好ましくは、第１サブセットは、すべてのその自由度に沿ってターゲット物体の運動を計測するべく、十分な数のトランスデューサを有する。３つの直交方向に沿った平行運動と、３つの直交軸の周りにおける回転運動と、を実行しうるように支持されたターゲット物体の場合には、トランスデューサの数は、６つの自由度、即ち、３つの平行運動自由度及び３つの回転自由度、に沿ってターゲット物体の運動を計測するのに十分なものとなるように選択することができる。必要とされるトランスデューサの正確な数は、トランスデューサの能力と、ターゲット物体との関係におけるその相対的位置と、ターゲット物体の運動の自由度の数と、に依存している。例えば、トランスデューサが３軸加速度計である際には、３つのトランスデューサがまっすぐなラインに沿って位置決めされていない際には、第１サブセットを形成するべく、３つのトランスデューサで十分である。又、好ましくは、３つのトランスデューサは、互いに十分に離れるように、或いは、ターゲット物体の形状が許容する範囲において、少なくとも互いに離れるように、位置決めする必要がある。

いくつかの実施形態において、方法は、第１サブセットが、十分な数のトランスデューサを含んでいるかどうか、並びに／或いは、トランスデューサが十分に位置決めされているかどうか、を検証するステップを有する。例えば、それぞれのトランスデューサが３軸加速度計を有するいくつかの実施形態において、それぞれのトランスデューサは、トランスデューサの計測座標系の軸に沿って前記トランスデューサによって計測される加速度に対するターゲット物体のローカル座標との関係におけるトランスデューサの加速度、具体的には、平行運動及び回転運動加速度、の間の変換を表すトランスデューサ伝達率行列と関連付けることができる。一般に、ターゲット物体に装着されたトランスデューサのトランスデューサ伝達率行列は、ターゲット物体のローカル座標系との関係におけるトランスデューサの計測座標系の位置及び向きを表している。従って、ターゲット物体との関係におけるトランスデューサの既知の位置及び向きの場合には、前記既知の位置及び向きから、トランスデューサ伝達率行列を決定することができる。その結果、第１サブセットが、十分な数のトランスデューサを含んでいるかどうかを検証するステップは、第１サブセットのトランスデューサのトランスデューサ変換行列から構築された変換行列が、フルランクを有しているかどうかを検証するステップを有しうる。変換行列がフルランクを有している場合には、プロセスは、第１サブセットが十分な数のトランスデューサを含んでいると決定することができる。

いくつかの実施形態において、プロセスは、誘発された運動が、トランスデューサ、即ち、トランスデューサの第２サブセットのトランスデューサ、の空間的構成を決定するために十分なものであるかどうかを決定するステップを有する。これを目的として、それぞれのトランスデューサが３軸加速度計を有するいくつかの実施形態において、生成された運動が、トランスデューサの第２サブセットのトランスデューサの空間的構成を決定するために十分であるかどうかを決定するステップは、
－誘発された運動の際のトランスデューサの計測された加速度を示すセンサ信号を受け取るステップと、
－任意選択により、センサ信号の低域通過フィルタリングを実行するステップと、
－計測された加速度を表すローカル加速度行列が、既定の数の主要特異値、即ち、制限されていない運動に対応する６つの主要特異値、を有するかどうを決定するステップと、
を有する。この決定は、最大の特異値（例えば、最大の６つの主要値）を残りの特異値のうちの最大のものと比較することにより、実施することができる。

ローカル加速度行列が、十分な数（例えば、６つ）の主要特異値を有している場合に、生成された運動は、トランスデューサの第２サブセットのトランスデューサの空間的構成を決定するために十分なものであると決定される。

一般に、いくつかの実施形態においては、特に、ターゲット物体がすべての６つの自由度に沿って運動しうる特定の実施形態においては、方法は、誘発された運動の撓みモードと関連する周波数を抑制するように、記録された信号をフィルタリングするステップを有する。具体的には、いくつかの実施形態において、方法は、センサ信号の低域通過フィルタリングステップを有する。これを目的として、プロセスは、最低撓みモードの固有周波数と、最高剛体モードの固有周波数と、を検出するステップと、撓みモードの少なくとも過半数を抑圧するように、且つ、剛体モードの少なくとも過半数を維持するように、低域通過フィルタのカットオフ周波数を選択するステップと、を有しうる。これを目的として、カットオフ周波数は、検出された最高剛体モードの固有周波数よりも大きくなるように、且つ、検出された最低撓みモードの固有周波数よりも小さくなるように、選択することができる。

トランスデューサの第２サブセットのトランスデューサの空間的構成を決定するステップは、
－第１サブセットのトランスデューサの受け取られたセンサ信号及び第１サブセットのトランスデューサの取得された空間的構成から基準座標系との関係におけるターゲット物体の加速度を示す物体加速度を演算するステップと、
－演算された物体加速度から、且つ、第２サブセットのそれぞれのトランスデューサの受け取られたセンサ信号から、第２サブセットのトランスデューサの前記空間的構成を演算するステップと、
を有することができる。

具体的には、物体加速度は、慣性基準座標系との関係におけるターゲット物体のローカル座標系の、時間との関係において２回にわたって微分された空間的構成を表しうる。

物体加速度を演算するステップは、一次方程式の系に対する最小二乗解を演算するステップを有していてもよく、この場合に、それぞれの一次方程式は、基準加速度計のうちの１つによって計測された、計測された加速度と未知の物体加速度の間の変換を表している。

上述のように、いくつかの実施形態において、プロセスは、第１サブセットのトランスデューサの間の空間的関係の組を示す情報を受け取っており、且つ、プロセスは、受け取られた情報から、且つ、受け取られたセンサ信号から、第１サブセットのトランスデューサの空間的構成を演算している。第１サブセットのトランスデューサの間の空間的関係に関する情報は、第１サブセットのトランスデューサのそれぞれのペアの間の距離を含みうる。或いは、この代わりに、又はこれに加えて、第１サブセットのトランスデューサの間の空間的関係に関する情報は、相対距離及び向きの組合せ、或いは、空間的パラメータの別の組合せ、を含みうる。

第１サブセットのトランスデューサの空間的構成を演算するステップは、
－トランスデューサの第１サブセットからの受け取られたセンサ信号からのセンサ信号の第１の組及びトランスデューサの第１のサブセットからの前記受け取られたセンサ信号からのセンサ信号の第２の組を定義するステップと、
－センサ信号の第１の組のセンサ信号から、且つ、第１サブセットのトランスデューサの間の空間的関係を示す受け取られた情報から、センサ信号の第２の組のセンサ信号の推定値を演算するステップであって、センサ信号の推定値は、第１サブセットのトランスデューサの推定された空間的構成に依存している、ステップと、
－演算された推定値とセンサ信号の第１の組の間の誤差を低減する、具体的には、極小化する、ように、推定された空間的構成を調節するステップと、
を含みうる。

従って、第１の組のトランスデューサの空間的構成は、演算された空間的構成が、受け取られたセンサ信号の第１の組が付与された場合に、受け取られたセンサ信号の第２の組をどれだけ良好に表しているかを示す極小化問題から、演算することができる。いくつかの実施形態において、第１サブセットのトランスデューサの空間的構成の演算は、例えば、重力の方向に沿って、などのように、既知の方向に沿って指し示す１つの軸に伴って、第１サブセットのトランスデューサを取り付けることにより、更に改善することができる。

一般に、本明細書において記述されている方法のいくつかの実施形態は、すべての６つの自由度におけるターゲット物体の運動を許容するように、吊り下げられた、或いは、その他の方法で支持された、ターゲット物体との関係において実行することができる。

本明細書において開示されている方法のその他の実施形態は、物体支持部により、（例えば、１つ又は複数の自由度に沿った）そのモーションが制限されている、変形可能なターゲット物体との関係において実行することができる。これを目的として、方法のいくつかの実施形態は、
－第１サブセットのトランスデューサの空間的構成から、且つ、受け取られたセンサ信号から、ターゲット物体の時間依存加速度係数を演算するステップと、
－演算された時間依存加速度係数から、且つ、受け取られたセンサ信号から、例えば、推定されたセンサ信号と受け取られたセンサ信号の間の差を示すフィット誤りを極小化することにより、第２サブセットのトランスデューの空間的構成を演算するステップと、
を有する。

時間依存加速度係数は、ターゲット物体のモード座標と関連する係数であってよい。このケースにおいて、係数を演算するための手順は、例えば、モード分解におけるモード座標の取得に類似しうる。或いは、この代わりに、時間依存加速度係数は、基底関数の適切な組に基づいた別の拡張の係数であってもよい。基底関数の適切な組は、好ましくは、例えば、所謂基本境界条件などの、境界条件の組を充足することを要し、且つ、互いに直交することを要する。

上述の、且つ、以下の、コンピュータ実装された方法の特徴は、ソフトウェア又はファームウェアにおいて少なくとも部分的に実装することができると共に、コンピュータ実行可能命令などのプログラムコード手段の実行によってもたらされるデータ処理システム又はその他の処理手段上において実行することができることに留意されたい。ここでは、且つ、以下においては、処理手段という用語は、上述の機能を実行するように適切に適合された任意の回路及び／又は装置を有する。具体的には、上述の用語は、汎用又は特殊目的プログラム可能マイクロプロセッサ、デジタル信号プロセッサ（ＤＳＰ：Ｄｉｇｉｔａｌ Ｓｉｇｎａｌ Ｐｒｏｃｅｓｓｏｒ）、用途固有の集積回路（ＡＳＩＣ：Ａｐｐｌｉｃａｔｉｏｎ Ｓｐｅｃｉｆｉｃ Ｉｎｔｅｇｒａｔｅｄ Ｃｉｒｃｕｉｔ）、プログラム可能なロジックアレイ（ＰＬＡ：Ｐｒｏｇｒａｍｍａｂｌｅ Ｌｏｇｉｃ Ａｒｒａｙ）、フィールドプログラム可能なゲートアレイ（ＦＰＧＡ：Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）、グラフィカル処理ユニット（ＧＰＵ：Ｇｒａｐｈｉｃａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、特殊目的電子回路など、或いは、これらの組合せを有する。

本開示は、上述の且つ以下における方法、それぞれがその他の態様の１つ又は複数との関連において記述されている利益及び利点のうちの１つ又は複数をもたらす、且つ、それぞれが本明細書において記述されているその他の態様の１つ又は複数との関連において記述された、且つ／又は、添付の請求項において開示された、実施形態に対応する１つ又は複数の実施形態を有する、更なる方法、システム、装置、及び製品手段を含む様々な態様に関する。

具体的には、本明細書において開示されている別の態様は、
－上述の、且つ、以下における、第１の態様による方法のステップを実行することにより、ターゲット物体に装着された複数のトランスデューサの空間的構成を決定するステップと、
－複数のトランスデューサから更なるセンサ信号を受け取るステップであって、受け取られた更なるセンサ信号は、ターゲット物体の誘発された運動の際のそれぞれのトランスデューサのそれぞれの運動を示している、ステップと、
－受け取られた更なるセンサ信号に基づいて構造解析計算を実行するステップと、
を有する、コンピュータ実装された計測プロセスの実施形態に関する。

更なるセンサ信号は、意図された構造解析計算を実行するのに適した運動に基づいて取得することができる。一般に、更なるセンサ信号は、トランスデューサの空間的構成が決定された誘発された運動とは異なる、ターゲット物体の更なる誘発された運動の際のそれぞれのトランスデューサの運動を示しうる。更なる誘発された運動は、剛体運動に加えて、或いは、このかわりに、撓み運動を含みうる。実行されるべく意図されている構造解析に応じて、更なる誘発された運動は、すべての６つの自由度に沿った剛体運動を含んでいてもよく、或いは、これを含んでいなくてもよい。

本明細書において開示されている更に別の態様は、ターゲット物体に装着された複数のトランスデューサの空間的構成を決定する方法の実施形態に関し、この場合に、複数のトランスデューサは、トランスデューサの第１の組及びトランスデューサの第２の組を有しており、この方法は、
－ターゲット物体の運動を誘発するステップと、
－トランスデューサの第１及び第２サブセットからセンサ信号を受け取るステップであって、受け取られたセンサ信号は、それぞれのトランスデューサのそれぞれの運動を示している、ステップと、
－第１サブセットのトランスデューサの空間的構成を取得するステップと、
－受け取られたセンサ信号及び第１サブセットのトランスデューサの取得された空間的構成から、第２サブセットのトランスデューの空間的構成を決定するステップと、
を有する。

いくつかの実施形態において、第１サブセットのトランスデューサの空間的構成は、ターゲット物体の運動を誘発する前に、且つ／又は、センサ信号を受け取る前に、取得することができる。但し、第１サブセットのトランスデューサの空間的構成が、受け取られたセンサ信号に基づいて決定される実施形態においては、第１サブセットのトランスデューサの空間的構成は、前記センサ信号の受け取りの後に取得されていることを理解されたい。

ターゲット物体の運動を誘発するステップは、例えば、ターゲット物体を手動で振り動かすことにより、且つ／又は、励起ハンマーによってターゲット物体を打撃することにより、ターゲット物体上に１つ又は複数の力を作用させるステップを有することができる。いくつかの実施形態において、運動は、すべての利用可能な自由度に沿った成分を含むことを要する。すべての６つの自由度に沿って運動可能であるターゲット物体の場合に、誘発された運動は、３つの直交方向に沿った成分と、３つの線形独立軸の周りの回転運動成分と、を含むことを要する。これを目的として、ターゲット物体は、物体が、３つの直交方向に沿った少なくとも小さな運動と、３つの線形独立軸の周りにおける少なくと小さな回転と、を実行しうるように、弾性バンドから吊り下げられていてもよく、空気クッション上において支持されていてもよく、或いは、その他の方法で支持されていてもよい。その他の実施形態において、ターゲット物体は、更に詳細に後述するように、その運動が１つ又は複数の自由度に沿って制約されるように、支持することができる。

本明細書において開示されている更に別の態様は、
－上述の、且つ、以下における、ターゲット物体に装着された複数のトランスデューサの空間的構成を決定する方法のステップを実行することにより、ターゲット物体に装着された複数のトランスデューサの空間的構成を決定するステップと、
－ターゲット物体の更なる運動を誘発するステップと、
－トランスデューサの第１及び第２サブセットから更なるセンサ信号を受け取るステップであって、受け取られた更なるセンサ信号は、前記誘発された更なる運動の際のそれぞれのトランスデューサのそれぞれの運動を示している、ステップと、
－受け取られた更なるセンサ信号に基づいて構造解析計算を実行するステップと、
を有する計測プロセスの実施形態に関する。

具体的には、トランスデューサの空間的構成を決定する、という目的のために、撓みモードを回避又は少なくとも抑圧することが有利でありうるが、構造解析がそれに基づいて実行される誘発された更なる運動は、抑圧される必要はないと共に実際には構造解析の目的のために興味深いものでありうる、大きな撓み運動を含むこともできる。

本明細書において開示されている更なる別の態様は、データ処理システムが、上述の、且つ、以下における、コンピュータ実装された方法のステップを実行するようにするべく構成されたコンピュータプログラムの実施形態に関する。コンピュータプログラムは、プログラムコードがデータ処理システム上において実行された際に、データ処理システムが、上述の、且つ、以下における、コンピュータ実装された方法のステップを実行するようにするべく適合されたプログラムコードを有しうる。コンピュータプログラムは、コンピュータ可読ストレージ媒体上において保存されていてもよく、或いは、データ信号として実施されてもよい。ストレージ媒体は、ＲＡＭ、ＲＯＭ、ＥＰＲＯＭ、ＥＥＰＲＯＭ、フラッシュメモリ、ＣＤ－ＲＯＭ、ＤＶＤ、ハードディスクなどの、磁気又は光学ストレージ装置、及び／又はこれらに類似したものなどの、データを保存するための任意の適切な回路又は装置を有しうる。

本明細書において開示されている、更に別の態様は、本明細書において開示されているコンピュータ実装された方法のステップを実行するように構成されたデータ処理システムを有する計測システムの実施形態に関する。計測システムは、
－それぞれの無線又は有線信号接続を介して、ターゲット物体に装着された複数のトランスデューサのそれぞれのトランスデューサ出力信号に結合されるように構成された複数の入力チャネルと、
－例えば、手動的なデータ入力による、或いは、無線又は有線データ通信リンクを介した別の装置からの、第１サブセットのトランスデューサのそれぞれの空間的座標データ又は少なくとも第１サブセットのトランスデューサの間の空間的関係に関する情報の受け取りのために構成されたデータ入力インターフェイスと、
を更に有しうる。

計測システムは、複数の信号入力チャネルに結合された複数のトランスデューサの決定された空間的構成を示すディスプレイを有することができる。

本明細書において開示されている更なる別の態様は、ターゲット物体に跨って分散された複数の既定の計測場所において取り付け可能である複数のトランスデューサを有する計測組立体の実施形態に関する。計測組立体は、上述の、且つ、以下における、計測システムを更に有する。

上述の且つその他の態様については、以下の図面を参照して以下において記述されている実施形態から、明らかとなり、且つ、解明されることになろう。

計測組立体の概略ブロック図を示す。 計測組立体のトランスデューサの概略図を示す。 本説明を目的として使用される座標系の例を示す。 本明細書において開示されている計測プロセスの一例のフロー図を示す。 本明細書において開示されているトランスデューサの空間的構成を決定するコンピュータ実装されたプロセスの一例のフロー図を示す。 例示用のターゲット物体の計測された周波数応答関数を示す。 例示用のターゲット物体の計測された周波数応答関数を示す。 ターゲット物体の例示用の運動について決定された正規化された特異値の例を示す。 それぞれの加速度計の決定された空間的構成を示すターゲット物体の表示された３Ｄデジタルモデルを概略的に示す。 それぞれの加速度計の決定された空間的構成を示すターゲット物体の表示された３Ｄデジタルモデルを概略的に示す。 本明細書において開示されている、第２サブセットのトランスデューサ、具体的には、３軸加速度計、の空間的構成を決定するコンピュータ実装されたプロセスの一例のフロー図を示す。 そのモーションが制約されている、変形可能なターゲット物体の一例を示す。 その上部において取り付けられたトランスデューサを有するターゲット物体の一例の図を示す。 本明細書において開示されているプロセスの一例によって決定される、決定された空間的構成の位置及び向きの誤りを示す。 本明細書において開示されているプロセスの一例によって決定される、決定された空間的構成の位置及び向きの誤りを示す。

図１Ａ～図１Ｃは、本明細書において開示されている計測組立体の一例を概略的に示している。具体的には、図１Ａは、計測組立体の概略ブロック図を示しており、図１Ｂは、計測組立体のトランスデューサの概略図を示しており、且つ、図１Ｃは、本説明を目的として使用される座標系の例を示している。

計測組立体は、計測システム１１０と、ターゲット物体１２０に装着された、且つ、計測装置１１０に通信自在に結合された、複数の３軸加速度計１００と、を有する。

３軸加速度計１００は、例えば、いずれも、Ｂｒｕｅｌ ａｎｄ Ｋｊｅｒ Ｓｏｕｎｄ ａｎｄ Ｖｉｂｒａｔｉｏｎ Ｍｅａｓｕｒｅｍｅｎｔ， Ｎｅｒｕｍ， Ｄｅｎｍａｒｋによって製造された４５２４、４５２４－Ｂ、又は４５０４型の３軸圧電加速度計を有することができる。当業者は、例えば、その他のタイプの加速度計などの、その他のタイプのトランスデューサが代替実施形態において使用されうることを理解するであろう。例えば、当然のことながら、代替肢においては、本発明を実装するべく、多数のその他の製造者からのその他のタイプの単軸又は多軸加速度計を使用することができる。加速度計の特定の構造的な詳細に応じて、１つ又は複数の圧電トランスデューサ要素１０１は、トランスデューサハウジング１０７の内側において取り付けられてもよく、この場合に、その特定の方向における加速度に比例した電圧又は電流などの、センサ信号を生成することにより、１つ又は複数のトランスデューサ要素は、加速度計との関係において３つの直交する空間的方向（本説明を目的として、ｉ、ｊ、ｋ方向と表記されている）の特定の方向における加速度に応答している。従って、方向ｉ、ｊ、ｋは、加速度計の計測座標系を定義している。同様に、低ノイズ前置増幅器、フィルタ、Ａ／Ｄコンバータ、電源などのような、様々なタイプの電子信号コンディショニング回路及び／又はメモリ装置が、適切な担持体上において、トランスデューサハウジング１０７の内側に取り付けられうる。電子信号コンディショニング回路は、それぞれの感度方向について、加速度計１００の低インピーダンス及び恐らくは周波数成形された出力信号を提供するべく、１つ又は複数の圧電トランスデューサ要素のそれぞれの出力端子に結合されうる。加速度計ハウジング１０７は、例えば、湿度、機械的衝撃、埃、光、及び熱などのような様々な有害な環境汚染物質から、１つ又は複数のトランスデューサ要素を保護するべく、チタニウムなどの金属性の組成又は材料を有しうる。又、金属性の加速度計ハウジング１０７は、ＥＭＩ遮蔽目的のためにも、有用でありうる。加速度計ハウジング１０７は、いくつかの異なるターゲット物体への加速度計１００の装着又はフィッティングのために、取付クリップを支持する、いくつかの任意選択のスロット１０９ａ、１０９ｂを有する。加速度計ハウジング１０７の、下部の実質的に平坦な外側ハウジング表面（図においては、隠蔽されており―外側ハウジング表面１０１とは反対側に配置されている）は、直接又は間接的に、（例えば、適切な接着剤又は接着層を介して）ターゲット物体との物理的接触状態となっている。従って、加速度計ハウジング１０７の、下部の実質的に平坦な外側ハウジング表面は、ターゲット物体に対する接合又は結合表面として機能している。

３軸加速度計１００は、共通接地端子と、それぞれ、３軸加速度計１００の１つ又は複数のトランスデューサ要素の感度の対応する第１、第２、及び第３直交軸を表す加速度計ｉ、ｊ、ｋ成分出力信号を搬送する３つの別個の出力信号端子と、を有する４ピンコネクタを有しうる電気コネクタ１０２を更に有する。

ｉ成分出力信号は、第１平坦外側ハウジング表面１０１を通じて直交するように突出する加速度計の予め定義されたｉ方向における加速度計１００の加速度を示している。同一の内容が、それぞれ、関連するターゲット物体の振動の真の３軸加速度計測値を提供するべく、ｊ軸及びｋ軸に沿った加速度を示すｊ及びｋ成分出力信号についても適用される。ターゲット物体１２０は、例えば、自動車ボディ、航空機構造、列車構造、風力タービンブレード、衛星構造、エンジン、トランスミッションなど、或いは、これらの部品を有しうる。電気コネクタ１０２は、例えば、低ノイズの遮蔽ケーブルなどの適切なケーブルを介して、ｉ、ｊ、及びｋ成分出力信号を計測システム１１０に結合するべく、使用することができる。当業者は、実際の計測セットアップは、適切な電気配線を介して計測システムに結合された、例えば、２０個超、或いは、場合によっては、１００個超、などの、多数の、それぞれの３軸加速度計１００を含みうることを理解するであろう。

計測システム１１０は、パーソナルコンピュータ又はその他の演算ハードウェアプラットフォーム上において実行される、様々なタイプのデータ取得ソフトウェアとの組合せにおいて、適切な加速度計計装システムを有しうる。具体的には、計測システムは、データ取得インターフェイス１１３、データ処理ユニット１１１、メモリ１１２、又はその他のデータストレージ装置、並びに、ユーザーインターフェイス１１４を有することができる。

データ取得インターフェイス１１３は、多数の３軸加速度計のそれぞれが、計測システムの特定の計測チャネルに結合されることを許容するべく、コネクタを含みうる。

データ処理ユニット１１１は、例えば、コンピュータ又はその他のデータ処理システムの適切にプログラミングされた中央処理ユニットであってよい。データ処理ユニット１１１は、ターゲット物体との関係において、且つ、加速度計１００から受け取られたセンサ信号に基づいて、構造解析計算を実行する、構造解析ソフトウェアアプリケーションを実行するように更に構成することもできる。

ユーザーインターフェイス１１４は、計測チャネルのそれぞれのものの状態及び識別子のみならず、例えば、連番などの、それぞれの加速度計の識別及び状態情報をグラフィカルに描くように構成されたディスプレイを含むことができる。ディスプレイは、加速度計の決定された位置及び向きを示すターゲット物体のデジタル３Ｄモデルの表示を更に許容することもできる。

計測システムは、３軸加速度計１００のそれぞれのものの決定された空間的６次元位置及び向きデータ、並びに、任意選択により、その連番、及び／又は、ターゲット物体上における３軸加速度計１００の加速度の、それぞれ、ｉ、ｊ、ｋ成分を表す、それぞれの加速度計のｉ、ｊ、ｋ出力信号のそれぞれのものに（４つのピンコネクタ１０２を通じて）接続された特定の計測チャネルを有するそれぞれの加速度計に属するその他の情報を関連付け又はリンクしている。計測システムは、複数の別個のコンポーネントを含みうることを理解されたい。コンポーネントは、データ転送を許容するように、互いに通信自在に結合されていてもよい。例えば、計測システムは、トランスデューサからのセンサ信号を記録するように構成されたデータ取得システムと、記録されたセンタ信号を分析するように構成されたデータ処理システムと、を含みうる。

ターゲット物体との関係におけるトランスデューサの空間的構成を決定するプロセスのいくつかの実施形態に関する以下の説明を目的として、ターゲット物体１２０は、３Ｄ空間内において制限されてはいない状態において運動しうる剛体Ωとして見なすことができる。ボディの位置及び向きは、慣性グローバル座標系（ＧＣＳ：Ｇｌｏｂａｌ Ｃｏｏｒｄｉｎａｔｅ Ｓｙｓｔｅｍ）において定義されている。ボディ上において固定されたローカル座標系（ＬＣＳ：Ｌｏｃａｌ Ｃｏｏｒｄｉｎａｔｅ Ｓｙｓｔｅｍ）を定義しよう。ポイントＰ_nにおいてボディに装着されたＮ個の３軸加速度計１００を検討すれば、ＬＣＳにおけるこれらの座標は、

である。図示の容易性を目的として、図１Ｃにおいては、単一の加速度計のみが明示的に示されているが、加速度計の数Ｎは、通常、１よりも大きく、具体的には、３よりも大きく、且つ、しばしば、１０よりも大きく、或いは、場合によっては、１０００よりも大きくなることを理解されたい。

それぞれの加速度計１００は、（慣性）ＧＣＳとの関係において適切な加速度を計測しており、加速度ベクトルは、加速度計の計測軸に沿って計測された、その３つのスカラー成分を介して提供される。次いで、それぞれの加速度計の計測軸と関連する、Ｎ個の計測座標系（ＭＣＳ：Ｍｅａｓｕｒｅｍｅｎｔ Ｃｏｏｒｄｉｎａｔｅ Ｓｙｓｔｅｍ）を定義しよう。ＭＣＳの位置及び向きは、ＬＣＳ内において固定されている。

ＬＣＳがターゲット物体上において固定されていることから、ＬＣＳの位置及び向きは、ターゲット物体が剛体として振る舞う際には、ターゲット物体の位置及び向きを十分に特徴付けている。任意の時点ｔにおいて、ＧＣＳにおいて付与された状態において、ベクトル

であり、

は、ＧＣＳ内のＬＣＳの位置及び向きを完全に記述しており、即ち、このベクトルは、ターゲット物体と関連するローカル座標系の空間的構成を記述している。ここで、

は、ＬＣＳの原点の半径ベクトルであり、且つ、θ₁（ｔ）；θ₂（ｔ）；θ₃（ｔ）は、ＧＣＳをＬＣＳに移行させる連続的な回転を定義する、３つのオイラー角である。

オイラー角の異なる定式化が存在しており、本明細書においては、我々は、Ｔａｉｔ－Ｂｒｙａｎのｘ－ｙ’－ｘ’’表記法（航海又はカルダン角とも呼称される：ヨー、ピッチ、及びロール）を使用している。ＧＣＳ内のポイントＰの座標が、｛ｒ｝_Gであり、且つ、ＬＣＳ内においては、｛ｒ_L｝であり、且つ、ＧＳＣ及びＬＣＳ原点が一致しているとしよう。その結果、座標｛ｒ_L｝は、まずは、θ₁によるＧＣＳのｚ軸の周り、次いで、θ₂による新しいｙ軸の周り、そして、最後に、θ₃による新しいｘ軸の周り、という３つの連続的な回転として取得することが可能であり、

ここで、それぞれの回転は、回転行列による半径ベクトルの事前乗算により、表されている。組み合わせられた３つの回転行列は、単一の回転行列｛Ｒ_GL｝をもたらす。ＬＣＳからＧＳＣへの逆座標変換は、

として提示することが可能であり、この場合に、以下のとおりである。

この場合の第２の式の理由は、回転行列が直交行列であるからである。

回転行列を使用することにより、ＧＳＣ内のポイントＰ_nにおいて配置された加速度計の座標は、その両方が時間依存性を有する、ＧＣＳ内におけるその平行運動

と回転の合計として提示することが可能である。

このポイントにおいて計測される加速度ベクトルは、以下のようになる。

ポイントＰ_nに装着された３軸加速度計は、加速度計と関連するＭＣＳの軸に投射された加速度ベクトル（５）に対応する、３つの信号

、

、及び

を提供しており、

ここで、｛ｉ_n（ｔ）｝_G、｛ｊ_n（ｔ）｝_G、及び｛ｋ_n（ｔ）｝_Gは、ＧＣＳ内において提示されているＭＣＳの場所（ｏｒｔ）である。ＧＣＳ内における場所の座標は、一連の、まずは、ＬＣＳへの、且つ、次いで、ＧＳＣへの、ＭＣＳの連続的な回転として取得することが可能であり、

ここで、

であり、且つ、

は、ポイントＰ_nにおいて取り付けられた加速度計と関連するＬＣＳからＭＣＳへの回転を表す行列である。

これは、オイラー角の組

、

、及び

を使用することにより、（２）と類似した方式により、構築することが可能であり、

ここで、

であり、且つ、

である。

（５）及び（７）を（６）に代入し、且つ、行列の形態を書き換えることにより、以下が得られる。

一般性の喪失を伴うことなしに、我々は、まず、ボディが休止状態にある際に、ＬＣＳは、ＧＣＳと一致している、と仮定することが可能であり（さもなければ、ＧＣＳをＬＣＳの初期場所まで運動させるべく、余分な平行運動及び回転を適用することが可能であり、且つ、次いで、逆変換を結果に適用することができる）、

である。

次に、その初期位置に関するボディの運動が小さいものと仮定しよう。

ここで、εは、付加された項の大きさのレベルを示すためのブックキーピング装置として使用されている。

（２）を（９）に代入し、且つ、２回にわたって微分することにより、結果がテイラー級数に展開され、その結果、以下が得られうる。

ε²以上の大きさのレベルを有する項を省略し、且つ、再構成することにより、以下が得られる。

（１２）における２つの行列は、時間不変３×６行列

に組み合わせることが可能であり、その要素は、ＬＣＳからＭＣＳへの回転を提供する３つのオイラー角を使用することにより、提示される、ＬＣＳ内において付与されたポイントＰ_nの座標の非線形関数と、加速度計の向きと、である。

これらのパラメータのいずれもが、小さいと仮定されてはいないことに留意されたい。

要すれば、小さなモーションのケースにおいて、ｎ番目の３軸加速度計によって計測される信号は、時間において２回にわたって微分された、行列［β_n］及びベクトル｛Ｃ（ｔ）｝の積として近似することができる。

式を短縮するべく、ベクトル

を導入すれば、これにより、式（１４）は、以下のようにコンパクトになる。

従って、一般的に、加速度計による、計測された加速度は、前記トランスデューサと関連するトランスデューサ伝達率行列から、且つ、ターゲット物体のローカル座標系の、時間との関係において２回にわたって微分された、空間的構成から、演算することができる。トランスデューサ伝達率行列は、ターゲット物体のローカル座標系との関係における前記加速度計の計測座標系の空間的構成を表している。

図２は、本明細書において開示されている計測プロセスの一例のフロー図を示している。計測プロセスは、例えば、図１Ａ～図１Ｃとの関連において記述されている組立体を使用することにより、実行することができる。

初期ステップＳ１０において、例えば、図１Ａ及び図１Ｃにおいて示されている物体１２０などの、構造解析の実行を要するターゲット物体が、試験サイトにおいて取り付けられている。取付のタイプは、その容積、質量、及び構造的完全性などの、ターゲット物体の特性に、且つ、例えば、実行されるべき構造解析のタイプなどの取得対象の計測値の目的に、依存していることを理解されたい。例えば、ターゲット物体が、ゴムバンジー及び／又は空気まくら上において取り付けられている際には、任意の方向における、その剛体モーションは、実質的に制限されてはいない状態において留まっている。但し、いくつかの実施形態においては、ターゲット物体の剛体モーションを部分的に制限する取付が望ましい場合もあることを理解されたい。

後続のステップＳ２０において、ターゲット物体の運動を記録しうるように、例えば、図１Ｂに示されている、３軸加速度計などの、加速度計などの複数のトランスデューサがターゲット物体に取り付けられている。例えば、トランスデューサは、例えば、適切な取付要素、接着剤、又はこれらに類似したものにより、ターゲット物体の表面に跨って分散されてもよく、且つ、これに装着されてもよい。トランスデューサは、ターゲット物体との関係におけるトランスデューサの第１サブセットの、空間的構成、即ち、位置及び向き、が、判明する、或いは、例えば、計測される、などのように、手動的に決定されるように、取り付けられている。例えば、第１サブセットのトランスデューサは、ターゲット物体の容易に認識可能な位置において位置決めされてもよく、且つ、例えば、それぞれのエッジ、コーナー、又は類似の表面特徴におけるもの、などのように、ターゲット物体の容易に認識可能な方向的特徴とアライメントされてもよい。基準トランスデューサの位置は、ＬＣＳとの関係における、その座標として提供されてもよい。ＬＣＳとの関係における基準トランスデューサの向きは、向きを定義するための任意の適切な表記法を使用することにより、提供されてもよい。例えば、これらは、ｚ－ｙ’－ｘ’’表記法に準拠した３つのＴａｉｔ－Ｂｒｙａｎ角として提供されてもよい。これらの角度は、基準トランスデューサのＬＣＳからＭＣＳへの連続的な回転に対応している。一般に、基準トランスデューサによって取り囲まれたエリアが、ターゲット物体の大部分をカバーするように、基準トランスデューサを位置決めすることが好ましくありうる。更に詳細に後述するように、いくつかの実施形態においては、その完全な空間的構成ではなく、基準トランスデューサの間の空間的関係に関する情報のみ（例えば、互いからのその相互距離）の入力が必要とされうる。

トランスデューサが３軸加速度計を有する際には、第１サブセットは、好ましくは、まっすぐなラインに沿って位置決めされてはいない、少なくとも３つの３軸加速度計を有する。本説明を目的として、第１のサブセットの加速度計は、基準速度計とも呼称することとする。

ステップＳ３０において、例えば、図１の実施形態の計測システム１１０などの、計測システムが初期化されている。これを目的として、トランスデューサは、計測システムのデータ取得インターフェイスに通信自在に接続されている。更には、基準トランスデューサの既知の空間的構成（或いは、少なくとも、その空間的関係に関する情報）が、計測システムに入力されてもよい。又、例えば、ＣＡＤモデル及びＦＥモデル、或いは、これらに類似したものなどの、ターゲット物体のデジタルモデルが、計測システムのメモリに読み込まれてもよい。最後に、例えば、調節可能な減衰、フィルタパラメータ、及び／又はこれらに類似したものなどの、計測システムによって実行されるべき信号処理の１つ又は複数の動作パラメータが設定されてもよい。例えば、まずは、取得されたセンサ信号に適用されるべき高域通過フィルタの値が、例えば、０．７Ｈｚなどの、１Ｈｚ未満のような、小さな値に設定されてもよい。適切な値は、特定のデータ取得システム及び／又は試験対象のターゲット物体の特性に依存しうることを理解されたい。

いくつかの実施形態において、パラメータのいくつかは、初期計測値に基づいて設定されてもよい。具体的には、基準トランスデューサではない、トランスデューサの空間的構成を決定することを目的として、ターゲット物体の運動の剛体モードと関連する周波数を維持しつつ、ターゲット物体の運動の撓みモードに関連する周波数を抑圧するように、センサ信号をフィルタリングすることが望ましい場合がある。

これを目的として、ターゲット物体は、例えば、モードハンマー、ハンマー様の物体を使用することにより、或いは、別の適切な方式により、励起することができる。計測システムは、例えば、モードハンマー又はハンマー様の物体の衝撃に応答して、などのように、励起に応答して、トランスデューサからセンサ信号を取得することができる。このステップは、例えば、それぞれの励起ポイントにおいて、複数の励起について、反復されてもよい。プロセスは、すべてのトランスデューサについて周波数スペクトル又は周波数応答関数を演算することができると共に、最低撓みモードの固有周波数及び最高剛体モードの固有周波数を推定することができる。次いで、プロセスは、トランスデューサの空間的構成を決定することを目的として、後続のセンサ信号に適用されるべき低域通過フィルタのカットオフ周波数を推定することができる。具体的には、カットオフ周波数は、低域通過フィルタが剛体モーションに関する情報を維持することになるが、ターゲット物体の撓みモーションを排除することになるように、選択することができる。更には、次いで、計測システムのデータ取得部分は、既定の閾値超であるサンプリング周波数によって記録を実施するように構成することができる。閾値は、例えば、低域通過フィルタのカットオフ周波数の２．５６倍などのように、低域通過フィルタの決定されたカットオフ周波数に応じて選択することができる。

ステップＳ４０において、プロセスは、トランスデューサからセンサ信号を取得している。これを目的として、ターゲット物体の運動、具体的には、剛体運動を誘発することができると共に、誘発された運動に応答するセンサ信号を取得することができる。

具体的には、データ取得が開始されてもよく、且つ、すべての６つの自由度（３つの平行運動及び３つの回転運動自由度）を経験するような方式により、ターゲット物体が振れ動くようにすることができる。これは、例えば、手で吊り下げられた物体を振り動かすことにより、或いは、非常にソフトな先端を有する重いハンマーを使用することにより、などのように、手動的に実行することができる。剛体モードの共振周波数は、通常、非常に低いことから、加速度信号内において十分な信号対ノイズ比を得るように剛体モードを励起することは、困難でありうるであろう。従って、例えば、６つの方向のそれぞれにおいて少なくとも２０～４０回の振れに基づいて、などのように、６つの自由度のそれぞれにおいて複数の振れに基づいて十分なデータを記録することが有利でありうる。取得されたセンサ信号は、時系列の計測値として表すことができる。いくつかの実施形態において、第１サブセットのトランスデューサの空間的関係のみが、まず、入力された場合には、このステップは、第１セットのトランスデューサの完全な空間的構成の演算によって後続されうる。

後続のステップＳ５０において、プロセスは、基準トランスデューサの既知の空間的構成から、且つ、取得されたセンサ信号から、基準トランスデューサではない、トランスデューサの空間的構成を演算している。トランスデューサの空間的構成を演算するための演算的プロセスの一例については、図３を参照して、以下において詳述することとする。以下の説明から明らかとなるように、演算の一部分として、プロセスは、望ましい精度によってすべてのトランスデューサの空間的構成を演算するべく、取得されたデータが、適している、或いは、十分である、程度を検出することができる。プロセスは、このような演算が可能ではない場合には、ユーザーに警告することができる。次いで、ユーザーは、例えば、信号対ノイズ比を改善するように、更に多くのデータを取得することにより、例えば、低域通過フィルタのカットオフ周波数などの、パラメータのいくつかを再調節することにより、或いは、（基準加速度計を定位置において維持しつつ）加速度計をその他の場所に移動させ、且つ、計測を反復することにより、データ取得を反復することができる。

すべてのトランスデューサの空間的構成が決定されたら、プロセスは、意図された構造解析のために、実際のデータ取得に進むことができる。これを目的として、ステップＳ６０において、更なるデータが取得されうる。具体的には、ターゲット物体は、撓み及び剛体モードの両方を励起するように運動させられてもよく、且つ、すべてのトランスデューサからの更なるセンサ信号が取得されてもよい。このステップにおいては、撓みモードを抑圧するための以前の低域通過フィルタリングを省略することができる。

ステップＳ７０において、プロセスは、取得された更なるデータに基づいて、且つ、すべてのトランスデューサの既知の且つ決定された空間的構成に基づいて、望ましい構造解析演算を実行することができる。一般的に、構造解析という用語は、力が、印加された際の、且つ、特に、動的に印加された際の、ターゲット物体の振る舞いに関する情報を結果的にもたらす、記録されたトランスデューサ信号に基づいた任意の演算方法を包含するものと解釈されたい。構造解析の例は、動的変位、時間履歴の分析、及びモード分析、振動分析などを含みうる。いくつかの実施形態において、構造解析は、例えば、モード形状及び周波数を算出するように、有限要素分析を含みうる。

図３は、本明細書において開示されている、トランスデューサ、特に、３軸加速度計、の空間的構成を決定するためのコンピュータ実装されたプロセスの一例のフロー図を示している。プロセスは、例えば、図１Ａとの関連において記述されている、計測システム１１０によって実行することができる。

初期ステップＳ５１において、プロセスが、入力データを受け取っている。具体的には、プロセスは、それぞれの基準加速度計の位置及び向きを受け取っている。例えば、位置は、ＬＣＳに関連する基準加速度計の座標として受け取られてもよい。それぞれの基準加速度計の向きは、ｚ－ｙ’－ｘ’’表記法に準拠した、３つのＴａｉｔ－Ｂｒｙａｎ角として受け取られてもよい。角度は、基準加速度計の、ＬＣＳからＭＣＳへの連続的な回転に対応している。位置及び向きは、手動的に入力されていてもよく、別の自動化された手順によって決定されていてもよく、或いは、基準トランスデューサの間の空間的関係の入力された情報に戻づいて、或いは、センサ信号に基づいて、演算されていてもよい。

プロセスは、例えば、図２のプロセスのステップＳ４０において取得された記録などの、ターゲット物体の誘発されたモーションに応答して記録された、すべての加速度計のセンサ信号を更に受け取っている。プロセスは、加速度計が計測の際に接続される、データ取得インターフェイスのチャネルテーブルと、基準加速度計に対応するチャネル、即ち、第１サブセットの加速度計に対応するチャネル及び第２サブセットの加速度計に対応するチャネル、に関する通知と、を更に受け取ることができる。

プロセスは、図２のプロセスのステップＳ３０との関連において記述されている、低域通過フィルタのカットオフ周波数などの、更なる情報を更に受け取ることができる。

後続のステップにおいて、プロセスは、第２サブセットの加速度計の空間的構成を演算している。それぞれのステップについて説明する前に、第２サブセットの加速度計の位置及び向きの計算の全体的な方法について、更に詳細に説明することとする。トランスデューサの第１サブセットが、組

からのＲ個の基準加速度計を有していると仮定し、我々は、これらを基準セットと呼称する。基準セットのすべての加速度計の位置及び向きは、既知である。組内のすべての加速度計用の式（１５）を積層することにより、以下が得られ、

この場合に、行列

は、トランスデューサ伝達率行列とも呼称することとする。基準加速度計用のトランスデューサ伝達率行列は、（１３）から判明する。それぞれのトランスデューサ伝達率行列は、対応するトランスデューサのセンサ信号を、ターゲット物体を表す剛体（即ち、加速度計）の（即ち、ターゲット物体と関連するローカル座標系ＬＣＳの）位置及び向きの第２導関数のベクトルに関連付けている。加速度計信号｛ａ_r（ｔ）｝は、基準加速度計について入手可能である。この結果、ベクトル

を計測された加速度の最小二乗平均として推定することが可能であり、

ここで、シンボル

は、疑似逆行列を表記している。

行列

は、基準加速度計のトランスデューサ伝達率行列から構築されており、これは、３Ｒ行と、６列と、を有する。式（１７）が解決可能となるには、行列

は、フルランクであることを要し、従って、これは、少なくとも６つの行を有することを要する。但し、２つの３軸基準加速度計（Ｒ＝２）は、フルランクの

行列を構築するには、十分でなく、且つ、少なくとも３つの基準加速度計（Ｒ＝３）が必要とされる（且つ、これらは、まっすぐなライン上において配置されることにならない）。

上述の条件が充足され、且つ、

が、（１７）を使用することにより、取得されうると仮定しよう。同一のボディに装着された、且つ、組

からのものと同期した状態において計測している、加速度計の別の組

を検討しよう。加速度計読取は、入手可能であるが、これらの位置及び向きは、未知である。それぞれの加速度

ごとに、その読取は、

に対する線形依存性を有しており、（１５）を参照されたい。

（１２）から、行列［β_q］は、

の関数であり、且つ、この３×６行列のそれぞれの要素について、分析（非線形）式が入手可能である、ということになる。

未知の

について解決することにより、組

からの任意の加速度計の位置及び向きを取得することができる。

この（非線形）問題を解決する別の方法が可能である。以下において、非線形問題を解決する方法の２つの例について開示する。

方式１
上述のように、（１８）内の行列［β_q］は、伝達率行列として見なすことが可能であり、これは、複数入力複数出力（ＭＩＭＯ：Ｍｕｌｔｉｐｌｅ Ｉｎｐｕｔ Ｍｕｌｔｉｐｌｅ Ｏｕｔｐｕｔ）モード分析から既知のエスティメータ、即ち、Ｈ₁、Ｈ_V、又はＨ_Sのうちの１つを使用することにより、計測された｛ａ_q（ｔ）｝及び推定された

から推定することができる（例えば、Ｊ． Ｓ． Ｂｅｎｄａｔ， Ａ． Ｇ． Ｐｉｅｒｓｏｌ， “Ｅｎｇｉｎｅｅｒｉｎｇ ａｐｐｌｉｃａｔｉｏｎｓ ｏｆ ｃｏｒｒｅｌａｔｉｏｎ ａｎｄ ｓｐｅｃｔｒａｌ ａｎａｌｙｓｉｓ”， Ｊｏｈｎ Ｗｉｌｅｙ ＆ Ｓｏｎｓ， Ｉｎｃ．， １９８０を参照されたい）。

においてサンプリングされた信号の場合に、加速度計読取及び

を行列として提示すれば、以下のとおりであり、

従って、（１７）は、以下のようになる。

最も単純なＨ₁エスティメータに従うことにより、以下が得られる。

次式を要素ごとに同等と見なすことにより、

６つの未知数について１８個の一次方程式を解くことができる。加速度計読取は、ノイズによって汚染されていることから、これは、［ε］の要素

或いは、その基準のいくつかを極小化する、行列方程式（２２）内の要素のベストフィットを生成する、このような

を見出す最適化問題として定式化することができる。

方式２
伝達率行列の１８個の要素をフィットする代わりに、計測された時間履歴をフィットすることができる。以下を検討してみれば、

ここで、

は、変数のベクトル

について演算された加速度計ｑの推定された読取である。計測された信号と推定値の間の不一致は、以下のとおりであり、

且つ、左手側の行列の行は、加速度計のｘ－、ｙ－、及びｚ－計測軸の予測誤差として見なすことができる。それぞれの計測軸のスカラーフィット誤差は、二乗平均として特徴付けすることが可能であり、

この場合に、式（２６）の左手側は、変数のベクトル｛ｖ｝の正のスカラー関数であり、これは、以下の正のスカラー関数に組み合わせることができる。

この結果、極小化問題は、次式として定式化することができる。

これは、例えば、加速度計の位置がターゲット物体の外側にはならない、且つ、オイラー角が通常は次式によって制約されている、などの、可能な制約を提示している。

上述の２つの方式のそれぞれは、スタンドアロン方式として、個別に利用することができる。或いは、この代わりに、両方の方式は、組合せにおいて利用することもできる。上述の方式１は、相対的に高速に収束しうるが、方式２は、相対的に正確な推定値を結果的にもたらしうる。従って、いくつかの実施形態において、方式１は、初期推定値を演算するべく使用されうる一方で、方式２は、方式１の出力を高度化するべく使用されうる。

一般に、上述の方法は、いくつかの仮定を実施している。具体的には、上述の方法は、ターゲット物体が剛体として運動している、誘発された運動の変位が小さい、且つ、ターゲット物体がすべての６つの自由度において運動している、と仮定している。

本明細書において開示されている方法の実施形態の精度は、上述の仮定が実際の計測状況において充足される程度に依存している。

以下、これについて更に詳細に説明することとする。

通常の計測状況において、ターゲット物体は、理想的な、制限されてはいない剛体として運動してはいない。但し、ターゲット物体が、例えば、自由－自由状態における物体の力学をモデル化するべく、ソフトバンジー及び／又は空気まくらによって支持されている際には、多くの計測シナリオが存在している。これは、物体のＦＥモデルの検証／チューニングのために実施される、実験的なモード分析（ＥＭＡ：Ｅｘｐｅｒｉｍｅｎｔａｌ Ｍｏｄａｌ Ａｎａｌｙｓｉｓ）用の通常のシナリオである。支持の剛性の選択は、剛体モードの周波数（物体の慣性プロパティ及び支持の剛性によって定義される）と撓みモード共振（物体の剛性に依存している）の最良の分離に基づいたものであってもよく、例えば、図４Ａ～図４Ｂに示されているように、これをチェックするべく、サンプルの周波数応答関数（ＦＲＦ：Ｆｒｅｑｕｅｎｃｙ Ｒｅｓｐｏｎｓｅ Ｆｕｎｃｔｉｏｎ）又は加速度計の自己スペクトルを使用することができる。図４Ａは、０～９００Ｈｚの周波数範囲におけるターゲット物体のＦＲＦを示している一方で、図４Ｂは、０～５０Ｈｚの周波数範囲におけるＦＲＦの拡大図を示している。図４Ａ～図４Ｂからわかるように、最低撓みモードは、約４６２ＨＺにおいて発生している一方で、最高剛体モードは、約６．５ＨＺにおいて発生している。従って、撓みモードと剛性モードの間の良好な分離が実際に可能である。従って、低域通過フィルタを適用し、且つ、最高剛体モード超において、且つ、最低撓みモード未満において、そのカットオフ周波数を設定することにより、応答の撓み成分を大幅に減衰させることが可能であり、且つ、剛体モーション用の仮定を充足することができる。

すべての６つの自由度におけるターゲット物体の運動に関連し、（１８）を解決する、すべての６つの未知数を見出すには、行列

は、フルランクであることを要する（即ち、そのランクは、６に等しいことを要する）。式（１７）は、（１９）と同様に、サンプリングされた加速度の信号について、以下のように書き換えられる。

Ｓｙｌｖｅｓｔｅｒの不等式によれば（Ｍａｔｒｉｘ Ｃｏｏｋｂｏｏｋを参照されたい）、以下のとおりである

のランクについては、既に説明済みであり、適切な方式によって基準加速度計を配置することにより、この行列のフルランクが容易に実現可能である（これは、まっすぐなラインに沿って取り付けられてはいない、３つ以上の基準加速度計であることを要することを思い出して頂きたい）。そして、

がフルランクである場合には、その疑似逆行列も、フルランクであり、

である。

従って、

は、

のフルランク用の必要な条件である。これを実現するべく、ターゲット物体は、すべての６つの自由度においてモーションを経験することを要し、且つ、データ取得の際にすべての６つの方向においてモーションを経験するような方式により、特別な注意が、ターゲット物体の励起に対して払われることを要する。

低域通過フィルタを適用した後に、

の特異値をチェックすることにより、この仮定が充足される方式をチェックすることが可能であり、これについては、後述することとする。

式（１１）を得るべく、我々は、平行運動及び回転運動変位が小さいものと仮定している。実際には、この仮定は、回転変位θが、

、

を充足するべく十分に小さい際に、充足され、これは、

である場合に、最大で小数点以下２桁まで、充足される、ことが見出された。同一の大きさを得るべく、平行運動変位ｄは、

である範囲内にあることを要し、この場合に、ｌは、ターゲット物体のなんらかの特徴的な長さである。

その一方で、加速度計信号の良好な信号対ノイズ比を保証する、いくつかの妥当な変位が励起されることを要し、これは、非常に低い周波数においては、実現が困難でありうることに留意されたい。

本発明者らは、ターゲット物体がすべての６つの方向において手で運動させられた際にターゲット物体を励起することで、或いは、これをソフトな先端を有する相対的に重いハンマーによって励起することで十分であることを見出した。ハンマーを使用する際には、ランダム方向における多くのポイントにおいて物体を打撃することが有用であることが見出された。この結果、構造は、その初期位置の周りにおいて発振を経験することになり、これは、この仮定を充足するための非常に有用な方式であると考えられる。

再度図３を参照すれば、ステップＳ５２において、プロセスは、決定されたカットオフ周波数を使用することにより、受け取られたセンサ信号の低域通過フィルタリングを実行している。

ステップＳ５３において、プロセスは、ターゲット物体の誘発された運動が、剛体モーションとして表現されうるかどうか、並びに、モーションがすべての６つの自由度において発生しているかどうか、を検証している。これを目的として、プロセスは、すべての加速度計信号との関係において［ａ（ｔ_i）］の特異値を演算することができると共に、特異値＃６及び＃７の間の差及び／又は比率が既定の閾値超であるかどうかを決定することができる。図５は、ターゲット物体の例示用の運動について決定された正規化済みの特異値の例を示している。特異値＃６及び＃７の間のギャップは、矢印５５１によって示されている。

誘発された運動が上述の基準を充足している場合には、プロセスは、Ｓ５４に進む。さもなければ、プロセスは、例えば、運動を誘発するステップを反復する（例えば、図２のプロセスのステップＳ４０に戻る）ようにユーザーに指示するなどの、適切なエラーメッセージと共に、戻っている。

ステップＳ５４において、プロセスは、（例えば、特異値分解と、これに続く切り捨てを実行することにより）６個超の特異値に対応するデータを除去し、基準の間の最大距離が１となるように座標をスケーリングし、且つ、式（１６）を使用することにより、行列

を演算している。

ステップＳ５５において、プロセスは、演算された行列

の状態数を既定の閾値と比較している。状態数が閾値超である場合には、プロセスは、基準加速度計の異なる配置、即ち、まっすぐなラインから相対的に大きく変位した配置、を選択するように（例えば、図２のプロセスのステップＳ２０に戻るように）、且つ、データ取得を反復するように、ユーザーに指示する適切なメッセージと共に、戻っている。さもなければ、プロセスは、ステップＳ５６に進む。

ステップＳ５６において、プロセスは、式（１７）を利用することにより、ターゲット物体のローカル座標系の加速度を演算している、即ち、プロセスは、

を演算している。上述のように、このステップは、一次方程式の系の最小二乗解を演算するステップを含みうる。

ステップＳ５７において、プロセスは、第２サブセットの加速度計のそれぞれのものの空間的構成、即ち、それぞれの非基準加速度計の位置及び向き、を演算している。第２サブセットの加速度計の空間的構成を演算する方法の一例については、図７を参照して更に詳細に後述することとする。

すべての基準加速度計の空間的構成が演算された際に、プロセスは、ステップＳ５８に進み、ここで、すべての座標が、再スケーリングにより、オリジナルのスケーリングに戻されている。

ステップＳ５９において、プロセスは、後続の計測の際に使用される、決定済みの空間的構成を保存している。又、プロセスは、例えば、図６Ａ～図６Ｂにおいて示されているように、例えば、グラフィカルユーザーインターフェイスを介して、結果をユーザーに提示することもできる。図６Ａは、それぞれの加速度計の決定された空間的構成を示すターゲット物体の表示された３Ｄデジタルモデルを概略的に示している。図６Ｂは、表示された３Ｄデジタルモデルの一部分の拡大図を示している。

図７は、本明細書において開示されている、第２サブセットのトランスデューサの、具体的には、３軸加速度計の、空間的構成を決定するコンピュータ実装されたプロセスの一例のフロー図を示している。プロセスは、例えば、図１Ａとの関連において記述されている計測システム１１０によって実行することができる。

初期ステップＳ５７１において、プロセスは、式（２１）を使用することにより、伝達率行列

を演算している。

ステップＳ５７２において、プロセスは、例えば、基準加速度計の位置の平均として、開始ポイントを演算している。

ステップＳ５７３において、プロセスは、例えば、ニュートンの方法を使用することにより、未知の位置及び向きの初期推定を演算している。

ニュートンの方法が収束した場合には、プロセスは、ステップＳ５７５に進み、そうではなくて、方法が収束しない場合には、プロセスは、ステップＳ５７４に進み、ここで、プロセスは、ランダムな変位を開始点に追加し、且つ、ステップＳ５７３に戻っている。アルゴリズムが、既定数の試みの後に（例えば、５～１０回の試みの後に）収束しない場合には、プロセスは、誤りによって現時点のアクセレレータにマーキングし、且つ、次の加速度計に進む。誤りによってマーキングされた加速度計は、再位置決めされなければならない場合もあり、或いは、それらの位置が、手動で決定されなければならない場合もある。

ステップＳ５７５において、プロセスは、現時点の加速度計の演算済みの位置及び向きを保存している。

任意選択により、プロセスが更なる極小化ステップＳ５７６を実行する際には、結果の精度が更に改善されうることが見出された。具体的には、追加の極小化ステップは、開始ポイントとして、ニュートンの方法の結果を使用してもよく、例えば、式（２９）を使用することにより、制約の組を生成してもよく、且つ、例えば、式（２８）において定義されている適切な目的関数を極小化する極小化アルゴリズムを実行してもよい。最適化問題（２８）は、６つの変数

の強力な線形関数であると共に、後の３つの変数において周期的である、関数

を極小化するべく試みている。極小化は、線形制約によって拘束されている。極小化は、例えば、ＭＡＴＬＡＢ（登録商標）のＦＭＩＮＣＯＮの使用などの、任意の適切なルーチンにより、実行することができる。このルーチンは、関数及び制約のタイプに応じて、最適化方法を自動的に選択している。この特定の最適化問題の場合には、Ｊａｃｏｂｉａｎ及びＨｅｓｓｉａｎ（変数との関係における第１及び第２導関数）を提供するニーズは、存在していない。但し、感度を必要としている別の最適化アルゴリズムを使用するケースにおいては、これらは、容易に取得されうる。目的関数の値が、極小化の後に、相対的に小さくなった場合には、プロセスは、最適化結果を位置及び向きの最終的な推定値として保存することができる。プロセスは、例えば、フィットの相対誤差などの、結果の品質インジケータを更に保存することができる。

上述の方法に対する多数の変形が可能であることを理解されたい。以下においては、トランスデューサの空間的構成を決定するプロセスの一代替実施形態について説明することとする。上述の実施形態は、基準加速度計の位置及び向きの知識に依存しており、少なくとも３つ加速度計が、既知の座標を有する位置において取り付けられなければならず、且つ、例えば、ＧＣＳに対して、慎重にアライメントされなければならない。基準加速度計の場所／向きにおける精度が、その他の加速度計の結果的に得られる位置／向きに直接的に影響を及ぼしている。以下において提示されているプロセスの実施形態においては、基準加速度計の位置及び向きを決定する要件が大幅に緩和されており、これにより、その間の距離により、正確な基準加速度計の位置及び向きのニーズが置換されている。従って、この実施形態においては、図２のプロセスのステップＳ２０において、ターゲット物体上において取り付けられた基準加速度計の間の相対距離のみが決定される。この実施形態は、基準なし環境とも呼称することとする。その他の基準なし実施形態に対する入力としては、相対距離の代わりに、トランスデューサの間のその他の空間的関係が使用されうることを理解されたい。

３つの基準加速度計の間の距離が提供された状態において、方法は、これらのうちの１つに装着されたＬＣＳ内におけるその位置及び向きを出力している。方法の出力は、上述の方法に対する入力として見なすことが可能であり、従って、同一のデータセットに基づいて、すべてのその他の加速度計の位置及び向きを取得することができる。

これを目的として、３Ｄ空間内において制限されていない状態において運動しうる変形不能なターゲット物体上に取り付けられた３つの３軸加速度計を検討しよう。ターゲット物体が、すべての６つの自由度において小さな変位を経験するものと仮定しよう。更には、加速度計を文字

、

、及び

によってマーキングしよう。以前と同一の表記法を使用することにより、加速度計のＬＣＳをＭＣＳに回転させる、ＬＣＳ内におけるその座標及び（ｚ－ｙ’－ｘ’’表記法を使用する）３つのＴａｉｔ－Ｂｒｙａｎ角を使用して、それぞれの加速度計の位置及び向きを定義しよう。利便を目的として、ベクトル

を使用しよう。ここで、

は、

、

、又は

のいずれかである。次に、加速度計

の初期位置をＬＣＳの原点として選択し、且つ、その初期向きにおける加速度計

のＭＣＳをＬＣＳの向きとして選択しよう。この結果、

である。

加速度計の読取が入手可能であって、

であり、且つ、我々は、この９つの信号の組を２つの組に分割することが可能であり、

この場合に、

であり、且つ、

であり、この説明を目的として、我々は、第１の組を基準セットと呼称し、且つ、第２の組を試験セットと呼称することとする。

（２９）によれば、両方の組内の加速度を

ベクトルから推定することができる。

第２の式から

を隔離し、且つ、第１のものに代入することにより、以下が得られる。

行列［β_r］及び行列［β_q］の要素は、加速度計（１３）の位置及び向きに依存しており、且つ、

及び

として記述することができる。任意の

及び

の場合に、試験セット信号は、次式として、基準セット信号から推定することが可能であり、

ここで、行列

は、基準セットからの信号と試験セットからの信号の間の伝達率行列として解釈することができる。

試験セットからの予測されたｊ番目の信号と対応する計測された信号の間の差は、以下のとおりである。

サンプリングされた形態において、

であり、この場合に、時間サンプルは、

である。この結果、差は、その二乗平均により、特徴付けることが可能であり、

且つ、すべての３つの試験チャネルについて、以下のとおりである。

いまや、極小化問題を次式として定式化することができる。

（４０）における制約は、加速度計位置がターゲット物体寸法の外側において求められることを許容していない幾何学的制約と、オイラー角に関する制約（２９）と、である。

本発明者らは、積

は、線形スケールに対して不変である、即ち、

及び

におけるすべての線形寸法が同一のスカラーによって乗算された場合に、結果的に得られる行列積は、同一である、と考えている。直観的に、線形距離を計測する際には、例えば、インチ又はミリメートルなどの異なる線形単位を使用するものとして理解されうるであろう。これは、解が最大で乗算子までの距離を含むことになることを意味しており、極小化を使用するケースにおいては、極小化ルーチンは、収束において困難を有することになる。これを回避するための１つの方法は、例えば、最適化ルーチンに３つの加速度計の間の距離

、

、及び

を提供するなどの、余分な制約を導入する、というものである。これらの制約は、以下のとおりとなり、

ここで、

は、

及び

について演算された加速度計

及び

の間の推定距離であり、且つ、ε_lは、許容された線形許容範囲である。

基準セットは、すべての３つの加速度計からの６つの信号の組合せから構成されていることに留意されたい。便利な方法（但し、必須ではない）は、これを加速度計

のすべての３つの信号、

からの１つ又は２つの信号、並びに、

からの残りのもの、から構成するというものであり、例えば、可能な基準セットは、

であってよく、この場合に、添え字は、所与の加速度計の計測方向を示している。この結果、試験セットは、残りの信号、即ち、

を含むことを要する。最適な信号を選択する方式については、後述することとする。

［β_r］及び［β_q］行列を以下のように構築することが便利であり、まず、以下の完全な９×６行列を構築し、且つ、

次いで、基準及び試験セットに対応する行を抽出する。

同一の内容は、時間履歴についても当て嵌まる。

が、３つの加速度計からのすべての時間履歴の集合体である場合には、以下のとおりである

式（３５）は、［β_r］の逆元を必要としていることから、｛ａ_r（ｔ）｝内の６つの時間履歴が可能な限り独立していることが好ましい場合がある。サンプリングされた形態

において記述された場合に、行列｛ａ_r（ｔ_i）｝の状態値が最小になるような方式により、組Ｒの内容を選択することが好ましい場合がある。Ｒ内の信号の可能な組合せの数が

であり、且つ、ＳＶＤ動作は、非常に高速である、ことから、すべての可能な組合せを反復し、且つ、最小の状態数を付与するものを選び出すことが可能である。

要すれば、（４１）によって拡張された制約の組を有する極小化問題（４０）を解決することにより、加速度計

のＭＣＳによって定義されたＬＣＳ内における加速度計

及び

の位置及び向き

及び

が得られることになる。

（４０）における目的関数

は、最適化問題が非線形制約（４１）に拘束された状態において、その１２個のパラメータの強力な非線形関数である。いくつかの状況において、最適化問題の収束が、良好な開始ポイントの選択の影響を受けやすい場合がある。

この問題を解決するには、最適化問題の緩和された定式化を探すことが有益でありうる。このような緩和された定式化は、依然として、実際的であることを要することに留意することが重要である。

以下、我々は、変数の数が１２から８に低減される、可能な緩和された定式化の１つについて示唆している。これは、そのＸ０Ｙプレーンが水平方向であるような方式により、ＬＣＳを定義すことにより、実現することができる。ＬＣＳは、加速度計

のＭＣＳによって定義されていることから、実際には、これは、そのＺ軸が垂直方向であることを要することを意味している。同一の内容は、加速度計

及び

についても必要とされることになろう。この要件が充足された場合に、後の２つの加速度計のＬＣＳからＭＣＳへの回転は、１つの回転、即ち、角度

及び

による、且つ、その他の４つの角度を０に設定する、即ち、

である、ＬＣＳのＺ軸を中心とした回転、のみから構成されている。

この要件の実際的な実装形態は、充足が特に困難ではない。スピリットレベルを使用することにより、そのＸ０Ｙプレーン（又は、ＭＣＳの別のプレーン）が水平方向となるように、３つの加速度計をアライメントすることは、相対的に容易である。Ｘ０Ｙプレーン内の加速度計の向き（即ち、その向き、

及び

角度）に関する特別な要件は、存在していないが、実際には、極小化ルーチンの開始ポイントとしてこれらの値を使用するべく、加速度計の概略的な向きを知ることが有利であることに留意されたい。

実際には、緩和されていない最適化が、後述するように、緩和されたものに後続することができる。
１．３つの基準加速度計が、緩和された定式化に従って取り付けられる。
２．開始ポイントとして、基準加速度計の概略位置と、その向き

及び

と、を使用することにより、緩和された最適化が実行される。
３．緩和された定式化の収束の後に、結果的に得られた値は、高度化された位置及び向きを取得するべく、緩和されていない最適化用の開始ポイントとして使用することができる。
４．最後に、同一のデータセットに基づいて、すべてのその他の加速度計の座標及び向きを取得するべく、基準加速度計の座標及び向きが使用される。

このようにして、例えば、加速度計が重力の方向と完全にアライメントされていない状況において、結果の精度を向上させることができる。

最適化問題（４０）は、１２個の変数（緩和されたケースの場合には、８つの変数）の強力な線形関数であると共に、そのいくつかにおいて周期的である、実際的な正のスカラー目的関数

の極小値を求めている。この極小化は、非線形制約の組（４１）によって拘束されている。

ＭＡＴＬＡＢの最適化ツールボックスの一部分である、ＭＡＴＬＡＢのＦＭＩＮＣＯＮは、極小化問題を解決するための便利なルーチンである。この特定の定式化において、ＦＭＩＮＣＯＮは、変数との関係において、目的関数及び制約の感度を提供するステップを必要としている。これは、目的関数及び制約の１次導関数（Ｊａｃｏｂｉａｎ）及び２次感度（Ｈｅｓｓｉａｎ）の両方を必要としている。ＦＭＩＮＣＯＮは、有限差分を使用することにより、導関数を自動的に演算することができる。但し、可能な場合には常に、（正確な）分析感度を提供することが好ましくありうる。感度の導出については、後述する。

目的関数（３９）は、緩和された又は緩和されていない定式化が検討されるかどうかに応じて、８個又は１２個の変数の関数である。すべての変数が単一のベクトル

として組み合わせられる場合には、目的関数のＪａｃｏｂｉａｎは、そのｉ番目の要素が以下である、ベクトル

である。

目的関数（３９）の定義から、以下のとおりである。

（３５）～（３８）を利用することにより、以下が得られ、

この場合に、予測された且つ計測された信号の間の差

は、サンプリングされた信号の式

及び

を使用することにより、以下のとおりである。

｛｝又は［］括弧の後の添え字ｊは、それぞれ、ベクトルのｊ番目の要素又は行列のｊ番目の行を表記していることに留意されたい。ここで、予測された且つ計測された信号の間の差は、

である。

以下を微分し、

（４８）を代入し、

且つ、（３６）を利用することにより、以下が得られる。

行列の導関数が、かなり簡単である一方で、逆元の導関数は、以下をより必要としている。

上述のすべてのものを組み合わせ、且つ、（４５）に代入することにより、Ｊａｃｏｂｉａｎの分析式を得ることができる（数値方法を使用して演算することが相対的に便利である逆行列を除く）。

Ｈｅｓｓｉａｎは、その要素が以下のとおりである、正方行列である。

Ｈｅｓｓｉａｎの要素は、Ｊａｃｏｂｉａｎにおけると同一の考慮事項を使用することにより、演算することができる。

制約の感度
最適化問題の標準的な制約式を踏襲することにより、制約（４１）は、以下のように再定式化することが可能であり、

この場合に、所与の可変値

の二乗距離の推定値は、以下のとおりであり、

且つ、制約（５４）の勾配は、非常に単純である。制約勾配は、３×８又は３×１２行列であり、一例として、わずかな数の要素が以下において提供されている。

オイラー角との関係における導関数は、ゼロである。

上述の方法の実施形態は、ターゲット物体が剛性であり、且つ、そのモーションが制約されていない、と仮定している。以下においては、そのモーションが制約されうる、曲がりやすいターゲット物体に適用されうる、一実施形態について説明することとする。多くの用途において、このような物体は、提示された方法が有用でありうる、ターゲット物体の大きな部分を表している。例は、試験リグ上において取り付けられた風力タービンブレード、橋、タワー、煙突のような土木構造、及びこれらに類似したものを含む。

図８は、そのモーションが制約されている、変形可能なターゲット物体の一例を示しており、これは、いくつかの境界条件によって記述されている。ターゲット物体８１０の変形されていない状態において、ターゲット物体のそれぞれのポイントＰは、ＧＣＳにおいては｛ｒ｝_G∈Ωである、半径ベクトルによって記述された座標を有しており、この場合に、Ωは、ターゲット物体が定義されているドメインである。任意の時点ｔにおいて、ポイントＰの変位は、ベクトル

によって定義されている。この結果、その加速度は、以下のとおりである。

ここでは、且つ、更には、本説明においては、ＧＣＳを表記する添え字が省略されている。

変形されていない状態８１０におけるターゲット物体を検討しよう。３軸加速度計８１１は、ＣＧＳ内の座標｛ｒ_n｝∈Ωを有するポイントＰ_nにおいて物体上において取り付けられており、且つ、その向きは、そのＭＣＳが回転行列

によってＧＣＳから得られうるようになっている。

ターゲット物体が振動した際に、加速度計によって検知される３つの（スカラー）信号は、以下のとおりであり、

ここで、｛ｉ_n（ｔ）｝、｛ｊ_n（ｔ）｝、及び｛ｋ_n（ｔ）｝は、ＧＣＳ内において提示されたＭＣＳの場所である。ＧＣＳ内における場所座標は、以下のように取得することが可能であり、

ここで、次式

は、ＭＣＳからＧＣＳへの回転行列である。（５９）、（６０）を利用し、且つ、（５８）内のスカラーをベクトル

として組み合わせることにより、以下が得られる。

Ｇａｌｅｒｋｉｎ方法を使用することにより、以下を拡張することが可能であり、

この場合に、

は、境界条件（所謂、基本境界条件）を充足する、且つ、互いに直交する、時間独立ベクトルの組である。Ｄ_k（ｔ）は、時間依存スカラーである。

の便利な（但し、必須ではない）選択肢の１つは、物体のモード形状であり、その結果、Ｄ_k（ｔ）は、所謂モード座標である。

任意の所与の制限された周波数範囲において、（６２）は、次式のとして提示することが可能であり、

この場合に、第１合計は、ターゲット物体の全体的な変位に大きく寄与しており、且つ、後者は、そうではなく、従って、無視されうる。時間との関係において２回にわたって微分することにより、以下が得られ、

ここで、

である

（６４）を（６１）に代入ことにより、以下が得られる。

時間依存回転行列

は、変形されていない状態における回転行列

と回転行列

の組合せとして提示することが可能であり、この場合に、添え字

は、物体の変形した状態におけるＭＣＳを表記している。

これは、変形されていない状態におけるＭＣＳから変形した状態におけるＭＣＳへの回転のシーケンスを定義する、３つのオイラー角θ₁（ｔ）、θ₂（ｔ）、及びθ₃（ｔ）に依存している。

振動の大きさが、小さいものと仮定し、且つ、式中内の小さな項をマーキングするためのブックキーパーとしてεを利用し、

且つ、

、

、及び

であるものと仮定しよう。回転行列の構築は、（８）に類似しており、小さなオイラー角

及び

の場合に、行列は、以下のようになる。

（６８）を（６６）に代入し、且つ、次いで、（６７）と共に、（６５）に代入することにより、以下が得られる

レベルε²以上の項を省略し、且つ、ブックキーパーを除去することにより、以下が得られる。

この合計は、以下の積として書き換えることが可能であり、

ここで、ｋは、組

におけるインデックスであり、且つ、Ｋは、そのサイズである。従って、以下のとおりである。

いまや、組

からのＲ個の加速度計の位置及び向きは、既知であり、我々は、これらを基準セットと呼称するものと仮定しよう。従って、

用の

及び

は、既知である。基準セットからのすべての加速度計用の式を１つに積層することにより、以下が得られる。

又、これらの加速度計の読取（時間履歴）が任意の時点ｔにおいて入手可能であるものと仮定しよう。従って、（７３）の右手側における時間依存係数は、次式として推定することができる。

（７４）の解決可能性のための必要な条件は、（７３）内の行列が、正方又は垂直方向である、即ち、その行の数が列の数以上である、即ち、

である、というものである。

記述されている、係数Ａ_kを取得するための手順は、例えば、（Ｄ． Ｊ． Ｅｗｉｎｓ， “Ｍｏｄａｌ ｔｅｓｔｉｎｇ： Ｔｈｅｏｒｙ， ｐｒａｃｔｉｃｅ ａｎｄ ａｐｐｌｉｃａｔｉｏｎｓ”， Ｒｅｓｅａｒｃｈ Ｓｔｕｄｉｅｓ Ｐｒｅｓｓ Ｌｔｄ．， ２０００）において詳細に付与されている、モード座標を算出する手順に似ている。このケースにおいては、ベクトル

は、ターゲット物体のモード形状である。

いまや、演算された係数Ａ_k（ｔ）から、基準セット内にはない加速度計（即ち、第２サブセットの加速度計）の位置及び向きを以下のように決定することができる。

利便のために、以下においては、書き換えられる、式（７２）は、いくつかの周波数範囲において、ポイント｛ｒ_n｝において配置された、且つ、その回転行列が

となるように方向付けされた、ｎ番目の加速度計の読取が、次式のように推定されうるものと主張している。

ここで、我々は、加速度計読取｛ａ_n（ｔ）｝（３つの時間履歴）が入手可能であり、ベクトル

が、任意のポイント｛ｒ｝∈Ωにおいて、既知であり（或いは、近似されうる）、且つ、時間依存重みＡ_k（ｔ）が、以前の節において記述されているように、見出されうる、ものと見なしている。

においてサンプリングされた信号の場合に、加速度読取｛ａ_n（ｔ）｝及び

は、以下の行列として提示することが可能であり、

従って、（７５）は、次式となる。

式（７７）は、（２０）に似ており、これは、制約されていない剛体のケースについて定式化されたものであり、且つ、２つのケースの間の類似性を示している。

回転行列は、オイラー角

の関数であることから、３つのオイラー角

、

、及び

によって記述された、加速度計座標｛ｒ_n｝及びその向き、という６つの未知数が存在している。すべての６つの未知数を単一のベクトル

として組み合わせることが便利である。

ベクトル｛ｖ_q｝によって特徴付けられたトライアル位置／向きにおいて、推定された加速度は、以下のとおりであり、

且つ、それぞれの時点サンプルにおいて、計測された且つ推定された加速度の間の誤差は、以下のとおりである。

それぞれの計測軸の誤りは、二乗平均により、特徴付けすることが可能であり、

この場合に、式（８０）の左手側は、変数のベクトル｛ｖ_q｝の正のスカラー関数であり、これは、以下のように、スカラー関数として組み合わせることができる。

この結果、極小化問題は、以下のように定式化される。

これは、例えば、加速度計の位置がターゲット物体の外側にはならない

などの、可能な制約を提示しており、且つ、オイラー角は、通常、以下によって拘束されている。

以上において提示されている実施形態は、加速度計の位置／向きを見出す方法が、そのモーションがいくつかの（十分に剛性の）支持部によって制約されている、変形可能なターゲット物体に適用されうる方式を実証している。このようなターゲット物体の例は、「原位置」において試験されている多くの物体を含み、これらは、通常、例えば、風力タービン、タワー、橋などのような、大きな土木構造である。

方法は、（ｉ）ターゲット物体の基本境界条件を充足する、且つ、（ｉｉ）相互に直交する、空間的関数の組

を見出しうるものと仮定している。このような関数の最良の組は、例えば、ターゲット物体のＦＥ分析から判明する、設定されたモード形状である。上述の条件を充足する任意の別の組が有効であるが、これらは、相対的に多くの形状の検討を必要としうると共に、その結果、相対的に多くの基準加速度計を必要としうる。

方法は、
－時間依存係数

を演算するステップであって、この手順は、例えば、モード分解においてモード座標を取得するステップに類似している、ステップと、
－フィット誤りを極小化することにより、加速度計の位置及び向きを算出するステップと、
を有する。

例
図９は、その上に取り付けられたトランスデューサを有するターゲット物体の一例の異なる図を示している。トランスミッションハウジング９０１は、図９において示されているように、ゴムバンド９０２によって支持されている。Ｂｒｕｅｌ ＆ Ｋｊｅｒ Ｓｏｕｎｄ ＆ Ｖｉｂｒａｔｉｏｎ Ｍｅａｓｕｒｅｍｅｎｔ Ａ／Ｓ， Ｎｅｒｕｍ， Ｄｅｎｍａｒｋから入手可能である４５０６Ｂ型の３つの加速度計９０３が、基準加速度計として選択され、且つ、図９において示されているように、取り付けられている。いずれも、Ｂｒｕｅｌ ＆ Ｋｊｅｒ Ｓｏｕｎｄ ＆ Ｖｉｂｒａｔｉｏｎ Ｍｅａｓｕｒｅｍｅｎｔ Ａ／Ｓ， Ｎｅｒｕｍ， Ｄｅｎｍａｒｋから入手可能である、スピリットレベルＵＡ－１４０及びスイベルベースＵＡ－１４７３を使用することにより、加速度計は、そのＹ軸が、垂直であり、且つ、上方を指し示すように、アライメントされている。基準加速度計の座標は、Ｂｏｕｌｄｅｒ Ｉｎｎｏｖａｔｉｏｎ Ｇｒｏｕｐ， Ｉｎｃ．， Ｂｏｕｌｄｅｒ， ＣＯ， ＵＳＡから入手可能である、ＢＩＧ ＦＰ５５００ ３Ｄ Ｃｒｅａｔｏｒを使用することにより、サブミリメートルの精度により、計測された。
・基準

（赤色）： ２２７．７； ２４８．５； ６６５．８ ｍｍ
・基準

（青色）： １７９．２； ３８０．６； ２２．９ ｍｍ
・基準

（緑色）： ３７６．９； ８６．４； ５９．８ ｍｍ

座標は、ＧＣＳにおいて付与されており、ＧＣＳのＺ軸が上方を指し示している状態において、トランスミッションが垂直方向において直立している際に計測されている。それぞれのセンサの８つのコーナーがデジタル化されており、且つ、提供されている座標は、８つのコーナーの幾何学的中心に対応している。

構造の主要な特徴（大部分の別個のエッジ及びリブ）が、非常に正確なワイヤフレームモデルを提供するＡＴＣソフトウェアを使用することにより、ＢＩＧ ＦＰ５５００ ３Ｄ Ｃｒｅａｔｏｒにより、ポイントごとに、デジタル化されている。又、視覚化を改善するべく、トランスミッションハウジングは、Ｓｏｎｙ Ｘｐｅｒｉａ ＸＺ１携帯電話機上において稼働するＳｏｎｙ ３Ｄ Ｃｒｅａｔｏｒアプリケーションを使用することにより、デジタル化されており、これは、Ｗａｖｅｆｒｏｎｔ ＯＢＪファイルを介してＭＡＴＬＡＢにエクスポートされている。この表現は、ノード及び三角面から構成されている。

更に詳しくは、以下のプロセスが利用されている。
１．基準加速度計が、上述のように、取り付けられ、且つ、アライメントされた。
２．構造が、ＢＩＧ ＦＰ５５００を使用することにより、デジタル化された。
ａ．（そのコーナーに基づいた）基準加速度の中心の座標が取得された。
ｂ．十字形によって注記された、且つ、付番された、３１個のポイントが、構造上において選択され、且つ、ＢＩＧ ＦＰ５５００によって取得された十字形の座標が表記された。
ｃ．構造に対する加速度計の視覚的なアライメントを許容するべく、３１個のポイントの１２個が、向き精度をチェックするべく、選択された、即ち、（コーナーのような）構造の特徴の隣に配置された。
３．６つのセットアップ用の計測が実行された（セットアップ［１、２、４、５、６、７］と命名された）。
ａ．それぞれのセットアップは、３つの基準加速度計と
ｂ．５つの更なる加速度計９０４と、を含み、これらは、既知の場所においてトランスミッションハウジングに装着された。
ｃ．セットアップ＃＃４、５、６は、向きの精度を計測するべく、特別に設計された。
４．それぞれのセットアップごと、
ａ．すべての８つの加速度計の位置及び向きが、本明細書において記述されている基準なし実施形態を使用することにより、演算された。
ｂ．加速度計の取得された位置が、その既知の座標と比較された。これは、表面に垂直である、表記されている計測軸を使用することにより、加速度計の中心を構造の表面に「投射」することにより、実行された。
ｃ．セットアップ＃＃４、５、６の場合に、ＭＣＳの場所が取得され、且つ、場所の間の角度が演算された。これは、場所のスキャッタに関する統計の生成を許容し、これにより、向きの誤差について結論付けることを許容している。
５．３５個の推定値に基づいて、平均位置決め誤差は、１１．２ｍｍである（σ＝４．８）。

図１０Ａは、位置決め誤差のヒストグラムを示している。このヒストグラムは、非基準加速度計に基づいている。又、基準加速度計は、試験加速度計に属するものと見なされうることから、これらの位置も、既知の位置と比較されうるが、これらは、格段に小さな誤差を結果的にもたらす。従って、これらは、統計には含めなかった。図６Ａ～図６Ｂには、センサの位置特定の一例が示されている。具体的には、図６Ａは、トランスミッションハウジング９０１のデジタル化された３Ｄモデルの図を示しており、この場合に、加速度計の決定された位置及び向きが、例示用のセットアップ＃２について示されている。トランスミッションハウジングの拡大された一部分である、図６Ｂには、黒色ドットによって示された加速度計の既知の場所を含む、２つの加速度計位置が示されている。位置特定の誤りは、黒色ドット（３Ｄ Ｃｒｅａｔｏｒによって取得されたもの）と、赤色矢印によって図６Ｂにおいて示されている、取付表面に対する推定された加速度計の中心の投射である、第２のドットの間の距離として、推定されている。

図１０Ｂは、向きの誤差のヒストグラムを示している。それぞれのセットアップごとに、５つの試験加速度計、従って、５つの場所／方向、が存在していた。場所のそれぞれのペアごとに、その間の角度が（場所のスカラー積のアークコサインとして）演算され、これは、１０個の角度／方向、即ち、３つの方向について３０個の角度、を結果的にもたらした。３つのセットアップ（＃４、５、６）の場合に、統計は、以下のとおりである。
ｉ．セットアップ＃４： 平均 ４．３⁰， σ ＝２．３⁰
ｉｉ．セットアップ＃５： 平均 ４．４⁰， σ ＝２．０⁰
ｉｉｉ．セットアップ＃６： 平均 ６．０⁰， σ ＝２．５⁰

図１０Ｂには、セットアップ＃４の場合のヒストグラムが示されている。

従って、以上に基づいて、本明細書において開示されている方法は、現時点の技法によって実現されているものに匹敵する、且つ、通常の用途には十分である、妥当な品質の推定値を提供している、と結論付けることができる。

本明細書において記述されている方法の実施形態は、試験セットアップの際に、特に、ターゲット物体上のトランスデューサの位置及び向きの決定のために、利用することができる。

又、本明細書において記述されている方法の実施形態は、例えば、加速度計の位置／向きが試験の際に変化しなかったことを検証するため、などのように、加速度計の取付を検証するべく、利用することもできる。

又、同様に、本明細書において記述されている方法の実施形態は、ターゲット物体の単純化された形状の、例えば、デジタル３Ｄモデルなどの、デジタル表現を生成するべく、利用することもできる。いくつかの計測シナリオにおいては、例えば、取得された周波数応答関数又はモード形状を視覚化するために、ターゲット物体の単純化された形状を生成することが便利である。示唆されている技法は、加速度計の位置の決定を許容していることから、これは、取得されたノードを接続することによる、ターゲット物体の迅速な且つ概略的な「デジタル化」のために使用することが可能であり、この場合に、加速度計は、ラインと共に、配置されている。

本明細書において記述されている方法の実施形態は、いくつかの別個の要素を有するハードウェアを利用することにより、且つ／又は、少なくとも部分的に適切にプログラミングされたマイクロプロセッサを利用することにより、実装することができる。いくつかの手段を列挙している装置請求項において、これらの手段のいくつかは、１つの且つ同一の要素、コンポーネント、又はハードウェア品目により、実施することができる。特定の手段が、相互に異なる従属請求項において記述されている、或いは、異なる実施形態において記述されている、という単なる事実は、これらの手段の組合せが有利に使用されえないことを示すものではない。

本明細書において使用されている「有する（ｃｏｍｐｒｉｓｅｓ／ｃｏｍｐｒｉｓｉｎｇ）」という用語は、記述されている特徴、要素、ステップ、又はコンポーネントの存在を規定するものと解釈されるが、１つ又は複数のその他の特徴、要素、ステップ、コンポーネント、又はこれらのグループの存在又は追加を除外するものではないことを強調しておきたい。
上述の実施形態は下記のように記載され得るが、下記に限定されるものではない。
［構成１］
ターゲット物体に装着された複数のトランスデューサの、それぞれの位置及び向きを示すそれぞれの空間的構成を決定するコンピュータ実装された方法であって、
－前記複数のトランスデューサのそれぞれから、センサ信号を受け取るステップであって、前記センサ信号は、前記ターゲット物体の誘発された運動の際の前記それぞれのトランスデューサのそれぞれの運動を示している、ステップと、
－前記複数のトランスデューサの前記空間的構成を決定するステップであって、前記トランスデューサの少なくともサブセットの、前記ターゲット物体との関係における前記位置及び向きを示す前記空間的構成は、前記受け取られたセンサ信号に基づいて決定される、ステップと、
を有する方法。
［構成２］
前記トランスデューサは、それぞれの加速度計、具体的には、３軸加速度計、を有する、構成１に記載の方法。
［構成３］
入力として、前記複数のトランスデューサの前記空間的構成の一部分のみを表す情報を受け取るステップを有し、且つ、前記複数のトランスデューサの前記空間的構成を決定するステップは、前記複数のトランスデューサのすべてのトランスデューサの完全な空間的構成を確立するべく、前記空間的構成の残りの部分を決定するステップを有する、構成１又は２に記載の方法。
［構成４］
前記複数のトランスデューサは、トランスデューサの第１サブセットと、トランスデューサの第２サブセットと、を有し、且つ、前記複数のトランスデューサの前記空間的構成を決定するステップは、
－前記第１サブセットの前記トランスデューサの空間的構成を取得するステップと、
－前記受け取られたセンサ信号から、且つ、前記第１サブセットの前記トランスデューサの前記取得された空間的構成から、前記第２サブセットの前記トランスデューサの前記空間的構成を決定するステップと、
を有する、構成１から３のいずれか１項に記載の方法。
［構成５］
前記第１サブセットは、３つの３軸加速度計を有する、構成４に記載の方法。
［構成６］
第１サブセットは、３つの直交方向に沿った前記ターゲット物体の平行運動と、３つの直交軸の周りにおける回転運動と、を計測するのに十分な数のトランスデューサを有する、構成４又は５に記載の方法。
［構成７］
前記第１サブセットが十分な数のトランスデューサを含んでいるかどうかを検証するステップを有する、構成６に記載の方法。
［構成８］
前記トランスデューサは、それぞれの３軸トランスデューサを有し、且つ、検証ステップは、前記第１サブセットの前記トランスデューサのそれぞれのトランスデューサ変換行列から構築された変換行列がフルランクを有するかどうかを検証するステップを有し、それぞれのトランスデューサ変換行列は、前記ターゲット物体のローカル座標系との関係における前記対応するトランスデューサの加速度と前記トランスデューサの計測座標系の軸に沿った前記トランスデューサによって計測される加速度の間の変換を表している、構成７に記載の方法。
［構成９］
前記誘発された運動が、トランスデューサの前記第２サブセットの前記トランスデューサの前記空間的構成を決定するために十分なものであるかどうかを決定するステップを有する、構成４から８のいずれか１項に記載の方法。
［構成１０］
前記トランスデューサは、それぞれの３軸加速度計を有し、前記受け取られたセンサ信号は、前記誘発された運動の際の前記トランスデューサの計測された加速度を示しており、且つ、前記生成された運動が、前記第２サブセットの前記トランスデューサの前記空間的構成を決定するために十分なものであるかどうかを決定するステップは、
－任意選択により、前記センサ信号の低域通過フィルタリングを実行するステップと、
－前記計測された加速度を表すローカル加速度行列が既定数の主要特異値を有するかどうかを決定するステップと、
を有する、構成９に記載の方法。
［構成１１］
前記第２サブセットの前記トランスデューサの空間的構成を決定するステップは、
－前記第１サブセットの前記トランスデューサの前記受け取られたセンサ信号及び前記受け取られた空間的構成から、基準座標系との関係における前記ターゲット物体の加速度を示す物体加速度を演算するステップと、
－前記演算された物体加速度から、且つ、前記第２サブセットのそれぞれのトランスデューサの前記受け取られたセンサ信号から、前記第２サブセットの前記トランスデューサの前記空間的構成を演算するステップと、
を有する、構成４から１０のいずれか１項に記載の方法。
［構成１２］
前記第１サブセットの前記トランスデューサの前記空間的構成を取得するステップは、
－前記第１サブセットの前記トランスデューサの間の空間的関係の組を示す情報を受け取るステップと、
－前記受け取られた情報から、且つ、前記受け取られたセンサ信号から、前記第１サブセットの前記トランスデューサの前記空間的構成を演算するステップと、
を有する、構成４から１１のいずれか１項に記載の方法。
［構成１３］
前記ターゲット物体は、そのモーションが物体支持部によって制限されている、変形可能な物体であり、且つ、前記方法は、
－前記第１サブセットの前記トランスデューサの前記空間的構成から、且つ、前記受け取られたセンサ信号から、前記ターゲット物体の時間依存加速度係数を演算するステップと、
－前記演算された時間依存加速度係数から、且つ、前記受け取られたセンサ信号から、前記第２サブセットの前記トランスデューサの前記空間的構成を演算するステップと、
を有する、構成４から１２のいずれか１項に記載の方法。
［構成１４］
前記誘発された運動の撓みモードと関連する周波数を抑圧するように、前記記録された信号をフィルタリングするステップを有する、構成１から１３のいずれか１項に記載の方法。
［構成１５］
コンピュータ実装された計測プロセスであって、
－構成１から１４のいずれか１項に記載の前記方法の前記ステップを実行することにより、ターゲット物体に装着された複数のトランスデューサの空間的構成を決定するステップと、
－前記複数のトランスデューサから更なるセンサ信号を受け取るステップであって、前記受け取られた更なるセンサ信号は、前記ターゲット物体の誘発された運動の際の前記それぞれのトランスデューサのそれぞれの運動を示している、ステップと、
－前記受け取られた更なるセンサ信号に基づいて構造解析計算を実行するステップと、
を有する、プロセス。
［構成１６］
構成１から１５のいずれか１項に記載の前記方法の前記ステップを実行するように構成された、データ処理システム。
［構成１７］
コンピュータプログラムであって、前記コンピュータプログラムがデータ処理システムによって実行された際に、前記データ処理システムが、構成１から１５のいずれか１項に記載の前記方法の前記ステップを実行するようにするべく構成されている、コンピュータプログラム。
［構成１８］
計測システムであって、
－構成１７に記載のデータ処理システムと、
－それぞれの無線又は有線信号接続を介して、ターゲット物体に装着された複数のトランスデューサのそれぞれのトランスデューサ出力信号に結合されるように構成された複数の入力チャネルと、
－例えば、手動的データ入力による、或いは、無線又は有線データ通信リンクを介した別の装置からの、前記第１サブセットの前記トランスデューサのそれぞれの空間的座標データの受け取りのために構成されたデータ入力インターフェイスと、
を有する、システム。
［構成１９］
トランスデューサ組立体であって、
－構成１８に記載の計測システムと、
－ターゲット物体に跨って分散された複数の既定の計測場所において取り付け可能である複数のトランスデューサと、
を有する、組立体。

Claims (18)

1. ターゲット物体に装着された複数のトランスデューサの、それぞれの位置及び向きを示すそれぞれの空間的構成を決定するコンピュータ実装された方法であって、
   －前記複数のトランスデューサのそれぞれから、センサ信号を受け取るステップであって、前記センサ信号は、前記ターゲット物体の誘発された運動の際の前記それぞれのトランスデューサのそれぞれの運動を示している、ステップと、
   －前記複数のトランスデューサの前記空間的構成を決定するステップであって、前記トランスデューサの少なくともサブセットの、前記ターゲット物体との関係における前記位置及び向きを示す前記空間的構成は、前記受け取られたセンサ信号に基づいて決定される、ステップと、
   を有し、
   前記複数のトランスデューサは、トランスデューサの第１サブセットと、トランスデューサの第２サブセットと、を有し、且つ、前記複数のトランスデューサの前記空間的構成を決定するステップは、
   －前記第１サブセットの前記トランスデューサの空間的構成を取得するステップと、
   －前記受け取られたセンサ信号から、且つ、前記第１サブセットの前記トランスデューサの前記取得された空間的構成から、前記第２サブセットの前記トランスデューサの前記空間的構成を決定するステップと、
   を有する、方法。
2. 前記トランスデューサは、それぞれの加速度計を有する、請求項１に記載の方法。
3. 入力として、前記複数のトランスデューサの前記空間的構成の一部分のみを表す情報を受け取るステップを有し、且つ、前記複数のトランスデューサの前記空間的構成を決定するステップは、前記複数のトランスデューサのすべてのトランスデューサの完全な空間的構成を確立するべく、前記空間的構成の残りの部分を決定するステップを有する、請求項１又は２に記載の方法。
4. 前記第１サブセットは、３つの３軸加速度計を有する、請求項１から３のいずれか１項に記載の方法。
5. 第１サブセットは、３つの直交方向に沿った前記ターゲット物体の平行運動と、３つの直交軸の周りにおける回転運動と、を計測するのに十分な数のトランスデューサを有する、請求項１から４のいずれか１項に記載の方法。
6. 前記第１サブセットが十分な数のトランスデューサを含んでいるかどうかを検証するステップを有する、請求項５に記載の方法。
7. 前記トランスデューサは、それぞれの３軸トランスデューサを有し、且つ、検証ステップは、前記第１サブセットの前記トランスデューサのそれぞれのトランスデューサ変換行列から構築された変換行列がフルランクを有するかどうかを検証するステップを有し、それぞれのトランスデューサ変換行列は、前記ターゲット物体のローカル座標系との関係における前記対応するトランスデューサの加速度と前記トランスデューサの計測座標系の軸に沿った前記トランスデューサによって計測される加速度の間の変換を表している、請求項６に記載の方法。
8. 前記誘発された運動が、トランスデューサの前記第２サブセットの前記トランスデューサの前記空間的構成を決定するために十分なものであるかどうかを決定するステップを有する、請求項１から７のいずれか１項に記載の方法。
9. 前記トランスデューサは、それぞれの３軸加速度計を有し、前記受け取られたセンサ信号は、前記誘発された運動の際の前記トランスデューサの計測された加速度を示しており、且つ、前記誘発された運動が、前記第２サブセットの前記トランスデューサの前記空間的構成を決定するために十分なものであるかどうかを決定するステップは、
   －前記センサ信号の低域通過フィルタリングを実行するステップと、
   －前記計測された加速度を表すローカル加速度行列が既定数の主要特異値を有するかどうかを決定するステップと、
   を有する、請求項８に記載の方法。
10. 前記第２サブセットの前記トランスデューサの空間的構成を決定するステップは、
    －前記第１サブセットの前記トランスデューサの前記受け取られたセンサ信号及び前記受け取られた空間的構成から、基準座標系との関係における前記ターゲット物体の加速度を示す物体加速度を演算するステップと、
    －前記演算された物体加速度から、且つ、前記第２サブセットのそれぞれのトランスデューサの前記受け取られたセンサ信号から、前記第２サブセットの前記トランスデューサの前記空間的構成を演算するステップと、
    を有する、請求項１から９のいずれか１項に記載の方法。
11. 前記第１サブセットの前記トランスデューサの前記空間的構成を取得するステップは、
    －前記第１サブセットの前記トランスデューサの間の空間的関係の組を示す情報を受け取るステップと、
    －前記受け取られた情報から、且つ、前記受け取られたセンサ信号から、前記第１サブセットの前記トランスデューサの前記空間的構成を演算するステップと、
    を有する、請求項１から１０のいずれか１項に記載の方法。
12. 前記ターゲット物体は、そのモーションが物体支持部によって制限されている、変形可能な物体であり、且つ、前記方法は、
    －前記第１サブセットの前記トランスデューサの前記空間的構成から、且つ、前記受け取られたセンサ信号から、前記ターゲット物体の時間依存加速度係数を演算するステップと、
    －前記演算された時間依存加速度係数から、且つ、前記受け取られたセンサ信号から、前記第２サブセットの前記トランスデューサの前記空間的構成を演算するステップと、
    を有する、請求項１から１１のいずれか１項に記載の方法。
13. 前記誘発された運動の撓みモードと関連する周波数を抑圧するように、前記受け取られたセンサ信号をフィルタリングするステップを有する、請求項１から１２のいずれか１項に記載の方法。
14. コンピュータ実装された計測プロセスであって、
    －請求項１から１３のいずれか１項に記載の前記方法の前記ステップを実行することにより、ターゲット物体に装着された複数のトランスデューサの空間的構成を決定するステップと、
    －前記複数のトランスデューサから更なるセンサ信号を受け取るステップであって、前記受け取られた更なるセンサ信号は、前記ターゲット物体の誘発された運動の際の前記それぞれのトランスデューサのそれぞれの運動を示している、ステップと、
    －前記受け取られた更なるセンサ信号に基づいて構造解析計算を実行するステップと、
    を有する、プロセス。
15. 請求項１から１４のいずれか１項に記載の前記方法の前記ステップを実行するように構成された、データ処理システム。
16. コンピュータプログラムであって、前記コンピュータプログラムがデータ処理システムによって実行された際に、前記データ処理システムが、請求項１から１４のいずれか１項に記載の前記方法の前記ステップを実行するようにするべく構成されている、コンピュータプログラム。
17. 計測システムであって、
    －請求項１６に記載のデータ処理システムと、
    －それぞれの無線又は有線信号接続を介して、ターゲット物体に装着された複数のトランスデューサのそれぞれのトランスデューサ出力信号に結合されるように構成された複数の入力チャネルと、
    －手動的データ入力による、或いは、無線又は有線データ通信リンクを介した別の装置からの、前記第１サブセットの前記トランスデューサのそれぞれの空間的座標データの受け取りのために構成されたデータ入力インターフェイスと、
    を有する、システム。
18. トランスデューサ組立体であって、
    －請求項１７に記載の計測システムと、
    －ターゲット物体に跨って分散された複数の既定の計測場所において取り付け可能である複数のトランスデューサと、
    を有する、組立体。

Images