JP3662682B2 - 航空機主翼に作用する荷重の推定方法 - Google Patents
航空機主翼に作用する荷重の推定方法 Download PDFInfo
- Publication number
- JP3662682B2 JP3662682B2 JP24983896A JP24983896A JP3662682B2 JP 3662682 B2 JP3662682 B2 JP 3662682B2 JP 24983896 A JP24983896 A JP 24983896A JP 24983896 A JP24983896 A JP 24983896A JP 3662682 B2 JP3662682 B2 JP 3662682B2
- Authority
- JP
- Japan
- Prior art keywords
- strain
- load
- deflection
- stage
- main wing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0041—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress
- G01M5/005—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress by means of external apparatus, e.g. test benches or portable test systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
- G01B7/18—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge using change in resistance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/20—Measuring force or stress, in general by measuring variations in ohmic resistance of solid materials or of electrically-conductive fluids; by making use of electrokinetic cells, i.e. liquid-containing cells wherein an electrical potential is produced or varied upon the application of stress
- G01L1/22—Measuring force or stress, in general by measuring variations in ohmic resistance of solid materials or of electrically-conductive fluids; by making use of electrokinetic cells, i.e. liquid-containing cells wherein an electrical potential is produced or varied upon the application of stress using resistance strain gauges
- G01L1/2206—Special supports with preselected places to mount the resistance strain gauges; Mounting of supports
- G01L1/2218—Special supports with preselected places to mount the resistance strain gauges; Mounting of supports the supports being of the column type, e.g. cylindric, adapted for measuring a force along a single direction
- G01L1/2225—Special supports with preselected places to mount the resistance strain gauges; Mounting of supports the supports being of the column type, e.g. cylindric, adapted for measuring a force along a single direction the direction being perpendicular to the central axis
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Measurement Of Force In General (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は歪又はたわみを計測して、この測定値から航空機主翼に作用する荷重を推定する技術に関する。
【０００２】
【従来の技術】
歪測定手段（例えば、歪ゲージ）を構造物に貼りつけ、その点における歪を直接計測し、この歪データから応力を求め、これらから作用荷重を推定する方法がある。
例えば、単純形状である棒に軸力を作用させたときに次の式が成立する。
Ｅ＝σ／ε＝（Ｐ／Ａ）／（λ／ｌ）
ただし、Ｅは縦弾性係数、σは単純な垂直応力、εは縦歪、Ｐは軸荷重、Ａは横断面積、λは部材の伸び又は縮み、ｌは部材の原長である。
【０００３】
【発明が解決しようとする課題】
しかし、実際の構造物は、断面が一様でなく、荷重の作用方向が複雑であり且つ作用点が複数であることから、歪は実測できるが、この歪から作用荷重を換算することは不可能に近い。
そこで、歪を測ることを断念して、コンピュータ・シミュレーションを用いると構造物のモデル化の仕方によって計算結果が著しく異なったり、荷重条件の仮定が難しいなどの理由から、作用荷重を計算により推定することは信頼性が乏しい。そして、コンピュータ・シミュレーションは設計段階で荷重を推定する手法であって、構造物を現実に使用しているときの作用荷重を求めるには限界がある。
【０００４】
例えば、航空機の主翼の各部に、実際の飛行状態において、いくらの荷重が作用しているかを知ることができれば、安全管理上有益であるが、この技術は極めて難しい。
すなわち、主翼の各部に歪ゲージを多数貼り付けたとしても、それらの点における歪が分かるだけで、これらの歪が許容値以内にあるか否かが分かるだけであって、機体の特定部位にいくらの荷重が作用しているのかを知ることは困難である。
【０００５】
【課題を解決するための手段】
しかし、本発明者は、研究を進めた結果、複数の歪ゲージ又はたわみゲージを巧みに利用して作用荷重を精度よく求めることのできる技術を開発することに成功した。
即ち、請求項１は、航空機主翼上のｎ個の点に各々作用する荷重を推定する方法を次に述べる第１段階と第２段階で実施するものであり、
第１段階は、ｎ個の箇所に各々歪測定手段又はたわみ測定手段を設ける工程と、前記ｎ個の箇所のうちの１箇所に基準荷重を集中的に掛け、そのときに１番からｎ番までに発生する歪又はたわみを測定し、得られた測定値をその時の基準荷重で割って単位荷重当りの歪又はたわみにする工程と、この工程を順次繰り返して下記の実用式を作成する工程とからなり、
【０００６】
【数２】
【０００７】
第２段階は、航空機主翼を実使用状態に保持して歪測定手段又はたわみ測定手段によりｎ個の歪又はたわみ生データを集録する工程と、これらの歪又はたわみ生データを前記実用式に代入して推定荷重を計算する。
【０００８】
例えば、第１段階を実験場で実施し、第２段階は航空機主翼を使用状態に維持して実施する。
第２段階では歪又はたわみ生データを集録し、簡単な計算のみで各点に作用する荷重が直接推定できる。
従来は高速な電子計算機にて膨大な量の繰り返し計算を行なっても、比較的精度の低い推定値しか得られなかった。これに対して、本発明によれば多元連立方程式を解くだけで精度の高い推定値が得られ、高速な電子計算機は不必要となった。
【０００９】
【発明の実施の形態】
本発明の実施の形態を添付図に基づいて以下に説明する。
図１（ａ）〜（ｃ）は本発明方法を導くためのサンプルモデル図であり、梁状の構造物１に、基準荷重Ｆ1〜Ｆ3を適当な間隔で、個別に作用させたときに発生する各部における歪を考える。前記基準荷重Ｆ1〜Ｆ3は集中荷重である。
（ａ）は、１個の基準荷重Ｆ1のみを掛けたときの状態を示し、このときに各部に発生する歪を順に、ε11，ε12及びε13とする。
（ｂ）は、１個の基準荷重Ｆ2のみを掛けたときの状態を示し、このときに各部に発生する歪を順に、ε21，ε22及びε23とする。
（ｃ）は、１個の基準荷重Ｆ3のみを掛けたときの状態を示し、このときに各部に発生する歪を順に、ε31，ε32及びε33とする。
【００１０】
前記歪は例えば次に示す要領で計測する。
図２は本発明方法の第１段階に係る装置の原理図であり、構造物１の荷重作用点に歪測定手段としての歪ゲージ２ａ，２ｂ，２ｃを貼り付け、ケーブル３・・・（・・・は複数本を示す。）でアンプ４に結合し、このアンプ４をＡ／Ｄ変換機構を備えた記録計５に繋き、更に記録計５を計算機６に繋ぐ。
上記構造物１は、例えば航空機の主要部材（主翼、胴体、尾翼など）や車両の主要部材（車体フレーム）である。
実験では１０ｍｍ×１０ｍｍの正方形断面のアルミニウム梁を用い、左の支持点から１／４スパン（２００ｍｍ）毎に荷重点を設定した。
【００１１】
歪ゲージ２ａ，２ｂ，２ｃは、周知のホイストンブリッヂや圧電素子などであり、歪量に応じた電気信号を発するものである。
実験場において、構造物１に、既知の１個の荷重Ｆ1を掛けて歪ε11，ε12及びε13を実測し、同様に１個の荷重Ｆ2を掛けて歪ε21，ε22及びε23を、１個の荷重Ｆ3を掛けて歪ε31，ε32及びε33を実測する。
【００１２】
表１は基準荷重Ｆ1を掛けたときの歪ε11，ε12及びε13を実測したものである。Ｆ1は８００ｇ、１６００ｇ、２１００ｇの３種類とした。
【００１３】
【表１】
【００１４】
（ａ）に示す通り、ε11については１１２×１０-6，２２６×１０-6，３０４×１０-6の３個の測定値を得た。
（ｂ）横軸をＦ1、縦軸をεとしたグラフに次の要領で測定値をプロットする。（８００，１１２×１０-6）、（１６００，２２６×１０-6）、（２１００，３０４×１０-6）の３点をプロットする。これらの３点と原点（０，０）との４点で最小二乗法にて近似式ε＝０．１４３２×１０-6・Ｆ1を得た。近似式の傾きは単位荷重当りの歪に相当する。
（ｃ）前記近似式の傾き、即ち単位荷重当りの歪をκ11とする。
同様に（ａ）のε12の３個の測定値から、（ｂ）のグラフを経て単位荷重当りの歪κ12を求め、更に、（ａ）のε13の３個の測定値から、（ｂ）のグラフを経て単位荷重当りの歪κ13を求め、（ｃ）にまとめる。
【００１５】
同様にＦ2，ε21，ε22及びε23について処理したものを表２、Ｆ3，ε31，ε32及びε33について処理したものを表３に示す。
【００１６】
【表２】
【００１７】
【表３】
【００１８】
次に示す数３は本発明の連立方程式を３元式にしたものである。
【００１９】
【数３】
【００２０】
この連立方程式の▲１▼式に表１（ｃ）に示したκ11，κ12，κ13の値を代入し、▲２▼式に表２に示したκ21，κ22及びκ23の値を代入し、▲３▼式に表３に示したκ31，κ32，κ33の値を代入する。ただし、両辺に１０6を乗じてκ11〜κ33における１０-6は消去する。
この結果が次の数４である。
【００２１】
【数４】
【００２２】
この連立方程式が、このときの本発明における「実用式」である。
この実用式の使用方法を次に説明する。
図３は本発明方法の第２段階に係る装置の原理図であり、この原理図の構成は前記図２と同じである。
歪ゲージ２ａ〜２ｃの箇所に実際に作用荷重Ａ1＝８００ｇ、Ａ2＝２１００ｇ、Ａ3＝１５００ｇを加える。
このときの各作用点における歪を計測したところ以下の通りであった。歪ゲージ２ａで３８６×１０-6、歪ゲージ２ｂで６２７×１０-6、歪ゲージ２ｃで４５６×１０-6の歪を計測した。これらは歪生データε1〜ε3である。これらを前記数４の右辺に代入したものが次の数５である。
【００２３】
【数５】
【００２４】
この連立方程式を解くと、Ａ1は８１９ｇ、Ａ2は２０４７ｇそしてＡ3は１５４１ｇである。
即ち、構造物１に発生した歪ε1〜ε3を実測することにより、各点に作用する荷重Ａ1〜Ａ3が推定できたわけである。
【００２５】
この推定値の精度を検討する。
推定値と実際に作用した荷重（Ａ1＝８００ｇ、Ａ2＝２１００ｇ、Ａ3＝２１００ｇ）とを比較して、推定精度を計算すると、Ａ1は２．４％、Ａ2は２．４％、Ａ3は２．７％であり、荷重の推定精度は良好であることが確認できた。
【００２６】
本発明の実施例を整理すると次の通りである。
図４は第１段階における計算フロー図であり、ＳＴ××はステップ番号を示す。なお、第１段階は設備の整った実験場で実施すると良い。
ＳＴ０１；構造物にｎ個の歪ゲージを取付ける。
ＳＴ０２；まず、１番歪ゲージの箇所に基準荷重（例えば８００ｇ）を集中的に掛ける。
ＳＴ０３；１番〜ｎ番までの各点に発生した歪を測定し、これらの歪を前記基準荷重で割って、単位荷重当りの歪とする。
ＳＴ０４；次に、２番歪ゲージの箇所に基準荷重（例えば８００ｇ）を集中的に掛ける。
ＳＴ０５；１番〜ｎ番までの各点に発生した歪を測定し、これらの歪を前記基準荷重で割って、単位荷重当りの歪とする。
【００２７】
ＳＴ０６；以上を繰り返すことにより、ｎ番歪ゲージの箇所に基準荷重（例えば８００ｇ）を集中的に掛ける。
ＳＴ０７；１番〜ｎ番までの各点に発生した歪を測定し、これらの歪を前記基準荷重で割って、単位荷重当りの歪とする。
ＳＴ０８；実用式（ｎ元連立方程式）を記憶する。
【００２８】
図５は第２段階における計算フロー図であり、ＳＴ××はステップ番号を示す。なお、第２段階は飛行中若しくは走行中に実施する。
ＳＴ１０；構造物に取付けたｎ個の歪ゲージでｎ個の歪データを集録する。
ＳＴ１１；歪データを実用式（ｎ元連立方程式）に代入して、Ａ1，Ａ2，・・・ＡNを算出する。これらＡ1，Ａ2，・・・ＡNが１番〜ｎ番に作用したと推定できる推定荷重である。
【００２９】
以上は歪で説明したが、歪をたわみに置き換えることで同様に荷重を推定することができる。たわみ測定手段としては、ダイヤルゲージなどの接触型ゲージや、空気マイクロメータ，オプチメータなどの非接触型ゲージが適当である。
【００３０】
また、「得られた測定値をその時の基準荷重で割って単位荷重当りの歪又はたわみにする」処理は、上記説明では最小二乗法で実施したが、これに準ずる手法であれば処理法は任意である。
【００３１】
尚、本実施の形態では、簡単な両端単純支持はりを例に説明したが、本発明方法は多点動的荷重の計測が可能であり、具体的には、航空機の飛行試験時の翼荷重の計測、海洋構造物、船体等に作用する流体力の計測が可能であり、このことから、本発明方法は、静的荷重、振動，衝撃等の動的荷重の計測に応用でき、センサや演算器を搭載することで、リアルタイムに荷重を推定して、過負荷信号を発信させることもできる。従って、本発明の適用分野は多岐にわたる。
【００３２】
【発明の効果】
本発明は上記構成により次の効果を発揮する。
請求項１は、航空機主翼上のｎ個の点に各々作用する荷重を推定する方法を次に述べる第１段階と第２段階で実施するものであり、
第１段階は、ｎ個の箇所に各々歪測定手段又はたわみ測定手段を設ける工程と、前記ｎ個の箇所のうちの１箇所に基準荷重を集中的に掛け、そのときに１番からｎ番までに発生する歪又はたわみを測定し、得られた測定値をその時の基準荷重で割って単位荷重当りの歪又はたわみにする工程と、この工程を順次繰り返して下記の実用式を作成する工程とからなり、
【００３３】
【数６】
【００３４】
第２段階は、航空機主翼を実使用状態に保持して歪測定手段又はたわみ測定手段によりｎ個の歪又はたわみ生データを集録する工程と、これらの歪又はたわみ生データを前記実用式に代入して推定荷重を計算する。
【００３５】
例えば、第１段階を実験場で実施し、第２段階は航空機主翼を使用状態に維持して実施する。
第２段階では歪又はたわみ生データを集録し、簡単な計算のみで各点に作用する荷重が直接推定できる。
従来は高速な電子計算機にて膨大な量の繰り返し計算を行なっても、比較的精度の低い推定値しか得られなかった。これに対して、本発明によれば多元連立方程式を解くだけで精度の高い推定値が得られ、高速な電子計算機は不必要となった。
【図面の簡単な説明】
【図１】本発明方法を導くためのサンプルモデル図
【図２】本発明方法の第１段階に係る装置の原理図
【図３】本発明方法の第２段階に係る装置の原理図
【図４】第１段階における計算フロー図
【図５】第２段階における計算フロー図
【符号の説明】
１…構造物、２ａ〜２ｄ…歪測定手段（歪ゲージ）、５…記録計、６…計算機。
【発明の属する技術分野】
本発明は歪又はたわみを計測して、この測定値から航空機主翼に作用する荷重を推定する技術に関する。
【０００２】
【従来の技術】
歪測定手段（例えば、歪ゲージ）を構造物に貼りつけ、その点における歪を直接計測し、この歪データから応力を求め、これらから作用荷重を推定する方法がある。
例えば、単純形状である棒に軸力を作用させたときに次の式が成立する。
Ｅ＝σ／ε＝（Ｐ／Ａ）／（λ／ｌ）
ただし、Ｅは縦弾性係数、σは単純な垂直応力、εは縦歪、Ｐは軸荷重、Ａは横断面積、λは部材の伸び又は縮み、ｌは部材の原長である。
【０００３】
【発明が解決しようとする課題】
しかし、実際の構造物は、断面が一様でなく、荷重の作用方向が複雑であり且つ作用点が複数であることから、歪は実測できるが、この歪から作用荷重を換算することは不可能に近い。
そこで、歪を測ることを断念して、コンピュータ・シミュレーションを用いると構造物のモデル化の仕方によって計算結果が著しく異なったり、荷重条件の仮定が難しいなどの理由から、作用荷重を計算により推定することは信頼性が乏しい。そして、コンピュータ・シミュレーションは設計段階で荷重を推定する手法であって、構造物を現実に使用しているときの作用荷重を求めるには限界がある。
【０００４】
例えば、航空機の主翼の各部に、実際の飛行状態において、いくらの荷重が作用しているかを知ることができれば、安全管理上有益であるが、この技術は極めて難しい。
すなわち、主翼の各部に歪ゲージを多数貼り付けたとしても、それらの点における歪が分かるだけで、これらの歪が許容値以内にあるか否かが分かるだけであって、機体の特定部位にいくらの荷重が作用しているのかを知ることは困難である。
【０００５】
【課題を解決するための手段】
しかし、本発明者は、研究を進めた結果、複数の歪ゲージ又はたわみゲージを巧みに利用して作用荷重を精度よく求めることのできる技術を開発することに成功した。
即ち、請求項１は、航空機主翼上のｎ個の点に各々作用する荷重を推定する方法を次に述べる第１段階と第２段階で実施するものであり、
第１段階は、ｎ個の箇所に各々歪測定手段又はたわみ測定手段を設ける工程と、前記ｎ個の箇所のうちの１箇所に基準荷重を集中的に掛け、そのときに１番からｎ番までに発生する歪又はたわみを測定し、得られた測定値をその時の基準荷重で割って単位荷重当りの歪又はたわみにする工程と、この工程を順次繰り返して下記の実用式を作成する工程とからなり、
【０００６】
【数２】
第２段階は、航空機主翼を実使用状態に保持して歪測定手段又はたわみ測定手段によりｎ個の歪又はたわみ生データを集録する工程と、これらの歪又はたわみ生データを前記実用式に代入して推定荷重を計算する。
【０００８】
例えば、第１段階を実験場で実施し、第２段階は航空機主翼を使用状態に維持して実施する。
第２段階では歪又はたわみ生データを集録し、簡単な計算のみで各点に作用する荷重が直接推定できる。
従来は高速な電子計算機にて膨大な量の繰り返し計算を行なっても、比較的精度の低い推定値しか得られなかった。これに対して、本発明によれば多元連立方程式を解くだけで精度の高い推定値が得られ、高速な電子計算機は不必要となった。
【０００９】
【発明の実施の形態】
本発明の実施の形態を添付図に基づいて以下に説明する。
図１（ａ）〜（ｃ）は本発明方法を導くためのサンプルモデル図であり、梁状の構造物１に、基準荷重Ｆ1〜Ｆ3を適当な間隔で、個別に作用させたときに発生する各部における歪を考える。前記基準荷重Ｆ1〜Ｆ3は集中荷重である。
（ａ）は、１個の基準荷重Ｆ1のみを掛けたときの状態を示し、このときに各部に発生する歪を順に、ε11，ε12及びε13とする。
（ｂ）は、１個の基準荷重Ｆ2のみを掛けたときの状態を示し、このときに各部に発生する歪を順に、ε21，ε22及びε23とする。
（ｃ）は、１個の基準荷重Ｆ3のみを掛けたときの状態を示し、このときに各部に発生する歪を順に、ε31，ε32及びε33とする。
【００１０】
前記歪は例えば次に示す要領で計測する。
図２は本発明方法の第１段階に係る装置の原理図であり、構造物１の荷重作用点に歪測定手段としての歪ゲージ２ａ，２ｂ，２ｃを貼り付け、ケーブル３・・・（・・・は複数本を示す。）でアンプ４に結合し、このアンプ４をＡ／Ｄ変換機構を備えた記録計５に繋き、更に記録計５を計算機６に繋ぐ。
上記構造物１は、例えば航空機の主要部材（主翼、胴体、尾翼など）や車両の主要部材（車体フレーム）である。
実験では１０ｍｍ×１０ｍｍの正方形断面のアルミニウム梁を用い、左の支持点から１／４スパン（２００ｍｍ）毎に荷重点を設定した。
【００１１】
歪ゲージ２ａ，２ｂ，２ｃは、周知のホイストンブリッヂや圧電素子などであり、歪量に応じた電気信号を発するものである。
実験場において、構造物１に、既知の１個の荷重Ｆ1を掛けて歪ε11，ε12及びε13を実測し、同様に１個の荷重Ｆ2を掛けて歪ε21，ε22及びε23を、１個の荷重Ｆ3を掛けて歪ε31，ε32及びε33を実測する。
【００１２】
表１は基準荷重Ｆ1を掛けたときの歪ε11，ε12及びε13を実測したものである。Ｆ1は８００ｇ、１６００ｇ、２１００ｇの３種類とした。
【００１３】
【表１】
（ａ）に示す通り、ε11については１１２×１０-6，２２６×１０-6，３０４×１０-6の３個の測定値を得た。
（ｂ）横軸をＦ1、縦軸をεとしたグラフに次の要領で測定値をプロットする。（８００，１１２×１０-6）、（１６００，２２６×１０-6）、（２１００，３０４×１０-6）の３点をプロットする。これらの３点と原点（０，０）との４点で最小二乗法にて近似式ε＝０．１４３２×１０-6・Ｆ1を得た。近似式の傾きは単位荷重当りの歪に相当する。
（ｃ）前記近似式の傾き、即ち単位荷重当りの歪をκ11とする。
同様に（ａ）のε12の３個の測定値から、（ｂ）のグラフを経て単位荷重当りの歪κ12を求め、更に、（ａ）のε13の３個の測定値から、（ｂ）のグラフを経て単位荷重当りの歪κ13を求め、（ｃ）にまとめる。
【００１５】
同様にＦ2，ε21，ε22及びε23について処理したものを表２、Ｆ3，ε31，ε32及びε33について処理したものを表３に示す。
【００１６】
【表２】
【表３】
次に示す数３は本発明の連立方程式を３元式にしたものである。
【００１９】
【数３】
この連立方程式の▲１▼式に表１（ｃ）に示したκ11，κ12，κ13の値を代入し、▲２▼式に表２に示したκ21，κ22及びκ23の値を代入し、▲３▼式に表３に示したκ31，κ32，κ33の値を代入する。ただし、両辺に１０6を乗じてκ11〜κ33における１０-6は消去する。
この結果が次の数４である。
【００２１】
【数４】
この連立方程式が、このときの本発明における「実用式」である。
この実用式の使用方法を次に説明する。
図３は本発明方法の第２段階に係る装置の原理図であり、この原理図の構成は前記図２と同じである。
歪ゲージ２ａ〜２ｃの箇所に実際に作用荷重Ａ1＝８００ｇ、Ａ2＝２１００ｇ、Ａ3＝１５００ｇを加える。
このときの各作用点における歪を計測したところ以下の通りであった。歪ゲージ２ａで３８６×１０-6、歪ゲージ２ｂで６２７×１０-6、歪ゲージ２ｃで４５６×１０-6の歪を計測した。これらは歪生データε1〜ε3である。これらを前記数４の右辺に代入したものが次の数５である。
【００２３】
【数５】
この連立方程式を解くと、Ａ1は８１９ｇ、Ａ2は２０４７ｇそしてＡ3は１５４１ｇである。
即ち、構造物１に発生した歪ε1〜ε3を実測することにより、各点に作用する荷重Ａ1〜Ａ3が推定できたわけである。
【００２５】
この推定値の精度を検討する。
推定値と実際に作用した荷重（Ａ1＝８００ｇ、Ａ2＝２１００ｇ、Ａ3＝２１００ｇ）とを比較して、推定精度を計算すると、Ａ1は２．４％、Ａ2は２．４％、Ａ3は２．７％であり、荷重の推定精度は良好であることが確認できた。
【００２６】
本発明の実施例を整理すると次の通りである。
図４は第１段階における計算フロー図であり、ＳＴ××はステップ番号を示す。なお、第１段階は設備の整った実験場で実施すると良い。
ＳＴ０１；構造物にｎ個の歪ゲージを取付ける。
ＳＴ０２；まず、１番歪ゲージの箇所に基準荷重（例えば８００ｇ）を集中的に掛ける。
ＳＴ０３；１番〜ｎ番までの各点に発生した歪を測定し、これらの歪を前記基準荷重で割って、単位荷重当りの歪とする。
ＳＴ０４；次に、２番歪ゲージの箇所に基準荷重（例えば８００ｇ）を集中的に掛ける。
ＳＴ０５；１番〜ｎ番までの各点に発生した歪を測定し、これらの歪を前記基準荷重で割って、単位荷重当りの歪とする。
【００２７】
ＳＴ０６；以上を繰り返すことにより、ｎ番歪ゲージの箇所に基準荷重（例えば８００ｇ）を集中的に掛ける。
ＳＴ０７；１番〜ｎ番までの各点に発生した歪を測定し、これらの歪を前記基準荷重で割って、単位荷重当りの歪とする。
ＳＴ０８；実用式（ｎ元連立方程式）を記憶する。
【００２８】
図５は第２段階における計算フロー図であり、ＳＴ××はステップ番号を示す。なお、第２段階は飛行中若しくは走行中に実施する。
ＳＴ１０；構造物に取付けたｎ個の歪ゲージでｎ個の歪データを集録する。
ＳＴ１１；歪データを実用式（ｎ元連立方程式）に代入して、Ａ1，Ａ2，・・・ＡNを算出する。これらＡ1，Ａ2，・・・ＡNが１番〜ｎ番に作用したと推定できる推定荷重である。
【００２９】
以上は歪で説明したが、歪をたわみに置き換えることで同様に荷重を推定することができる。たわみ測定手段としては、ダイヤルゲージなどの接触型ゲージや、空気マイクロメータ，オプチメータなどの非接触型ゲージが適当である。
【００３０】
また、「得られた測定値をその時の基準荷重で割って単位荷重当りの歪又はたわみにする」処理は、上記説明では最小二乗法で実施したが、これに準ずる手法であれば処理法は任意である。
【００３１】
尚、本実施の形態では、簡単な両端単純支持はりを例に説明したが、本発明方法は多点動的荷重の計測が可能であり、具体的には、航空機の飛行試験時の翼荷重の計測、海洋構造物、船体等に作用する流体力の計測が可能であり、このことから、本発明方法は、静的荷重、振動，衝撃等の動的荷重の計測に応用でき、センサや演算器を搭載することで、リアルタイムに荷重を推定して、過負荷信号を発信させることもできる。従って、本発明の適用分野は多岐にわたる。
【００３２】
【発明の効果】
本発明は上記構成により次の効果を発揮する。
請求項１は、航空機主翼上のｎ個の点に各々作用する荷重を推定する方法を次に述べる第１段階と第２段階で実施するものであり、
第１段階は、ｎ個の箇所に各々歪測定手段又はたわみ測定手段を設ける工程と、前記ｎ個の箇所のうちの１箇所に基準荷重を集中的に掛け、そのときに１番からｎ番までに発生する歪又はたわみを測定し、得られた測定値をその時の基準荷重で割って単位荷重当りの歪又はたわみにする工程と、この工程を順次繰り返して下記の実用式を作成する工程とからなり、
【００３３】
【数６】
第２段階は、航空機主翼を実使用状態に保持して歪測定手段又はたわみ測定手段によりｎ個の歪又はたわみ生データを集録する工程と、これらの歪又はたわみ生データを前記実用式に代入して推定荷重を計算する。
【００３５】
例えば、第１段階を実験場で実施し、第２段階は航空機主翼を使用状態に維持して実施する。
第２段階では歪又はたわみ生データを集録し、簡単な計算のみで各点に作用する荷重が直接推定できる。
従来は高速な電子計算機にて膨大な量の繰り返し計算を行なっても、比較的精度の低い推定値しか得られなかった。これに対して、本発明によれば多元連立方程式を解くだけで精度の高い推定値が得られ、高速な電子計算機は不必要となった。
【図面の簡単な説明】
【図１】本発明方法を導くためのサンプルモデル図
【図２】本発明方法の第１段階に係る装置の原理図
【図３】本発明方法の第２段階に係る装置の原理図
【図４】第１段階における計算フロー図
【図５】第２段階における計算フロー図
【符号の説明】
１…構造物、２ａ〜２ｄ…歪測定手段（歪ゲージ）、５…記録計、６…計算機。
Claims (1)
- 航空機主翼上のｎ個の点に各々作用する荷重を推定する方法を次に述べる第１段階と第２段階で実施するものであり、
第１段階は、ｎ個の箇所に各々歪測定手段又はたわみ測定手段を設ける工程と、前記ｎ個の箇所のうちの１箇所に基準荷重を集中的に掛け、そのときに１番からｎ番までに発生する歪又はたわみを測定し、得られた測定値をその時の基準荷重で割って単位荷重当りの歪又はたわみにする工程と、この工程を順次繰り返して下記の実用式を作成する工程とからなり、
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP24983896A JP3662682B2 (ja) | 1996-09-20 | 1996-09-20 | 航空機主翼に作用する荷重の推定方法 |
| US08/934,419 US5922967A (en) | 1996-09-20 | 1997-09-19 | Method and apparatus for estimating loads imposed on structural body |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP24983896A JP3662682B2 (ja) | 1996-09-20 | 1996-09-20 | 航空機主翼に作用する荷重の推定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH1090087A JPH1090087A (ja) | 1998-04-10 |
| JP3662682B2 true JP3662682B2 (ja) | 2005-06-22 |
Family
ID=17198943
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP24983896A Expired - Fee Related JP3662682B2 (ja) | 1996-09-20 | 1996-09-20 | 航空機主翼に作用する荷重の推定方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5922967A (ja) |
| JP (1) | JP3662682B2 (ja) |
Families Citing this family (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6354155B1 (en) * | 1999-06-02 | 2002-03-12 | Bertec Corporation | Multi-component force and moment measuring platform and load transducer |
| US6173618B1 (en) * | 1999-07-27 | 2001-01-16 | The United States Of America As Represented By The Department Of Health And Human Services | Ore pass level and blockage locator device |
| US6393898B1 (en) | 2000-05-25 | 2002-05-28 | Symyx Technologies, Inc. | High throughput viscometer and method of using same |
| US6664067B1 (en) | 2000-05-26 | 2003-12-16 | Symyx Technologies, Inc. | Instrument for high throughput measurement of material physical properties and method of using same |
| US6857309B2 (en) | 2001-08-24 | 2005-02-22 | Symyx Technologies, Inc. | High throughput mechanical rapid serial property testing of materials libraries |
| US6690179B2 (en) * | 2001-08-24 | 2004-02-10 | Symyx Technologies, Inc. | High throughput mechanical property testing of materials libraries using capacitance |
| US6860148B2 (en) | 2001-08-24 | 2005-03-01 | Symyx Technologies, Inc. | High throughput fabric handle screening |
| US6650102B2 (en) | 2001-08-24 | 2003-11-18 | Symyx Technologies, Inc. | High throughput mechanical property testing of materials libraries using a piezoelectric |
| US6736017B2 (en) | 2001-08-24 | 2004-05-18 | Symyx Technologies, Inc. | High throughput mechanical rapid serial property testing of materials libraries |
| US6772642B2 (en) | 2001-08-24 | 2004-08-10 | Damian A. Hajduk | High throughput mechanical property and bulge testing of materials libraries |
| US6769292B2 (en) * | 2001-08-24 | 2004-08-03 | Symyx Technologies, Inc | High throughput rheological testing of materials |
| US6837115B2 (en) | 2001-08-24 | 2005-01-04 | Symyx Technologies, Inc. | High throughput mechanical rapid serial property testing of materials libraries |
| US20030055587A1 (en) * | 2001-09-17 | 2003-03-20 | Symyx Technologies, Inc. | Rapid throughput surface topographical analysis |
| US7324007B2 (en) * | 2001-12-31 | 2008-01-29 | The United States Of America As Represented By The Secretary Of The Department Of Health And Human Services, Centers For Disease Control And Prevention | Instrumented rock bolt, data logger and user interface system |
| US7013709B2 (en) * | 2002-01-31 | 2006-03-21 | Symyx Technologies, Inc. | High throughput preparation and analysis of plastically shaped material samples |
| US7112443B2 (en) * | 2002-10-18 | 2006-09-26 | Symyx Technologies, Inc. | High throughput permeability testing of materials libraries |
| US8596135B2 (en) * | 2006-12-07 | 2013-12-03 | Technion Research and Dvelopment Foundation Ltd. | System and method for monitoring health of structural joints |
| KR101041332B1 (ko) | 2008-10-31 | 2011-06-14 | 현대건설주식회사 | 변형률 측정에 의한 구조물의 변위측정방법 및 이를 이용한변위측정장치 |
| WO2015047360A1 (en) * | 2013-09-29 | 2015-04-02 | Rinand Solutions Llc | Force sensing compliant enclosure |
| KR101645546B1 (ko) * | 2014-12-05 | 2016-08-04 | 연세대학교 산학협력단 | 변형률 측정을 이용한 임의 지점의 하중 추정 장치 및 방법 |
| US10013081B1 (en) | 2017-04-04 | 2018-07-03 | Google Llc | Electronic circuit and method to account for strain gauge variation |
| US10635255B2 (en) | 2017-04-18 | 2020-04-28 | Google Llc | Electronic device response to force-sensitive interface |
| US10514797B2 (en) * | 2017-04-18 | 2019-12-24 | Google Llc | Force-sensitive user input interface for an electronic device |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4493220A (en) * | 1982-11-23 | 1985-01-15 | Advanced Mechanical Technology, Inc. | Force measuring platform and load cell therefor using strain gages to measure shear forces |
| US4488441A (en) * | 1983-04-15 | 1984-12-18 | Jr3, Inc. | Apparatus for simultaneous measurement of mutually perpendicular forces and moments |
| DE3406059A1 (de) * | 1983-07-27 | 1985-02-07 | Siemens AG, 1000 Berlin und 8000 München | Mehrkomponentenkraft- und -momentenmesskoerper mit dehnungsmessstreifen |
| US4671118A (en) * | 1984-08-04 | 1987-06-09 | Yotaro Hatamura | Load sensor |
| US5461933A (en) * | 1994-02-04 | 1995-10-31 | Acutus Industries, Inc. | Shear web load cell having thermal compensation |
| US5500509A (en) * | 1994-03-01 | 1996-03-19 | Borg-Warner Automotive, Inc. | Microprocessor based universal digital pressure sensor |
| JP3131642B2 (ja) * | 1994-09-14 | 2001-02-05 | 日本電子工業株式会社 | 応力複合化センサ及びこれを用いた構造体の応力測定装置 |
| US5668325A (en) * | 1996-03-27 | 1997-09-16 | Cyprus Amax Coal Company | Method and apparatus for determining compressive stress in pillars |
-
1996
- 1996-09-20 JP JP24983896A patent/JP3662682B2/ja not_active Expired - Fee Related
-
1997
- 1997-09-19 US US08/934,419 patent/US5922967A/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH1090087A (ja) | 1998-04-10 |
| US5922967A (en) | 1999-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3662682B2 (ja) | 航空機主翼に作用する荷重の推定方法 | |
| CN109870134B (zh) | 一种旋转叶片非接触式动应变场测量方法及其*** | |
| CN108896230B (zh) | 一种基于有限元的螺栓紧固力超声检测及关键检测参数确定方法 | |
| Cornwell et al. | Application of the strain energy damage detection method to plate-like structures | |
| CN108254153B (zh) | 光纤气动力测量天平温度补偿方法 | |
| CN113237661B (zh) | 一种翼吊发动机的动力学试验载荷测量装置和方法 | |
| CN109933872B (zh) | 基于增强稀疏的复合材料结构冲击载荷识别方法及装置 | |
| CN113049217A (zh) | 一种大型风洞柔板结构多状态信息的动态监测方法 | |
| GB2189615A (en) | Improvements in dynamometers | |
| Ostachowicz et al. | Wave propagation in delaminated beam | |
| US3460383A (en) | Control surface having flexure areas for force measurements | |
| US3899922A (en) | Double lever strain multiplier apparatus and method | |
| US2785569A (en) | Wind tunnel force and moment measuring device | |
| Heaney et al. | Distributed sensing of a cantilever beam and plate using a fiber optic sensing system | |
| CN110823104B (zh) | 基于应变监测的工装构件关键区域位移场预测方法 | |
| CN109492323B (zh) | 基于变形测试的毛坯初始残余应力高精度反求方法 | |
| CN110823122A (zh) | 一种基于光纤光栅与弹性体的弯曲曲率测量装置及方法 | |
| JP2005172589A (ja) | 材料の試験方法、材料試験機 | |
| CN108195532B (zh) | 一种梁结构裂纹等效刚度的测量方法 | |
| Miller et al. | Evaluation of wing load calibration and sensing methods using conventional strain gages and a fiber optic sensing system installed on a straight tapered wing | |
| Wagner et al. | Analysis of several loading methods for simultaneous determination of Young's and shear moduli in composites | |
| Kalyukh | Applicability of the Quasistatic Approach to the Calculation of the Characteristics of an Accelerating Towed System | |
| KR100577349B1 (ko) | 교량 행어 케이블의 장력 측정장치 | |
| CN114964667B (zh) | 一种用于识别拱桥吊杆损伤的应变影响线方法 | |
| Nakamura | Estimation of dynamic load on a beam using central-difference scheme and FEM |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040416 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20040706 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040902 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050322 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050324 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080401 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090401 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090401 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100401 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110401 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |