Next Article in Journal
Conjugated Polymer/Graphene Oxide Nanocomposites—State-of-the-Art
Next Article in Special Issue
Wrinkle Formation and Initial Defect Sensitivity of Steered Tow in Automated Fiber Placement
Previous Article in Journal
State-of-the-Art Review on Experimental Investigations of Textile-Reinforced Concrete Exposed to High Temperatures
Previous Article in Special Issue
Low-Frequency Magnetoelectric Effects in Magnetostrictive–Piezoelectric Bilayers: Longitudinal and Bending Deformations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Application of Effective Thickness Approach to Box Aluminium Sections

Department of Civil Engineering, University of Salerno, 84084 Fisciano, Italy
*
Author to whom correspondence should be addressed.
J. Compos. Sci. 2021, 5(11), 291; https://0-doi-org.brum.beds.ac.uk/10.3390/jcs5110291
Submission received: 18 October 2021 / Revised: 29 October 2021 / Accepted: 3 November 2021 / Published: 5 November 2021
(This article belongs to the Special Issue Feature Papers in Journal of Composites Science in 2021)

Abstract

:
The ultimate behaviour of aluminium members subjected to uniform compression or bending is strongly influenced by local buckling effects which occur in the portions of the section during compression. In the current codes, the effective thickness method (ETM) is applied to evaluate the ultimate resistance of slender cross-sections affected by elastic local buckling. In this paper, a recent extension of ETM is presented to consider the local buckling effects in the elastic-plastic range and the interaction between the plate elements constituting the cross-section. The theoretical results obtained with this approach, applied to box-shaped aluminium members during compression or in bending, are compared with the experimental tests provided in the scientific literature. It is observed that the ETM is a valid and accurate tool for predicting the maximum resistance of box-shaped aluminium members during compression or in bending.

1. Introduction

In recent decades, the use of aluminium alloys has spread into many engineering fields. In fact, this material is successfully applied in transportation, such as in the rail industry (subway coaches, sleeping cars), automotive industry (car and moto components, moving cranes), and shipping industry (motorboats, sailboats). Aluminium alloys are also successfully used in the civil engineering field where they are employed both as a main element and as a secondary element within civil structures [1,2]. The main advantages of aluminium alloys are related to their lightness, corrosion resistance, and high functionality of structural shapes due to the extrusion process. For these reasons, aluminium material is competitive with steel in the application of structural design. However, used as structural material, it is affected by the same problems concerning the steel members, particularly the instability phenomena which can occur at different levels (Figure 1) and either in the elastic or plastic range. These phenomena are even more relevant in the case of aluminium structures, considering that the Young modulus of aluminium is about one third of that of the steel. In this paper, the attention is focused on studying the effect of local buckling effects on the ultimate response of aluminium members subjected to uniform compression or non-uniform flexural bending. Clearly, local and distortional buckling occurs when the portions of a single aluminium member are subjected to uniform or non-uniform compression. For this reason, the real prediction of the ultimate response of aluminium members becomes very complex and, generally, the proposed methodologies by the current codes underestimate their maximum capacity during compression or bending. Recently, Krishanun Roy et al. [3,4] highlighted that, in the case of channel section members, the design codes (Aluminium Design Manual [5] Australian/New Zealand Standards [6], Eurocode 9 [7] and Eurocode 3 [8]) are conservative within the 20–30% range, with respect to the experimental tests. The same problem was encountered by Feng et al. [9,10] in the case of the circular hollow sections of the aluminium alloy with holes during compression or bending. In the case of compression, the authors provided a new design equation based on the current design formulae that use the effective area method (EAM).
The excessive underestimation of the maximum capacity prediction is related to the current codes which consider only the elastic local buckling, which affects the ultimate resistance of the slender aluminium sections and neglects the interaction of the plate elements constituting the cross-sections and the strain-hardening behaviour of the aluminium alloys. Instead, for thicker sections, a plastic analysis is adopted in a similar way to the steel members. This approach leads to underestimated maximum capacity values with respect to the real behaviour of the aluminium material. Recently, in order to obtain more accurate predictions of the ultimate resistance of the aluminium sections, different approaches were presented within the scientific community. Su et al. [11] evaluated the applicability of the Continuous Strength Method (CSM), which was originally developed for stainless steel stocky (i.e., small width-to-thickness ratio) cross-sections. Zhu and Young [12] modified the current Direct Strength Method (DSM) to achieve more accurate design provisions for flexural SHS members, while Kim and Peköz [13] presented a simplified design approach called the Numerical Slenderness Approach (NSA) for determining the nominal stresses of each plate element of a complex section under flexure. In the case of aluminium box and H-shaped members, the same authors proposed simple mathematical formulations to predict the ultimate flexural resistance and rotation capacity as a function of the main nondimensional slenderness parameters [14,15,16,17,18,19]. These results are currently proposed in the Annex L of new Eurocode 9 draft [7].
Recently, starting from the effective thickness method (ETM) currently adopted by Eurocode 9, an extension of this approach is presented in [17], in order to compute the ultimate resistance of box-shaped members during compression or bending. The new version of ETM considers the plastic interactive buckling of cross-sections and the strain-hardening behaviour of aluminium material. In order to evaluate the accuracy of this simplified approach, a numerical application is provided. In particular, with reference to the box sections, the experimental results presented in [20,21,22] are compared with those obtained by the extension of the effective thickness method.

2. Different Approaches in the Design of Aluminium Section

In the technical literature, there are different simplified methods to evaluate the behaviour of aluminium members under uniform or non-uniform compression, taking into account that the response of thin-walled sections is strongly affected by local instability phenomena which arise in the compressed section. The main approaches, codified and adopted by the European Eurocodes, are: (1) the width effective approach, (2) the reduced strength approach, and (3) the effective thickness approach.
The first approach is very well known, being codified for many years with particular reference to steel structures [8]. It was first introduced by Von Kármán in 1932 [23]. He stated that, for a fixed thickness, a fictitious plate with the width of b e f f would have the critical stress equal to the yield stress. If the actual plate has a larger width, the capacity is the same as that of the fictitious plate. In a plate, the real stress distribution is approximated, or replaced, with two strips which describe the load-carrying, effective width of the plate. Consequently, this method is based on reducing the cross-section area in the parts affected by plate buckling.
The reduced stress method, used in the past in the Aluminium Associated Code [24], verifies the stress level at which a plate part buckles; if a cross-section is built up from multiple plate parts, the lowest stress governs the entire cross-section. Therefore, this method evaluates the capacities of a slender section by considering a reduced value of the limiting stress acting on the full section.
Finally, the effective thickness approach, first introduced in the British Code of Practice for Structural Aluminium [25], was recently introduced in EN1999-1-1, which concerns aluminium alloy structures [7]. This approach is based on replacing the true section with an effective section, obtained to reduce the actual thickness of the compressed parts. The main advantage of the effective thickness approach is that the influence of the heat-affected zones is more easily accounted for in the case of the cross-sections composed by welding. This is very important in the case of aluminium alloy structures where a significant reduction in the material properties arises in the heat-affected zones.
Figure 2 shows a qualitative comparison between the different design approaches for a generic box section subjected to non-uniform compression.
Nowadays, the current European Code Provisions suggest the use of the effective width approach and the effective thickness method for the steel and aluminium cross-section, respectively. In particular, these methods are indirectly used to evaluate the behaviour of the slender sections, belonging to class 4, affected by the local buckling occurring in the elastic region.

3. Extension of Effective Thickness Method (ETM)

3.1. Stability of a Single Plate in the Elastic-Plastic Region

According to a single plate under uniform compression (Figure 3), the elastic critical stress is equal to:
σ c r . e = k π 2 E 12 ( 1 ν e 2 ) ( b / t ) 2
where E is the elastic modulus, ν e is the Poisson’s ratio in the elastic range, b is the plate width, a is the length plate, t is the plate thickness and k is the buckling factor. The factor k   accounts for the edge-restraining conditions and the stress distribution along the loaded edges. In the case of the plate elements constituting the section of a structural member, the occurrence of elastic buckling is also affected by the interaction between the plate elements constituting the member section and by the longitudinal stress gradient. These effects can be considered by modifying Equation (1) using two factors: ζ and α L S . The ζ factor accounts for the interaction between the plate elements that constitute the section, while the α L S factor accounts for the influence of the longitudinal stress gradient occurring in the structural members under non-uniform bending. Clearly, in the case of uniform compression, α L S = 1.00 . Therefore, including the effects of interactive buckling, the longitudinal stress gradient and the non-linear behaviour of material, the elastic buckling stress can be expressed, by rearranging Equation (1) in the following form:
σ c r . e = ζ   α L S k π 2 E 12 ( 1 ν e 2 ) ( b / t ) 2
The correction factor ζ   for the interactive buckling can be evaluated by taking into account that it represents the ratio between the buckling factor k , accounting for interactive buckling, and the buckling factor k 0 , evaluated for the isolated plate element, i.e., ζ = k k 0 . The main expressions are provided by BS5950-5 [26], while, regarding the factor accounting for the influence of the longitudinal stress gradient α L S , the relationships are presented in [27,28].
The occurrence of buckling in the plastic range can be accounted for using a correction factor that depends on the non-linear behaviour of the material. By denoting with ξ such a correction factor, the buckling stress in the plastic range σ c r . p is given by:
σ c r . p = ζ ξ   α L S k π 2 E 12 ( 1 ν 2 ) ( b / t ) 2
Concerning the Poisson’s ratio in the yield region, Gerard and Wildhorn (1952) [29] studied the problem in the case of several aluminium alloys and showed that they were seriously affected by the anisotropy of the material. Therefore, they proposed the following expression:
ν = ν p   E s E ( ν p ν e )
where E s is the secant modulus and ν p represents the Poisson’s ratio in the plastic range, which is equal to 0.50 . Regarding the plastic coefficient ξ , there are many expressions in the scientific literature [30]. The main relationships are reported in the following:
• tangent modulus theory ξ 1 : ξ 1 = E t E
• secant modulus theory ξ 2 : ξ 2 = E s E
• Pearson (1950), Bleich (1952), Vol’Mir (1965) ξ 3 : ξ 3 = E t E
• Gerard (1957) ξ 4 : ξ 4 = E s E t E 2
• Weingarten et al. (1960) ξ 5 : ξ 5 = E s E E t E
•  Stowell (1948), Bijlaard (1949) ξ 6 : ξ 6 = E s E [ 1 3 + 2 3 1 4 + 3 4 E t E s ]
•  Li and Reid (1992) ξ 7 : ξ 7 = E s E [ 1 2 + 1 2 1 4 + 3 4 E t E s ]
where E t is the tangent modulus and n represents the exponent of the Ramberg–Osgood law that describes the behaviour of aluminium alloys. Recently, the use of a different formula was suggested by the same authors for the proposed Annex L of EN 1999-1-1:
ξ 8 = E s E [ n 8 8 + 8 n E t E s ]
It is easy to recognize that Equation (5) is a combination of the secant modulus theory with the Gerard formula. In particular, for the small values of the Ramberg–Osgood exponent n , Equation (5) tends to provide values close to those given by Gerard. Conversely, for high values of n, Equation (5) tends to provide values close to the secant modulus theory.

3.2. Mathematical Steps for the Definition of the New Version of ETM

In this section, the mathematical steps, needed to compute the effective thickness in the non-linear range as a function of the strain level, are reported. In particular, Equation (3) can be rewritten as:
E s ε = ζ ξ   α L S k π 2 E 12 ( 1 ν 2 ) ( b / t ) 2
where ε is the strain corresponding to buckling. So, the effective ratio b / t   can be defined as a function of the strain level as:
( b t ) e f f = ξ E   E s   ζ   α L S k π 2 12 ( 1 ν 2 )     1 ε
The previous equation can be written as:
( t b ) e f f =     E s ξ E 1 ζ   α L S k 12 ( 1 ν 2 ) π 2   ε
by introducing the parameter μ , which accounts the nonlinear behaviour of aluminium material:
μ = ξ E   E s 1 1 ν 2
Therefore, Equation (8) can be rewritten as:
t e f f t = b t   12 π 2   1 μ   1 ζ   α L S k   ε
remembering that:
λ ¯ p = t e f f t
where λ ¯ p is the normalised slenderness in the elastic-plastic region. To use the buckling curves of Eurocode 9 [7], with the normalised slenderness corrected to account for the non-linearity depending on the strain level, it must be considered that:
λ ¯ p = 0.03143   β ϵ 0
where ϵ 0 = ( 250 f 0.2 ) and f 0.2 represents the conventional yield stress of the aluminium material. Therefore:
β ϵ 0 = λ ¯ p 0.03143 = 1 0.03143   3 π 2     ε μ   ζ   α L S   2 k   b t
The factor accounting for the stress distribution along the loaded edge is:
η = 2 k
The final expression used to compute the slenderness parameter of the single-plate element to use with the buckling curves of EN1999-1-1 [7] is:
β ϵ 0 = 17.54   η   b t   ε μ   ζ   α L S    
The last equation computes the effective thickness in the non-linear range as a function of the strain level ε . In fact, according to EN1999-1-1 [7], the reduction factor accounting for local buckling is computed as:
ρ c = 1                     i f             β ϵ 0 1 2 ( C 1 + C 1 2 C 2 ( 3 + ψ ) ) ρ c = C 1 β / ϵ 0 C 2 ( 3 + ψ ) 4   ( β / ϵ 0 ) 2                     i f             β ϵ 0 > 1 2 ( C 1 + C 1 2 C 2 ( 3 + ψ ) )
The parameter ψ accounts for the strain distribution along the loaded edge of the plate. It is given by the ratio between the maximum compression strain at one end of the plate and the strain at the second end of the plate element. In the case of uniform compression, it is ψ = 1 while ψ < 0 when the second end of the plate element is subject to tension. Meanwhile, the values of coefficients C 1 and C 2 are reported in Eurocode 9 and they depend on the class of the buckling curves.

4. Numerical Application to the Box Sections

In this section, the numerical application is provided with reference to the box-shaped members under uniform compression or bending. The aluminium alloy, considered in this procedure, is EN-AW 6082 T6 and, according to the nominal mechanical proprieties presented in Eurocode 9, it is characterized by the conventional yield stress f 0.2 = 260   MPa and the ultimate stress f u = 310   MPa ; the exponent of Ramberg–Osgood is n = 25 . The nominal Young modulus is assumed to be equal to E = 70   GPa while the geometric properties of the box section are the plate width equal to 100   mm and the plate thickness equal to 10   mm .

4.1. Box Sections under Uniform Compression

The application of the effective thickness approach, as previously described, requires a procedure under displacement control. Regarding the stub columns subjected to uniform compression, by increasing the values of the axial displacement δ , the corresponding average strain ε = δ / h is derived; h represents the height of members. Therefore, the slenderness parameter is given by Equation (14) for the increasing values of the strain level and computed for the plate elements constituting the member section. It increases with the increasing values of ε . The appropriate buckling curve is used, according to EN1999-1-1, so that the effective thickness is computed for all the plate elements and the effective cross-section area A e f f is computed. The axial force corresponding to the axial displacement δ is obtained as N = σ   A e f f , where σ is the stress level corresponding to ε , evaluated according to the constitutive stress–strain curve of the material, i.e., the Ramberg–Osgood model:
ε = σ E + 0.002 ( σ f 0.2 ) n
It is important to underline that the aluminium material could be expressed with other mathematical models, for example a two stage R-O model, recently proposed by Yun et al. [31], which appears to be more suitable than the simple R-O law for use in analytical modelling, numerical simulations, and the advanced design of aluminium alloy structures.
The numerical procedure is performed step-by-step by means of the MATLAB program [32] and it is repeated by varying the plastic coefficient ξ according to the relationships presented in Section 3.1.
The generic steps of the procedure are shown in Figure 4a. At the initial step, a generic box section is subjected to the strain level ε i and, consequently, to the stress level σ i . At step i + 1 , the strain level increases, and it is equal to ε i + 1 , where the corresponding stress state is σ i + 1 . By applying Equations (14) and (15), the effective thickness of the box section decreases (red line) with the definition of the reduction coefficient ρ ( ε i + 1 ) .
In Figure 4b, the results of the numerical procedure are depicted in terms of the load–strain curves. It is possible to observe that the tangent modulus theory ( ξ 1 ) provides lower load values, while the higher values are obtained by means of the secant modulus theory ( ξ 2 ). Instead, the results provided by the new formulation, i.e., Equation (5), are located between the Stowell-Bijlaard ( ξ 6 ) and Li-Reid ( ξ 7 ) results.

4.2. Box Sections under Non-Uniform Bending Moment

The same procedure was repeated in the case of the box-shaped aluminium beam under the non-uniform bending moment. The mechanical and geometric properties are the same as those reported in Section 4.1 and, moreover, the beam length is assumed as equal to 2000   mm .
In this case, the iterative procedure is more complicated than the previous case. In fact, the strain level and, consequently, the corresponding stress σ , are not constant, but vary along the section height, as shown in Figure 5a. For this reason, it is not possible to apply a continuous relation, as depicted for the uniform compression, but a fiber model is used to evaluate the curvature χ as a function of the strain level ε . Finally, the bending moment M is evaluated by means of the rotational equilibrium equation between the compression parts and the tension parts. Moreover, the rotation θ corresponding to the attainment of the flexural resistance M was calculated by integrating the curvature diagram χ along the shear length L s of the structural member:
θ = 0 L s χ [ M ( x ) ] d x
Obviously, the rotation θ u , can be estimated by means of Equation (18) for χ [ M ( x ) ] = χ u , i.e., the curvature corresponding to the maximum bending moment M u . Meanwhile, the maximum rotation θ m , corresponding to the point where the moment resistance drops below the conventional yield value M 0.2 , can be computed as:
θ m = θ u + ( χ m χ u ) L p
where χ m is the ultimate curvature and L p is the length of the plastic hinge evaluated as the distance between the points where the conventional yield curvature χ 0.2 and χ m occurs.
The generic steps of the procedure are shown in Figure 5a, while in Figure 5b, the results of the numerical procedure are depicted in terms of the moment–curvature curves. The considerations, made in the case of uniform compression, are still valid for the box-shaped beam under a non-uniform bending moment.

5. Comparison with Experimental Tests

The accuracy of the previous methodology is evaluated by comparing the theoretical results with those provided in the scientific literature. The following analyses were performed by varying the plastic coefficient ξ according to the relationships provided in Section 3. The numerical results are reported in Appendix B. In particular, in Table A3 and Table A5, the experimental and numerical values are shown in Table A3 and Table A5, while the ratio values are shown in Table A4 and Table A6. In the case of uniform compression, the theoretical results are compared with the stub column tests performed by Su et al. [21], while, in the case of bending moment, the theoretical values are compared with the results of the three-point bending tests performed by Moen et al. [20] and Su et al. [22]. For the sake of completeness, the experimental results were compared with those obtained by Eurocode 9, taking into account the local buckling effects for the fourth-class sections. The mean values, the standard deviations σ , and the variation coefficients COV of the ratios between the theoretical and experimental results are reported in Table A4 and Table A6.
By the comparison between the experimental tests and the theoretical methodologies, it is possible to observe that the Eurocode 9 provides underestimated resistance values with respect to the experimental tests and, consequently, the new methodology results are more accurate. Moreover, by observing the theoretical values obtained by varying the plastic coefficient ξ , it is evident that the new expression reported in Equation (5) provides the values of the maximum compressive load and ultimate flexural resistance, which are more similar to the real behaviour of the aluminium members. In fact, in the case of uniform compression, the mean value of the ratio N u . E T M / N u . e x p is about 1.00, while the standard deviation is equal to 0.09. In the case of bending moment, the mean value of the ratio M u . E T M / M u . e x p   is about 0.96 with the standard deviation equal to 0.06. In Figure 6a,b, the comparison between the results of ETM with the plastic coefficient ξ 8 and the experimental results are reported, respectively, in the case of uniform compression and bending moment. Instead, as seen in the numerical application (Section 4), it is evident that the tangent modulus theory ( ξ 1 ) underestimates the resistance values, while the secant modulus theory ( ξ 2 ) overestimates them. The theories of Stowell-Bijlaard ( ξ 6 ) and Li-Reid ( ξ 7 ) provide the values very similar to the experimental one; however, from a stochastic point of view, they are characterized by a greater dispersion of the results.

6. Conclusions

In this work, a numerical application of the extension of the effective thickness method (ETM) is presented. This approach was developed to define a simplified method to predict the ultimate behaviour of aluminium members in compression or in bending. The extension of ETM considers the strain-hardening behaviour of aluminium alloys, the interaction between plate elements constituting a generic cross-section, and, finally, the local buckling in the elastic-plastic region. In order to evaluate the influence of the plastic coefficient on the ultimate resistance of box-shaped aluminium members, an iterative procedure was performed by means of the MATLAB program. It was clear that the tangent modulus theory provided the lower values of the ultimate resistance, while the higher values were obtained by means of the secant modulus theory.
In order to consider the influence of the strain-hardening behaviour of aluminium alloys, a new formulation of the plastic coefficient was proposed. This relationship is a function of the Ramberg–Osgood exponent n ; in this sense, it is the different inelastic behaviour which characterizes each aluminium alloy. The results provided by the new formulation are located between the Stowell-Bijlaard and Li-Reid results.
For evaluating the accuracy of this procedure, a comparison between the numerical results with those presented in the scientific literature is provided. The theoretical results were computed according to Eurocode 9 and the new version of ETM by varying the plastic coefficients.
The values of the ultimate resistance computed by the new methodology are more similar to the corresponding experimental values. Specifically, by comparing the results, it is evident that the Eurocode 9 approach underestimates the maximum resistance of the aluminium members both during compression and bending. Consequently, the current code has a safety advantage, but the aluminium alloys may be uncompetitive from the economic and design point of view. Conversely, the new approach provides the maximum resistance values similar to the experimental results. In fact, in the case of uniform compression, the mean values of the ratios N u . E T M / N u . e x p are nearly 1.00, except when the plastic coefficient is computed by means of the tangent modulus theory ( ξ 1 ) or Weingarten theory ( ξ 5 ) . In case of the prediction of the flexural resistance, the mean values of the ratios M u . E T M / M u . e x p are almost greater than 0.90. Additionally, in this case, when the plastic coefficient is computed by means of the tangent modulus theory ( ξ 1 ) or Weingarten theory ( ξ 5 ) , the theoretical values are underestimated with respect to the experimental results.
Additionally, the theoretical values obtained by applying the new formulation of the plastic coefficient ( ξ 8 ) are more accurate. In fact, the mean value, and the standard deviation of the ratio between the experimental and theoretical results are equal to 0.99 and 0.09 in the case of uniform compression and 0.96 and 0.06 for the beams during bending. Therefore, the new formulation (Equation (5)) provides the theoretical results with a very low dispersion index both during compression and bending.
Consequently, we can conclude that the new version of the effective thickness approach is valid and more accurate for predicting the maximum resistance of box-shaped aluminium members during compression or bending.
In a future development, for evaluating the accuracy of this method to predict the maximum deformation capacity of the members during compression or bending, stub column tests and three-point bending tests could be performed to compare the theoretical load–strain and bending-curvature curves with directly experimental ones. Moreover, in order to define a general procedure to be adopted within the design codes, this new version of ETM and its numerical applications could be extended to the other cross-sections (H-shaped, channel and angle sections) and other aluminium alloys with different heat-treatments.

Author Contributions

Conceptualization, V.P. and A.P.; methodology, V.P. and A.P.; software, A.P.; validation, A.P., V.P. and E.N.; formal analysis, A.P.; investigation, A.P.; resources, A.P.; data curation, A.P.; writing—original draft preparation, A.P., writing—review and editing, E.N.; visualization, A.P.; supervision, E.N. and V.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

According to Figure A1, the main mechanical and geometrical properties of the specimens, considered in the comparison with the theoretical results, are collected in the following Table A1 and Table A2:
Figure A1. Geometric scheme of box section.
Figure A1. Geometric scheme of box section.
Jcs 05 00291 g0a1
Table A1. Experimental properties of stub column tests.
Table A1. Experimental properties of stub column tests.
SpecimenAlloyB1
[mm]
t1
[mm]
B2
[mm]
t2
[mm]
a
[mm]
A
[mm2]
E
[MPa]
f0.2
[MPa]
fu
[MPa]
n
[-]
H64 × 64 × 3A6061 T663.902.8163.902.81191.10686.6566,00023424812
B6061 T663.902.8563.902.85191.50695.9766,00023424812
H70 × 55 × 4.2A6061 T669.904.0854.904.08209.80951.7865,00019320722
B6061 T669.904.0954.904.09209.90953.9565,00019320722
H95 × 50 × 10.5A6061 T694.8010.3649.7010.36284.902564.7271,00022924211
H120 × 70 × 10.5A6061 T6119.9010.3969.9010.39360.003512.2469,00022623810
H120 × 120 × 9A6061 T6120.008.91120.008.91360.203959.2565,00022523413
N95 × 50 × 10.5A6063 T594.9010.3749.7010.37285.202568.8669,00017922010
N120 × 70 × 10.5A6063 T5119.9010.4569.8010.45360.903527.9271,0001391949
N120 × 120 × 9A6063 T5120.008.92120.008.92361.303963.3369,0001812289
Table A2. Experimental properties of three-point bending tests.
Table A2. Experimental properties of three-point bending tests.
SpecimenAlloyB1
[mm]
t1
[mm]
B2
[mm]
t2
[mm]
L
[mm]
A
[mm2]
E
[MPa]
f0.2
[MPa]
fu
[MPa]
n
[-]
Q1-1m-16082 T699.605.94100.305.891000222468,886315.50323.5064.0
Q1-1m-26082 T699.605.94100.305.891000222468,886315.50323.5064.0
Q1-2m-16082 T699.605.94100.305.892000222468,886315.50323.5064.0
Q1-2m-36082 T699.605.94100.305.892000222468,886315.50323.5064.0
Q2-1m-16082 T4100.005.95100.005.881000222566,868176.60283.4038.0
Q2-1m-26082 T4100.005.95100.005.882000222566,868176.60283.4038.0
Q2-2m-16082 T4100.005.95100.005.882000222566,868176.60283.4038.0
Q2-2m-26082 T4100.002.8999.702.831000111066,853120.10221.0026.0
Q3-1m-16082 T4100.002.8999.702.831000111066,853120.10221.0026.0
Q3-1m-36082 T4100.002.8999.702.832000111066,853120.10221.0026.0
Q3-2m-16082 T4100.002.8999.702.832000111066,853120.10221.0026.0
Q3-2m-26082 T4100.105.94100.005.982000224366,880314.00333.4065.0
Q4-2m-17108-T7100.105.94100.005.982000224366,880314.00333.4065.0
Q4-2m-27108-T760.002.29119.402.58100086766,577288.50302.3051.0
R1-1m-16082 T660.002.29119.402.58200086766,577288.50302.3051.0
R1-2m-16082 T660.002.29119.402.58200086766,577288.50302.3051.0
R1-2m-26082 T660.002.29119.402.58300086766,577288.50302.3051.0
R1-3m-16082 T660.002.29119.402.58300086766,577288.50302.3051.0
R1-3m-26082 T660.102.94100.002.94100090666,225281.40290.4045.0
R2-1m-16082 T660.102.94100.002.94100090666,225281.40290.4045.0
R2-1m-26082 T660.102.94100.002.94200090666,225281.40290.4045.0
R2-2m-16082 T660.102.94100.002.94200090666,225281.40290.4045.0
R2-2m-26082 T660.102.94100.002.94300090666,225281.40290.4045.0
R2-3m-16082 T660.102.94100.002.94300090666,225281.40290.4045.0
R2-3m-26082 T669.804.0955.204.0969590667,000207.00222.0016.0
H55 × 70 × 4B36061 T654.704.0969.804.0969395667,000207.00222.0016.0
H95 × 50 × 10B36061 T694.7010.3449.6010.3469595168,000229.00242.0011.0
H50 × 95 × 10B36061 T649.5010.3494.6010.34693255668,000229.00242.0011.0
H64 × 64 × 3B36061 T663.902.8963.802.89693255267,000232.00245.0010.0
H120 × 120 × 9xB36061 T6120.008.90119.908.9069170565,000225.00234.0013.0
H120 × 70 × 10xB36061 T6119.8010.2869.8010.28691395368,000226.00238.0010.0
H70 × 120 × 10xB46061 T669.8010.26119.8010.26692347568,000226.00238.0010.0
H70 × 55 × 4B3-R6061 T669.804.0754.804.07694347065,000193.00207.0022.0
H50 × 95 × 10B3-R6061 T649.5010.3394.7010.3369394868,000229.00242.0011.0
H64 × 64 × 3B3-R6061 T663.902.8363.902.83696255267,000232.00245.0010.0
N120 × 70 × 10B36063 T5120.0010.4069.9010.4068969171,000139.00194.009.0
N70 × 120 × 10B36063 T569.9010.40119.9010.40688351771,000139.00194.009.0
N120 × 120 × 9B36063 T5119.908.90119.908.90693351569,000181.00228.009.0

Appendix B

In this section, the numerical comparison between the theoretical and experimental results are collected in the following tables. In particular, the experimental results presented in the scientific literature are reported in the first column; the theoretical results computed according to Eurocode 9 are listed in the second column. In the last columns, the values derived by the new version of ETM are reported in the last columns. Specifically, they are computed by varying the plastic coefficients according to the relationships provided in Section 3. Moreover, the mean value [ μ ] , the standard deviation [ σ ] and the coefficient of variation [ COV ] of the ratio between the theoretical and experimental results are reported in the final part of Table A4 and Table A6.
Table A3. Comparison between the theoretical and experimental results in the case of uniform compression.
Table A3. Comparison between the theoretical and experimental results in the case of uniform compression.
SpecimenNu.exp [kN]Nu.EC9 [kN]Nu.ETM [kN]
ξ1ξ2ξ3ξ4ξ5ξ6ξ7ξ8
H64 × 64 × 3A164.2150.68129.07154.95138.10136.50135.18145.12147.57158.13
B165.4152.86131.54158.33141.07139.34137.92148.37150.87160.54
H70 × 55 × 4.2A196.2183.69167.99193.69179.37176.82175.13188.61190.26191.31
B196.9184.11168.37194.30179.82177.37175.53189.21190.79191.74
H95 × 50 × 10.5A626.2587.32620.66776.37783.85690.55651.70750.18757.69722.86
H120 × 70 × 10.5A862.5793.77795.711008.01980.70886.35843.13968.12979.62924.23
H120 × 120 × 9A981.5890.83839.361027.42959.40914.06888.85991.791001.69987.83
N95 × 50 × 10.5A609.8459.83494.87626.53649.24554.78516.96603.68610.10574.76
N120 × 70 × 10.5A736.9490.38504.22646.89667.38568.00528.23620.91627.97603.69
N120 × 120 × 9A811.1717.36693.58889.77857.29774.95735.97850.14861.62812.48
Table A4. Ratio values between the theoretical and experimental results in the case of uniform compression.
Table A4. Ratio values between the theoretical and experimental results in the case of uniform compression.
SpecimenNu.EC9/Nu.expNu.ETM/Nu.exp
ξ1ξ2ξ3ξ4ξ5ξ6ξ7ξ8
H64 × 64 × 3A0.920.790.940.840.830.820.880.900.96
B0.920.800.960.850.840.830.900.910.97
H70 × 55 × 4.2A0.940.860.990.910.900.890.960.970.98
B0.940.860.990.910.900.890.960.970.97
H95 × 50 × 10.5A0.940.991.241.251.101.041.201.211.15
H120 × 70 × 10.5A0.920.921.171.141.030.981.121.141.07
H120 × 120 × 9A0.910.861.050.980.930.911.011.021.01
N95 × 50 × 10.5A0.750.811.031.060.910.850.991.000.94
N120 × 70 × 10.5A0.670.680.880.910.770.720.840.850.82
N120 × 120 × 9A0.880.861.101.060.960.911.051.061.00
Mean value [μ]:0.880.841.030.990.920.880.991.000.99
Standard deviation [σ]:0.090.080.110.130.100.090.110.110.09
Coefficient of variation [COV]:0.110.100.110.130.100.100.110.110.09
Table A5. Comparison between the theoretical and experimental results in the case of bending moment.
Table A5. Comparison between the theoretical and experimental results in the case of bending moment.
SpecimenMu.exp [kNm]Mu.EC9 [kNm]Mu.ETM [kNm]
ξ1ξ2ξ3ξ4ξ5ξ6ξ7ξ8
H70 × 55 × 4.2B34.754.164.766.095.645.184.944.945.484.75
H55 × 70 × 4.2B36.764.634.705.104.834.814.644.644.726.56
H95 × 50 × 10.5B312.0912.6812.9113.9112.9112.9112.8412.8412.9111.33
H50 × 95 × 10.5B321.0917.0523.6325.0424.0723.7423.5723.5723.9519.25
H64 × 64 × 3.0B34.13.724.104.564.494.304.124.124.354.17
H120 × 120 × 9.0xB344.4238.4745.3646.3541.4242.2742.6445.6446.3441.79
H120 × 70 × 10.50xB323.5922.9225.3626.3526.4226.2725.6425.6426.3422.41
H70 × 120 × 10.50xB437.8629.5337.9337.9137.9337.9337.5937.5937.9232.63
H70 × 55 × 4.2B3-R4.823.834.294.494.294.294.184.184.294.75
H64 × 64 × 3.0B3-R10.283.653.554.053.983.663.493.493.8611.33
N120 × 70 × 10.5B320.7214.316.4320.3621.3917.9216.4016.4019.7918.17
N70 × 120 × 10.5B337.318.4316.3831.0231.3028.9626.0126.0130.8029.36
N120 × 120 × 9.0B340.5330.9326.6042.4043.0338.2035.1141.1141.4542.01
Q1-1m-127.8722.4223.4127.7424.5924.2124.1126.1127.0625.44
Q1-1m-227.8722.4223.4127.7424.5924.2124.1124.1127.0625.44
Q1-2m-127.1522.4223.4127.6324.4824.0723.9724.9726.9225.32
Q1-2m-327.3922.4223.4127.6324.4824.0723.9724.9726.9225.32
Q2-1m-118.4312.5516.5319.2717.4716.9916.7617.7618.9516.78
Q2-1m-218.312.5516.5319.2717.4716.9916.7617.7618.9516.78
Q2-2m-117.7612.5516.5319.2717.4716.9916.7617.7618.9516.78
Q2-2m-218.312.5516.5319.2717.4716.9916.7617.7618.9516.78
Q3-1m-14.784.33.764.764.404.154.054.054.434.67
Q3-1m-34.694.33.754.744.384.134.034.034.434.67
Q3-2m-14.644.33.754.744.384.134.034.034.434.65
Q3-2m-24.874.33.754.744.384.134.034.034.434.65
Q4-2m-128.8122.3922.9027.5524.1824.1124.0425.0426.8225.22
Q4-2m-227.8522.3922.9027.5524.1824.1124.0425.0426.8225.22
R1-1m-18.17.537.629.148.288.138.018.318.987.90
R1-2m-18.17.537.629.148.288.138.018.318.987.86
R1-2m-28.17.537.629.138.278.138.008.318.977.86
R1-3m-17.627.537.629.138.278.138.008.318.977.85
R1-3m-27.547.536.968.657.857.607.407.808.457.85
R2-1m-18.657.226.968.657.857.607.407.808.458.35
R2-1m-28.657.226.968.627.837.587.387.788.428.35
R2-2m-18.587.226.968.627.837.587.387.788.428.33
R2-2m-28.87.226.968.617.827.577.387.888.418.33
R2-3m-19.037.226.968.617.827.577.387.888.418.33
R2-3m-28.957.226.968.617.827.577.387.888.418.33
Table A6. Ratio values between the theoretical and experimental results in the case of bending moment.
Table A6. Ratio values between the theoretical and experimental results in the case of bending moment.
SpecimenMu.EC9/Mu.expMu.ETM/Mu.exp
ξ1ξ2ξ3ξ4ξ5ξ6ξ7ξ8
H70 × 55 × 4.2B30.881.001.281.191.091.041.041.151.00
H55 × 70 × 4.2B30.680.700.750.710.710.690.690.700.97
H95 × 50 × 10.5B31.051.071.151.071.071.061.061.070.94
H50 × 95 × 10.5B30.811.121.191.141.131.121.121.140.91
H64 × 64 × 3.0B30.911.001.111.101.051.011.011.061.02
H120 × 120 × 9.0xB30.871.021.040.930.950.961.031.040.94
H120 × 70 × 10.50xB30.971.081.121.121.111.091.091.120.95
H70 × 120 × 10.50xB40.781.001.001.001.000.990.991.000.86
H70 × 55 × 4.2B3-R0.790.890.930.890.890.870.870.890.99
H64 × 64 × 3.0B3-R0.360.350.390.390.360.340.340.381.10
N120 × 70 × 10.5B30.690.790.981.030.860.790.790.950.88
N70 × 120 × 10.5B30.490.440.830.840.780.700.700.830.79
N120 × 120 × 9.0B30.760.661.051.060.940.871.011.021.04
Q1-1m-10.800.841.000.880.870.870.940.970.91
Q1-1m-20.800.841.000.880.870.870.870.970.91
Q1-2m-10.830.861.020.900.890.880.920.990.93
Q1-2m-30.820.851.010.890.880.880.910.980.92
Q2-1m-10.680.901.050.950.920.910.961.030.91
Q2-1m-20.690.901.050.950.930.920.971.040.92
Q2-2m-10.710.931.090.980.960.941.001.070.94
Q2-2m-20.690.901.050.950.930.920.971.040.92
Q3-1m-10.900.791.000.920.870.850.850.930.98
Q3-1m-30.920.801.010.930.880.860.860.951.00
Q3-2m-10.930.811.020.940.890.870.870.961.00
Q3-2m-20.880.770.970.900.850.830.830.910.95
Q4-2m-10.780.790.960.840.840.830.870.930.88
Q4-2m-20.800.820.990.870.870.860.900.960.91
R1-1m-10.930.941.131.021.000.991.031.110.98
R1-2m-10.930.941.131.021.000.991.031.110.97
R1-2m-20.930.941.131.021.000.991.031.110.97
R1-3m-10.991.001.201.091.071.051.091.181.03
R1-3m-21.000.921.151.041.010.981.031.121.04
R2-1m-10.830.801.000.910.880.860.900.980.97
R2-1m-20.830.801.000.910.880.850.900.970.97
R2-2m-10.840.811.000.910.880.860.910.980.97
R2-2m-20.820.790.980.890.860.840.900.960.95
R2-3m-10.800.770.950.870.840.820.870.930.92
R2-3m-20.810.780.960.870.850.820.880.940.93
Mean value [μ]:0.820.851.020.940.910.890.920.990.95
Standard deviation [σ]:0.130.150.140.130.130.130.140.140.06
Coefficient of variation [COV]:0.160.180.140.140.140.150.150.140.06

References

  1. Mazzolani, F.M. 3D aluminium structures. Thin-Walled Struct. 2012, 61, 258–266. [Google Scholar] [CrossRef]
  2. Georgantzia, E.; Gkantou, M.; Kamaris, G.S. Aluminium alloys as structural material: A review of research. Eng. Struct. 2021, 227, 111372. [Google Scholar] [CrossRef]
  3. Roy, K.; Chen, B.; Fang, Z.; Uzzaman, A.; Chen, X.; Lim, J.B.P. Local and distortional buckling behaviour of back-to-back built-up aluminium alloy channel section columns. Thin-Walled Struct. 2021, 163, 107713. [Google Scholar] [CrossRef]
  4. Roy, K.; Chen, B.; Fang, Z.; Uzzaman, A. Axial capacity of back-to-back built-up aluminium alloy channel section column. J. Struct. Eng. 2022. [Google Scholar] [CrossRef]
  5. The Aluminium Association. Aluminium Design Manual; The Aluminium Association: Washington, DC, USA, 2020. [Google Scholar]
  6. Australian/New Zealand Standard. Aluminium Structures Part 1: Limit Stade Design; Australian/New Zealand Standard: Sydney, Australia, 1997.
  7. EN1999-1-1. Eurocode 9: Design of Aluminium Structures—Part 1-1: General Structural Rules; European Committee for Standardization: Brussels, Belgium, 2007. [Google Scholar]
  8. EN 1993-1-1. Eurocode 3: Design of Steel Structures—Part 1-1: General Rules and Rules for Buildings; European Committee for Standardization: Brussels, Belgium, 2005. [Google Scholar]
  9. Feng, R.; Mou, X.; Chen, Z.; Roy, K.; Chen, B.; Lim, J.B.P. Finite-element modelling and design guidelines for axial compressive capacity of aluminium alloy circular hollow sections with holes. Thin-Walled Struct. 2020, 157, 107027. [Google Scholar] [CrossRef]
  10. Feng, R.; Chen, Z.; Shen, C.; Roy, K.; Chen, B.; Lim, J.B.P. Flexural capacity of perforated aluminium CHS tubes–An experimental study. Structures 2020, 25, 463–480. [Google Scholar] [CrossRef]
  11. Su, M.N.; Young, B.; Gardner, L. The continuous strength method for the design of aluminium alloy structural elements. Eng. Struct. 2016, 122, 338–348. [Google Scholar] [CrossRef]
  12. Zhu, J.H.; Young, B. Design of aluminum alloy flexural members using direct strength method. J. Struct. Eng. 2009, 135, 558–566. [Google Scholar] [CrossRef]
  13. Kim, Y.; Peköz, T. Numerical Slenderness Approach for design of complex aluminium extrusion subjected to flexural loading. Thin-Walled Struct. 2018, 127, 62–75. [Google Scholar] [CrossRef]
  14. Piluso, V.; Pisapia, A.; Nastri, E.; Montuori, R. Ultimate resistance and rotation capacity of low yielding high hardening aluminium alloy beams under non-uniform bending. Thin-Walled Struct. 2019, 135, 123–136. [Google Scholar] [CrossRef]
  15. Pisapia, A.; Piluso, V.; Nastri, E.; Montuori, R. Ultimate behaviour of high yielding low hardening aluminium alloys I-beams. Thin-Walled Struct. 2020, 146, 106463. [Google Scholar]
  16. Nastri, E.; Piluso, V. The influence of strain-hardening on the ultimate behaviour of aluminium RHS-beams under moment gradient. Thin-Walled Struct. 2020, 157, 107091. [Google Scholar] [CrossRef]
  17. Piluso, V.; Pisapia, A. Interactive Plastic Local Buckling of Box-shaped Aluminium Members under Uniform Compression. Thin-Walled Struct. 2021, 164, 107828. [Google Scholar] [CrossRef]
  18. Castaldo, P.; Nastri, E.; Piluso, V. FEM simulations and rotation capacity evaluation for RHS temper T4 aluminium alloy beams. Comp. Part B Eng. 2017, 115, 124–137. [Google Scholar] [CrossRef]
  19. Castaldo, P.; Nastri, E.; Piluso, V. Ultimate behaviour of RHS temper T6 aluminium alloy beams subjected to non-uniform bending: Parametric analysis. Thin-Walled Struct. 2017, 115, 129–141. [Google Scholar] [CrossRef]
  20. Moen, L.A.; Hopperstad, O.S.; Langseth, M. Rotational capacity of aluminum beams under moment gradient. I: Experiments. J. Struct. Eng. 1999, 125, 910–920. [Google Scholar] [CrossRef]
  21. Su, M.N.; Young, B.; Gardner, L. Deformation-based design of aluminium alloy beams. Eng. Struct. 2014, 80, 339–349. [Google Scholar] [CrossRef] [Green Version]
  22. Su, M.N.; Young, B.; Gardner, L. Testing and Design of Aluminium Alloy Cross Sections in Compression. J. Struct. Eng. 2014, 140, 04014047. [Google Scholar] [CrossRef] [Green Version]
  23. Von Kármán, T.; Sechler, E.E.; Donnell, L.H. The Strength of Thin Plate in Compression. Trans. Am. Soc. Mech. Eng. 1932, 54, 53–58. [Google Scholar]
  24. Landolfo, R.; Mazzolani, F.M. Different Approaches in the Design of Slender Aluminium Alloy Sections. Thin-Walled Struct. 1997, 27, 85–102. [Google Scholar] [CrossRef]
  25. BS 8118-1. The Structural use of Aluminium; British Standard Institute: London, UK, 1991. [Google Scholar]
  26. BS 5950-5. Structural Use of Steelwork in Building—Formed Thin Gauge Sections Part 5: Code of Practice for Design of Cold Formed Thin Gauge Sections; Steel Construction Institute: London, UK, 1998. [Google Scholar]
  27. Rebiano, R.; Silvestre, N.; Camotim, D. GBT Formulation to Analyze the Buckling Behavior of Thin-Walled Members subjected to Non-Uniform Bending. Int. J. Struct. Stab. Dyn. 2007, 7, 23–54. [Google Scholar]
  28. Yu, C.; Schafer, B.W. Effect of Longitudinal Stress Gradients on Elastic Buckling of Thin Plates. J. Eng. Mech. 2007, 133, 452–463. [Google Scholar] [CrossRef]
  29. Gerard, G.; Wildhorn, S. A Study of Poisson’s Ratio in the Yield Region; NACA Tech.: Washington, DC, USA, 1952; No. 2561. [Google Scholar]
  30. Stowell, E.Z. A Unified Theory of Plastic Buckling of Columns and Plates; NACA Tech.: Washington, DC, USA, 1948; No.1556. [Google Scholar]
  31. Yun, X.; Wang, Z.; Gardner, L. Full-range stress-strain curves for aluminium alloys. J. Struct. Eng. 2021, 147, 04021060. [Google Scholar] [CrossRef]
  32. MathWorks Inc. MATLAB-High Performance Numeric Computation and Visualization Software. In User’s Guide; MathWorks Inc.: Natick, MA, USA, 1997. [Google Scholar]
Figure 1. Instability phenomena on the metal structures.
Figure 1. Instability phenomena on the metal structures.
Jcs 05 00291 g001
Figure 2. Comparison of the different design approaches for a rectangular section in bending.
Figure 2. Comparison of the different design approaches for a rectangular section in bending.
Jcs 05 00291 g002
Figure 3. Single plate under uniform compression.
Figure 3. Single plate under uniform compression.
Jcs 05 00291 g003
Figure 4. Scheme of a generic box section in compression (a). The load-strain curves by varying the plastic coefficient ξ   (b).
Figure 4. Scheme of a generic box section in compression (a). The load-strain curves by varying the plastic coefficient ξ   (b).
Jcs 05 00291 g004
Figure 5. Scheme of a generic box section in bending (a). The moment-curvature curves by varying the plastic coefficient ξ   (b).
Figure 5. Scheme of a generic box section in bending (a). The moment-curvature curves by varying the plastic coefficient ξ   (b).
Jcs 05 00291 g005
Figure 6. Comparison between theoretical and experimental results in the case of uniform compression (a). Comparison between theoretical and experimental results in the case of non-uniform bending (b).
Figure 6. Comparison between theoretical and experimental results in the case of uniform compression (a). Comparison between theoretical and experimental results in the case of non-uniform bending (b).
Jcs 05 00291 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nastri, E.; Piluso, V.; Pisapia, A. Numerical Application of Effective Thickness Approach to Box Aluminium Sections. J. Compos. Sci. 2021, 5, 291. https://0-doi-org.brum.beds.ac.uk/10.3390/jcs5110291

AMA Style

Nastri E, Piluso V, Pisapia A. Numerical Application of Effective Thickness Approach to Box Aluminium Sections. Journal of Composites Science. 2021; 5(11):291. https://0-doi-org.brum.beds.ac.uk/10.3390/jcs5110291

Chicago/Turabian Style

Nastri, Elide, Vincenzo Piluso, and Alessandro Pisapia. 2021. "Numerical Application of Effective Thickness Approach to Box Aluminium Sections" Journal of Composites Science 5, no. 11: 291. https://0-doi-org.brum.beds.ac.uk/10.3390/jcs5110291

Article Metrics

Back to TopTop