低周波振動環境におけるエネルギー収集技術は、従来の線形構造では効率が悪く、特に振幅が小さい場合に限界があります。一方で、非線形な双稳態構造はこの課題を解決する可能性を持っています。COMSOL Multiphysicsにおいてそのモデル構築を試みる際には、構造力学・圧電効果・回路解析のマルチフィジックス結合を活用するのが一般的です。

まず、材料特性の定義より前に、双稳態機構の中心となる「予備曲げ梁」の形状を決定します。幾何構成として、互いに逆向きの弧状梁を2本配置し、これを「背中合わせ」の構造とします。パラメトリックカーブを使って、以下のような式で梁の形状を定義します：

curve1 = model.geom.create("curve1", "ParametricCurve");
curve1.set("x", "L/2*(1-cos(s*pi))");  
curve1.set("y", "-h0*sin(s*pi)");
curve2 = model.geom.create("curve2", "ParametricCurve");
curve2.set("x", "L/2*(1+cos(s*pi))");
curve2.set("y", "h0*sin(s*pi)");

このようにして作成された2つの弧形梁が、初期変位パラメータ h0 によってポテンシャルエネルギーの井戸深さに影響を与えます。また、圧電層の極化方向は誤って設定すると出力電圧の符号が逆になるため注意が必要です。過去にこのミスを経験したことがあるため、非常に重要なポイントです。

境界条件では、固定端にすべての自由度を拘束するのではなく、面内の滑りを許容する条件を適用することで、実際の支持状態に近いモデルを作成できます。さらに、初期変位の導入も重要です。ここでは初期の小さな変位を与えることで、システムを「刺激」し、双稳態遷移を誘導します。以下のようにガウス分布に基づいた初期変位を設定します：

model.physics("solid").feature("init").set("u0", "0.5*h0*exp(-(x^2+y^2)/0.1L^2)");

この初期条件により、より早く双稳態遷移が起こるようになります。解析手法としては、時間応答解析（瞬態解析）を採用し、時間ステップは初期は1e-4秒程度から試行的に設定します。変位が大きくなると自動的に可変ステップサイズへ切り替えるようにします。また、非線形変形を考慮するためには、幾何学的非線形性を有効にする必要があります。

結果の可視化では、単なる電圧の出力だけでなく、ポテンシャルエネルギーの分布図を確認することが重要です。典型的な二重ポテンシャル井戸構造が観察できれば、モデルの妥当性が確認できます。また、1次元プロットで梁の中央部の変位を時間軸に沿って描画すると、通常の正弦波のような波形が「圧縮された」様子になることが見られます。これは双稳態系の特徴的な弛緩振動であることを示しています。

出力電圧の変化が急激に増加するタイミングは、ポテンシャル井戸間の遷移が起きたことを意味します。このとき、負荷抵抗と圧電素子のインピーダンスを適切にマッチングさせることで最大出力を得ることができます。具体的には、圧電層の内部抵抗の2〜3倍の値を負荷として選択することが一般的です。ただし、外部振動周波数と構造体の固有周波数比が0.6〜1.2の範囲にある場合にのみ、双稳態構造の利点が顕著に現れます。

最後に、低周波の環境下でのエネルギー収集において、線形共振構造に頼るだけでは限界があることを再認識しましょう。COMSOL上で双稳態構造を用いたシミュレーションにより、5Hz程度の低周波でも効率的にエネルギーを捕らえることが可能になるのです。これはまさに、構造の非線形性を活かした高性能発電機の実現です。