"構造物の強度の変化を検出する方法であって、構造物の振動を検出して、検出した振動を表す振動波形を発生させ、振動波形を所定長さの短区間の波形に分割し、各短区間波形を分析して、含まれる周波数成分のうちパワーが最大となる最大周波数成分の周波数及びパワーが最小となる最小周波数成分の周波数を求め、パワーが最大となる最大周波数成分の周波数の検出頻度を示す最大周波数成分周波数頻度分布及びパワーが最小となる最小周波数成分の周波数の検出頻度を示す最小周波数成分周波数頻度分布を求め、最大周波数成分周波数頻度分布及び最小周波数成分周波数頻度分布の少なくとも1つの時間変化から構造物の強度の変化を検知するステップより成る構造物の強度変化検出方法。短区間波形の所定長さをTとすると、短区間波形の分析は、kを1より小さい正の係数としてk/Tの整数倍の周波数で与えられる周波数成分を求める非調和周波数分析であることを特徴とする請求項1に記載の方法。パワーが最大または最小となる最大または最小周波数成分の周波数の検出頻度を示す頻度分布の標本の数の大きさは、最大周波数成分周波数頻度分布と最小周波数成分周波数頻度分布の積が所望の観測周波数帯域で周波数に関してほぼ一定の値になる時の数で与えられることを特徴とする請求項1または2に記載の方法。建造物の点検または監視は従来の検査法によるのであれば専門技術を必要とし、時間も費用も嵩む。本発明の目的は、専門技術を必要とせず、客観的且つ定量的に構造物の強度変化を検出し、構造物の強度変化の点検または常時監視を可能にする方法を提供することにある。上記目的を達成するために、本発明は、構造物の強度の変化を検出する方法であって、構造物の振動を検出して、検出した振動を表す振動波形を発生させ、振動波形を所定長さの短区間の波形に分割し、各短区間波形を分析して、含まれる周波数成分のうちパワーが最大となる最大周波数成分の周波数及びパワーが最小となる最小周波数成分の周波数を求め、パワーが最大となる最大周波数成分の周波数の検出頻度を示す最大周波数成分周波数頻度分布及びパワーが最小となる最小周波数成分の周波数の検出頻度を示す最小周波数成分周波数頻度分布を求め、最大周波数成分周波数頻度分布及び最小周波数成分周波数頻度分布の少なくとも1つの時間変化から構造物の強度の変化を検知するステップより成る構造物の強度変化検出方法を提供する。また、本発明は、パワーが最大または最小となる最大または最小周波数成分の周波数の検出頻度を示す頻度分布の標本の数の大きさが、最大周波数成分周波数頻度分布と最小周波数成分周波数頻度分布の積が所望の観測周波数帯域で周波数に関してほぼ一定の値になる時の数で与えられることを特徴とする構造物の強度の変化の検出方法を提供する。橋、ビル、ダムなどの構造物は、その大きさや重さに関係なく風や水流あるいは地盤の振動によって常に微小な振動をしている。本発明による構造物の強度変化検出方法によると、上記のような自然の力によって振動している構造物の雑多な振動波形から構造物の固有振動特性を推定することが可能となり、その特性の変化から構造物の強度変化を検知できるという効果がある。本発明による構造物の強度変化検出方法を従来の検査法と併用すると、点検する建造物を選択して重点的な検査を行うことが可能になり、点検と検査に要する時間と費用を大幅に削減できるという効果がある。本発明による構造物の強度変化検出方法は、建造物の振動を検出する振動検出器、振動波形をデジタル信号に変換するAD変換器、信号処理用コンピューター及び出力表示装置から成るシステムにより実施可能である。構造物の振動特性を推定する周波数帯域の下限周波数をGとする時、短区間波形の時間長さTを1/G、mを10以上で150以下の正の整数、k=0.1として、mk/Tで与えられる周波数を分析周波数とし、短区間波形の数N頻度分布の標本数を最大周波数成分周波数頻度分布と最小周波数成分周波数頻度分布の積が所望する観測周波数帯域で周波数に関してほぼ一定の値になる時の数とする。構造物の基本周期が10秒なら1/Gは10秒、1秒なら1/Gは1秒程度とするのが適当である。Nは周囲の振動環境により異なるが、2000乃至4000が典型的な数である。図1は、本発明の一実施例の構成を示すブロック図である。構造物の振動を検出する振動検出器101の出力電気信号は前置増幅器102で増幅され、ローパスフィルター103を経て、AD変換器104によりデジタル信号に変換され、コンピューター105により処理される。コンピューター105による信号の処理は、図2に示すように、デジタル信号の短区間波形短いサンプル列への分割201、分割された短区間波形ごとの非調和周波数分析非調和フーリエ解析202、分析により求められた周波数成分の中からの最大パワーと最小パワーを有する周波数成分の検出203、検出された周波数成分の各周波数の記録204、こうして記録された周波数に基づく頻度分布の確定の順序で行なわれ、その結果が表示部106に表示される。上記の頻度分布は建物を例とすると、一度だけ最大周波数成分周波数頻度分布と最小周波数成分周波数頻度分布からNの数を決めれば、それ以降は同じ標本数Nで頻度分布を求めて、頻度分布の変化を監視するようにしてもよい。また、最大周波数成分周波数頻度分布と最小周波数成分周波数頻度分布は互いに逆特性の関係を有するため、どちらか一方だけを用いて頻度分布の変化を監視してもよい。コンピューター105で行なわれる信号処理を数式により説明すると、以下のようになるが、これらはソフトウェアで行なう。AD変換器の出力デジタル信号である構造物の振動波形から得た短区間波形のサンプル列をWr=1,2,,K、分析周波数をfmとして、係数aおよびbをそれぞれ以下の式で求める。KKa=ΣWsin/Σsin2r1r1KKb=ΣWcos/Σcos2r1r1ただし、x=r-1/2-K/2、Fは標本化周波数とする。次に、Q=a2b2として、Qが最大となる周波数をfp、最小となる周波数とをfqとすると、構造物の振動波形から得た全短区間波形よりfpの頻度分布最大周波数成分周波数頻度分布およびfqの頻度分布最小周波数成分周波数頻度分布が得られる。図3は、最大周波数成分周波数頻度分布の時間変化を示す例である。これは構造物の模型に雑多な振動非定常振動を加え、構造物の強度を変えて求めたものである。各頻度分布の標本の数Nは2400である。図4は、図3の頻度分布を求めた振動波形を、従来法の周波数分析フーリエ解析、DFTを使って分析して、図3の頻度分布を得るために分割した短区間波形の長さの10倍の長さの区間波形240個のスペクトルを求めて平均した特性曲線の変化を示したものである。図3と図4の比較から、雑多な振動の周波数特性の影響は、従来法による振動特性の推定に比べて本発明による方法の方が少ないことが分かる。図3において、たとえば時刻n=1における頻度分布をy1、時刻nにおける頻度分布をynm=10,11,,150とし、y1を基準の分布としたynの変化の度合いEを150150E=Σyn-y12/Σyn2y12m10m10で与えるか、y1とynの一致の度合いVを150150V=2Σy1yn/Σy12yn2m10m10で与えれば、頻度分布の変化を定量的に示すことができる。同図において、n=2でE=0.02、V=0.98、n=15でE=0.38、Vが所定の値を超えた時、あるいはVが所定の値以下になった時、建造物の強度に有意または重大な劣化が生じたとして警報を自動的に発生することができる。また、飛行中の機体の振動を検出して、その振動波形から本発明の方法で頻度分布を求めることにより、機体の強度変化の有無を航空機の飛行中に常時監視することが可能である。雑多な振動下にあるパワープラントにおいて、様々な構造物の振動を検出し、その振動波形から本発明の方法で頻度分布を求めるなら、それぞれの部分の強度変化の有無をプラントの稼動中に常時監視することができる。"