Next Article in Journal
Mesenchymal Stromal/Stem Cells Isolated by Explant Culture Method from Wharton’s Jelly and Subamnion Possess Similar Biological Characteristics
Next Article in Special Issue
Investigations of the Windage Losses of a High-Speed Shrouded Gear via the Lattice Boltzmann Method
Previous Article in Journal
Research on Intraparticle to Interparticle Entanglement Swapping Protocols
Previous Article in Special Issue
Numerical Investigation of Low-Velocity Ice Impact on a Composite Ship Hull Using an FEM/SPH Formulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Yield Behaviour of Welded I-Shaped Steel Cross-Sections

by
Luigi Palizzolo
and
Salvatore Benfratello
*
Department of Engineering, University of Palermo, Viale delle Scienze, I-90128 Palermo, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(17), 8037; https://doi.org/10.3390/app14178037
Submission received: 5 August 2024 / Revised: 27 August 2024 / Accepted: 6 September 2024 / Published: 8 September 2024
(This article belongs to the Special Issue Mathematical Methods and Simulations in Mechanics and Engineering)

Abstract

:
The limit behaviour of I-shaped welded steel cross-sections subjected to axial force, shear, and bending moment is a crucial matter to ascertain the reliability of framed structures constituted by non-standard beam elements. International standards provide an approximate solution to the problem, and other studies have proposed improved approximate formulations to ascertain the real features of the relevant cross-sections. The present paper is devoted to enhancing the problem of the limit behaviour of plane I-shaped welded steel cross-sections subjected to axial force N , shear T and bending moment M ; therefore, new appropriate formulations are proposed in order to define suitable new domains, both in planes N , T , N , M , and M , T and in the space N , T , M . The material is assumed as elastic–perfectly plastic and the Von Mises limit condition is adopted as the resistance criterion. The elastic stresses are described by the Navier formula and the Jourawski formula. The limit stress condition related to the contemporaneous presence of the acting forces is defined as the one that, at each point of the cross-section, fulfils the Von Mises limit condition as equality. The formulation is rigorously devoted to factory-made welded I-shaped steel cross-sections. Some numerical examples are reported in the application stage and useful comparison are carried out, with the results being obtainable by the application of the classical known standard formulae, proving the reliability and effectiveness of the determined domains.

1. Introduction

In many cases of practical interest, in the design of steel frame structures, it is very useful to adopt factory-welded I-shaped beams instead of the standard rolled ones. One of these cases is that of seismic retrofitting of existing masonry structures. In such a case, to obtain the required seismic behaviour, it is necessary to modify the stiffness of the masonry panels. For this goal, steel hooping is considered the most adequate, and the adoption of suitably designed and factory-welded steel elements allows us to obtain optimal structures with reduced size and, consequently, with little effect on usability. Another task is that of new openings in masonry panels, which requires the steel hooping to restore both strength and stiffness of those panels where the opening is created. In such cases, usually, the required stiffness is often very high and, consequently, the chance of adopting steel frames realized with suitably designed and factory-welded I-shaped beam elements with assigned geometric characteristics is a fundamental strategy which also allows for a reduction in the opening dimensions.
On the other hand, the above referenced chance of adopting welded structural elements is also often advantageous in new frame design. As is well known in such structures, the connections between beams and columns are very critical, since their design plays a crucial role in the safety and ductility of the overall structure. This topic is relevant when elastoplastic analysis is performed both statically [1,2,3,4] or dynamically [5,6,7,8,9], as well as in design problems for frame and truss structures [10,11,12], both taking into account buckling and slenderness [13,14,15,16] and in the case of seismic and dynamic actions [17,18,19,20].
In practical applications, an important strategy is to plan interventions able to reduce the strength of the element in suitably selected portions of the element close to the connections between beams and columns. In this way, the onset of plastic strain also allows for the protection of the welded connections between beams and columns. The most common approach belonging to this strategy is the so-called dogbone, while a recently proposed approach is referred to as the Limited Resistance Plastic Device (LRPD). Another approach is based on dissipative approaches leading, e.g., to devices for seismic-resistant steel frames developed by the FUSEIS EC project (RFSR-CT-2008-00032) or other devices [21,22,23,24].
The main idea of the dogbone [25,26,27,28,29] is that of suitably reducing the base of the I-shaped beam in a portion close to the column to obtain a reduction in the limit bending moment. Such a reduction in the case of high loads causes the required yielding and decreases in the actions on the welded connections between beams and columns. The seismic performances of structures equipped with dogbones is the topic of many papers available in the literature (see, e.g., [30,31,32]), and developments of the main idea (i.e., reducing the web) have been also proposed (see, e.g., [33,34,35,36]).
In recent years, a number of innovative systems based on energy dissipation and damping have been invented as a result of national and European research projects, such as dissipative connections assigned to braced steel frames (developed in “Two INnovations for Earthquake Resistant Design”, INERD project) [24] and dissipative links (see, e.g., [37,38]). These innovative systems consist of small and dismountable dissipative parts where seismic damage is concentrated; hence, they have increased repairability while displaying comparable stiffness and ductility to the conventional ones. The structural behaviour of structures equipped with these devices is assessed in many papers available in the literature (see, e.g., [39]).
The LRPDs represent the most recent proposal in this topic; these innovative devices can be positioned either at the ends of beams and at the bases of columns, and their peculiar design not only incites the onset of plastic strains on the selected portions of the equipped beam elements, but, furthermore, also guarantees the invariance of the bending stiffness of the overall involved elements.
The first version of LRPD was referred to as limited-resistance rigid perfectly plastic flexural hinges (LRPH). Later, the design problem was enriched by taking into account the simultaneous presence of both normal and tangential stresses and adopting a special analysis of the beam elements by means of smart beam elements specialized to the design of a frame with assigned capacity curve [40], as well as performing both numerical and experimental validation. Other formulations, among which one is approximated to be useful in practical applications, have been proposed to compare the results with those achieved by different typologies of reduced beam section (RBS) connections [41]; some new formulations have also been proposed to consider the local buckling and to protect the welded connections between beams and columns [42].
The actual version, referred to as LRPD, is a beam element with multistep geometry, constituted by three subsequent I-shaped portions whose web and flange thicknesses are suitably designed.
In this last version of LRPD, the inner portion with reduced dimensions shows limited strength and allows the onset of plastic strain, while the outer portions are equal to each other and show greater web and flange thicknesses to guarantee that the bending stiffness of the overall involved element is unaltered. As can be easily deduced, the LRPD design requires the solution to appropriate optimal design problems, and their production has to be performed in a factory by welding steel plates of appropriate thickness.
Expanding the focus to welded I-shaped sections of any size and referring to beam elements belonging to frame structures, these sections are subjected to the simultaneous action of axial force N , shear force T , and bending moment M . It is evident that the knowledge of both the elastic and plastic limit domains of such sections is essential to evaluate the overall structural safety. The actual international standards [43] provide approximated expressions to identify such domains, and, therefore, it is appropriate to examine this topic in depth to obtain more rigorous expressions for practical applications. In the literature, some papers have recently been proposed for either elastic [41] or plastic [44], domains but the formulations for plastic domains still present approximations, and they involve standard sections.
Therefore, the aim of the present paper is to propose a rigorous procedure to identify the plastic domains of welded I-shaped sections subjected to the simultaneous action of axial force N , shear force T , and bending moment M , assuming an elastic–perfectly plastic material behaviour and imposing that the section belongs to Class 1 of the EC3 standard [43]. The expressions of the plastic domain contours in the ( N , T ,   M ) space, as well as those in the ( N , M ) , ( N , T ) , and ( M , T ) planes, will be provided. Some numerical applications for the case of welded sections and a comparison with the results from the EC3 standard conclude the paper.

2. Materials and Methods

The mechanical behaviour of the I-shaped welded section sketched in Figure 1, subjected to the simultaneous presence of axial force N , shear force T z , and bending moment M y , is presented. An elastic–perfectly plastic behaviour of the material is considered, and the adopted yielding criterion is the Von Mises one:
σ 0 2 = σ x 2 + 3 τ x y 2 + τ x z 2
σ 0 is the yielding stress of the material. Due to the symmetry of the domains to be determined with respect to the N , M , N , T , T , M planes as well as in the space N , T , M , in the following section, reference will be made only to positive quadrants or octants.
Referring to the geometrical characteristics of the section reported in Appendix A, the following limit internal forces can be determined:
N E =   N p = A σ 0
M E = W E y σ 0
M p = W p y σ 0
T E = σ 0 3 I y t w S y , G
T p = σ 0 3 t w h 2 t f + 2 σ 0 t w 3 · S y h 2 t f h 2 t f h 2 S y z d z
N E , N p , M E , M p , T E , and T p are the elastic and plastic limit axial force, bending moment, and shear force, respectively. In Equation (5), τ 0 = σ 0 / 3 is assumed. In Equation (6) T p is derived referring to the limit condition sketched in Figure 2, that is, the web is assumed to be fully plasticized and the presence of τ x z in the flanges is taken into account.
The following final remarks can be made: when N = N p , then T = 0 and M = 0 result; when M = M p , then N = 0 and T = 0 result; when T =   T p , then N     0 and M     0 can be.

2.1. Simultaneous Presence of T = T p and N   M = 0

The presence of T p as the acting shear force does not exclude the simultaneous presence of a limited axial force. In such a case, all along the web, τ x z = τ 0 =   σ 0 / 3 occurs, so the presence of normal stresses is excluded. Instead, in the flanges, both τ x y and τ x z are present, assuming values lesser than τ 0 , so that the presence of normal stresses is possible. Referring to Figure 3a, in the ranges h 2 t f     z     h 2 and t w 2     y   b 2   , the expression of τ x y is given by the following relation:
τ x y y = T p b 2 y h t f 2 I y
The expression of τ x z , in the ranges h 2 t f z h 2 and 0 y b 2 (Figure 3b), is given by:
τ x z z = T p h 2 4 z 2 2 I y
By imposing the full yielding condition on all the points of the cross-section in the ranges h 2 t f z h 2 and t w 2 y b 2 , the normal stresses acting in the flange must respect Equation (1) as equality; that is, taking into account Equations (1), (7), and (8), the following relation is obtained:
σ x y , z = σ 0 2 3 T p 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2
By imposing the full yielding condition on all the points of the cross-section where τ x y = 0 (i.e., in the ranges h 2 t f z h 2 and 0   y t w 2 ), the normal stresses acting in the flange must respect Equation (1) as equality; that is, taking into account Equations (1) and (8), the following relation is obtained:
σ x z = σ 0 2 3 T p 2 h 2 4 z 2 2 4 I y 2
It follows that the limit axial force acting when T = T p is given by:
N l i m T p = 4 t w 2 b 2 h 2 t f h 2 σ 0 2 3 T p 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 d z d y + 2 t w h 2 t f h 2 σ 0 2 3 T 2 h 2 4 z 2 2 4 I y 2 d z
The pair of values N l i m T p , T p represents the internal force condition, which determines the full yielding of the whole cross-section. It is worth noting that all the pairs of values N ,   T p , with 0 N < N l i m T p , always represent limit conditions, but only for shear; indeed, in this latter condition, the shear force reaches its limit yield value, but the web suffers plastic strains due to the shear, while the flanges behave elastically. It follows that the segment 0 , T p ; N l i m T p , T p in the plane N , T represents the upper plateau of the relevant yield domain boundary. The entire yield domain boundary is defined in Section 2.4 and is represented in the application stage for two representative cross-sections.

2.2. Simultaneous Presence of T = T p and M   N = 0

Analogously, on the grounds of Equations (7)–(10), the limit bending moment acting when T = T p is given by:
M l i m T p = 4 t w 2 b 2 h 2 t f h 2 σ 0 2 3 T p 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 z d z d y + 2 t w h 2 t f h 2 σ 0 2 3 T 2 h 2 4 z 2 2 4 I y 2 z d z
As previously noted for the axial force, the pair of values M l i m T p , T p represents the internal force condition, which determines the full yielding of the whole cross-section. It is worth noting that all the pairs of values M , T p , with 0 M < M l i m T p , always represent limit conditions, but only for shear; indeed, in the latter condition, the shear force reaches its limit yield value, but the web suffers plastic strains due to the shear, while the flanges behave elastically. It follows that the segment 0 , T p ; M l i m T p , T p in the plane M , T represents the upper plateau of the relevant yield domain boundary. The entire yield domain boundary is defined in Section 2.5 and represented in the application stage for two representative cross-sections.

2.3. Simultaneous Presence of T = T p   N and M

A successive step is considering the simultaneous presence of T p as the acting shear force together with both the axial force N < N l i m T p and bending moment M < M l i m T p .
With the aim of taking into account the influence of both the axial force and bending moment, reference can be made to the usual procedure utilized for the definition of the yield domain in the plane N , M . Therefore, each flange is divided in two parts by a segment parallel to the y -axis; the inner ones are subjected to normal stress related to the axial force, while the outer ones to normal stresses related to the bending moment. Due to these remarks, the distance z = z N dividing the above-described parts lies in the range h 2 t f z N h 2 . The infinite pairs of N , M values which determine the yielding conditions for the whole cross-section are given by the following relations:
N T p = 4 t w 2 b 2 h 2 t f z N σ 0 2 3 T p 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 d z d y + 2 t w h 2 t f z N σ 0 2 3 T p 2 h 2 4 z 2 2 4 I y 2 d z
M T p = 4 t w 2 b 2 z N h 2 σ 0 2 3 T p 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 z d z d y + 2 t w z N h 2 σ 0 2 3 T p 2 h 2 4 z 2 2 4 I y 2 z d z
It is worth noting that the deduced values define a plane boundary for T = T p ; the surface within this boundary is the locus of the points for which the yield condition is reached referring just the shear force, while the flanges still behave in an elastic condition.

2.4. Yield Domain in N , T Plane ( M = 0 )

The next step is to define the boundary of the yield domain N , T . To this goal the range 0 T T P is divided in the following subranges 0 T T E and T E T T p .
In the range 0 T T E , for an assigned value of the shear force T , the expression of τ x y , for h 2 t f z h 2 and t w 2 y b 2 , is given by:
τ x y y = T b 2 y h t f 2 I y
In the flange identified by the ranges h 2 t f z h 2 and 0 y b 2 , the expression of τ x z is given by:
τ x z z = T h 2 4 z 2 2 I y
In the web, for the range 0 z h 2 t f , τ x z is given by:
τ x z z = T b t f h t f + t w h 2 t f 2 z 2 2 I y t w
By imposing the Von Mises full yielding condition in all points of the cross-section where both tangential stresses are present, the following expression for σ x y , z is obtained:
σ x y , z = σ 0 2 3 T 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2
By performing analogous steps for the flange portion where only τ x z z acts, the following relation is obtained:
σ x z = σ 0 2 3 T 2 h 2 4 z 2 2 4 I y 2
By performing analogous steps for the web where only τ x z z acts, the following relation is obtained:
σ x z = σ x z = σ 0 2 3 T 2 b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2
Finally, the limit axial force N l i m T is given by
N l i m T = 4 t w 2 b 2 h 2 t f h 2 σ 0 2 3 T 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 d z d y + 2 t w h 2 t f h 2 σ 0 2 3 T 2 h 2 4 z 2 2 4 I y 2 d z + 2 t w 0 h 2 t f σ 0 2 3 T 2 b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 d z
In the range T E T T P a portion of the web is fully plasticized with τ x z = τ 0 = σ 0 / 3 due to the acting shear. In the other portions of the web, τ x z shows the usual parabolic expression, while in the flanges, τ x y acts with a linear expression.
Let the fully plasticized portion of the web be identified by a value of z = z ¯ in the range of 0 ,   h 2 t f . Clearly, in this portion:
τ x z z ¯ = σ 0 3
In the remaining portion of the web, i.e., for z ¯ z h 2 t f ,
τ x z z = σ 0 3 S y z S y z ¯
In the flanges, for h 2 t f z h 2 ,
τ x z z = σ 0 3 S y z S y z ¯ t w b
It follows that, for an assigned 0 z ¯ h 2 t f , the shear force is given by:
T z ¯ = σ 0 3 2 t w z ¯ + 2 z ¯ h 2 t f σ 0 3 S y z S y z ¯ t w d z + 2 h 2 t f h 2 σ 0 3 S y z S y z ¯ t w d z
In the web portion outside the fully plasticized one, the acting normal stress must fulfil the Von Mises yielding condition, leading to the following relation:
σ x z = σ x z = σ 0 1 S y z S y z ¯ 2
and, as a consequence, the axial force N l i m T z ¯ reads
N l i m T z ¯ = 4 t w 2 b 2 h 2 t f h 2 σ 0 2 3 T 2 z ¯ 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 d z d y + 2 t w h 2 t f h 2 σ 0 2 3 T 2 z ¯ h 2 4 z 2 2 4 I y 2 d z + 2 t w z ¯ h 2 t f σ 0 2 3 T 2 z ¯ b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 d z

2.5. Yield Domain in M , T Plane ( N = 0 )

The boundary of the yield domain in the plane M , T can be defined by performing analogous steps to those in the foregoing section. Therefore, the following relations are obtained.
In the range 0 T T E , M l i m T is given by:
M l i m T = 4 t w 2 b 2 h 2 t f h 2 σ 0 2 3 T 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 z d z d y + 2 t w h 2 t f h 2 σ 0 2 3 T 2 h 2 4 z 2 2 4 I y 2 z d z + 2 t w 0 h 2 t f σ 0 2 3 T 2 b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 z d z
In the range T E T T p , M l i m T is given by:
M l i m T z ¯ = 4 t w 2 b 2 h 2 t f h 2 σ 0 2 3 T 2 z ¯ 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 z d z d y + 2 t w h 2 t f h 2 σ 0 2 3 T 2 z ¯ h 2 4 z 2 2 4 I y 2 z d z + 2 t w z ¯ h 2 t f σ 0 σ 0 2 3 T 2 z ¯ b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 z d z

2.6. Yield Domain in N , T , M Space

For T = T p , the yield conditions for N and M have been already determined in the foregoing sections by Equations (13) and (14). In the following section, the relations characterizing the 3D domain N , T , M will be determined separately by considering the two subsequent ranges: 0 T T E and T E T T p .
In the range 0 T T E , the 3D domain N , T , M is obtained by assigning a set of values to T and determining the consequent N , M plane domain boundaries. To this aim, as usual, increasing values to z N are assigned, with 0 z N h 2 , 2 z N being the portion of the cross-section area (symmetric with respect to the y -axis) that reaches the full plasticization as a result of the presence of shear force and axial force. For the typical value of T , in the flanges, τ x y shows the expression reported in Equation (15), while τ x z shows the expression reported in Equation (16). In the web, the expression of τ x z shows the expression reported in Equation (17). It is important to note that, due to the geometric discontinuities characterizing the shape of the section under consideration for z N = h 2 t f , the relations which define the limit values for N and M are defined for each range: 0 z N h 2 t f and h 2 t f z N h 2 .
In the range 0 z N h 2 t f , the normal stress acting in the range 0 z z N is given by Equation (20). It follows that the axial force N z N is given by:
N z N = 2 t w 0 z N σ 0 2 3 T 2 b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 d z
and the bending moment M z N is given by:
M z N = 4 t w 2 b 2 h 2 t f h 2 σ 0 2 3 T 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 z d z d y + 2 t w h 2 t f h 2 σ 0 2 3 T 2 h 2 4 z 2 2 4 I y 2 z d z + 2 t w z N h 2 t f σ 0 2 3 T 2 b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 z d z
In the range h 2 t f z N h 2 , the axial force N z N is given by:
N z N = 2 t w 0 h 2 t f σ 0 2 3 T 2 b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 d z + 4 t w 2 b 2 h 2 t f z N σ 0 2 3 T 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 d z d y + 2 t w h 2 t f z N σ 0 2 3 T 2 h 2 4 z 2 2 4 I y 2 d z
The bending moment M z N is given by:
M z N = 4 t w 2 b 2 z N h 2 σ 0 2 3 T 2 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 z d z d y + 2 t w z N h 2 σ 0 2 3 T 2 h 2 4 z 2 2 4 I y 2 z d z
In the range T E T T p , for 0 z ¯ h 2 t f the shear force is given by
T z ¯ = σ 0 3 2 t w z ¯ + z ¯ h 2 t f S y z S y z ¯ d z + h 2 t f h 2 S y z S y z ¯ d z
Once the shear force is given, the corresponding N , M domain boundary is obtained as described in the foregoing sections. By selecting different values of z N in the range of z ¯ z N h 2 and identifying the portion of the web where the yield condition is reached due to the presence of normal stresses due to the axial force and shear stresses, the relevant pairs of axial force and bending moment defining the searched boundary can be obtained by considering the two subranges z ¯ z N h 2 t f and h 2 t f z N h 2   s e p a r a t e l y .
In the range z ¯ z N h 2 t f , the axial force N T z ¯ , z N and the bending moment M T z ¯ , z N are given by:
N T z ¯ , z N = 2 t w z ¯ z N σ 0 2 3 T 2 z ¯ b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 d z
M T z ¯ , z N = 4 t w 2 b 2 h 2 t f h 2 σ 0 2 3 T 2 z ¯ 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 z d z d y + 2 t w h 2 t f h 2 σ 0 2 3 T 2 z ¯ h 2 4 z 2 2 4 I y 2 z d z + 2 t w z N h 2 t f σ 0 2 3 T 2 z ¯ b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 z d z
In the range h 2 t f z N h 2 , the axial force N T z ¯ , z N and the bending moment M T z ¯ , z N are given by:
N T z ¯ , z N = 4 t w 2 b 2 h 2 t f z N σ 0 2 3 T 2 z ¯ 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 d z d y + 2 t w h 2 t f z N σ 0 2 3 T 2 z ¯ h 2 4 z 2 2 4 I y 2 d z + 2 t w z ¯ h 2 t f σ 0 2 3 T 2 z ¯ b t f h t f + t w h 2 t f 2 z 2 2 4 I y 2 t w 2 d z
M T z ¯ , z N = 4 t w 2 b 2 z N h 2 σ 0 2 3 T 2 z ¯ 4 I y 2 b 2 y 2 h t f 2 + h 2 4 z 2 2 z d z d y + 2 t w z N h 2 σ 0 2 3 T 2 z ¯ h 2 4 z 2 2 4 I y 2 z d z

3. Numerical Application

In order to test the reliability of the domain formulations defined in the foregoing section, in the following section, numerical applications to two welded sections are presented and compared with the domain defined in the EC3 standard. The selected sections possess the geometrical characteristics reported in Table 1. As can be easily observed, with the aim of performing the necessary comparison, the geometry of section W1 is chosen with width, depth, and thickness equal to those of HEA300, while the geometry of section W2 is chosen with width, depth, and thickness equal to those of IPE360.
The equations defining the proposed domains were written in Mathematica® 14.0 software to obtain the corresponding graphs. For the plain domains, the investigated range for z was divided into 20 intervals, leading to satisfactory results, which are reported in the following section. For 3D domains, the investigated range for T was divided into 10 intervals, and for each of these intervals, 100 points were calculated. The overall computing time for each section was equal to 5 min on HP I7-32 GB RAM workstation.
As the first application, the case of the W1 section was considered. In Figure 4, the N , M domain obtained by Equations (13) and (14), the N , T domain obtained by Equations (21) and (27), and the M , T domain obtained by Equations (28) and (29) are sketched, respectively.
In Figure 5, the N , T , M domain is sketched both with and without a bounding surface.
In Figure 5a, the plane domains reported in Figure 4 are also reported to enhance the overall coherence of the proposed approach.
It is worth noting that the boundary of the upper plateau of the domain in Figure 5b represents the yield conditions for the whole cross-section, while all the points within the relevant boundary define the limit conditions related to the shear forces; in other words, for this plane domain, just the web of the cross sections is fully plasticized, while the flanges remain in the elastic regime. On the contrary, all the points of the lateral surface represent yield conditions for the whole cross-section.
To check the affordability of the proposed approach, the plane domains reported in Figure 4 are compared with those arising from the EC3 international standard and sketched in Figure 6. It is important to emphasize that the EC3 code does not provide any information about the interaction between axial and shear force. Two outcomes can be inferred from the last remark: the first is that no comparison can be made for the N , T domain sketched in Figure 4b; the second is that no limit condition is derived from the EC3 for the N , T domain, and, consequently, the 3D N , T , M is that reported in Figure 7. An examination of Figure 5, Figure 6 and Figure 7 allows us to draw the following conclusions: the proposed domains are fully coherent with those proposed by the international standards, but they represent a rigorous definition of the real yield conditions of the cross-section subjected to the presence of axial force, shear, and bending moment; furthermore, the proposed approach can easily be extended to the case of cross-sections subjected to more complex combinations of internal forces ( N , T y , T z , M x , M y , M z ), but this point will be considered in future studies.
As a second application, the case of the W2 section is considered. In Figure 8, the N , M domain obtained by Equations (13) and (14), the N , T domain obtained by Equations (21) and (27), and the M , T domain obtained by Equations (28) and (29) are sketched, respectively.
In Figure 9, the N , T , M domain is sketched both with and without a bounding surface. In Figure 9a, the plane domains reported in Figure 8 are also reported to enhance the overall coherence of the proposed approach.
As previously stated, the boundary of the upper plateau of the domain in Figure 9b represents yield conditions for the whole cross-section, while all the points within the relevant boundary define limit conditions related to the shear forces, and all the points of the lateral surface represent yield conditions for the whole cross-section.
To check the affordability of the proposed approach, the plane domains reported in Figure 8 are compared with those arising from the EC3 international standard and sketched in Figure 10. As previously stated, the EC3 code does not provide any information about the interaction between axial and shear force, so no comparison can be made for the N , T domain sketched in Figure 8b, and, consequently, the 3D N , T , M related to EC3 is reported in Figure 11. The examination of Figure 9, Figure 10 and Figure 11 allows us to make analogous remarks, as previously stated with reference to Figure 5, Figure 6 and Figure 7.

4. Conclusions

The proposed approach for the identification of the yield domains of I-shaped welded steel sections is rigorous and takes full advantage of the well-known classical formulas usually adopted in beam theory. The comparison with the corresponding domains obtained in agreement with the EC3 standard confirms the affordability of the proposed approach, and it allows the fact that EC3 does not propose the formulation related to the determination of the yield domain for the interaction between axial and shear forces to be overcome. The proposed approach can therefore be adopted for the analysis of I-shaped welded cross-sections of any geometry, allowing for a more reliable design. Possible developments of the proposed study are represented by an extension to the case of the presence of different combinations of internal forces and by the specialization of the approach to the case of rolled sections, as well as to the case of welded cross-sections of any geometry. Moreover, due to the necessary manufacturing production of the welded cross-section beam elements, further in-depth analyses must be performed, consisting of appropriate experimental tests and the related FE computations. All these detailed topics will be the subjects of future papers.

Author Contributions

Conceptualization, S.B. and L.P.; validation, S.B. and L.P.; writing—original draft preparation, L.P.; writing—review and editing, S.B. and L.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Referring to Figure A1,
Figure A1. Geometrical sketch of the cross-section.
Figure A1. Geometrical sketch of the cross-section.
Applsci 14 08037 g0a1
The following geometrical properties of the cross-section can be easily determined:
A = 2 b t f + t w h 2 t f = b h b t w h 2 t f
I y = b t f 3 6 + b t f h t f 2 2 + t w h 2 t f 3 12 = b h 3 12 b t w h 2 t f 3 12
W E y = 2 I y h
W p y = b t f h t f + t w h 2 t f 2 4 = b h 2 4 b t w h 2 t f 2 4
A , I y , W E y and W P y represent the area, the second-order moment of area (with respect to the y -axis), the elastic modulus (with respect to the y -axis) and the plastic modulus (with respect to the y -axis), of the given section, respectively.
Finally, the first-order moment of area of the portion of the area of the given section below the y -axis evaluated with respect to the y -axis, is given by:
S y , G = b t f h t f 2 + t w 2 h 2 t f 2 = b h 2 8 b t w h 2 t f 2 2

References

  1. Gokhfeld, D.A.; Cherniavsky, D.F. Limit Analysis of Structures at Thermal Cycling; Springer: Dordrecht, The Netherland, 1980; ISBN 978-90-286-0455-1. [Google Scholar]
  2. Marti, K. Limit load and shakedown analysis of plastic structures under stochastic uncertainty. Com. Meth. App. Mech. Eng. 2008, 198, 42–51. [Google Scholar] [CrossRef]
  3. Palizzolo, L.; Benfratello, S.; Tabbuso, P.; Vazzano, S. Reliable measures of plastic deformations for elastic plastic structures in shakedown conditions. In Proceedings of the 24th Conference of the Italian Association of Theoretical and Applied Mechanics, AIMETA 2019, Rome, Italy, 15–19 September 2019; ISBN 978-303041056-8. [Google Scholar] [CrossRef]
  4. Chen, W.; Duan, L. Plasticity, Limit Analysis, Stability and Structural Design: An Academic Life Journey from Theory to Practice; World Scientific: Singapore, 2021; ISBN 978-981122974-9/978-981122973-2. [Google Scholar]
  5. Huang, M.F.; Li, Q.; Chan, C.M.; Lou, W.J.; Kwok, K.C.S.; Li, G. Performance-based design optimization of tall concrete framed structures subject to wind excitations. J. Wind. Eng. Ind. Aerodyn. 2015, 139, 70–81. [Google Scholar] [CrossRef]
  6. Tabbuso, P.; Spence, S.M.J.; Palizzolo, L.; Pirrotta, A.; Kareem, A. An efficient framework for the elasto-plastic reliability assessment of uncertain wind excited systems. Struct. Saf. 2016, 58, 69–78. [Google Scholar] [CrossRef]
  7. Benfratello, S.; Di Paola, M.; Palizzolo, L.; Tabbuso, P. Evaluation of the shakedown limit load multiplier for stochastic seismic actions. Meccanica 2017, 52, 2735–2750. [Google Scholar] [CrossRef]
  8. Zhang, R.; Wang, W.; Alam, M.S. Seismic evaluation of friction spring-based self-centering braced frames based on life-cycle cost. Earthq. Eng. Struct. Dyn. 2022, 51, 3393–3415. [Google Scholar] [CrossRef]
  9. Hu, S.; Zhu, S. Life-cycle benefits estimation for hybrid seismic-resistant self-centering braced frames. Earthq. Eng. Struct. Dyn. 2023, 52, 3097–3119. [Google Scholar] [CrossRef]
  10. Banichuk, N.V. Introduction to Optimization of Structures; Springer: New York, NY, USA, 1990; ISBN 978-1-4612-7988-4. [Google Scholar] [CrossRef]
  11. Se-Hyu, C.; Seung-Eock, K. Optimal design of steel frame using practical nonlinear inelastic analysis. Eng. Struct. 2002, 24, 1189–1201. [Google Scholar] [CrossRef]
  12. Degertekin, S.O. Improved harmony search algorithms for sizing optimization of truss structures. Comp. Struct. 2012, 92–93, 229–241. [Google Scholar] [CrossRef]
  13. Benfratello, S.; Giambanco, F.; Palizzolo, L.; Tabbuso, P. Structural design of frames able to prevent element buckling. In Proceedings of the 11th International Conference on Computational Structures Technology, CST 2012, Civil-Comp Proceedings, Dubrovnik, Croatia, 4–7 September 2012; Volume 99, ISBN 978-190508854-6. [Google Scholar]
  14. Benfratello, S.; Giambanco, F.; Palizzolo, L.; Tabbuso, P. Optimal design of steel frames accounting for buckling. Meccanica 2013, 48, 2281–2298. [Google Scholar] [CrossRef]
  15. Palizzolo, L.; Caffarelli, A.; Tabbuso, P. Minimum volume design of structures with constraints on ductility and stability. Eng. Struct. 2014, 68, 47–56. [Google Scholar] [CrossRef]
  16. Palizzolo, L.; Benfratello, S.; Tabbuso, P. Discrete variable design of frames subjected to seismic actions accounting for element slenderness. Comput. Struct. 2015, 147, 147–158. [Google Scholar] [CrossRef]
  17. Kaveh, H.; Nasrollahi, A. Performed-based seismic design of steel frames utilizing charged system search optimization. Appl. Soft Comput. 2014, 22, 213–221. [Google Scholar] [CrossRef]
  18. Benfratello, S.; Palizzolo, L.; Tabbuso, P. Optimal design of elastic plastic frames accounting for seismic protection devices. Struct. Multidiscip. Optim. 2014, 49, 93–106. [Google Scholar] [CrossRef]
  19. Zhang, R.; Hu, S. Optimal design of self-centering braced frames with limited self-centering braces. J. Build. Eng. 2024, 88, 109201. [Google Scholar] [CrossRef]
  20. Benfratello, S.; Palizzolo, L.; Tabbuso, P. Optimization of structures with unrestricted dynamic shakedown constraints. Struct. Multidiscip. Optim. 2015, 52, 431–445. [Google Scholar] [CrossRef]
  21. Colajanni, P.; La Mendola, L.; Monaco, A.; Pagnotta, S. Seismic Performance of Earthquake-Resilient RC Frames Made with HSTC Beams and Friction Damper Devices. J. Earth. Eng. 2022, 26, 7787–7813. [Google Scholar] [CrossRef]
  22. Colajanni, P.; La Mendola, L.; Monaco, A.; Pagnotta, S. Design of RC joints equipped with hybrid trussed beams and friction dampers. Eng. Struct. 2021, 227, 111442. [Google Scholar] [CrossRef]
  23. Colajanni, P.; Pagnotta, S. Friction-based beam-to-column connection for low-damage RC frames with hybrid trussed beams. Struct. Eng. Mech. 2022, 45, 231–248. [Google Scholar] [CrossRef]
  24. Plumier, A.; Doneux, C.; Castiglioni, C.; Brescianini, J.; Crespi, A.; Dell’ Anna, S.; Lazzarotto, L.; Calado, L.; Ferreira, J.; Feligioni, S.; et al. Two Innovations for Earthquake Resistant Design—The INERD Project—Final Report; Research Programme of the Research Fund for Coal and Steel. Steel RTD, Report EUR22044EN; 2006; ISBN 92-79-01694-6. Available online: https://op.europa.eu/en/publication-detail/-/publication/ccd50305-1c81-4eea-b937-fa05f81e2b18/language-en (accessed on 5 July 2024).
  25. Iwankiw, N.R.; Carter, C.J. The dogbone: A new idea to chew on. Mod. Steel Constr. 1996, 36, 18–23. [Google Scholar]
  26. Plumier, A. The dogbone: Back to the future. Eng. J. 1997, 34, 61–67. [Google Scholar] [CrossRef]
  27. Engelhardt, M.D.; Winneberger, T.; Zekany, A.J.; Potyray, T.J. Experimental investigation of dogbone moment connections. Eng. J. 1998, 35, 128–139. [Google Scholar] [CrossRef]
  28. Miller, D.K. Lessons learned from the Northridge earthquake. Eng. Struc. 1998, 20, 249–260. [Google Scholar] [CrossRef]
  29. Mahin, S.T. Lessons from damage to steel buildings during the Northridge earthquake. Eng. Struct. 1998, 20, 261–270. [Google Scholar] [CrossRef]
  30. Shen, J.; Kitjasateanphun, T.; Srivanich, W. Seismic performance of steel moment frames with reduced beam sections. Eng. Struct. 2000, 22, 968–983. [Google Scholar] [CrossRef]
  31. Horton, T.A.; Hajirasouliha, I.; Davison, B.; Ozdemir, Z. More efficient design of reduced beam sections (RBS) for maximum seismic performance. J Constr Steel Res. 2021, 183, 106728. [Google Scholar] [CrossRef]
  32. Shakeri, K.; Akrami, V.; Moradpour, S.; Khedmati, S. Post-earthquake Behavior of Steel Moment Resisting Frames with Connections Modified by Introducing Reduced Beam Section (RBS). Int. J. Steel Struct. 2024, 24, 462–476. [Google Scholar] [CrossRef]
  33. Mirghaderi, S.R.; Shahabeddin, T.; Imanpour, A. Seismic performance of the accordion-web RBS connection. J. Constr. Steel Res. 2010, 66, 277–288. [Google Scholar] [CrossRef]
  34. Saleh, A.; Mirghaderi, S.R.; Zahrai, S.M. Cyclic testing of tubular web RBS connections in deep beams. J. Constr. Steel Res. 2016, 117, 214–226. [Google Scholar] [CrossRef]
  35. Momenzadeh, S.; Kazemi, M.T.; Asl, M.H. Seismic performance of reduced web section moment connections. Int. J. Steel Struct. 2017, 17, 413–425. [Google Scholar] [CrossRef]
  36. Tabar, A.M.; Alonso-Rodriguez, A.; Tsavdaridis, K.D. Building retrofit with reduced web (RWS) and beam (RBS) section limited-ductility connections. J. Constr. Steel Res. 2022, 197, 107459. [Google Scholar] [CrossRef]
  37. Dimakogianni, D.; Dougka, G.; Vayas, I. Seismic behavior of frames with innovative energy dissipation systems (FUSEIS1-2). Eng. Struct. 2015, 90, 83–95. [Google Scholar] [CrossRef]
  38. Valente, M.; Castiglioni, C.A.; Kanyilmaz, A. Welded fuses for dissipative beam-to-column connections of composite steel frames: Numerical analyses. J. Constr. Steel Res. 2017, 128, 498–511. [Google Scholar] [CrossRef]
  39. Avgerinou, S.; Thanopoulos, P.; Hoffmeister, B.; Vayas, I. Seismic resistant buildings with dissipative elements made of high strength steel [Erdbebensichere Bauten mit dissipativen Elementen aus hochfestem Stahl]. Stahlbau 2022, 91, 326–337. [Google Scholar] [CrossRef]
  40. Benfratello, S.; Caddemi, S.; Palizzolo, L.; Pantò, B.; Rapicavoli, D.; Vazzano, S. Targeted steel frames by means of innovative moment resisting connections. J. Constr. Steel Res. 2021, 183, 106695. [Google Scholar] [CrossRef]
  41. Benfratello, S.; Palizzolo, L. Prevention of brittle failure for steel connections utilizing special devices. Structures 2024, 62, 106153. [Google Scholar] [CrossRef]
  42. Benfratello, S.; Palizzolo, L.; Vazzano, S. A New Design Problem in the Formulation of a Special Moment Resisting Connection Device for Preventing Local Buckling. J. Appl. Sci. 2022, 12, 202. [Google Scholar] [CrossRef]
  43. European Committee for Standardization. EN 1993-1-8:2006; Eurocode 3: Design of Steel Structures Part 1–8: Design of Joints; European Committee for Standardization: Bruxelles, Belgium, 2006. [Google Scholar]
  44. Goczek, J.; Supeł, Ł. Resistance of steel cross-sections subjected to bending, shear and axial forces. Eng. Struct. 2014, 70, 271–277. [Google Scholar] [CrossRef]
Figure 1. Geometrical sketch of the cross-section with the acting internal forces.
Figure 1. Geometrical sketch of the cross-section with the acting internal forces.
Applsci 14 08037 g001
Figure 2. Tangential stress ( τ x y , τ x z ) diagrams for T = T p .
Figure 2. Tangential stress ( τ x y , τ x z ) diagrams for T = T p .
Applsci 14 08037 g002
Figure 3. Portion of the cross-section considered: (a) in Equation (7); (b) in Equation (8).
Figure 3. Portion of the cross-section considered: (a) in Equation (7); (b) in Equation (8).
Applsci 14 08037 g003
Figure 4. Yield domains for W1 section: (a) N , M domain; (b) N , T domain; (c) M , T domain.
Figure 4. Yield domains for W1 section: (a) N , M domain; (b) N , T domain; (c) M , T domain.
Applsci 14 08037 g004
Figure 5. Proposed yield domain in N , M , T space for section W1: (a) 3D domain with plane domains evidenced; (b) 3D domain with bounding surface.
Figure 5. Proposed yield domain in N , M , T space for section W1: (a) 3D domain with plane domains evidenced; (b) 3D domain with bounding surface.
Applsci 14 08037 g005
Figure 6. Comparison between the proposed plane domains and those obtained by EC3 for W1 section: (a) N , M domain: (black) EC3, (red) proposed one; (b) M , T domain: (black) EC3, (green) proposed one.
Figure 6. Comparison between the proposed plane domains and those obtained by EC3 for W1 section: (a) N , M domain: (black) EC3, (red) proposed one; (b) M , T domain: (black) EC3, (green) proposed one.
Applsci 14 08037 g006
Figure 7. EC3 domain in N , M , T space for W1 section.
Figure 7. EC3 domain in N , M , T space for W1 section.
Applsci 14 08037 g007
Figure 8. Plane domains for W2 section: (a) N , M domain; (b) N , T domain; (c) M , T domain.
Figure 8. Plane domains for W2 section: (a) N , M domain; (b) N , T domain; (c) M , T domain.
Applsci 14 08037 g008
Figure 9. Proposed domain in N , M , T space for section W1: (a) 3D domain with plane domains evidenced; (b) 3D domain with bounding surface.
Figure 9. Proposed domain in N , M , T space for section W1: (a) 3D domain with plane domains evidenced; (b) 3D domain with bounding surface.
Applsci 14 08037 g009
Figure 10. Comparison between the proposed plane domains and those obtained by EC3 for W2 section: (a) N , M domain: (black) EC3, (red) proposed one; (b) M , T domain: (black) EC3, (green) proposed one.
Figure 10. Comparison between the proposed plane domains and those obtained by EC3 for W2 section: (a) N , M domain: (black) EC3, (red) proposed one; (b) M , T domain: (black) EC3, (green) proposed one.
Applsci 14 08037 g010
Figure 11. EC3 domain in N , M , T space for W2 section.
Figure 11. EC3 domain in N , M , T space for W2 section.
Applsci 14 08037 g011
Table 1. Geometrical characteristics of the welded sections under examination.
Table 1. Geometrical characteristics of the welded sections under examination.
Section Width   b
(mm)
Depth   h
(mm)
Web   Thickness   t w (mm)Flange
Thickness  t f (mm)
W13002908.514
W2170360812.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Palizzolo, L.; Benfratello, S. Yield Behaviour of Welded I-Shaped Steel Cross-Sections. Appl. Sci. 2024, 14, 8037. https://doi.org/10.3390/app14178037

AMA Style

Palizzolo L, Benfratello S. Yield Behaviour of Welded I-Shaped Steel Cross-Sections. Applied Sciences. 2024; 14(17):8037. https://doi.org/10.3390/app14178037

Chicago/Turabian Style

Palizzolo, Luigi, and Salvatore Benfratello. 2024. "Yield Behaviour of Welded I-Shaped Steel Cross-Sections" Applied Sciences 14, no. 17: 8037. https://doi.org/10.3390/app14178037

APA Style

Palizzolo, L., & Benfratello, S. (2024). Yield Behaviour of Welded I-Shaped Steel Cross-Sections. Applied Sciences, 14(17), 8037. https://doi.org/10.3390/app14178037

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop