JP2004132901A - Component force measuring apparatus - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To accurately determine the component forces of an unknown load acting on an object to be measured in consideration of the nonlinear relation between deformation and load. <P>SOLUTION: A distortion gauge is mounted to the object to be measured (S301). A load acting only on the direction of each component force is provided for the object to be measured as increasing its magnitude. A computing device displays a graph of the detection output of the distortion gauge and load values on a display device (S302). A workman performs the processing of dividing the graph into a plurality of approximately linear regions by viewing the graph (S303). A load acting on the direction of each component force and corresponding to a unit vector is provided for the object to be measured in each linear region (S304-S306). The computing device computes a deformation-load transformation determinant in each of the linear regions α, β, and γ (S307). The computing device computes each unknown component force on the basis of the determinant. <P>COPYRIGHT: (C)2004,JPO

Description

【０００１】
【発明の属する技術分野】
本発明は、被測定物に作用する未知な分力を求める分力測定装置に関する。
【０００２】
【従来の技術】
従来より、被測定物に作用する未知な荷重による分力を計測するために、分力測定装置が用いられている。このような分力測定装置の一例が、特開平１１−１０８７８０号公報（特許文献１）に記載されている。なお、被測定物としては一般に存在する様々な物を対象とすることができる。例えば、自動車であればエンジン、トランスミッションなどの内部部品などである。
【０００３】
この文献の分力測定装置では、複数の変形検出器を被測定物の任意の位置に取り付け、その取付状態で較正を行って変形−荷重伝達係数行列を記憶し、被測定物に荷重が加わると、変形−荷重伝達係数行列を用いて各分力を求めている。
【０００４】
この分力測定装置では、次のように較正を行っている。まず、変形検出器の数に対応した数の異なる既知の分力を被測定物に与える。そして、各分力毎の較正変形データを入力して荷重−変形伝達係数行列を作成し、荷重−変形伝達係数行列を変形−荷重伝達係数行列に逆変換する。
【０００５】
【特許文献１】
特開平１１−１０８７８０号公報
【０００６】
【発明が解決しようとする課題】
上述した従来技術の分力測定装置では、較正により得られる変形−荷重伝達係数行列は一つだけであった。しかし、被測定物の変形は、荷重に対して非線形の特性を持つため、変形−荷重伝達係数行列を一つだけ求め、この変形−荷重伝達係数行列に基づいて被測定物に作用する分力を演算しても、正確な変形量を計算できないといった問題があった。
【０００７】
なお、変形と荷重の関係は、あらゆる系において一般的に非線形となる。但し、被測定物がゴムなどを含む場合には、特に非線形関係が顕著となる。自動車であれば、ゴムを材料とするブッシュ、ショックアブソーバを含む場合である。
【０００８】
本発明は、上述の課題に鑑みてなされたものであり、被測定物に作用する未知な荷重の分力を、変形と荷重の非線形関係を考慮して、より正確に求めることができる分力測定装置を提供することを目的とする。
【０００９】
【課題を解決するための手段】
本発明に係る分力測定装置は、被測定物に作用する未知な分力の大きさを求める分力測定装置であって、被測定物の変形を検出する複数の変形検出器と、変形と分力の非線形関係に関するデータを保存する記憶装置と、前記変形検出器の検出出力と、前記記憶装置に保存された非線形関係データに基づき、前記被測定物に作用する分力の大きさを求める演算装置と、を備えるものである。
【００１０】
これにより、変形と荷重の非線形関係を考慮して、被測定物に作用する未知な分力をより正確に求めることができる。ここで、変形検出器とは、被測定物の変形を検出する様々なものを用いることができる。例えば、一般的な荷重検出器の他、歪みゲージなどを用いることができ、これらの組み合わせでもよい。また、記憶装置に保存される変形と分力の非線形関係に関するデータとは、例えば行列や関数などの形式で保存されるデータであり、所定の変形に対応した分力を一義的に算出可能なデータである。
【００１１】
上述の分力測定装置は、より具体的には、次のような態様とすることができる。
【００１２】
つまり、前記記憶装置は、非線形な変形−荷重曲線における各線形領域に対応付けられた複数の変形−荷重変換行列式を保存し、前記演算装置は、前記変形検出器の検出出力に基づき、前記複数の変形−荷重変換行列式からいずれか一つを選択し、選択された変形−荷重変換行列式に、前記検出出力を代入して、前記被測定物に作用する分力の大きさを求める。このような態様においては、全体として非線形な変形−荷重曲線を複数の線形領域に分割し、各線形領域に対応した変形−荷重変換行列式に基づいて、分力を演算するため、変形−荷重変換行列式と検出出力の列ベクトルの積を求めるだけの単純な演算するだけで済む。つまり、演算処理が簡易であり、迅速に分力を求めることができる。
【００１３】
また、分力測定装置は、次のような態様とすることもできる。
【００１４】
前記記憶装置は、変形量を変数とする非線形関係式を成分として含む変形−荷重変換行列式を保存し、前記演算装置は、前記変形−荷重変換行列式に、前記検出出力を代入して、前記被測定物に作用する分力の大きさを求める。ここで、前記記憶装置が保存する関係式は、最小自乗法などの近似方法により求められた近似式としてもよい。この態様においては、変形−荷重伝達行列が非線形関係式を含むため、より実際に近い分力の値を求めることができる。
【００１５】
【発明の実施の形態】
以下に、図面を参照して、本発明の実施形態について説明する。
【００１６】
第１の実施形態．
（１）分力測定装置の構成
第１の実施形態に係る分力測定装置について説明する。図１は、第１の実施形態に係る分力測定装置１０の概略構成を示すブロック図である。
【００１７】
第１の実施形態に係る分力測定装置１０は、被測定物に発生する変形量を検出する変形検出器１２、変形検出器１２の出力を増幅する増幅器１４、増幅された出力をアナログ値からデジタル値に変換するＡ／Ｄ変換器１６、変形−荷重変換行列式の保存する記憶装置１８、変形−荷重変換行列式の演算処理、分力の演算処理などを行う演算装置（ＣＰＵ）２０、及び表示装置２２が接続されて構成されている。
【００１８】
変形検出器１２は、本実施形態では歪みゲージであり、被測定物の表面に取り付けられる。荷重を受ける被測定物がエンジン２４である場合には、図２に示されるように、エンジンマウント２６に歪みゲージ１２Ｘ，Ｙ，Ｚを取り付ければよい。
【００１９】
なお、変形検出器１２としては、歪みゲージの他、一般的な様々な荷重検出器を用いることができる。また、複数の種類の荷重検出器を組み合わせて用いてもよい。一般的に変形検出器１２の個数ｍは、測定する分力の数ｎの２乗以上の数が必要である。但し、以下の説明では、説明を簡略化するために歪みゲージの個数を、実際に必要な個数より少なくしている。
【００２０】
（２）変形−荷重変換行列式を求める方法
次に、図３のフローチャートを参照して、荷重−変形伝達係数行列を求める方法について説明する。
【００２１】
まず、被測定物に歪みゲージ１２Ｘ，Ｙ，Ｚを取り付ける（Ｓ３０１）。次に、被測定物に対して各分力方向にのみ作用する荷重を、その大きさを徐々に増加させて与える。このような荷重は、荷重を負荷するための一般的な装置を用いることにより与えることができる。本実施形態では、被測定物の前後方向（Ｘ軸方向）、左右方向（Ｙ軸方向）、上下方向（Ｚ軸方向）を分力方向としているので、始めにＸ軸方向に負荷Ｆｘを徐々にその大きさを増加させて与え、その後Ｙ軸方向、Ｚ軸方向についても、負荷Ｆｘ，Ｆｙを徐々にその大きさを増加させて与える。このとき、この荷重により歪みゲージ１２Ｘ，Ｙ，Ｚの取り付け位置が変形し、変形が抵抗値変化として各歪みゲージ１２Ｘ，Ｙ，Ｚで検出される。演算装置２０は、各歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力を取得すると共に、荷重負荷装置から荷重値を入力し、歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力と荷重値を関連づけて記憶装置１８に保存する。そして、この保存されたデータから、図４（ａ）〜（ｃ）に示されるような、各分力Ｆｘ，Ｆｙ，Ｆｚと歪みゲージ検出出力の関係を示すグラフを作成し、表示装置２２に表示する（Ｓ３０２）。
【００２２】
次に、作業者は、これらのグラフを見て、グラフを複数のほぼ線形な領域に分割する処理を行う（Ｓ３０３）。図４の場合、具体的には、作業者がグラフに３つの線形領域α，β，γがあることを目視で認識し、グラフが屈曲し線形性が失われるポイントαｍａｘ，βｍａｘをキーボード、マウス等の入力手段を用いて、演算装置２０に入力することにより行う。入力が行われると、演算装置２０は、ポイントの名称（αｍａｘ，βｍａｘなど）と、そのポイントにおける荷重値Ｆｘ（αｍａｘ），Ｆｙ（αｍａｘ），Ｆｚ（αｍａｘ），Ｆｘ（βｍａｘ），Ｆｙ（βｍａｘ），Ｆｚ（βｍａｘ）と、そのポイントにおける各歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力、つまりポイントαｍａｘについては検出出力｛ε１１（αｍａｘ），ε１２（αｍａｘ），ε１３（αｍａｘ）｝，｛ε２１（αｍａｘ），ε２２（αｍａｘ），ε３２（αｍａｘ）｝，｛ε３１（αｍａｘ），ε３２（αｍａｘ），ε３３（αｍａｘ）｝を関連付けて、またポイントβｍａｘについては｛ε１１（βｍａｘ），ε１２（βｍａｘ），ε１３（βｍａｘ）｝，｛ε２１（βｍａｘ），ε２２（βｍａｘ），ε３２（βｍａｘ）｝，｛ε３１（βｍａｘ），ε３２（βｍａｘ），ε３３（βｍａｘ）｝を関連付けて、記憶装置１８に保存する。なお、図４の例では、各グラフの屈曲ポイントが互いに対応しているため、Ｘ軸についてのみグラフの線形領域の設定を行っているが、各軸についてグラフを線形領域に分割する処理を行ってもよい。
【００２３】
次に、演算装置２０を用いて、各線形領域α，β，γについて、変形−荷重伝達係数行列を演算する処理を行う。このために先ず、各線形領域において被測定物に各分力方向にのみ作用し、単位ベクトルに相当する荷重を与える。すなわち、線形領域αについては、
（ｘ，ｙ，ｚ）＝（１，０，０）、（０，１，０）、（０，０，１）
線形領域βについては、
（ｘ，ｙ，ｚ）＝（Ｆｘ（αｍａｘ）＋１，０，０）、（０，Ｆｙ（αｍａｘ）＋１，０）、（０，０，Ｆｙ（αｍａｘ）＋１）
線形領域γについては、
（ｘ，ｙ，ｚ）＝（Ｆｘ（βｍａｘ）＋１，０，０）、（０，Ｆｙ（βｍａｘ）＋１，０）、（０，０，Ｆｙ（βｍａｘ）＋１）
の荷重を与える。そして、これらの荷重を被測定物に与えたときの、歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力のデータを得る（Ｓ３０４〜Ｓ３０６）。
【００２４】
次に、演算装置２０が、得られた検出出力に基づき、各線形領域α，β，γにおける変形−荷重変換行列式を演算するために、以下に説明する処理を行う（Ｓ３０７）。
【００２５】
線形領域αにおいて、上記３つの荷重を与えたときの歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力がそれぞれ｛ａ１１（α），ａ１２（α），ａ３１（α）｝、｛ａ２１（α），ａ２２（α），ａ３２（α）｝、｛ａ３１（α），ａ３２（α），ａ３３（α）｝である場合には、荷重−変形伝達係数行列［Ａ（α）］は次の式１で表される。
【数１】

【００２６】
この荷重−変形伝達係数行列［Ａ（α）］を逆変換することにより、変形−荷重伝達係数行列［Ｂ（α）］を求めることができる。
【数２】

【００２７】
線形領域αにおいては、この変形−荷重伝達係数行列［Ｂ（α）］に各歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力の列ベクトルを積算することにより、各分力（Ｆｘ，Ｆｙ，Ｆｚ）を求めることができる。つまり、この積算した行列式が、線形領域αの変形−荷重変換行列式である。
【００２８】
また、線形領域βにおいて、上述の３つの荷重を与えたときの歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力がそれぞれ｛ａ１１（β），ａ１２（β），ａ１３（β）｝、｛ａ２１（β），ａ２２（β），ａ３２（β）｝、｛ａ３１（β），ａ３２（β），ａ３３（β）｝である場合には、荷重−変形伝達係数行列［Ａ（β）］は次の式３で表される。
【数３】

【００２９】
この荷重−変形伝達係数行列［Ａ（β）］を逆変換することにより、変形−荷重伝達係数行列［Ｂ（β）］を求めることができる。
【数４】

【００３０】
線形領域βにおいては、次式の変形−荷重変換行列式に示すように、変形−荷重伝達係数行列［Ｂ（β）］を用いた次式に各歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力（ε１，ε２，ε３）を代入することにより、荷重値（Ｆｘ，Ｆｙ，Ｆｚ）を求めることができる。
【数５】

【００３１】
また、線形領域γについても、線形領域βと同様に、変形−荷重変換行列式を求めることができる。線形領域が多数あるときには、一般的に線形領域βと同様の方法で変形−荷重変換行列式を求めればよい。
【００３２】
（３）実働荷重の演算処理
次に、被測定物に未知の外力が作用したときに、演算装置２０が上述の処理により求められた変形−荷重変換行列式を用いて外力の各分力を演算する処理について説明する。図５は、各分力を演算する処理のフローチャートである。
【００３３】
まず、演算装置２０が各歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力を入力する（Ｓ５０１）。
【００３４】
この検出値の出力に基づき、被測定物に作用する外力がいずれの線形領域に属するかを判別する。本実施形態の場合、歪みゲージＡの検出出力を、記憶装置１８に保存された屈曲ポイントαｍａｘ，βｍａｘ等における検出値と比較し、線形領域を選択する（Ｓ５０２）。
【００３５】
次に、演算装置２０は、選択された線形領域の変形−荷重変換行列式を、記憶装置１８から読み出し、この変形−荷重変換行列式（式５）に、ステップＳ５０１で入力された歪みゲージ検出出力の列ベクトル｛ε１，ε２，ε３｝を代入する（Ｓ５０３）。この結果、各分力が求められる。演算装置２０は、これらの各分力｛Ｆｘ，Ｆｙ，Ｆｚ｝を表示装置２２にリアルタイムで表示する（Ｓ５０４）。
【００３６】
本実施形態では、上述したように、各線形領域ごとに保存された変形−荷重変換行列式を、記憶装置１８から読み出し、この行列式に基づき、被測定物に作用する未知の外力の各分力を演算するため、荷重−変形特性が非線形であっても、真の値に近い分力を求めることができる。また、本実施形態では、変形−荷重変換行列式が簡易な行列式であるため、演算装置２０は歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力を取得後、速やかに各分力を演算することができる。
【００３７】
第２の実施形態．
次に、第２の実施形態に係る分力測定装置について説明する。第２の実施形態に係る分力測定装置の装置構成は、第１の実施形態の分力測定装置と同じであるため、説明を省略する。第２の実施形態では、「変形−荷重変換行列式を求める方法」が、第１の実施形態と異なる。以下に、図６のフローチャートを参照して、第２の実施形態に係る「変形−荷重変換行列式を求める方法」について説明する。
【００３８】
まず、被測定物に歪みゲージ１２Ｘ，Ｙ，Ｚを取り付ける（Ｓ６０１）。次に、被測定物に対して各分力方向にのみ作用する荷重、すなわち上下方向にのみ作用する荷重Ｆｘ、前後方向にのみ作用する荷重Ｆｙ、左右方向にのみ作用する荷重Ｆｚを、その大きさを所定値ごと徐々に増加させて与え、負荷検定を行う（Ｓ６０２〜Ｓ６０４）。各方向の荷重が与えられると、各歪みゲージ１２Ｘ，Ｙ，Ｚから検出出力｛ａ１１，ａ１２，ａ１３｝、｛ａ２１，ａ２２，ａ２３｝、｛ａ３１，ａ３２，ａ３３｝が得られる。なお、これらの各検出出力ａ１１〜ａ３３の値は、荷重Ｆｘ，Ｆｙ，Ｆｚの値により変化する値である。演算装置２０は、各歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力と共に、荷重負荷装置から荷重値を入力し、歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力と荷重値を関連づけて記憶装置１８に保存する。
【００３９】
次に、演算装置２０は、各検出出力ａ１１〜ａ３３について、最小自乗法を用いて、荷重値Ｆｘ，Ｆｙ，Ｆｚを変数とする検出出力値の非線形な関係式ａ１１（Ｆｘ），ａ１２（Ｆｘ），ａ１３（Ｆｘ），ａ２１（Ｆｙ），ａ２２（Ｆｙ），ａ２３（Ｆｙ），ａ３１（Ｆｚ），ａ３２（Ｆｚ），ａ３３（Ｆｚ）を求める。これらの関係式は、荷重値Ｆｘ，Ｆｙ，Ｆｚの多次式とすればよい。このように本実施形態においては、検出出力値の関係式として、最小自乗法により近似式を求めているが、近似を用いない多次式としてもよい。各関係式ａ１１（Ｆｘ）〜ａ３３（Ｆｚ）を求めることにより、次式の荷重−変形伝達係数行列［Ａ］を求めることができる。
【数６】

【００４０】
さらに、この荷重−変形伝達係数行列［Ａ］を逆変換して、変形−荷重変換行列式［Ｂ］、すなわち変形−分力変換行列式を求めることができる（Ｓ６０５）。
【００４１】
被測定物に作用する未知の荷重の分力を測定するときには、演算装置２０は未知な荷重により生じた各歪みゲージ１２Ｘ，Ｙ，Ｚの検出出力を取得し、変形−荷重伝達係数行列［Ｂ］の各成分に検出出力を代入し、さらに変形−荷重伝達係数行列［Ｂ］と検出出力の列ベクトルを積算する。これにより、未知の荷重の各分力を求めることができる。
【００４２】
以上に説明した第２の実施形態に係る分力測定装置では、変形−荷重変換行列式が、変形量を変数とする非線形関係式を成分として含む行列であるため、第１の実施形態と同様に、変形と分力の非線形関係を補正して、真の値に近い分力を求めることができる。
【００４３】
以上、本発明の好適な実施形態について説明したが、本発明は上記実施形態に限定されるものではなく、等価な範囲で様々な変形が可能である。例えば、第１の実施形態に係るグラフを複数の線形領域に分割する方法と、第２の実施形態に係る最小自乗法を用いた方法を組み合わせた変形−荷重変換行列式を求め、未知な外力の分力を演算してもよい。
【００４４】
【発明の効果】
本発明では、被測定物に作用する荷重と、この荷重による変形の特性が非線形である場合でも、実際の値に近い分力を求めることができる。
【図面の簡単な説明】
【図１】第１の実施形態に係る分力測定装置１０の構成を示すブロック図である。
【図２】被測定物２４，２６に歪みゲージ１２Ｘ，Ｙ，Ｚが取り付けられた様子を示す斜視図である。
【図３】第１の実施形態に係る分力測定装置１０の較正処理を示すフローチャートである。
【図４】第１の実施形態に係る分力測定装置１０の較正処理を説明するための説明図である。
【図５】第１の実施形態に係る分力測定装置１０の分力演算処理を示すフローチャートである。
【図６】第２の実施形態に係る分力測定装置１０の較正処理を示すフローチャートである。
【符号の説明】
１０　分力測定装置、１２，１２Ｘ，１２Ｙ，１２Ｚ　変形検出器（歪みゲージ）、１４　増幅器、１６　Ａ／Ｄ変換器、１８　記憶装置、２０　演算装置、２２　表示装置、２４　エンジン、２６　マウント。[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a component force measuring device for determining an unknown component force acting on an object to be measured.
[0002]
[Prior art]
Conventionally, a component force measuring device has been used to measure a component force due to an unknown load acting on an object to be measured. An example of such a component force measuring device is described in Japanese Patent Application Laid-Open No. H11-108780 (Patent Document 1). In addition, various objects that are generally present can be used as the object to be measured. For example, in the case of an automobile, internal components such as an engine and a transmission are used.
[0003]
In the component force measuring device of this document, a plurality of deformation detectors are attached to an arbitrary position of an object to be measured, calibration is performed in the attached state, a deformation-load transfer coefficient matrix is stored, and a load is applied to the object to be measured. And each component force is obtained using the deformation-load transfer coefficient matrix.
[0004]
In this component force measuring device, calibration is performed as follows. First, a known number of different component forces corresponding to the number of deformation detectors are applied to the device under test. Then, the load-deformation transfer coefficient matrix is created by inputting the calibration deformation data for each component force, and the load-deformation transfer coefficient matrix is inversely transformed into the deformation-load transfer coefficient matrix.
[0005]
[Patent Document 1]
JP-A-11-108780
[Problems to be solved by the invention]
In the above-described prior-art component force measuring device, only one deformation-load transfer coefficient matrix is obtained by calibration. However, since the deformation of the DUT has a non-linear characteristic with respect to the load, only one deformation-load transfer coefficient matrix is obtained, and the component force acting on the DUT based on the deformation-load transfer coefficient matrix is determined. , There is a problem that an accurate deformation amount cannot be calculated.
[0007]
Note that the relationship between deformation and load is generally nonlinear in any system. However, when the object to be measured includes rubber or the like, the non-linear relationship is particularly significant. In the case of an automobile, the case includes a bush and a shock absorber made of rubber.
[0008]
The present invention has been made in view of the above-described problems, and a component force of an unknown load acting on an object to be measured can be more accurately determined in consideration of a nonlinear relationship between deformation and load. It is an object to provide a measuring device.
[0009]
[Means for Solving the Problems]
A component force measuring device according to the present invention is a component force measuring device that determines the magnitude of an unknown component force acting on an object to be measured, and includes a plurality of deformation detectors that detect deformation of the object to be measured, A storage device for storing data relating to the non-linear relationship of component forces, a detection output of the deformation detector, and a magnitude of a component force acting on the device under test based on the non-linear relationship data stored in the storage device. And an arithmetic device.
[0010]
Thus, the unknown component force acting on the object to be measured can be more accurately determined in consideration of the nonlinear relationship between the deformation and the load. Here, various types of deformation detectors that detect deformation of an object to be measured can be used. For example, in addition to a general load detector, a strain gauge or the like can be used, and a combination thereof may be used. Further, the data on the nonlinear relationship between the deformation and the component force stored in the storage device is, for example, data stored in a form such as a matrix or a function, and the component force corresponding to a predetermined deformation can be uniquely calculated. Data.
[0011]
More specifically, the above-described component force measuring device can be configured as follows.
[0012]
That is, the storage device stores a plurality of deformation-load conversion determinants associated with each linear region in the non-linear deformation-load curve, and the arithmetic device is configured to execute the processing based on a detection output of the deformation detector. Any one of a plurality of deformation-load conversion determinants is selected, and the detected output is substituted into the selected deformation-load conversion determinant to determine the magnitude of the component force acting on the object to be measured. . In such an embodiment, the deformation-load curve is divided into a plurality of linear regions as a whole, and the component force is calculated based on the deformation-load conversion determinant corresponding to each linear region. All that is required is a simple operation of finding the product of the transformation determinant and the column vector of the detection output. That is, the arithmetic processing is simple, and the component force can be quickly obtained.
[0013]
Further, the component force measuring device may have the following modes.
[0014]
The storage device stores a deformation-load conversion determinant including, as a component, a non-linear relational expression having a deformation amount as a variable, and the arithmetic unit substitutes the detection output into the deformation-load conversion determinant, The magnitude of the component force acting on the object to be measured is determined. Here, the relational expression stored in the storage device may be an approximate expression obtained by an approximation method such as a least square method. In this embodiment, since the deformation-load transfer matrix includes a non-linear relational expression, it is possible to obtain a more realistic component force value.
[0015]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, embodiments of the present invention will be described with reference to the drawings.
[0016]
First embodiment.
(1) Configuration of Component Force Measurement Apparatus A component force measurement apparatus according to the first embodiment will be described. FIG. 1 is a block diagram illustrating a schematic configuration of a component force measuring device 10 according to the first embodiment.
[0017]
The component force measuring device 10 according to the first embodiment includes a deformation detector 12 that detects a deformation amount generated in an object to be measured, an amplifier 14 that amplifies an output of the deformation detector 12, and an amplified output obtained from an analog value. An A / D converter 16 for converting to a digital value, a storage device 18 for storing a deformation-load conversion determinant, a calculation device (CPU) 20 for performing a calculation process of a deformation-load conversion determinant, a component force calculation process, And the display device 22 are connected.
[0018]
The deformation detector 12 is a strain gauge in the present embodiment, and is attached to the surface of the object to be measured. When the object to be loaded is the engine 24, the strain gauges 12X, 12Y, 12Z may be attached to the engine mount 26 as shown in FIG.
[0019]
As the deformation detector 12, various general load detectors can be used in addition to the strain gauge. Further, a plurality of types of load detectors may be used in combination. Generally, the number m of the deformation detectors 12 needs to be equal to or larger than the square of the number n of the component forces to be measured. However, in the following description, the number of strain gauges is set to be smaller than the actually required number for the sake of simplicity.
[0020]
(2) Method for Obtaining Deformation-Load Conversion Matrix Equation Next, a method for obtaining the load-deformation transfer coefficient matrix will be described with reference to the flowchart of FIG.
[0021]
First, the strain gauges 12X, 12Y, 12Z are attached to the measured object (S301). Next, a load acting only on the component to be measured in each component force direction is applied while gradually increasing its magnitude. Such a load can be applied by using a general device for applying a load. In the present embodiment, since the longitudinal direction (X-axis direction), the left-right direction (Y-axis direction), and the up-down direction (Z-axis direction) of the DUT are the component force directions, the load Fx is gradually increased in the X-axis direction first. The loads Fx and Fy in the Y-axis direction and the Z-axis direction are also gradually increased in size. At this time, the mounting position of the strain gauges 12X, Y, and Z is deformed by the load, and the deformation is detected by each of the strain gauges 12X, Y, and Z as a change in resistance. The arithmetic unit 20 acquires the detection output of each of the strain gauges 12X, Y, and Z, inputs a load value from the load application device, and associates the detection output of the strain gauges 12X, Y, and Z with the load value and stores the storage device 18 therein. To save. Then, a graph showing the relationship between each component force Fx, Fy, Fz and the strain gauge detection output as shown in FIGS. It is displayed (S302).
[0022]
Next, the operator looks at these graphs and performs a process of dividing the graph into a plurality of substantially linear regions (S303). In the case of FIG. 4, specifically, the operator visually recognizes that there are three linear regions α, β, and γ in the graph, and points αmax and βmax at which the graph is bent and linearity is lost are defined by a keyboard and a mouse. This is performed by inputting to the arithmetic unit 20 using input means such as. When the input is made, the arithmetic unit 20 calculates the name of the point (αmax, βmax, etc.) and the load values Fx (αmax), Fy (αmax), Fz (αmax), Fx (βmax), Fy (βmax) at the point. ), Fz (βmax) and the detection output of each strain gauge 12X, Y, Z at that point, ie, the detection output {ε11 (αmax), ε12 (αmax), ε13 (αmax)}, ｛ε21 ( αmax), ε22 (αmax), ε32 (αmax)｝, {ε31 (αmax), ε32 (αmax), ε33 (αmax)}, and for the point βmax, {ε11 (βmax), ε12 (βmax), ε13 (βmax)｝, {ε21 (βmax), ε22 (βmax), ε32 (βmax)}, ｛ε31 (βmax) ε32 (βmax), in association with ε33 (βmax)}, stored in the storage device 18. In the example of FIG. 4, since the inflection points of each graph correspond to each other, a linear region of the graph is set only for the X axis. However, a process of dividing the graph into linear regions for each axis is performed. You may.
[0023]
Next, a process of calculating a deformation-load transfer coefficient matrix for each of the linear regions α, β, and γ is performed using the arithmetic device 20. For this purpose, first, a load corresponding to a unit vector is applied to the object to be measured only in each component force direction in each linear region. That is, for the linear region α,
(X, y, z) = (1, 0, 0), (0, 1, 0), (0, 0, 1)
For the linear region β,
(X, y, z) = (Fx (αmax) +1, 0, 0), (0, Fy (αmax) +1, 0), (0, 0, Fy (αmax) +1)
For the linear region γ,
(X, y, z) = (Fx (βmax) +1, 0, 0), (0, Fy (βmax) +1, 0), (0, 0, Fy (βmax) +1)
Apply a load of Then, data on the detection output of the strain gauges 12X, Y, Z when these loads are applied to the object to be measured are obtained (S304 to S306).
[0024]
Next, the arithmetic unit 20 performs a process described below in order to calculate a deformation-load conversion determinant in each of the linear regions α, β, and γ based on the obtained detection output (S307).
[0025]
In the linear region α, the detection outputs of the strain gauges 12X, Y, and Z when the above three loads are applied are {a11 (α), a12 (α), a31 (α)}, {a21 (α), and a22, respectively. (Α), a32 (α)}, {a31 (α), a32 (α), a33 (α)}, the load-deformation transfer coefficient matrix [A (α)] is given by the following equation 1. expressed.
(Equation 1)

[0026]
By inversely transforming the load-deformation transfer coefficient matrix [A (α)], a deformation-load transfer coefficient matrix [B (α)] can be obtained.
(Equation 2)

[0027]
In the linear region α, each component (Fx, Fy, Fz) is obtained by integrating the column vector of the detection output of each strain gauge 12X, Y, Z into this deformation-load transfer coefficient matrix [B (α)]. Can be requested. That is, this integrated determinant is the deformation-load conversion determinant of the linear region α.
[0028]
In the linear region β, the detection outputs of the strain gauges 12X, Y, and Z when the above three loads are applied are {a11 (β), a12 (β), a13 (β)}, and ｛a21 (β ), A22 (β), a32 (β)}, {a31 (β), a32 (β), a33 (β)}, the load-deformation transfer coefficient matrix [A (β)] is It is expressed by Equation 3.
[Equation 3]

[0029]
By inversely transforming the load-deformation transfer coefficient matrix [A (β)], a deformation-load transfer coefficient matrix [B (β)] can be obtained.
(Equation 4)

[0030]
In the linear region β, as shown in the deformation-load conversion determinant of the following expression, the detection output of each of the strain gauges 12X, Y, and Z is expressed by the following expression using the deformation-load transfer coefficient matrix [B (β)]. By substituting ε1, ε2, ε3), the load value (Fx, Fy, Fz) can be obtained.
(Equation 5)

[0031]
Also, for the linear region γ, the deformation-load conversion determinant can be obtained in the same manner as the linear region β. When there are many linear regions, the deformation-load conversion determinant may be generally obtained by the same method as in the linear region β.
[0032]
(3) Calculation Processing of Actual Working Load Next, when an unknown external force acts on the object to be measured, the arithmetic unit 20 calculates each component of the external force using the deformation-load conversion matrix obtained by the above processing. The calculation process will be described. FIG. 5 is a flowchart of a process for calculating each component force.
[0033]
First, the arithmetic unit 20 inputs detection outputs of the respective strain gauges 12X, Y, Z (S501).
[0034]
Based on the output of the detected value, it is determined to which linear region the external force acting on the measured object belongs. In the case of the present embodiment, the detection output of the strain gauge A is compared with the detection values at the bending points αmax, βmax and the like stored in the storage device 18 to select a linear region (S502).
[0035]
Next, the arithmetic unit 20 reads out the deformation-load conversion determinant of the selected linear region from the storage device 18 and adds the strain-gauge detection determinant input in step S501 to this deformation-load conversion determinant (expression 5). The output column vector {ε1, ε2, ε3} is substituted (S503). As a result, each component is determined. The arithmetic device 20 displays these component forces {Fx, Fy, Fz} on the display device 22 in real time (S504).
[0036]
In the present embodiment, as described above, the deformation-load conversion determinant stored for each linear region is read out from the storage device 18, and based on this determinant, each of the unknown external forces acting on the DUT is calculated. Since the force is calculated, a component force close to the true value can be obtained even if the load-deformation characteristics are non-linear. Further, in the present embodiment, since the deformation-load conversion determinant is a simple determinant, the arithmetic device 20 can quickly calculate each component force after acquiring the detection output of the strain gauges 12X, Y, Z. it can.
[0037]
Second embodiment.
Next, a component force measuring device according to a second embodiment will be described. The device configuration of the component force measuring device according to the second embodiment is the same as that of the component force measuring device of the first embodiment, and a description thereof will be omitted. The second embodiment differs from the first embodiment in the “method of obtaining a transformation-load conversion determinant”. Hereinafter, with reference to the flowchart of FIG. 6, a “method of obtaining a deformation-load conversion matrix” according to the second embodiment will be described.
[0038]
First, the strain gauges 12X, 12Y, 12Z are attached to the measured object (S601). Next, the load acting on the object to be measured only in each component force direction, that is, the load Fx acting only in the vertical direction, the load Fy acting only in the front-rear direction, and the load Fz acting only in the left-right direction, The load is gradually increased by a predetermined value, and a load test is performed (S602 to S604). When a load is applied in each direction, detection outputs {a11, a12, a13}, {a21, a22, a23}, and {a31, a32, a33} are obtained from the strain gauges 12X, Y, and Z. Note that the values of these detection outputs a11 to a33 are values that change depending on the values of the loads Fx, Fy, and Fz. The arithmetic unit 20 inputs the load value from the load load device together with the detection output of each of the strain gauges 12X, Y, and Z, and stores the detected output of the strain gauges 12X, Y, and Z in the storage device 18 in association with the load value. .
[0039]
Next, the arithmetic unit 20 uses the least squares method for each of the detection outputs a11 to a33, and uses the load values Fx, Fy, and Fz as variables to obtain nonlinear relational expressions a11 (Fx) and a12 (Fx ), A13 (Fx), a21 (Fy), a22 (Fy), a23 (Fy), a31 (Fz), a32 (Fz), and a33 (Fz). These relational expressions may be polynomials of the load values Fx, Fy, Fz. As described above, in the present embodiment, as the relational expression of the detection output value, the approximate expression is obtained by the least squares method, but may be a polynomial expression without using approximation. By obtaining the relational expressions a11 (Fx) to a33 (Fz), the load-deformation transfer coefficient matrix [A] of the following expression can be obtained.
(Equation 6)

[0040]
Further, the load-deformation transfer coefficient matrix [A] is inversely transformed to obtain a deformation-load transformation determinant [B], that is, a transformation-component force transformation determinant (S605).
[0041]
When measuring the component force of the unknown load acting on the object to be measured, the arithmetic unit 20 acquires the detection output of each of the strain gauges 12X, Y, Z generated by the unknown load, and obtains the deformation-load transfer coefficient matrix [B ] Is substituted for each component, and the deformation-load transfer coefficient matrix [B] and the column vector of the detected output are integrated. Thereby, each component of the unknown load can be obtained.
[0042]
In the component force measuring device according to the second embodiment described above, the deformation-load conversion determinant is a matrix including, as components, a non-linear relational expression having the deformation amount as a variable, and thus is the same as in the first embodiment. In addition, by correcting the nonlinear relationship between the deformation and the component force, a component force close to the true value can be obtained.
[0043]
The preferred embodiment of the present invention has been described above, but the present invention is not limited to the above embodiment, and various modifications can be made within an equivalent range. For example, a transformation-load conversion determinant obtained by combining the method of dividing the graph according to the first embodiment into a plurality of linear regions and the method using the least square method according to the second embodiment is obtained, and the unknown external force is obtained. May be calculated.
[0044]
【The invention's effect】
According to the present invention, a component force close to an actual value can be obtained even when the load acting on the object to be measured and the characteristics of deformation due to the load are nonlinear.
[Brief description of the drawings]
FIG. 1 is a block diagram illustrating a configuration of a component force measuring device 10 according to a first embodiment.
FIG. 2 is a perspective view showing a state where strain gauges 12X, Y, and Z are attached to DUTs 24 and 26;
FIG. 3 is a flowchart illustrating a calibration process of the component force measuring device 10 according to the first embodiment.
FIG. 4 is an explanatory diagram for explaining a calibration process of the component force measuring device 10 according to the first embodiment.
FIG. 5 is a flowchart showing a component force calculation process of the component force measuring device 10 according to the first embodiment.
FIG. 6 is a flowchart illustrating a calibration process of the component force measuring device 10 according to the second embodiment.
[Explanation of symbols]
10 component force measuring device, 12, 12X, 12Y, 12Z deformation detector (strain gauge), 14 amplifier, 16 A / D converter, 18 storage device, 20 computing device, 22 display device, 24 engine, 26 mount.

Claims (4)

被測定物に作用する未知な分力の大きさを求める分力測定装置であって、
被測定物の変形を検出する複数の変形検出器と、
変形と分力の非線形関係に関するデータを保存する記憶装置と、
前記変形検出器の検出出力と、前記記憶装置に保存された非線形関係データに基づき、前記被測定物に作用する分力の大きさを求める演算装置と、
を備えることを特徴とする分力測定装置。A component force measuring device for determining the magnitude of an unknown component force acting on an object to be measured,
A plurality of deformation detectors for detecting deformation of the device under test;
A storage device for storing data related to the non-linear relationship between deformation and component force,
An arithmetic unit that calculates the magnitude of the component force acting on the device under test based on the detection output of the deformation detector and the non-linear relationship data stored in the storage device;
A component force measuring device comprising: 請求項１に記載の分力測定装置であって、
前記記憶装置は、非線形な変形−荷重曲線における各線形領域に対応付けられた複数の変形−荷重変換行列式を保存し、
前記演算装置は、前記変形検出器の検出出力に基づき、前記複数の変形−荷重変換行列式からいずれか一つを選択し、選択された変形−荷重変換行列式に、前記検出出力を代入して、前記被測定物に作用する分力の大きさを求めることを特徴とする分力測定装置。The component force measuring device according to claim 1,
The storage device stores a plurality of deformation-load conversion determinants associated with each linear region in the nonlinear deformation-load curve,
The arithmetic unit selects one of the plurality of deformation-load conversion determinants based on the detection output of the deformation detector, and substitutes the detection output for the selected deformation-load conversion determinant. A force component acting on the object to be measured. 請求項１に記載の分力測定装置であって、
前記記憶装置は、変形量を変数とする非線形関係式を成分として含む変形−荷重変換行列式を保存し、
前記演算装置は、前記変形−荷重変換行列式に、前記検出出力を代入して、前記被測定物に作用する分力の大きさを求めることを特徴とする分力測定装置。The component force measuring device according to claim 1,
The storage device stores a deformation-load conversion determinant including, as a component, a nonlinear relational expression having a deformation amount as a variable,
The component force measuring device, wherein the arithmetic device substitutes the detection output into the deformation-load conversion determinant to obtain a magnitude of a component force acting on the object to be measured. 請求項３に記載の分力測定装置であって、
前記記憶装置が保存する関係式は、最小自乗法により求められた近似式であることを特徴とする分力測定装置。The component force measuring device according to claim 3,
The component force measuring device, wherein the relational expression stored in the storage device is an approximate expression obtained by a least square method.

Images