Numerical Analysis for U-Shaped Thin-Walled Structure Reinforced Timber Beam Based on Thin-Layer Beam Theory

/is paper presents a theoretical model, taking into account the shear deformation subjected to the influence of U-shaped member by geometric parameters as flange height based on thin-layer beam theory, to analyze the structural bending behavior of U-shaped member reinforced timber composite beams, and the feasible design forms of U-section have been pointed out. /e algorithm for this composite beam is the most practical and effective method to meet the accurate solution. /e formulas for the common forms of U-section are presented. It aims to develop a rational engineering approach. /e proposed model has been validated by comparing the results obtained in the present analysis with experimental results and finite element analysis. Furthermore, the results suggested that the value of flange height can be one-fifth the beam height based on the present analysis by comparison of two types of beams. And it is shown that the model provided here correlates consistently and satisfactorily with a wide range of timber beams reinforced by a thin-walled structure such as steel or aluminum alloy sheet bonded to their tension faces.


Introduction
Timber structures' natural features are attractive in applied engineering.As timber structures are brittle materials, different reinforcement techniques can be adopted for avoiding easily damageable parts of the building, such as concrete, steel, and FRP used for timber beams and sandwich beams [1][2][3][4][5][6][7].However, the reinforcements often make beam confront stiffness degradation and bonded layers separation due to interlay slip associated with flexibility and deformation of connectors [8][9][10].Shear stress plays a significant role in bending strength of composite beam as discussed in many research studies, and some numerical methods to estimate it have been established [5,9,[11][12][13][14][15][16].Malek et al. proposed that shear stress in the interface of concrete and plate beams leads to flexural cracks compared to other items including tensile stress, analytically given by unique definition from high-order ordinary differential equations applied to the concrete beam [17].
A system of U-shaped member provided to reinforce timber beam shows a ductile behavior with respect to the plate [18][19][20].e materials such as steel, aluminum, cement asbestos, and GFRP adopted are all available [21].And, a U-shaped member using GFRP fabric was first discussed by eakston [22].
Many publications [11,12,23,24] reported that using the U-shaped member can effectively improve the maximum bending capacity.Rescalvo et al. [19,20] presented data to demonstrate the benefits in flexural strength, and stiffness of timber improved due to the U-shaped member [19,20].e same results were obtained with applications of the U-shaped member reinforced LVL beam in a study by Subhani et al. [25].e improvements of 15.94% and 16.10% in flexural strength and stiffness, respectively, were observed from the experiment comparison [25].
Until now, various multilayered composite beam theories have been developed by many authors [26][27][28][29][30][31].Several exact mathematical solutions were established based on the concepts of Timoshenko beam theory, including displacementbased theory, first-order shear deformation theory, and higher-order theory [32].Alternatively, Edward and omas [2] developed a plate theory to reflect the mechanical behavior of thin-layer composite structures shown by the model without externally imposed shear force in the cross section of thin layer [2].Based on the plate theory, Hussein has completed analytical and experimental studies that give a calculated simulation focussed on the influence of the adhesive [9].Recently, the plate theory was used in the thin-walled steel to strengthen timber composite structures [33].
e current study for simulating the flexural behavior of the U-shaped member reinforced composite beam mainly adopted the method calculating the strain relation based on the Bernoulli's principle [6,18].Subhani et al. [25] presented an algorithm to be more accurate with a nonlinear model [25].However, both simulations related to the experiment were given by available literature, and this method works by neglecting the transverse shear strains completely defined by the model considering the balance of compressive and tensile forces.
In this paper, the mechanical behavior of the timber beam was reinforced with the U-shaped member as the thinwalled plate made in the steel plate or aluminum alloy strip unlike the composite beam reinforced by two parts fully glued.e method estimating the bending capacity cannot adequate to fit the principle, and so the research was presented to analysis based on thin-layer beam theory.e steel plate or aluminum alloy is widely used in engineering; notwithstanding that FRP has advantages of high tensile strength with easily moldable properties, but the premature catastrophic failure of FPR could easily occur by wood splinters when used in the composite beam [23]. is paper will show the bending capacity of U-shaped thin-walled member reinforced timber beam by the plate theory.e approach will be prescribed and compared with the thicker composite beams [34][35][36].It will be applied to reflect correctly the potential function of U shape in terms of bending capacity and stiffness with respect to the layer slip or shear displacement.

Analytical Models
Several exact mathematical analytical solutions of displacement and stress for composite beams with different combinations have been presented in the professional literature.
In general, two different approaches have been adopted to study composite beams: thick-layer beam theories and thin-layer beam theories.In thicker layer beam theory approaches, researchers wanted to overcome connect issues and explain the behaviors of composite structures more accurately, which is known by the beam being inherently subjected to transverse shear between the sections.Girhammar et al. [36] presented the internal actions and established the first-and second-order models with respect to the slip effect deflection.Nie et al. [37], based on the uniform loading, have proposed the stiffness reduction method in deflection analysis.Girhammar and Gopu [38] presented an exact procedure using two second-order ordinary differential equations in terms of interlayer slip for calculating deflection of composite beams [37,38].In thinlayer beam theory approaches, mainly known as sandwich construction in the form of thin face panels, the analytical model and solution to estimate the effective mechanical properties of adhesives are unavoidable.eoretical studies about the bending capacity subjected to various working loads have been given, and the research on nonrigid adhesive applied the plate theory developed by Hussein [1,2,9,39,40].In fact, the approaches were initiated by the Kirchhoff-Love plate theory.
2.1.Mechanical Model.Due to the highly nonlinear behavior of timber joints, fasteners such as screws are potentially ductile in nature compared with glued joints and carpenter joints.Under compressive stresses, steel or aluminum alloy plate used for reinforced timber can be loaded far over its elastic limit and be able to deform in a distinctly plastic manner, especially the aluminum alloy performs corresponding with timber in terms of mechanism properties behavior.
For a thin-layer beam, the basic structural principle is that the facings act together to resisting moment counteracting the external imposed bending moment.Figure 1 depicts the isolation unit which shows that moment is constant through the thickness, and the core resists the shear stresses set up by the external loads.In Figure 1, q(x) are external loads, P(x) represents the bonding force between two sections, and V i is the value of shear owing to the facings act.en, M i and N i are bending moment and tensile with isolation units.
is kind of theory's model defines a structural sandwich covered by a thin sheet as a major contributor to tensile strength [41].
e isolation unit of thick-layer beam model as two isolated beams connected together by the adhesive is shown in Figure 2. Especially, the moment M 1 is different from M 2 by each part, but in Figure 1, the value of moment M 1 is about equal to M 2 with the thin layer.
For verifying the reliability of the design models of Ushaped plate reinforced timber beam, De la Rosa et al. have analyzed them using the transformed section method and the section conditions of equilibrium.De la Rosa et al. have presented theory based on the formula derived from stressstrain relationship of materials and the conditions of equilibrium with Bernoulli's principle.However, the estimate method was checked with the results obtained experimentally on reinforcement of FPR [18].A different calculation method to simulate the U-shaped plate as reinforcement based on Shen's work was presented, which is translated from the thick-layer beam theory and displayed based on adding the contribution of U-section strength to the element without strengthening, while the reinforcement material is steel plate in concrete beam [42].U section is nonrigid adhesive on the tension face of the timber beam and owed the property of self-supporting due to the unwrapped form and is also weaker in free edge under compressive force caused by the nonbonded adhesive flange of the U section in the lateral sides of the beam.
us, the thin-layer beam may be necessary and sufficient to reflect the mechanical simulation of timber beam assembled by a U-shaped member in such kinds of material properties.
e formulation of the planar two-layer composite beam model uses the following assumptions: (1) the materials are orthotropic and linear elastic; (2) the U-shaped layer is symmetrical with respect to the plane of deformation and forces, and the layers are continuously connected only on the bottom of beam; and (3) the faces are thin in comparison with the core depth.is implies that the exural rigidity of each face about its own middle surface is negligible, and consequently, the inplane stresses resisted by each skin are uniformly distributed across its thickness.
A distribution of U-shaped reinforcement is shown in Figure 3, which is fabricated by attaching timber to the web part of the U-section.e tensile force obviously generates from the connection between the surfaces of two parts.ere are loads which result in bending with the e ect of U-section warping, such as the shear stress between interlayers and the friction both on web and ange between the timbers and due to the shear displacement lead by the distortion that induces in uence of the extent to which the plane sections remain oriented after bending.

Shear Displacement
Analysis.An interlayer slip often develops if it is up to su cient magnitude and therefore a ects signi cantly the deformation and the stress distribution of the composite system.e plate assembled to timber beam having nite sti ness, because of distortion, contributes most to any unequal displacement of the plate and the lower face of the timber with interlayer slip under edgewise loads, as assumed in the model of Figure 4. e U section can not only hold o vertical tensile force P by connectors but also prevent the shear slip by the strengthening function of the ange.
e friction horizontally distributed between plane ange and timber functionally resists partial interaction especially contributed to the force induced by edgewise warping, works together with reinforcement of connector at the bottom of the beam and with extrusion forces generated by the U-section, and plays the role of supplementary stress as well as friction vertically transferred by lateral interfaces.

U-Shape Simulation Based on the in-Layer Beam eory.
e elemental longitudinal section of a prismatic rectangular beam integrates with the U-shaped member under uniformly distributed load along a beam of length L and the total load given by Q. e distance between centroids of principal moment-carrying member h is shown in Figure 3. Considering the interlayer e ects of shear bonding sti ness, de ned as α, the second-order di erential equation [2] for the shear stress τ and de ection ω can be expressed as follows: Figure 1: e isolation unit of the thin-layer beam.

Advances in Civil Engineering
where EI 0 is the sti ness of all beam parts as if unglued and EI ∞ is the sti ness of composite beam if parts were glued together; the bending moment and shear load at a point located at distance x from a reaction are calculated as follows: Substituting the expression for V given in equation ( 4) into equation ( 2), Constants of integration A and B are evaluated using the following boundary conditions; the boundary conditions are evaluated at x 0 and x L/2 , where the length of 0 represents the ends of the beam and half of span represents the midspan of the beam.At points, dτ/dx 0 or τ 0, can be shown by Equation ( 6): Using the above boundary conditions, the following expressions for A and B can be obtained: By substituting equations ( 7) and (8) into equation ( 5), the shear stress is expressed as Due to the strain obtained in equation ( 9), the maximum de ection ω max which under the total compressive Q with integration and evaluation of integration constants leads to where where K is de ned as shear modulus of adhesive or shear load per unit span length to cause unit slip.In a similar manner, an expression for the midspan de ection ful lls other load position conditions that can be obtained, and the function W obtained to represent the constant of integration is given by ) where k is the load position.
Based on the graphs [2,9], and by taking α into numerical evaluation for the parametric study, ω max always increase inducing values of αL/2 decrease when EI ∞ /EI 0 remains increasing proportionally.
A parameter χ which represented constant coe cient under di erent forms of load distribution is used in the various expressions in the interests of brevity.e results can be simpli ed according to equations (10), (16), and (17) with a function ∆κ which re ected reactions between the interfaces of the elements.
e thin-layer beams in equation ( 10) describe the result that can be represented through a function W as e function Δκ is denoted as equation ( 17) due to di erent models in terms of the di erent loading modes: 4

Advances in Civil Engineering
For the first case in equation ( 16), where Δκ is calculated in the term under the offset point load.
It can be seen that parameter α in equations ( 17), (22), and ( 23) significantly determines the function Δκ contributing to midspan deflection.

3.2.
ick-Layer Beam. is method is to simulate the composite beam by second-order shear deformation theory presented [43].It assumes that the displacement u has a cubic variation through the thickness of each section, with moment differential equations of equilibrium for beams considering corrected shear force with respect to slip effects.
e resulting equations for midspan deflection of the beam under uniformly distributed load are given as follows: Likewise, under offset point, the load is given as in all of which and the thick-layer beams in equations ( 18) and ( 19) can be simplified by For the second case in equation ( 21), where Δκ2 is calculated in the term under the uniformly distributed load: In the two cases, the general expressions for displacement are presented, all for a simply supported beam.It can be seen that the algorithms formed in two cases are different shown by the relationships of Δκ with fundamental function χ(QL 3 /EI ∞ ).e thin-layer theory's simplified expression of W contains the Δκ as argument in χ(QL 3 /EI ∞ ) function calls.
e thick-layer theory's expression contains the ∆κ representing the intercept with respect to χ(QL 3 /EI ∞ ).Equation ( 26) can be obtained by where M is the bending moment, N is the axial force, V is the external shear force, i represents each part, and H is the height of the beam.When neglecting external shear force, the method assumes to be applied to simulate the thin-layer model by equations ( 24)-( 26); the results to be calculated are It can be seen that the both sides of equation ( 27), a set of axial interface force, commensurate with extensional stiffness of each parts leading to balancing themselves; that is, the equation obtained a solution that has no value to concern in this way.It means that the equation for two cases could not be interchanged for calculation of composite beam in two forms.Numerical and mechanical analysis are done to illustrate the effectiveness and efficiency of adopting thin layer beam theory for numerical analysis of the U-shaped member composite timber beam without losing any generality, there are insignificant effects to interlayer shears, but the face normal stress presents greater sensitivity to the variation of bond stiffness value as researches [9,17] have shown.In principle, increased bond stiffness induces increased face normal stress.e same function implies that adopting expressions is reasonable due to the absence of loads which were induced by flange assuming an advantage in the remaining composite beam bearing the load together herein.

Finite Element Analysis.
is work was not aimed at identifying the relationship of load and slip deflection for the composite structure but rather to analyze the flexural behavior strengthened by the U-shaped member provided.Here, the load-displacement plot for models of bending capacity simulation with regard to the height of flange was obtained with two types of the composite beam containing two kinds of materials for U-shaped member: one is steel and the other is aluminum alloy.e finite element analysis based on experiments shown in Figure 6 will be discussed further.
e U-shaped member in steel is an SU-OT beam according to data of L-1 [44], while the AU-T beam is an aluminum alloy U-shaped sheet-timber composite beam; the main mechanical properties of aluminum alloy are f t � 206 MPa, E � 0.7 × 10 5 MPa, and c � 0.34; the values of timber according to the test of Pseudotsuga menziesii wood are f c � 28.32 MPa, E � 11003 MPa, and c � 0.3, where E is Advances in Civil Engineering the elastic modulus, c is the Poisson ratio, f t is the yield strength, and f c is the compressive strength.Timber simulation based on the trilinear strain-softening constitutive model [45] and the damaged plasticity model was used to describe the behavior of the aluminum alloy.e nite element analysis carried out the AU-T and SU-OT; it also found the modeling method to be e ective with which the models were calibrated against the experimental load-de ection data.
en, simulations have done to keep the mechanical properties same and change the ange height of AU-T to 0, 5 mm, 13 mm, 26 mm, 39 mm, 52 mm, and 65 mm and 0, 5 mm, 15 mm, 30 mm, and 45 mm for ange height of SU-OT.
It can be seen that the load-displacement response of AU-T beam simulated by ANSYS in Figure 7 suggested that increasing the ange height of U-shaped member has a signi cant in uence on the sti ness as shown by the decrease trend of displacement at midspan beam with the same load condition.e tted curves (e.g., Figures 8 and  9) slightly overestimate both AU-T beam and SU-OT beam, and the reduction in maximum displacement and compressive stress indicates providing higher bearing capacity.
e curves indicate the ange height of the U-shape member adopted 1/5 beam height that is available, and it is obvious that the carrying capacity would not improve when the height is above 1/5 beam height.
Comparing the carbon ber-reinforced timber beam, an increase in its carrying capacity is likely to be produced which proved the good behavior of a U-shaped member with ber-reinforced plastic as well [46][47][48].
e study [18] does not show the relationship of ange height and carrying capacity, but it can be found that reducing the ange height would be superior for carrying capacity according to the equilibrium and the structure analysis sketch.

Experiment of Timber Beams Assembled by U-Shaped Member
To check its e ectiveness, the thin-layer beam calculation model is compared with the results of experiments on a U-shaped thin-walled strip-timber composite beam.SU-OT [44] and AU-T were analyzed by using both the methods described in this paper and the nite element method.From the data, results of two kinds of timber beam are used for  6 Advances in Civil Engineering discussion here.e general view of SU-OT and AU-T is shown in Figure 6.Specimens made of timber with identical characteristics to the beams in terms of timber origin, quality, drying time, and moisture content were tested.In this case, loads and displacements were obtained, and with this information, the corresponding stress/strain diagrams were traced.e values of the elasticity modulus, maximum compression strength, and Poisson's ratio were presented in Section 3.3.e tests were conducted at two points of load application, and the displacement of the beams at midspan and at the supports was measured by dial gauges to an accuracy of 0.01 mm.Investigating the details of behavior of the e ective load-de ection response of the components was carrying out by several full-scale bending tests in experiments; there was relatively ductile behavior exhibited in load-de ection curves, and the simulation showed very similar results compared to test results.With regard to the load bearing capacity observed in the results for identical models, the e ectiveness of parameters used to predict the behavior of beams with di erent heights of U-shaped ange provided an adequate estimate.
e experimental results are shown in Figures 10 and 11, and comparison of displacement at midspan between calculation and test or simulation using thin layer theory is listed in Table 1.e data obtained correspond to the loads applied and to the vertical midspan displacements experienced by the beams where the load is applied, up to fracture.With these data, load-displacement graphs of the beam were obtained as shown in Figure 7, and the graphs of mean values of the strain under ultimate fracture load are shown in Figures 10 and 11. Figure 10 shows the deformation behavior and strain distribution along the beam, in which the values obtained for the U-shaped member ultimate tensile stress corresponded to coinciding with the maximum compressive stress of timber reached by the reinforcing U shape member.When ultimate load is reached, the tensile stress coinciding with the maximum compressive stress would reach the ultimate point.e structure exhibited a ductile behavior.e load-strain relationship in Figure 11 presented the longitudinal strain along the loading procedure.
e results show that there is no occurrence of separation.
It can be seen from Table 1 that the results have shown the ultimate load and corresponding ultimate midspan de ection with calculations by equation ( 12) and the exact value of simulations.e values of ultimate load indicate fracture load, which was adopted from the experimental data in [49], and conclusion by Pellicane [50] can be applied to the material of aluminum alloy when using the algorithm.In this paper, the K value is 0.3482 for SU-OT and 0.4527 for AU-T.e di erences in the results from comparison are less than 5%.e validation is used to check the e ectiveness of the numerical model in calculation and simulation of U-shaped member reinforced timber beam.It still needs various experiments for further research and construction application.

Conclusions
Accordingly, great e ort has been devoted to the development of prototypes of plate strengthening timber hybrid components to be used in composite structure.e U-shaped member combined with timber provides support not only to vertical forces but also to horizontal forces as fractions between interfaces for the failure are triggered by shear slip.Particularly, analysis of timber beams nds that the mild softening part leads the bolts or screws separate before braking by not carrying the vertical load.
Analysis of the mechanical model of interface bond-slip behavior between U-shaped member and timber, and shown in the results of nite element calculation, indicates the Ushape member contributes to the bending capacity, and the available height of ange has been found.e calculation method based on thin-layer beam theory is validated by numerical analysis that can be used to predict the bending capacity of U-shaped sheet reinforced timber.
e results should be treated as qualitative rather than quantitative due to the small range of data.However, the research for the composite system and the U-shaped member reinforced timber beam still needs to be developed.

Figure 2 :Figure 3 :
Figure 2: e isolation unit of the thick-layer beam.

Figure 4 :Figure 5 :
Figure 4: e mechanics model of bond-slip behavior between interface of steel and timber.

Table 1 :
Comparison of displacement at midspan between calculation and test or simulation.