首页/文章/ 详情

LS-DYNA电池行业应用介绍(二)

8小时前浏览17
   

LS-DYNA使用同一模型可以同时求解结构-热-电等多方面的多物理场问题,可以应用在电池的挤压和针 刺方面,可一次性得到结构变形信息、热信息、电流电压及SOC剩余载荷等信息。多物理场电池挤压和针 刺采用分布式等效电路模型,可以模拟电池的局部短路,模型中电池电的相关输入参数和热的输入参数可以与Fluent共用。


本系列将介绍LS-DYNA中电池结构分析相关的推荐材料模型,以及多物理场挤压、跌落和针 刺案例。


     

LS-DYNA等效电路模型


等效电路模型可用于电池的正常充放电、挤压、针 刺、跌落等滥用工况。



上图展示了常规的电化学模型,锂离子从隔膜的一极流向另一极,从而产生电流。想要模拟该微观电化学过程非常耗时,在工程上想要求解模组级别或pack级别的工况极为复杂。而我们可以将它简化为等效电路。LS-DYNA最高可支持三阶等效电路,目前使用较多的是一阶等效电路,将每一个正极和负极节点之间简化成等效电路。用户只需输入参数(或曲线) R0,R10,C10等,LS-DYNA可自动建立等效电路进行求解。



随后利用电路的发热r0乘以i2,考虑可逆放热以及不可逆的焦耳热,同时考虑将热量反馈给电路,以及结构进行求解。



如何获得电路的参数?用户需进行不同SOC、温度与电流方向下的HPPC实验,拟合成相应的参数进行求解。



上图展示了某个短路工况,在结构变形最剧烈的区域,自动调用短路判距*EM_RANDLES_SHORT,发生局部短路时,使用阻值比较小的r0替换原先的等效电路,并开始发生急剧的放电。



用户可以自行在K文件中编写C语言函数定义内短路判据。



LS-DYNA提供四种电路搭建模型,如实体单元模型(包含电芯所有层数的模型),可以使用*EM_RANDLES_SOLID关键字定义,进行正常的充放电。



同时,如果事先已定义了局部短路判定准则,还可以进行挤压仿真,挤压到一定程度,电流从另一侧流出,后处理软件中可以查看电流矢量图。



Ford进行的电池实验和仿真对比,实验包括了两种类型的电池A和B,材料与尺寸略有差别。分别进行HPPC放电实验后与LS-DYNA仿真结果对比。



上图可以看到实验与仿真的电压数据基本吻合。



后续可进行正常的充放电实验。上图展示了在不同倍率下进行放电实验(黑色曲线)与仿真(其他颜色曲线)对比,温度变化、电压变化、整体趋势,以及峰值等关键数据与实验数据吻合较好。



外短路模型设置。



假设该电池模型发生外短路现象,使用*EM_ISOPOTENTIAL_CONNECT关键字短路连接bus bar,发生短路和放电,热量从温度较高的bus bar逐渐传向内部。



在电池包里设置多个监测点,不同位置的电芯的电流、变化、SOC变化均可以在后处理软件中画出。右侧图表为实验(实线)中不同位置的温度与仿真(虚线)结果对比。



上图案例展示了单电芯的圆柱挤压工况下造成短路放电的实验及数据分析。



建立LS-DYNA实体单元模型,为每一层的电芯进行建模。



上图展示了仿真与实验的电压及温度对比。



若要使用*EM_RANDLE_SOLID进行Pack模组模型(包含上百层电芯)的挤压仿真,巨大的网格数需要消耗大量的计算资源。此时模型中可使用厚壳单元加速计算,每个厚壳单元仅包含1~2层电芯即可满足,且可考虑厚度方向的积分点上电的信息。



通过*PART_COMPOSITE_TSHELL在积分点上定义电材料,以及结构材料,包括比热容和电导率。



小球撞击厚壳单元建立的电芯模型,由于结构单元数量较少,计算速度也相对较快。



圆柱撞击厚壳单元建立的电芯模型仿真案例,动画模型中可以观测到电压变化,电流密度变化,SOC变化,以及相关云图。



除了方型和软包电池外,LS-DYNA电池模拟同样适用于圆柱型电池。上图为厚壳单元建立的圆柱电池模型挤压模拟。

- END -

       


来源:摩尔芯创
LS-DYNAFluent化学电路半导体光学电力电子材料
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-12-05
最近编辑:8小时前
摩尔芯创
光学仿真、光学培训、硅基光电子
获赞 6粉丝 5文章 42课程 0
点赞
收藏
作者推荐

LS-DYNA | 如何将隐式计算转为显式计算?

01.摘要本文介绍了4种隐式计算转化显式计算的方法,分别为重启动方式、*INTERFACE_SPRINGBACK_LSDYNA关键字、载荷曲线设置和隐式动态松弛方法。 02.正文隐式分析后再进行显式计算(隐式-显式顺序求解),有以下4种方法: 1、重启动 隐式分析结束后,以完全重启动方式提交显式计算。该方法应用广泛,可以在转换计算时调整模型,如添加或删除零件、增加或修改边界条件、载荷等。 2、*INTERFACE_SPRINGBACK_LSDYNA关键字 在隐式分析k文件中加入 *INTERFACE_SPRINGBACK_LSDYNA关键字。在计算完成后,将创建一个名为“dynain”的输出文件,文件中含有变形几何信息与初始应力信息。Dynain文件与keyword格式一致,将此文件与初始的输入文件(替换原始*NODE和*ELEMENT_SHELL命令),去掉隐式选项,重新进行应力初始化显式计算。无需重启动;显式计算是一个全新的分析,该分析是根据之前隐式分析的结果初始化的。请注意,dynain文件不能初始化梁单元或离散元,也不能初始化接触力。 3、载荷曲线设置 可以设定一条载荷曲线,该曲线设定哪些时间段进行隐式计算,哪些进行显式计算。曲线的横坐标是时间,隐式计算时纵坐标是1.0,而在显式计算时(曲线是阶跃函数)纵坐标是0.0。将*CONTROL_IMPLICIT_GENERAL中的IMFLAG设置为-|curve ID|。无缝从一种分析切换到另一种分析,并且没有CPU或I/O开销。 4、隐式动态松弛方法 进行隐式动态松弛计算以实现应力初始化。来源:摩尔芯创

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