基于Abaqus的塑料缸盖罩耐钢球撞击的仿真分析.docxVIP

论文所属行业：汽车

基于Abaqus的塑料缸盖罩耐钢球撞击的仿真分析

徐熹，曾庆强，贾正锋

(重庆长安汽车股份有限公司动力研究院重庆市401120)

摘要:钢球撞击试验是塑料缸盖罩开发过程中一项重要的单体试验，主要目的是模拟不明外物撞击塑料缸盖罩后，

缸盖罩是否能够不产生裂纹而继续正常工作。本文通过有限元的显示动力学分析方法，进行塑料缸盖罩耐钢球撞击

的仿真分析。

关键词:缸盖罩，碰撞，有限元，显示动力学，Abaqus

1.简介

在Abaqus中，Abaqus/Standard和Abaqus/Explicit都能分析多种类型的问题，应根据问题的特点和求解效率来选

择合适的分析类型。一般来说，对于光滑的非线性问题，Abaqus/Standard更有效，而Abaqus/Explicit适于求解复杂非

线性动力学问题，特别是用于模拟短暂、瞬时的动态问题，如冲击和爆炸问题。而且对于复杂接触问题，Abaqus/Explicit的求解效率优于Abaqus/Standard，包括模拟加工成形过程中接触条件改变的问题。

对于钢球撞击塑料缸盖罩仿真分析，属于碰撞分析，接触非线性、几何非线性和材料非线性三大非线性同时存

在，因此采用基于直接解法的显示动力学分析(explicitdynamicanalysis)，通过显式直接积分来求解该非线性动态问题。

2.理论基础

钢球撞击塑料缸盖罩仿真分析的一个重要特点是接触和冲击载荷影响碰撞的全过程。碰撞分析涉及到一个含未

知边界条件的偏微分方程求解问题，即动态接触问题。采用拉格朗日描述增量法，根据连续介质力学的质量、动量和能量守恒方程，碰撞体系的基本方程可以分别写成：

动量方程：

ρxi?=σij,j+ρfi公式（1）

质量方程：

ρ?+ρxi,j=0公式（2）

能量守恒方程：

E?=VS

ijεi?j?(p+q)V?公式（3）

上式中：xi表示质点坐标；xi?表示加速度；σij，εi?j，Sij分别表示柯西应力张量、应变率张量、偏应力张量；fi表

1

示单位质量体积力；V，ρ分别为当前构形的体积和质量密度；偏应力Sij=σij+(p+q)δij，p=σkk?q，q为体积

3

黏性阻力。根据最小势能原理，由上述基本方程，并运用d’Alembert原理和散度理论，可得接触系统的运动方程必须满足如下方程：

δH=∫Vρxi?δxidv+∫Vσijδi′jdv?∫Vρfiδxidv?∫Sgiδxids

1

=0公式（4）

式中：δxi满足位移边界条件；gi=gi(t)为作用在边界S1上的力；各项积分分别表示单位时间内体系的惯性力、内力和外力的虚功。在碰撞中惯性力起主导作用，体系的体力和表面力相对而言可以忽略不计。

3.钢球碰撞过程的模拟

3.1有限元模型

塑料缸盖罩的钢球撞击试验中，用质量为0.5kg的钢球在高度0.3m的位置处，自由落体撞击缸盖罩表面任意位置，可多次重复试验，要求撞击后缸盖罩不得出现裂纹。

根据试验过程，设计如下有限元模型：钢球初始位置位于缸盖罩上方0.005mm处；初始移动2424.9mm/s(根据

自由落体公式计算得到)，方向竖直向下。见图1。

图1.有限元模型示例图.

缸盖罩采用修正的二阶四面体C3D10M单元划分，钢球简化为离散刚体，其中缸盖罩和钢球可能接触区域网格需要细化。缸盖罩螺栓安装孔处全约束。

设置仿真分析总时间为4ms。由于缸盖罩和钢球的接触区域不可精确预测，接触设置为通用接触(general

contact)。

缸盖罩的主要材料参数见表1。在Abaqus的求解文件中，通过*Material下的*Damping命令定义材料的瑞利阻

尼(Rayleighdamping)，即阻尼矩阵C是质量矩阵M和刚度矩阵K的线性组合：

?=α?+β?公式（5）

表1.材料参数表.

部件材料弹性模量(GPa)泊松比(-)密度(g/cm3)

缸盖罩PA66-GF3570.361.4

3.2碰撞仿真的可信性评估

通过查看分析结果中的能量平衡，可以评估Abaqus/Explicit分析过程的正确性。整体模型的能量平衡可以表示

为公式（6）：

EI+EV+EFD+EKE?EW=Etotal公式（6）

式中：EI为内能；EV为粘性耗散能；EFD为摩擦耗散能；EKE为动能；EW为外载荷所做的功；Etotal为上述能量的总和，它必须保持为常数。

图2是整个碰撞过程中的能量变化曲线，其中ALLIE(InternalEnergy)为系统内能，ALLKE(K