Contribution702 A Final
Contribution702 A Final
1 INTRODUCTION 123
Over the last few years, CLT structures have been
increasingly spreading in the construction sector, mainly
in Europe and North America. CLT panels are made up of
several layers of boards stacked crosswise and glued
together on their faces. Among the main advantages of
CLT, we can highlight its fast and efficient on-site
installation, its favourable seismic performance and its
excellent strength. As CLT panels are light-weight
structural elements with high stiffness and strength to
bending, compression and shear, they are an economically
competitive building system when compared to traditional
options and therefore, are a suitable candidate for some
applications in buildings which currently use concrete,
masonry and steel [1].
One important issue in the design of CLT structures which
still requires further investigation is the progressive
damage that may occur across the thickness of CLT
layers, or between adjacent layers [2]. In particular, the
design of CLT floor systems with low span-to-depth ratios
is often governed by the rolling shear capacity of CLT
plates, and therefore, its full understanding is of
paramount importance to provide the required structural
strength and prevent damage in CLT structures. This type
of damage can also be potentially present in CLT walls
subject simultaneously to eccentric compressive loads and
out-of-plane bending or shear loads.
Rolling shear failure consists of inter-fibre cracking due
to shear strains in the plane perpendicular to the
longitudinal axis of wood fibres. Figure 1 shows a typical
2 THEORETICAL FRAMEWORK
Over the last two decades, homogenisation-based multiscale computational modelling techniques have been
developed to predict accurately the mechanical response
of
heterogeneous
materials
with
hierarchical
microstructures [4]. Successful use of this methodology in
the context of wood materials and structures can be found
in [5-9]. In this constitutive theory, each material scale is
associated with a microstructure whose most statistically
relevant features are incorporated within a representative
volume element (RVE). This RVE is assumed to have a
(microscopic) characteristic length much smaller than the
macro-continuum, and at the same time, a size large
enough to capture the microscopic heterogeneities in an
averaged sense.
In this theory it is also assumed that the macroscopic or
homogenised strain tensor component ij at any arbitrary
point of the macroscopic continuum is the volume average
of the microscopic strain tensor component ij over the
domain of the RVE.
Similarly, the macroscopic or homogenised stress tensor
component ij is assumed to be the volume average of the
microscopic stress tensor component ij over .
In addition, we can introduce a convenient decomposition
in the total displacements field over the RVE domain as a
sum of a linear displacement component ui and a
displacement fluctuation i which represents local
variations about the linear displacement component and
does not contribute to the macroscopic scale strain.
By taking into account the Hill-Mandel Principle of
Macro-homogeneity [10, 11], which establishes that the
macroscopic stress power must equal the volume average
of the microscopic stress power over , the virtual work
equation for the RVE can be reduced to
1
= 2 ( + ) ,
(1)
3 WOOD MICROMECHANICS
At the scale of some few nanometres, wood contains three
fundamental phases: cellulose, hemicellulose and lignin.
These three fundamental constituents form the wood cellwall composite material whose basic unit building block
is called microfibril. This unit can be considered as a
periodic block of rectangular cross-section with infinite
length. This composite comprises reinforcing microfibres which are oriented mainly in a single direction (in
almost the whole cell-walls volume) and are embedded
periodically in a softer matrix. The reinforcing microfibre is made up of periodic alternations of crystalline and
amorphous cellulose fractions. The degree of cellulose
crystallinity is defined as the volume fraction of the
crystalline portion of cellulose with respect to the total
volume of cellulose. The high stiffness of the cellulose is
due to its crystalline fraction whereas its flexibility is
provided by the amorphous part. Only its non-crystalline
fraction may absorb moisture and consequently, change
its mechanical properties.
The matrix of the cell-wall composite is made up of
hemicellulose and lignin polymers. Contrary to crystalline
cellulose, hemicellulose is a polymer with little strength,
built up of sugar units. Its structure is partially random,
with mechanical properties highly sensitive to moisture
changes, softening with increasing water content.
Lignin is an amorphous and hydrophobic polymer and its
main purpose is to cement the individual wood fibres
together and to provide inter-fibre shear strength. It is the
most hydrophobic component in the cell-wall, with
relatively stable mechanical properties under moisture
changes. Figure 3 shows a schematic representation of the
cross section of the microfibril and its fundamental
constituents.
4 COMPUTATIONAL MODELLING
In this section, we describe the multi-scale finite element
modelling of CLT. The type of wood chosen for this study
is radiata pine, which has several applications in building
and engineering structures.
As commented in the previous Section 3, three material
scales (at the level of nanometres, micrometres and
millimetres) form the base material for the analysis of the
CLT plate at the structural scale. Therefore, four spatial
scales are interlinked within the present multi-scale
modelling framework in order to predict the structural
behaviour.
A typical finite element mesh of the RVE associated with
the microfibril scale is shown in Figure 6.
5 NUMERICAL RESULTS
In this section we present the numerical results obtained
by the present multi-scale computational approach
enriched with macroscopic cohesive laws.
The computational loading programme consists of a
displacement-controlled load at mid-span where the
deflection increases monotonically from zero up to reach
extensive failure.
Figure 1 shows the computational simulation of the
progressive damage process in the CLT. The
corresponding contour plot shows the damage variable d
that ranges from 0, when an interface point is fully
healthy, to 1, when the interface point is completely
damaged. In our simulation, two parallel cracks are
predicted consistently with experiments (not detailed
here), along with some delaminated regions, mainly
between the bottom and central layer. In our numerical
predictions, we note that after a cohesive interface has
become fully damaged, it is converted into a contact
interface as well.
6 CONCLUSIONS
In this paper, we have investigated the progressive
damage process in CLT plates by a computational
homogenisation approach in conjunction with a CZM.
Four spatial scales have been interlinked within a multiscale framework in order to predict the undamaged
structural behaviour. The CZM has been adopted at the
macroscopic scale to simulate the formation of
interlaminar and interfibre cracks in the material. The
results obtained in this study demonstrate that the present
approach may represent a promising framework to
describe multiple crack fronts or to describe not wellknown crack paths in complex materials and structures.
We believe that the present multi-scale modelling strategy
combined with CZMs, and with the support of
experimental works, can provide a robust platform for
further investigations on the dissipative response of mass
timber structures with a particular view to clarify features
of the microscopic behaviour, and their influence on the
macroscopic response, which are not well-understood at
present. For example, the investigation of non-linear
buckling combined with delamination in CLT walls
subject to compressive loads is relatively straightforward
and is likely to shed more light into the complex
behaviour of such structures. Finally, we remark that
studies are currently under way to predict numerically inplane shear strength of CLT plates based on the present
framework. This will the subject of a future publication.
ACKNOWLEDGEMENT
E.I. Saavedra Flores acknowledges the financial support
from the Chilean National Commission for Scientific and
Technological Research (CONICYT), FONDECYT
REGULAR research project No 1140245.
K. Saavedra acknowledges the financial support from
CONICYT, FONDECYT Initiation into Research project
No 11130623.
REFERENCES
[1] M. Mohammad, S. Gagnon, B.K. Douglas and L.
Podesto (2012). Introduction to cross laminated
timber. Wood Design Focus 22(2), 3-12.
[2] E.I. Saavedra Flores, K. Saavedra, J. Hinojosa, Y.
Chandra and R. Das. Multi-scale modelling of rolling
shear failure in cross-laminated timber structures by
homogenisation and cohesive zone models.
International Journal of Solids and Structures, 81
(2016), pp. 219232.
[3] Q. Zhou, M. Gong, Y. Hei Chui and M. Mohammad.
Measurement of rolling shear modulus and strength
of cross laminated timber fabricated with black
spruce. Construction and Building Materials 64
(2014) 379-386.
[4] D. Peric, E. A. de Souza Neto, R. A. Feijoo, M.
Partovi, A. J. C. Molina, On micro-to-macro
transitions for multi-scale analysis of non-linear
heterogeneous materials: unified variational basis