(57)【要約】  本発明は、非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方を提供し、以下のステップを含む。フラッター時刻ｔ０時のモデル入力パラメータ取得し、モデル入力パラメータは２つの質量体の質量比、各質量体の位置座標、速度、加速度及び反発係数を含み、前記質量比及び反発係数はいずれも定数である。第１解モデルを確立し、フラッター終了時刻ｔ∞のパラメータを取得する。第２解モデルを確立し、フラッター終了時刻ｔ∞時の２つの質量体の位置パラメータを取得する。第３解モデルを確立し、第３解モデル及びモデル入力パラメータに基づき、フラッター終了時刻ｔ∞時の２つの質量体の速度パラメータを取得する。本発明は、あるフラッター時刻から数値シミュレーションを直接継続して終了し、フラッター過程をスキップしてフラッター終了の時刻とこの時刻の２つの衝突質量体の位置及び速度とを直接取得することができ、それにより計算精度を向上させ、且つ大量の計算時間を節約する。
2021525379.tif 000002
【特許請求の範囲】
【請求項１】
非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方であって、
フラッター時刻ｔ０時のモデル入力パラメータを取得し、前記モデル入力パラメータは２つの質量体の質量比、各質量体の位置座標、速度、加速度及び反発係数を含み、前記質量比及び反発係数はいずれも定数であるステップＳ１と、
第１解モデルを確立し、第１解モデル及びモデル入力パラメータに基づき、フラッター終了時刻ｔ∞のパラメータを取得するステップＳ２と、
第２解モデルを確立し、第２解モデル及びモデル入力パラメータを基づき、フラッター終了時刻ｔ∞時の２つの質量体の位置パラメータを取得するステップＳ３と、
第３解モデルを確立し、第３解モデル及びモデル入力パラメータを基づき、フラッター終了時刻ｔ∞時の２つの質量体の速度パラメータを取得するステップＳ４と、を含む
ことを特徴とする非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方。
【請求項２】
ステップＳ２において、第１解モデルは
2021525379.tif 000012
である
ことを特徴とする請求項１に記載の非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方。
【請求項３】
ステップＳ３において、第２解モデルは
2021525379.tif 000013
である
ことを特徴とする請求項１に記載の非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方。
【請求項４】
ステップＳ４において、第３解モデルは
2021525379.tif 000014
である
ことを特徴とする請求項１に記載の非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方。
【請求項５】
前記反発係数の値は０〜１の範囲をとる
ことを特徴とする請求項１に記載の非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、機械的運動数値シミュレーションの分野に属し、特に非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方に関する。
【背景技術】
【０００２】
衝突は、機械的な運動の過程で避けられない現象である。復元力の作用により、繰り返しの小さな振幅の衝突が発生することがよくあり、一般的にフラッターと呼ばれる。フラッター過程は、２つの衝突質量体が完全に同期運動して粘性を生じるまで続き、完全フラッターと呼ばれる。完全フラッターは、数え切れない回数衝突することを意味し、数値シミュレーションに大きな困難をもたらし、主に以下のように表現する。（１）フラッター振幅が小さく、数値シミュレーション丸め誤差が大きい。（２）衝突の回数が多く、エンドレスループに入りやすい。（３）小さい振幅の衝突を無視すると、発生する誤差は無視できなくなる。現在、固定拘束付きの衝突、すなわち衝突に関与する両質量体のうちの一方が固定されて動かない場合に対し、予備的な解決策がある。
【０００３】
しかし、最も一般的な二重衝突質量体はいずれも固定的に拘束されない状況に対し、その運動過程が複雑であるため、非常に簡略化された方法でシミュレーションすることがよく使用され、以下の欠陥をもたらす。（１）小さいフラッター運動を直接無視し、最終結果の誤差が大きくなる。（２）上記の比較的小さいフラッター運動のために計算時間を増加させ、計算時間が長く、効率が低い。そのため、現在、シミュレーション過程における欠陥を解決する効果的な方法がなく、非固定拘束フラッターシミュレーションの問題は、既に衝突研究の分野における技術的なネックとなっている。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明の目的は、フラッター過程をスキップしてフラッター終了の時刻とこの時刻での２つの衝突質量体の位置及び速度とを直接取得することができ、それによりシミュレーション計算がエンドレスループを超えて計算精度を大幅に向上させることができる非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方を提供することである。上記目的を達成するために、本発明が採用した技術的解決手段は以下のとおりである。
【課題を解決するための手段】
【０００５】
非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方であって、
フラッター時刻ｔ０時のモデル入力パラメータを取得し、前記モデル入力パラメータは２つの質量体の質量比、各質量体の位置座標、速度、加速度及び反発係数を含み、前記質量比及び反発係数はいずれも定数であるステップＳ１と、
第１解モデルを確立し、第１解モデル及びモデル入力パラメータに基づき、フラッター終了時刻ｔ∞のパラメータを取得するステップＳ２と、
第２解モデルを確立し、第２解モデル及びモデル入力パラメータを基づき、フラッター終了時刻ｔ∞時の２つの質量体の位置パラメータを取得するステップＳ３と、
第３解モデルを確立し、第３解モデル及びモデル入力パラメータを基づき、フラッター終了時刻ｔ∞時の２つの質量体の速度パラメータを取得するステップＳ４と、を含む
ことを特徴とする非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方。
【０００６】
好ましくは、ステップＳ２において、第１解モデルは
2021525379.tif 000003
である。
【０００７】
好ましくは、ステップＳ３において、第２解モデルは
2021525379.tif 000004
である。
【０００８】
好ましくは、ステップＳ４において、第３解モデルは
2021525379.tif 000005
である。
【０００９】
好ましくは、前記反発係数の値は０〜１の範囲をとる。
【発明の効果】
【００１０】
従来技術と比べ、本発明の利点は、あるフラッター時刻から数値シミュレーションを直接継続して終了し、フラッター過程をスキップしてフラッター終了の時刻とこの時刻での２つの衝突質量体の位置及び速度とを直接取得することができることであり、それにより計算精度を向上させ、且つ大量の計算時間を節約する。
【図面の簡単な説明】
【００１１】
【図１】本発明の一実施例の非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方のフローチャートである。
【図２】本発明の一実施例のバネプロトンシステム図である。
【図３】本発明の一実施例における両質量体の相対変位の時間変化曲線である。
【図４】本発明の一実施例における両質量体の相対速度の時間変化曲線である。
【発明を実施するための形態】
【００１２】
以下、模式図を参照して本発明の非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方についてより詳細に説明し、ここで、本発明の好ましい実施例を示しており、本発明の有利な効果を実現する限り、当業者であれば、ここで説明された本発明を修正することができることを理解されるべきである。そのため、以下の説明は、当業者であれば広く理解されるべきであり、本発明を限定するものではない。
【００１３】
図１に示すように、非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方であって、ステップＳ１〜Ｓ４を含み、具体的には以下のとおりである。
【００１４】
2021525379.tif 000006
【００１５】
ステップＳ２：第１解モデルを確立し、第１解モデル及びモデル入力パラメータに基づき、フラッター終了時刻ｔ∞のパラメータを取得する。第１解モデルは式（１）に示すようになる。
2021525379.tif 000007
【００１６】
ステップＳ３：第２解モデルを確立し、第２解モデル及びモデル入力パラメータに基づき、フラッター終了時刻ｔ∞時の２つの質量体の位置パラメータを取得する。第２解モデルは式（２）〜（３）に示すようになる。
2021525379.tif 000008
【００１７】
ステップＳ４：第３解モデルを確立し、第３解モデル及びモデル入力パラメータに基づき、フラッター終了時刻ｔ∞時の２つの質量体の速度パラメータを取得する。第３解モデルは以下のようになる。
2021525379.tif 000009
【００１８】
本実施形態では、反発係数ｒの値は０〜１の範囲をとる。
【００１９】
図２に示すように、バネ１、台車２、２つの質量体で１つのバネプロトンシステムを構成する。台車２の質量はＭであり、ボール３は台車２の内部のキャビティ４に置かれる。ボール３は水平方向に自由移動可能であり、ボール３の質量はｍである。台車２の変位をｘとし、ボール３の変位をｙとし、いずれも水平右方向を正方向とし、バネ１が平衡位置にある時に、ｘ座標及びｙ座標の原点は台車２の中心と一致する。バネ１の剛性Ｋ＝１０００Ｎ・ｍ−１、減衰係数Ｃ＝１０Ｎ・ｓ・ｍ−１、台車２の質量Ｍ＝１０ｋｇ、ボール３の質量ｍ＝０．５ｋｇ、台車キャビティ４の自由長５ｍｍ、加振力Ｆ０＝１００Ｎ、加振周波数Ω＝１０Ｈｚである。
【００２０】
従来技術における数値シミュレーション方法を用い、ｔ＝０．１１８１０３１ｓである時にボール３が台車２に衝突したことがわかり、この時刻をｔ０時刻とし、この時刻の台車２及びボール３の運動パラメータを表１に示す。それぞれ数値シミュレーション方法及び本発明の方法を用い、フラッター終了時刻ｔ∞＝０．１３７９１２７ｓを得て、この時刻の台車２及びボール３の運動パラメータも表１に示す。
【００２１】
2021525379.tif 000010
【００２２】
表１から明らかなように、本発明の方法を用いて得られる計算結果の精度は高く、変位及び速度の相対誤差はいずれも０．０７％未満である。計算時間では、数値シミュレーションを用いる時、１回のフラッターの計算時間は約５０ｓであるが、本発明の方法を用いると、フラッター終了時の運動パラメータを直接得ることができる。上記のフラッター過程は、数値シミュレーションにおいて数百〜数千万回繰り返し発生するため、全体的に見ると、計算効率を大幅に向上させることができる。
【００２３】
2021525379.tif 000011
【００２４】
本実施例では、この方法は以下の利点を有する。（１）この方法は一般性があり、利用範囲が広い。この方法は、２つの衝突質量体がいずれも非固定拘束であって、２つの衝突質量体がいずれも可変加速度を有することを許容する場合に適用可能であり、そのため、一般的なフラッター衝突の場合に属する。また、以下のような特別な場合もある。一般的なボールが落下して跳ね返って衝突するのは、一方が固定拘束（大地）であり、他方が固定加速度（重力加速度）である場合に属し、本方法の応用の特殊な場合であり、一般的なバネプロトンシステムと固定バッフルとの衝突であれば、一方が固定拘束（バッフル）であり、他方が可変加速度である場合に属し、本方法の応用の特殊な場合でもある。以上の２つの場合は、いずれも質量比μｍを＋∞とし、固定拘束の初期速度、初期加速度はいずれも０に設定することができる。本発明の方法を用いて計算し、フラッター終了時、２つの衝突質量体の速度はいずれも０であり、固定拘束の解に収束する。従って、固定拘束の場合は本発明の特殊な場合であると考えられる。本発明は固定拘束の制限がなく、最も一般的な非固定拘束の場合に適用することができるため、広い使用範囲を有する。
【００２５】
（２）この方法は、解を求める速度が速く、効率が高い。
式（１）〜（４）は、ｔ０時刻の位置、速度と加速度及びｒとμｍの２つの定数のみを使用し、従って、ｔ０時刻において衝突によるフラッター終了時の時刻とその時２つの衝突質量体の位置及び速度とを計算することができる。従来の方法では、計算精度を向上させるために、反復ステップは一般的に小さい値をとり、多くの逐次反復をもたらし、計算時間が長い。本発明は、１つのステップで最終結果を得ることを実現すると同時に、従来の計算において克服することが困難な無限回数の衝突によるエンドレスループになる問題を回避する。
【００２６】
3.  この方法は計算精度が高く、誤差が小さい。
【００２７】
加速度が一定のシステムに対して、本発明は何ら誤差を持たない。加速度可変のシステムに対して、ｔ０時刻の２つの衝突質量体の相対速度が十分に小さいため、フラッター過程全体が継続する時間は十分に小さく、この期間は加速度が一定であると仮定してもよい。一般的には、フラッター期間の加速度の変化量が１％未満であれば、本方法による計算結果の精度を保証することができる。このことは、非線形衝突等の初期値に対する敏感な問題の数値シミュレーションにとって、非常に必要である。
【００２８】
以上のように、本発明の実施例が提供した非固定拘束両質量体に基づく完全フラッター終了パラメータの求め方において、あるフラッター時刻から数値シミュレーションを直接継続して終了し、フラッター過程をスキップしてフラッター終了の時刻とこの時刻での２つの衝突質量体の位置及び速度とを直接取得することができ、それにより計算精度を向上させ、且つ大量の計算時間を節約する。
【００２９】
上記は本発明の好ましい実施例にすぎず、本発明を何ら制限するものではない。当業者であれば、本発明の技術的解決手段から逸脱しない範囲で、本発明に開示された技術的解決手段及び技術的内容をいずれかの形で置換又は修正するなどの変更は、いずれも本発明の技術的解決手段の内容から逸脱することなく、本発明の技術的範囲に含まれる。
【符号の説明】
【００３０】
１  バネ、２  台車、３  ボール、４  キャビティ。
【図１】
2021525379.tif 000015
【図２】
2021525379.tif 000016
【図３】
2021525379.tif 000017
【図４】
2021525379.tif 000018
【国際調査報告】
2021525379.tif 000019
2021525379.tif 000020
2021525379.tif 000021
2021525379.tif 000022
2021525379.tif 000023