建立了风栅中玻璃的冷却模型，数值模拟玻璃冷却的温度和应力变化规律，反演了不同季节风温时的合理匹配风压。结果表明，在玻璃淬冷过程，约3 s时玻璃表面拉应力达到最大，若该应力大于玻璃此时的抗拉强度，玻璃将破裂。此后玻璃从外到内降温速率逐渐减小，在约15~17 s时玻璃表层受内部影响减弱，表面应力趋于稳定。与钢化玻璃表面应力测试结果相比，数值模拟结果略小，但相对误差不超过5%。随冷却风温降低，玻璃钢化所需的风压逐渐减小。在玻璃钢化程度接近的情况下，风压随风温降低近似线性减小，钢化风压调节量与环境温度变化量的相关系数为0.103 kPa/K。
Wind temperature and wind pressure are important factors affecting tempering degree of glass. And the reasonable matching relationship should be sought between wind pressure and wind temperature due to seasonal variation, which satisfied the quality requirements of tempered glass and achieve the effect of energy conservation. In this work, the air-grid model of glass tempering was established to numerically simulate the change rules of glass temperature and stress during glass specimen process, and at the same time, the reasonable wind pressure was numerically calculated to match the wind temperatures in different season. The results showed that in the process of glass cooling, the glass surface cooling rate was fastest, and the internal cooling rate was slower. And the cooling rate of glass decreased gradually from outside to inside, and the influences of glass inner to surface weakened. But the curves of internal cooling rate and surface cooling rate reversed at about 15~17 s which meant that the internal cooling rate was larger, while the external cooling rate was smaller. The glass surface stress weakened by the internal temperature difference, in other words, the surface stress tended to be stable from 15 to 17 s in cooling process. At about 3 s, the tensile stress on the glass surface reached the maximum and was greater than the stress at all internal positions. When the stress was greater than the tensile strength of the glass at this time, the glass started to crack from the surface. Compared with the stress test results of the tempered glass, the numerical results were slightly smaller, but the relative error was not more than 5%. With the decrease of cooling air temperature, the wind pressure required for glass tempering was also decreasing. When the tempering degree of glass was close, the wind pressure decreased approximately linearly with the decrease of wind temperature.
[1]许伟光.钢化玻璃工艺参数的设定[J].玻璃与搪瓷, 2011, 39(4):25-29
[2]郭凯.从专利申请角度分析我国玻璃制品钢化技术发展[J].新型建筑材料, 2013, 40(7):91-94
[3]程金树，朱连英，楼贤春.玻璃钢化方法的探讨[J].材料导报, 2012, 26(20):135-137
[4]NIELSEN J H, OLESEN J F, POULSEN PN，et al.Simulation of residual stresses at holes in tempered glass: a parametric study[J].Materials and Structures, 2010, 43(7):947-961
[5]李会平，李冬生，王云正.玻璃钢化加热过程的初步分析[J].材料导报, 2012, 26(12):137-144
[6]赵德清，唐仪，黄志刚.基于的钢化玻璃淬冷系统流场的建模与仿真[J].机电工程技术, 2016, 45(4):1-6
[7]郭卫，童树庭，朱雷波.玻璃钢化工艺过程与钢化应力的研究[J].建筑材料学报, 2005, 8(1):100-104
[8]朱雷波，童树庭，郭卫，等.玻璃钢化工艺的新进展[J].玻璃与搪瓷, 2002, 30(1):39-42
[9]Abena H, Antona J, Paemurrub M, et al.A new method for tempering stress measurement in glass panels[J].Estonian Journal of Engineering, 2013, 19(4):292-297
[10]封福明，冯悦冲，许世清，等.浮法玻璃在退火窑 - 区内退火过程的数值模拟研究[J].燕山大学学报, 2017, 41(4):335-342
[11]张景超，李立，丁瑞，等.玻璃中永久应力的理论分析与数值模拟[J].材料科学与工程学报, 2011, 29(5):790-793
[12]张范斌，温治，苏福永，等.浮法玻璃退火过程中应力变化的数值模拟[J].硅酸盐通报, 2012, 31(5):1057-1061
[13] 李潮洋.基于 ANSYS 玻璃退火过程中应力的模拟分析[D]. 燕山大学，2013.
[14] 唐秋华，吴再豪，余震，等.玻璃钢化冷却过程的建模与仿真[J][J].机械设计与制造, 2014, (12):152-158
[15] 翟守元.物理钢化玻璃内应力的分析[J][J].山东轻工业学院学报, 1987, (1):17-20
[16] 吴再豪.玻璃钢化冷却过程预控模型与仿真[D]. 武汉科技大学，2015.
[17]夏国华，朱锦杰，童树庭.平板玻璃钢化有限元数值模拟[J].玻璃与搪瓷, 2008, 36(5):1-5
[18] Bernard F, Gy R, Daudeville L.Finite element computation of residual stresses near holes in tempered glass plates[J].[J].Glass Technology, 2002, (43c):290-295
[19] Aronen A.Modelling of Deformations and Stresses in Glass Tempering[D]. Tampere University of Technology, 2012.
[20] Benbahouche1 S, Roumili1 F, Sangleboeuf J C.Modelling of Internal Stresses in Sheet Glass during the Cooling[J].[J].World Journal of Nano Science and Engineering, 2011, (1):21-26
[21] Fan BaoWei, Zhu KeQian, Shi Qiang, et al.Effect of glass thickness on temperature gradient and stress distribution during glass tempering[J][J].Journal of Non-Crystalline Solids, 2016, (437):72-79
[22]黄谦，桂业伟，耿湘人.层状平板内热应力的计算研究[J].工程热物理学报, 2005, 26(3):492-494
[23]ZHOU Huamin，XI Guodong，LI Dequn.Modeling and simulation of residual stresses during glass bulb pressing process[J].Science in China Series E: Technological Sciences, 2007, 50(1):103-117
[24]IRZA E M.Mathematical model of the formation of residual stresses in glass bodies of revolution in the process of cooling[J].Materials Science, 2008, 44(2):156-162
[25]AJOVALASIT A，PETRUCCI G，SCAFIDI M.Photo elastic analysis of edge residual stresses in glass by automated “test fringes” methods[J].Experimental Mechanics, 2012, 52(8):1057-1066
