摘要
玻璃由液态向固态转变的过程中，内部的应力变化情况非常复杂，而对形成玻璃暂时应力和永久应力的影响因素进行分析，有利于减少成型过程中玻璃的炸裂情况和改善成型后玻璃内的应力分布情况，实现玻璃高效率、高质量的生产。本文模拟计算了自由边界和约束边界两种情况下，加热平板玻璃一边时应力的分布情况，有助于找到在发生火灾的情况下玻璃应力最大处的位置，即最容易炸裂的位置。
首先，基于ANSYS模拟软件，通过建模，参数设置，建立玻璃的二维模型。设置玻璃退火的曲线，模拟计算玻璃的表面层和中间层在退火过程中的应力，用MATLAB进行数据处理，得到它们随退火时间及随表面层与中间层温差的变化规律。
其次，通过改变B区的退火速度，得到不同速度下所形成的永久应力，用MATLAB对得到的结果进行处理，得到永久应力与B区退火速度的关系，并与理论计算值比较是否相符。进一步用相同方法得到C区的暂时应力随温差的变化关系。
最后，建立自由边界和约束边界两种条件下玻璃被加热的模型，分析它们X、Y方向应力随加热过程应力的变化。分析自由边界和约束边界时玻璃内X方向应力的性质，得到X方向应力随玻璃位置的变化规律。相同方法分析自由边界和约束边界条件下Y方向应力情况，得到Y方向应力随玻璃位置的变化规律。
The internal stress is very complex, when glass changes into solid-state from liquid.The analysis of the factors which affect the temporary and permanent stress is work toreduce the burst of glass in modeling process and improve the stress’distribution in glass.It helps to improve efficiency and quality. The stress’distribution have been simulatedunder Free and constraint boundaries conditions in this essay, which helps to find themaximum stress at the glass location where is most likely to burst.
Firstly, the two-dimensional glass model has been established based on ANSYS,through modeling and parameter setting. Setting the glass annealing progress, middle andsurface layers’stress have been simulated. Processing the data with MATLAB, The stressvariation changing with the temperature difference between middle and surface layers, hasbeen given.
Secondly, by changing the annealing rates of B area, the permanent stress at differentrates has been obtained. Processing the data with MATLAB, the relationship betweenpermanent stress and B area annealing rates is gotten, which is then compared with thetheoretical calculations. Also the same is C area.
Finally, the heat models of free and constraints boundaries have been established.X-direction and Y-directions stress, changing with the heating process, have been analyzed.The nature of stress in the X-direction, under the free and constraint boundaries wascompared. Also, the X-direction stress changing with location of the glass is gotten. Thenature of stress in the Y-direction, under the free and constraint boundaries was alsocompared. The Y-direction stress changing with the location of the glass is gotten, too.
引文
[1]罗漪青.浮法玻璃退火窑的研究与设计[D].成都:成都理工大学, 2009: 1-2.
[2] T. A. Egorova, V. Y. Matveenko, V. I. Kondrashov, et al. Production of Floatglass of More thanEquilibrium Thickness [J]. Glass and Ceramics. 2002, 59(2): 265-366.
[3] M. Pawluk-Ko c, J. Zieba-Palus, A. Parczewski. The effect of re-annealing on the distribution ofrefractive index in a windscreen and a windowpane Classification of glass samples [J]. ForensicScience International, 2008, (174): 222-228.
[4] Oriel Goodman, Brian Derby. The mechanical properties of float glass surfaces measured bynanOindentation and acoustic microscopy [J]. Acta Materialia, 2011, (59): 1790-1799.
[5] K. P. Rushton, S. A. Coulson, A. W. N. Newton. The effect of annealing on the variation of glassrefractive index [J]. Forensic Science International, 2011, (209): 102-107.
[6]冯跃冲.浮法玻璃退火过程及其应力的数值模拟[D].秦皇岛:燕山大学, 2006: 1-2.
[7]孙晋涛.硅酸盐工业热工基础[M].武汉:武汉工业出版社, 1994:123-141.
[8] J. Klouzek. Modeling of the refining space working under reduced pressure glasstech ber [J].Glass science Technol, 2000, 73(11): 329-336.
[9]李立.浮法玻璃退火过程中应力的理论研究[D].秦皇岛:燕山大学, 2011: 1-2.
[10] Shen Junwu, D. J. Green, R. E. Tressler. Stress relaxation of a soda lime silicate glass below theglass transition temperature [J]. Journal of Non-Crystalline Solids, 2003, (324): 277-288.
[11]叶先磊. ANSYS工程分析软件应用实例[M].北京:科学出版社, 2003: 76-84.
[12]刘涛,杨凤鹏.精通ANSYS[M].北京:清华出版社, 2002: 2-10, 406-418.
[13] J. Wang. Validation of an improved batch model in a coupled combustion space/melt tank/batchmelting glass furnace simulation glasstech.ber [J]. Glass science Technol, 2000, 73(10): 299-307.
[14] S. D. Cesar, J. M. A. Numerical modeling of glass forming processes [J]. Engineering computerpine ridge press, 1986, 86(3): 226-274.
[15]白葳,喻海良.通用有限元分析ANSYS8.0基础教程[M].北京:清华大学出版社, 2005: 1-6.
[16]邵宏根.浮法玻璃成型的数值模拟与实验研究[D].秦皇岛:燕山大学, 2004: 7-8.
[17] LAB Pinkington. The float glass process [J]. ProcRoyLond, 1969, (314): 1-25.
[18] H. Fischer, M. Hemelik, R. Telle. Influence of annealing temperature on the strength of dentalglass ceramic materials [J]. Dental Materials, 2005, (21): 671-677.
[19] M. Aljerf, K. Georgarakis, A. R. Yavari. Shaping of metallic glasses by tress-annealing withoutthermal embrittlement [J]. Acta Materialia, 2011, (59): 3817-3824.
[20]许超,程继健等.玻璃工艺原理[M].北京:中国建筑工业出版社, 1984: 42-45.
[21]林亢.普通工业玻璃退火与应力的几个问题[J].玻璃,2000,27(1): 11-33.
[22]林亢.普通工业玻璃退火与应力的几个问题[J].玻璃,2000,27(2): 10-14.
[23] J. Malzbender. G. De. With cracking and residual stress in hybrid coatings on floating glass [J].Thin solid films, 2000, 359(4): 210-214.
[24] F. W. Prestion. Annealing of flat glass [J]. Glass Tech, 1952, (36): 287.
[25]林亢.普通工业玻璃退火与应力的几个问题[J].玻璃,2000, 27(3): 7-10.
[26] M. Vaniseghem. Simulation of material processing: theory, methods, and application [J]. BalkemaRotterdam, 1998, 113(5): 547-552.
[27]干福喜.光学玻璃[M].北京:科学出版社, 1964:354.
[28]作花济夫,等.玻璃手册[M].北京:建筑工业出版社, 1985:345.
[29]陈恭源.玻璃的退火[J].玻璃, 2005, (2): 56-58.
[30]王传辉.浮法玻璃退火的热应力及其分析[J].洛阳工业高等专科学校校报, 1996, 6(3): 6-10.
[31]韩彤.浮法玻璃的应力与退火[J].玻璃, 2011, (4): 11-14.
[32]杜永明.平板玻璃表层应力在线检测[J].中国玻璃, 1998, (5): 35-37.
[33]陈恭源.玻璃退火和浮法玻璃退火理论[J].玻璃, 2010, (2): 3-15.
[34]冯毅.真空绝热玻璃应力理论分析及数值计算[J].低温与超导, 2009, 37(7): 78-80.
[35]周天辉.玻璃应力的测定方法[J].玻璃与搪瓷, 2001, 29(6): 44-48.
[36]吕岑.玻璃应力图像处理测量方法研究[J].陕西科技大学学报, 2005, 6(23): 65-68.
[37]高佳栋.应力分布的光学图像处理研究[J].玻璃与搪瓷, 2007, 35(4): 34-38.
[38]吴晓.浮法玻璃在线应力检测及相关退火概念[J].玻璃, 2005, (1): 18-22.
[39]王桂英.浮法玻璃应力检测系统及其应用[J].玻璃, 2011, (5): 25-27.
[40] Zheng Hongzhi. Temperature adaptability of magnetooptic modulation in a disconnectmechanically azinuth measurement [J]. Acta Photonica Sinica, 2004, 33(5): 638-640.
[41] Gilbert Lebrun, Johanne Denault. Effect of annealing on the thermal expansion and residualstresses of bidirectional thermoplastic composite laminates [J]. Composites, 2010, (41): 101-107.
[42]吴易明.一种玻璃材料内应力精密测定的方法[J].光子学报, 2010, 39(3): 490-493.
[43]张志勇.几种波片位相延迟测量方法的比较[J].光学精密工程，2007, 15(11): 1678-1686.
[44] W. Feng. A Spectroscopic Method for Determining Thickness of Quartz Wave Plate[J]. ChineseOptics Letters, 2006，4(12): 705-708.
[45] J. Klouzek. Modeling of the Refining Space Working under Reduced Pressure [J]. Glastech. Ber.GlassSci. Technol, 2000, 73(11): 329-336.
[46] J. Bauer. Validation of a Mathematical Glass Tank Model [J]. Glastech. Ber. GlassSci. Technol,1999, 72(6): 71-181.
[47] W. Wang. Validation of an Improved Btch Model in a Coupled Combustion Space /Melt Tank/Batch Melting Glass Furnace Simulation [J]. Glastech. Ber. GlassSci. Technol, 2000, 73(10):301-307.
[48]韩文梅.基于ANSYS航空层合玻璃的热应力分析[J].兵器材料科学与工程, 2010, 33(6): 1-3.
[49]满金仓.浮法在线Low-E玻璃退火窑的设计与控制[J].玻璃与搪瓷, 2008, 36(6): 20-23.
[50] Wang Qingsong, Zhang Yi, Sun Jinhua. Temperature and thermal stress simulation of windowglass exposed to fire [J]. Procedia Engineering, 2011, (11): 452-460.
