Hygromorphic response dynamics of oak: towards accelerated material characterization

When a wood board is exposed to a change in relative humidity on only one of its surfaces, e.g. in case of flooring or a panel painting, the resulting asymmetric moisture content profile induces differential expansion over the thickness. Consequently a bending moment causes the board to curve. A theory is presented to describe the bending of a wood board due to a step change in relative humidity. The board is assumed to be homogeneous, isotropic, and linearly elastic. Moisture transport is presumed to obey the diffusion equation with constant coefficients, such that moisture transport can be directly related to the bending of the board. It is shown that the transient deflective behavior provides the diffusion coefficient and the final length change yields the linear hygroscopic expansion coefficient. Derived diffusion coefficients are in good agreement with values in literature. Furthermore, a scaling law for the deflection of the board is proposed, which is seen to be followed qualitatively but not quantitatively by experiments. Finally, by assuming the deflection of the board to be the response of a linear system, the deflective frequency response of the board can be predicted from its step response. The results allow upscaling of deflection and expansion, such that behavior of thick boards can be determined from an experiment using a thin board.


Introduction
The shape change due to moisture sorption is a familiar phenomenon in natural systems. A wellknown example of such a hygromorphic effect is the opening and closing of pine cones at different environmental conditions [1]. Due to the bilayered structure of the individual scales, a change in the relative humidity of the ambient air causes differential expansion of the two layers. As these two layers are firmly attached to each other this will result in a bending moment, thereby opening the cone. Another well-known example of a hygromorphic material is wood. Upon a change in the relative humidity of the ambient air, wood will exchange moisture with the ambient air [2], resulting in moisture transport in the material in order to reach an equilibrium. The resulting change in moisture content of the wood leads to expansion or shrinkage, as the fibers which constitute the wood are pushed apart by the entering moisture or approach when moisture is released. This property of wood has been used e.g. in mechanical hygrometers [3], or in barrels ensuring the sealing of the barrel. The anisotropic nature of wood, with differences in expansion in its three principal directions, can also induce out-of-plane deformation. The warping of wood boards has been applied in the construction of 19th century boathouses in Nordmøre, Norway [4]. The walls of these boathouses are cladded with pine boards, nailed towards their upper edges. The anisotropy in expansion causes the boards of the boathouse to bend upwards during dry conditions, allowing natural ventilation through the façade. During humid conditions, the boards bend downwards, ensuring weathertight sealing. More recently, a prototype façade has been developed which reacts to environmental changes by the warping of individual thin wood slices constituting the façade [5]. The effect may be amplified by adding a passive layer with different expansion to the active wood veneer, e.g. a polymer [6], thereby mimicking a pine cone or a bimetallic strip. Bending in wood can also be caused by differential expansion resulting from an asymmetric moisture content distribution. This phenomenon has been studied during kiln-drying [7][8][9][10][11], and in the case of a board simulating a painted panel [12].
The expansion and shrinkage of wood are, however, more commonly associated with stresses in the material [13]. When these stresses exceed the strength of the material, cracks can be formed. These cracks are aesthetically undesired when a wooden art object is concerned [14], or in the case of a wooden parquet floor [15]. Cracks can also have consequences in a structural component, which can ultimately give rise to structural failure. A wide variety of experimental methods has been used to study the hygromorphic behavior of wood, e.g. using displacement measurement gauges [16], a linear variable differential transformer mounted to a micrometer [17], phase contrast synchrotron X-ray tomography [18], digital image correlation [19,20], a combination of nuclear magnetic resonance (NMR) and a fiber Bragg grating optical sensor [21], or a combination of NMR and X-ray microtomography [22]. All these methods are either complex or require sophisticated post-processing to derive useful, generalizable knowledge about the relation between moisture transport and expansion. Alternatively one could relate moisture transport to hygromorphic macroscopic behavior by analyzing the transient morphological or mechanical behavior during non-equilibrium conditions. Scherer and co-workers [23][24][25][26] have derived transport properties of different porous media from a bending experiment. The method relies on the redistribution of liquid in the pores of a liquid-saturated porous beam when subjected to three-point bending. Its mechanical response at constant deflection is influenced by the liquid flow and thus provides information about the permeability of the material. Moisture-induced bending of a beam of clay-bearing stone exposed to water on one side was thoroughly analyzed by Scherer and Jimenez-Gonzalez [27]. Similarly, Holmes et al. [28] considered the bending of gel beams when wetted unilaterally by hexane.
The basic idea in this study, shown in Fig. 1, is analogous. Here we will consider a wood board with a thickness much smaller than its length and width, clamped at one end. The board has an impermeable layer on five surfaces, leaving one surface exposed to the atmosphere surrounding it. A semi one-dimensional moisture transport experiment is thus created. Here we will assume that the impermeable layer has no influence on the mechanical behavior of the board. The board itself is initially in equilibrium with the surrounding air, and therefore has a flat moisture content profile (Fig. 1a) when suddenly the relative humidity of the surrounding air is changed. The moisture content at the exposed surface consequently rises, and moisture is transported into the material resulting in a moisture content gradient. The consequent asymmetric moisture profile causes unequal expansion over the thickness; the expansion near the exposed surface is the highest and decreases towards the coated surface. These differences in expansion result in a constant moment, and as a consequence the board bends (Fig. 1b). Due to moisture transport the moisture content profile continuously changes and consequently the bending changes continuously over time. The deflection has a maximum in time when the moment induced by the differences in expansion is the largest. Progressing in time, the moisture starts to equilibrate throughout the sample and the moisture content becomes homogeneous again. In the final situation the board is straight again, but elongated due to the increase in moisture content (Fig. 1c). Since the process is governed by moisture transport, one can expect that the dynamics of the bending is dependent on both object geometry, in this case thickness, and material properties such as diffusivity and moisture expansion.
The aim of this study is the analysis of the dynamic moisture-induced bending of a wooden board. A theoretical model is introduced by coupling moisture transport to a linear elastic model for the moistureinduced bending. We analyze the bending in time after a sudden increase in relative humidity and couple it to internal moisture transport. We will first discuss the theory. Experimental verification is provided with boards with different thicknesses. The postulated step change in relative humidity is, however, not very realistic since changes of environmental conditions are continuous, generally manifested in daily or seasonally repeating cycles [29]. Nevertheless, the hygromorphic step response of the board can provide information on the frequency response characteristics. The behavior at very low frequencies, which is timeconsuming to determine experimentally, can thus be retrieved from a short experiment. The experimental study of moisture transport in wood during cycles is often limited to block waves; examples investigating sinusoidal fluctuations are scarce [16,17,30]. We will discuss the theory to derive the frequency response, and the experimental verification with sinusoidal relative humidity fluctuations.

Hygro-bending theory
In this section, the bending of a clamped hygromorphic board as shown schematically in Fig. 2, sealed on five of its six sides, is described by introducing a socalled hygro-bending theory. First the moisture transport arising from a sudden change of the ambient relative humidity will be described analytically. A simple, linear elastic model is then proposed to express the consequent bending of the board. x Fig. 1 The conceptual idea of the bending experiment used in this study. A board of hygroscopic material is exposed to a step change in the ambient relative humidity (RH) at one side. The board is sealed on all sides except the one exposed to the environment resulting in a one-dimensional moisture transport. a Initially the sample is in equilibrium and the moisture content (MC) is constant throughout the sample. b A sudden increase in the ambient relative humidity causes a moisture content gradient in the sample resulting in a bending moment due to differential expansion. c In the final situation the moisture has equilibrated throughout the sample and the board is straight again, but elongated due to a higher moisture content dimensional moisture transport. Here we will assume that moisture transport is characterized by a constant diffusion coefficient D throughout the thickness d of the board. The other surface is sealed and therefore has a no-flux boundary condition. Mathematically this can be formulated as: where c is the moisture content in kg water/kg dry wood. An analytical solution can be obtained by separation of variables, see e.g. [31]. The result can be formulated in the form of the Fourier series A second solution can be formulated when the medium can be considered semi-infinite, i.e. when the penetration depth of the moisture l is small compared to the thickness of the board (l/d ( 1). In this regime, the moisture concentration can be described as with the penetration depth of the moisture 2.2 Transient bending

Coupling with moisture distribution
A change in moisture content in a hygromorphic material results in expansion or shrinkage. In case of a uniform moisture distribution c and no lateral contraction, the expansion of the board e is given by where L is the length of the board at time t, L 0 the length of the board at t 0 , a the linear hygroscopic expansion coefficient [3,29], and c 0 the uniform moisture content at time t 0 . After a step change in the moisture content at the open surface, the moisture content will not be uniform. As a result the timeevolution of the expansion of the board can no longer be described by Eq. (5). Furthermore, the board bends Fig. 2 Schematic representation of a board as used in the theory of the bending. The board has an initial length L 0 and thickness d, and is exposed on one surface to a change in moisture content from c 0 to c d due to the uneven expansion across its thickness. In order to describe the bending of the board, the assumption is made that the material is linearly elastic with no lateral contraction (for a more in-depth discussion see the strictly analogous problem of thermal expansion in chapter 10 of Boley and Weiner [32]). By doing so, the expansion of the board at halfthickness e can be expressed as Note that for a constant moisture content, Eq. (6) simplifies to Eq. (5). The bending of the board can be characterized by the angle h of its cross-section (see Fig. 2 Here the integral represents the moment exerted on the board as a consequence of an unsymmetrical moisture content distribution. It can be shown that for the general case of a symmetric moisture content profile, i.e. c(x, t) = c(d-x, t), the integral in Eq. (7) vanishes and the angle of the cross-section of the board equals zero. Hence in case of a two-sided exposure to a change in relative humidity there will be no deflection, only an expansion in the y-direction. In case of a onesided step change, the moisture profile will approach the constant value c d for t ? ?, and consequently the angle h approaches zero again.
The bending of the board causes deflection, which is maximal at the free end of the board. The xdeflection at the free end at half-thickness is given by and the corresponding y-deflection can be expressed as where u = tan -1 h. Under the assumption that the deflection of the beam remains small, i.e. u & h, a compact expression for the deflection can be derived by taking into account only the first two terms of the Taylor series expansion of Eq. (8) and assuming

Explicit expressions for expansion and deflection
We have now obtained analytical expressions to describe the expansion of the board at half-thickness and the deflection of its free end following a step change in the moisture content at one surface. We also introduced two expressions for the moisture content evolution. Three expressions for the expansion and the deflection can be constructed with the two descriptions of the moisture content evolution. First of all, a total solution for the expansion e tot and the deflection w tot is obtained by substitution of Eq. (2) into the expressions for e and w in Eqs. (6) and (10), which results in and w tot t ð Þ ¼ 24 Note that the expansion and deflection in Eqs. (11) and (12) can be rewritten as where e* and w* are the scaled expansion and deflection respectively, as a function of the dimensionless time t* = Dt/d 2 . The scaled expansion and deflection for the total solution can accordingly be expressed as Note that, as in Eq. (2), the expressions contain an infinite sum of modes, which makes a concise analysis impossible. A second set of solutions for the expansion e M and deflection w M can be obtained by taking into account only the first two modes: and This results in a 'long-term' solution, since higher order modes are neglected and the terms with the smallest decay constants remain. At small times, the long-term solution deviates largely from the total solution. At large times, however, it will converge to the total solution, due to the dominance of the low modes. Similar to total solution, the long-term solution can also be rewritten according to Eq. (13), which yields for the scaled expansion and deflection: A third set of expressions for the expansion e S and deflection w S can be obtained by substitution of the short-term solution, i.e. Eq. (3), in the expressions for e and w to yield and The scaled expansion and deflection for the shortterm solution are then obtained as: The scaled expansions, i.e. Eqs. (14), (18), and (22), are shown in Fig. 3a as a function of the square root of t * . As can be seen, the short-term solution e S * is linear with a slope 2/Hp. This linear behavior is initially also seen in the total solution e tot * , up to a value for Ht * of *0.5. The long-term solution deviates at low values of t * , but converges to the total solution from a value for Ht * of *0.2. The scaled deflections from Eqs. (15), (19), and (23) are plotted in Fig. 3b as a function of t * . All three solutions rise up to a point of maximum deflection, after which the deflection decreases again. This is also apparent by comparing the long-term and the short-term solution for the scaled deflection. We observe that both solutions are composed of a rising component which is dominant at short times, and a negative component, which gains dominance at large times. In the long-term solution, these components are exponentials with different decay constants; in the short-term solution, these components are a square root and a linear term respectively. Nevertheless, the qualitative behavior is similar. Furthermore it is seen that the short-term solution follows the complete solution perfectly up to the point of maximum deflection, and starts to deviate considerably hereafter. The opposite is observed for the long-term solution, which deviates significantly for short times, but converges to the complete solution even before the point of maximum deflection. A hybrid solution can therefore be proposed, which consists of the short-term solution for times up to the maximum deflection, and the long-term solution for the deflection after the maximum.

Derivation of parameters from transient data
To find the time t m at which the maximum deflection w m occurs we can take the zero of the derivative with respect to time for the long-term and short-term solution, which, for both solutions, has the form with Hence the time at which the board is deflected maximally is proportional to the typical diffusion time d 2 /D. Inserting t m into Eqs. (17) or (21) yields the maximum deflection of the board with The maximum deflection is thus proportional to the linear expansion coefficient, the length of the board squared, and the step in moisture content at the surface, and inversely proportional to the thickness of the board. Furthermore, the point of maximum deflection is similar for the long-term and short-term solution, with / M // S = 1.05 and w M /w S = 1.01. Moreover, it appears that the maximum deflection is independent of the diffusion coefficient. The diffusion coefficient only plays a role in the time at which the maximum occurs. By performing an experiment with a board of thickness d, the diffusion coefficient can be derived from the results by rewriting Eq. (24) to arrive at To determine the time of maximum deflection t m , the x-deflection is first fitted with the double exponential function where b 0 , c 1 , and c 2 are fitting parameters. Note that Eq. (29) has a similar form as Eq. (17). Equation (29) has a maximum at the time Substitution of Eq. (30) in Eq. (29) yields the maximum deflection. Moreover, the diffusion coefficient can be calculated by substitution of t m in Eq. (28). Analogous to the x-deflection, the y-deflection can be fitted by: where b 3 , c 3 , b 4 and c 4 are fitting parameters. For t ? ? the x-deflection approaches zero, and consequently the y-deflection approaches the length change of the board. Hence by analyzing the bending response of a board after a step change in the ambient relative humidity the diffusivity can be obtained. The linear hygroscopic expansion coefficient a can be derived from Eq. (16) as the asymptote for t ? ?, in case the step in moisture content is known. Furthermore, Eq. (26) provides us a scaling law with which experimental results using boards of different thicknesses can be compared.

From step response to frequency response
To derive frequency response characteristics from a step response, the deflective behavior of the wood  Fig. 3 a The scaled expansion e* and b the scaled deflection w* of the three different solutions (total, long-term, short-term, subscripts tot, M, and S respectively) as a function of the square root of dimensionless time t* board is considered as the dynamic response of a linear system [33]. The frequency response for the xdeflection can be obtained by first determining the Laplace transform of Eq. (29) in terms of the complex variable s: Equation (29) is the response of the linear system to a step change in the moisture content at the surface at t = 0. Since the board is approximated by a linear system, only the difference in initial and final value of the moisture content is of interest, not the actual values. The step change u(t) can therefore be written as where a is the step in moisture content. The Laplace transform of Eq. (33) is: In a linear system, the Laplace transforms X and U of the output and input function respectively, are related by: where H(s) is the transfer function of the system. Note that this relation holds for all input functions u(t) for which a Laplace transform exists. As a consequence, if the transfer function of a system is known, the inverse Laplace transform of Eq. (35) will yield the transient response x(t) of the system to an arbitrary input u(t).
Using Eqs. (32) and (34) in Eq. (35), the transfer function of the x-deflection can be expressed as If the transfer function of the system is known, its frequency response can be determined by considering the real and imaginary parts of H(ix), with x the angular frequency. H(ix) is a vector in the complex plane; its magnitude is the amplitude of the steady-state output signal to a sine input signal with unit amplitude, and its angle with the positive real axis determines the phase shift of the output signal. Using these results, Bode plots can be constructed which show the amplitude on a logarithmic scale and phase on a linear scale as a function of frequency. A similar analysis can be performed to determine the frequency dependence of the y-deflection of the board. The transfer function of the y-deflection can be obtained using Eq. (34) and the Laplace transformation of Eq. (31) as In a similar fashion, the frequency behavior of the y-deflection can be calculated by considering the real and imaginary part of G(ix).

Materials and methods
Experiments were conducted on oak boards where the bending due to one-sided exposure to a change in relative humidity is assessed. First, the material of interest and the preparation of samples are described, after which the experimental procedure is explained with the description of the set-up.

Material
Oak boards with a length of 100 mm, a width of 30 mm, and different thicknesses (2, 4, 6, 8 mm) were sawn from a large oak wood board. Because wood is an anisotropic material with properties dependent on the characteristic directions, attention is paid to the grain direction in the samples. All samples were prepared such that the grain direction is along the width of the board. Perpendicular to that, the tangential direction is in the length of the board, and the radial direction along the thickness. Transport hence occurs in the radial direction of the wood, whereas the expansion along the length of the board is the expansion in tangential direction. Since all experiments start at a relative humidity of 50%, the boards were placed in a desiccator containing a saturated Mg(NO 3 ) 2 Á6H 2 O solution, ensuring a relative humidity of 53%. Prior to the experiment, Bison silicone kitÓ was applied on five sides of the board, leaving one of the two largest surfaces open. The silicone kit was allowed to dry for at least 1 day.

Experimental set-up
The experimental set-up designed to measure the deflection of a wood board due to humidity fluctuations is shown in Fig. 4. One end of the wood board is clamped between two PVC strips. The other end of the sample is also clamped between two strips, on which a pointer is mounted. This setup is placed in a plastic container in which the relative humidity can be controlled by the help of a humidifier, which operates in the relative humidity range 0-95%. A Dino-LiteÓ digital microscope records images to measure the deflection. Using a Matlab optical recognition program the position of the pointer is recorded automatically.
In the humidifier, compressed air is first dried by a membrane dryer. The dehumidified stream is then divided over two BronkhorstÓ flow controllers, with which the flow can be set precisely. One of the air streams is conducted through tubing with semipermeable sides, embedded in a temperature controlled bath, to moisten the air. After leaving the bath, the air stream is conducted through a condenser to remove excess moisture. Both streams are mixed to produce a uniform flow with a relative humidity, the mixing proportions then determine the relative humidity of the air: where V wet and V dry are the wet and dry flow respectively [m 3 /s].

Step response
We have first measured the step response of boards with various thicknesses. In Fig. 5a the deflection in both directions and the calculated length change at half-thickness of the free end of a 4 mm thick board is shown as a function of time after a step change in the relative humidity from 50 to 90%. The characteristic shape as predicted from the theory is seen. After the step change in relative humidity the board immediately starts to bend up to a maximum deflection, subsequently it more slowly bends back towards its initial x-position. The deflective response in the y-direction is dependent on the thickness of the board, since it is caused by two processes. Purely geometrically, the bending of an initially straight board causes a change in the y-position of its end. In case of a hygro-expansive material, the length of the board also changes during moisture sorption. This double effect causes the y-position of the end of the board to first decrease before increasing. In a thicker board, which deflects less in the x-direction, expansion dominates and no initial dip in the y-position is recorded.
The deflections in x-and y-direction are fitted with Eqs. (29) and (31) respectively, as shown in Fig. 5b. The resulting time of maximum deflection t m for the various boards is shown in Fig. 6a as a function of the thickness d of the different boards. A fit with Eq. (24) is added to illustrate the agreement between theory and experiments. The fitted diffusion coefficient has a value of 4.5 9 10 -11 m 2 /s, which is in accordance with values reported in literature for oak [34][35][36]. The derived diffusion coefficients for the individual boards are shown in Fig. 6b. The resulting range of 2-6 9 10 -11 m 2 /s is relatively narrow for a highly inhomogeneous material such as oak.   -(b 3 ? b 4 ). Since the step in ambient relative humidity is known instead of the step in moisture content, an additional operation needs to be performed to obtain the expansion coefficient a. The step in moisture content is measured on a sample of the same board using dynamic vapor sorption (DVS). A step in relative humidity from 50% to 90% results in c d -c 0 = 0.11. This result is used to retrieve the expansion coefficient. The expansion coefficients derived from the asymptote of the expansion after the step change are shown in Fig. 7. The values are between 0.2 and 0.5, with an average of the derived expansion coefficients (0.35) being in accordance [35] or slightly higher than values reported in literature for the tangential expansion of oak (*0.3) [29]. The scatter in values is most likely caused by differences in initial and final moisture content of the different boards. This can be due to structural inhomogeneities or sorption hysteresis.
With the diffusion and expansion coefficients of the boards determined, the maximum deflection can be scaled according to Eq. (26). The result is shown in Fig. 8, where the theoretical curve 1/d should ideally overlay the experimental results. However, a theoretical overestimation of the maximum deflection is seen. Scaling the theoretical line by a factor results in a more accurate description of the experimental results. The qualitative scaling behavior of 1/d is thus followed by experiments with different thicknesses, but the absolute values are overestimated, which may have several causes. In the present theory, it is x-deflection y-deflection Fit (a) (b) Fig. 5 a The x-and y-deflections w and v of the free end and the length change DL of a clamped, 4 mm thick board at halfthickness due to a step change in the relative humidity of the ambient air from 50 to 90%, and b the same experiment fitted with Eqs. (29) and ( assumed that the surface of the board reaches equilibrium with the surrounding air immediately. A stationary boundary layer at the surface can, however, increase resistance to moisture sorption. If the transport of moisture is limited by the resistive boundary layer, less asymmetry in the moisture profile is expected over time. As a consequence, the deflective response of the board is less intense. In the experimental set-up used in this study, humid air is blown into the plastic container with a flow such that the air in the container is refreshed at least once a minute. It is therefore presumed that no stagnant boundary layer is formed which attenuates moisture sorption. A more probable cause of the overestimated maximum deflection values is the assumption of wood as a homogeneous, isotropic, and linearly elastic material. Local structural fluctuations can e.g. lower the moisture content asymmetry over time, decreasing the deflection. Also expansion is presumed to follow the change in moisture content immediately, whereas in literature viscoelastic effects are assigned to wood [37]. These time-dependent effects can attenuate the bending of the board. Additional experiments should therefore be performed with a hygro-expansive model-material, which is preferably highly homogeneous and isotropic.
The scaled expansion e* as a function of the square root of dimensionless time t* is shown in Fig. 9a for different boards. The theoretically predicted time evolution of e tot * from Fig. 3 is plotted in the same graph. As can be seen, the expansion evolution of the experiments is in reasonable agreement with the predicted course. The experimental slope, however, is higher than the predicted slope. The dimensionless x-deflections w* of the different boards are shown in Fig. 9b. As already observed in Fig. 8, the predicted deflection is approximately twice as large as the scaled experimental results. Similar results were obtained by Holmes et al. [28], which was explained with the use of equilibrium data for scaling kinetic results. Nevertheless, the scaled experiments are all on one main master curve. The general behavior is similar for all experiments, with an initial linear deflection on the square root scale. After the maximum is reached, the boards bend back towards their initial x-position, i.e. towards zero deflection. To compare the experiments to the theory, the theoretical curve for w tot * scaled by a factor is shown too. It is seen that most boards bend back faster towards their initial x-position after the maximum has been reached. The opposite was found in experiments with clay-bearing stone conducted by Scherer and Jimenez Gonzalez [27]; a delay in the post-peak flattening resulted from the dry side being compressed by bending stresses, making it more difficult for water to cause swelling. The stresses which develop due to the bending can have an influence on the moisture transport due to their impact on the chemical potential of water and thus the moisture absorption capacity of the matrix [38]. Since a gradient in chemical potential is the driving force for moisture transport, stresses arising from swelling and shrinkage can have an influence on the moisture distribution over time [39], which again influence the bending.
In the analysis introduced by Scherer and Jimenez Gonzalez [27], the moisture-dependence of the stiffness was taken into account. An estimate of this effect on the bending of the wooden boards is shown in Fig. 10a, in which the E-modulus is varied as E = 10 b-kc , with b = 10 and k = 1, 2, 5, 10, 20, and 50 (see e.g. Scherer [40] for the influence on the bending). A step in the moisture content c from 0 to 0.2 is considered. As can be seen, the maximum deflection is approximately halved in case of k = 50. This, however, means that the E-modulus decreases with ten orders of magnitude over the considered moisture content range, which is far from realistic. Saft and Kaliske [35] proposed an E-modulus in the tangential direction which decreases *50% over the moisture content range considered here. This would correspond  Fig. 10, which has a negligible influence on the deflection of the board over time.
We have so far assumed the diffusion coefficient to be constant. In literature, moisture-dependent values have been measured. The effect on the scaled deflection is shown in Fig. 10b, where a diffusion coefficient of the form D = D 0 e kc is assumed. As can be seen, the moisture-dependent diffusion coefficient affects the maximum deflection. A positive k, i.e. an increasing diffusion coefficient with increasing moisture content, results in a higher scaled maximum deflection. Conversely, the maximum deflection is attenuated for negative k. Also the post-peak flattening is influenced by the moisture-dependent diffusion coefficient. A positive k advances the post-peak flattening, a negative k delays the flattening. A similarly pronounced influence is not found in the experiments presented in Fig. 9b. Furthermore, literature values for the diffusion coefficient generally increase with moisture content [3,34,41], resulting in a higher peak. The moisture-dependent diffusion coefficient alone can thus not explain the mismatch between predicted and experimental scaled deflection.
The discrepancy between predicted and experimental scaled deflection may be an indication that transport processes take place other than described by the diffusion equation, as proposed e.g. by Cunningham [42] and Krabbenhoft and Damkilde [43]. Cunningham [42] assumes vapor diffusion takes place in the cell cavities and moisture sorption occurs from the cavities to the cell wall, mathematically described as moisture content to a driving potential for moisture sorption. It is assumed here that the ratio l/m is the slope of the sorption curve relating cell wall moisture content to vapor content, and that the water vapor and cell wall moisture are initially and finally in equilibrium. This yields u 0 = l/m v 0 and u ? = l/m v ? , with the subscripts 0 and ? indicating initial and final values respectively. Equations (39) and (40) can then be scaled by introducing dimensionless variables: to arrive at The effect of sorption from cavities to cell walls on the bending of a board is estimated using Eqs. (42) and (43), complemented with a step change in the vapor content on one surface and a no-flux boundary condition on the other surface. A value for m * of 3.3 is used, which is the reciprocal of a linear sorption curve with a slope of 0.3. Equations are discretized using a backward Euler scheme and implemented in Matlab. Numerical simulations with different values of h * are performed.
Scaled deflections as a function of the square root of t * are shown in Fig. 11a for different values of h * . It can be seen that for slower sorption the maximum deflection is delayed and attenuated compared to the case of immediate equilibrium (h * = ?). The delay and attenuation of the maximum scaled deflection are a function of h * , as shown in Fig. 11b. Since it can be assumed that l, h, and D are no functions of the thickness, h * scales with d 2 . The scaled maximum deflection is then expected to be dependent on board thickness. Experimental results in Fig. 8, however, show that the scaled deflection follows Eq. (26), so that w m * is no function of d. Furthermore, the scaled time of maximum deflection would be influenced by board thickness. This would be reflected in a thickness-dependent diffusion coefficient, which is not observed in Fig. 6. A stationary boundary layer at the macroscopic surface attenuating the increase in moisture content at the surface has a similar effect, i.e. a thickness-dependent scaled maximum deflection. It is therefore not likely that moisture sorption to the cell wall or a stationary boundary layer can explain the observed discrepancy between the simple theory and experimental results. Experiments in which the moisture content distribution and bending are measured simultaneously, e.g. using nuclear magnetic resonance, are expected to provide more insight.

Frequency response
To arrive at a rapid means of characterizing its frequency behavior, the deflective step response of a 2 mm thick board was used according to the procedure introduced in Sect. 2.4. The resulting Bode amplitude plots of the deflection of the board in x-and y-direction are shown as a function of the frequency in Fig. 12. Three regimes can be distinguished. At low frequencies, the board is seen to deflect to a smaller extent in the x-direction, but its amplitude in the y-direction reaches a constant value. During these slow fluctuations, the asymmetry in the moisture content distribution is minor, resulting in small x-deflections. The moisture penetration is high, resulting in a high amplitude in the y-deflection. In an intermediate frequency interval, the board reacts heavily to the imposed fluctuations in relative humidity. For increasing frequencies, the amplitude in both x-and y-deflection is seen to decline; the variations in relative humidity are too fast to follow, i.e. the moisture does not penetrate far enough into the material to cause deflections.
In order to verify the frequency response predicted from its step response, we have conducted experiments with the same board where the relative humidity was varied sinusoidally at different frequencies. The equilibrium relative humidity was 50% and the amplitude 40% in all experiments. The result of one of these measurements for a board with a thickness of 2 mm is given in Fig. 13, and shown in the animation (Online Resource 1). The board is seen to deflect cyclically to the imposed relative humidity variations. Furthermore it is seen that the deflection during moisture absorption is slightly larger than during desorption. This can be explained with the sorption curve of wood, which has an S-shape with a relatively large increase at high relative humidity [3]. Although the input, i.e. the relative humidity, has a pure sinusoidal shape, the non-linearity of the sorption curve causes the moisture content at the exposed surface to follow a transformed course over time. In Fig. 12 we have plotted the amplitudes from individual experiments at different frequencies (periods of 134 h to 2 min). As can be seen these results overlap well with the frequency response as calculated from a single step experiment. Some quantitative deviations between predicted and experimentally found values are encountered, especially with the y-deflection. As mentioned earlier, the y-deflection has a dual source: it arises as a geometrical consequence of the bending and is caused by expansion/shrinkage. This explains the peak in the predicted curve, which is not seen to be confirmed by experiments. Nevertheless, the agreement between prediction based on the assumption of linear system behavior and actual experiments is striking. Despite non-linear effects affecting moisture transport and expansion in wood which are often found in literature, the behavior in these experiments is seen to be fairly linear. This result enables simple and insightful prediction of hygromorphic behavior of wood. Instead of measuring the response at various frequencies, which is time-consuming, the frequency behavior can be predicted with a short experiment.

Conclusions
A simple theory is proposed based on the diffusion equation and linear elasticity to analyze experiments x-deflection from step response y-deflection from step response x-deflection from direct measurement y-deflection from direct measurement  with oak boards exposed unilaterally to moisture. The diffusion coefficient and linear hygroscopic expansion coefficient of wood can be derived accordingly, which are in agreement with literature values. The scaling law concerning the maximum deflection of the board during a step change in the relative humidity is seen to describe the experiments well qualitatively, but overestimates the deflection quantitatively. Possible causes are the oversimplification of the involved processes, neglecting non-linear, non-elastic effects, local inhomogeneities and the influence of stresses on moisture transport. From a simple parameter study we can conclude that moisture-dependent values for the E-modulus and the diffusion coefficient cannot explain the mismatch in scaled deflection between the simple theory and experimental results. A different moisture transport model, with vapor diffusion and sorption to the cell wall, cannot explain the observed discrepancy either due to a thickness-dependence of the scaled deflection, which is not seen in experiments.
Approximation of the bending behavior of the wood board as linear system behavior allows the derivation of the frequency response from its step response. The behavior is well captured, reflected in the agreement in shape of the Bode plots predicted from a step response and direct measurement with sinusoidal relative humidity fluctuations. The ability to predict the frequency behavior of the board from a relatively short experiment is an encouraging result, especially in combination with the scalable responses for different boards.
Step responses of thick boards and frequency behavior at slow fluctuations, both highly time-consuming experiments, can thus be predicted with a short experiment with a thin board. This allows a fast characterization of the hygromorphic response of a material.
A promising extension of the experimental part of the study should be sought in the simultaneous measurement of the deflection and the local moisture content. The deflection can be assessed optically, as shown in this study, whereas the moisture content could be measured non-destructively by nuclear magnetic resonance. In this way the moisture content profile can be related directly to the hygromorphic response at any time. This combination provides more information about the moisture transport during the bending experiment, and possibly elucidates the interrelation between the stress distribution and moisture transport.