電流・電圧・抵抗 — オームの法則からジュール発熱・CAE電磁解析まで
1. 電気とCAE — 現代設計に欠かせない連成解析
電気自動車（EV）・燃料電池車・産業用パワーエレクトロニクス・データセンターの冷却設計——現代の最重要技術領域はどれも電気と熱が密接に絡み合う。「電気解析だけ」「熱解析だけ」では解けない問題が急増しており、電気-熱連成解析（Electro-Thermal Coupling）のニーズが爆発的に高まっている。
2. 電荷・電流・電圧・電力の定義
電荷（Electric Charge） $Q$ [クーロン: C]：電気の量。電子1個の電荷量は $-e = -1.602\times10^{-19}$ C。
電流（Electric Current） $I$ [アンペア: A]：単位時間あたりに流れる電荷量。
$$I = \frac{dQ}{dt}$$電圧（Voltage）（電位差）$V$ [ボルト: V]：単位電荷を移動させるのに必要なエネルギー（電位のポテンシャル差）。
電力（Power） $P$ [ワット: W]：単位時間あたりのエネルギー変換量。
$$P = IV \quad \text{[W]}$$エネルギー $E = Pt = IVt$ [J]（充電量は Wh = 3600 J で表すことも多い）。
3. オームの法則と電気抵抗率
線形抵抗素子では電圧と電流が比例する——これがオームの法則（Ohm's Law）だ。
$$V = IR \quad \text{（} R \text{: 電気抵抗 [Ω]}）$$物体の抵抗値は形状と材料で決まる:
$$R = \frac{\rho_e L}{A}$$$\rho_e$: 電気抵抗率 [Ω·m]、$L$: 長さ [m]、$A$: 断面積 [m²]。
| 材料 | 電気抵抗率 ρ_e [Ω·m] at 20℃ | 主な用途 |
|---|---|---|
| 銀（Ag） | 1.59 × 10⁻⁸ | 高周波回路・接点 |
| 銅（Cu） | 1.72 × 10⁻⁸ | 電線・バスバー・コイル |
| アルミ（Al） | 2.82 × 10⁻⁸ | 送電線・バスバー（軽量化） |
| 炭素鋼 | 1.0 × 10⁻⁷ | 構造材（電気配線には使わない） |
| ステンレス（SUS304） | 7.2 × 10⁻⁷ | 発熱体・抵抗素子 |
| シリコン（半導体） | ~6.4 × 10² | トランジスタ・太陽電池 |
| CFRP（繊維方向） | ~1.5 × 10⁻⁵ | 航空機構造（落雷対策が必要） |
電気抵抗率の温度依存性
金属の抵抗率は温度上昇とともに増加する（格子振動による電子散乱）。
$$\rho_e(T) = \rho_0\left[1 + \alpha_{elec}(T - T_0)\right]$$銅の場合 $\alpha_{elec} \approx 3.9\times10^{-3}$ /K——つまり温度が100℃上昇すると抵抗が約39%増加する。この温度依存性がジュール発熱の自己増強効果（電流→発熱→抵抗増大→さらに発熱）を引き起こすため、熱解析と電気解析の連成が必須になる。
4. ジュール発熱
電流が抵抗を流れると電気エネルギーが熱エネルギーに変換される。これがジュール発熱（Joule Heating）だ。
$$P = I^2 R = \frac{V^2}{R} = IV \quad \text{[W]}$$連続体（FEM）での単位体積あたりのジュール発熱密度:
$$q_J = \rho_e J^2 = \vec{J}\cdot\vec{E} \quad \text{[W/m}^3\text{]}$$$J$ [A/m²]: 電流密度ベクトル、$E$ [V/m]: 電界ベクトル。
5. キルヒホッフの法則
複雑な電気回路を解析するための基本法則がキルヒホッフの法則（Kirchhoff's Laws）だ。
電流則（KCL: Kirchhoff's Current Law）
任意の節点（接続点）では、流入する電流の和 = 流出する電流の和。
$$\sum I_{in} = \sum I_{out} \quad \left(\text{等価: } \sum_k I_k = 0\right)$$電荷保存則（電荷が節点に蓄積しない）から導かれる。
電圧則（KVL: Kirchhoff's Voltage Law）
任意の閉ループ（閉回路）を一周したときの電圧の代数和はゼロ。
$$\sum_k V_k = 0 \quad \text{（閉ループ内）}$$エネルギー保存則から導かれる。
FEMネットワーク解析との類似性
回路解析の KCL+KVL と構造FEMの「力のつり合い+変位適合性」は数学的に全く同じ構造を持つ。
| 電気回路 | 構造FEM |
|---|---|
| 電圧 $V$ | 変位 $u$ |
| 電流 $I$ | 力 $F$ |
| 抵抗 $R$ | 剛性の逆数 $1/k$ |
| KCL（節点の電流和） | 力のつり合い（節点の力の和） |
| KVL（閉ループの電圧和） | 変位適合性（閉じた変形ループ） |
| $V = RI$ | $F = ku$ |
6. RC/RL回路の過渡応答
抵抗とキャパシタ（RC）または抵抗とインダクタ（RL）を組み合わせた回路では、スイッチをオンにした瞬間から「定常状態」に達するまで時間的な過渡現象が起きる。
RC回路の充電（ステップ応答）
$$v_C(t) = V_s\left(1 - e^{-t/\tau}\right), \quad \tau = RC$$$\tau$ [s] が時定数で、$t = \tau$ のとき $v_C = 0.632\,V_s$（最終値の63.2%）。$t = 5\tau$ でほぼ定常値（99.3%）に達する。
RL回路の電流立ち上がり
$$i(t) = \frac{V_s}{R}\left(1 - e^{-t/\tau}\right), \quad \tau = \frac{L}{R}$$インダクタは電流変化を妨げる（$v_L = L\,di/dt$）ため、スイッチ投入直後は電流がゼロで、徐々に $V_s/R$ に近づく。
CAEとの接続：電磁過渡解析
電磁解析（Transient Electromagnetic Analysis）では、RC/RL の1次応答が有限要素でも支配方程式の基礎になる。モーター始動時の過渡電流、変圧器のインラッシュ電流、雷サージの伝播解析などが対象だ。
7. 表皮効果（Skin Effect）
高周波電流は導体の表面に集中し、内部にはほとんど流れない。この現象を表皮効果（Skin Effect）と呼び、電流が内部に到達できる深さが表皮深さ（Skin Depth） $\delta$ で表される。
$$\delta = \sqrt{\frac{2\rho_e}{\omega\mu}} = \sqrt{\frac{\rho_e}{\pi f \mu}}$$$\omega = 2\pi f$ [rad/s]: 角周波数、$\mu$: 透磁率 [H/m]（真空中 $\mu_0 = 4\pi\times10^{-7}$）
銅ケーブルの表皮深さ（数値例）
| 周波数 | 銅の δ [mm] | 実用上の意味 |
|---|---|---|
| 60 Hz（商用電源） | 8.5 | ~17 mm 以上の丸線では内部が無駄 |
| 10 kHz | 0.66 | インバーターの高調波問題 |
| 1 MHz | 0.066 | RF回路・コイル設計に影響 |
| 1 GHz | 0.0021 | マイクロ波——表面めっきが支配的 |
モーター・変圧器コアでの渦電流損失
電磁鋼板（ケイ素鋼板）は積層構造にして各層を薄くすることで渦電流の流れる断面積を小さくし、渦電流損失を低減する。板厚を $t$ にすると渦電流損失は $t^2$ に比例して減少する。FEM渦電流解析（2D axisymmetric または3D）でコア損失分布を計算してモーター効率を最適化する。
8. FEM電磁解析への架け橋
電磁現象はMaxwell方程式によって記述される。FEMではこれを弱形式化して解く。
静電場（直流電気解析）
静電位 $V$（電圧）のポアソン方程式:
$$\nabla^2 V = -\frac{\rho_{charge}}{\varepsilon}$$導体内の定常電流解析では $\rho_{charge} = 0$（電荷分布なし）でラプラス方程式 $\nabla^2 V = 0$ になる。
定電流解析（DC Conduction Analysis）
$$\nabla\cdot\vec{J} = 0, \quad \vec{J} = \sigma_e\vec{E} = -\sigma_e\nabla V$$$\sigma_e = 1/\rho_e$: 電気伝導率 [S/m]。FEMでは電気伝導率を材料プロパティとして入力し、電位 $V$ を未知数として解く。
渦電流解析（Eddy Current Analysis）
時間変化する磁場による誘導電流（渦電流）はMaxwell方程式から:
$$\nabla\times\vec{H} = \vec{J} + \frac{\partial\vec{D}}{\partial t}, \quad \nabla\times\vec{E} = -\frac{\partial\vec{B}}{\partial t}$$FEMでは磁気ベクトルポテンシャル $\vec{A}$（$\vec{B} = \nabla\times\vec{A}$）を未知数として変位電流項を無視した近似で解く。Ansys Maxwell、COMSOL Multiphysics、GetDP（オープンソース）が主要ソルバーだ。
9. 実践：電気自動車バスバーの最適化
EV（電気自動車）のバッテリーパックとモーターインバーターを接続するバスバー（Bus Bar）は、数百アンペアの直流電流を流す高電流導体だ。銅またはアルミ製の平板状導体で、電流密度の不均一分布とジュール発熱が設計の核心問題になる。
問題設定（典型例）
- 電流: 500 A（急速充電時）
- 材料: 銅（$\rho_e = 1.72\times10^{-8}$ Ω·m, $\lambda = 400$ W/mK）
- 最大温度: 90℃ 以下（インバーター要求）
- 目標: 電流密度均一化 + 重量最小化
電気-熱連成FEM解析フロー
- 電気解析：入力端・出力端の電位境界条件で電位分布→電流密度分布→ジュール発熱密度 $q_J = \rho_e J^2$ を計算
- 熱解析：$q_J$ を熱源として温度分布を計算（対流冷却境界条件: $h = 10〜20\,\text{W/m}^2\text{K}$ の空冷）
- 反復収束：温度依存の電気抵抗率 $\rho_e(T)$ で両解析を収束するまで繰り返す
電流密度均一化の設計手法
電流は最短経路（角の内側）に集中するため、バスバーの曲がり部では内側に電流が集中して局所過熱しやすい。
- 角部の内側を切り欠いて電流集中を緩和（電流均一化）
- トポロジー最適化で最適な穴パターンを自動探索
- 積層バスバー（複数枚の薄板を絶縁層で積層）でインダクタンス低減と電流均一化
10. 圧電効果・電歪
圧電材料（PZT, 水晶など）では機械的応力と電気的変位が相互に変換される。CAEでは電気-構造連成解析の対象だ。
直接圧電効果（センサー：応力→電圧）
$$D_i = d_{ijk}\sigma_{jk} + \varepsilon_{ij}^T E_j$$力が加わると電荷が発生する。スマートフォンの加速度センサー、超音波探傷子（UT）、自動車エンジンのノックセンサーに使用。
逆圧電効果（アクチュエータ：電圧→ひずみ）
$$\varepsilon_{ij} = d_{kij}E_k + S_{ijkl}^E\sigma_{kl}$$電圧印加で微小変位を発生させる。超音波発生素子（医療用エコー）、インクジェットプリンターのノズル、半導体製造用ナノポジショナーに使用。
FEM圧電解析
圧電FEM要素では電位と変位が連成自由度となる。Abaqusの PIEZOELECTRIC 材料、COMSOL の「Structural Mechanics + Electrostatics」モジュールが代表的な実装だ。