搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于电化学-应力耦合模型的锂离子电池硅/碳核壳结构的模拟与优化

柳小伟 宋辉 郭美卿 王根伟 迟青卓

引用本文:
Citation:

基于电化学-应力耦合模型的锂离子电池硅/碳核壳结构的模拟与优化

柳小伟, 宋辉, 郭美卿, 王根伟, 迟青卓

Simulation and optimization of silicon/carbon core-shell structures in lithium-ion batteries based on electrochemical-mechanical coupling model

Liu Xiao-Wei, Song Hui, Guo Mei-Qing, Wang Gen-Wei, Chi Qing-Zhuo
PDF
HTML
导出引用
  • 硅基电极材料在应用中的一个主要问题是巨大的体积膨胀, 以及由此带来的电极材料破裂、粉化. 本文在有限变形假设前提下, 基于电化学-力学耦合理论, 研究球形Si/C核壳结构在嵌锂过程中的浓度、应力场的演化, 并在此基础上讨论了核壳结构的优化设计. 计算结果显示: 壳层可以很好地保护硅颗粒的膨胀; 然而核内产生的较大的径向压缩应力可能导致核壳界面的脱黏, 而核壳界面处的切向拉伸应力可能会导致壳层的断裂. 进一步为有效提高核壳结构的电化学与力学性能, 从而实现锂离子电池更长的循环寿命, 考虑了两种结构的优化: 1)单层核壳结构; 2)双层核壳结构. 结果表明对于单层核壳结构应使用更软的包覆层材料; 而双层核壳结构中优化的材料布置方案为内软外硬, 对双层核壳结构的硬度分析表明, 内层材料的杨氏模量应低于10 GPa, 而外层材料的应不高于70 GPa. 本文的结论对球形材料颗粒电极的设计及优化具有一定的指导意义.
    Silicon is considered as a first candidate for ideal anode material of the next-generation lithium-ion battery due to its high theoretical capacity to meet the demand for higher energy density. On the other hand, high theoretical capacity is accompanied by massive volume expansion, which gives arise to high stress and crack and pulverization of anode particles. Finally, the capacity of the battery fades gradually. While some kinds of factors contribute to the failure of silicon-based electrodes, the most important one is the diffusion-induced stress generated in silicon-based electrode particles. The cyclic processes of lithiation and delithiation are accomplished by the intercalation into and deintercalation from the silicon particles of lithium ions. During the cycle, physical processes and chemical processes, such as diffusion of lithium ions, phase transition, and volume expansion, take place simultaneously, making the cyclic process a strong-coupling problem to be addressed. For example, the intercalation of lithium ions into the electrode results in volume expansion and phase transition of anodes, thereby inducing stress; in turn, stress affects the diffusion process of lithium ions. Aiming to probe this problem, with the finite deformation hypothesis, an electrochemical-mechanical coupling model is used to study the variation and distribution of concentration and stress of core-shell structure during lithiation. And more importantly, great emphasis is put on the optimal design of core-shell structure. The numerical results show that the shell is useful in prohibiting the volume expansion of silicon core, but large compressive radial stress in silicon core may cause the core and shell to be detached, while the tangential tensile stress at the core-shell interface leads the shell to fracture. To improve the electrochemical and mechanical performance and hence lengthen the cycle life of lithium-ion batteries, two kinds of optimal designs are considered: 1) single-layered core-shell structure and 2) double-layered core-shell structure. The numerical results suggest that the softer shell material is suitable for a single-layered core-shell structure and the inner-soft & outer-hard design is optimal for the double-layered core-shell structure. Furthermore, the effects of Young's modulus of the inner and outer carbon layer materials on the chemical and mechanical performance of anode are explored. The simulation shows that the optimal Young's modulus of the inner shell is less than 10 GPa, and that of the outer shell is not higher than 70 GPa. This research is helpful in designing and optimizing the silicon-based anode electrodes of lithium-ion batteries.
      通信作者: 宋辉, songhui@tyut.edu.cn ; 郭美卿, guomeiqing@tyut.edu.cn
    • 基金项目: 国家自然科学基金 (批准号: 11872265)和山西省自然科学基金 (批准号: 201901D111087, 201801D121281)资助的课题
      Corresponding author: Song Hui, songhui@tyut.edu.cn ; Guo Mei-Qing, guomeiqing@tyut.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 11872265) and the Natural Science Foundation of Shanxi Province of China (Grant Nos. 201901D111087, 201801D121281)
    [1]

    程昀, 李劼, 贾明, 汤依伟, 杜双龙, 艾立华, 殷宝华, 艾亮 2015 物理学报 64 210202Google Scholar

    Cheng Y, Li J, Jia M, Tang Y W, Du S L, Ai L H, Yin B H, Ai L 2015 Acta Phys. Sin. 64 210202Google Scholar

    [2]

    Li M, Lu J, Chen Z, Khalil A 2018 Adv. Mater. 30 1800561Google Scholar

    [3]

    Xie Z, Ma Z, Wang Y, Zhou Y, Lu C 2016 RSC Adv. 6 22383Google Scholar

    [4]

    Vetter J, Novák P, Wagner M R, Veit C, Möller K C, Besenhard J O, Winter M, Wohlfahrt-Mehrens M, Vogler C, Hammouche A 2005 J. Power Sources 147 269Google Scholar

    [5]

    孙凤楠, 冯露, 卜家贺, 张静, 李林安, 王世斌 2019 物理学报 68 120201Google Scholar

    Sun F N, Feng L, Bu J H, Zhang J, Li L A, Wang S B 2019 Acta Phys. Sin. 68 120201Google Scholar

    [6]

    Baggetto L, Niessen R A H, Roozeboom F, Notten P H L 2008 Adv. Funct. Mater. 18 1057Google Scholar

    [7]

    Panat R 2015 Thin Solid Films 596 174Google Scholar

    [8]

    Zhao Y, Stein P, Bai Y, Al-Siraj M, Yang Y, Xu B X 2019 J. Power Sources 413 259Google Scholar

    [9]

    Prussin S 1961 J. Appl. Phys. 32 1876Google Scholar

    [10]

    Zhang X, Shyy W, Sastry A M 2007 J. Electrochem. Soc. 154 A910Google Scholar

    [11]

    Song Y, Lu B, Ji X, Zhang J 2012 J. Electrochem. Soc. 159 A2060Google Scholar

    [12]

    Bagheri A, Arghavani J, Naghdabadi R 2019 Mech. Mater. 137 103134Google Scholar

    [13]

    Lu Y, Zhang P, Wang F, Zhang K, Zhao X 2018 Electrochim. Acta 274 359Google Scholar

    [14]

    Christensen J 2010 J. Electrochem. Soc. 157 A366Google Scholar

    [15]

    Zhang X Y, Hao F, Chen H S, Fang D N 2014 J. Electrochem. Soc. 161 A2243Google Scholar

    [16]

    Cheng Y T, Verbrugge M W 2008 J. Appl. Phys. 104 1876

    [17]

    Deluca C M, Maute K, Dunn M L 2011 J. Power Sources 196 9672Google Scholar

    [18]

    Li Y, Zhang K, Zheng B, Zhang X Q, Wang Q 2015 J. Appl. Phys. 117 245103

    [19]

    Hu B, Ma Z S, Lei W X, Zou Y L, Lu C S 2017 Theor. Appl. Mech. Lett. 7 199Google Scholar

    [20]

    Zhao Y F, Lu B, Zhang J 2018 Acta. Mech. Solida. Sin. 31 290Google Scholar

    [21]

    Gao X, Lu W, Xu J 2020 J. Power Sources 449 227501Google Scholar

    [22]

    Hao F, Fang D 2013 J. Electrochem. Soc. 160 A595Google Scholar

    [23]

    Zhao K J, Pharr M, Hartle L, Vlassak J J, Suo Z G 2012 J. Power Sources 218 6Google Scholar

    [24]

    Wu B, Lu W 2017 J. Phys. Chem. C 121 19022Google Scholar

    [25]

    Zhang Y, Zhu Z 2020 Int. J. Adv. Manuf. Technol. 108 499Google Scholar

    [26]

    Bohn E, Eckl T, Kamlah M, McMeekingc R 2013 J. Electrochem. Soc. 160 A1638Google Scholar

    [27]

    Zhang X, Chen H S, Fang D 2019 Int. J. Mech. Sci. 169 105323

    [28]

    Xu C, Weng L, Chen B, Zhou J Q, Cai R 2019 Int. J. Mech. Sci. 157 87

    [29]

    Zhao K, Pharr M, Cai S, Vlassak J J, Suo Z G 2011 J. Am. Ceram. Soc. 94 s226Google Scholar

    [30]

    Zhang K, Zheng B, Yang F, Li Y 2020 Int. J. Mech. Sci. 177 105602Google Scholar

    [31]

    Liu B H, Wang X, Chen H S, Chen S, Yang H X, Xu J, Jiang H Q, Fang D N 2019 J. Appl. Mech. 86 041005Google Scholar

    [32]

    Ding N, Xu J, Yao Y X, Wegner G, Fang X, . Chen C.H, Lieberwirth I 2009 Solid State Ionics 180 222Google Scholar

    [33]

    Sethuraman V A, Chon M J, Shimshak M, Winkle N V, Guduru P R 2010 Electrochem. Commun. 12 1614Google Scholar

    [34]

    Wu B, Lu W 2017 J. Power Sources 360 360Google Scholar

    [35]

    Deng Q, Hu R, Xu C, Chen B, Zhou J 2019 J. Solid State Electrochem. 23 2999Google Scholar

  • 图 1  锂离子电池的用途

    Fig. 1.  Usage of lithium-ion batteries.

    图 2  电化学-力学耦合机理

    Fig. 2.  Electrochemistry-mechanics coupling mechanism.

    图 3  (a) 锂离子电池核壳结构示意图; (b) 径向和切向示意图

    Fig. 3.  (a) Schematic diagram of core-shell structure of lithium-ion battery; (b) radial and tangential schematic diagram.

    图 4  不同假设条件下充电过程中沿着半径方向的径向应力和切向应力 (a)小变形下的径向应力; (b)有限变形下的径向应力; (c)小变形下的切向应力; (d)有限变形下的切向应力

    Fig. 4.  Radial stress and tangential stress along the radial direction during charging under different hypothesis: (a) Radial stress under small deformation; (b) radial stress under finite deformation; (c) tangential stress under small deformation; (d) tangential stress under finite deformation.

    图 5  嵌锂过程中沿着半径方向的(a)浓度和(b) von Mises 应力分布示意图

    Fig. 5.  Schematic diagram of (a) concentration; (b) von Mises stress along the radius during lithium insertion.

    图 6  嵌锂过程中(a)浓度, (b) von Mises应力, (c)径向应力, (d)切向应力随着时间的分布示意

    Fig. 6.  Schematic diagram of the distribution of (a) concentration, (b) von Mises stress, (c) radial stress, and (d) tangential stress over time during lithium insertion.

    图 7  双层核壳结构在嵌锂过程(a)浓度, (b) von Mises应力, (c)径向应力, (d)切向应力的分布示意图

    Fig. 7.  Schematic diagram of the distribution of the double-layer core-shell structure during lithium insertion process: (a) concentration; (b) von Mises stress; (c) radial stress; and (d) tangential stress.

    图 8  Si/C1/C3双层核壳结构和Si/C3/C1双层核壳结构的(a)浓度, (b) von Mises应力, (c)径向应力, (d)切向应力分布示意图

    Fig. 8.  (a) Concentration, (b) von Mises stress, (c) radial stress, and (d) tangential stress distribution schematic diagram of Si/C1/C3 double-layer core-shell structure and Si/C3/C1 double-layer core-shell structure.

    图 9  不同的内外层壳杨氏模量情况下 (a)浓度, (b)浓度, (c) von Mises应力, (d)径向应力分布示意图

    Fig. 9.  (a) Concentration, (b) concentration, (c) von Mises stress, (d) schematic diagram of radial stress distribution under different Young's modulus of inner and outer shells.

    表 1  核壳结构材料参数[16,32,33]

    Table 1.  Material parameters of core-shell structure[16,32,33]

    参数
    饱和浓度 Cmax/(mol·m–3)2.95 × 1052.4 × 104
    初始浓度C0/(mol·m–3)00
    扩散系数 D/(m2·s–1)1.0 × 10–161.45 × 10–13
    内径/外径 R/m4 × 10–85 × 10–8
    泊松比 ν0.230.30
    杨氏模量E/Pa8 × 10106 × 1010
    偏摩尔体积 Ω/(m3·mol–1)3/Cmax3.497 × 10–6
    室温 T/K298
    下载: 导出CSV

    表 2  包覆层碳的材料参数

    Table 2.  Material parameters of the outer shell.

    参数杨氏模量E/Pa饱和浓度Cmax/(mol·m–3)偏摩尔体积Ω/(m3·mol–1)泊松比ν扩散系数D/(m2·s–1)
    碳层[34]1 × 1010241613.497 × 10–60.31.0 × 10–14
    外包覆层[35]3 × 10102.5 × 1043.497 × 10–60.31.0 × 10–14
    下载: 导出CSV
  • [1]

    程昀, 李劼, 贾明, 汤依伟, 杜双龙, 艾立华, 殷宝华, 艾亮 2015 物理学报 64 210202Google Scholar

    Cheng Y, Li J, Jia M, Tang Y W, Du S L, Ai L H, Yin B H, Ai L 2015 Acta Phys. Sin. 64 210202Google Scholar

    [2]

    Li M, Lu J, Chen Z, Khalil A 2018 Adv. Mater. 30 1800561Google Scholar

    [3]

    Xie Z, Ma Z, Wang Y, Zhou Y, Lu C 2016 RSC Adv. 6 22383Google Scholar

    [4]

    Vetter J, Novák P, Wagner M R, Veit C, Möller K C, Besenhard J O, Winter M, Wohlfahrt-Mehrens M, Vogler C, Hammouche A 2005 J. Power Sources 147 269Google Scholar

    [5]

    孙凤楠, 冯露, 卜家贺, 张静, 李林安, 王世斌 2019 物理学报 68 120201Google Scholar

    Sun F N, Feng L, Bu J H, Zhang J, Li L A, Wang S B 2019 Acta Phys. Sin. 68 120201Google Scholar

    [6]

    Baggetto L, Niessen R A H, Roozeboom F, Notten P H L 2008 Adv. Funct. Mater. 18 1057Google Scholar

    [7]

    Panat R 2015 Thin Solid Films 596 174Google Scholar

    [8]

    Zhao Y, Stein P, Bai Y, Al-Siraj M, Yang Y, Xu B X 2019 J. Power Sources 413 259Google Scholar

    [9]

    Prussin S 1961 J. Appl. Phys. 32 1876Google Scholar

    [10]

    Zhang X, Shyy W, Sastry A M 2007 J. Electrochem. Soc. 154 A910Google Scholar

    [11]

    Song Y, Lu B, Ji X, Zhang J 2012 J. Electrochem. Soc. 159 A2060Google Scholar

    [12]

    Bagheri A, Arghavani J, Naghdabadi R 2019 Mech. Mater. 137 103134Google Scholar

    [13]

    Lu Y, Zhang P, Wang F, Zhang K, Zhao X 2018 Electrochim. Acta 274 359Google Scholar

    [14]

    Christensen J 2010 J. Electrochem. Soc. 157 A366Google Scholar

    [15]

    Zhang X Y, Hao F, Chen H S, Fang D N 2014 J. Electrochem. Soc. 161 A2243Google Scholar

    [16]

    Cheng Y T, Verbrugge M W 2008 J. Appl. Phys. 104 1876

    [17]

    Deluca C M, Maute K, Dunn M L 2011 J. Power Sources 196 9672Google Scholar

    [18]

    Li Y, Zhang K, Zheng B, Zhang X Q, Wang Q 2015 J. Appl. Phys. 117 245103

    [19]

    Hu B, Ma Z S, Lei W X, Zou Y L, Lu C S 2017 Theor. Appl. Mech. Lett. 7 199Google Scholar

    [20]

    Zhao Y F, Lu B, Zhang J 2018 Acta. Mech. Solida. Sin. 31 290Google Scholar

    [21]

    Gao X, Lu W, Xu J 2020 J. Power Sources 449 227501Google Scholar

    [22]

    Hao F, Fang D 2013 J. Electrochem. Soc. 160 A595Google Scholar

    [23]

    Zhao K J, Pharr M, Hartle L, Vlassak J J, Suo Z G 2012 J. Power Sources 218 6Google Scholar

    [24]

    Wu B, Lu W 2017 J. Phys. Chem. C 121 19022Google Scholar

    [25]

    Zhang Y, Zhu Z 2020 Int. J. Adv. Manuf. Technol. 108 499Google Scholar

    [26]

    Bohn E, Eckl T, Kamlah M, McMeekingc R 2013 J. Electrochem. Soc. 160 A1638Google Scholar

    [27]

    Zhang X, Chen H S, Fang D 2019 Int. J. Mech. Sci. 169 105323

    [28]

    Xu C, Weng L, Chen B, Zhou J Q, Cai R 2019 Int. J. Mech. Sci. 157 87

    [29]

    Zhao K, Pharr M, Cai S, Vlassak J J, Suo Z G 2011 J. Am. Ceram. Soc. 94 s226Google Scholar

    [30]

    Zhang K, Zheng B, Yang F, Li Y 2020 Int. J. Mech. Sci. 177 105602Google Scholar

    [31]

    Liu B H, Wang X, Chen H S, Chen S, Yang H X, Xu J, Jiang H Q, Fang D N 2019 J. Appl. Mech. 86 041005Google Scholar

    [32]

    Ding N, Xu J, Yao Y X, Wegner G, Fang X, . Chen C.H, Lieberwirth I 2009 Solid State Ionics 180 222Google Scholar

    [33]

    Sethuraman V A, Chon M J, Shimshak M, Winkle N V, Guduru P R 2010 Electrochem. Commun. 12 1614Google Scholar

    [34]

    Wu B, Lu W 2017 J. Power Sources 360 360Google Scholar

    [35]

    Deng Q, Hu R, Xu C, Chen B, Zhou J 2019 J. Solid State Electrochem. 23 2999Google Scholar

  • [1] 李晓杰, 喻云泰, 张志文, 董小瑞. 基于电化学老化衰退模型的锂离子动力电池外特性. 物理学报, 2022, 71(3): 038803. doi: 10.7498/aps.71.20211401
    [2] 董军, 张晨雪, 程小同, 邢宇, 韩庆艳, 严学文, 祁建霞, 刘继红, 杨祎, 高伟. 构建NaYF4:Yb3+/Ho3+/Ce3+@NaYF4:Yb3+/Nd3+纳米核壳结构增强Ho3+离子的上转换红光发射. 物理学报, 2021, 70(15): 154208. doi: 10.7498/aps.70.20210118
    [3] 洪文鹏, 兰景瑞, 李浩然, 李博宇, 牛晓娟, 李艳. 基于时域有限差分法的核壳双金属纳米颗粒光吸收率反转行为. 物理学报, 2021, 70(20): 207801. doi: 10.7498/aps.70.20210602
    [4] 李涛, 程夕明, 胡晨华. 锂离子电池电化学降阶模型性能对比. 物理学报, 2021, 70(13): 138801. doi: 10.7498/aps.70.20201894
    [5] 张佳晨, 鱼卫星, 肖发俊, 赵建林. 金薄膜衬底上介质-金属核壳结构的光学力调控. 物理学报, 2020, 69(18): 184206. doi: 10.7498/aps.69.20200214
    [6] 严学文, 王朝晋, 王博扬, 孙泽煜, 张晨雪, 韩庆艳, 祁建霞, 董军, 高伟. 构建核壳结构增强Ho3+离子在镥基纳米晶中的红光上转换发射. 物理学报, 2019, 68(17): 174204. doi: 10.7498/aps.68.20190441
    [7] 曾建邦, 郭雪莹, 刘立超, 沈祖英, 单丰武, 罗玉峰. 基于电化学-热耦合模型研究隔膜孔隙结构对锂离子电池性能的影响机制. 物理学报, 2019, 68(1): 018201. doi: 10.7498/aps.68.20181726
    [8] 庞辉. 基于扩展单粒子模型的锂离子电池参数识别策略. 物理学报, 2018, 67(5): 058201. doi: 10.7498/aps.67.20172171
    [9] 宋旭, 陆勇俊, 石明亮, 赵翔, 王峰会. 集流体塑性变形对锂离子电池双层电极中锂扩散和应力的影响. 物理学报, 2018, 67(14): 140201. doi: 10.7498/aps.67.20180148
    [10] 庞辉. 基于电化学模型的锂离子电池多尺度建模及其简化方法. 物理学报, 2017, 66(23): 238801. doi: 10.7498/aps.66.238801
    [11] 林莹莹, 李葵英, 单青松, 尹华, 朱瑞苹. ZnSe/ZnS/L-Cys核壳结构量子点光声与表面光伏特性. 物理学报, 2016, 65(3): 038101. doi: 10.7498/aps.65.038101
    [12] 马昊, 刘磊, 路雪森, 刘素平, 师建英. 锂离子电池正极材料Li2FeSiO4的电子结构与输运特性. 物理学报, 2015, 64(24): 248201. doi: 10.7498/aps.64.248201
    [13] 李娟, 汝强, 胡社军, 郭凌云. 锂离子电池SnSb/C复合负极材料的热碳还原法制备及电化学性能研究. 物理学报, 2014, 63(16): 168201. doi: 10.7498/aps.63.168201
    [14] 李娟, 汝强, 孙大伟, 张贝贝, 胡社军, 侯贤华. 锂离子电池SnSb/MCMB核壳结构负极材料嵌锂性能研究. 物理学报, 2013, 62(9): 098201. doi: 10.7498/aps.62.098201
    [15] 黄乐旭, 陈远富, 李萍剑, 黄然, 贺加瑞, 王泽高, 郝昕, 刘竞博, 张万里, 李言荣. 氧化石墨制备温度对石墨烯结构及其锂离子电池性能的影响. 物理学报, 2012, 61(15): 156103. doi: 10.7498/aps.61.156103
    [16] 彭薇, 岳敏, 梁奇, 胡社军, 侯贤华. 锂离子电池LiMn1-xFexPO4(0x<1)正极材料的制备及性能研究. 物理学报, 2011, 60(3): 038202. doi: 10.7498/aps.60.038202
    [17] 白莹, 王蓓, 张伟风. 熔融盐法合成锂离子电池正极材料纳米LiNiO2. 物理学报, 2011, 60(6): 068202. doi: 10.7498/aps.60.068202
    [18] 侯贤华, 胡社军, 石璐. 锂离子电池Sn-Ti合金负极材料的制备及性能研究. 物理学报, 2010, 59(3): 2109-2113. doi: 10.7498/aps.59.2109
    [19] 侯贤华, 余洪文, 胡社军. 锂离子电池Sn-Al薄膜电极的制备及电化学性能研究. 物理学报, 2010, 59(11): 8226-8230. doi: 10.7498/aps.59.8226
    [20] 侯柱锋, 刘慧英, 朱梓忠, 黄美纯, 杨 勇. 锂离子电池负极材料CuSn的Li嵌入性质的研究. 物理学报, 2003, 52(4): 952-957. doi: 10.7498/aps.52.952
计量
  • 文章访问数:  1403
  • PDF下载量:  74
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-03-10
  • 修回日期:  2021-04-23
  • 上网日期:  2021-06-07
  • 刊出日期:  2021-09-05

基于电化学-应力耦合模型的锂离子电池硅/碳核壳结构的模拟与优化

  • 1. 太原理工大学机械与运载工程学院, 太原 030024
  • 2. 太原理工大学, 材料强度与结构冲击山西省重点实验室, 太原 030024
  • 3. 太原理工大学机械与运载工程学院, 应用力学研究所, 太原 030024
  • 4. 大连理工大学能源与动力学院, 大连 116024
  • 通信作者: 宋辉, songhui@tyut.edu.cn ; 郭美卿, guomeiqing@tyut.edu.cn
    基金项目: 国家自然科学基金 (批准号: 11872265)和山西省自然科学基金 (批准号: 201901D111087, 201801D121281)资助的课题

摘要: 硅基电极材料在应用中的一个主要问题是巨大的体积膨胀, 以及由此带来的电极材料破裂、粉化. 本文在有限变形假设前提下, 基于电化学-力学耦合理论, 研究球形Si/C核壳结构在嵌锂过程中的浓度、应力场的演化, 并在此基础上讨论了核壳结构的优化设计. 计算结果显示: 壳层可以很好地保护硅颗粒的膨胀; 然而核内产生的较大的径向压缩应力可能导致核壳界面的脱黏, 而核壳界面处的切向拉伸应力可能会导致壳层的断裂. 进一步为有效提高核壳结构的电化学与力学性能, 从而实现锂离子电池更长的循环寿命, 考虑了两种结构的优化: 1)单层核壳结构; 2)双层核壳结构. 结果表明对于单层核壳结构应使用更软的包覆层材料; 而双层核壳结构中优化的材料布置方案为内软外硬, 对双层核壳结构的硬度分析表明, 内层材料的杨氏模量应低于10 GPa, 而外层材料的应不高于70 GPa. 本文的结论对球形材料颗粒电极的设计及优化具有一定的指导意义.

English Abstract

    • 锂离子电池是一种可反复充放电的二次电池, 因具有工作温度范围宽、电压高、荷电保持能力强、比能量高、循环寿命长和无污染等诸多优点, 不但成为了便携式储能设备的首选, 而且在电动汽车、卫星、航空航天等众多领域发挥重要的作用 (图1)[1,2]. 科技的发展与社会的进步对锂离子电池比能量和使用寿命的需求日益提高, 对具有更高能量密度的电极材料的研究日益增加. 对于负极材料而言, 硅、锡、锗等新型材料有望代替几乎没有提升空间的石墨, 尤其是硅, 其理论比容量[3](${\rm{L}}{{\rm{i}}_{{\rm{22}}}}{\rm{S}}{{\rm{i}}_{\rm{5}}}$, 4200 mA·h/g)比传统的石墨负极材料的理论比容量(372 mA·h/g)[4]大10倍多, 因此, 硅被认为是最具发展前景的负极材料之一[5]. 然而硅作为锂离子电池的负极材料, 在完全锂化后有高达300%—400%的体积变化[6], 从而导致电极发生断裂、粉碎与电池寿命衰减等[7,8]. 这种失效情况制约着硅基锂离子电池的进一步发展. 因此, 研究硅基锂离子电池的失效机理具有重要意义.

      图  1  锂离子电池的用途

      Figure 1.  Usage of lithium-ion batteries.

      造成锂离子电池容量衰减甚至失效的原因是多方面的, 其中最主要的是电池内部产生的扩散诱导应力. 锂离子电池的工作原理是锂离子在正负极之间的移动来完成充电、放电过程. 在充电过程中, 锂离子在外加电场的驱动下从正极材料脱嵌后透过电解质嵌入到负极中, 放电时则相反. 由于锂离子电池频繁地循环充放电, 锂不断嵌入和脱出, 使得电极颗粒发生(晶格大小、晶格结构和晶相转变等)变化, 产生扩散诱导应力. 在其作用下电极颗粒可能发生破碎、粉化, 最终造成电池内阻增加, 循环性能下降, 容量衰减, 电池失效. 由此可见, 研究电极颗粒中扩散诱导应力的分布与演化, 有助于认识锂离子电池的失效机理.

      锂离子电池的充放电过程, 涉及到锂离子的扩散、电极结构的变形与内部应力场演化, 是一个物理扩散、电化学与力学的多物理场强耦合问题. 锂离子在电极中的嵌入和脱嵌会产生应力; 反之, 应力也会影响锂离子的扩散过程. Prussin[9]首次类比于温度应力推导扩散应力, 从此热应力类比法广泛应用于锂离子电池的模拟分析中 (图2). 基于Prussin[9]的热应力类比法, Zhang等[10]通过建立数学模型研究了球形颗粒和椭球颗粒的扩散诱导应力的产生与演变, 在线弹性、小变形的情况下确定了颗粒的扩散诱导应力表达式, 而且对椭球粒子的模拟结果表明, 较大的长宽比有利于减小扩散诱导应力. Song等[11]通过研究由于浓度不均匀引起的扩散诱导应力, 得到了圆柱形多层电极中的扩散诱导应力的解析解.

      图  2  电化学-力学耦合机理

      Figure 2.  Electrochemistry-mechanics coupling mechanism.

      由于耦合作用, 扩散诱导应力与扩散过程之间相互影响, 众多研究表明, 拉应力会促进锂离子的快速扩散, 而压应力则在一定程度上抑制扩散过程. 如Bagheri等[12]研究了静水应力对锂离子电池球形活性颗粒中锂离子扩散动力学的影响, 结果表明, 静水应力可以加快锂离子的扩散, 并且在恒电流状态下静水应力会导致浓度梯度减小, 从而降低应力. Lu等[13]也研究了表面张力对纳米球形颗粒中力学与电化学的作用, 数值结果表明, 在嵌入过程中电极颗粒表面的压应力抑制了锂离子的吸收.

      由于扩散诱导应力直接关系到电极颗粒的膨胀与破裂、粉化, 研究其影响因素以有效降低应力水平, 从而提高电极的机械性能就成为一个重要研究课题. Christensen[14]阐明了体积膨胀因子、非线性晶格膨胀和压力扩散对扩散诱导应力的作用, 发现压力扩散和非线性晶格膨胀在决定大体积膨胀材料的电压和应力响应方面扮演着重要的角色. Zhang等[15]建立了恒电流作用下横向各向同性圆柱状电极的扩散诱导应力分析模型, 结果表明径向和切向应力基本不依赖于杨氏模量比和扩散膨胀系数比, 而轴向应力在很大程度上取决于杨氏模量比和扩散膨胀系数比. Cheng等[16]研究了表面张力对球形纳米颗粒内扩散诱导应力的影响. 结果表明, 当颗粒直径在纳米范围内时, 应力的大小及分布对颗粒表面的断裂行为具有重要影响, 随着颗粒半径的减小, 拉伸应力的大小显著降低, 甚至可以恢复到压缩应力状态. Deluca等[17]使用电化学-力学耦合模型研究了锂离子嵌入硅颗粒过程中颗粒形态(空心颗粒与实心颗粒)和尺寸对应力产生的影响, 发现产生的应力随着颗粒尺寸减小而减小, 在纳米级时达到最小值, 并且结果表明在相同的几何尺度下, 空心球形颗粒中的扩散应力水平远小于实心球形颗粒, 表现出更好的力学性能. Li等[18]为了研究复合梯度球状电极颗粒中的力学机制, 采用了有限变形理论和可逆化学反应, 得出可逆的电化学反应导致了容量的衰退和力学性能的下降, 并提出3种非均匀梯度材料参数有利于减轻扩散诱导应力带来的负面影响. Hu等[19]建立了电极颗粒的二维电化学-力-热耦合模型, 研究了电极颗粒的形状对颗粒锂的浓度和应力的影响, 结果表明电极的几何形状在锂离子的扩散动力学中起着重要的作用. Zhao等[20]考虑到在晶界处的锂离子偏析, 建立了用于锂离子在多晶膜电极的晶界和晶界中扩散的二维解析解. 通过有限元模拟表明晶界可以帮助锂离子快速扩散并有效地降低应力.

      另一方面, 针对锂离子电池充放电过程中伴随的巨大体积变化和容量衰减问题, 研究者们提出了核壳结构的设计. 核壳结构的核心为电极材料, 壳层为起辅助改善作用的包覆层材料. 包覆层材料的作用主要有以下几个方面: 包覆层可以降低电极材料被电解液腐蚀的可能性、强化或者是带来新的物理化学性能、限制体积膨胀保持结构稳定性等. 其中碳材料具有成本较低、良好的导电性和力学性能等优点, 是目前最常用的包覆材料之一[21]. Hao等[22]建立了碳包覆正极材料锰酸锂的模型来研究核壳纳米结构的扩散应力, 结果表明, 碳层能够有效缓冲材料的体积膨胀, 降低核活性材料的拉应力. Zhao等[23]研究了包覆硅纳米管颗粒的氧化铝外壳中的扩散应力, 发现对于厚壳颗粒, 核壳界面容易发生脱黏, 而对于核较大、壳相对较薄的颗粒, 则更容易发生壳破裂. 并且借助能量释放率量化了应力和破坏模式之间的关系. 基于Zhao等[23]关于核壳结构中能量释放率的研究, Wu和Lu[24]建立了核为LiNi0.8Mn0.1Co0.1O2、包覆层为LiNi0.4Mn0.4Co0.2O2的核壳结构, 研究发现嵌锂过程中, 在壳内产生了拉伸切向应力, 可能导致壳体断裂, 而在脱锂过程中, 在核壳界面处产生了拉伸径向应力, 可能导致核心和外壳之间的脱黏, 从力学的角度解释了可能发生的两种失效模式, 并且采用能量释放率量化了核半径、壳层厚度、外加电流密度对壳层断裂和核壳界面处脱黏两种失效模式的影响, 构建了设计图用于指导核壳结构的设计. Zhang等[25]研究了恒压充放电条件下锂离子初始浓度和电极尺寸对空心碳纳米球壳电极中扩散诱导应力大小的影响, 结果表明, 初始锂离子浓度越低, 电极尺寸越小, 电极颗粒中的扩散诱导应力越小.

      核壳结构在缓释应力的同时有效地限制了电极颗粒的膨胀, 是电极材料设计的一个发展方向. 为了进一步深化对于核壳结构的认识, 本文针对电极颗粒在充放电过程中的几何大变形问题, 采用有限变形假设建立了核壳颗粒模型来描述锂离子的扩散过程和由此产生的应力, 分析研究锂离子电池的失效机理, 进而重点研究核壳结构对硅负极材料的改善作用. 对于可能出现的失效情况, 计算了目前电极材料设计领域采用的两种方案, 探索了单、双层核壳结构对应力的缓冲效果, 为锂离子电池核壳结构电极的设计提供了理论依据和优化策略.

    • 本文利用有限元软件COMSOL建立了球形颗粒模型. 如图3所示, 图3(a)是锂离子嵌入核壳结构的示意图, 整个结构分为核心1和壳层2两个域, 其中核心和壳层的半径分别为ab. 图3(b)为核壳结构沿着半径方向的二维截面, 图中rθ代表了核壳结构中的径向和切向.

      图  3  (a) 锂离子电池核壳结构示意图; (b) 径向和切向示意图

      Figure 3.  (a) Schematic diagram of core-shell structure of lithium-ion battery; (b) radial and tangential schematic diagram.

    • 核壳结构中锂离子的扩散用Fick第二定律[26]来描述. 根据动力学理论, 活性粒子在化学势μ的梯度的驱动下扩散[27], 表达式为

      $\frac{{\partial c}}{{\partial t}} = - \nabla \cdot J = - \nabla \cdot ( - Mc\nabla \mu ),$

      式中M是锂在固体颗粒中的迁移率, 根据能斯特-爱因斯坦关系, 其与扩散系数在稀溶液中存在以下关系: $D = MRT$.

      在化学势μ包含静水应力项${\sigma _{\rm{h}}}$, 体现应力对扩散过程的作用[28]:

      $\mu = {\mu _0} + RT\ln X - \varOmega {\sigma _{\rm{h}}},$

      其中, ${\mu _0}$为常数, 表示无应力状态下的化学势, R为标准气体常数, T是绝对温度, X是锂离子的摩尔分数, Ω是锂离子与电极材料发生化学反应后生成的化合物(如硅锂合金)的偏摩尔体积. ${\sigma _{\rm{h}}}$为静水应力, 定义为${\sigma _{\rm{h}}} = ({\sigma _{11}} + {\sigma _{22}} + {\sigma _{33}})/3$.

      将方程(2)代入(1)式中, 可以得到:

      $\nabla (RT\ln X) = RT\frac{1}{X}\nabla X = RT\frac{1}{c}\nabla c.$

      根据方程(2)与方程(3), 方程(1)可以化简为

      $J = - D\left( {\nabla c - \frac{{\varOmega c}}{{RT}}\nabla {\sigma _{\rm{h}}}} \right).$

      综上所述, 将简化后的扩散通量代回到方程(1), 可以得到扩散过程中浓度的控制方程为

      $\frac{{\partial c}}{{\partial t}} = D\left( {{\nabla ^2}c - \frac{\varOmega }{{RT}}\nabla c \cdot \nabla {\sigma _{\rm{h}}} - \frac{{\varOmega c}}{{RT}}{\nabla ^2}{\sigma _{\rm{h}}}} \right){{.}} $

      锂在电极颗粒中扩散时, 在电极颗粒与电解液接触表面的锂通量由颗粒表面电流密度${i_n}$决定, 根据法拉第定律二者相互转换. 因此初始条件和边界条件如下:

      $ {c}_{1}{|}_{t=0}=0, $

      ${J_2}{|_{r = b}} = \frac{{{i_n}}}{F}, $

      式中下标1, 2表示其所在的域.

    • 首先认为力学平衡过程远快于扩散过程, 因而是一个静态问题. 考虑锂化过程发生的几何大变形, 控制方程组表述如下.

      1)几何方程: 对于球对称的各向同性球形颗粒, 若球体的径向位移为u, 应变分量与径向位移u有关[29]. 在有限变形问题中, 区分参考构型与当前构型, 以下方程都在当前构型中建立. 基于连续介质力学大变形理论, 应变-位移关系可以表示为

      ${\varepsilon _{ij}} = \frac{1}{2}\left( {\frac{{\partial {u_i}}}{{\partial {x_j}}} + \frac{{\partial {u_j}}}{{\partial {x_i}}} + \frac{{\partial {u_k}}}{{\partial {x_j}}}\frac{{\partial {u_k}}}{{\partial {x_j}}}} \right),$

      式中保留了高阶项来表征真实的变形. 而小变形假设中的应变张量是线性的, 去掉了高阶项:

      ${\varepsilon _{ij}} = \frac{1}{2}\left( {\frac{{\partial {u_i}}}{{\partial {x_j}}} + \frac{{\partial {u_j}}}{{\partial {x_i}}}} \right).$

      2)平衡方程: 对于球体结构, 不考虑任何体力时, 球体的静力平衡方程为[30]

      ${\sigma _{ij,j}} = 0.$

      3)物理方程: 通过扩散诱导应力与热应力的类比, 以浓度c变量代替温度T, 偏摩尔体积Ω的作用类似于热膨胀系数α, 则本构关系可以表示为[31]

      ${\varepsilon _{ij}} = \frac{1}{E}\left[ {\left( {1 + v} \right){\sigma _{ij}} - v{\sigma _{ij}}{\delta _{ij}}} \right] + \frac{{\tilde c\varOmega }}{3}{\delta _{ij}}, $

      式中, $\tilde c = c - {c_0}$为在当前状态浓度c与初始无应力状态浓度${c_0}$之差, 表示核壳结构中锂离子浓度之间的差异. Eν分别代表杨氏模量和泊松比.

      本文考虑球形颗粒不受外在束缚的理想状态. 则有如下的边界条件:

      ${u_1}{|_{r = 0}} = 0, $

      ${\sigma _{r,2}}{|_{r = b}} = 0.$

    • 将具有高导电性、充放电过程中体积变形小等优良特点的碳作为硅负极电极颗粒的包覆层材料, 核与壳的材料参数如表1所列. 其中对于硅电极颗粒, 完全锂化状态的体积膨胀约300%, 因此其在有限变形条件下的偏摩尔体积是通过${C_{{\rm{max}}}} \cdot \varOmega = 3$来确定的[29].

      参数
      饱和浓度 Cmax/(mol·m–3)2.95 × 1052.4 × 104
      初始浓度C0/(mol·m–3)00
      扩散系数 D/(m2·s–1)1.0 × 10–161.45 × 10–13
      内径/外径 R/m4 × 10–85 × 10–8
      泊松比 ν0.230.30
      杨氏模量E/Pa8 × 10106 × 1010
      偏摩尔体积 Ω/(m3·mol–1)3/Cmax3.497 × 10–6
      室温 T/K298

      表 1  核壳结构材料参数[16,32,33]

      Table 1.  Material parameters of core-shell structure[16,32,33]

    • 本文的瞬态数值模拟经过网格与时间步长的无关性验证, 得到了如下结果.

    • 针对高能量密度所带来的体积膨胀, 计算了小变形理论的径向应力与切向应力, 与有限变形理论的数值结果进行了对比, 如图4所示. 随后模拟计算得到浓度场和应力场, 结果如图5所示.

      图  4  不同假设条件下充电过程中沿着半径方向的径向应力和切向应力 (a)小变形下的径向应力; (b)有限变形下的径向应力; (c)小变形下的切向应力; (d)有限变形下的切向应力

      Figure 4.  Radial stress and tangential stress along the radial direction during charging under different hypothesis: (a) Radial stress under small deformation; (b) radial stress under finite deformation; (c) tangential stress under small deformation; (d) tangential stress under finite deformation.

      图  5  嵌锂过程中沿着半径方向的(a)浓度和(b) von Mises 应力分布示意图

      Figure 5.  Schematic diagram of (a) concentration; (b) von Mises stress along the radius during lithium insertion.

      对比图4(a)图4(b), 可以看出, 有限变形假设下的应力在各个时刻均大于小变形假设下的应力, 而且小变形条件下的半径随着锂化的进行, 并没有表现出增大的趋势, 而在有限变形的基础下, 核壳结构的半径随着充电过程的进行越来越大, 表明有限变形假设能正确地描述材料的膨胀. 在颗粒的外表面, 径向应力为0, 这是由方程(13)所定义的. 外表面不受约束的情况下, 在锂化过程中随着锂原子的扩散, 硅颗粒开始膨胀, 内部承受拉应力. 然而, 在图4(a)图4(b)中核壳结构表现出径向压缩应力, 这是由于包覆层材料对硅体积膨胀的限制.

      图4(c)图4(d)为材料半径方向的切向应力分布图, 可以看出, 不同时刻下, 锂浓度的变化对切向应力的影响较大, 有限变形条件下的应力水平明显大于小变形. 为了限制硅核心的膨胀, 在核壳界面处产生了较大的切向应力, 这种切向拉伸应力可能会导致壳层材料的断裂.

      锂离子电池有限变形理论下的浓度场与von Mises应力场结果如图5所示.

      图5(a)给出了核壳结构在锂化过程中, 锂离子浓度沿径向的分布图. 嵌锂开始前, 整个核壳结构的浓度为0, 随着锂不断嵌入至外部的碳壳层饱和(停止条件), 硅核心材料与壳层的锂离子浓度均在增大, 而且注意到核壳结构趋向于均匀分布, 这是因为在核壳结构中的扩散应力会促进锂离子的快速扩散. 图5(b)为不同时刻碳包覆硅电极内的von Mises应力沿着径向的分布云图. 随着锂离子的不断嵌入, 颗粒会发生膨胀而产生应力, 可以看出核壳两相之间的von Mises应力差异较大, 核心材料相比于壳包覆层材料应力水平明显很低, 壳层为了保持结构稳定性以及限制硅颗粒的膨胀, 承受着巨大的应力.

      针对有限变形条件下核壳结构中可能出现的失效情形, 考虑材料和结构效应带来的影响, 模拟了以下两种结构的优化: 1)单层核壳结构, 2)双层核壳结构. 模拟所使用的参数如表2所示.

      参数杨氏模量E/Pa饱和浓度Cmax/(mol·m–3)偏摩尔体积Ω/(m3·mol–1)泊松比ν扩散系数D/(m2·s–1)
      碳层[34]1 × 1010241613.497 × 10–60.31.0 × 10–14
      外包覆层[35]3 × 10102.5 × 1043.497 × 10–60.31.0 × 10–14

      表 2  包覆层碳的材料参数

      Table 2.  Material parameters of the outer shell.

    • 为了降低核壳结构有限变形模型下的应力, 本文数值模拟了不同的包覆层碳, 以比较包覆层材料对活性颗粒结构的电化学和力学性能的影响. 图6所示为嵌锂过程中浓度、von Mises应力、径向应力和切向应力随着时间的分布示意.

      图  6  嵌锂过程中(a)浓度, (b) von Mises应力, (c)径向应力, (d)切向应力随着时间的分布示意

      Figure 6.  Schematic diagram of the distribution of (a) concentration, (b) von Mises stress, (c) radial stress, and (d) tangential stress over time during lithium insertion.

      由于锂离子浓度沿着颗粒径向基本保持均匀而且核壳界面处的应力数值最大, 因此选取了核壳界面处的一点作为随时间变化曲线图的基点. 图中黑线代表了3.1节讨论部分的核壳结构, 表示为Si/C1, 红线代表使用了不同碳包覆材料之后的结果图, 用Si/C2来表示.

      图6(a)可看出在达到相同的饱和浓度时, 使用表2参数的包覆层所用的时间更少, 表明选取合适的包覆材料可以更快地传输锂离子, 提高锂离子电池的倍率性能.

      图6(b)图6(d)可以看出: 随着时间的延长, Si/C2在核壳界面处的Von Mises应力、径向应力、切向应力都明显降低. 应力水平的整体降低体现了使用更软的包覆层材料对于减少活性颗粒的力学损伤、限制硅核心的体积变化、提高结构的稳定性有着重要的作用.

      因此, 使用更软的包覆层材料, 可以有效地提高核壳结构的电化学与力学性能, 从而实现锂离子电池更长的循环寿命.

    • 另一种抑制体积膨胀、提高电极机械性能的结构是双层壳结构, 该结构既可以缓解电极颗粒的体积膨胀, 也可以维持内部核壳结构的稳定性.

      本文在Si/C1核壳结构的外表面, 再包覆一层机械和导电性能良好的纳米级别的碳材料, 形成Si/C1/C3双层核壳结构, 在下图中用Si/C1/C3表示. 最外层碳的厚度与第一层碳的厚度一致, 均为10 nm.

      表2表1参数下构建的核壳结构和双层核壳结构进行多物理场仿真计算, 得到结果如图7所示. 图7(a)是Si/C1与Si/C1/C3双层核壳结构的浓度随着时间的分布图, 可以看出, Si/C1/C3双层核壳结构更快地达到了壳的饱和浓度, 这表明双层核壳结构有着更加优秀的扩散性能, 加快了锂离子嵌入核壳结构, 是改善负极颗粒电化学性能的有效方法.

      图  7  双层核壳结构在嵌锂过程(a)浓度, (b) von Mises应力, (c)径向应力, (d)切向应力的分布示意图

      Figure 7.  Schematic diagram of the distribution of the double-layer core-shell structure during lithium insertion process: (a) concentration; (b) von Mises stress; (c) radial stress; and (d) tangential stress.

      图7(b)为达到停止条件时单、双层核壳的von Mises应力沿着径向的分布图. 对比显示, 双层核壳的应力水平整体低于只有一层碳包覆的核壳结构, 这是内层壳与外层壳之间, 壳与核之间相互制约和相互作用所致.

      图7(c)是Si/C1核壳结构与Si/C1/C3双层核壳结构的径向应力沿着半径方向的分布图, 在有限变形的情况下, 双层核壳结构表现出了更大的径向压缩应力. 这种径向应力是核壳脱黏的原因, 这也导致了双层核壳结构相较于核壳结构有更大的脱黏风险.

      图7(d)为达到停止条件时单、双层核壳的切向应力沿着径向的分布图. 相对于单层核壳结构, 双层核壳结构的切向应力则大幅降低, 且在核壳界面处变化较为平缓, 显示了更好的结构稳定性. 根据Wu和Lu[24]的发现, 切向拉伸应力可能导致壳结构的断裂, 因此双层核壳结构有效地避免了切向应力过大而造成的包覆硅结构的碳材料的失效, 进而提高了电池的循环稳定性.

      从以上分析可以看出, 双层核壳结构有着更加优越的电化学性能, 有利于避免包覆硅材料的碳结构的断裂, 随之而来的是增加了核与壳之间的脱黏概率.

    • 双层壳为电极设计提供了一个新的自由度: 壳层材料的选取与组合. 本文考虑两层壳之间杨氏模量的影响, 将内外层壳的材料互换, 即形成Si/C3/C1双层核壳结构, 来比较对于硅电极材料的包覆效果, 结果如图8所示. 其中两种双层核壳结构分别用Si/C1/C3和Si/C3/C1来表示.

      图  8  Si/C1/C3双层核壳结构和Si/C3/C1双层核壳结构的(a)浓度, (b) von Mises应力, (c)径向应力, (d)切向应力分布示意图

      Figure 8.  (a) Concentration, (b) von Mises stress, (c) radial stress, and (d) tangential stress distribution schematic diagram of Si/C1/C3 double-layer core-shell structure and Si/C3/C1 double-layer core-shell structure.

      图8(a)以核壳界面处的一点为研究对象, 对比了Si/C1/C3双层核壳结构和Si/C3/C1双层核壳结构对于嵌锂过程中锂离子浓度扩散的影响. 由图8可以看出, Si/C1/C3双层核壳结构更快地达到了结构中锂离子的饱和状态. 与图7(a)对比, 这表明了内外层壳材料的不同, 影响锂离子扩散的反应动力学, 从而影响倍率性能.

      双层碳包覆层不仅防止了硅颗粒与电解液之间的直接接触, 也充当了硅颗粒在充电过程中体积膨胀的缓冲介质. 图8(b)图8(c)分别显示了当达到饱和浓度时, Si/C1/C3双层核壳结构和Si/C3/C1双层核壳结构的von Mises应力、径向应力沿着径向的分布. 由图可得, 双层核壳的应力在嵌锂过程中不断增加, 沿着径向Si/C3/C1双层核壳结构的von Mises应力、径向应力均低于Si/C1/C3双层核壳结构. 而且对比von Mises应力可以发现, 互换壳层材料之后的von Mises应力在核壳界面处更加平缓. 结合图7可以看出, Si/C3/C1双层核壳结构的稳定性更加优越, 因此内软外硬的双层核壳结构有利于缓解结构中的von Mises应力和径向应力. 图8(d)为当达到饱和浓度时, 内外层壳材料不同的双层核壳结构的切向应力沿着径向的分布示意图. 虽然Si/C3/C1双层核壳结构的切向应力大于Si/C1/C3双层核壳结构的切向应力, 但是对比图7(d), 仍然低于核壳结构的应力水平.

      在本文的理论框架范围内, 为了探究双层核壳结构优化, 对内外碳层材料的杨氏模量进行了进一步讨论. 图中E1E2分别表示内层壳、外层壳的杨氏模量, 且E1 = 60 GPa与E2 = 30 GPa代表了Si/C1/C3双层核壳结构. 图9(a)表示了4组不同杨氏模量下, 核与内层壳界面处一点的锂离子浓度随时间分布图.

      图  9  不同的内外层壳杨氏模量情况下 (a)浓度, (b)浓度, (c) von Mises应力, (d)径向应力分布示意图

      Figure 9.  (a) Concentration, (b) concentration, (c) von Mises stress, (d) schematic diagram of radial stress distribution under different Young's modulus of inner and outer shells.

      图9结果表明, E2 = 70 GPa时, E1 = 10 GPa与E1 = 30 GPa锂离子浓度均在3600 s内达到饱和浓度. 而当E2 = 90 GPa时, 锂离子浓度在3600 s内并没有到规定的饱和浓度(2.4 × 104 mol/m3), 说明活性颗粒的电化学性能变弱. 随着外层杨氏模量增大到300 GPa时, 双层核壳结构的嵌锂能力更差. 进一步缩小外层壳杨氏模量的范围, 以E2 = 80 GPa为研究对象, 由图9(b)可得3600 s时仍未达到壳层的饱和浓度. 因此综合图9(a)图9(b), 从锂离子的电化学性能考虑外层的硬度不能高于70 GPa, 而内层壳的杨氏模量最小可到10 GPa(文献中可查到的碳最小杨氏模量).

      本文将此优化范围内的扩散应力与Si/C1/C3双层核壳结构中的应力水平进行了对比, 得图9(c)图9(d). 图9(c)图9(d)分别显示了结构中的von Mises应力、径向应力沿着径向的变化规律, 其中红线与绿线代表了内软外硬的双层核壳结构. 数据显示在合理范围内的内软外硬双层壳结构应力水平显著低于Si/C1/C3双层核壳结构.

      综上所述, E1不低于10 GPa, E2不高于70 GPa的内软外硬的双层壳, 可对硅基双层壳结构中的应力起到很好的缓解作用.

    • 本文以Si/C核壳结构为研究对象, 建立了一个基于有限变形理论的电化学-应力耦合理论模型, 探索了锂离子浓度和应力应变场的分布. 研究表明, 随着充电过程的进行应力逐渐增大, 有限变形条件下的应力水平明显大于小变形假设下的应力水平, 同时发现在核壳界面了产生了较大的切向应力, 这种应力可能导致核壳结构的失效.

      针对这种失效情况, 研究了两种解决方案: 1)单层核壳结构, 2)双层核壳结构. 通过分析可以得出, 使用不同的包覆层材料, 能有效地改善材料内部的应力水平, 从而提高核壳结构的循环性能. 使用两种不同的材料包覆硅电极材料, 得到双层核壳结构, 可以有效降低充电过程中的扩散诱导应力. 进一步对双层核壳结构的硬度进行优化, 发现在锂化过程中内软外硬的双层壳的应力水平更低, 有助于更好地减少活性材料的失效, 并且经过计算发现优化范围为: 内层壳杨氏模量不低于10 GPa, 外层壳杨氏模量不高于70 GPa.

参考文献 (35)

目录

    /

    返回文章
    返回