Mechanical properties of wood such as stiffness and strength vary locally especially due to heterogeneities and anisotropy. Analytical models and numerical simulations of wooden boards are able to represent varying material orientation e.g. with local fiber directions from laser scanning as input for the prediction of strength. Current Finite Element Models reconstructed the grain orientation by means of computationally demanding fluid analysis around obstacles like knots; whereas the available fiber pattern, captured by means of laser scanning, was passed solely into the detection of knots, but not directly processed for the inclusion of material fiber orientation. Therefore, the goal of this paper was the development of a numerical approach to directly include locally varying measured fiber orientation with orthotropic material properties and to predict the tensile strength of boards with reduced computational effort. Therefore, the stiffness was transformed element-wise according to the measured fiber deviations and the local fiber stress components were computed for the specific tensile load case. For the virtual strength prediction, numerical maximum stress values were compared to experimental tensile strength. Good agreements were observed with reduced computational effort compared to existing approaches between numerical and experimental results.@en

When curved glulam beams are subjected to bending, tension perpendicular to grain is introduced. In this work, Monte-Carlo analysis is conducted to study the influence of the annual-ring orientation of individual board on stress concentration in the radial direction. A virtual cutting program is developed which simulates the sawing processes of boards and yields a database of sawn boards, each with a characterized material direction with respect to the pith location. Lamellas are randomly selected from the database to construct curved glulam beams in FE models for Monte-Carlo simulation. Parameters including lumber diameter distribution, taper function, orthotropic elastic properties are considered for four common wood genera in Europe: spruce, pine, larch, and beech. Statistical analysis shows that when the annual-ring effect is considered, the average radial stress in the middle height area of the glulam beams can reach up 1.56 times (even higher on the middle width and height area) compared to the common practice, where wood is usually simplified as transverse isotropic material. Such stress concentration, which can trigger early damage of the structure, indicates an underestimation of stress perpendicular to grain in common practice. Moreover, parameter studies show that the stress redistribution depends not only on the mechanical properties of wood species but also on the width of the board layer.@en

This contribution aims for an enhanced numerical representation for strength prediction of timber. This implies a validated elasto-plastic continuum damage model which considers orthotropy and heterogeneity of the material, and represents the ductile behavior under compression and the brittle material behavior under tension dependent on the three-dimensional orthogonal fiber directions. The behavior under compression is captured by Hill (1948) plasticity and an exponential hardening law enhanced by the loading direction dependency. The same model covers also the brittle damaging behavior by means of continuum damage mechanics (CDM). In this study, a separated damage mode (SDM) criterion with simultaneously evolving damage variables is investigated. After the experimental validation of the model for axially loaded clear wood samples, the developed numerical model is implemented to a sawn timber with fiber deviation, where homogenization of the material and simplification to transverse isotropy is not anymore possible. The 3D orthotropic material behavior is experimentally validated for this application example with bi-axial loading and aims for further numerical investigation of wood with heterogeneities as occurring in sawn (hard)wood for its efficient use in engineered wood products such as glued laminated timber.@en

This contribution aims for an enhanced numerical representation for strength prediction of timber. This implies a validated elasto-plastic continuum damage model which considers orthotropy and heterogeneity of the material, and represents the ductile behavior under compression and the brittle material behavior under tension dependent on the three-dimensional orthogonal fiber directions. The behavior under compression is captured by Hill (1948) plasticity and an exponential hardening law enhanced by the loading direction dependency. The same model covers also the brittle damaging behavior by means of continuum damage mechanics (CDM). In this study, a separated damage mode (SDM) criterion with simultaneously evolving damage variables is investigated. After the experimental validation of the model for axially loaded clear wood samples, the developed numerical model is implemented to a sawn timber with fiber deviation, where homogenization of the material and simplification to transverse isotropy is not anymore possible. The 3D orthotropic material behavior is experimentally validated for this application example with bi-axial loading and aims for further numerical investigation of wood with heterogeneities as occurring in sawn (hard)wood for its efficient use in engineered wood products such as glued laminated timber.@en