1. 项目概述：电容器FEM仿真的核心价值

电容器作为电子电路中的基础元件，其内部电场分布直接影响着器件性能。传统解析法在处理复杂边界条件时存在明显局限，而有限元方法（FEM）通过区域离散化和插值函数构建，能够精确求解任意形状电容器的电位分布问题。这个项目将展示如何用Matlab实现二维轴对称电容器模型的完整FEM仿真流程，包含从网格划分到后处理可视化的全链条操作。

在电力电子、传感器设计等领域，工程师经常需要评估不同结构电容器的场强分布、储能密度等关键参数。通过自主编写FEM代码（而非依赖现成商业软件），我们可以深度掌控算法细节，方便后续定制化修改。本文代码已实测通过Matlab R2021a及以上版本，完整工程文件可在文末获取。

2. 有限元方法理论基础

2.1 泊松方程与变分原理

电容器内部静电场满足泊松方程： ∇²φ = -ρ/ε 其中φ为电位，ρ为电荷密度，ε为介电常数。通过伽辽金加权残差法将其转化为弱形式：

∫Ω (∇w·∇φ) dΩ = ∫Γ w (∂φ/∂n) dΓ + ∫Ω w (ρ/ε) dΩ

式中w为权函数，Γ为边界区域。采用三角形线性单元离散后，可得到全局刚度矩阵K和载荷向量F，最终形成线性方程组： KΦ = F

2.2 轴对称模型简化

对于圆柱形电容器，采用柱坐标(r,z)下的轴对称建模可将三维问题降维为二维处理。此时梯度算子变为： ∇ = (∂/∂r, ∂/∂z)

相应的单元刚度矩阵需要进行修正，包含r的积分权重。这种处理既能保证计算精度，又可大幅减少网格数量。

3. Matlab实现详解

3.1 前处理：几何建模与网格生成

% 定义电容器几何参数 inner_radius = 0.01; % 内电极半径(m) outer_radius = 0.02; % 外电极半径 height = 0.03; % 电容器高度 dielectric_const = 4; % 介电常数 % 使用PDETool创建几何模型 model = createpde(); g = decsg([3 4 0 outer_radius outer_radius 0 0 0 height height]'); geometryFromEdges(model,g);

提示：对于复杂几何，建议先用CAD软件绘制后导入。Matlab的importGeometry函数支持STEP格式文件。

3.2 有限元矩阵组装

核心刚度矩阵计算代码如下：

function [K,F] = assemble_system(nodes,elements,epsilon) n_nodes = size(nodes,1); K = sparse(n_nodes,n_nodes); F = zeros(n_nodes,1); % 高斯积分点与权重 gauss_points = [1/6 1/6; 2/3 1/6; 1/6 2/3]; weights = [1/6 1/6 1/6]; for el = 1:size(elements,1) idx = elements(el,:); verts = nodes(idx,:); % 计算雅可比矩阵 J = [verts(2,1)-verts(1,1) verts(3,1)-verts(1,1); verts(2,2)-verts(1,2) verts(3,2)-verts(1,2)]; detJ = abs(det(J)); % 形函数导数 dN = [-1 -1; 1 0; 0 1] / J; % 单元刚度矩阵 Ke = zeros(3,3); for gp = 1:3 r = sum(verts(:,1).*[1-gauss_points(gp,1)-gauss_points(gp,2); gauss_points(gp,1); gauss_points(gp,2)]); Ke = Ke + (dN'*dN) * (weights(gp)*detJ*r); end K(idx,idx) = K(idx,idx) + epsilon*Ke; end end

3.3 边界条件处理

典型平行板电容器边界设置：

% 内电极设为1V inner_edge = find(abs(nodes(:,1)-inner_radius)<1e-6); fixed_dofs = inner_edge; fixed_values = ones(size(inner_edge)); % 外电极接地 outer_edge = find(abs(nodes(:,1)-outer_radius)<1e-6); fixed_dofs = [fixed_dofs; outer_edge]; fixed_values = [fixed_values; zeros(size(outer_edge))]; % 处理Dirichlet边界条件 free_dofs = setdiff(1:size(nodes,1), fixed_dofs); K_ff = K(free_dofs,free_dofs); F_f = F(free_dofs) - K(free_dofs,fixed_dofs)*fixed_values;

4. 后处理与结果分析

4.1 电位分布可视化

phi = zeros(size(nodes,1),1); phi(free_dofs) = K_ff\F_f; phi(fixed_dofs) = fixed_values; trisurf(elements,nodes(:,1),nodes(:,2),phi); xlabel('径向距离(m)'); ylabel('轴向距离(m)'); zlabel('电位(V)'); title('电容器电位分布'); colormap jet; shading interp;

4.2 电场强度计算

[Ex,Ey] = pdegrad(nodes,elements,phi); E_mag = sqrt(Ex.^2 + Ey.^2); % 绘制电场矢量图 pdeplot(model,'XYData',phi,'FlowData',[Ex;Ey]); title('电场强度分布');

5. 性能优化技巧

5.1 稀疏矩阵处理

% 创建稀疏模式加速组装 sparsity_pattern = sparse(n_nodes,n_nodes); for el = 1:size(elements,1) idx = elements(el,:); sparsity_pattern(idx,idx) = 1; end K = spalloc(n_nodes,n_nodes,nnz(sparsity_pattern));

5.2 自适应网格加密

基于电场梯度实现局部网格细化：

error_indicator = sqrt(Ex.^2 + Ey.^2); [new_nodes,new_elements] = refinemesh(... nodes,elements,error_indicator>0.7*max(error_indicator));

6. 常见问题排查

1. 矩阵奇异警告

   • 检查边界条件是否施加充分
   • 确认网格中无重复节点
   • 验证材料参数是否为零值

2. 电场计算结果异常

   • 检查单位制统一性（米/毫米）
   • 确认介电常数输入正确
   • 验证高斯积分阶数是否足够

3. 内存不足错误

   • 改用稀疏存储格式
   • 采用分块组装策略
   • 对于大型模型考虑使用PCG迭代求解器

7. 工程应用扩展

7.1 多层介质电容器

通过修改单元属性矩阵实现：

epsilon_array = ones(size(elements,1),1)*epsilon0; epsilon_array(elements_region2) = epsilon0*4;

7.2 瞬态场计算

引入时间离散项：

C = assemble_mass_matrix(nodes,elements); % 质量矩阵 [M,K,F] = transient_solver(C,K,F,dt,total_time);

完整项目代码包含：

• 基础FEM求解器（fem_solver.m）
• 后处理工具包（postprocess.m）
• 示例模型库（parallel_plate.m, cylindrical.m）
• 可视化脚本（plot_results.m）