圧電解析（静的）
理論と物理
圧電効果
先生、圧電効果って機械的な力を電圧に変換することですか？
正圧電効果（力→電圧）と逆圧電効果（電圧→変形）の両方がある。
先生、圧電効果って機械的な力を電圧に変換することですか？
正圧電効果（力→電圧）と逆圧電効果（電圧→変形）の両方がある。
構成方程式:
$\{T\}$: 応力、$\{S\}$: ひずみ、$\{E\}$: 電界、$\{D\}$: 電束密度、$[c^E]$: 弾性定数（定電界）、$[e]$: 圧電定数、$[\varepsilon^S]$: 誘電率（定ひずみ）。
構造と電磁場が連成するんですね。
そう。FEMでは変位$u$と電位$\phi$を同時に未知数として解く圧電連成解析になる。
主要な圧電材料
| 材料 | $d_{33}$ [pC/N] | 用途 |
|---|---|---|
| PZT (チタン酸ジルコン酸鉛) | 300〜600 | アクチュエータ、センサ |
| BaTiO₃ | 190 | セラミックコンデンサ |
| PVDF | -33 | フレキシブルセンサ |
| AlN | 5 | MEMS共振器 |
| LiNbO₃ | 6 | SAWフィルタ |
まとめ
圧電効果の発見——ピエール・キュリーとポール・キュリーの1880年の実験
圧電効果（Piezoelectric Effect）は1880年にフランスのピエール・キュリー（Marie Curieの夫）と兄ポール・ジャック・キュリーが石英結晶で発見した。彼らはクォーツに力を加えると表面に電荷が現れる「正圧電効果」を実証し、翌1881年にガブリエル・リップマンが逆の「逆圧電効果（電場を加えると変形する）」を理論予測した。圧電効果の工業応用は20世紀初頭のソナー（水中音波探知）から始まり、第一次世界大戦でPaul Langevinが圧電クォーツを使った水中探信儀を開発した。現代のスマートフォンのカメラ手振れ補正アクチュエータ、医療用超音波診断装置、インクジェットプリンタのヘッドはすべてキュリー兄弟の発見の子孫だ。
各項の物理的意味
- 電場項 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：ファラデーの電磁誘導法則。時間変動する磁束密度が起電力を生じさせる。【日常の例】自転車のダイナモ（発電機）は、磁石を回転させることで近くのコイルに電圧が発生する——磁場が時間的に変化すると電場が誘起されるというこの法則の直接的応用。IHクッキングヒーターも同じ原理で、高周波磁場の変化が鍋底に渦電流を誘起し、ジュール熱で加熱する。
- 磁場項 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：アンペア-マクスウェルの法則。電流と変位電流が磁場を生成する。【日常の例】電線に電流を流すと周囲に磁場が生じる——これがアンペアの法則。電磁石はこの原理で動作し、コイルに電流を流して強力な磁場を作る。スマートフォンのスピーカーも、電流→磁場→振動板の力というこの法則の応用。高周波（GHz帯のアンテナ等）では変位電流 $\partial D/\partial t$ が無視できなくなり、電磁波の放射を記述する。
- ガウスの法則 $\nabla \cdot \mathbf{D} = \rho_v$：電荷が電束の発散源であることを示す。【日常の例】下敷きで髪の毛をこすると静電気で髪が逆立つ——帯電した下敷き（電荷）から電気力線が放射状に広がり、軽い髪の毛に力を及ぼす。コンデンサ（キャパシタ）の設計では、電極間の電場分布をこの法則で計算する。ESD（静電気放電）対策もガウスの法則に基づく電場解析が基盤。
- 磁束保存 $\nabla \cdot \mathbf{B} = 0$：磁気単極子が存在しないことを表す。【日常の例】棒磁石を半分に割っても、N極だけ・S極だけの磁石は作れない——必ずN極とS極がペアで存在する。これは磁力線が「始点も終点もない閉じたループ」を描くことを意味する。数値解析では、この条件を満たすためにベクトルポテンシャル $\mathbf{B} = \nabla \times \mathbf{A}$ という定式化を用い、磁束保存を自動的に保証する。
仮定条件と適用限界
- 線形材料仮定：透磁率・誘電率が磁場・電場強度に依存しない（飽和領域では非線形B-Hカーブが必要）
- 準静的近似（低周波）：変位電流項を無視可能（$\omega \varepsilon \ll \sigma$）。渦電流解析で一般的
- 2D仮定（断面解析）：電流方向が一様で、端部効果を無視できる場合に有効
- 等方性仮定：異方性材料（珪素鋼板の圧延方向等）では方向別の特性定義が必要
- 適用外ケース：プラズマ（電離気体）、超伝導体、非線形光学材料では追加の構成則が必要
数値解法と実装
圧電FEMの定式化
離散化すると:
$[K_{uu}]$: 機械剛性、$[K_{\phi\phi}]$: 誘電剛性、$[K_{u\phi}]$: 圧電連成項。
構造と電気の自由度が1つの行列に統合されるんですね。
各節点に変位DOF（$u_x, u_y, u_z$）と電位DOF（$\phi$）の合計4つ（3D）。Abaqusの圧電要素C3D8Eなどが対応。
まとめ
圧電FEMの連成設定——機械-電気連成方程式の弱定式化と材料テンソルの扱い
圧電材料のFEM解析は「機械-電気連成問題」で、変位u（機械場）と電気ポテンシャルphi（電界場）を連立して解く。弱定式化では圧電構成方程式（応力テンソル＝弾性定数×歪み－圧電定数×電界）を仮想仕事原理に組み込み、剛性行列K（機械）、誘電行列（電気）、連成行列（圧電）の3ブロック行列で表す。実装上の落とし穴は圧電テンソルeの座標変換——結晶軸方向と解析座標系のミスマッチがあると連成が弱くなり電圧出力が過小評価される。COMSOLのPiezoelectric Devices PhysicsとANSYS MechanicalのPIEZOモジュールはこの変換を自動処理するが、ユーザー定義材料では変換行列を手動で設定する必要がある。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実務
超音波トランスデューサ、圧電アクチュエータ（インクジェットヘッド）、振動センサ、エナジーハーベスタの設計。
チェックリスト
超音波洗浄機の圧電トランスデューサ設計——キャビテーション分布とFEM最適化
超音波洗浄機は圧電トランスデューサが発生する超音波振動で液中にキャビテーション気泡を生成し、その崩壊衝撃で汚れを除去する。トランスデューサの共振周波数設計（20〜40kHz）と洗浄槽内の音圧分布が洗浄均一性を決定する。FEM（圧電連成）でトランスデューサの振動モードと音響場（Helmholtz方程式）を解析し、音圧が均一になる配置と形状を最適化する。実際の設計事例では、4個の円板型トランスデューサの配置を非対称にすることでキャビテーション分布の均一性を30%改善し、精密洗浄ラインでの洗浄不良率を0.8%から0.1%以下に低減した。FEMによる設計最適化が半導体ウェーハ洗浄・医療機器洗浄の品質管理に直結している。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| COMSOL (Piezoelectric) | 最も使いやすい。GUIで分極方向設定 |
| Abaqus | C3D8E等の圧電要素。研究標準 |
| Ansys Mechanical | SOLID226/227圧電要素 |
| OnScale (旧PZFlex) | 超音波トランスデューサ専用 |
圧電FEMツール比較——COMSOL Multiphysics vs ANSYS Mechanical PZEの機能差
圧電解析（機械-電気連成FEM）の主要ツールを比較する。COMSOL Multiphysicsの「Piezoelectric Devices」モジュールは連成問題の設定が最も直感的で、周波数応答・固有値・時間依存解析をGUIで切り替えられる。材料データベースにPZT-4/5H/5A等の標準圧電材料が事前登録されているため、初期設定の手間が少ない。ANSYS Mechanicalの圧電要素（SOLID226, SOLID227）は大規模問題の並列計算性能が高く、複数の圧電素子を含むシステムレベルのFSI（流体-構造-電気連成）に向く。OpenFOAMは圧電連成ソルバーを標準では持たないため、FEniCSやCode_Asterなどの学術用FEMコードが圧電研究には使われる。
選定で最も重要な3つの問い
- 「何を解くか」：圧電解析（静的）に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端
圧電MEMSの最前線——ナノ発電素子とAIoTセンサへの応用
圧電効果を利用した「エネルギーハーベスティング（振動発電）」が、バッテリーレスのAIoTセンサ電源として注目されている。PZT（チタン酸ジルコン酸鉛）やAlN（窒化アルミニウム）薄膜を用いた片持ち梁型MEMSが振動エネルギーを電力に変換し、環境振動（橋梁・機械の振動）から数uW〜数mWの電力を収穫する。FEM解析（機械-電気連成）で梁の共振周波数と電力出力を最適化するが、圧電材料の異方性と機械境界条件（固定端の残留応力）の適切なモデル化が精度の鍵だ。ETH ZurichとImecが共同開発した広帯域圧電ハーベスター（複数固有振動数チューニング）では、FEM最適化により環境振動の周波数変動（25〜70Hz）に対して電力収穫効率を従来比2.5倍に向上させた。
トラブルシューティング
トラブル
圧電FEMの固有振動数が実測と5%ズレる——境界条件と残留応力の影響
圧電FEMで計算した共振周波数が実測より5%程度高い場合、境界条件の設定ミスと残留応力の無視が主原因だ。「Free-Free（自由端）境界条件」を想定したFEMは試験台での「完全自由」状態（ワイヤだけで吊るした状態）と対応する。しかし実製品では接着剤・ネジ・はんだで固定されており、これらが付加剛性（Added Stiffness）として共振周波数を下げる効果がある。接着層の厚さ10〜20umでも弾性率が低い（エポキシのEr = 2〜5 GPa）接着剤は、PZTの弾性率（Er = 60〜70 GPa）に対して10〜30倍柔らかく、共振周波数を1〜3%シフトさせる。また焼成プロセスで生じる残留応力（圧縮〜引っ張り）が圧電係数d33を5〜10%変化させる効果も無視できない。実測との1%以内の一致を目標とするなら、製造プロセス条件を反映した詳細モデルが必要だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——圧電解析（静的）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告