Hybrid Dynamics of Rocking Blocks Beyond Overturning: Saltation Analysis, Bifurcations, and Stability Characterization

要約

地面の上に置かれた、重くて長方形をした石の塊を想像してみてください。もし地面を前後に揺らすと、その石の塊は角の部分で、まるでシーソーに乗っている子供のように、前後に揺れ動くようになります。エンジニアたちは、地震の際に建物や彫像、給水塔などが倒れないようにする方法を理解するために、数十年にわたりこの現象を研究してきました。

長い間、科学者たちは、ブロックが地面に当たるたびにどれだけのエネルギーを失うかを正確に予測するために、「経験則」（ハウズナーという人物による公式）を使用してきました。これはバスケットボールのゲームのようなものです。ボールを落とすと跳ね返りますが、床にエネルギーを奪われるため、最初はそれほど高くは跳ね返りません。ハウズナーのルールは、この「跳ね返り」（あるいは揺れ）がどの程度になるかを計算するための特定のやり方です。

しかし、近年の実験により、現実の世界は必ずしもハウズナーのルール通りにはいかないことが示されました。より現実に即した、より複雑な公式（マオとその同僚らによるもの）が作成されました。

この論文は、大きな問いを投げかけています。「実際に、どちらのルールを使うかが重要なのでしょうか？」 もし古いルールと、より正確な新しいルールを使い分けた場合、ブロックの将来の挙動は変わってしまうのでしょうか？

研究者たちの発見を、簡単に説明します：

1. 「跳ね返り」が物語を完全に変える

研究者たちは、地面をさまざまな強さで揺らす数千通りのシナリオをシミュレーションしました。その結果、「跳ね返りのルール」の選択が、ブロックに何が起こるかという物語を完全に変えてしまうことがわかりました。

• 古いルール（ハウズナー）： ブロックが地面に当たるとき、多くのエネルギーを失うと予測します。それは重く、鈍い衝撃音のように振る舞います。ブロックは安定した、予測可能なリズムの中に落ち着く傾向があります。
• 新しいルール（マオ）： ブロックが失うエネルギーは「より少ない」と予測します。それは、より活発な跳ね返りのように振る舞います。ブロックはより多くのエネルギーを保持するため、より早く「興奮（激しく揺れ動く状態）」状態になります。

2. 「背が高いか、低いか」によるブロックの違い

論文では、ブロックの形状が非常に重要であることが判明しました。

• 低くて太いブロック（低スレンダーネス）： ルール間の食い違いが最も大きいのは、これらのブロックです。もし、短くて太いブロックに対して、より正確な新しいルールを使用した場合、そのブロックはより早く混沌とした状態になります。複雑で予測不可能なパターンで揺れ始め、古いルールが予測していたよりもはるかに転倒・落下しやすくなります。古いルールは、実態よりもブロックが安全であると告げる「嘘」をついていたことになります。
• 背が高くて細いブロック（高スレンダーネス）： ブロックがより高く、細くなるにつれて、2つのルールは一致し始めます。「鈍い衝撃」と「活発な跳ね返り」の違いは、これらの高いブロックにおいては重要性が低くなります。どちらのルールも、同様の挙動を予測します。

3. 「壊れやすいバランス」の比喩

トランプの束を積み上げる場面を想像してみてください。

• 古いルールの場合： 積み上げられたカードは安定しているように見えます。少し突っついても、ゆらゆらと揺れるだけで、そのまま立っています。「安全圏（倒れない範囲）」は広く、簡単に見つけることができます。
• 新しいルールの場合： 同じ積み上げられたカードが、より脆弱になります。その「安全圏」は、細く、うねった筋のような形になります。もしブロックを少しでも異なる方向に突っつくと（ほんのわずかな差であっても）、突然崩壊してしまうかもしれません。

研究者たちは、短いブロックの場合、新しいルールを用いると「安全圏」が実際には非常に小さく、断片化していることを明らかにしました。これは、現実の世界において、短いブロックはどのように動き始めるかに対して非常に敏感であることを意味します。ブロックへの「突き方」のわずかな違いが、安全な揺れと、壊滅的な転倒の分かれ目となるのです。

4. 主な教訓

この論文は、「跳ね返り」を単なる小さな局所的な詳細として扱うべきではないと結論付けています。それは、システムの「個性」の根本的な部分なのです。

• 短いブロックに対して： 古い単純なルールを使うことは、偽りの安心感を与えます。新しいルールは、それらがより混沌としており、より転倒しやすいことを示しています。
• 高いブロックに対して： 2つのルールが一致するため、古いルールでも「十分」です。

要するに、もしあなたが重い物体が地震を生き延래られるかどうかを予測しようとしているなら、短くて幅の広い物体に対して、単に古い単純な数学を用いることはできません。より複雑で正確な数学を用いる必要があります。さもなければ、物体が転倒の瀬戸際にあるにもかかわらず、安全であると思い込んでしまうことになるでしょう。物体が地面にどのように当たるかは、その後の運命を決定づけるのです。

技術概要

技術要約：転倒を超えた領域におけるロッキング・ブロックのハイブリッド力学

問題提起
ロッキング・ブロックは、石造りの壁、彫像、高架水槽などの剛体構造物の地震応答を理解するための基礎的なモデルである。従来、転倒を防ぐことが主要な工学的目的であったが、崩壊に至る前の非線形振動は、多周期的なレジームやカオス的なレジームを含む複雑な挙動を示す。このダイナミクスを支配する決定的な要因は、衝突時のエネルギー散逸であり、これは反発係数（COR）を介してモデル化される。

古典的なハウスナー（Housner）モデルは、角運動量保存則から導出されたものであり、数十年にわたり標準として用いられてきたが、近年の研究では、その予測と実験的な衝突挙動との間に乖離があることが指摘されている。Maoらによって提案された運動エネルギー再分配モデルのような代替の定式化は、特に細長比の低いブロックにおいて、実験データにより適合することが示唆されている。しかし、これらの代替的な反発定式化を採用することが、ロッキング系のグローバルなダイナミクス（具体的には、アトラクタの安定性、分岐構造、および有界運動へのアクセシビリティ）に与える影響については、ほとんど解明されていない。

手法
著者らは、調和基盤励起を受けるロッキング・ブロックを、ハイブリッド非平滑力学系として定式化している。系は、質量 $M$、幅 $2b$、高さ $2h$ を持つ矩形の剛体で構成され、剛体基礎上で滑ることなくロッキングを行う。ダイナYNAMICSは、活動的なピボットに基づいて切り替わる、区分的な自律微分方程式によって支配される。

解析にあたり、著者らは2つの異なる反発モデルを比較している：

1. ハウスナー・モデル ($e_H$): 角運動量保存に由来し、ブロックの幾何学的形状（臨界角 $\theta_c$）のみに依存する。
2. Maoらによるモデル ($e_M$): エネルギー散逸と接触界面における弾性回復を考慮した、速度場のベクトルベースの解析から導出される。

本研究では、以下の手法を用いた包括的な数値的特性評価を行っている：

• 分岐図: 励起振幅 ($\alpha$) の変化に伴う漸近応答（周期、多周期、カオスの、転倒）の定性的な変化をマッピングする。
• 最大リアプノフ指数 ($\lambda_{max}$): 不連続点での塩殺行列（saltation matrices）を用いたベンベティン・ウルフ（Benettin–Wolf）アルゴリズムを用いて計算し、カオス的レジームと安定な周期レジームを識別する。
• 吸引盆（Basins of Attraction）: 初期条件空間を離散化することで、異なる力学的帰結（減衰運動、持続振動、転倒）へのグローバルなアクセシビリティを構築する。

シミュレーションは、同一の幾何学的および励起パラメータの下で、2つの細長比（$h/b = 1$ および $h/b = 2$）に対して実施され、反発則がシステムの進化に与える影響を分離して抽出している。

主な貢献および結果
本研究は、反発モデルの選択が、特に細長比の低いブロックにおいて、予測されるダイナミクスの景観を著しく変化させることを示している。

• 低細長比ブロック ($h/b = 1$) における乖離:

  • Maoらの定式化は、ハウスナー・モデルと比較して高い反発係数（低いエネルギー散逸）を予測し、その対称的なパーセンテージ差は約79.85%に達する。
  • この散逸の減少により、複雑な振動挙動の早期の発現と、カオス力学のためのより広いパラメータ範囲がもたらされる。
  • 吸引盆の解析: Maoモデルは、安定な対称振動の吸引盆を大幅に収縮させ（初期条件の約55%から約20%へ減少）、転倒の結果が生じる割合を著しく増加させる（約44%から約80%へ）。振動の吸引盆の境界は断片化し不規則になり、初期条件に対する感受性の高まりを示している。
  • 両モデルによるダイナミカルな帰結の全体的な一致度は、$h/b=1$ においては**64.09%**のみであり、これは初期条件の3分の1以上が、使用する衝突モデルのみによって定性的に異なる運命を辿ることを意味している。

• 高細長比ブロック ($h/b = 2$) における収束:

  • 細長比が増加するにつれ、2つの反発係数の間の乖離は減少する（対称差は約17%に低下）。
  • その結果、分岐構造およびリアプノフ指数の分布は、両モデルで漸次的に収束する。
  • Maoモデルは依然として有界領域内に多周期アトラクタを導入するものの、モデル間のグローバルな一致度は**86.93%**に達しており、背の高いブロックは反発係数の特定の定式化に対して感度が低いことを示唆している。

意義と主張
本論文は、反発係数を単なる単一衝突におけるエネルギー損失を制御する局所的なパラメータとしてではなく、位相空間のグローバルな組織化を根本的に変え得るモデリング上の仮定として捉えるべきであると論じている。

著者らは、実験的な衝突データに適合させるために採用されることの多い代替の反発定式化が、古典的なハウスナー・モデルと比較して、著しく異なる長期予測を生み出す可能性があると主張している。具体的には、モデルの選択は以下に影響を与える：

1. 安定性の閾値および力学的遷移（分岐）の発生位置。
2. 安全な有界ロッキングのための吸引盆のサイズと堅牢性。
3. 同一の励起条件下における転倒への感受性。

本研究は、低細長比のブロックにおいては、衝突モデルがシステムの安全性と挙動の決定的な要因となる一方で、背の高いブロックにおいては、ダイナミクスの特徴が反発モデリングの変動に対してよりロバストであることを結論付けている。著者らは、本解析が理想的な剛体近似に基づいていること、および材料の変形、滑り、あるいは粘弾性散逸を考慮していないことを述べ、これらが将来の改良に向けた領域であることを明記し、その範囲を限定的に設定している。