Analysis of poplar timber finger joints by means of Digital Image Correlation (DIC) and finite element simulation subjected to tension loading

Due to the strong demand for wood, and the need to reduce the carbon footprint in construction, fast-growing low-graded planted species like poplar are promising wood for the supply chain in the context of engineered wood products, EWPs. Finger joints constitute a key technology for the manufacture of long structural elements of EWPs. Thus, evaluation of the mechanical properties of finger joints is important when designing these structural members. This paper shows the results of mechanical behaviour in tension of poplar timber of the I-214 cultivar, for specimens with and without finger joints, using experiments monitored by DIC (Digital Image Correlation) and simulated by the Finite Element Method. For moderate and intermediate loads, the samples with finger joints showed behaviour similar to those without joints. However, the presence of the fingers decreased the mechanical properties of the modulus of elasticity in the longitudinal direction anywhere from 7.7 to 23.7%; and there was a decrease of around 27.5% in tensile strength. An agreeable correlation between the DIC results and FEM simulations is obtained for the longitudinal, transversal, and shear strain fields, thus demonstrating the high potential of both methodologies.


Introduction
Wood is a natural organic, anisotropic, and renewable material with excellent ecological properties, acting as a carbon sink and characterized by a low embodied energy. Timber is one of the most important materials used over centuries within the construction sector (Harte 2016). Nowadays, because the timber market is continually increasing, the need for elements with special shapes and higher spans has become overwhelming. Complex connection-mechanically or adhesively bonded-is a common demand. Finger joints are currently used to produce engineered wood products like glued laminated timber beams (glulam) for different construction sectors, including buildings, bridges, or sport halls.
A main benefit/advantage of implementing finger joints is the possibility to avoid certain types of flaws (e.g., grain deviation, knots, cracks, etc.) and geometrical imperfections.
Splicing two timber pieces together has been a challenging and difficult task since the early stage of finger joint applications, and the manufacturing process is highly industrialized nowadays. Finger joint parameters are, therefore, observed to determine which are more closely related to the strength of a given specimen. Many compressive research works (Özçifçi and Yapıcı 2008;Hernández et al. 2011;Rao et al. 2012;Yeh and Lin 2012;Morin-Bernard et al. 2021) have been carried out, employing different finger joint lengths and pitch lengths, resulting in finger joints with different slopes. Experimental results have shown that the length of the finger is one key parameter affecting finger joint strength (Ayarkwa et al. 2000): fingers with relatively flat slopes and sharp tips could achieve a high strength finger joint, even for finger lengths shorter than the usual length in current practice. The applied end pressure is likewise a key parameter for the evaluation of a specimen's ultimate tensile strength, according to Bustos et al. (2011); under lower end pressure, numerous air bubbles could trigger a significant reduction in the tensile strength. It is also important to remark that higher pressure affords the greatest strength of the finger. However, Ayarkwa et al. (2000) noted that after a certain value (12 MPa) for the applied end pressure, the tensile strength of the finger joint does not significantly increase. There is moreover a risk of failure of the specimen by splitting.
In addition to these parameters, the stiffness and thickness of the adhesive would be highly relevant for finger strength.
As Groom and Leichti (1994) stated, a stiff adhesive is recommended to reduce longitudinal and radial stress concentration at the finger base, while significantly reducing shear stress within the adhesive layer. Other important aspects regarding the behaviour of finger joints bonded with different adhesives and temperature-dependent behaviour were underlined by Frangi et al. (2012).
Numerical finite element simulations (Serrano and Gustafsson 1999) were developed based on nonlinear fracture mechanics to predict the behaviour of finger joints and to assess the influence of the finger joint strength when defects are located in the bond-line. Various authors (Oudjene and Khelifa 2009;Khennane et al. 2014;Tran et al. 2014;Khelifa et al. 2016) carried out further finite element models for timber members dealing with failure prediction, fracture modes, different damage definition criteria, tensile strength, and geometrical optimization of the finger joint.
The Digital Image Correlation technique (DIC) holds great possibilities for exploring full-field displacements and strain measurements in the realm of experimental mechanics. As a non-contact, optical method, DIC has proven to be an ideal tool for the study of material behaviour for a large range of products (Lava et al. 2009). DIC is traditionally used to analyse steel, concrete and thermoplastic behaviour, but not wood material (Quanjin et al. 2020).
In recent years, however, considerable and comprehensive numerical and experimental approaches have been used to address timber fracture characterization by means of DIC and the finite element method (FEM) (Dubois et al. 2012;Milch et al. 2017;Ostapska and Malo 2021a, b). As is well known, the joints are the most critical part of most of the timber structures, being designed as mechanical or adhesively bonded joints. For that reason, joint behaviour is an interesting research area for the application of DIC (Angelidi et al. 2018). Specific applications in determining the Poisson ratio of a ductile adhesive in tension and compression have been undertaken through DIC (Angelidi et al. 2017).
The present study is focused on experimental and numerical investigations of the I-214 poplar cultivar behaviour in tension, with and without finger joints. The study is centered on the linear elastic part of the analysis, since the modulus of elasticity is the main variable of interest for characterization of the specimens. However, to take into account the influence of the contact and cohesion at the finger joints, nonlinear contact and cohesive models that behave almost linearly until failure were used. The experimental campaign involved the optical evaluation of the specimen during testing by means of DIC. The main contribution of the paper is to characterize the mechanical behaviour of the poplar timber with and without finger joints through an experimental analysis supported by DIC monitoring and FEM simulation.

Poplar wood and adhesive
The wood used for this work was extracted from a 9-yearold-poplar plantation of the cultivar I-214 (Populus x Euroamericana [Dode] Guinier) located nearby the city of Granada, Spain. From the logs, planks of 2000 mm length and a section of 35 × 75 mm were sawn and dried for 6 months in natural conditions, always ensuring proper ventilation and avoiding direct exposure to the sun or rain. The final moisture content (MC) of wood was 10 ± 2%. To manufacture the finger joint samples, a monocomponent polyurethane adhesive was used. The density of the adhesive was 1.12 g/cm 3 and its viscosity (Brookfield RVT) was 7000 mPa.s, spindle-6, with a rotational speed of 20 rpm and with a temperature of 20 °C, measured between 16 and 36 h after production (ISO 2555:2018).

Sample description and manufacturing
A set of 20 poplar specimens were prepared to be subjected to tension in the grain direction. The samples were divided into two batches of 10 specimens each, without finger joints (TT) and with finger joints (TF), respectively, without being grouped into strength classes. All the specimens were sawn to a length of 586 mm, with a general cross-section of 50 × 17 mm. So as to promote failure in the central area of the specimen and to avoid undesired breakage, the width of this central part was narrowed to 20 mm, resulting in a mid-section of the specimen of 20 × 17 mm (Fig. 1). For the TF samples, the finger joint was located at the geometrical centre. The adhesive for the finger joint was applied at a temperature of 20 ± 2 °C, with a wood MC of 9.8 ± 1.5%, complying with the requirements of the manufacturer for its application. The adhesive was applied at a rate of 250 g/ m 2 . After that, the specimens were prepared and introduced in a climatic chamber at 20 °C and 65% relative humidity, according to the EN 408:2010 + A1:2012 (2012) standard.
As Jokerst (1981) stated, the geometry of the finger joint dictates the strength of a joint, and consequently the strength of the specimen. Among the research findings reported to date, the finger length is found to mostly affect the strength of the finger joint, while the pitch has a negligible effect on joint strength. In addition, tip thickness and slope of the finger bear a noticeable impact on joint strength (Strickler 1 3 1980). According to EN 14080:2013 (2013)-Annex I, the geometric parameters that define a joint include the finger length (l j ), finger pitch (p), finger angle (α), tip gap (l t ), and tip width (b t ) (see Fig. 2). These parameters are interrelated, meaning that a change in one parameter will influence all the others. This complicates research of the effect of any single parameter on strength (Jokerst 1981). According to studies on the performance of finger joints (Hu et al. 2012;Habipi and Ajdinaj 2015), good performance in terms of tensile strength and modulus of elasticity can be attained for finger lengths of 13-14 mm; hence, a length of 14 mm was selected for the specimens used here (Fig. 2). Afterwards, the geometric parameters of the finger joint were set to fulfil the requirements stated in Annex I from EN 14080:2013 (2013) (see Fig. 2). Taking into account the length of the finger joint, a final pressure of 10 N/mm 2 was applied, in view of the graph from Annex I (EN 14,080:2013).

Experimental test
Tensile tests were carried out in a multi-testing machine, specially manufactured (ad hoc) by Microtest S.A Company, with an electrical actuator having a maximum capacity of 200 kN, with an accuracy of 0.01 kN, and a speed test rate of 0.5 mm/s. To avoid slipping during the test, the samples were clamped 90 mm on each side. Figure 3 shows the setup of the experiment with the DIC equipment during a test. The where F max is maximum axial load, and A is the cross-section of the specimen at the mid-length.

Digital image correlation
For monitoring the strain field in the central part of the specimens, a three-dimensional image correlation non-contact optical measurement system was used, Aramis 3D® (GOM mbH, Braunschweig, Germany). The ARAMIS 3D System is a powerful non-contact measuring instrument that is suitable for assessing three-dimensional deformation and strain distributions of real components under static or dynamic loads.
The surface of the specimen should be smooth and clean. Then, a black-on-white random speckle pattern was applied as uniformly as possible, using an opaque spray paint on the face of each sample.
The DIC system was calibrated using the specific Aramis 3D calibration protocol, which includes camera positions and lens distortion parameters. With this calibration information, the system ensures measurement accuracy. The DIC monitoring area is defined by a face size with 19 pixels, as well as 16 pixels for point distance.
The optical module of Aramis (made up of two 12 M industrial cameras) was oriented towards the specimen The placement of the Aramis optical module was also conditioned by the recommended calibration distance to the specimen, in this case 295 mm, with an angle of 25º between cameras. By means of this technique and the capabilities of GOM Correlate software, developed by GOM GmbH (2020), the modulus of elasticity in tension was calculated using two methods: (I) Based on the relationship established in the standard EN 408:2010 + A1:2012 (2012); (II) based on the average stress-strain relation, within the defined area.
For method I, according to st andard EN 408:2010 + A1:2012 (2012), the following relation was used: where F 2 − F 1 are the increment of load on the straight line portion of the load deformation curve at 20-30% of the maximum load, w 2 − w 1 are the increment of the corresponding displacement, l 1 is the length of the extensometer and A is the cross-section of the specimen at the mid-length. To derive accurate and reliable results, three pairs of virtual orthogonal extensometers with different lengths were placed in three distinct locations (Fig. 4-left). The first extensometer, Ext.1, was placed along the evaluation area (100 mm length) in accordance with EN 408:2010 + A1:2012 (2012) for E Ext.100 t,0 computation. The second extensometer, Ext.2, was placed at the top end of the gauge length at a length of ) for control of the results. To evaluate the influence of the finger, a third extensometer (Ext. 3) was placed within the middle of the finger joint at a length of 14 mm to calculate E Ext.14 t,0 . In all cases, to calculate Poisson's ratio, horizontal virtual extensometers 15 mm in length (Transversal and Transversal 1) were placed at the midlength of the vertical ones. Thus, the relationship between transversal and longitudinal strains could be computed.
In the case of method II, since the DIC technique allows one to determine the mean strains of a defined area with different geometric sizes (Fig. 4-center and right), several areas were considered, starting from the mid-gauge length to the ends of the specimens. To determine the mean longitudinal strain corresponding to these areas, for the specimen with no finger joint, a gauge area corresponding to the entire gauge length was defined, allowing for computation of the corresponding modulus of elasticity (E ga t,0 ) (see Fig. 4-center). For the specimens with finger joints, the modulus of elasticity is calculated in several distinct areas: E c.a t,0 corresponding to the central area, E a1 t,0 to "Area 1", and so on (see Fig. 4-right). These moduli are computed as the slope of corresponding stress-strain relation between the interval of 20%-30% of the maximum load (EN 408:2010 + A1:2012).
For both methods, Poisson's ratio of the specimen was computed. To this end, two different values of Poisson's ratio were obtained employing methods I ( ) using the extensometer, Ext. 3, and II ( * ) using the defined central area, respectively, as follows, where Δ transv is the transversal strain increment (method I) or the transversal average strain (method II), and Δ axial is the longitudinal strain increment (method I) or the longitudinal average strain (method II).

Finite element method: numerical simulations
A 3D parametric finite element model to investigate the behaviour in tension of the specimens with and without finger joints was developed employing Abaqus software (Abaqus CAE 2020) and the programming language Python (Python Core Team 2015). Figure 5 depicts the finite element model of the sample with finger joints. This 3D solid linear material elastic anisotropic model is composed of two glued parts that interact along the finger interface by contact. The mechanical behaviour of the interface is governed by contact and damage cohesive traction laws from Abaqus (Camanho et al. 2003;Abaqus CAE 2020). The failure criterion used for this study is the quadratic traction law.
The boundary conditions are defined according to the experimental test set-up; hence, they were modelled as clamped conditions applied to the all-lateral area of the bottom flange of the specimen. For the top flange, however, the boundary conditions allowed the imposed longitudinal displacement of the specimen with a constant rate of 0.4 mm/s, similar to the experimental part. In the case of the TF samples, the failure of the model is brittle due to the opening of the contacts between fingers. Owing to the sharp geometry of the fingers, a four-node tetrahedral element (C3D4) was used for the modelling of the specimen. The meshing process was refined in the surrounding finger area, with a mean element size of 0.8 mm for the fingers. For the rest of the gauge length, the mean size was 2.4 mm, and the coarse size corresponds to the tapered part and clamping area with mean sizes of 6.4 mm and 7.4 mm, respectively. The size geometry of the meshing process was changed based on converge analysis until the results were stable. The constitutive parameters as moduli of elasticity, Poisson's ratio, and the contact parameters are defined in Tables 1 and 2. The longitudinal input elastic parameter (E 1 ) was taken from the performed experiments, while the rest of the elastic input parameters were considered as a ratio of 1/8 for E 2 /E 1 , and 1/30 for E 3 /E 1 , taken close to the ratios corresponding to "Yellow Poplar" wood (Table 5-1 from Wood Handbook 2010). In turn, the input parameters for the contact interface were defined in view of previous, similar numerical simulations (Tran et al. 2014). An iterative calibration process modifying the initial parameters from the reference was performed by trial and successive verification until the numerical model had fitted the experimental relation.

Experimental results
All the manufactured specimens, TT and TF, were subjected to the static axial tensile test. Figure 6 shows the stress-strain relations based on method I using extensometer with 100 mm length (top), based on method II, using mean longitudinal strains from the central area ( Fig. 6 bottomleft), and method I, using the extensometer of 14 mm length (Fig. 6 bottom-right).
All the relations highlighted linear-elastic behaviour without any evidence of plastic deformation until failure of the specimens, so that the failure is characterized by brittle failure. Clearly the TF specimens exhibit lower stiffness than the TT specimens, identified by a slight curvature within the stress-strain relation induced by the non-linear behaviour of the finger joint.
For the TF specimens with an extensometer or area located at the finger joint, the relations are slightly lower due to the strain concentration captured by DIC at the base of the fingers. Therefore, this local effect produces a slight reduction in the moduli of elasticity. Tables 3 and 4 show the longitudinal moduli of elasticity using methods I and II, for the extensometer of 100 mm length ( E Ext.100 t,0 ) according to EN 408:2010 + A1:2012, the corresponding area ( E ga t,0 ) of the extensometer located in the finger joint or mid-gauge length for the specimen without finger ( E Ext.14 t,0 ) and its corresponding central area ( E ca t,0 ). The third extensometer is located away from the finger joint at 20 mm length ( E Ext.20 t,0 ) and a corresponding area that coincides with the extensometer location ( E a2 t,0 ). It is important to remark that in Table 4, the results computed based on the gauge area with method II are not meaningful, given the reduced size of the finger area (about 14 mm length) compared to the entire gauge area (area of finger joints is 11% of the gauge area). Thus, the mean longitudinal strain using the gauge area is representative only for the TT specimens, since the strain is homogeneous, giving values that match the results provided by method I using Ext.1 ( E Ext.100 t,0 ). For the TT specimens, the moduli of elasticity are considerably higher than for the TF specimens, considering that the finger joint introduces a weakness in the longitudinal behaviour of the specimen.
For the TT specimens, Table 3 shows that the mean values obtained by means of method I and method II are comparable, ranging between 8495 N/mm 2 and 9120 N/mm 2 . For the specimens with finger joints (TF), Table 4 shows that the mean values employing methods I and II have a wider range, 6942 N/mm 2 -7736 N/mm 2 , possibly induced by different strength classes of the top and bottom timber pieces. These higher differences related to the finger specimens are directly linked with the strain acquisition using method II, which is able to capture the high strain concentration. Consequently, the mean longitudinal strains are high, resulting in low longitudinal moduli values. Under method I, the displacements are computed by integration of the strains, for which reason the influence of the concentration of strains is averaged.
The last row in Table 4 shows the moduli reduction between TT and TF, which varies from 11.0 to 19.0%, indicating a relatively low dispersion of the results that is mainly induced by the finger joint and the lack of timber part classification in strength classes. The maximum reduction (19.0%) is found for the measurements located in the central area (finger joint) using method II. In this local measurement, high strains are captured by DIC, increasing the mean longitudinal strain. Hence, the computation of the modulus of elasticity in this area provides reduced values. For the same location, the reduction is very similar to the one provided by method I, with a variation of 0.5%. This trend is visible for all of the cases since the virtual extensometers do not take into account the local nonlinear strain concentrations developed at the fingers. With regard to the obtained results using methods I and II for the top part of the specimen (Ext. 2 and Area 2), the maximum mean reduction is 16.4% and 11.0%, which highlighted the highest variation (5.4%) between methods. Consequently, the mean variation value of modulus of elasticity in tension, evaluated according to the   TT2  7611  7562  7793  7827  7247  7218  TT3  7533  7421  7926  7481  7590  7602  TT4  7531  7468  7740  7494  7396  7853  TT5  8627  8545  8124  8294  8275  8368  TT6  8661  8695  8990  8741  8844  9557  TT7 10 , (2012) standard, is around 14.6%, mainly due to the finger joint presence. Table 5 summarizes the average value of the experimental results and the mean difference related to the Poisson ratios obtained using methods I-II (ν and ν*)-using the extensometer with 10 mm length (Ext. 3) and central area, tensile strength (f t,0 ) and density (ρ) of the TT and TF specimens. As the current results revealed, the finger joints do not have any influence on Poisson's ratio value, because the fingers do not significantly affect the transversal strains. Yet a significant reduction (27.5%) in the tensile strength is observed for the TF specimens, because the weakness induced in the longitudinal direction triggered a failure mechanism. According to EN 14080:2013EN 14080: (2013, the tensile strength obtained for the TF samples is higher than that required for the maximum resistance class, GL32h / GL32c. An important aspect deserving mention is that, for both specimens, the mean Poisson ratios computed using both approaches (Methods I and II) are nearly identical (0.0% and 5.0% of variation). Therefore, the proportionality between transversal and longitudinal deformation is similar for both computational methods. Subsequently, the presence of the finger joint does not exert any influence related to density (− 2.6%).
For a better understanding of the longitudinal moduli variation along the gauge length, the mean strains computed based on the defined areas (Fig. 4) -Method II, are  /mm 2 )   TF1  8079  -8449  7817  8082  8925  TF2  7884  -7164  6952  9489  7455  TF3  9036  -7882  8971  10,453  8191  TF4  7182  -7008  7015  8093  6613  TF5  7297  -7741  6835  7088  8193  TF6  9356  -9215  8910  9470  10,840  TF7  5669  -5302  5526  5194  7006  TF8  6468  -5403  5982  7739  7028  TF9  4944  -5486  4923  4246  4391  TF10  6647  -5772  6611 Fig. 7. It is worth mentioning that for the TF specimens, a significant decrease is emphasized in the central area, corresponding to the finger joint location. However, the TT specimens revealed a homogeneous distribution of longitudinal moduli along the gauge length, mainly due to the absence of the finger and therefore to the one-piece manufacturing process of the specimen. It is moreover important to remark that for the TF specimens, different moduli values are observed for the top and bottom part of the specimen, indicating different timber strength classes used in the manufacturing process. In general, for the TT specimens, the failure is mainly caused by crack propagation along the matrix, while initiation of the failure is triggered within the tapered part (transition between gauge and top flange) (see Fig. 8a). Yet for the TF specimens, the failure is brittle and the crack initiation is located within the finger joint, developing in the direction perpendicular to the fingers (see Fig. 8b). Figure 9 shows the normal distributions of the mean moduli of elasticity presented in Tables 3 and 4, based on the considered computation methods (I and II). As Fig. 9 shows, the mean longitudinal strain distributions for specimen TT, using method I (left), are similar and identical from the pattern variation point of view, while for TF specimens, a visible decrease and a higher dispersion are revealed by the distribution (Fig. 9left).
This heterogenic aspect of the results is likewise observed in the results based on method II (Fig. 9right). Figure 10 shows the experimental longitudinal strain evolution for one representative TT (top) and TF (bottom) specimen with intermediate behaviour. In addition, the corresponding DIC strain map is shown at four loading stages: 25%, 50%, 75% and 100%, respectively.
As the DIC strain map shows, for low and intermediate load step levels (up to 50% of the maximum load), the strain field is homogeneous for both specimens. For the TF specimen, the force-elongation relation also captures the stiffness reduction, characterized by a reduced slope compared to the TT relation.
For higher load levels, several inhomogeneities start to appear for TF around the fingers due to the non-linear behaviour of the finger joint, while for TT the field is still homogeneous in the mid-section of the specimen. For loads closer to the fracture point, the strain field remains homogenous for TT, but for the TF specimen, a significant concentration

Finite element simulation and DIC results
For the discussion of this subsection, a representative experimental sample and a calibrated numerical model were selected. Figure 11 presents the force-elongation relations for the TF and TT specimens obtained from the experimental and numerical simulation part. The relations from the experimental part are plotted in a representative sample with an intermediate behaviour based on the relations provided in Fig. 6-top. The numerical simulations through the calibration process matched the experimental relation and the failure point of the specimens. The calibration process started with the experimental parameters; then, through successive steps, the final values were obtained. These parameters of the FEM calibrated models are presented in Tables 2 and 3. Figure 12 shows a comparison of the longitudinal (Ɛ x ), transverse (Ɛ y ) and shear (Ɛ xy ) strain fields between DIC and FEM before failure. Both techniques capture the high strain concentration around the tip and base of the fingers due to the initiation of detachment of the fingers and the shear distortion along the lateral faces. The maximum longitudinal strain concentration occurs within the base of the fingers, giving similar values for DIC and FEM (0.8%). It is noteworthy that longitudinal strains within the tip are almost null due to the gap produced during the manufacturing process. The maximum longitudinal strain (Ɛ x ) occurs at the base owing to the action of the cohesive forces along the lateral sides of the fingers, the strain flowing through the finger base. Regarding the transversal strains and shear strains, distributions by FEM and DIC are in accordance, showing low magnitude in comparison to the longitudinal one, with a maximum transversal strain of 0.3%, and 0.12% for shear strain. For the transversal (Ɛ y ), the flow is symmetric and is mainly caused by the Poisson ratio, producing a small necking due to a high Poisson's ratio (value of 0.4 in Table 5). For the shear (Ɛ xy ), the strains at the fingers are equal and opposite, and the maximum strains occur at the tips of fingers because of the reduced size of the section of the finger close to the tip. While the DIC shear strain distribution suggested small negative shear strains, which can be produced by the imperfections between the top and bottom clamps of the machine.

Conclusion
This study analyzed and characterized the mechanical behaviour of the poplar timber of I-214 cultivar for specimens with and without finger joints. Based on the obtained results, the poplar wood behaviour subjected to axial tensile stress was analyzed in parallel with the numerical simulations.
The experimental results obtained employing different methodologies using DIC technique allow the possibility to analyze the samples along their entire length. As the results show, for the specimens with fingers joint, the reduction in the stiffness and the tensile strength becomes relatively small. The results show that the presence of the fingers decreases the mechanical properties between 11.0 and 19.0% for moduli of elasticity, and around 27.5% of the tensile strength, respectively. However, the tensile strength remains higher than the limits specified in EN 14080 standard.
One of the most important experimental results obtained by means of DIC technique is the evaluation of the longitudinal modulus of elasticity in tension along the gauge length of the specimen. Consequently, two particular behaviours are identified, a homogeneous moduli of elasticity distribution along the gauge length for the TT specimens, while for the TF specimens, the results were highly dispersed, with a minimum value located within finger joint. It is worth mentioning that both methods provide similar results for the strain gauges placed within finger length and the corresponding area of the finger surface. Nevertheless, the central area is averaging the strains of the entire length of the finger, which revealed high parts with high strains concentration. Therefore, the low variation of the mean reduction obtained employing both methods, confirms the validity of both investigation methods. Fig. 12 Strain field distribution for TF specimen. a Longitudinal strain (Ɛ x ); b Transversal strain (Ɛ y ); c Shear strain (Ɛ xy ) from DIC. d Longitudinal strain (Ɛ x ); e Transversal strain (Ɛ y ); f Shear strain (Ɛ xy ) from FEM Besides, the DIC technique is a powerful tool that allows evaluating the strain state at every rate load, and the strain concentrations and flow within the finger geometry.
Furthermore, numerical simulations were carried out to obtain a suitable FE model, able to capture, verify and complement the structural behaviour provided by DIC, and showing good agreement. Through these numerical simulations, the main aim was to obtain the parameters related to the contact interface, which can be extended to a further study related to the optimization of fingers joints. The results provided through this study, are useful for further research tasks, which involve glued poplar beams or other engineered wood products, where the tension behaviour of the raw material has yet to be known.