JP2013057691A - 荷重推定方法 - Google Patents

Abstract

【課題】部材に衝撃力を加えたときに、その部材の任意の位置にどの程度の荷重が加わったかを推定する荷重推定方法を提供する。
【解決手段】被計測部材１０の表面の任意の位置に歪ゲージ１〜６を複数貼り付け、被計測部材に実際に荷重を加えたときの歪ゲージにより計測された値と、あらかじめコンピュータシミュレーションにより被計測部材の形状モデルを用いて荷重を加えたときの歪ゲージを貼り付けた位置に相当する位置での歪と当該位置における荷重値との関係を表す関係式を求めておき、シミュレーションによって得られた関係式に歪ゲージにより計測された歪の値を代入して、被計測部材に実際に荷重を加えたときの歪ゲージを貼り付けた位置における荷重値を求めることを特徴とする荷重推定方法である。
【選択図】図２

Description

本発明は、荷重推定方法に関し、詳しくは、たとえば車両の衝突試験のように部材に荷重を加えた際にその部材内部にかかる荷重値を推定する荷重推定方法に関する。

車両衝突試験は、バリアに実際の車両を衝突させることにより行われる。従来、この試験において車両を構成する部材の変形量を知るための方法として、車両の部材に加速度計を取り付け、衝突試験時の加速度から加速度計を取り付けた部位の変位量を算出する方法がある。この方法では、加速度計を変形する部位と変形しない部位に取り付け、それら各部位の加速度から変形部位の加速度を補正することにより、各部位の変形量を算出することとなっている（特許文献１）。

特開２００２−０９０３８３号公報

しかしながら、上記従来の方法では、加速度計で得られた値を部材各部へ分配することになるが、実際に部材にかかる荷重値は計測できないため変形量が予定と違う場合などに原因が特定できない。たとえば、あらかじめＦＥＭ解析などで変形量を求めて実験と比較する場合、たとえば車両減速度や生存空間（ダッシュ後退量）の解析値が加速度計で求めた実験値と合わない時などに、真の原因解明ができないという問題がある。

このようなことは、たとえば、変形部位と変形しない部位の加速度（衝突試験においては減速度となる）の差が大きいため、つぶれが大きく生存空間が足らない事態が発生した場合などに、その生存空間を形成している部材が受け持つ質量配分が違っており、解析で予定した以上の荷重がその部材にかかりその部材が座屈して生存空間が足りなくなったのか、その部材の接合強度が不足しており、荷重が大きくないにもかかわらず変形部位と変形しない部位の加速度の差が大きく生存空間が足らなくなったのか、原因を見極める必要がある。

しかし、従来の荷重推定法では車体の部材内部の荷重を直接測定しているものではないため、生存空間不足などが生じても真の原因は知り得ないことになる。

また、ある条件で実験結果と衝突解析結果とが合っている場合に、条件を違えると合わなくなるような場合に、実験結果と衝突解析結果の条件の違いによる符合性を実車を使って何度も衝突試験を行い条件の違いによる結果の違いを見つけることは事実上不可能であるため、真の原因解明ができないという問題がある。これもひとえに、実験力学的な意味のある車体の部材内部における荷重値の実験データが取得できていない、という点に起因する問題である。

このように車体を構成する部材内部の荷重値は、部材にかかる現実の荷重値を測るべきものである。しかし、現状で部材にかかる現実の荷重値を計測するためには、車体を切断してロードセルを入れるしかない。しかし、このようにしてロードセルを入れると、車体の構造を変えることとなってしまい計測目的に照らして許されないことである。

そこで本発明の目的は、部材に衝撃力を加えたときに、その部材の任意の位置にどの程度の荷重が加わったかを、部材を切断するなどしなくても推定することのできる荷重推定方法を提供することである。

上記課題を解決するための本発明は、被計測部材に力が加えられたときに当該被計測部材の断面に加わった荷重を推定する荷重推定方法であって、前記被計測部材の表面の任意の位置に歪ゲージを複数貼り付け、当該被計測部材に実際に荷重を加えたときの前記歪ゲージにより計測された値と、あらかじめコンピュータシミュレーションにより前記被計測部材の形状モデルを用いて荷重を加えたときの歪ゲージを貼り付けた位置に相当する位置での歪と当該位置で切り取られる断面における荷重値との関係を表す関係式を求めておき、前記コンピュータシミュレーションによって得られた前記関係式に前記歪ゲージにより計測された歪の値を代入して、前記被計測部材に実際に荷重を加えたときの前記歪ゲージを貼り付けた位置における荷重値を求め、当該求めた荷重値を、前記被計測部材の表面に貼り付けた前記複数の歪みゲージの貼り付け位置で切り取られる断面に加わった荷重とすることを特徴とする荷重推定方法である。

本発明によれば、荷重が加えられて変形する部材の任意の位置にかかる荷重値を高精度で推定することができる。

本発明を適用した実施形態における歪ゲージを貼り付ける位置の一例を説明するための説明図である。 実施形態のコンピュータシミュレーションにおいてディスプレイに映し出されたモデルを示す図である。 実施形態の歪ゲージの位置ずらしを説明するための図である。 実施形態の試験例１において荷重測定を行った車両の骨格部材を説明するための車両ボディ構造概略を示す図面である。 試験例１においてロードセルによって得られた荷重測定値と本実施形態１による通過力（荷重値）推定値の時間変化を示す図面である。 歪と応力の関係を示した図である。 実施形態２の試験例２によるロードセルの測定値、歪ゲージを貼り付けた位置における荷重推定値を示す図面である。

以下、図面を参照して本発明を実施するための形態について説明する。

［実施形態１］
まず、本実施形態１による荷重推定方法の概要について説明する。

本実施形態１による荷重推定方法は、車両衝突試験において、衝突時に加わる荷重を知りたい部位に歪ゲージを貼り付けて衝突試験を行う一方、当該部位にかかる応力と荷重の関係式をあらかじめ有限要素法（ＦＥＭ）によるコンピュータシミュレーション（以下、ＦＥＭ解析と称する）によって求めておく。そして、歪ゲージの値をＦＥＭ解析によって得られた関係式（行列式）に当てはめて、衝突時に部材の歪ゲージを貼り付けた位置にかかった荷重を推定するものである。

歪ゲージを貼り付ける位置は、車両内の部材の任意の位置でよい。部材としてはたとえば車両を構成するための充当部材である骨格部材などである。

図１は、歪ゲージを貼り付ける位置の一例を説明するための説明図である。

歪ゲージを貼り付ける位置は、たとえば、図示するように、骨格部材１０の周囲面に、合わせて６点歪ゲージ１〜６を貼り付ける。基本的には部材の断面線１１（図示一点鎖線で示す）の長い面の対向する面に２枚ずつと、この面に直交する部材の断面線１１の短い対向面に１枚ずつの合計６枚である。図示するものでは、上面１０ａが部材の断面線が長いので、この上面１０ａと下面１０ｃにそれぞれ２つ（符号１、２、４および５）、両側面１０ｂおよび１０ｄが断面線が短いのでそれぞれ１つ歪ゲージ（符号３および６）を貼り付けている。歪ゲージを貼り付ける方向は、断面にかかる荷重、すなわち部材を通過する力を測定する場合は、部材に衝撃が加わる方向（ここでは部材の長手方向）の歪を計れるように貼り付ける。

一方、ＦＥＭ解析は、市販されているＦＥＭ解析用プログラムを用いればよく、物体（解析モデル）を変形させたときに、その物体にかかる荷重をメッシュの節点ごとに求めることのできるシミュレーションプログラムを用いればよい。たとえば、ＮＡＳＴＲＡＮなどの有限要素解析ソルバーを用いることができる。

ＦＥＭ解析においては、解析モデルとして歪ゲージを貼り付けた骨格部材の有限要素モデル（メッシュモデル）を設定して、歪ゲージを貼り付けた部分に相当する位置のメッシュ節点をエレメントとする。そして、この解析モデル全体に荷重（力）を加えて、各エレメントに生じる応力と荷重値を求める。図１に示したように６枚の歪ゲージを貼り付けた場合、エレメントの数も６個となる。

ＦＥＭ解析条件は、モデル形状、節点数（またはメッシュ間隔）、加える荷重の値などを設定して解析すればよい。具体的には、モデル形状は解析を実行する部材形状に合わせればよい。節点数（またはメッシュ間隔）は、たとえば荷重値を知りたい断面位置と部材の大きさにより決めればよい。たとえば車体骨格部材の場合、５ｍｍ間隔程度のメッシュとなるように設定すればよい。部材全体に加える荷重（力）は、シミュレーションによって部材内部の各エレメントにおける応力とその位置の荷重値がわかればよいので、たとえば１Ｎ（ニュートン）など分かりやすい値を設定すればよい。部材全体に加える荷重（力）は、衝突試験を想定して、部材のある一点から力が加わるようにシミュレーションする。

図２は、コンピュータシミュレーションにおけるディスプレイに映し出されたモデルを示す図である。

ＦＥＭ解析では、図示するようなメッシュモデルが用いられ、このモデルに、図示（力）Ｆで示される位置、方向で荷重が加えられ、荷重値を求めたい断面Ｄにおけるエレメントに相当する部分のメッシュ節点の応力および荷重値を求めることになる。

衝突試験では、６枚の歪ゲージから得られた衝突試験時の骨格部材の歪の値εｍｘとする。一方、ＦＥＭ解析から算出した荷重値のうち、歪ゲージを貼り付けた位置で切り取られる断面での値としてその部分のエレメントの値を抽出し、これを荷重値を求めたい断面として、Ｈ行列（以下、行列を［Ｈ］とあらわす。下記式参照）とする。このＨ行列を特異値分解して得た直交基底（行列［Ｕ］）およびその逆行列（行列［Ｕ］^＃）と、ＦＥＭ解析から求めた断面上の全エレメントの断面積を表す行ベクトル｛ｂ｝^Ｔと、全エレメント（ｎ個）から計測に使用するｍ個（ｍ＝６程度）のエレメントを選択する行列［Ｓ］とを用いて、下記の式で車体骨格部材の通過力Ｔｘ（その断面にかかる荷重値）を推定する。ここで車体骨格部材の通過力は、図２に示すように、ある断面Ｄに生じている合力である。

骨格通過力Ｔｘ＝｛ｂ｝^Ｔ＊Ｅ＊［Ｕ］［Ｓ］［Ｕ］^＃＊｛εｘｍ｝ …（１）
なお、Ｅはヤング率である。

この式を具体的に行列式として示すと、（２）式のとおりである。

この（２）式の応力σ^ｉ _ｘは、シミュレーションにより求めた歪ゲージを貼り付けた部位におけるエレメントの応力の値である。また、ここでのＴｘはシミュレーションによって得られた歪ゲージを貼り付けた部位の断面にかかる荷重値である。

ここで応力と歪ゲージで計測される歪の関係は（３）式のとおりである。

また、歪は曲げ方向との関係で（４）式のように表すことができる。

この（５）式に対し、ｂ１，ｂ２，…ｂｎにＦＥＭ解析モデルにおける断面上の全エレメントの断面積を表す行ベクトル、ε^ｉ _ｘｉｎ，ε^ｉ _ｘｏｕｔに衝突試験で得られた各歪ゲージの値を入れることで、所望する断面（図１で一点差線で示す歪ゲージを貼り付けた部分の断面）の骨格通過力Ｔｘ、すなわち、その断面に加わった荷重値を得ることができる。なお、歪ゲージの値は内側か外側かの曲げ方向を考慮して使用する。

このようにＦＥＭ解析結果と歪ゲージの値から、車両を構成する部材の任意の位置における断面に加わる荷重値を推定することができるようになる。

次に、この方法において、最適な歪ゲージの貼り付け位置について説明する。

まず、部材に加わる荷重をその方向により分解すると、Ｆｘ、Ｆｙ、Ｆｚとなる。また、回転方向の力としてモーメントがあり、これを回転方向で分割すると、Ｍｘｙ、Ｍｙｚ、Ｍｚｘとなる。これら荷重のかかる方向とモーメントの方向を合わせて６分力と称する。

そこで、これら分解した６分力を種々変化させてＦＥＭ解析を行い、下記（６）式を用いて、行列［Ｈ］を決定する。ここで、Ｆｘ，Ｆｙ，Ｆｚ，Ｍｘｙ，Ｍｙｚ，Ｍｚｘは単位荷重とする。たとえば，１Ｎ（ニュートン）など分かりやすい値を設定すればよい。そして歪分布εｘを正規直交基底［Ｕ］を使って推定する。

特異値分解は（７）式のとおりである。

この特異値分解の式を用いて、下記（８）〜（１２）式のとおり、骨格通過力Ｔｘを求める式を完成される。

（１２）式のとおり骨格通過力Ｔｘを求める式ができる。

ここで、ＦＥＭ解析モデル上で、少しずつ歪ゲージの位置をずらし、その位置の違いによる可能な全ての組み合わせで歪ゲージを選択する行列［Ｓ］を作成する。

そしてＦＥＭ解析モデル上の付与荷重値と同定された骨格通過力Ｔｘの残差を全ての［Ｓ］について計算する。残差の最も小さくなる［Ｓ］が表すエレメントの組み合わせが最も同定精度の良い歪ゲージ配置である。

すなわち、歪ゲージ貼り付け位置をさまざまに変えてコンピュータシミュレーションによって歪と荷重との関係を求めて、あらかじめコンピュータシミュレーションにより被計測部材の形状モデルを用いて荷重を加えたときの歪ゲージを貼り付けた位置に相当する位置での歪に誤差を加えた量の歪を、コンピュータシミュレーションによって得られた歪と荷重の関係式である、上記（１２）式に代入してコンピュータシミュレーションによる誤差のある荷重を求め、コンピュータシミュレーションによる歪ゲージを貼り付けた位置に相当する位置での荷重との誤差が最小となる前記歪ゲージ貼り付け位置を実際の部材に対する歪ゲージ貼り付け位置とするのである。

ここで少しずつ歪ゲージの位置をずらすとは、たとえば、図３に示すように、解析モデル上において歪ゲージを貼り付ける位置として、たとえば１メッシュごとに、歪ゲージ１、２、４、および５をその面に沿うｚ方向にずらし、歪ゲージ３および６をその面に沿うｙ方向にずらすことである。なお、図３は解析モデルを示す図でありがメッシュは省略した。

ｎ個の全エレメントからｍ個の計測エレメントを選択する行列［Ｓ］を設定するにあたり、可能な全ての組み合わせの［Ｓ］行列で荷重Ｆｘ同定値の残差を計算し、残差の最も小さな計測エレメントの組み合わせで骨格部材の歪ε_ｘｍを計測する。

このように、あらかじめ部材の荷重値を得たい目的とする位置における歪と力との関係ＦＥＭ解析により求めておいて、歪ゲージ貼り付け位置をさまざまに変えて同様にＦＥＭ解析を行った結果と対比し、それらの差が最も少なくなる位置を実際の部材に対する歪ゲージ貼り付け位置とすることで、最も良好な精度で車体骨格通過力Ｔｘを計測することができる。

（試験例１）
実際に車両の側面衝突試験を行った際の車両の骨格部材における荷重実測値と、本実施形態に基づいた荷重推定結果について説明する。

図４は、荷重測定を行った車両の骨格部材を説明するための車両ボディ構造概略を示す図面である。

この試験では、図４に示した車両内における横方向の骨格部材１０において側面衝突試験を行った際の荷重を実測した。なお、図４は、あくまでも測定した骨格部材のボディ内での概略位置を説明するためのものであり、実際に側面衝突試験に使用した車両そのものを表すものではない。

荷重値の実測は、この骨格部材１０の中程に、骨格部材１０の強度ができるだけ低下しないように注意して切り込みを付けて、そこにロードセルを設置して測定した。ロードセルの設置位置は、骨格部材の端（その部材に衝撃が加わる点）から４４０ｍｍの位置に設置した。

そして、この骨格部材のロードセル設置位置から衝撃が加わる点から遠い側に６０ｍｍ離れた位置に本実施形態に基づき歪ゲージを６枚貼り付けた。

通過力（荷重値）Ｔｘを求めるための行列式（上記（５）式）および歪ゲージ貼り付け位置はあらかじめＦＥＭ解析により求めておいた。

このように設定された車両の側面衝突試験を行った。このときロードセルの値および歪ゲージの値はリアルタイムで計測した。

図５は、ロードセルによって得られた荷重測定値と本実施形態１による通過力（荷重値）推定値の時間変化を示す図面である。図５においては、ロードセルによる荷重測定値、歪ゲージを６枚使用して推定した場合の推定値、歪ゲージを４枚さまざまな組み合わせで使用して推定した場合の推定値について示した。図において、各グラフを指示している「ロードセル」ロードセルによる荷重測定値を示し、「１，２，３，４，５，６」、「２，３，４，６」、「１，３，４，６」、「２，３，５，６」、「１，３，５，６」の数字は図１に示した符号の歪ゲージに対応した歪ゲージを選択して推定した荷重値に対応している。

図からわかるように、歪ゲージを６枚使用したものは、ロードセルで計測された荷重値からわずかに低い値を示している。これら２つの値が違うのはロードセルの設置位置と歪ゲージ貼り付け位置に差があるためと考えられる。したがって、本実施形態１による荷重値の推定精度が非常に高いことがわかる。

また、荷重値推定に使用する歪ゲージの選択を変えた場合、６枚使用した場合が最も良い値となっていることがわかる。

次に、部材に生じる歪が弾性限界を越えている場合について説明する。

部材に貼り付けた歪ゲージによって計測される歪の値（全歪）εは、塑性歪相当分ε^（ｐ）と弾性歪相当分ε^（ｅ）の合計値である。すなわち、ε＝ε^（ｐ）＋ε^（ｅ）となる。

歪と応力の関係は、上述した（３）式のとおりである。図６は、この歪と応力の関係を示した図である。

上記のとおり全歪εは塑性歪相当分ε^（ｐ）と弾性歪相当分ε^（ｅ）での合計値であるから全歪εから塑性歪相当分ε^（ｐ）を引けば弾性歪相当分ε^（ｅ）が得られる。

そこで、衝突試験においては、衝突終了後衝突荷重が０（ゼロ）となった時点での歪ゲージの値を塑性歪相当分ε^（ｐ）と仮定すれば、最終的に得られた歪ゲージの値からこの衝突荷重が０（ゼロ）となった時点での歪ゲージの値を引けば、弾性歪相当分ε^（ｅ）の値が得られる。衝突荷重が０（ゼロ）となった時点は、たとえば、衝突試験において、車両をバリアなどに衝突させた後０．１ｍｓｅｃなど、衝撃が加わった時点から極短い時間を規定して、その時点の歪ゲージの値を衝突荷重が０（ゼロ）となった時点での歪ゲージの値とすればよい。

このようにして得られる弾性歪ε^（ｅ）の値を前述した（５）式におけるε^ｉ _ｘｉｎ，ε^ｉ _ｘｏｕｔの値として用いれば、弾性限界を超える骨格部材の任意の位置における荷重値（通過力）を求めることができる。下記（１３）式は、（５）式におけるε^ｉ _ｘｉｎ，ε^ｉ _ｘｏｕｔを弾性歪ε^（ｅ）に変えた式である。

［実施形態２］
本実施形態２は、荷重値推定に使用する行列式をより簡単に求めるものである。

下記（１４）式は、部材のある断面での６分力Ｆｘ，Ｆｙ，Ｆｚ，Ｍｘｙ，Ｍｘｚ，Ｍｙｚと歪εｚとの関係を示す式である。

εｚ＝α１Ｆｘ＋α２Ｆｙ＋α３Ｆｚ＋α４Ｍｘｙ＋α５Ｍｘｚ＋α６Ｍｙｚ …（１４）
骨格部材の任意の位置における荷重値を求めるために、図１に示したと同様に６枚の歪ゲージを貼り付けた場合、各歪ゲージで測定される歪εとの関係を（１４）式をもとに行列式に展開すれば、（１５）式のとおりとなる。

この（１５）式の行列式のαの行列部分を［Ａ］行列と称する。この（１５）式を簡略化して書き表せば、｛εｉ｝＝［Ａ］｛ｆ｝となる（ただし、ｉはひずみゲージの数であり、６枚の歪ゲージを用いた場合１〜６となる。以下同様）。

本実施形態２では、この（１５）式の行列式の［Ａ］行列部分をＦＥＭ解析によって求める。ＦＥＭ解析は実施形態１と同様に、荷重値を求めたい骨格部材の解析モデルを作成する。そして、部材に力からを加えたときに、歪ゲージを貼り付けた位置における６分力と応力を求める。応力と歪の関係は、前述した（３）式のとおりである。

ＦＥＭ解析によって上記（１５）式の行列式の［Ａ］行列が求まるので、この（１５）式から得られた［Ａ］行列の逆行列［Ａ］^♯を使用して、骨格部材の任意の位置における荷重値（通過力）を求める。このための逆行列式は、下記（１６）式のとおりとなる。

｛ｆ｝＝［Ａ］^♯｛εｉ｝ …（１６）
そして、所望する荷重値を求めたい骨格部材の任意の位置に歪ゲージを貼り付けて、衝突試験を行い、得られた各歪ゲージの値を上記（１６）式の｛εｉ｝に入れれば、６分力の行列｛ｆ｝が得られる。そして、通過力は骨格部材の長手方向へ通る力であるから６分力の行列｛ｆ｝のうちＦｘの値が、歪ゲージを貼り付けた位置における荷重値（通過力）となる。

（試験例２）
上述した実施形態１における試験例１の車両側面衝突試験を行ったときの歪ゲージの値を用いて、本実施形態２による推定方法を用いて骨格部材１０の任意の位置における荷重値を求めた。なお、ロードセルによる荷重値の実測値は実施形態１の試験例１と同じである。したがって、ロードセルの設置位置は骨格部材の端（その部材に衝撃が加わる点）から４４０ｍｍの位置に設置した。

歪ゲージは、この骨格部材１０のロードセル設置位置から衝撃が加わる点（衝撃点）に近い側に６０ｍｍ離れた位置に６枚、ロードセル設置位置から衝撃点から遠い側に６０ｍｍ離れた位置に６枚、それぞれ貼り付けた。ロードセルの値、および歪ゲージの値はリアルタイムに測定した。

図７に、ロードセルの測定値、本実施形態２に基づいて推定された歪ゲージを貼り付けたい値における荷重値（通過力）を示す。

図７に示すように、衝撃点に近い側における荷重推定値がロードセルの計測値より高く、逆に、衝撃点から遠い側における荷重推定値がロードセルの計測値より低い。したがって、衝撃が加わった部材を通った力が本実施形態２によって正しく推定されていることを示している。

以上説明した実施形態１および２によれば、シミュレーション（ＦＥＭ解析）によって得られる部材内部に加わる力と応力との関係と、実際の衝突試験の際に部材に貼り付けた歪ゲージの値から部材内部を通る荷重値を求めることとしたので、部材内部の任意の位置における荷重値を精度よく推定することが可能となる。しかも、歪ゲージは部材の外側に貼り付けるだけでよいため、これらから衝突試験を行う部材に切れ込みを入れてロードセルを設置するなどといった操作が不要であり、部材強度になんら変化を与えることなく、部材内部にかかる荷重を任意の位置で知ることができる。

また、荷重値を得たい位置とその位置で最も誤差が少なくなる歪ゲージ貼り付け位置をＦＥＭ解析により求めたので、計測断面全周に歪ゲージを貼り付けることなく、わずかに６枚の歪ゲージで精度のよい荷重値の推定を行うことができる。このため歪ゲージの貼り付け工数も大幅に省略することができる。

また、歪ゲージで測定された値を弾性歪相当分と塑性歪相当分に分け、弾性歪相当分を用いることとすれば、部位の応力が弾性限界を越えている場合でも、歪ゲージのデータを用いて、高精度で部材の任意の位置にかかる荷重値を求めることができる。

以上本発明を適用した実施形態について説明したが、本発明は、これらの実施形態に限定されるものではない。上述した実施形態では車両衝突試験において、車両の骨格部材にかかる荷重を測定対象としているが、本発明は車両の骨格部材に限らず、さまざまな剛体（部材）について実施可能である。たとえば、車両のホイールなどについても同じように衝撃が加わった際の任意の位置における荷重値を推定することができる。そのほか、車両にかかわりのないような部材であっても同様に衝撃が加わった際に、その部材内部における任意の位置で荷重を推定することができる。

１、２、３、４、５、６ 歪ゲージ、
１０ 骨格部材。

Claims (5)

1. 被計測部材に力が加えられたときに当該被計測部材の断面に加わった荷重を推定する荷重推定方法であって、
   前記被計測部材の表面の任意の位置に歪ゲージを複数貼り付け、当該被計測部材に実際に荷重を加えたときの前記歪ゲージにより計測された値と、
   あらかじめコンピュータシミュレーションにより前記被計測部材の形状モデルを用いて荷重を加えたときの歪ゲージを貼り付けた位置に相当する位置での歪と当該位置で切り取られる断面における荷重値との関係を表す関係式を求めておき、
   前記コンピュータシミュレーションによって得られた前記関係式に前記歪ゲージにより計測された歪の値を代入して、前記被計測部材に実際に荷重を加えたときの前記歪ゲージを貼り付けた位置における荷重値を求め、
   当該求めた荷重値を、前記被計測部材の表面に貼り付けた前記複数の歪みゲージの貼り付け位置で切り取られる断面に加わった荷重とすることを特徴とする荷重推定方法。
2. 前記関係式は、有限要素法によるコンピュータシミュレーションによって、前記被計測部材のある一点から力が加えられたときに、歪みゲージ貼り付け位置で切り取られる前記被計測部材の断面におけるメッシュ節点のうち歪みゲージ貼り付け位置における節点の応力と荷重の関係から前記断面に加わる荷重値を推定する関係式として求めておくものであり、
   前記歪ゲージを貼り付けた位置における荷重値は、当該関係式中の応力を歪みに置換したうえで、当該関係式に歪みゲージによる計測値を代入して被計測部材断面に加わる荷重値として推定することを特徴とする請求項１記載の荷重推定方法。
3. 前記歪ゲージの貼り付け位置は、
   前記被計測部材の荷重値を得たい目的とする位置における前記歪と前記荷重との関係式をコンピュータシミュレーションにより求め、
   前記歪ゲージ貼り付け位置をさまざまに変えて同様にコンピュータシミュレーションによって前記歪と前記荷重との関係を求めて、あらかじめコンピュータシミュレーションにより前記被計測部材の形状モデルを用いて荷重を加えたときの歪ゲージを貼り付けた位置に相当する位置での歪に誤差を加えた量の歪を、前記コンピュータシミュレーションによって得られた前記歪と荷重の関係式に代入してコンピュータシミュレーションによる誤差のある荷重を求め、前記コンピュータシミュレーションによる歪ゲージを貼り付けた位置に相当する位置での荷重との誤差が最小となる前記歪ゲージ貼り付け位置を実際の部材に対する前記歪ゲージ貼り付け位置とすることを特徴とする請求項１または２記載の荷重推定方法。
4. 前記歪ゲージにより計測された値を、弾性歪相当分と塑性歪相当分に分け、当該弾性歪相当分を前記コンピュータシミュレーションによって得られた歪と荷重の関係式に代入して、前記被計測部材に実際に荷重を加えたときの前記歪ゲージを貼り付けた位置における荷重値を求める、ことを特徴とする請求項１〜３のいずれか一つに記載の荷重推定方法。
5. 前記歪ゲージは前記被計測部材の荷重値を得たい位置の周囲に６枚、それぞれの歪みゲージが部材に力が加わる方向の歪を測れるように貼り付けることを特徴とする請求項１〜４のいずれか一つに記載の荷重推定方法。