Abstract

With respect to the collapse behavior of joints in a steel pipe truss structure, in the present study, we carry out numerical analyses using a general-purpose FEM software package on the yield strength and the stress distribution near the connection in order to obtain information on the collapse behavior and the yield strength of joints. 𝑇-joint and 𝑌-joint are investigated. It is found that the joint strength is proportional to the 𝜆th power of thickness 𝑇 of a main tube for various joints with a fixed branch-to-chord diameter ratio 𝑑/𝐷. The index 𝜆 changes with 𝑑/𝐷 and in the range of 𝑑/𝐷=0.1 to 0.91, 𝜆1.93 to 1.55 holds with smaller values corresponding to larger 𝑑/𝐷. A compressive axial force along the main tube will reduce the joint strength, and it is also found that the strength-reduction ratio for 𝑌-joint is almost the same as that for 𝑇-joint having the same 𝑇 and 𝑑/𝐷. Also, the joint strength is influenced by the span length of the main tube. The joint strength decreases as the span length increases. However, the influence is small for small 𝑑/𝐷.

1. Introduction

Steel pipe truss structures are widely used in marine structures and large space structures, and the yield strength of joints in a steel pipe truss structure has been extensively studied [15]. A book entitled the “Recommendations for the Design and Fabrication of Tubular Truss Structures in Steel,” published by Architectural Institute of Japan, is available for practical design [1]. However, although there are some studies based on numerical FEM analysis among these past works [6], most studies are purely experimental, and the proposed design equations are based on limited experimental conditions. Therefore, with respect to the collapse behavior of joints in a steel pipe truss structure, in the present study, we carry out numerical analyses using a general-purpose FEM software package on the yield strength and the stress distribution near the connection in order to obtain information on the collapse behavior and the yield strength of joints.

2. Numerical Analysis Method

In this study, assuming an analysis model for the geometry of the 𝑇- and 𝑌-joints as shown in Figure 1, large-deformation elastic-plastic FEM analysis was carried out by using the general-purpose FEM software package, MSC. Marc on the deformation behavior of a joint assuming that compressive load is applied to the branch tube in the axial direction. The boundary condition is such that the compressive load is applied from above by a moving rigid wall in contact with the branch tube, while both ends of the main tube are simply supported.

Furthermore, structural steel pipes are the assumed material, and the stress-strain relation is assumed to obey the bilinear hardening rule. Unless otherwise noted, the material constants are set as follows: yield stress 𝜎𝑌= 360 MPa, Young's modulus 𝐸 = 205 GPa, Poisson's ratio 𝜈= 0.3, and strain-hardening coefficient 𝐸=0.01𝐸. The geometry parameters of the joint are span length 𝐿/𝐷=5, main tube diameter 𝐷=318.5mm, and common thickness of branch and main tubes 𝑡=𝑇=4.5mm.

In the finite element modeling, the size of the element is assumed to be about 20×20mm2 so that the deformation proceeds smoothly, and discretization is performed using 4-node quadrilateral thick shell finite elements. In this analysis, to model nonlinear behavior, the updated Lagrange method is used to consider geometrical nonlinearity, and the algorithm based on the Newton-Raphson method and the return-mapping method is used to solve the nonlinear equation.

3. Analytical Results and Discussions

3.1. Crush Behavior of a 𝑇-joint

Figure 2 shows the relation between the compressive load 𝑃 of a branch tube and displacement 𝑈, as evaluated for a 𝑇-joint with a branch-to-main diameter ratio of 𝑑/𝐷=0.44. The figure shows that load 𝑃 initially increases with compression displacement 𝑈 and shows a peak at a certain level of displacement. As known from the deformation behavior before and after the load peak shown in the figure, we have a compressive load peak because of the denting deformation on the surface of a cylindrical main tube. In addition, Figure 3 shows the contour diagram of von Mises stress distribution near the branch joint, where the brighter color denotes higher stress (the yellow color denotes the highest stress). The schematic view in Figure 3 shows the relation between each contour diagram and joint. Figure 3 also shows the contour diagram of von Mises stress distribution of the main tube at peak load for different branch-to-main diameter ratios (for the ratios of 𝑑/𝐷=0.1 and 𝑑/𝐷=0.26). As the figure shows, if 𝑑/𝐷 is small (e.g., 𝑑/𝐷=0.1), the largest stress appears at both ends of the branch tube (both ends of the branch tube as viewed from the axial direction of the main tube; since they are at the side surface of the main tube, we call them the “side” area). If 𝑑/𝐷 is large (e.g., 𝑑/𝐷=0.44), the largest stress appears before and after the branch tube (before and after the branch tube as viewed from the axial direction of the main tube; since they are at the uppermost parts of the main tube, we call them the “top” area).

Figure 2 also shows the experimental values obtained by other researchers [2]. As the figure shows, the experimental values and the numerical results of FEM analysis exhibit similar trends; therefore, it is confirmed that experimental results can be reproduced by the numerical simulation.

In the following discussion, the collapse load corresponding to the peak in the load-displacement curve shall be defined as the yield strength 𝑃cr of the said joint.

Figure 4 shows the relation between the yield strength of joint 𝑃cr and 𝑑/𝐷, with the thickness 𝑇 of the main tube as a parameter. In this figure, the vertical axis is the ratio of yield strength 𝑃cr and (𝜎𝑌𝑇2): 𝑃cr/𝜎𝑌𝑇2. As shown in the figure, the ratio 𝑃cr/𝜎𝑌𝑇2 is mainly a function of 𝑑/𝐷, increasing with 𝑑/𝐷 but also depends on the thickness 𝑇 of the main tube.

Various empirical formulas have been proposed to date for the yield strength of a 𝑇-joint [13].

In [3] 𝑃cr𝜎𝑌𝑇2𝑑=6.431+4.60𝐷2,(1)

In [2] 𝑃cr𝜎𝑌𝑇2=7.310.833(𝑑/𝐷),(2)

In [1] 𝑃cr𝜎𝑌𝑇2=2.14×𝑓𝛽𝑓𝛾𝑓𝑛.(3)

In (3), 𝑓𝛽𝑑=0.96+4.7𝐷2,𝑓𝛾=𝑑(2𝑇)0.2.(4)

Furthermore, 𝑓𝑛 is a coefficient expressing the effect of the axial force on the main tube, 𝑁, where the following equation holds:𝑓𝑛=1+0.3𝑛0.3𝑛2𝑁,𝑛=𝐴𝜎𝑌,(5) where A is the cross-sectional area of main tube. In Figure 4, 𝑓𝑛=1 because the axial force is 𝑁=0.

Figure 4 also shows the values predicted by the experimental formulas. In many cases, good agreement is found between the experimental formulas and the result of FEM analysis, but there are cases of large error, depending on the thickness 𝑇.

In (1) and (2), 𝑃cr𝑇2 is used for yield strength 𝑃cr, while (3) uses 𝑃cr𝑇1.8. To probe the relation between yield strength 𝑃𝑐𝑟 and main tube thickness 𝑇, the log-log graph of Figure 5 plots yield strength 𝑃cr against main tube thickness 𝑇 for 𝑇-joints with 𝑑/𝐷=0.1, 0.27, 0.44, 0.60, 0.76, and 0.91. As the figure shows, for any constant value of 𝑑/𝐷, the yield strength 𝑃cr is expressed by the exponent function in the form𝑃cr𝑇𝜆,(6) where the exponent 𝜆 depends on 𝑑/𝐷, decreasing against 𝑑/𝐷. In Figure 5, 𝜆1.93 holds for 𝑑/𝐷=0.1 and 𝜆1.55 for 𝑑/𝐷=0.91.

In this study, common thickness of branch and main tubes 𝑡=𝑇=4.5 mm was used. Then, axial compressive buckling for branch tube did not occur. In this case, the yield strength 𝑃cr does not appreciably depend on the branch tube thickness 𝑡.

3.2. Crushing Load Behavior of 𝑌-Joints

The angle between the branch and main tubes is 90 in a 𝑇-joint, while it is expressed by 𝜃 in a 𝑌-joint as shown in Figure 1(b).

As shown in Figure 6, the yield strength 𝑃𝑐𝑟|𝑌 of a 𝑌-joint is greater than that of a 𝑇-joint assuming the same values of 𝑑/𝐷 and 𝑇. To estimate the yield strength of a 𝑌-joint from that of a 𝑇-joint, we normally use the following approximation equation:𝑃cr||𝑌𝑃cr=1sin𝜃.(7)

Figure 7 shows a comparison of the yield strength for a 𝑌-joint predicted by the approximation equation (7) and the values obtained by FEM analysis. The figure shows that the prediction by (7) always lies in the safe side and agrees well with the data obtained by FEM analysis. In particular, in the case of a 𝑌-joint with 𝜃=45, the values estimated by (7) agree well with the values from FEM analysis. However, in a 𝑌-joint with 𝜃<45, the results of FEM analysis vary greatly with 𝑑/𝐷, and the error of using (7) is larger. One reason for this result is thought to be the existence of a force component 𝑃cr|𝑌cos𝜃 in the axial direction of the main tube due to the load 𝑃cr|𝑌. As discussed later in Section 3.3, the yield strength of a joint depends also on the axial force of the main tube, decreasing if the axial force is compressive and increasing if it is tensile. However, as has been pointed out by other researchers [4], in a 𝑌-joint, the force component 𝑃cr|𝑌cos𝜃 in the axial direction of a main tube seems to act differently from the axial force acting on the main tube. Figure 8 shows the von Mises stress distribution at peak load for a 𝑌-joint of 𝑑/𝐷=0.44. (The schematic diagram in the left side of Figure 8 shows the positional relation between the contour diagram and the joint.) Because of the existence of a force component 𝑃cr|𝑌cos𝜃 in the axial direction, we have different stress values at the two top positions, with stress increasing in one and decreasing in the other. Furthermore, the joint geometry of branch and main tubes, which depends on the angle 𝜃 of the 𝑌-joint, also affects the denting deformation, that is, the yield strength of the joint. Therefore, it is difficult to predict accurately the yield strength of a 𝑌-joint from the yield strength of a 𝑇-joint. Note, however, that if 𝑑/𝐷 is small, the effect of a force component 𝑃cr|𝑌cos𝜃 of an inclined load 𝑃cr|𝑌 in the axial direction of the main tube is considered small because the denting part of the main tube is mainly subject to the stress in the side part of the main tube as shown in Figure 3.

3.3. Effect of Axial Force along Main Tube

There is an axial force along the main tube in the actual structural object. Although the effect of the axial force along the main tube on the yield strength of a joint has been studied and design formulas have been proposed [1], the number of reports is few because the condition setting by experiment is difficult. Therefore, we examined a wide range of conditions through numerical analyses.

Figures 9 and 10 show the ratio of yield strength 𝑃cr/𝑃cr|𝑁=0 for a 𝑇-joint and the ratio 𝑃cr|𝑌/𝑃cr|𝑌,𝑁=0 for a 𝑌-joint with and without axial force along the main tube. As the figures show, the yield strength decreases if the axial force is compressive and increases if the axial force is tensile. The decrease and increase rates relate to the diameter ratio 𝑑/𝐷 between branch and main tubes of the joint, increasing with 𝑑/𝐷. Furthermore, from Figures 9 and 10, we can assume that the change rate of yield strength 𝑃cr|𝑌/𝑃cr|𝑌,𝑁=0 at a 𝑌-joint and the change rate of yield strength 𝑃cr/𝑃cr|𝑁=0 at a 𝑇-joint are approximately the same for the same axial force and the same 𝑑/𝐷, and the effect of the axial force along the main tube on yield strength 𝑃cr is of similar magnitude in 𝑇- and 𝑌-joints.

Figure 11 shows the relation between compressive load 𝑃 and branch tube displacement 𝑈 obtained from the numerical FEM analysis for a 𝑇-joint with and without axial force along the main tube; it is found that in comparison with the case without axial force along the main tube, the peak of compressive load 𝑃 in the branch tube appears later if the axial force along the main tube is tensile. Figure 12 shows the variation of the stress 𝜎𝑥 in the axial direction at the top part of the main tube, owing to the indentation by the branch tube, as evaluated by the numerical FEM analysis. We focus on the top part of the main tube because a large stress appears in this part, as shown earlier in Figure 3, and this stress governs the yield strength 𝑃cr of the joint. As Figure 12 shows, if the axial force is zero, the stress 𝜎𝑥 in the direction of axial force, appearing at the top part due to the load 𝑃, is compressive, of which the absolute value increases with deformation and reaches yield stress as the compressive load reaches a peak (mark A in the figure). Furthermore, if the axial force is compressive, the compressive stress due to the axial force in the main tube is added to the compressive stress in the top part; the yield strength becomes small because the stress is high from the very early stage of deformation. On the other hand, if the axial force is tensile, the whole body of main tube is initially in the tensile state; if load is applied to the branch tube, compressive stress appears in the top part but is canceled by the initial tensile stress; a larger amount of load is needed for the deformation to reach the peak corresponding to yield strength 𝑃cr, and it would be natural to assume that the yield strength 𝑃cr increases.

Equation (5) is proposed in [1] to predict the effect of axial force, and Figures 9 and 10 also show the values obtained from (5). It is found that the result of numerical analysis, which shows the same trend as (5), is always on the danger side and the difference is great especially if 𝑑/𝐷 is large. For example, in a 𝑇-joint of 𝑑/𝐷=0.76, the ratio between the yield strength obtained by FEM analysis and the yield strength with zero axial force 𝑃cr/𝑃cr|𝑁=0 is 0.600 if the axial force is 𝑁/(𝐴𝜎𝑌)=0.6 while (5) gives 0.712 with the error 19% on the danger side. Therefore, Figures 9 and 10 also show the calculation using a modified formula:𝑓𝑛=1+0.35𝑛0.35𝑛2𝑁𝑛=𝐴𝜎𝑌.(8) As the figure shows, the value calculated using (8) is nearly equal to the change rate of yield strength 𝑃cr/𝑃cr|𝑁=0 due to the compressive axial force for a joint of 𝑑/𝐷=0.44.

3.4. Effect of Span Length of Main Tube

Figure 13 shows the variation of yield strength 𝑃cr against the variation of span length 𝐿, which is obtained for the 𝑇-joints with various values of 𝑑/𝐷 taking 𝐿0=1592.5mm as a standard span length. Most approximation equations used to date do not specifically take into account the effect of span length 𝐿. However, the yield strength of a joint decreases with increasing span length, as shown in Figure 13. Although the effect is small if 𝑑/𝐷 is small, the effect of span length 𝐿 on the yield strength of a joint is measurably great if 𝑑/𝐷 is large. For example, against the joint of 𝑑/𝐷=0.91, the yield strength decreases to about 73% for the joint of 𝐿/𝐿0=2.

The effect of span length on the yield strength of a joint can be discussed considering the bending stress in the main tube due to the load 𝑃 of the branch tube. First, it would be qualitatively understood that the effect of span is small when 𝑑/𝐷 is small if one considers that the dent of the main tube is spatially limited to the side or the top area, respectively, for small or large 𝑑/𝐷 and that the bending stress of the main tube is the greatest at the top area. Then, in the following, the effect of span length is quantitatively examined from the viewpoint of bending stress in the main tube due to the load 𝑃.

Here, it is assumed that yield strength is reached in the joint when the stress 𝜎𝑥 in the axial direction at the top area of the main tube reaches a certain critical value 𝜎cr of stress. Furthermore, it is assumed that the stress in the axial direction 𝜎𝑥 is the sum of local stress 𝜎loc due to the indentation caused by the branch tube and the bending stress 𝜎glo at the top part of the main tube due to the load 𝑃:𝜎cr=𝜎loc+𝜎glo.(9)

Since the stress 𝜎loc is not related to span length and is proportional to the yield strength 𝑃cr, it can be evaluated by using𝜎loc=𝑘1𝑃𝐷𝑇cr.(10)

Meanwhile, since the stress 𝜎glo is bending stress due to the load 𝑃, it is proportional to span length 𝐿 and can be evaluated by using𝜎glo=𝑘2𝐷2𝑇𝑃cr𝐿.(11) In (10) and (11), 𝑘1 and 𝑘2 are unknown coefficients.

Using (9)–(11), we can estimate the yield strength for any span from the yield strength for the two different spans. For example, assuming the yield strength 𝑃cr|𝑎 and 𝑃cr|𝑏 for the joints with span lengths 𝐿𝑎 and 𝐿𝑏, the yield strength 𝑃cr for a joint with an arbitrary span length 𝐿 is estimated as follows. From (9) to (11), the following equations hold for the joints with the span lengths 𝐿𝑎, 𝐿𝑏, and 𝐿:𝜎cr=𝑘1𝑃𝐷𝑇cr||||𝑎+𝑘2𝐷2𝑇𝑃cr||||𝑎𝐿𝑎𝜎cr=𝑘1𝑃𝐷𝑇cr||||𝑏+𝑘2𝐷2𝑇𝑃cr||||𝑏𝐿𝑏𝜎cr=𝑘1𝑃𝐷𝑇cr+𝑘2𝐷2𝑇𝑃cr𝐿.(12)

From (12), from 𝑃cr|𝑎 and 𝑃cr|𝑏, we have an equation to estimate the yield strength 𝑃cr of a joint of length 𝐿 as follows:𝑃cr=𝑃cr||𝑎𝑃cr||𝑏𝐿𝑏𝐿𝑎𝑃cr||𝑎𝑃cr𝐿𝐿𝑎+𝑃cr||𝑏𝐿𝑏𝐿.(13) Assuming 𝐿𝑎=5𝐷 and 𝐿𝑏=10D, Figure 14 shows a comparison of the approximation by (13) and the numerical FEM analysis for the yield strength 𝑃cr of joints of length 𝐿=2.5D to 12.5D. The figure shows that the two mostly agree with each other for a long span that satisfies, for example, 𝐿/𝐷>4.

We can further examine the effect of span on yield strength through the effect of axial force. That is, we assume that the effect of the yield strength due to the bending stress 𝜎glo in (11) is analogous to the effect of the axial force of the main tube; the ratio of axial force 𝑛 as used in (8) can be evaluated by using𝑘𝑛=0𝜎glo𝜎glo||0𝜎𝑌.(14) Here, 𝜎glo|0 is the bending stress at the standard span length 𝐿0. By substituting (11) into (14), we obtain𝑘𝑛=0𝑘2𝐷2𝑇𝑃cr𝐿𝑃cr||0𝐿0𝜎𝑌.(15) Here, 𝑃cr|0 is the yield strength of joint corresponding to the standard span length 𝐿0. Furthermore, the coefficient 𝑘0𝑘2 is assumed to be approximately 𝑘0𝑘21/1.5 based on investigating the application of the equation to various span lengths.

Therefore, using (8), we have an equation for 𝑃cr as follows:𝑃cr𝑃cr||0𝑃=10.35cr𝐿𝑃cr||0𝐿01.5𝐷2𝑇𝜎𝑌𝑃0.35cr𝐿𝑃cr||0𝐿01.5𝐷2𝑇𝜎𝑌2.(16) By solving this quadratic equation, we can estimate the yield strength 𝑃cr of a joint with a span length 𝐿 from the yield strength 𝑃cr|0 for the standard length 𝐿0.

Figure 15 shows a comparison of the approximation by (16) and the numerical FEM analysis for the yield strength of a joint with length 𝐿=2.5𝐷 to 12.5𝐷 assuming 𝐿0=5𝐷. It is shown that the two agree roughly with each other for any long span satisfying, for example, 𝐿/𝐷>4.

Meanwhile, although the span length clearly influences yield strength, the issue of how to define the standard span length for design of an actual structure and apply it in the design equations remains. These issues are left for future study.

4. Conclusions

Although the yield strength of steel joints depends not only on buckling and plastic yield but also on the brittle failure of welded parts and punching shear, we focused on only the yield strength originating from buckling and plastic yield and performed analysis with a general-purpose analysis software using the finite element method on the yield strength of 𝑇- and 𝑌-joints of steel pipes in the range of diameter ratio 𝑑/𝐷=0.1 to 0.91 and span length 𝐿/𝐷=2.5 to 12.5. As a result, the following conclusions were obtained.(1)When the diameter ratio 𝑑/𝐷 is kept constant between branch and main tubes, the yield strength 𝑃cr depends on the power function of the thickness 𝑇 of the main tube, 𝑃cr𝑇𝜆, where the power 𝜆 depends on 𝑑/𝐷, decreasing against 𝑑/𝐷. 𝜆1.931.55 holds in the range 𝑑/𝐷=0.10.91.(2)The yield strength of a joint decreases as span length increases. Furthermore, although the effect is small if 𝑑/𝐷 is small, the effect of span length 𝐿 on the yield strength of a joint is large if 𝑑/𝐷 is large. We can examine the effect of span length on the yield strength of a joint by considering the bending stress of the main tube due to the load 𝑃.(3)The yield strength of a joint decreases due to the axial compressive force of the main tube, and the magnitudes of the effect are similar in 𝑇- and 𝑌-joints. The effect of axial force can be understood considering that the dent of the main tube is governed mainly by the axial stress 𝜎𝑥 due to the load 𝑃 and the axial stress 𝜎𝑥 is compressive.(4)Although the relation 𝑃cr|𝑌<𝑃cr/sin𝜃 always holds between the yield strength 𝑃cr|𝑌 of a 𝑌-joint and that of a 𝑇-joint for the same value of 𝑑/𝐷, the equation 𝑃cr|𝑌𝑃cr/sin𝜃 approximately holds. In particular, the error is small at 𝜃=45 or when 𝑑/𝐷 is small.