汽车玻璃钢化分析.docx

汽车玻璃钢化分析

上海耀皮康桥汽车玻璃有限公司王旸翁洋

摘要：应用有限元分析软件Abaqus对某车型汽车的车门玻璃钢化效果进行分析，直观而且准确的预测各种工艺参数的变化

对钢化玻璃产生的影响，并且通过更改工艺条件来消除其对产品产生的不良影响。

关键词：Abaqus，钢化，热应力

1.引言

现在随着汽车行业的发展，汽车玻璃的外形变得越来越复杂，对玻璃的成型，质量以及成本都提出了更高的

要求。由于模拟技术应用于产品设计开发阶段，可以在设计初期发现并更改设计和加工中的缺陷，在生产首次试

件时回避掉绝大部分缺陷，极大的节省了工艺开发的时间以及费用。

本文所研究的模拟分析是针对玻璃快速冷却至室温时表面形成压应力层分布，在模拟过程中通过应用流体力

学来模拟空气与玻璃的热转换数值，应用固体力学来模拟分析玻璃的应力变化。数值模型的建立基于玻璃的粘弹

性理论，玻璃结构的松弛性也被考虑在内。模拟可以使我们预知在指定条件下玻璃的温度分布、瞬时应力，以及

变形情况。

淬冷时玻璃表面的压应力与中心张应力可按下式计算：

σt=αE/(1-υ)×1/3(ΔTmax)=αE/(1-υ)×bQ/24K

σc=2σt=αE/(1-υ)×bQ/12K

式中σt－－玻璃中心张应力；

σc－－玻璃表面压应力；

α－－玻璃的膨胀系数；

E－－玻璃的弹性模量；

υ－－玻璃的泊松比；

ΔTmax－－玻璃厚度中心与表面最大温差；

b－－玻璃厚度；

Q－－玻璃单位面积上均匀放出的热量；

K－－玻璃的热导率。

从这两个公式我们也不难得出理论上表面压应力的最大值近似为中心张应力最大值的2倍。都与玻璃的材料

参数，温差，玻璃厚度以及玻璃的热学参数密切相关。

2.计算模型和边界条件的处理

2.1有限元模型

玻璃的有限元模型分别采用DC3D20（热传分析）及C3D20R（应力分析）网格划分。由于比较关心玻璃厚

度方向上的残余应力分布情况，所以在厚度方向上划分较密的网格，以便于输出结果。

玻璃的有限元模型

2.2材料数据

在低温和常温下，玻璃基本上是服从虎克定律的弹性体，而在高温下则表现出粘弹性力学性能（松弛和蠕变）。

下面这个表格列出了部分在分析中定义的性能参数。

杨氏模量E（Pa）7.138E10

泊松比ν0.218

密度ρ（kg/m3）2520

2.3边界条件

淬冷过程是将加热后的玻璃迅速冷却使其达到室温时表面形成高度均匀的压应力层。整个过程的传热方式包

括冷空气与玻璃上下表面的对流换热，玻璃内部的热传导以及玻璃向外界环境的热辐射。这里起主导作用的方式

为对流换热，因此忽略辐射的影响。

模拟分析分两个步骤进行，分别为温度响应分析和应力响应分析。

3计算结果

3.1温度响应分析结果

下图分别为玻璃表面及厚度方向上的温度变化。

玻璃表面温度分布

玻璃厚度方向上温度变化分布

从结果数值上可以看出在玻璃表面与中心最大的温差达到80度左右，此时的温度分布，我们认为是理想应

力生成的主要区间。

3.2应力响应分析结果

钢化处理在冷却阶段会产生残余应力。理论上在厚度方向上残余应力的分布如下：

表面压应力可以改善弯曲强度以及抗撞击强度，在破碎测试中中间面的张应力直接关系到碎片的密度。

下图为玻璃中间部分及侧边在厚度方向上残余正应力分布。

应力在整个厚度方向为抛物线分布，这也验证了前面提到的理论。理论与实验都可得出中间面的张应力近似

为表面压应力的一半并且在破碎测试中得出小颗粒的密度直接关系到中间面的张应力。另外通过分析对比得出残

余应力的大小主要取决于玻璃的对流换热系数以及初始温度。

下图分别为玻璃表面及厚度方向上残余正应力的分布。

玻璃表面残余应力分布玻璃厚度方向残余应力分布

计算的结果显示在玻璃边缘部分存在压应力的峰值，并且在玻璃的圆角处有小的张应力。针对这个问题我们

可以将玻璃圆角处模拟更加详细的边界条件。因为汽车玻璃的边部与车壳接触，因此边部应力尤为重要。为了降

低边部应力，可以通过在钢化过程中减少边部地区的热传递来解决。

4总结

汽车玻璃的钢化模拟可以评估钢化方式以及得出玻璃的残余应力值。它以一种理想方式来分析各种工艺参数

对钢化结果的灵敏度。在这里我们得到了温度以及应力的分布情况，然而边部应力部分以及各种工艺参数的优化

将是我们研究的方向。

通过使用数值模拟，我们能够执行一个瞬时分析以及建立一个虚拟的钢化实验来缩短产品设计开发阶段并且

为产品质量的改进提供了评鉴的信息。

参考文献

1.AbaqusVersion6.7Documentation[M].Abaqus,INC,2007

2.张义同.热粘弹性理论.天津大学出版社