首页/文章/ 详情

弱形式到底是什么?为什么有限元离不开它?

1天前浏览17

在学习有限元方法(FEM)时,我们总会遇到一个看似抽象又略显神秘的概念——微分方程的弱形式(Weak Form)

但教材往往“一上来就讲推导”,很多人没弄明白它到底是什么,就已经陷入了公式的旋涡中。

本篇文章,我们从头梳理弱形式的本质,力图回答三个关键问题:

  • 什么是“弱形式”?它和我们熟悉的微分方程形式有何不同?
  • 为什么有限元方法必须依赖弱形式才能成立?
  • 弱形式到底怎么推导出来?它的结构、意义和数值实现方式是什么?

什么是“弱形式”?

所谓弱形式,是对微分方程的一种积分等价重写
在“强形式”中,微分方程需要在每一个点都被严格满足,同时要求解函数具有较高的光滑性(例如二阶连续可导)。但在实际工程问题中,物理条件复杂,解往往不够光滑,这就给解析方法和数值求解带来了巨大挑战。

弱形式的出现,就是降低对解的光滑性(可微性)要求,它是许多强大的数值方法(尤其是有限元法 Finite Element Method, FEM)的理论基础。有限元方法的本质,就是将一个偏微分方程问题转化为一个代数系统进行求解。这个转化过程的核心步骤,正是将强形式变成弱形式。

强形式 vs 弱形式

为了说明“强形式”和“弱形式”的区别,我们从一个最简单的一维边值问题出发,这样更容易理解其背后的逻辑结构:

例子:一维泊松方程

考虑以下边值问题:

这是一个非常基础的微分方程,描述了例如一根受力杆的稳态变形问题,其中 u(x) 是位移, f(x) 是分布载荷。

什么是“强形式”?

强形式就是上面直接写出来的这个微分方程,它要求:

  1. 函数 u(x) 至少是二阶可导的(即 u ∈ C^2(0,1));
  2. 在定义域内的每一个点都严格满足该微分方程;
  3. 边界条件 u(0)=u(1)=0 也必须严格满足。

这样的要求非常“严格”,因此称为“强形式”。

如何得到“弱形式”?——四个步骤

我们现在来将这个强形式转换为弱形式,过程如下:

第一步:引入试函数(test function)

选择一个试函数 v(x) ∈ V,其中 V = { v ∈ H^1(0,1) ∣ v(0)=v(1)=0 },即满足同样边界条件的“光滑但不一定二阶可导”的函数。

第二步:乘以试函数并积分(加权残差法)

我们将方程两边都乘以 v(x),然后对整个定义域积分:

第三步:使用分部积分(积分转移导数)

我们不希望对 u'' 求值,因为太“强”。于是使用积分分部公式:

 

因为试函数 v(0) = v(1) = 0 ,边界项为零。

第四步:得到弱形式表达

最终我们得到弱形式:

 

🔍 强 vs 弱:真正的差别在哪里?

项目      
强形式      
弱形式      
函数要求      
u ∈ C^2      
u ∈ H^1(只需一阶导数)      
成立方式      
每一点都要满足      
对所有试函数 v 满足积分等式      
数学处理      
偏微分方程      
积分方程      
是否利于数值解      
     
是(可直接构造有限维子空间)      

弱形式更“宽容”,允许函数更粗糙;同时也更适合数值方法处理,如有限元、Galerkin方法等。

小结

这个例子展示了一个重要思想:

我们不是“简化”微分方程,而是“换个角度理解”它。
从点的约束 → 积分意义的约束,正是有限元能够成立的关键。


来源:锂电芯动

Comsolpython理论控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-06-10
最近编辑:1天前
锂电芯动
博士 中科院博士,电芯仿真高级工程师
获赞 2粉丝 0文章 38课程 0
点赞
收藏
作者推荐

如何测量极片的迂曲度

锂离子电池极片是多孔结构的。极片的孔隙率、孔径大小与分布以及迂曲度等微观结构参数是决定锂离子传输效率的关键因素。本文将重点介绍一下迂曲度的基本概念和常用的测量方法。极片迂曲度代表了多孔电极锂离子传输路径的曲折程度,即锂离子在电极涂层的实际传输路径∆L与涂层厚度∆x的比值(图1)。可以看出有效离子电导率与迂曲度成反比,因此需要设计低迂曲度的电极结构来提高电池的性能。图1. 多孔电极锂离子传输路径示意图 目前主要有以下几种方法可以用来测量或评估极片的迂曲度: 1.压汞法 这是一种传统的测量孔隙结构的方法,它的基本原理是利用液态汞的表面张力非常大,与各类固体材料的接触角都大于90度,因此不会润湿材料。需要通过施加外力才能使汞渗透到材料内部的孔隙中,根据压入汞的体积及样品的表观体积就可计算样品的孔隙率。压力越大,汞能够进入孔的半径越小。通过测量不同压力下进入孔中汞的量,可以得到相应孔大小的体积,从而得到多孔材料孔径分布的情况。结合孔径、孔径分布和孔隙率等数据可以计算出极片的迂曲度。需要注意的是由于汞无法进入极片的封闭孔中,所以它只能测量开口孔隙,即联通孔隙和半通孔。压汞法最小能测的孔径为2 nm,最大能测的孔径能到几百个微米。测试过程中,若加压过高时,还会破坏电极的涂层,导致出现裂痕从而产生新的孔隙。2.电化学阻抗谱(EIS)通过测量电池的电化学阻抗谱,可以得到电极的离子传输电阻,进而评估迂曲度。这种方法非破坏性,可以快速且准确地评估电极的迂曲度。此外,EIS作为一种无损技术,可以用于原位和操作研究,尤其在监控电池循环期间的状态方面发挥重要作用。这种无损性和实时性进一步增强了EIS在表征迂曲度方面的应用价值。利用电化学工作站测试对称电池的阻抗,得到的EIS如图2所示。此时电化学阻抗谱的Nyquist图具有低频区域线段和高频区域线段相交的形状特点,这是无电化学反应的典型Nyquist图。将Nyquist图中低频线段延长,直至与X轴相交,该交点与高频线段和X轴的交点的差值的3倍即为该极片涂层的离子阻抗Rion。图2. 对称电池的电化学阻抗谱将拟合得到的离子阻抗Rion代入下面的公式中计算可得到极片曲折度:A是电极横截面积,k是电解液电导率,ε是孔隙率,d是电极涂层厚度。3.X射线计算机断层扫描(XCT)XCT技术能够无损检测样品的内部结构, 具有高穿透性、友好的成像环境和丰富的衬度机制等特点, 能够获取大尺寸电极样品的三维微观结构,从而评估迂曲度。XCT测试电池极片微观结构时, 以X射线照射样品, 部分X射线被吸收, 透射的光束通过闪烁体探测器转换成可见光,经过放大处理最后在电荷耦合元件(CCD)图像传感器上形成可见光图像。旋转样品或者调节X射线强度, 形成了一系列投影图像。然后将系列图像通过计算机软件重构形成3D图像,测试原理示意图见图3。测试时为获取全面详细的统计学结构特征, 极片样品尺寸至少需要几百微米。随着XCT技术的进步, 空间分辨率也不断提高, 已经从微米级达到纳米级分辨率, 使用Zernike相位衬度也可以辨别碳胶相(导电碳和粘结剂)的分布。但设备成本较高,不是一般研发机构能够拥有的。图3. XCT测试原理图 4.聚焦离子束-扫描电子显微镜(FIB-SEM)FIB-SEM技术是一种高空间分辨破坏性成像技术。与XCT类似,也是通过对电极进行三维重构,从而评估迂极片曲度。但FIB-SEM技术重构的3D结构分辨率更高,一般为几十纳米,能够清晰地区分活性物质、导电剂和粘结剂的混合相以及孔洞。在使用该技术对电极进行3D研究中, 需要采用聚焦离子束(FIB)在电极基片上切出类似正方体形状,用SEM对正方体的侧面进行形貌扫描, 然后利用FIB把这个侧面切除几十纳米,再使用SEM进行形貌扫描,将反复切除和扫描成像后形成的系列图片进行3D重构, 如图4所示。但是由于样品需要层层切割, 样品信息采集耗时长, 测试成本高,而且样品尺寸不能太大, 一般20-30μm, 这种小的微观结构体积单元往往仅仅只包含了2-3个活性物质颗粒, 并不具备代表性, 无法对电极微观结构特征(如颗粒尺寸分布、比表面积、迂曲度等)进行统计分析。图4. FIB测试示意图 5.基于图像处理的迂曲度测量软件Bruggeman EstimatorEbner等基于微分等效介质理论, 开发了一个软件专门快速评估锂离子电池极片的孔隙迂曲度。只需输入电极表面和横截面两张扫描电子显微镜照片, BruggemanEstimator软件就可以通过分析每个活性物质颗粒在X, Y和Z三轴上的取向分布状态来计算迂曲度。下表展示的是该软件对球形NMC正极,片状石墨负极和菱形LCO正极的估算结果: 参考文献:WUXiangkun, ZHAN Qiushe, ZHANG Lan, et al. Progress on MicrostructuralOptimization and Controllable Preparation Technology for Lithium Ion BatteryElectrodes[J]. Chinese Journal of Applied Chemistry, 35(9): 1076-1092FroboeseL, Titscher P, Westphal B, et al. Mercury intrusion for ion and conversionbasedbattery electrodes-Structure and diffusion coefficient determination[J].Materials Characterization, 2017, 133.Landesfeind,J, et al. Tortuosity determination of battery electrodes and separators byimpedance spectroscopy. Journal of the Electrochemical Society, 2016,163(7): A1373-A1387.ChenWiegart, A, et al. Tortuosity characterization of 3D microstructure atnano-scale for energy storage and conversion materials - ScienceDirect. Journalof Power Sources, 2014, 249 (1):349-356. Dubeshter, T, et al. Measurement of Tortuosityand Porosity of Porous Battery Electrodes. Journal of the ElectrochemicalSociety, 2014, 161(4): A599-A605.EbnerM, Wood V. Tool for Tortuosity Estimation in Lithium Ion Battery PorousElectrodes[J]. Journal of the Electrochemical Society, 2014, 162(2): A3064-A3070汪晨阳,张安邦, 常增花, 等 .锂离子电池用多孔电极结构设计及制备技术进展[J].材料工程, 2022, 50 (1): 67-79.来源:锂电芯动

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈