【学术论文】白车身电泳工艺的仿真分析与实验研究廖毅1，刘强强1，许冰1，李奥2，杨磊3,4，杨锐4，肖刚4 (1.上汽通用五菱技术中心,广西柳州 545007; 2.湖南大学汽车工程学院,湖南长沙410082; 3.上海交通大学金属基复合材料国家重点实验室，上海 200240; 4.上海格麟倍信息科技有限公司, 上海 200120) 摘要：针对汽车电泳过程中，现场影响因素复杂，过程控制要求高，质量问题分析不足，根据流体N-S方程、连续方程、动量守恒方程、能量守恒方程及麦克斯韦方程建立全耦合数学计算模型，应用ECoatMaster有限元分析软件，对三维状态下的汽车电泳工艺进行全耦合场计算机仿真分析，得到整车电流密度分布及漆膜厚度分布。结果表明在全耦合计算机仿真分析过程中应用ECoatMaster软件能有效计算电泳液属性、阳极管排布、电压特性、行车轨迹特性对汽车电泳性能的影响，所得整车外板及B柱内腔膜厚的不同开孔位置膜厚差异，并与已有实验结果及研究数值吻合良好。 关键词：电泳工艺；全耦合；数值分析 中图分类号：TG401 文献标识码：A Numerical analysis of E-Coating process under the fully coupled conditions LIAO Yi1, LIUQiang-qiang1, XU Bing1, LI Ao2, Yang Lei3,4, YangRui4, Xiao Gang 4 (1. SGMW Automotive Company Technology Center, Liuzhou 545007; 2. Automotive College, Hunan University, Changsha 410082; 3. State Key LAB of Metal Matrix Composites, Shanghai Jiaotong University, Shanghai 200240; 4.Shanghai GLB Infor&Tech.Ltd, Shanghai 200120） Abatract: Aiming at the defects predict of Automotive E-coating process in difficult conditions and highquality control. According to the N-S equation, continuality equation, momentum conservation equation, energy conservation equation and Maxwell equations. Applied in ECoatMaster finite element software to analyze the E-Coating process 3D field under fully coupled conditions to predict electric current distribution and painting film distribution. The results show that the paint characters, anodes distribution, voltagesetting and vehicle trajectory will effect film distribution，The film distribution in in outer vehicle and B-pillar is different based on the sharp and product design. It is good agreement with experimental and research result. Key Words: E-Coating process; Fully Coupled; Numerical Analysis 0 前言 电泳过程是一个非常复杂的物理－化学过程，影响因素众多，而且当前由于对汽车安全性、NVH特性等要求越来越高，车身设计也越来越复杂，白车身内腔结构如A柱、B柱、门槛、横梁、前纵梁等部位往往存在涂层厚度较小或者根本电泳不到的问题。汽车电泳过程伴随流体力学、能量转换、能量传输、电场作用等一系列动态变化的过程，很难通过直接观测的方法得到电泳工艺的准确行为特征，而数值分析则可在一定的数学模型下提供完整的流场、电场及化学特性参数信息。由于汽车电泳过程中涂料特性及工艺的控制参数稳定对整体车辆防腐质量评价起到关键性作用因此国内外众多研究者对汽车的电泳工艺开展了大量的研究，并建立了相对成熟稳定的数学模型。XiaoHSong为计算电泳工艺建立了直角坐标系下的电场模型，并研究了电泳力大小和方向对电场的影响行为，使其更接近实际的结果。 Elsyca Paulo等根据流体守恒控制方程组和麦克斯韦方程组对整个电泳过程电场及流场区域建立了三维的数学模型，考虑了电流分布、流动速度、行车轨迹、涂料属性等之间的相互作用，对白车身的电泳工艺做了认真的分析和研究，同时对比二维的计算结果发现此种模型更加符合实际，同时也研究了不同阳极条件和外力作用下涂装性能的改变。 实际电泳过程中存在电场作用、流场作用的同时，温度场的变化也对涂料性能及电场强度及方向有影响，本文建立了不同方程作用下的全耦合数学模型，利用ECoaECoatMaster软件对上汽通用五菱汽车有限公司某款车型进行了全耦合场分析研究，是研究更符合实际，并对白车身电泳工艺进行指导和改进。 1 模型的数学描述及分析 1.1 控制方程 描述电泳工艺的控制方程组如下式（1～11）所示： 质量连续方程为
式中：ρ为流体的密度，u为流体在x方向的速度，v为流体在y方向的速度，w为流体在z方向上的速度。 动量守恒方程为
式（2～4）分别为流体在3个方向上的动量守恒方程。其中，μ为流体的粘度，fx为x方向流体所受的体积力，fy为y方向流体所受的体积力，fz为z方向流体所受的体积力，其中fx、fy和fz主要为洛伦兹力，而在z方向时还应考虑重力的作用。 能量守恒方程为
式中：k为电弧介质的热传导系数，Cp为电弧介质的比热容，St为热源项。
式中：j为电流密度，o为电导率，V为电势，B为磁感应强度，A为磁矢量式，u0为真空磁导率。 1.2 几何模型及网格划分 根据现场工况条件分别建立整车数模、电泳槽数模、阳极排布数模，使三维模型符合实际工况，如图1所示：
图1. 现场工况的有限元网格划分 在求解过程中，为了得到更加精确的结果，需在变量梯度大的区域划分细密的网格，故该模型采用的是自由化网格。由于电极附近温度、电流密度、速度等变量的梯度都很大，故该电极区域选用细密网格。模型内部由于是电场较难到达部位及电泳工艺关心位置，该故此区域也选用细密网格。在求解域的内部，中通道附近的电流密度梯度也很大，故该区域也需选用细密网格。除了上述需要特别注意的区域外，求解域的其他部分可以采用相对稀疏的网格，以达到在求解精度满足要求的前提下减小计算量，如图2：
图2. 白车身数模网格划分 1.3 边界条件 在对电泳所涉及的温度场、流场、电场进行求解时，需要对其边界条件进行设定。在温度场中，考虑电泳漆特性变化属性。在流场中，电极及工件表面设为无滑动的壁面，而电场求解域的四周则设为开边界。在电磁场中，电极设为向内的电流密度，而求解域的白车身设为接地，四周电泳槽及其他边界则设为电绝缘。白车身电流随电压分布，如图3所示：
图3. 电泳过程中的电流趋势 1.4 模拟参赛的选择 模拟过程中介质为涂料及气体，在计算过程中需要对求解域中的介质进行属性设定，分段工作电压选择、行车轨迹运行，其参数的选择如表一所示：
2 计算结果与分析 2.1 控制方程 漆膜厚度分布 外板膜厚结果如图4
(a)
(b) 图4. 整车外板膜厚模拟结果 内板膜厚分布如图5：
图5. 内板膜厚分布 2.2 仿真值与真实值对比
图6. 仿真值与实际测量值对比 3 结论 利用ECoatMaster有限元电泳分析软件对三维状态下,电泳过程电场、流动场及涂料电泳能行进行全耦合计算机仿真模拟。 1、在能量守恒、动量守恒质量守恒三大控制方程状态下利用MaxWell方程求解计算电泳模型是可行、可靠的。 2、在全耦合条件下求解得出的电流密度分布、电泳膜厚厚度分布，符合在大电流状态下有利于漆膜厚度的生长，随漆膜厚度增加，电泳阻抗增加，电流密度降低，向内板转移，使内板电泳膜厚增加。与实际测量值相吻合。 参考文献 [1] 李建 ,张新成 .电泳工艺电场数值模型研究 [J] 南京工业大学学报, 2006, 27（11）；25-28. [2] 李婷婷 ,金超 .汽车电泳工艺计算机仿真分析[M] 吉林大学, 2012,2； [3] ECoatMaster Manual [J] Elsyca ,2010,15-23； |
