Tomohiko Hachiya，Yoshihiro Obata

Department of Mechanical and Aerospace Engineering，Graduate School of Engineering，Tottori University，Tottori 680－8552，Japan

1 Introduction

The importance of energy saving against the global warming and for the sustainable development has become widely understood.The newlybuilt buildings and houses have had the higher performance of heat insulating property and air tightness than older ones as the demand for airconditioning increases.On the other hand，the area of windows has not become smaller than older ones because of day lighting and designing against the idea of energy saving.As a result，the energy loss through windows is estimated to be 70%of total energy loss from houses.It is demanded strongly to reduce energy loss from windows with satisfying the demand for the lighting and designing.

The solar light reaching at the window glass surface reflects，is absorbed in the glass，or transmits through the glass into the room.Those optical properties of glass is defined as the trans－mittance，the absorptivity and reflectance，respectively.The summation of these three optical properties is constant of 1.Therefore increasing the reflectance and／or the absorptivity is needed to decrease the amount of heat－ray transmitting into the room.The heat－ray reflecting glass enhances a function of mirror at surface and is one of the effective candidates to decrease the transmitting heat－ray.But it has defects that the reflected light may blind the driver′s eye and that the mirror may raise radio disturbance for mobile phone because the glass reflects not only heat－ray but also light and radio wave.The heat－ray absorbing glass is another candidate to decrease the incoming heat－ray amount into the room.It contains dispersed metal particles which are absorbing heat－ray.But the heat－ray absorbing glass sometimes yields heat cracks by thermal stress.It is important to know the thermal stresses in the glass to use the heat－ray absorbing glass effectively.

We treat thermal stress in the heat－ray absorbing sheet glass in this paper.Although the design manual to avoid heat crack is open for builders and／or clients［1］，it refers to steady state and assumes that the absorbing of heat－ray occurs at only the glass surface.We introduce the idea of local absorptivity per unit length and propose the new model that the heat－ray is absorbed at not only surface but also all over the thickness.We also analyze the problem as the unsteady problem.

In addition，we consider the sudden weather changes such as sudden heavy rain and／or sudden strong wind after the glass is heated by the absorbed heat－ray and is reached to be steady state.The thick cloud may decrease the amount of absorbed heat－ray.The surface of the glass is suddenly cooled by rain and wind.Such sudden temperature change may cause heat cracking in the glass.These sudden weather changes are treated with the change of amount of absorbed heat－ray and／or the change of heat transfer coefficient between the glass surface and the outside atmosphere.

2 Analysis

2.1 Heat－ray absorbing process

We define that the heat－ray passing in the glass is the energy per unit time and per unit area，i.e.，it has the same units to the heat flux.We may treat the internal heat absorbing process as the steady state because the passing time in the sheet glass of the heat－ray with light speed is much smaller than the time that the heat transfers in the glass by the heat conduction.Therefore we assume that the heat－ray passing at the arbitrary position of xdepends on only the position xand not on the time.

Now we consider the energy balance in a tiny element as shown in Fig.1.The incoming heatray at the position xis I（x）；the outgoing heat flux at x＋dxis I（x）＋dI（x）.We introduce

the local absorptivity is constant because metal particles considered to be dispersed uniformly in sheet glass.Then the absorbed heat－ray in the tiny element is aIdx.

The energy balance gives the following equation

Integrating Eq.（1）gives Eq.（2）

Boundary condition is given as follows

where I0is the amount of the heat－ray reaching the glass surface.Applying Eq.（3）into Eq.（2）gives

We also define the total absorptivity for sheet glass with thickness b

Analyzing Eq.（5）gives the relationship between total absorptivity A and that per unit length aas follows

2.2 Temperature field

We consider a heat－ray absorbing sheet glass with thickness bas shown in Fig.2.The fine metal particles are dispersed uniformly and absorb heat－ray I0incomes at the surface of x＝0.Now we assume that the quantity of heat－ray passing through the cross－section at the arbitrary position of xis I（x），which is the function of position x.We treat the heat conduction and thermal stress in which the absorbed heat－ray becomes the internal heat generation as the unsteady problems.We also assume that the thermal boundary conditions

Fig.2 Coordinate system and thermal boundary conditions

are that heat convention occurs at x＝0and x＝b with heat transfer coefficient h1and h2between around air，respectively.

We assume that the heat－ray absorbed at dispersed metal becomes into the heat generation.Then the governing equation，boundary conditions and the initial condition are given as follows

舉報(bào)者為中行揚(yáng)州支行原員工王燕茹。她自稱(chēng)與黃宇曾交往多年，打算結(jié)婚。去年7月，黃宇因懷疑她與其他男人交往而對(duì)其進(jìn)行毆打?！拔覉?bào)警后才發(fā)現(xiàn)他已婚，后來(lái)去婚姻登記處詢(xún)問(wèn)，才知道他早就結(jié)婚了”。

Basic equation

whereρ，Candκis the density，the specifice heat and the thermal diffusivity，respectively.

Boundary conditions

Initial condition：heating

Initial condition：sudden weather change

where T0is the initial temperature，TMAXthe maximum temperature of the glass.

By solving Eq.（7）under the conditions Eqs.（8－11），we obtain the following temperature equation

whereβis coefficient of solar radiation and pmthe eigenvalues which satisfy the following equation

and

2.3 Thermal stress field

Let apply the result of temperature to the thermal stress of a plate.At the time t，assuming that the distribution of temperature is T（t，x），the strain and curvature at x＝0areε0and r0，respectively，the thermal stressσ（t，x）is given by

where Gis the modulus of transverse elasticity，α the coefficient of linear thermal expansion andν the Poisson′s ratio.The strain and curvature at x＝0，i.e.，ε0and r0，are determined to satisfy the mechanical boundary conditions if the elongation（E）and bending（B）are free（f）or clamped（c）.

3 Analysis Result and Consideration

3.1 Absorptivity

Table 1shows the total absorptivity of commercial heat－ray absorbing sheet glass［2］and the local absorptivity per unit length calculated by Eq.（6）.The total absorptivity depends on the thickness of sheet glasses and the thicker sheet glass has the larger total absorptivity.This result suggests that our modeling of internal heat absorbing process is reasonable.We use these absorptivities for the numerical calculation and refer other physical properties needed for the calculation to Ref.［2］.

Table 1 Absorptivity of commercial heat－ray absorbing sheet glass

3.2 Sudden change of weather

Fig.3shows unsteady temperature distribution in cases of a strong wind and rain.They are calculated for the sheet glass with 12mm nominal thickness，which is the maximum thickness among commercial heat－ray absorbing sheet glass.

In case of a strong wind，the temperature shows distribution to gradually fall.In case of rain，the temperature shows distribution to suddenly fall.However，the room side becomes the moderate temperature change.

Fig.3 Unsteady temperature distribution

Fig.4 Unsteady thermal stress distribution

Fig.4shows the distributions of unsteady thermal stress for temperature distributions in Fig.3.In case of a condition of EcBf，compress stress occurs in the room outside and tensile stress occurs in the room inside.The stress oc－curs in the tensile side because of a sudden temperature change.However，the stress becomes small gently because the room side is a moderate temperature change.In addition，it is unlikely to be the heat crack because the thermal stress is not exceeding allowable stress.

3.3 Soda glass and heat－ray absorbing glass

Next，we compare the common soda glass with heat－ray absorbing glass.Fig.5shows the distributions of maximum temperature of common soda glass and heat－ray absorbing glass.The absorption factor of common soda glass is a＝0.018mm－1.The glass does not absorb much heat，and the temperature does not rise，too.Therefore the thickness of soda glass is up to b＝19mm，but the thickness of heat－ray absorbing glass is up to b＝12mm.

Fig.5 Steady temperature distributions of different glasses

Fig.6 Unsteady thermal stress distributions of different glasses（EfBf）

Fig.6shows the distribution of unsteady thermal stress.Under the influence of rain，big tensile stress occurs temporarily.Big stress occurs in the heat－ray absorbing glass which is higher in initial temperature than common soda glass.However，it is unlikely to be the heat crack because the thermal stress is not exceeding allowable stress with both glasses either.

［1］ Nippon Sheet Grass Co Ltd.Heat crack in sheet glasses［EB／OL］.http：∥glass－catalog.jp／pdf／g12－010.pdf，2013.（in Japanese）

［2］ Asahi Glass Co Ltd.Thermal and optical properties of heat－ray reflecting sheet glasses［EB／OL］.http：∥www.asahiglassplaza.net／cata－logue／sougo＿sho2012／spdfdata／00342＿26s.pdf ＆00342＿27s.pdf，2013.（in Japanese）