Issue 
Mechanics & Industry
Volume 18, Number 8, 2017
Experimental Vibration Analysis



Article Number  804  
Number of page(s)  8  
DOI  https://doi.org/10.1051/meca/2017042  
Published online  21 March 2018 
Regular Article
Using experimental modal analysis to assess the behaviour of timber elements
^{1}
University of Mons, Faculty of Engineering, Department of the Theoretical Mechanics, Dynamics and Vibrations,
Place du Parc 20,
7000
Mons, Belgium
^{2}
University of Mons, Faculty of Engineering, Department of Civil Engineering,
Place du Parc 20,
7000
Mons, Belgium
^{*} email: georges.kouroussis@umons.ac.be
Received:
12
December
2016
Accepted:
6
November
2017
Timber frameworks are one of the most important and widespread types of structures. Their configurations and joints are usually complex and require a high level of craftsmanship to assemble. In the field of restoration, a good understanding of the structural behaviour is necessary and is often based on assessment techniques dedicated to wood characterisation. This paper presents the use of experimental modal analysis for finite element updating. To do this, several timber beams in a free supported condition were analysed in order to extract their bending natural characteristics (frequency, damping and mode shapes). Corresponding ABAQUS finite element models were derived which included the effects of local defects (holes, cracks and wood nodes), moisture and structural decay. To achieve the modal updating, additional simulations were performed in order to study the sensitivity of the mechanical parameters. With the intent to estimate their mechanical properties, a procedure of modal updating was carried out in MatLab with a Python script. This was created to extract the modal information from the ABAQUS modal analysis results to be compared with the experimental results. The updating was based on a minimum of unconstrained multivariable function using a derivativefree method. The objective function was selected from the conventional comparison tools (absolute or relative frequency difference, and/or modal assurance criterion). This testing technique was used to determine the dynamic mechanical properties of timber beams, such as the anisotropic Young's Moduli and damping ratio. To verify the modulus, a series of static 4point bending tests and STS04 classifications were conducted. The results also revealed that local defects have a negligible influence on natural frequencies. The results demonstrate that this assessment tool offers an effective method to obtain the mechanical properties of timber elements, especially when onsite and nondestructive techniques are needed, for example when retrofitting an existing structure.
Key words: Structural behaviour assessment / modal analysis / timber beam / NDT / FRF measurement / modal updating
© AFM, EDP Sciences 2017
1 Introduction
Wood is widely used as main material in building construction due to its good thermal and acoustic properties, and to its aesthetic. Contrary to other building materials, wood is heterogeneous, anisotropic and hygroscopic. It is mainly composed of 50% carbon, 43% oxygen, 6% hydrogen and 1% nitrogen in mass, against the dry timber mass (all the carbon contained in the wood comes from the carbon dioxide used for photosynthesis). Beginning with the inside layer of a tree (Fig. 1), the wood is composed of:

Heartwood or duramen: This part is called “perfect wood” and only contains dead cells which ensure the role of support for the tree. It can be distinguished by its colour, as it is much darker than living wood.

Sapwood constituted of living cells: This is the younger, outermost wood.

Based generator or cambium: Located directly under the cork, the cambium produces the outside of the internal cork and the inside of the wood cells. In temperate zones, it works only in spring and summer and produces a layer of wood every year, called an annual growth ring, through which it is possible to distinguish the initial wood (spring) and the final wood (summer).

Internal cork is the inner part of cork and conducts the sap containing sugars produced by photosynthesis.

External cork, more commonly called cork, which is composed of dead cells and acts as a protective coating.
Wood properties are affected by the climate and the conditions the tree's growth. Wood is also a porous material which can lose or regain moisture according to the relative moisture, and it can eventually store water for its cellular cavities. The variations of humidity influence a piece's dimensions and its mechanical strength. In general, the mechanical properties are considered for 12% of moisture content. As wood is a highly variable natural material, susceptible to physical, chemical and biological deterioration, it is very difficult to assess it from a mechanical point of view. Before any intervention on existing structures, the first step is to fully understand how the structure works, meaning that the knowledge of its mechanical properties is essential.
Sounding methods are the oldest methods used to inspect in situ timber members [1]. It consists in exciting the timber beams with a hammer. This requires a trained inspector able to interpret the sound. This method allows a quick screening of a timber member but is highly qualitative and not quantitative. The diagnosis may also change according to the inspector. The ultrasonic method offers a way to control the quality of timber product production [2]. It consists of measuring the time difference between the arrival of 2 piezoelectrical transducers in a frequency range of 20–500 kHz with an accuracy of 1 μs. The transducers are located on the two opposite sides of the timber piece. The data acquisition system consists of a signal detector, an analyser and a recorder which collects and retains the wave information. The measured propagation velocity depends on numerous factors, such as the species, moisture content, grain direction and density, and it allows the determination of Young's modulus as well as any possible defects in the wood. New generations of devices deal with the analysis of the acoustoultrasonic response of the wood by measuring the speed of the transmitted wave and their maximal peak of energy. The speed of propagation is still correlated to Young's modulus but the energy is related to the local singularities, like knots or grain direction [3]. Other assessment techniques are available and their use must be discussed depending on the remaining level of doubt after visual inspection, the feasibility of repair/strengthening, and is always in combination with economical considerations.
Complementary to the existing nondestructive testing methods, the purpose of this work is to propose the first step in the development of an efficient technique for the mechanical assessment of old timber structures. The chosen method is the vibration technique, due to the simplicity of its implementation and the fact that it gives an overall view of the surveyed wood element. This paper presents the results investigated for the experimental assessment of the dynamic properties of timber elements, measured by forced vibration tests using modal identification methods.
Fig. 1 Constitution of wood. 
2 Existing vibration techniques for wood characteristics assessment
The vibration method consists of exciting the structure by hitting it with a hammer at several points of impact in order to induce a stress wave into the wood specimen [1]. The acceleration signal is measured using small accelerometers (Fig. 2). When applied to wood inspection, the high attenuation of wood, which worsen with high frequencies, must be considered. However, high frequency transmissions are more sensitive to internal defects such as knots or other singularities. The choice of the frequency range results from a compromise between the level of attenuation (small or great spans) on the studied length L and the sensitivity to defects. The frequency spectrum analysis can be used to assess and quantify timber decay. When fundamental modes are of interest, some formula exist to establish close relationships between Young's modulus and the studied fundamental mode.
The mechanical properties of a timber structure can be experimentally determined by forced or ambient vibration tests using dynamic properties or modal identification methods. These methods have been successfully used in several structural applications. Yang et al. [4] proposed to use the shape of the flexural vibration wave of wooden beams at the first mode as a primary indicator. The dynamic modulus of the elasticity of the studied beam sections was calculated from the ratio of the curvature of the corresponding mode shape. By using a cantilever beam configuration, Wang et al. [5] detected the first bending mode of the beam, using a simple technique based on the measurement of impulse response. Both Young's modulus and the damping ratio were estimated. As illustrated in Figure 3, the main difficulty is to isolate the first bending mode among all the recorded bending modes (and eventually torsional modes if the symmetry is not verified). It is clearly shown that a filtering becomes unavoidable, with the lack of awareness of its cutoff frequency. Moreover, the proposed method needs to be performed on wood species under laboratory testing conditions. A similar method was used by Wang et al. [6] and dedicated to in situ measurement. An analytical model was derived from the flexural beam theory to correlate the associated eigenfrequency. An alternative nondestructive impact test, based on sectional mode shapes, was presented by Alves et al. [7]. This allows the determination of transverse Young's modulus of wood species.
From this short literature review, it is seen that existing vibration techniques are mainly focused on one flexural mode, and provide one or two dynamic parameters. The effect of each mechanical parameter on the natural frequency is also disregarded. Another finding is that one specific mode shape is expected to be analysed without verifying that the correct mode is analysed. As suggested by Olsson and Källsner [8], a calibration with a finite element model reduces the uncertainty in the mode estimation. However, the experimental data needs to be completely mastered by identifying the various mode shapes well.
Fig. 2 Device generally used for the vibration method. 
Fig. 3 Example of a measured acceleration on a cantilever wood beam due to an impact. 
3 On the use of frequency response functions
In mechanical engineering, model updating techniques are widely used in order to calibrate models in structural dynamics. The potential of model updating techniques has close links to inverse parameter identification, as long as the obtained values are inside a realistic range. In this paper, it will be shown how finite element updating with wood beams can be used. The feasibility of experimental modal analysis on the characterisation of timber element dynamic properties will also be described.
A series of beams (more than twenty specimens) of the same species (Douglasfir) and dimensions (2.4 × 0.14 × 0.04 m) were analysed. Each timber beam was defined according to the fibre direction (along the beam length). Two configuration setups were retained (Fig. 4):
To avoid any influence related to the boundary conditions, a free condition was reproduced by supporting the beam elements with very soft rubber bands (Fig. 4a). To ensure a freebased condition, the natural frequency of the last rigid mode should be less than onetenth of the natural frequency of the first flexible mode [9].
A cantilever beam configuration was then adopted, with specific attention being paid to the oneside fixed condition. Since a perfect clamping was physically impossible without spoiling the wood surface (and therefore locally modifying the beam properties), neoprene pads were used between the beam and the support. This induced a unknown linear stiffness k_{s} of the pad which then had to be estimated.
The roving hammer technique was applied, by exciting the structure vertically on all the selected points with an instrumented impulse hammer. Three piezoelectric accelerometers were vertically placed on the beam surface. The influence of the mass of these added accelerometers on frequencies is negligible, since their total mass was less than 5% of the weight of the timber beam. The obtained result was a series of measured accelerations at point i for each excitation reference at point j, establishing the corresponding frequency response functions (FRFs) H_{ij} (f) as a function of the frequency f. The applied frequency range was from 0–1024 Hz with a resolution of 1 Hz. An exponential window was used to minimise the effects of leakage. According to the BettiMaxwell reciprocity theorem, the excitation points j and the response points i can be switched (H_{ij} (f) = H_{ji} (f)) so as to obtain a set of FRFs adapted to modal analysis. The leastsquare complex exponential (LSCE) method, implemented under the inhouse EasyMod toolbox [10,11], was used to determine the modal parameters (eigen frequencies, damping ratio, modal compliance) of the analysed structure.
Fig. 4 Test support mechanism with hammer in (a) freefree conditions and (b) cantilever installation. 
4 Implemented finite element updating procedure
The EasyMod toolbox was used as it presents an alternative to commercial software packages for applying experimental modal analysis. It consists of several MatLab functions, allowing the step by step establishment of a complete modal analysis from experimental data. This tool has the advantage of being open source. The structure of EasyMod functions is presented in Figure 5. Both identification and validation methods are implemented: circlefit, linefit and LSCE methods, as well as modal assurance criterion (MAC) and modal collinearity. By using the matrix abilities offered by MatLab, operations on FRF are made much easier.
The finite element analysis was performed using ABAQUS v.6.10 software. The model is a 3D beam element, defined in terms of bulk density ρ and 9 mechanical constants: Young's moduli E_{i}, the shear moduli and Poisson's ratios ν_{ij}. Indexes 1–3 correspond to the longitudinal, radial and tangential directions, respectively. Young's moduli and the Poisson's ratios are related by (1) to ensure the symmetry of the matrix. The direction of the axis of the anisotropic model was specified. In order to extract mode shapes, a set of nodes was defined, corresponding to the measurement setup. For the free boundary conditions, the hung rubber bands were modelled as springs whose stiffness k_{b} was to be optimised during the modal updating. For fixed conditions, the pad stiffness k_{s} was considered and included in the updating process. The goal of the procedure was to compare the modal parameters between the experimental measurements and the numerical simulations. To achieve the modal updating, the minimisation of one of several criteria (error functions) was considered [12]:
An algorithm used to achieve the modal updating was implemented in EasyMod. It is based on a Python script created to extract the modal information (frequencies, damping ratios and mode shapes) from ABAQUS modal analysis results to be compared to ones of counterparts obtained from experiment. The updating is based on a minimum of unconstrained multivariable functions using a derivativefree method. The objective function was selected from the three aforementioned criteria using a large number of n modes, so that the dynamic behaviour and the stiffness of the clamping or of the rubber bands could be estimated using an optimization algorithm [13].
With the intent to achieve modal updating, the mechanical behaviour of a timber beam relies on 10 parameters: E_{1}, E_{2}, E_{3}, ν_{12}, ν_{23}, ν_{13}, G_{12}, G_{13}, G_{23} and ρ, plus the boundary stiffness (k_{s} or k_{b}). If all of these parameters are considered for the modal updating, the number of variables is quite exhaustive and solving the problem would be badly conditioned. A numerical study was required to separate and study the influence of each parameter on the natural frequencies. Figure 6 shows the results from this sensitivity analysis, revealing that the effect of E_{1}, G_{12} and G_{13} on eigen frequencies f_{0,i} is mainly positive although the Young's moduli E_{2} and E_{3} are insignificant. For instance, a relative variation of 200% of the Young's modulus E_{2} only causes a change of 0.01% to the frequencies. With regards to the other mechanical parameters, they do not impact the frequencies. Obviously, since only bending modes are studied, Poisson's ratios do not impact the natural frequency results.
Another physical property of wood is the bulk density. This characteristic fluctuates severely according to the considered species, the position of the tree, and more precisely the position of the levied timber on the trunk. This is an important criterion which allows a physical description of wood. Moreover, the bulk density depends on the water content. Most commonly, the bulk density is considered for normal weather conditions (temperature of 20 °C and relative air humidity of 65%). This parameter should not be confused with the true density which corresponds to the ratio between the apparent density, and the porosity (practically identical whatever the considered species, 1530 kg/m^{3}).
Fig. 5 Schematic operating diagram of the EasyMod toolbox, including the finite element updating procedure. 
Fig. 6 Sensitivity study of the mechanical parameters: (a) all the moduli and (b) focusing on shear moduli. 
5 Results
The timber beams were quite large (14 cm) and, due to the effects of knots or other heterogeneities, it was important to check if the side of the measurement, or its opposite side, had any influence on the natural frequencies. After several tests, the differences never exceeded more than 1%, therefore the influence of the surface of measurement may be considered negligible. Other tests were performed, like the effect of local defects (holes or cracks) and it turns out that a small local modification on beam timber does not affect the modal properties.
Once the modal parameters were determined, several procedures that allowed the modal to be validated were used (visual verification, modal vector orthogonality,…). For instance, Figure 7 shows some FRF results obtained from impact measurements and the corresponding FRFs synthesised using the modal parameters extracted from the identification step. Clearly, both curves are accurately matched, serving as a check of the modal parameter estimation process.
The combination of the modal updating was tested for three aforementioned methods, defined by equations (2)–(4). Each of them provided good results with low relative errors. The initial characteristics of the timber beams to initiate the updating was fixed to values found in literature: ρ = 475 kg/m^{3}, E_{1} = 12000 MPa, E_{2} = 9067.4 MPa, E_{3} = 332 MPa, ν_{12} = 0.33, ν_{13} = 0.34, ν_{23} = 0.48, G_{12} = 600 MPa, G_{13} = 515.62 MPa and G_{23} = 53.34 MPa.
With the aim of estimating the advantages and drawbacks of considering E_{1} and G_{12} or E_{1}, G_{12} and G_{13}, the procedure was also carried out for each of these cases. Table 1 highlights the effectiveness of the procedure for the last mentioned case (simultaneously on the three parameters), showing the absolute and relative errors in natural frequencies and the MAC values between the experimental and numerical mode shapes. It is shown that the proposed identification method allows an accurate identification of both inplane and transverse elastic properties. Regarding the use of MAC as a measure of consistency between estimates of a modal vector, it provides an efficient way of determining the degree of causality between the estimates of different modal vectors from the same system. It is simply defined, according to this approach, as follows: (5)
Although other forms of assurance criteria generally offer a better modal parameter estimation consistency [14], the form defined by equation (5) gives a reasonable definition to expected modal vectors and their accurate comparison, as illustrated in Figure 8 showing MAC values of, almost always, greater than 90%.
To verify the obtained dynamic parameter values, two other tests were used in parallel to the vibration tests: a visual inspection of the timber beam based on the Belgian national recommendation STS04 for estimating the mechanical characteristics STS04 [15] (qualitative estimation of E_{STS04} based on a correspondence among various defined grating), and flexural tests carried out according to the standard NBN EN408 for the purpose of experimentally determining the Young's modulus E_{1} [16]. This standard defines a set of conditions to be respected in order to carry out the test and to quantitatively estimate E_{fl}. The range [E_{fl,min} − E_{fl,max}] defined by the last method was also provided and calculated by considering a tolerance of 0.5 mm on the measured deflection. Finally, the modal updating procedure was carried out for several timber beams and the results are presented in Table 2. At first sight, it seems that the results of the modal updating range from the values given by the STS04 guideline for the results of the flexural tests, except for some timbers (the bad results for beam #2 may be explained by the fact that during the flexural test, the timber had a clear trend to torsional buckling). With regards to these results, it appears that satisfactory results are provided by the modal updating, with a small difference between each method. This also pointed out the large discrepancies of the method provided by the standard NBN EN408. Both studied boundary conditions showed similar results and the assessment of the wood dynamic parameter with both cases was well deployed using the modal updating. In addition to the classical method, other parameters can be estimated, such as the damping loss factor measurable for each experimental test (a mean value around 0.003 was found through the different tests).
Fig. 7 Example of measured FRFs with the corresponding synthesised curve obtained from the estimated modal parameters. 
Example of results of the modal updating using method 3.
Fig. 8 Modal assurance criterion: (a) timber #4, (b) timber #5 and (c) timber #14. 
Modal updating − Final results for some beam specimens (free boundary conditions).
6 Conclusion
This paper analysed whether the procedure of modal updating ensures getting the mechanical characteristics of a timber beam using FRF measurements and finite element updating. This makes a link between the wellknow method in mechanical engineering and the material used in civil engineering, for which the structural behaviour is still poorly understood. Based upon the results and findings obtained, the following conclusions were made:

Free or fixed boundary conditions test can be used without any technical constraints and offer similar trends.

Using beam elements allows the estimation of the longitudinal Young's modulus E_{1} and, to a lesser extent, the corresponding shear moduli G_{12} and G_{13}. The other dynamic parameters have no effect on the modal characteristics.

Local defects, like holes or cracks, have a negligible influence on the natural frequencies so that the method could be generalised to use this technique on a large scale.

Satisfactory results were obtained as the existing methods do not provide accurate results. This provides proof of the concept for the modal updating of timber beam using FRF measurement.
To allow the use of this assessment technique on site, the introduction of realistic boundary conditions in the process should be mastered.
Acknowledgement
The authors would like to acknowledge Amandine Meersseman and Pauline Goffinet for their various tests and analyses on multiple timber beams.
References
 B. Kasal, T. Tannert (Eds.), In situ assessment of structural timber, volume 7 of RILEM state of the art reports, Springer, Netherlands, 2011 [CrossRef] [Google Scholar]
 B.M. Lempriere, Ultrasound and elastic waves: frequently asked questions, Academic Press, San Diego, USA, 2003 [Google Scholar]
 J.L. Sandoz, B. Yann, Timber grading machine using multivariate parameters based on ultrasonic and density measurement, in: COST E 53 Conference − Quality Control for Wood and Wood Products, Warsaw, Poland, 2007, pp. 167–173 [Google Scholar]
 X. Yang, T. Amano, Y. Ishimaru, I. Iida, Application of modal analysis by transfer function to nondestructive testing of wood II: modulus of elasticity evaluation of sections of differing quality in a wooden beam by the curvature of the flexural vibration wave, J. Wood Sci. 49 (2003) 140–144 [CrossRef] [Google Scholar]
 Z. Wang, L. Li, M. Gong, Measurement of dynamic modulus of elasticity and damping ratio of woodbased composites using the cantilever beam vibration technique, Constr Build. Mater. 28 (2012) 831–834 [Google Scholar]
 X. Wang, R.J. Ross, M.O. Hunt, J.R. Erickson, J.W. Forsman, Low frequency vibration approach for assessing performance of wood floor systems, Wood Fiber Sci. 37 (2005) 371–378 [Google Scholar]
 R.J. Alves, M.D.C. Magalhaes, E.V.M. Carrasco, Determination of the transverse Young's modulus (TYM) of wood by means of an input power technique, Constr. Build. Mater. 42 (2013) 11–21 [CrossRef] [Google Scholar]
 A. Olsson, B. Källsner, Shear modulus of structural timber evaluated by means of dynamic excitation and fe analysis, Mater. Struct. 48 (2015) 977–985 [CrossRef] [Google Scholar]
 A.G. Piersol, T.L. Paez, Harris' shock and vibration handbook (6th Edition), McGrawHill, New York, USA, 2009 [Google Scholar]
 G. Kouroussis, L. Ben Fekih, C. Conti, O. Verlinden, EasyMod: du développement d'un tool box sous MatLab vers l'enseignement des bases de l'analyse modale expérimentale, in: 3ième Colloque «Analyse vibratoire Expérimentale», Blois, France, 2012 [Google Scholar]
 G. Kouroussis, L. Ben Fekih, C. Conti, O. Verlinden, EasyMod: a MatLab/SciLab toolbox for teaching modal analysis, in: 19th International Congress on Sound and Vibration (ICSV19), Vilnius, Lithuania, 2012 [Google Scholar]
 J. Cugnoni, T. Gmür, A. Schorderet, Inverse method based on modal analysis for characterizing the constitutive properties of thick composite plates, Comput. Struct. 85 (2007) 1310–1320 [CrossRef] [Google Scholar]
 G. Kouroussis, L. Ben Fekih, T. Descamps, Assessment of timber element mechanical properties using experimental modal analysis, Constr. Build. Mater. 134 (2017) 254–261 [CrossRef] [Google Scholar]
 R.J. Allemang, The modal assurance criterion − twenty years of use and abuse, Sound Vib. 37 (2003) 14–21 [Google Scholar]
 Service Public fédéral Economie, STS 04: Bois et panneaux à base de bois, Brussels, Belgium, 2008 [Google Scholar]
 Institut belge de normalisation, NBN EN408: structures en bois − Bois de structure et bois lamellé collé − Détermination de certaines propriétés physiques et mécaniques, 2003 [Google Scholar]
Cite this article as: G. Kouroussis, L.B. Fekih, T. Descamps, Using experimental modal analysis to assess the behaviour of timber elements, Mechanics & Industry 18, 804 (2017)
All Tables
Modal updating − Final results for some beam specimens (free boundary conditions).
All Figures
Fig. 1 Constitution of wood. 

In the text 
Fig. 2 Device generally used for the vibration method. 

In the text 
Fig. 3 Example of a measured acceleration on a cantilever wood beam due to an impact. 

In the text 
Fig. 4 Test support mechanism with hammer in (a) freefree conditions and (b) cantilever installation. 

In the text 
Fig. 5 Schematic operating diagram of the EasyMod toolbox, including the finite element updating procedure. 

In the text 
Fig. 6 Sensitivity study of the mechanical parameters: (a) all the moduli and (b) focusing on shear moduli. 

In the text 
Fig. 7 Example of measured FRFs with the corresponding synthesised curve obtained from the estimated modal parameters. 

In the text 
Fig. 8 Modal assurance criterion: (a) timber #4, (b) timber #5 and (c) timber #14. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.