Determination of Shearing Properties for Tubular Pinewood under Torsion

In this study, tubular pinewood (Pinus sylvestris L.) specimens are tested and shear strain measurements are performed by applying torsion in z direction in the consideration of light weight aircraft engineering. The objective of this paper is to contribute and generate the nonlinear material model in terms of shear modulus presented with power functions under the consideration of nonlinear behavior of wood under torque. Strain gauge measurements are performed for the maximum shear stresses which develop on the tubular specimen, along the radial ( ) r r r in out , , circumferential ( ) in out , φ φ φ and z directions, in a point-wise manner. The data is gathered and examined for the determination of the local variations of empirical shear modulus functions on transversely isotropic surfaces of the specimens. The coordinate dependent shear modulus functions of ( ) z G r φ , ( ) z G φ φ , ( ) z G z φ are derived for ( ) , , z G r z φ φ as the function of r , φ and z , respectively, by analyzing the gathered data. It is proposed to represent the shear modulus functions, ( ) z G φ φ and ( ) z G z φ with the parabolic polynomials, and, to represent the shear modulus function ( ) z G r φ with a linear equation.


Introduction
Wood is modeled as a fiber composite material of transversely isotropic type.The elastic material constants of fiber composites vary along fiber directions.It is known that; the wood material behaves like a nonlinearly elas-tic composite structure under torsion loading.However, in addition to this phenomenon, it also exhibits viscoelastic properties.Transversely isotropic structure of wood is related to the natural growth of a tree.Annual rings of the wood form fairly regular concentric circles in nature, giving the wood cylindrically symmetric structure.Wood material enters the lightweight materials group with its composite structure.Therefore, the use of wood materials will take place between other light composite materials considered in aircraft engineering.
Coordinate dependent material elastic constants which were represented with mathematical equations, were studied analytically by scientists [1]- [4].An extensive body of experimental data according to the variation of in-plane shear modulus ( ) G φ , and the variation of Young's modulus ( ) E φ of plywood panels were gathered and analyzed, in this sense.The mathematical relations between elastic properties of fiber and matrix structure and the wood polymer elastic constants in micromechanical level were studied [5].In the literature we have found out that, the simulation of wood nonlinearity was formulated as equivalent stresses and equivalent plastic strain relationships by using power functions [6]- [9].
Fundamental works on wood compression were represented as continuous functions [10] [11].Beside these, these researchers have used bilinear or tri-linear functions for the wood nonlinear modeling [12]- [17].Generally, 2D analyses with plane-stress or plane-strain assumption were employed by most FE models for wood or wood composites and they were generally unable to completely describe the 3D stress fields.
Few attempts were found in the literature concerning with the 3D analysis of wood and associated constitutive models.In the literature survey, we haven't found an accepted nonlinear constitutive model for wood that describes the change of elastic constants in 3D as the stiffness becomes nonlinear [7] [16].
It is noteworthy that, "Japanese Iosipescu Shear Testing Method" standardized for fiber-reinforced plastics (FRP) was given in the American Society for Testing Material (ASTM) D5379.In this survey, radial and/or longitudinal plane shear behavior was mentioned.
The above mentioned test method was used in measuring the shearing properties of wood type materials [18].The nonlinear material model presents the power functions which represent the stiffness change based on initial stiffness or stiffness of previous recurrences.Interestingly, typical nonlinearity of wood was observed under the compression loading and along with the shearing directions.
According to another study, in the modeling of the lateral buckling strength of large cross-section beam-columns, the values of the shear modulii GLT and GLR should be known.Thereafter, their variation under creep can be easily defined.In the above mentioned study, power law and exponential law were adopted for creep testing on the wood-species and a linear-parabolic model was adopted for the short term behavior.Power and exponential laws have been chosen as creep functions due to their ease of handling properly in the numerical computation of the creep strain rates [19].
In the literature, it is possible to find specific studies on the simulation of wood nonlinearity under shear and compression.The objective of one of these studies was developing a "3D finite element modeling" that corresponds to the wood material behavior [20].
In the literature survey, researchers have focused on both theoretical and experimental studies that based on genetic optimization and optical metrology in order to obtain whole displacement fields of wood-based panels under bending loading.Tests were named as "Non-statically determined tests" because of the absence of unique analytical relations.In these studies, the developing complex stress fields were found as not only depending on boundary conditions, but also on constitutive equations [21].
In order to obtain the balanced modulii in torsion (GJ) of wood-based composite materials, two finite element method (FEM) Fortran approaches (in displacement and stress) were proposed under the Saint Venant assumption.In this study, finite element analysis results were compared with the experimental values which were obtained from the torsion tests on specimens [22].
The Norway spruce was analyzed numerically using the bilinear material law, and the experimental data are available in the literature.According to this study; as the rolling shear of RT plane has shown mostly nonlinearity, it is seldom observed on the LT plane.The results of analyses with modified bilinear parameters reflected good correspondence with the experimental findings.The parameters were found to be nearly equal to the Weibull statistical distributions [23].
In the another study, the test which was named as "Arcan shear test" was used to investigate orthotropic shear modulus values of softwood Norway spruce.Video extensometer method was used in measuring the strain on the critical specimen section, and the determined moduli were evaluated by means of FEM calculations.In this study it was stated that, the 3-parameter Weibull distribution was the most appropriate method for the probabil-istic description of the three orthotropic shear modulii of wood [24].
In another similar publication; in which the Norway spruce was investigated again by means of the Arcan shear test, the data were obtained until failure occurs in each of the three orthotropic shear planes.In this latter study, the stress-strain curves were adapted to the linear, bilinear and Voce models [25].
A new relationship to calculate the shear modulus 12 G in terms of off-axis modulus of elasticity o E of orthotropic specimens such as wood-based panels was proposed.The verification of shear modulus transformation rules to obtain basic engineering constants that agrees with reliable test data and the anisotropic elasticity theory was performed [26].The shear stress versus shear strain expressions of transversely isotropic cylindrical bar under torsion, having a finite length, was described in terms of partial differential equations [27] [28].In the current study, two coordinate systems were used to define the coordinate dependent strain variations of pinewood specimens under torque, applied in z direction.The results which were obtained by performing pointwise incremental measurements, were illustrated as curves of point-wise shear stress zφ τ vs. shear strain zφ γ , along the principal directions at the gauge locations.The torque loading was applied to the specimen in the ( ) r φ − plane at one end, while the other end was fixed.The related data which was gathered from all these points was shown as separate curves for each test case.The purpose was to show the variations of the

Experiment and Method
Tubular pinewood specimens were used in this study (Figures 1(a)-(c)).The average moisture content of the material was measured and tabulated around 6% -7% levels.Meanwhile, the geometry of the specimens were formed by taking fillet depths equal to the radius of the outer portion of the wood shaft (see Figure 1(c)).Thus, it became possible to minimize stress concentrations around specimen headings.The torsion load in Equation ( 1) was applied on the transverse plane (r φ − plane, Figure 1(b)) which had a unit normal vector k along the z direction.
The measured coordinate dependent point-wise shear strain data ( ) generated by the applied torque T .
All experiments were performed in room temperature by making use of uni-axial and tri-axial rosette gauges.These gauges were bonded onto the specimens at definite locations (Figure 1(a)).In order to minimize the deviations in the test results; these specimens were manufactured from a single tree trunk.The forming procedure was performed on the wood specimens in such a way that their fibers lay longitudinally parallel to the trunk.However, there were intrinsic structural differences stemming from the number and radii of annual rings in the part of the trunk from which the specimen was taken out.The experimental set-up was displayed as a schematic diagram in Figure 2. The torsion loading was applied by using TQSM-21 torsion testing machine [31].Data was collected by employing a specialized units [32]- [34].The set-up included torsion test machine and two different data acquisition systems.In this study, the data was simultaneously stored in terms of output voltage difference As a first step, the normal strain values were obtained by reading the voltage differences from the display of the unit, and as a second step, these normal strains were converted into shear strains by using Mohr's circle [35].

Measurements
In the consideration of the resultant curves ( )  [36].Only the forward loading case was considered in the empirical formulation of this study.Results of the experimental studies were analyzed mainly in two steps.Firstly, specific locations were defined along the r , φ , z directions on the specimens and shear stress versus strain ( ) (Each MA lies on the 1-3 principal plane which corresponds to r φ − plane (Figure 1(b)).Measurement Axis Angle MAA is relatively different from the geometrical central axis angle φ .MAA was defined on the 1-3 principal material plane which was attached on the surface of the cylinder.The first strain gauge was attached with 45 +  angle in the principal material axis-1.All other gauges were also located along this direction (Fig- ure 1(a)).Application of gauges in in r , out r , in φ , out φ locations were repeated in similar way.In the rest of the study, z-axis was rotated to be z′ such that it coincided with MA.Invariants 1 2 3 , , I I I remain constant with respect to coordinate transformations.Stress and strain compo- nents were uncoupled along the principal directions for Only along the principal directions, principal stress and strain values of the isotropic materials coincided with the principal values of fiber composite materials on the Mohr's circle [30] [37] [38].Thus, it became easier to understand the behavior of the wood material under torsion and to compare it with that of the isotropic material.Accordingly, the shear modulus functions were derived along the MA directions for the corresponding principal stresses and strains.In literature, the shear modulus, shear strength, and principal strain angle were measured from the shear stress/ shear strain relation obtained by Iosipescu shear tests.Tests were conducted, and the validity of the in-plane shear tests was examined for the thin strips of western hemlock (Tsuga Heterophylla Sarg.) [39].In this work, the torque value T as measured incrementally from the torsion testing machine.Variation of ( ) Experimental results were discussed in the following section in detail.

Results and Discussion
Distributions of shear modulus functions and related empirical formulas were determined for pinewood.Beside this, average and equivalent mechanical properties of shear modulus values avg G corresponding to the average shear stresses zφ τ for different types of wood (pine, oak, chestnut, hornbeam) were measured and results were listed in Table 1.These measurements were performed without using strain gauges and the related data was obtained from the second form of the nonlinear exact tri-linear curve ( ) In the application direction of torsion loading; since the resultant curves ( )  for the tubular specimen.This distribution is obtained through the specified coordinates of MAA at the outer and/or inner surfaces of specimens (see Figures 1(a)-(c)).
In Figure 3, the outer and inner surface ( ) z G φ φ reactions due to the shear stress-strain distributions were formed in reverse orders at the measurement point where z is equal to 200 mm .Through the circumferential direction, ( ) z G φ φ forms a convex curve at the inner surface and ( ) z G φ φ forms a concave curve at the outer surface.2) Shear stresses formed on the outer surface are smaller than the shear stresses formed on the inner surface because shear stresses on the inner surface accumulate densely and shear stresses on the outer surface do not accumulate that densely.

Conclusions
In this research study, variations of the shear modulus of transversely isotropic pinewood specimens under torsion were examined.The set of equations of shear modulus in z φ − and r φ − planes along the principal directions were formulated.All measurements were performed in terms of local coordinate system.As a continuation of this study for the future an extended experimental study is planned.In this future study it is suggested to calculate additional shear modulus functions and related coefficients by using different loading types on different wood species.As is known, the elastic constants' matrix ij C appearing in the Hooke's general law is used to define the elastic constants of materials.In this matrix, it is possible to identify these constants in the form of linear and nonlinear functions by replacing γ γ γ into the Hooke's general law (Equation ( 4)).The equations obtained by re- placing linear and nonlinear shear modulus functions with unknown elastic constants, can be used in theoretical and numerical researches [28].
2) According to our study it is concluded that,  measurements along the 0˚ oriented wood fibers, yield very small positive or negative shear strain values in general.Additionally, forward or backward loadings cause extension or contraction along the fibers and in turn, this generates positive or negative shear strains respectively.4) The negative shear strains which develop at the measurement points as a consequence of forward loading applying can be explained by contraction of the material along the fibers.Thus, behavior of the shear strains are dependent both upon arrangement of the wood fibers and the load applying direction (Table 3).
5) According to the results of measurements along the MA directions using tri-axial rosette gauges, we have concluded that, the relative shear strain changes between 0˚ -45˚ and 45˚ -90˚ angles are not equal to each other.The difference of shear strains between the 0˚ -45˚ directions is smaller than the directions between 45˚ -90˚.
We recommend further additional experiments, in order to enhance the understanding of complex wood behavior under torsion loading.The stress accumulation phenomenon at the critical points along the wood fibers [30].Consequently, these curves were formulated as empirical mathematical equations.In each test case, the coordinate dependent variations of the shear modulus functions ( ) r , φ and z directions, respectively, were generated on the circumferential surface of the tubular specimen.
stresses were established at the forward and backward rotation directions.Directions of forward and backward torque loadings,

−
curves were obtained for each location.

Figure 1 .
Figure 1.(a)The specimen was modeled by employing two different coordinate systems.The first coordinate system was the cylindrical coordinate system with , , r z φ coordinates (global system).The second coordinate system was fiber oriented principal axes system, 1-2-3 or (LRT); (b) MA stands for measurement axis which is taken at 45˚ from fiber direction on the tangent plane; (c) Dimensions of tubular pinewood specimen (mm).

Figure 2 .
Figure 2. Schematic diagram of the experimental torsion test set-up with data acquisition systems, cables, tubular specimen with gauges on it.Additionally, sequence of the experimental operations were also shown.Secondly, the variations of point-wise shear modulus values of each specimen were shown as shear modulii ( ) z G r φ using the inner or outer shear stresses zφ τ and their corresponding principal shear strains.Consequently, shear modulus functions for the principal shear strains and corresponding shear stresses were expressed in terms of , , the backward and forward loadings of torque, respectively.Forward loading applications yield higher failure values than the backward ones.The empirical equations of terms of various functions and were listed in the Table2.The descriptions of the results were explained according to the relative differences of slopes in shear stress-strain ( ) loading applications.
(GPa) versus φ (degree) distribution of data obtained from six gauges at 0 versus z (mm) distribution from four gauges at 65,157,177, 200 mm z = for the tubular specimen.Variation of function

(
GPa) versus r (mm) on the MA angles φ .In these

Figure 4 .
Figure 4. Two plots of

5 )
Fifth group coefficients of the parabolic functions ( )
pinewood (Pinus sylvestris L.) material.Results were summarized below as follows: 1) In our study, we have shown that zφ τ and rφ τ can be written as values can be calculated by substituting the values of the normal strain components , ,

φ
z is equal to 200 mm are greater than the ( ) z G z φ values at z is equal to 70 mm.Similarly, ( ) outer z G φ φ values are smaller than the values at z is equal to 200 mm are larger than ( ) z G r φ values at z is equal to 70 mm.3) In our study it is also defined that, obtained shear modulus values along the 0˚ fibers are larger than the other two 45˚ and 90  directed measurements performed with tri-axial rosette gauges.Results obtained from

Table 1 .
Average shear modulus values of typical wood types obtained without using strain gauges ( )

Table 2 .
Experimentally determined linear and nonlinear shear modulus functions of pinewood.

Table 3 .
Shear modulus linear/nonlinear functions including the statistical values of the coefficients (GPa).