Three coupled two-dimensional viscoelastic creep models for orthotropic material are analyzed. The models of different complexity are mathematically formulated and implemented in a finite element software. Required viscoelastic material parameters are determined by calibration procedure, where numerical results are compared against experimentally obtained viscoelastic strains caused by tensile or shear loading. Finally, a comparison method is used to evaluate the accuracy of strain predictions of each particular model. The analysis shows that all the models are able to accurately predict viscoelastic creep simultaneously in two perpendicular directions for various periods of time and wood species. Calculated numerical values of the viscoelastic material parameters suitable for the three models and wood species, i.e., Douglas fir (Pseudotsuga menziesii), Norway spruce (Picea abies), Japanese cypress (Chamaecyparis obtusa), and European beech (Fagus sylvatica L.), under constant tensile loading are also given.

Changes in relative humidity of the ambient air, RH (%), cause wetting and drying of wood material, which results in non-uniform moisture contents or moisture gradients, and consequently in moisture-induced stresses and strains in the glued-laminated timber (glulam) members. The aim of the present paper is to perform a hygro-mechanical analysis to predict the mechanical behavior of glulam specimens exposed to two RH regimes, causing wetting from 50% to 90% RH and drying from 90% to 50% RH, and compare the numerical to the experimental results. The aims are also to quantitatively analyze the influence of characteristic material parameters required in the multi-Fickian moisture transport model and the mechanical model on moisture-induced strains and stresses in glulam specimens and to determine the possibility of cracking of the material by analyzing the maximum tensile stresses perpendicular to the grain. Accurate numerical predictions of moisture contents and moisture-induced strains are obtained in the glulam specimens during wetting and drying as compared to the experimental results. The influence of a particular characteristic material parameter on moisture-induced strains and stresses is characterized as significant, but not crucial when a rough numerical estimation of the mechanical behavior of the glulam beam exposed to RH changes is required.

A hygro-mechanical (H-M) analysis of a wooden specimen sustaining a mechanical load while subjected to varying relative humidity was performed to predict the long-term rheological behavior of wood. The numerical analysis was based on the experimental results of total strains, monitored in two orthotropic material directions on oak wood specimens under constant uniaxial compression and with moisture content (MC) variation. For the moisture analysis, a multi-Fickian moisture transport model (MFMTM) was used to obtain temporal and spatial MC fields, which were the input data in the mechanical analysis. The presented mechanical model assumed a decomposition of the total strains into the elastic, viscoelastic and mechanosorptive strains and the strains due to shrinkage and swelling. The moisture and mechanical analyses required material parameters, which were taken from the literature or were empirically obtained by a fitting procedure. The performed H-M analysis gave accurate numerical predictions of the experimentally obtained total strains in two orthotropic directions simultaneously. Thus, the analysis developed has a high potential for predicting the long-term rheological behavior of timber structures, assuming that the material parameters are determined previously, based on specific, extensive, multidimensional experimental analyses.