汽车文摘
李鼎1,2
于保君1,2
刘创举1,2
朱学武1,2
曹正林1,2
（1.中国第一汽车股份有限公司研发总院，长春130013；2.汽车振动噪声与安全控制综合技术国家重点
实验室，长春130013）
Method ）与有限体积法（Finite
Volume
Method ）的理论基础、多相
流仿真分析方法及重叠网格技术理论及应用，分析了流固耦合仿真应用的国内外现状及电泳过程车身结构变形预期仿真
方法可能存在的挑战。
主题词：涂装电泳有限元法有限体积法流固耦合
中图分类号：U466
文献标识码：A
DOI:10.19822/jki.1671-6329.20200205
Overview on Simulation Method of Body Structure Deformation During
Electrophoresis
Li Ding 1,2,Yu Baojun 1,2,Liu Chuangju 1,2,Zhu Xuewu 1,2,Cao Zhenglin 1,2
（1.General Research and Development Institute,China FAW Corporation Limited,Changchun 130013;2.State Key Laboratory of Comprehensive Technology on Automobile Vibration and Noise &Safety Control,Changchun 130013）【Abstract 】Coating is an important process in the manufacture of car body products,and electrophoresis is an
important link in the coating process.In the process of electrophoresis,due to the impact of liquid and the movement of the body in white,there will be sheet metal deformation,structural degumming and insufficient strength of fixed points,which will lead to uneven size matching and large gap,which will seriously affect the user ’s visual experience.In this paper,the theoretical basis of finite element method and finite volume method,the theory and application of multiphase flow simulation analysis method and overlapping grid technology,the current situation of fluid structure coupling simulation analysis at home and abroad,and the possible challenges in the simulation method of body structure deformation in
electrophoresis process are summarized.
Key words:Coating electrophoresis,Finite element method,Finite volume method,Fluid
structure coupling
【欢迎引用】李鼎，于保君，刘创举，等.电泳过程车身结构变形仿真方法概述[J].汽车文摘,2021(4):39-45.
【Cite this paper 】Li D,Yu B,Liu C,et al.Overview on Simulation Method of Body Structure Deformation During Electrophoresis[J].
Automotive Digest (Chinese),2021(4):39-45.
电泳过程车身结构变形仿真方法概述
缩略语
FVM Finite Volume Method FSI Fluid-Structure Interaction VOF Volume Of Fluid
FEM Finite Element Method
FSI
Fluid-Structure Interaction
CFD
Computational Fluid Dynamics
CSD Computational Structural Dynamics
1前言
电泳涂装工艺一般由涂装前预处理、电泳涂装、电泳后清洗、电泳涂膜的烘干等多道工艺组成，如图1所示。为提高电泳质量和生产节拍，新型翻转线、穿梭机等电泳线逐步被采用，而与之带来的是白车身闭合件的变形、生产过程中车身倾斜和运动装置的脱钩
39
等现象在电泳过程中时有发生。
图1电泳涂装工艺示意
传统结构强度设计主要考虑整车的碰撞、强度、模态和刚度参数，并未充分考虑白车身在各大工艺过程的应力应变累积情况。而白车身在冲压、焊接工艺设计时，充分考虑了材料、模具、夹具因素的影
minicooper价格响，白车身焊接总成整体变形量不大。而当焊接总成进过涂装生产线，由冲压、焊接引入的残余应力，电泳泳池中的流体压力积累的塑性变形，烘烤引入的材料应力释放、热应力变形，最终造成总装时部分车身结构面差较大，无法顺利完成装配。
由于白车身经过涂装后，变形影响因素复杂，因此，在研究此类系统问题时，首先需要将问题进行分解，然后分析不同工序对结构变形的影响，出设计的薄弱环节，加以优化，再通过多学科仿真工具将全生命周期下的变形情况进行分析和验证。
白车身在电泳槽中的受力主要来自其在泳池中改变姿态而引起的阻力、浸入液面以下受到的浮力、自身重力、工装的约束力，因此可以将模型抽象为求解流体力学及结构力学问题。流体力学部分需要考虑气-液-固3相的交互作用，目前考虑多相流的方法主要是有限体积法（Volume Of Fluid,VOF），由于白车身在泳池中是运动的，因此若采用网格算法进行仿真，则需要采用重叠网格技术。考虑到流体-固体间的交互作用并不显著，因此计算时可以采用单向耦合法进行计算，即先计算不同位置、时间下白车身各零部件表面受到的流体压力情况，再将其作为压力边界条件导入到结构力学计算中。由于白车身钣金为标准固体模型，目前金属弹塑性变形数值模拟中的一种有效数值计算方法为有限元法，将求解对象离散为由各种单元组成的计算模型，计算流体载荷激励下各离散单元非线性位移响应，从而获取钣金残余变形。
当前对电泳过程车身结构变形仿真方法研究尚处于起步阶段。ESS斯太尔工程软件公司提出了一种基于有限差分法对工业EPD涂层过程进行准确有效的模拟方法[1]。埃尔朗根-纽伦堡大学的Dominik Bar⁃tuschat介绍了一种用于微流控电泳大规模并行直接数值模拟的耦合算法，这种多物理算法采用了流体和离子的欧拉描述，结合了移动带电粒子的拉格朗日表示[2]。卡尔斯鲁厄理工学院用数值方法证明了分段式对电极在不同电位差下可以有效地减少涂层材料的用量，将现有的模型扩展到三维空间中的任意非平面几何[3]。上汽通用五菱汽车股份有限公司应用电泳仿真CAE技术对车身电泳膜厚进行了仿真分析，提升了产品的强度和防腐性能[4]。Institut für Computerphysik 利用分子动力学模拟方法研究了聚电解质在平行非带电平板间的电泳拉伸，模拟结果显示，流体动力相互作用的影响，减少了限制板之间的距离[5]。
2流体侧仿真
电泳过程车身结构流固耦合仿真主要采用计算流体力学计算流体部分。白车身在电泳槽中运动，会产生复杂多变的流体流动特性，电泳液的流动符合湍流流动特性，而车身运动中包含进出槽过程，当中会卷入空气进入车身腔体内部，电泳液在车身表面存在浸润、物质交换的物理化学过程，携带的电泳液在出槽过程会有流淌现象。因此白车身电泳槽的运动过程实质上是一个多相流问题。根据具体工程应用，可以将电泳液定义为流体的液相，空气定义为气相，整个流动过程符合气液2相流计算模型。
2.1有限体积法网格化理论
有限体积法（Finite Volume Method，FVM）是近年发展非常迅速的一种离散化方法，其特点是计算效率高，目前在流体领域得到了广泛的应用。其基本思路是：将计算区域划分为网格，并使每个网格点周围有一个互不重复的控制体积；将待解的微分方程（控制方程）对每一个控制体积分，从而得到一组离散方程。其中的未知数是网格点上的因变量，为了求出控制体的积分，必须假定因变量值在网格点之间的变化规律。从积分区域的选取方法看来，有限体积法属于加权余量法中的子域法，从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。北京航空航天大学流体力学研究所的朱自强给出了一种生成分区对接网格和Euler方程分区解，以二维多段翼型为对象分别讨论了分区重叠结构网格方法、非结构网格方法和自适应的笛卡尔网格法[6]。杨建宏通过有限元空间和有限体积元空间的一种双射投影，得到了不可压缩流问题低次等阶稳定有限体积元方法[7]。D.
G.Roychowdhury讨论了计算流体力学中的各种控制方程及其推导，以及偏微分方程描述边界条件的物理和数学意义[8]
。40
2021年第4期
对于电泳过程车身高度复杂的几何，捕捉此几何可能比较困难。在这类情况下，包面可用于简化和捕捉复杂几何的水密表示。无论是否使用包面，通常都使用表面重构来重构几何的初始表面。几何的初始表面通常由许多三角表面组成。此三角形化也称为网格化，通常包含高度偏斜的三角形，不适用于生成高质量体网格。重构通过生成尺寸更均匀的三角形提高这些表面的整体质量，此方法最适合网格模型。然后，使用重构表面作为体网格的基础。
2.2多相流仿真方法
目前在多相流计算中，如果采用网格算法，通常采用欧拉-欧拉法和欧拉-拉格朗日法2种方法，其中欧拉-欧拉法常用流体体积模型（VOF）、混合物模型及欧拉模型，欧拉-拉格朗日法常用离散相模型。
VOF模型是一种在固定欧拉网格下表面追踪的方法，它通过研究网格单元中流体和网格体积分数来确定自由面，追踪流体的变化，而非追踪自由液面上质点的运动。在计算过程中，利用各流体的体积分数为计算参数，在一个动量方程中进行求解，把每个时刻各流体所占的体积比都记录下来。
体积追踪/交界面捕捉法非常适用于模拟具有大尺度交界面的分层多相流体。必须注意的是，计算网格
单元的尺寸在将交界面分类为大尺度还是小尺度时至关重要。由于体积追踪/交界面捕捉法依赖于完全求解交界面，因此在离散流态建模精度同样重要的多尺度流体模拟中使用这些方法会导致成本过高。流体体积(Volume Of Fluid,VOF)多相模型实施属于交界面捕捉方法系列，可预测不混溶相交界面的分布和移动。此建模方法假设网格分辨率足以求解相之间的交界面位置和形状。
交界面的相分布和位置由相体积分数∝i的场来描述。相i的体积分数定义如式（1）：
∝i=V i V（1）式中，V i为为网格单元中相i的体积，V为网格单元的体积。
网格单元中所有相的体积分数总和必须是式（2）：
∑i=1N∝i=1（2）式中，N为总相数；∝i=0表示网格单元完全没有相i；∝i=1表示网格单元完全由相i填充；0<∝i<1表示2个极限之间的值存在相间交界面i。
包含交界面的网格单元中计算的材料属性依赖于组成流体的材料属性。同一包含交界面的网格单元中的流体被视为混合物，其中混合物的密度、粘度及比热如式（3）（4）（5）：
ρ=∑iρi∝i（3）
μ=∑iμi∝i（4）
C
p
=∑i(C p)iρiρ∝i（5）式中，ρi为相i密度，μi为相i动力粘度，(C p)i相i 比热。
体积分数传输方程相i的分布由相质量守恒方程驱动，如式（6）：
∂∂
t∫V
0a
昌河铃木浪迪t d V+
∮A a i v d a=∫V0(S a i-a iρi Dρi D t)d V-∫V01ρi∇(a iρi v d,t)dV
（6）式中，a为表面积矢量，v为混合（质量平均）速
度，v d,t为扩散速度，S a
i
为相i的用户自定义源项，
D
ρ
i
D
t 为相密度ρi的材料或拉格朗日导数。
计算采用如下所示计算相体积分数：
1）存在2个VOF相时，仅对第1个相求解体积分数传输。在每个网格单元中，调整第2个相的体积分数，使这2个相的体积分数总和等于1。
2）存在3个或更多VOF相时，对所有相求解体积分数传输。然后，根据每个网格单元中的所有相的体积分数总和标准化每个相的体积分数。
连续性方程所有相的总质量守恒方程由式（7）给出：
∂∂
t
()
∫V0ρd V+∮Aρv d a=∫V0S d V（7）式中，S为与相源项相关的质量源项，如式（8）所示：
S=∑i Sαiρi（8）通过由式（3）给出的密度来解释对流体混合物组成相的体积分数的依赖性。
动量方程及能量方程如式（9）式（10）所示：
∂∂
t∫V
0ρv d V+∮0Aρv×v d a=-∮0A pI d a+∮0A T d a+∫V0ρg d V+∫V0f b d V-∑i∫A0αiρi v d,i×v d,t da
（9）式中，p为压力，I为单位张量，T为应力张量，
f
b为体积力矢量。
∂∂
t∫V
0ρE d V+∮0A[ρHv+p+∑iαi p i H i v d,i]d a=
-∮0A q″d a+∮0A Tv d a+∫V0f b v d V+∫V0S E d V（10）
41
汽车文摘
2021年第4期
工程中关注的大多数流体流都具有不规则的波动流量。此类波动通常尺度小但频率高，要在时间和空间上对其求解，需要高昂的计算成本。可以不求解湍流的精确控制方程（直接数值模拟），改为求解平均量或滤波量来近似小波动结构的影响，这样成本较低。各种湍流模型为此类结构的建模提供了不
同的方法。2.3
重叠网格技术
由于白车身在电泳槽中是运动的，如果采用网格重构技术，网格计算量很大而且会造成后期结构有限元计算网格节点不统一。因此需要使用重叠网格技术。重叠网格也称作“嵌套”网格，用于离散具有以任意方式相互重叠的多个不同网格的计算域。这些网格在处理涉及多个体或移动体的问题时以及优化研究中最为有用。Tapan K.Sengupta 通过计算二维Navier-Stokes 方程，使用单块结构和重叠网格来求解圆柱绕流，利用相对较少的网格点捕获了该布置的已
知流动特性[9]。
中国船舶科学研究中心的赵发明采用RANS 方法和重叠网格计算了带自由液面的船舶绕流问题，描述了重叠网格和单相Level Set 自由液面模拟方法的数学模型及求解过程[10]。上海交通大学船舶海洋与建筑工程学院的王建华采用重叠网格技术开发的多功能水动力学求解器naoe-FOAM-SJTU ，对标准船模DTMB5512裸船体在平面运动机构(PMM)控制作用下对纯摇首运动进行了数值模拟[11]。海军工程大学的周广礼以Suboff 标模为研究对象，基于RANS 方程、VOF 模型及重叠网格技术，建立了潜艇应急上浮运动的数值预报方法[12]。
在大多数情况下，使用重叠网格不需要在生成初始网格之后进行任何网格修改，比标准网格化方法更灵活。任何涉及重叠网格的研究都具有封闭整个求解域的背景区域，以及包含域中的体的一个或更多较小区域（图2）。
图2重叠网格示意
在重叠网格中，网格单元分组为活动、不活动或
郭美美 玛莎拉蒂受体网格单元。在活动网格单元中，将对离散控制方程进行求解。在不活动网格单元中，不会对任何方程进行求解，但是如果重叠区域在移动，这些网格单元可以变为活动状态。受体网格单元可分离背景区域中的活动和不活动网格单元，并与重叠区域中的重叠边界相连接。受体网格单元用于耦合2个重叠网格上的求解。一个网格的供体网格单元处的变量值，通过插值表示其他网格中受体网格单元处的变量值。供体网格单元是其他网格中最靠近受体网格单元的活动网格单元。
将同时为所有区域中的所有活动网格单元计算求解，即网格为隐式耦合。当引用离散方程中一个网格的受体网格单元中的变量值时，将使用其他网格中供体网格。
3固体侧仿真
电泳过程中，车身变形是由液体冲击及车身自身运动导致，因此需将流体侧计算获得的不同时刻下的
流体压力值，映射到结构计算的有限元网格上，作为压力边界条件计算整个过程中白车身的弹塑性变形情况。3.1
有限元法网格化理论雪铁龙价格
有限元法（Finite Element Method,FEM ）是广泛应用的固体力学数值计算方法。有限元法是将一个连续的求解域任意分成适当形状的许多微小单元，并于各小单元分片构造插值函数，然后根据极值原理（变分或加权余量法），将问题的控制方程转化为所有单元上的有限元方程，把总体的极值作为个单元极值之和，即将局部单元总体合成，形成嵌入了指定边界条件的代数方程组，求解该方程组就得到各节点上待求的函数值。有限元法的基础是极值原理和划分插值，它吸收了有限差分法中离散处理的内核，又采用了变分计算中选择逼近函数并对区域积分的合理方法，是这2类方法相互结合，取长补短发展的结果。它具有广泛的适应性，特别适用于几何及物理条件比较复杂的问题，而且便于程序的标准化。Nathan Mendes 在Basics in Practical Finite-Element Method 一书中详细介绍了有限元法及其应用[13]。庄
茁在基于ABAQUS 的有限元分析和应用中详细介绍了有限元的基本理论和数值计算方法，并对相关例题进行了讨论研究[14]。杜平安详细介绍了划分有限元网格时应考虑的一些基本原则，包括网格的数量、
疏密、阶次、质量、分界点、布局、编号等[15]。3.2
多体动力学仿真
42
汽车文摘
白车身在电泳槽中的运动模式与涂装工艺相关，如图3所示。在自身重力的作用下，不同的运动模式对车身变形有相应的影响，因此在模拟电泳过程固体侧变形时，需考虑动力学影响。
图3车身电泳过程示意
应用拉格朗日待定乘子法，动力学方程如式（11）：
新帕萨特2013款éëêùûúM ΦT q Φq 0æèçöø÷q λ=æè
çöø÷F （q,q ,t ）γ（11）式中，Φq 为约束方程雅克比矩阵；M 为质量矩阵；λ为拉格朗日乘子；q 为加速度；
F(q,q ,t)为作用力；γ为加速度公式中二阶导数项。
其中，系统约束方程为式（12）：
Φ()q,t =0
（12）
求约束方程广义坐标的偏微分，可得到约束方程
雅克比矩阵Φq ：
Φq =∂Φ
∂q
（13）
速度方程为式（14）
Φq q +Φt =0（14）加速度方程为式（15）
Φq q +(Φq q )q q +2Φqt q +Φtt =0
（15）
由加速度方程得到的拉格朗日乘子γ为式（16）
γ=-(Φq q )q q -2Φqt q -Φtt （16）3.3
结构力学非线性仿真
与流体方程不同，固体应力在拉格朗日架构中求解，网格依据材料中的材料点。这意味着特定网格单元中的质量将在变形时保持不变，因为它在物理上代表相同的材料。此网格单元质量由材料原始密度与原始体积的乘积来确定。因此，变形状态下的密度是原始质量除以网格单元的当前体积。
固体力学用于描述固体连续体对应用负载的响应行为。应用负载包括体积力、表面负载、点力或固体温度变化导致的热负载。应用负载会在结构中产
生应力场，并且可能导致结构位移从初始未变形配置到变形配置。
应用的负载可导致固体结构产生从初始配置到变形配置的位移。总位移由刚体运动和体中各点的相对位移（决定固体结构的变形）共同组成，如图4所示。
图4
各节点总位移示意[16]
如果未变形配置中材料点的位置为X ，且此点相对于变形配置的位移为u (X ,t )，则该材料点在变形配置中的位置为式（17）：
x ()X,t =X +u (X,t )
（17）采用分量形式，位移可以表示为式（18）：
u ={}
u x u y u z T
（18）
刚体的位移场完全由单个位移矢量定义。可变形体的位移场由其材料点的位移矢量集定义。
应力是作用于表面的单位面积力的度量。在固体结构中，外力或温度变化会产生可能导致固体结构运动和变形的应力场。而阻碍结构变形的内力也会引起可将结构恢复到原始未变形状态的内部弹性应力。一些材料还具有内置应力，这种应力在没有作用力和变形时就已存在。
一般情况下，作用于体的平面截面上的应力由矢量（称作应力矢量或拉力）定义如式（19）：
τ=τx τy τz =F x /A F y /A F z /A
（19）
式中，F 为作用于平面的力，A 为平面的面积。在任意点处，应力是单位面积上的力，应力状态完全
由与3个通过该点且相互垂直的平面关联的应力矢量定义。因此，任何平面上的任何点处的应力状态均通过式（20）形式的二阶张量定义：
σ=σxx σxy σxz
σyx σyy σyz
传祺gs6σzx σzy σzz
（20
）
43
发布评论