HTML
-
研究对象为锂离子电池负极极片材料,该材料为“三明治”结构,由铜箔及其上下表面涂覆的等厚度石墨层组成,材料总厚度为100μm,其中铜箔厚度为10μm,材料结构如图 1所示。切割过程中,激光束辐照至负极材料上表面,上层石墨材料吸收能量后温度迅速上升至气化点,同时热量传导至中间金属层与下层石墨,产生材料的去除[14]。通过扫描振镜控制激光光斑位置,调整其沿切割路径连续移动,由此完成切割过程[15]。
激光切割过程中,光斑所在区域的温度随时间急剧变化,激光切割温度场求解属于非线性的瞬态导热问题。温度场分析基于能量控制方程,利用傅里叶定律和能量守恒定律建立导热微分方程求解,求得内部各点的热流密度矢量,获得材料体内的温度场。为简化计算条件,提出3点假设[16]:(1)材料是各向同性的连续介质;(2)考虑材料物性参量(导热率、比热以及密度大小)的温度依赖性;(3)忽略石墨层与铜箔间的热损耗以及热阻的影响。
-
在激光与材料的作用过程中,温度t随时间τ变化,属于非稳态导热,直角坐标系中3维非稳态导热微分方程的一般形式为[17]:
式中,ρ(kg/m3)为材料的密度,c为材料的比热容,λ[W/(m·K)]为材料的热导率,q′为内热源。根据本节中对模型材料常物性的假设,对应该模型的导热过程,可将导热微分方程简化为物性参量为常数、无内热源的傅里叶方程,即:
式中,a=λ/(ρc), x, y, z为空间坐标。
-
初始条件下,设置初始温度t0=283.97K。由于材料在激光烧蚀过程中温度急剧升高,其温度变化还需考虑光斑辐照区域与周围介质间的对流换热与辐射换热。因此,在激光热源辐照区域内,激光热源以热流密度形式加载在材料上表面,除此之外,所有表面均设置为对流换热边界与辐射换热边界,热传递量q的对流边界条件[18]可表示为:
式中,Δt指极片材料表面与空气流体间的温度差,h表示表面换热系数。
材料的辐射边界条件,即极片与周围空气间的热传递量q可以表示为:
式中,ε为辐射系数,取ε=0.7,S表示极片材料的面积,tm为材料温度,ta为周围空气的温度,初始条件下,tm=ta=283.97K, 斯忒藩-玻尔兹曼常数σ=5.67×10-8W(m2·K4)。
激光切割过程中,温度的高低、温度梯度的大小取决于板材的受热面积、板材吸收的热流密度和加热时间[19-20],除(3)式、(4)式所表达的对流边界条件和辐射边界条件外,将激光热源视为施加于材料表面的一类边界条件。本文中的研究对象为厚度约100μm的负极极片材料,高斯热源模型可描述其切割温度场的分布。高斯热源模型是一种面热源,切割过程中通过激光热源的加载为材料去除提供热能,其中表面有效加载区域为光斑。加热区域内热流密度近似为高斯函数分布,热源S0的数学表达式[20]如下:
式中,γs为材料表面的反射系数,δ为材料沿厚度方向的吸收率,(x0, y0)为光斑中心坐标,r0为光斑半径,I0为光斑中心强度, P为激光功率。根据高斯分布理论,光斑中心强度可表示为:
入射激光波长为1065nm,该条件下石墨和铜的材料特性参量如表 1所示。为准确地模拟复合材料的特性,表面吸收性质与水平面传播方向上的热传导表现为石墨材料的性质。
properties copper graphite density/(kg·m-3) 8.96×103 1.73×103 thermal conductivity/(W·m-1·K-1) 317 18.1 specific heat/(J·kg-1·K-1) 385 2092.48 surface reflectance 0.97267 0.16815 melting temperature/K 1358 3652 boiling temperature/K 2835 4800 Table 1. Material properties of copper and graphite
-
仿真过程为3-D瞬态热分析,采用ANSYS solid 90单元,适用于3-D的稳态或瞬态热分析问题。如图 2所示,该单元为3-D 20个节点的热实体,每节点具有一个温度自由度。采用扫掠方式划分六面体网格,按照载荷特征,根据网格与切缝的距离决定划分精度,由外向内逐渐提高,对切缝部分的网格进一步细化,该模型一共具有5252个节点,3600个单元[21]。考虑到切缝两侧的对称性,故只取切缝一侧进行求解。由此建立材料3-D模型,其尺寸为1mm×0.5mm×0.1mm。