圧電解析（静的）
圧電解析（静的）の理論基礎
圧電効果
先生、圧電効果って機械的な力を電圧に変換することですか？
正圧電効果（力→電圧）と逆圧電効果（電圧→変形）の両方がある。
構成方程式:
$\{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が圧電クォーツを使った水中探信儀を開発した。現代のスマートフォンのカメラ手振れ補正アクチュエータ、医療用超音波診断装置、インクジェットプリンタのヘッドはすべてキュリー兄弟の発見の子孫だ。
圧電解析（静的）の数値計算手法
圧電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モジュールはこの変換を自動処理するが、ユーザー定義材料では変換行列を手動で設定する必要がある。
圧電解析（静的）の実務適用
実務
超音波トランスデューサ、圧電アクチュエータ（インクジェットヘッド）、振動センサ、エナジーハーベスタの設計。
チェックリスト
超音波洗浄機の圧電トランスデューサ設計——キャビテーション分布とFEM最適化
超音波洗浄機は圧電トランスデューサが発生する超音波振動で液中にキャビテーション気泡を生成し、その崩壊衝撃で汚れを除去する。トランスデューサの共振周波数設計（20〜40kHz）と洗浄槽内の音圧分布が洗浄均一性を決定する。FEM（圧電連成）でトランスデューサの振動モードと音響場（Helmholtz方程式）を解析し、音圧が均一になる配置と形状を最適化する。実際の設計事例では、4個の円板型トランスデューサの配置を非対称にすることでキャビテーション分布の均一性を30%改善し、精密洗浄ラインでの洗浄不良率を0.8%から0.1%以下に低減した。FEMによる設計最適化が半導体ウェーハ洗浄・医療機器洗浄の品質管理に直結している。
圧電解析（静的）のソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| 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コードが圧電研究には使われる。
圧電解析（静的）の先端研究
先端
圧電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%以内の一致を目標とするなら、製造プロセス条件を反映した詳細モデルが必要だ。
なった
詳しく
報告