电流、电压、电阻 — 从欧姆定律到焦耳热、CAE电磁解析
1. 电气与CAE — 现代设计中不可或缺的耦合解析
电动汽车（EV）、燃料电池车、工业电力电子设备、数据中心冷却设计——现代最重要的技术领域中，电气和热密切相关。"仅进行电气解析"或"仅进行热解析"已经无法解决新兴问题，电气-热耦合解析（Electro-Thermal Coupling）的需求呈爆发式增长。
2. 电荷、电流、电压、电力的定义
电荷（Electric Charge） $Q$ [库仑：C]：电的数量。单个电子的电荷量为 $-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]（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的一阶响应成为有限元支配方程的基础。电动机启动时的过渡电流、变压器的励磁涌流、雷击浪涌传播解析等都是其应用对象。
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 | 射频电路、线圈设计受影响 |
| 1 GHz | 0.0021 | 微波——表面镀层主导 |
电动机、变压器铁芯的涡流损耗
电磁钢板（硅钢板）采用积层结构，每层变薄来减少涡电流的流动截面积，从而降低涡电流损耗。将板厚设为 $t$ 时，涡电流损耗与 $t^2$ 成正比并减小。通过FEM涡电流解析（2D轴对称或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的"结构力学+静电学"模块是典型实现。