网站设计原则,建e全景官网,大学校园网站模板图片,东莞市网络seo推广平台上一期讲到Eq.1中所需要求解的参数#xff0c;我们已经获取了3个#xff0c;需要注意的是#xff0c;上期所获得的电化学反应过电势已经转化成时间的函数#xff0c;接下来只需要求正负极的固相电势随时间的变化就可以了。Eq.1 固相电势与材料的soc或体相锂离子浓度… 上一期讲到Eq.1中所需要求解的参数我们已经获取了3个需要注意的是上期所获得的电化学反应过电势已经转化成时间的函数接下来只需要求正负极的固相电势随时间的变化就可以了。Eq.1 固相电势与材料的soc或体相锂离子浓度有关其电位由材料颗粒表面锂离子浓度决定。即Eq.13 我们可以通过扣电或者三电极来获得正负极平衡电势随锂离子浓度的变化关系。为了方便对比MATLAB与COMSOL里两者单粒子模型的仿真结果我们使用与COMSOL中模型相同的参数其中正、负极电位-soc曲线如下 那么对于某一时刻t正极的电位E_pos我们可以通过MATLAB里的插值函数interp1获取soc_surfcs,surf(t)/csmax_pos %%表面soc求解E_posinterp1(soc,Eeq_pos,soc_surf) %%表面电势求解这样决定全电池电位E_cell的5个参数我们就都已经求出来了。下面来看下仿真结果的对比可以看到我们搭建的模型Modeling-spb倍率放电曲线与COMSOL结果spb吻合度较高。中间部分位置电位有偏差可能来源于我们的模型中为了方便求解使用了三参数抛物线法对锂离子扩散的偏微分方程进行了降阶近似处理。COMSOL中放电末端电压直线下降可能是软件内截至条件设置的关系。 为了进一步考察我们搭建模型的适用性我们将单粒子模型中两个重要的动力学参数参考交换电流密度i0ref与锂离子扩散系数Ds调低再与COMSOL结果对比可以看到模型仍然试用。关于模型的求解时间我们的模型对1C倍率放电的求解时间是7.3秒而COMSOL是4秒。在整个模型的搭建过程中中间卡住了几次但过一段时间回头看看对参数的含义又会有新的理解建模才得以继续。整个模型还存在可以改进的地方即是对固相内锂离子扩散的偏微分方程的直接求解。MATLAB中同样有函数可以直接求解偏微分方程这个我们或许得留待以后再讨论了。至此电化学单粒子模型的搭建与校验已经完成。前面没有提到的是单粒子模型对电池内部情况做了以下几个假设1.局部电势和电解质浓度差被忽略用一个溶液内阻Rsol来代替2.固相电势梯度被忽略3.不同位置颗粒等效4.单粒子模型满足固相扩散和嵌锂反应动力学5.一般满足薄电极高导电性的电极。为什么会有过电位为什么会有过电位电池欧姆极化电化学极化扩散极化随soc的变化是什么规律它们随电流大小的变化是什么样的传荷内阻的大小和什么参数有关正负极传荷内阻的规律是一样的吗正极传荷内阻随soc的变化为什么是船形/碗形的?下一期我们提取模型中的过程参数来看一看。系列文章【原创】干货用MATLAB搭建电化学单粒子模型(上)
