Wednesday, 16 March 2011 21:45

## The Physical Basis of Work in Heat

Rate this item

Thermal Exchanges

The human body exchanges heat with its environment by various pathways: conduction across the surfaces in contact with it, convection and evaporation with the ambient air, and radiation with the neighbouring surfaces.

Conduction

Conduction is the transmission of heat between two solids in contact. Such exchanges are observed between the skin and clothing, footwear, pressure points (seat, handles), tools and so on. In practice, in the mathematical calculation of thermal balance, this heat flow by conduction is approximated indirectly as a quantity equal to the heat flow by convection and radiation which would take place if these surfaces were not in contact with other materials.

Convection

Convection is the transfer of heat between the skin and the air surrounding it. If the skin temperature, tsk, in units of degrees Celsius (°C), is higher than the air temperature (ta), the air in contact with the skin is heated and consequently rises. Air circulation, known as natural convection, is thus established at the surface of the body. This exchange becomes greater if the ambient air passes over the skin at a certain speed: the convection becomes forced. The heat flow exchanged by convection, C, in units of watts per square metre (W/m2), can be estimated by:

C = hc FclC (tsk - ta)

where hc is the coefficient of convection (W/°C m2), which is a function of the difference between tsk and ta in the case of natural convection, and of the air velocity Va (in m/s) in forced convection; FclC is the factor by which clothing reduces convection heat exchange.

Every body emits electromagnetic radiation, the intensity of which is a function of the fourth power of its absolute temperature T (in degrees Kelvin—K). The skin, whose temperature may be between 30 and 35°C (303 and 308K), emits such radiation, which is in the infrared zone. Moreover, it receives the radiation emitted by neighbouring surfaces. The thermal flow exchanged by radiation, R (in W/m2), between the body and its surroundings may be described by the following expression:

where:

s is the universal constant of radiation (5.67 × 10-8 W/m2 K4)

e is the emissivity of the skin, which, for infrared radiation, is equal to 0.97 and independent of the wavelength, and for solar radiation is about 0.5 for the skin of a White subject and 0.85 for the skin of a Black subject

AR/AD is the fraction of the body surface taking part in the ex- changes, which is of the order of 0.66, 0.70 or 0.77, depending upon whether the subject is crouching, seated or standing

FclR is the factor by which clothing reduces radiation heat exchange

Tsk (in K) is the mean skin temperature

Tr (in K) is the mean radiant temperature of the environment —that is, the uniform temperature of a black mat sphere of large diameter that would surround the subject and would exchange with it the same quantity of heat as the real environment.

This expression may be replaced by a simplified equation of the same type as that for exchanges by convection:

R = hr (AR/AD) FclR (tsk - tr)

where hr is the coefficient of exchange by radiation (W/°C m2).

Evaporation

Every wet surface has on it a layer of air saturated with water vapour. If the atmosphere itself is not saturated, the vapour diffuses from this layer towards the atmosphere. The layer then tends to be regenerated by drawing on the heat of evaporation (0.674 Watt hour per gram of water) at the wet surface, which cools. If the skin is entirely covered with sweat, evaporation is maximal (Emax) and depends only on the ambient conditions, according to the following expression:

Emax = he Fpcl (Psk,s - Pa)

where:

he is the coefficient of exchange by evaporation (W/m2kPa)

Psk,s is the saturated pressure of water vapour at the temperature of the skin (expressed in kPa)

Pa is the ambient partial pressure of water vapour (expressed in kPa)

Fpcl is the factor of reduction of exchanges by evaporation due to clothing.

Thermal insulation of clothing

A correction factor operates in the calculation of heat flow by convection, radiation and evaporation so as to take account of clothing. In the case of cotton clothing, the two reduction factors FclC and FclR may be determined by:

Fcl = 1/(1+(hc+hr)Icl)

where:

hc is the coefficient of exchange by convection

hr is the coefficient of exchange by radiation

Icl is the effective thermal isolation (m2/W) of clothing.

As regards the reduction of heat transfer by evaporation, the correction factor Fpcl is given by the following expression:

Fpcl = 1/(1+2.22hc Icl)

The thermal insulation of the clothing Icl is expressed in m2/W or in clo. An insulation of 1 clo corresponds to 0.155 m2/W and is provided, for example, by normal town wear (shirt, tie, trousers, jacket, etc.).

ISO standard 9920 (1994) gives the thermal insulation provided by different combinations of clothing. In the case of special protective clothing that reflects heat or limits permeability to vapour under conditions of heat exposure, or absorbs and insulates under conditions of cold stress, individual correction factors must be used. To date, however, the problem remains poorly understood and the mathematical predictions remain very approximate.

Evaluation of the Basic Parameters of the Work Situation

As seen above, thermal exchanges by convection, radiation and evaporation are a function of four climatic parameters—the air temperature ta in °C, the humidity of the air expressed by its partial vapour pressure Pa in kPa, the mean radiant temperature tr in °C, and the air velocity Va in m/s. The appliances and methods for measuring these physical parameters of the environment are the subject of ISO standard 7726 (1985), which describes the different types of sensor to use, specifies their range of measurement and their accuracy, and recommends certain measurement procedures. This section summarizes part of the data of that standard, with particular reference to the conditions of use of the most common appliances and apparatus.

Air temperature

The air temperature (ta) must be measured independent of any thermal radiation; the accuracy of the measurement should be ±0.2ºC within the range of 10 to 30ºC, and ±0.5 °C outside that range.

There are numerous types of thermometers on the market. Mercury thermometers are the most common. Their advantage is accuracy, provided that they have been correctly calibrated originally. Their main disadvantages are their lengthy response time and lack of automatic recording ability. Electronic thermometers, on the other hand, generally have a very short response time (5 s to 1 min) but may have calibration problems.

Whatever the type of thermometer, the sensor must be protected against radiation. This is generally ensured by a hollow cylinder of shiny aluminium surrounding the sensor. Such protection is ensured by the psychrometer, which will be mentioned in the next section.

Partial pressure of water vapour

The humidity of the air may be characterized in four different ways:

1.      the dewpoint temperature: the temperature to which the air must be cooled to become saturated with humidity (td, °C)

2.      the partial pressure of water vapour: the fraction of atmospheric pressure due to water vapour (Pa, kPa)

3.      the relative humidity (RH), which is given by the expression:

RH = 100·Pa/PS,ta

where PS,ta is the saturated vapour pressure associated with the air temperature

4.      the wet bulb temperature (tw), which is the lowest temperature attained by a wet sleeve protected against radiation and ventilated at more than 2 m/s by the ambient air.

All these values are connected mathematically.

The saturated water vapour pressure PS,t at any temperature t is given by:

while the partial pressure of water vapour is connected to the temperature by:

Pa = PS,tw - (ta - tw)/15

where PS,tw is the saturated vapour pressure at the wet bulb temperature.

The psychrometric diagram (figure 1) allows all these values to be combined. It comprises:

Figure 1.  Psychrometric diagram.

• in the y axis, the scale of partial pressure of water vapour Pa, expressed in kPa
• in the x axis, the scale of air temperature
• the curves of constant relative humidity
• the oblique straight lines of constant wet bulb temperature.
• The parameters of humidity most often used in practice are:
• the relative humidity, measured by means of hygrometers or more specialized electronic appliances
• the wet bulb temperature, measured by means of the psychrometer; from this is derived the partial pressure of water vapour, which is the parameter most used in analysing thermal balance

The range of measurement and the accuracy recommended are 0.5 to 6 kPa and ±0.15 kPa. For measurement of the wet bulb temperature, the range extends from 0 to 36ºC, with an accuracy identical with that of the air temperature. As regards hygrometers for measuring relative humidity, the range extends from 0 to 100%, with an accuracy of ±5%.

The mean radiant temperature (tr) has been defined previously; it can be determined in three different ways:

1.      from the temperature measured by the black sphere thermometer

2.      from the plane radiant temperatures measured along three perpendicular axes

3.      by calculation, integrating the effects of the different sources of radiation.

Only the first technique will be reviewed here.

The black sphere thermometer consists of a thermal probe, the sensitive element of which is placed at the centre of a completely closed sphere, made of a metal that is a good conductor of heat (copper) and painted matt black so as to have a coefficient of absorption in the infrared zone close to 1.0. The sphere is positioned in the workplace and subjected to exchanges by convection and radiation. The temperature of the globe (tg) then depends on the mean radiant temperature, the air temperature and the air velocity.

For a standard black globe 15 cm in diameter, the mean temperature of radiation can be calculated from the temperature of the globe on the basis of the following expression:

In practice, the need must be stressed to maintain the emissivity of the globe close to 1.0 by carefully repainting it matt black.

The main limitation of this type of globe is its long response time (of the order of 20 to 30 min, depending on the type of globe used and the ambient conditions). The measurement is valid only if the conditions of radiation are constant during this period of time, and this is not always the case in an industrial setting; the measurement is then inaccurate. These response times apply to globes 15 cm in diameter, using ordinary mercury thermometers. They are shorter if sensors of smaller thermal capacity are used or if the diameter of the globe is reduced. The equation above must therefore be modified to take account of this difference in diameter.

The WBGT index makes direct use of the temperature of the black globe. It is then essential to use a globe 15 cm in diameter. On the other hand, other indices make use of the mean radiant temperature. A smaller globe can then be selected to reduce the response time, provided that the equation above is modified to take account of it. ISO standard 7726 (1985) allows for an accuracy of ±2ºC in the measurement of tr between 10 and 40ºC, and ±5ºC outside that range.

Air velocity

The air velocity must be measured disregarding the direction of air flow. Otherwise, the measurement must be undertaken in three perpendicular axes (x, y and z) and the global velocity calculated by vectorial summation:

The range of measurements recommended by ISO standard 7726 extends from 0.05 to 2 m/s The accuracy required is 5%. It should be measured as a 1- or 3-min average value.

There are two categories of appliances for measuring air velo-city: anemometers with vanes, and thermal anemometers.

Vane anemometers

The measurement is carried out by counting the number of turns made by the vanes during a certain period of time. In this way the mean velocity during that period of time is obtained in a discontinuous manner. These anemometers have two main disadvantages:

1. They are very directional and have to be oriented strictly in the direction of the air flow. When this is vague or unknown, measurements have to be taken in three directions at right angles.
2. The range of measurement extends from about 0.3 m/s to 10 m/s. This limitation to low velocities is important when, for instance, it is a question of analysing a thermal comfort situation where it is generally recommended that a velocity of 0.25 m/s should not be exceeded. Although the range of measurement can extend beyond 10 m/s, it hardly falls below 0.3 or even 0.5 m/s, which greatly limits the possibilities of use in environments near to comfort, where the maximum permitted velocities are 0.5 or even 0.25 m/s.

Hot-wire anemometers

These appliances are in fact complementary to vane anemometers in the sense that their dynamic range extends essentially from 0 to 1 m/s. They are appliances giving an instantaneous estimate of speed at one point of space: it is therefore necessary to use mean values in time and space. These appliances are also often very directional, and the remarks above also apply. Finally, the measurement is correct only from the moment when the temperature of the appliance has reached that of the environment to be evaluated.

Back

" DISCLAIMER: The ILO does not take responsibility for content presented on this web portal that is presented in any language other than English, which is the language used for the initial production and peer-review of original content. Certain statistics have not been updated since the production of the 4th edition of the Encyclopaedia (1998)."

## Heat and Cold References

ACGIH (American Conference of Governmental Industrial Hygienists). 1990. Threshold Limit Values and Biological Exposure Indices for 1989–1990. New York: ACGIH.

—. 1992. Cold stress. In Threshold Limit Values for Physical Agents in the Work Environment. New York: ACGIH.

Bedford, T. 1940. Environmental warmth and its measurement. Medical Research Memorandum No. 17. London: Her Majesty’s Stationery Office.

Belding, HS and TF Hatch. 1955. Index for evaluating heat stress in terms of resulting physiological strain. Heating Piping Air Condit 27:129–136.

Bittel, JHM. 1987. Heat debt as an index for cold adaptation in men. J Appl Physiol 62(4):1627–1634.

Bittel, JHM, C Nonotte-Varly, GH Livecchi-Gonnot, GLM Savourey and AM Hanniquet. 1988. Physical fitness and thermoregulatory reactions in a cold environment in men. J Appl Physiol 65:1984-1989.

Bittel, JHM, GH Livecchi-Gonnot, AM Hanniquet and JL Etienne. 1989. Thermal changes observed before and after J.L. Etienne’s journey to the North Pole. Eur J Appl Physiol 58:646–651.

Bligh, J and KG Johnson. 1973. Glossary of terms for thermal physiology. J Appl Physiol 35(6):941–961.

Botsford, JH. 1971. A wet globe thermometer for environmental heat measurement. Am Ind Hyg J 32:1–10.

Boutelier, C. 1979. Survie et protection des équipages en cas d’immersion accidentelle en eau froide. Neuilly-sur-Seine: AGARD A.G. 211.

Brouha, L. 1960. Physiology in Industry. New York: Pergamon Press.

Burton, AC and OG Edholm. 1955. Man in a Cold Environment. London: Edward Arnold.

Chen, F, H Nilsson and RI Holmér. 1994. Cooling responses of finger pad in contact with an aluminum surface. Am Ind Hyg Assoc J 55(3):218-22.

Comité Européen de Normalisation (CEN). 1992. EN 344. Protective Clothing Against Cold. Brussels: CEN.

—. 1993. EN 511. Protective Gloves Against Cold. Brussels: CEN.

Commission of the European Communities (CEC). 1988. Proceedings of a seminar on heat stress indices. Luxembourg: CEC, Health and Safety Directorate.

Daanen, HAM. 1993. Deterioration of manual performance in cold and windy conditions. AGARD, NATO, CP-540.

Dasler, AR. 1974. Ventilation and thermal stress, ashore and afloat. In Chapter 3, Manual of Naval Preventive Medicine. Washington, DC: Navy Department, Bureau of Medicine and Surgery.

—. 1977. Heat stress, work functions and physiological heat exposure limits in man. In Thermal Analysis—Human Comfort—Indoor Environments. NBS Special Publication 491. Washington, DC: US Department of Commerce.

Deutsches Institut für Normierung (DIN) 7943-2. 1992. Schlafsacke, Thermophysiologische Prufung. Berlin: DIN.

Dubois, D and EF Dubois. 1916. Clinical calorimetry X: A formula to estimate the appropiate surface area if height and weight be known. Arch Int Med 17:863–871.

Eagan, CJ. 1963. Introduction and terminology. Fed Proc 22:930–933.

Edwards, JSA, DE Roberts, and SH Mutter. 1992. Relations for use in a cold environment. J Wildlife Med 3:27–47.

Enander, A. 1987. Sensory reactions and performance in moderate cold. Doctoral thesis. Solna: National Institute of Occupational Health.

Fuller, FH and L Brouha. 1966. New engineering methods for evaluating the job environment. ASHRAE J 8(1):39–52.

Fuller, FH and PE Smith. 1980. The effectiveness of preventive work procedures in a hot workshop. In FN Dukes-Dobos and A Henschel (eds.). Proceedings of a NIOSH Workshop on Recommended Heat Stress Standards. Washington DC: DHSS (NIOSH) publication No. 81-108.

—. 1981. Evaluation of heat stress in a hot workshop by physiological measurements. Am Ind Hyg Assoc J 42:32–37.

Gagge, AP, AP Fobelets and LG Berglund. 1986. A standard predictive index of human response to the thermal environment. ASHRAE Trans 92:709–731.

Gisolfi, CV and CB Wenger. 1984. Temperature regulation during exercise: Old concepts, new ideas. Exercise Sport Sci Rev 12:339–372.

Givoni, B. 1963. A new method for evaluating industrial heat exposure and maximum permissible work load. Paper submitted to the International Biometeorological Congress in Paris, France, September 1963.

—. 1976. Man, Climate and Architecture, 2nd ed. London: Applied Science.

Givoni, B and RF Goldman. 1972. Predicting rectal temperature response to work, environment and clothing. J Appl Physiol 2(6):812–822.

—. 1973. Predicting heart rate response to work, environment and clothing. J Appl Physiol 34(2):201–204.

Goldman, RF. 1988. Standards for human exposure to heat. In Environmental Ergonomics, edited by IB Mekjavic, EW Banister and JB Morrison. London: Taylor & Francis.

Hales, JRS and DAB Richards. 1987. Heat Stress. Amsterdam, New York: Oxford Excerpta Medica.

Hammel, HT. 1963. Summary of comparative thermal patterns in man. Fed Proc 22:846–847.

Havenith, G, R Heus and WA Lotens. 1990. Clothing ventilation, vapour resistance and permeability index: Changes due to posture, movement and wind. Ergonomics 33:989–1005.

Hayes. 1988. In Environmental Ergonomics, edited by IB Mekjavic, EW Banister and JB Morrison. London: Taylor & Francis.

Holmér, I. 1988. Assessment of cold stress in terms of required clothing insulation—IREQ. Int J Ind Erg 3:159–166.

—. 1993. Work in the cold. Review of methods for assessment of cold stress. Int Arch Occ Env Health 65:147–155.

—. 1994. Cold stress: Part 1—Guidelines for the practitioner. Int J Ind Erg 14:1–10.

—. 1994. Cold stress: Part 2—The scientific basis (knowledge base) for the guide. Int J Ind Erg 14:1–9.

Houghton, FC and CP Yagoglou. 1923. Determining equal comfort lines. J ASHVE 29:165–176.

International Organization for Standardization (ISO). 1985. ISO 7726. Thermal Environments—Instruments and Methods for Measuring Physical Quantities. Geneva: ISO.

—. 1989a. ISO 7243. Hot Environments—Estimation of the Heat Stress on Working Man, Based on the WBGT Index (Wet Bulb Globe Temperature). Geneva: ISO.

—. 1989b. ISO 7933. Hot Environments—Analytical Determination and Interpretation of Thermal Stress using Calculation of Required Sweat Rate. Geneva: ISO.

—. 1989c. ISO DIS 9886. Ergonomics—Evaluation of Thermal Strain by Physiological Measurements. Geneva: ISO.

—. 1990. ISO 8996. Ergonomics—Determination of Metabolic Heat Production. Geneva: ISO.

—. 1992. ISO 9886. Evaluation of Thermal Strain by Physiological Measurements. Geneva: ISO.

—. 1993. Assessment of the Influence of the Thermal Environment using Subjective Judgement Scales. Geneva: ISO.

—. 1993. ISO CD 12894. Ergonomics of the Thermal Environment—Medical Supervision of Individuals Exposed to Hot or Cold Environments. Geneva: ISO.

—. 1993. ISO TR 11079 Evaluation of Cold Environments—Determination of Required Clothing Insulation, IREQ. Geneva: ISO. (Technical Report)

—. 1994. ISO 9920. Ergonomics—Estimation of the Thermal Characteristics of a Clothing Ensemble. Geneva: ISO.

—. 1994. ISO 7730. Moderate Thermal Environments—Determination of the PMV and PPD Indices and Specification of the Conditions for Thermal Comfort. Geneva: ISO.

—. 1995. ISO DIS 11933. Ergonomics of the Thermal Environment. Principles and Application of International Standards. Geneva: ISO.

Kenneth, W, P Sathasivam, AL Vallerand and TB Graham. 1990. Influence of caffeine on metabolic responses of men at rest in 28 and 5C. J Appl Physiol 68(5):1889–1895.

Kenney, WL and SR Fowler. 1988. Methylcholine-activated eccrine sweat gland density and output as a function of age. J Appl Physiol 65:1082–1086.

Kerslake, DMcK. 1972. The Stress of Hot Environments. Cambridge: Cambridge University Press.

LeBlanc, J. 1975. Man in the Cold. Springfield, IL, US: Charles C Thomas Publ.

Leithead, CA and AR Lind. 1964. Heat Stress and Head Disorders. London: Cassell.

Lind, AR. 1957. A physiological criterion for setting thermal environmental limits for everybody’s work. J Appl Physiol 18:51–56.

Lotens, WA. 1989. The actual insulation of multilayer clothing. Scand J Work Environ Health 15 Suppl. 1:66–75.

—. 1993. Heat transfer from humans wearing clothing. Thesis, Technical University. Delft, Netherlands. (ISBN 90-6743-231-8).

Lotens, WA and G Havenith. 1991. Calculation of clothing insulation and vapour resistance. Ergonomics 34:233–254.

Maclean, D and D Emslie-Smith. 1977. Accidental Hypothermia. Oxford, London, Edinburgh, Melbourne: Blackwell Scientific Publication.

Macpherson, RK. 1960. Physiological responses to hot environments. Medical Research Council Special Report Series No. 298. London: HMSO.

Martineau, L and I Jacob. 1988. Muscle glycogen utilization during shivering thermogenesis in humans. J Appl Physiol 56:2046–2050.

Maughan, RJ. 1991. Fluid and electrolyte loss and replacement in exercise. J Sport Sci 9:117–142.

McArdle, B, W Dunham, HE Halling, WSS Ladell, JW Scalt, ML Thomson and JS Weiner. 1947. The prediction of the physiological effects of warm and hot environments. Medical Research Council Rep 47/391. London: RNP.

McCullough, EA, BW Jones and PEJ Huck. 1985. A comprehensive database for estimating clothing insulation. ASHRAE Trans 91:29–47.

McCullough, EA, BW Jones and T Tamura. 1989. A database for determining the evaporative resistance of clothing. ASHRAE Trans 95:316–328.

McIntyre, DA. 1980. Indoor Climate. London: Applied Science Publishers Ltd.

Mekjavic, IB, EW Banister and JB Morrison (eds.). 1988. Environmental Ergonomics. Philadelphia: Taylor & Francis.

Nielsen, B. 1984. Dehydration, rehydration and thermoregulation. In E Jokl and M Hebbelinck (eds.). Medicine and Sports Science. Basel: S. Karger.

—. 1994. Heat stress and acclimation. Ergonomics 37(1):49–58.

Nielsen, R, BW Olesen and P-O Fanger. 1985. Effect of physical activity and air velocity on the thermal insulation of clothing. Ergonomics 28:1617–1632.

National Institute for Occupational Safety and Health (NIOSH). 1972. Occupational exposure to hot environments. HSM 72-10269. Washington, DC: US Department of Health Education and Welfare.

—. 1986. Occupational exposure to hot environments. NIOSH Publication No. 86-113. Washington, DC: NIOSH.

Nishi, Y and AP Gagge. 1977. Effective temperature scale used for hypo- and hyperbaric environments. Aviation Space and Envir Med 48:97–107.

Olesen, BW. 1985. Heat stress. In Bruel and Kjaer Technical Review No. 2. Denmark: Bruel and Kjaer.

Olesen, BW, E Sliwinska, TL Madsen and P-O Fanger. 1982. Effect of body posture and activity on the thermal insulation of clothing: Measurements by a movable thermal manikin. ASHRAE Trans 88:791–805.

Pandolf, KB, BS Cadarette, MN Sawka, AJ Young, RP Francesconi and RR Gonzales. 1988. J Appl Physiol 65(1):65–71.

Parsons, KC. 1993. Human Thermal Environments. Hampshire, UK: Taylor & Francis.

Reed, HL, D Brice, KMM Shakir, KD Burman, MM D’Alesandro and JT O’Brian. 1990. Decreased free fraction of thyroid hormones after prolonged Antarctic residence. J Appl Physiol 69:1467–1472.

Rowell, LB. 1983. Cardiovascular aspects of human thermoregulation. Circ Res 52:367–379.

—. 1986. Human Circulation Regulation During Physical Stress. Oxford: OUP.

Sato, K and F Sato. 1983. Individual variations in structure and function of human eccrine sweat gland. Am J Physiol 245:R203–R208.

Savourey, G, AL Vallerand and J Bittel. 1992. General and local adaptation after a ski journey in a severe arctic environment. Eur J Appl Physiol 64:99–105.

Savourey, G, JP Caravel, B Barnavol and J Bittel. 1994. Thyroid hormone changes in a cold air environment after local cold acclimation. J Appl Physiol 76(5):1963–1967.

Savourey, G, B Barnavol, JP Caravel, C Feuerstein and J Bittel. 1996. Hypothermic general cold adaptation induced by local cold acclimation. Eur J Appl Physiol 73:237–244.

Vallerand, AL, I Jacob and MF Kavanagh. 1989. Mechanism of enhanced cold tolerance by an ephedrine/caffeine mixture in humans. J Appl Physiol 67:438–444.

van Dilla, MA, R Day and PA Siple. 1949. Special problems of the hands. In Physiology of Heat Regulation, edited by R Newburgh. Philadelphia: Saunders.

Vellar, OD. 1969. Nutrient Losses through Sweating. Oslo: Universitetsforlaget.

Vogt, JJ, V Candas, JP Libert and F Daull. 1981. Required sweat rate as an index of thermal strain in industry. In Bioengineering, Thermal Physiology and Comfort, edited by K Cena and JA Clark. Amsterdam: Elsevier. 99–110.

Wang, LCH, SFP Man and AN Bel Castro. 1987. Metabolic and hormonal responses in theophylline-increased cold resistance in males. J Appl Physiol 63:589–596.

World Health Organization (WHO). 1969. Health factors involved in working under conditions of heat stress. Technical Report 412. Geneva: WHO.

Wissler, EH. 1988. A review of human thermal models. In Environmental Ergonomics, edited by IB Mekjavic, EW Banister and JB Morrison. London: Taylor & Francis.

Woodcock, AH. 1962. Moisture transfer in textile systems. Part I. Textile Res J 32:628–633.

Yaglou, CP and D Minard. 1957. Control of heat casualties at military training centers. Am Med Assoc Arch Ind Health 16:302–316 and 405.