Finite Element Analysis of Corroded Steel Angles Under Compression

Conference of Junior Researchers in Civil Engineering
163
Finite element analysis of corroded steel angles under compression

Katalin Oszvald
BME Department of Structural Engineering, e-mail: [email protected]
Abstract
Steel structures corrode almost in every environment. The measure of the corrosion can be especially significant if the
maintenance is insufficient. In practice it is necessary to determine the ultimate failure mode and estimate the resistance of the
corroded members to decide it should be replaced or it is enough to remove the corrosion and re-paint the structure. The subject of
this paper is a finite element study of corroded angle section members under centric compression. The joint influence of three
main corrosion parameters is analyzed on the buckling behaviour. The investigated basic elements have different cross-section
sizes and lengths. The behaviour is determined in the function of the parameters of the corrosion and the cross-section
characteristics.
Introduction
Single steel angles section members are used in many various structures as bridges, trusses and latticed
transmission towers. Every structure is exposed to the effects of different environmental influence. The
circumstances, which include the inadequate maintenance lead to corrosion on the structures. The rate of the
corrosion can be very significant and can cause failures in the structures. There are many various corrosion
types with different appearance forms like average, pitting, and crevice corrosion. Because of the diversity of
corrosion, it is difficult to describe it by only one parameter. Corrosion can occur anywhere along the
member length and has various size extensions and rarely extends to the whole member. The place of the
corrosion also can be various within the cross-section. The standards, as e.g. Eurocode [1], do not give any
suggestion [2] how to analyse the corroded members. In the practice the engineers apply an average crosssection decrease, what is supposed along the whole element. In some cases this assumption can be good
approach, but for example in case of localized corrosion it is not accurate enough to determine the behaviour.
In the current research a numerical model is developed and verified by experimental results. Geometrically
and material nonlinear analysis is used to follow the stability phenomena of the compressed angle section.
The resistance of the corroded members are determined, considering the cross-section classification of the
Eurocode. Three different corrosion parameters are studied: (i) cross-section reduction; (ii) extension of
corrosion and (iii) location of corrosion. The structural behaviour is determined in the function of the
corrosion parameters and the cross-section characteristics. In this paper the main focus is on the behaviour of
the corroded members, based on the finite element parametric study.
Previous research
The corrosion is a significant problem in the world, therefore many researcher analyze the effect of it on the
various members of the structures. All of the studies deal with the remaining capacity of the members and
give recommendation how can be assessed the influence of the corrosion. Numerical study was completed
on sheared plates with pitting corrosion by Paik at el. [3]. In the function of the corroded surface and the
plate thickness the ratio of the stress of the corroded and the non-corroded plates was determined. In other
study Rahgozar [4] developed residual capacity curves of corroded I-section on the basis of thickness
reduction. Corrosion on the lower flange and on the web was assumed. Effect of localized corrosion on
buckling plate was investigated by Sadovsy and Drdacky [5]. The influence of pitting corrosion on the hold
frames of bulk carrier was extensively analyzed by Nakai et al. [6, 7]. The research of Heinemeyer and
Feldman [8] focused on the influence of the corrosion on riveted connections. Angle section members, which
are investigated in this study too, were previously analyzed by Beaulieu et al. [2]. The specimens were
corroded by galvanic process and the tests results were compared to the analytical results which were
calculated according to Canadian and American standards. In the test set-up eccentric compression was
applied by gusset plate. Comprehensive study on corroded angle section members is not found in the
literature, where the effect of the main corrosion parameters on the buckling behaviour is analyzed.
Therefore, to determine the behaviour and the resistance decrease of the corroded members a numerical
study program is carried out.
Oszvald, K.: Finite element analysis of corroded steel angles under compression
164
Experimental and numerical background

In the first phase of the research compressive buckling tests were carried out on 24 specimens. The corrosion
was artificially created in the members; part of the leg thickness was eliminated by mechanical process.
During the test centric compression was applied and the obtained maximal load and the vertical and
horizontal displacements were measured. Details of the test can be found in [9]. In parallel with the tests, a
numerical model was developed using Ansys program [10]. The corroded angle section members and the
corrosion appearance in the model can be various by changing the different corrosion parameters. In the
linear and non-linear analyses the applied finite element is a brick element (SOLID 45 of [10]). The size of
the finite elements is half of the original thickness of the legs. The support is hinge connection in the center
of the gravity of the non-corroded elements. The model contains about 100000 elements. By the model the
experimental tests are simulated on the accuracy was verified. Based on the test results the equivalent
geometric imperfection was determined on every specimens of the test. The developed finite element model
is proved to be accurate and efficient to complete a parametric study on centric compression angle members.
Numerical study program

Design buckling resistance in the research is determined by numerical simulation; this process is a general
proposal in the Eurocode standard. In the finite element simulation the nominal yield stress and the
equivalent geometric imperfection are used. In the case of the corroded members the aim is to give a
recommendation based on the results of the numerical study, how can take into account the effect of the
corrosion in the design method. In the first step the behaviour must be determined. Therefore bifurcation
stability analysis is carried out on perfect corroded elements in order to find the critical compressive force
(geometric non-linear buckling analysis GNB). In the GNB analysis the first buckling mode (eigenvector)
is determined. The applied equivalent geometrical imperfection in the geometrically and material non-linear
analyses (GMNI) follows the shape of the determined first mode. In the simulations linear elastic plastic
material model is applied (steel grade of S235 with the nominal yield strength of fy = 235 N/mm2). In the
following parametric study three different cross-section sizes and three element lengths are applied, it means
nine different basic non-corroded elements are analyzed. Table 1 contains the cross-section sizes, the lengths
of the elements, the relative slendernesses and the legs initial b/t (width/thickness) ratios of the noncorroded elements.
ID
Cross-section
[mmmmmm]
Length
[mm]
Slenderness
b/t
RF-1
40404
510
0.7
10
RF-2
40404
840
1.15
10
RF-3
40404
1100
1.5
10
RF-4
60608
750
0.7
7.5
RF-5
60608
1250
1.15
7.5
RF-6
60608
1640
1.5
7.5
RF-7
10010012
1280
0.7
8.3
RF-8
10010012
2100
1.15
8.3
RF-9
10010012
2760
1.5
8.3
Table 1. Basic non-corroded elements
The thickness reduction (Tred), the extension of corrosion (Ext) and the position of the corrosion (pc) along
the whole element are the three main parameters which are considered in the analyses, as corrosion
parameters, defined in Eqs. (1) and (2). Two different basic corrosion patterns (A and B) are applied in the
analyses, as it is shown in Figs. 1 and 2. The crosshatched area marks the corrosion in the member.
165
L/2
0
Fig. 1. A pattern
Fig. 2. B pattern
pc =
Ext =
y
L/2
(1)
x
100
L
(2)
L [mm]: member length; x,y see Fig. 1.

Table 2 shows the applied parameters, in the table the corrosion is marked by black. The corrosion position
is calculated by the ratio of the centre of the corroded area and the half of the member length, Eq. (1). The
length of the corroded area is calculated by extension parameter Eq. (2).
The number of the studied members is about 2000; 238 different corroded cases are analyzed for each noncorroded members.
Pattern Cross - section Tred [%]Ext [%] Corrosion position (pc)
20
20
0.20, 0.47, 0.73, 1.00
30
30
0.30, 0.53, 0.77, 1.00
40
0.40, 0.70, 1.00
50
0.50, 0.75, 1.00
70
70
0.70, 1.00
80
100
1.00
40
50
60
B
Table 2. Corrosion parameters
GNB analyses equivalent geometric imperfections

The observed bucking around the weak axis is shown in Fig. 3. The width-to-thickness ratio of the legs is
changed due to thickness reduction by corrosion and it leads to change the buckling shape and failure mode.
For the relative slenderness of 1.15 and 1.5 the buckling modes follow the first buckling mode. This is also
valid if the thickness reduction is greater than 50%. If the relative slenderness is 0.7 the first mode is not the
same in every case as in the previous elements. Beside the global buckling shape a local buckling is
observed, as shown in Figs. 4 and 5. The length of the local buckling wave is approximately two times the
free width of the outstanding plates if the corrosion is on one leg (pattern A). In case of pattern B the length
of the buckling wave is the length of the corroded area, representing flexural-torsional buckling in the
corroded area, as it is shown in Fig. 6.
The buckling mode is very much dependent on the corrosion parameters. The reduced area can lead to
different modes of buckling due to the different b/t ratios. The reduced cross-section is in class 4 in all cases
if the thickness reduction is greater than 50%, but the b/t ratio is not the same on the studied elements. Table
3 presents the first mode of buckling of part of the studied elements. In RF-1 members with pattern A and
where the Tred is 80% the first buckling mode is local independently of the other two parameters. It is not
valid in case RF-4 and RF-7 members, because the p/t ratio is lower in these cases than in case of RF-1
elements. The position of corrosion also important parameter; if it is closer to the support (ID 171) local
buckling mode is observed by members with lower b/t ratio. But if the corrosion is in the middle (ID 174)
in the cases of RF-4 and RF-7 global bucking mode is observed. The extension of corrosion (Ext) is also an
166
important parameter. Having the same thickness reduction (Tred) and corrosion position (pc), but different
extension (Ext) causes different buckling mode, e.g. on RF-7 element in case of ID 175 and ID 182. This
is almost valid in cases of elements with corrosion pattern B. Local mode is observed already by 70%
thickness reduction, but the number of this type of mode is lower as in the cases of pattern A, beside the
same corrosion parameters.
The shape of the first mode of stability analysis is applied as equivalent geometrical imperfection in the
numerical simulation. The amplitude of the imperfection is L/200 in the analyses, following the
recommendation of the standard and the results of the model verification.
Fig. 3. Ext=20%, pc=0.47, Tred=80
Fig. 4. Ext=70%, pc=1.0, Tred=80
Fig. 5. Ext=70%, pc=1.0, Tred=80
Fig. 6. Ext=70%, pc=1.0, Tred=80
167
ID
Corrosion
pattern
Tred [%]
Ext [%]
pc
RF-1
RF-4
RF-7
171
80
20
0,20
Local
Local
Local
172
80
20
0,47
Local
Global
Local
173
80
20
0,73
Local
Global
Local
174
80
20
1,00
Local
Global
Global
175
80
30
0,30
Local
Global
Local
176
80
30
0,53
Local
Global
Local
177
80
30
0,77
Local
Global
Global
178
80
30
1,00
Local
Global
Global
179
80
40
0,40
Local
Global
Local
180
80
40
0,70
Local
Global
Global
181
80
40
1,00
Local
Global
Global
182
80
50
0,50
Local
Global
Global
183
80
50
0,75
Local
Global
Global
184
80
50
1,00
Local
Global
Global
185
80
70
0,70
Local
Global
Global
186
80
70
1,00
Local
Global
Global
187
80
100
1,00
Local
Global
Global
205
70
20
0,20
Local
Global
Global
209
70
30
0,30
Local
Global
Global
222
80
20
0,20
Local
Local
Local
223
80
20
0,47
Local
Global
Global
225
80
20
1,00
Local
Global
Global
226
80
30
0,30
Local
Global
Local
227
80
30
0,53
Local
Global
Global
228
80
30
0,77
Local
Global
Global
230
80
40
0,40
Local
Global
Global
233
80
50
0,50
Local
Global
Global
Table 3. Bucking modes in the function of corrosion parameters
GMNI analysis behaviour modes

The results of GMNI analyses are evaluated by the behaviour modes function of the different parameters.
Generally the observed failure mode is global flexural buckling about the weak axis on elements with 1.15
and 1.5 relative slenderness. The yield mechanism is developed at the half-length of the member according
to the expectation on the non-corroded members, but on corroded members it is shifted. In these cases the
buckling mode and the failure mode are same. In some cases, when the buckling mode is global, the
behaviour mode in the simulation is different. The observed behaviour is an elastic-plastic failure mode.
After an initial global behaviour the failure mode is local. A third type of behaviour mode is also observed,
in this case the first buckling mode and the failure mode are the same. This is valid if the corrosion is on one
or both of the legs. Figure 7 presents the deformed shape of three different behaviour modes. The first is the
global-global (RF-1-1), the second is global-local (RF-1-52) and the third is local-local (RF-1-222).The
second and the third type of failure modes are observed just on the element with 0.7 relative slenderness. In
two cases the behaviour is clearly global buckling around the y-y axis (the y-y axis is parallel with the leg).
168
The corrosion is on one leg and it is on the whole elements and the Tred is 60% and in the other case Tred is
70%.
RF-1-1
RF-1-52
RF-1-222
Fig. 7. Observed behaviour modes
The load and vertical displacement curve is similar in the case of the three modes. The curves are presented
on the right side of Fig. 8. On the left side of Fig. 8 the displacements of the middle point are shown (the
middle point is considered in the middle of the element length, in the corner edge of the cross-section). The
coordinate axes are the displacements in direction of x and y. The curve of RF-1-52 presents the behaviour;
in the initial phase the line goes together with RF-1-1 an in the moment of failure the direction changes,
because the reduced leg loses its stiffness.
60
4
3
50
2
1
30
RF-1-52
RF-1-1
RF-1-222
20
10
dy [mm]
Load [kN]
40
0
-2
-1 0
-2
RF-1-52
-3
RF-1-1
-4
0
0
0,5
1
Vertical displacement [mm]
1,5
-5
dx [mm]
Fig. 8. Behaviour curves
In GNB analyses the cases are determined when the first buckling is local, in GMNI analyses this statement
is also relevant. Table 4 contains some elements with the ID numbers and corrosion parameters, where the
first buckling mode and the failure mode are different in the GNB and GMNI analyses. The results show that
the b/t ratio is not negligible in the determination of the failure mode.
Changing the parameter b/t (modify Tred) and set to parameter of RF-1 the same behaviour occurs on the
examined RF-4 and RF-7 parameters. By more than 40% thickness reduction the second type of behaviour
mode appears if the b/t ratio is more or equal than 15. This is not a sufficient, but a necessary condition. The
position of the corrosion must be close to the support. As the Tred is getting greater the number of this type of
behaviour is increased; e.g. 70% thickness reduction just in one case the result is not the second mode. By
this scale of Tred the condition, corrosion close to the support, is not important. Necessary condition to local
failure mode is more than 80% thickness reduction and more or equal than 50 b/t ratio. If the thickness
169
reduction is on both of the legs, almost the same conditions are necessary to the various failure modes. In
these cases the second mode does not appear, because after the thickness reduction the symmetry of the
cross-section does not change. The local buckling failure mode is observed, however, by 70% thickness
reduction (more or equal than 35 b/t ratio). On the basis of the completed analyses the resistance of the
corroded members is determined, as a peak value of the applied compressive force.
ID
Cross-section
Tred
Ext
pc
b/t
RF-1-35
40404
40
20
0,2
16,67
RF-1-52
40404
50
20
0,2
20
RF-1-53
40404
50
20
0,47
20
RF-1-56
40404
50
30
0,3
20
RF-4-52
60608
50
20
0,2
15
RF-7-52
10010012
50
20
0,2
16,67
Table 4. Element parameters Global flexural buckling about strong axis
GMNI analysis resistances
1,20
1,20
1,00
1,00
0,80
0,80
Nb /Nb0
Nb /Nb0
The corrosion causes resistance reduction, but the measure of the decrease is different function of the
parameters. Results of RF-1 elements show the tendencies of the resistance decrease function of the different
parameters, as illustrated in Fig. 9. In the left side the horizontal axis is the whole cross-section reduction
(Cred [%]), what can be calculated from Tred. The tendency is almost linear, but there are some differences. In
the different cases the behaviour is the second and the third mode. If the results are plotted in the function of
the b/t ratio, the tendency is non-linear and there are no significant differences when the failure mode is
different. In Fig. 10 part of results are presented in the function of the b/t ratio. The corrosion parameters
what belongs to the plotted results are: Ext = 100% and pc = 1. In the case of the same b/t value the ratio of
the corroded and non-corroded element resistance can be different. The tendencies similar but the measure of
the decrease is greater if the initial b/t is lower.
0,60
0,60
0,40
0,40
0,20
0,20
0,00
0,00
0
50
C red [%]
100
20
40
60
b/t
Fig. 9. Resistances in the function of the Cred and b/t parameters
The development of a design method on the basis of the completed parametric study to determine the
buckling resistance of corroded angle section members is in process. In the method the application rule of
the Eurocode 3 is used with a cross-section dependent reduction factor.
Summary and conclusions

In the current research corroded compressed angle section elements are analysed. A numerical parametric
study is carried out on different corroded elements. In the completed numerical analyses, the effect of crosssection reduction, position and extension of corrosion are analyzed on the stability behaviour. On the basis of
the simulation results the first buckling modes and the failure modes are determined. On the basis of the
numerical study the following conclusions can be made:
170
1,2
RF-1
RF-4
Nb /Nbo
0,8
RF-7
0,6
0,4
0,2
0
5
10
15
b/t ratio
20
25
Fig. 10. Results of RF-1, RF-4 and RF-7 (Ext=100%, pc=1.0)
Three different failure modes are determined, as follows: first is global flexural buckling about weak
axis; second is an elastic-plastic failure mode (initial global behaviour with local failure mode); third is local
buckling. In the case of pattern A all of the modes are appeared, but in the case of pattern B just the first and
the third.
The knowledge of the corrosion parameters is not sufficient to determine the different buckling and
failure modes.
The width to remaining thickness ratio is important characteristic and it is a necessary but not
sufficient parameter to determine the first buckling and the ultimate failure modes.
The observed failure mode is global flexural buckling around the weak axis in the case of relative
slenderness greater than 1.15. Simulation results show that it is relevant in every corroded element,
irrespectively of the corrosion parameters.
The second behaviour mode is observed on the elements where, the b/t ratio is greater than 15 and the
corrosion is located close to the support, and the extension is around 20-30%.
Local failure mode can be observed on elements with greater than 35 b/t ratio. This is also a
necessary condition, but not sufficient.
Acknowledgement
This work is connected to the scientific program of the Development of quality-oriented and harmonized
R+D+I strategy and functional model at BME. This project is supported by the New Szchenyi Plan (Project
ID: TMOP-4.2.1/B-09/1/KMR-2010-0002).
References
[1] EN 1993-1-1:2005 Eurocode 3: Design of steel structures. Part 1-1: General rules and rules for buildings.
[2] Beaulieu LV, Legeron F, Langlois S, Compression strength of corroded steel angle members, Journal of Constructional
Steel Research, vol. 66, pp. 1366-1373, 2010.
[3] Paik JK, Lee JM, Ju Ko M, Ultimate shear strength of plate members with pit corrosion wastage, Thin-Walled Structures,
vol. 42, pp. 1161-1176, 2004.
[4] Rahgozar R, Remaining capacity assessment of corrosion damaged beams using minimum curves, Journal of
Constructional Steel Research, vol. 65, pp. 299-307, 2009.
[5] Sadovsky Z, Drdacky M, Buckling of plate strip subjected to localised corrosion a stochastic model, Thin-Walled
Structures, vol. 39, pp. 247-259, 2001.
[6] Nakai T, Matsushita H, Yamamoto N, Arai H, Effect of pitting corrosion on local strength of hold frames of bulk carrier
(1st report), Marine Structures, vol. 17, pp. 403-432, 2004.
[7] Nakai T, Matsushita H, Yamamoto N, Arai H, Effect of pitting corrosion on local strength of hold frames of bulk carrier
(2nd report) Lateral distortional buckling and local face buckling, Marine Structures, vol. 17, pp. 612-641, 2004.
171
[8] Heinemeyer C, Feldman M, The influence of rivet corrosion on the durability of riveted connections, 6th European
Conference on Steel and Composite Structures, Budapest, Hungary, vol. C, pp. 2217-222, 2011.
[9] Oszvald K, Dunai L, Effect of corrosion on the buckling of steel angle elements. Proc. 8th International PhD Symposium
in Civil Engineering, pp. 549-554, Lyngby, Denmark, 2010.
[10] ANSYS v11.0, Canonsburg, Pennsylvania, USA.

Finite Element Analysis of Corroded Steel Angles Under Compression

Uploaded by

Copyright:

Available Formats

Finite Element Analysis of Corroded Steel Angles Under Compression

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Finite Element Analysis of Corroded Steel Angles Under Compression

Uploaded by

Copyright:

Available Formats

Conference of Junior Researchers in Civil Engineering