AP05_3.vp 1 Introduction The mechanical properties of stainless steels are signifi- cantly different from those of carbon steel. Stainless steels display a pronounced response to cold working resulting in anisotropic, non-linear stress-strain behaviour and low pro- portional limits. The material properties of various stainless steels have been thoroughly investigated since the 1960s by a number of investigators, e.g. refs. [1, 2, 3, 4]. It has been generally concluded that the stress-strain behaviour of stain- less steels can be best described by the Ramberg-Osgood model [5], and Hill’s [6] modified form of the Ramberg- -Osgood equation is used in design specifications. The main design specification for cold formed stainless steel members in the USA is the ASCE specification [7] and in Europe, Eurocode 3: Part 1.4 [8] has been recently developed and is still under examination. The two codes use different approaches when dealing with the mechanical properties of the material. The ASCE code employs the modified form of the Ramberg-Osgood model to describe the stress-strain be- haviour of a material, whereas the Eurocode relies for most purposes on the specification of a linear stress-strain law, with the yield strength taken as the 0.2% proof stress. Refs. [9, 10] show a comparison of the Eurocode and ASCE code load ca- pacity predictions for lipped channel columns is illustrated. The simpler Eurocode analysis has been found to give rea- sonable estimates of concentrically loaded column strength without taking account of the non-linearity of the stress-strain curve. As part of a previous investigation ref. [11], a series of tensile tests were carried out on coupons cut from stainless steel lipped channel sections, and also on full sections, and the stress-strain characteristics are examined in this paper and incorporated into a non-linear finite element analysis. 2 Mechanical properties of stainless steel lipped channel members In the formation of a profiled section, the cold working occurs in localised areas, with the material at the bends being strain hardened. Therefore the properties of the material vary throughout the cross-section where at the formed bends, higher yield and tensile strengths exist, leading to a more complex stress-strain relationship for cold formed members, and in particular, for stainless steel members. The level of increase of both yield and tensile strength is dependent on the ratio of corner radius to material thickness (r/t). The cold formed lipped channels under investigation are of stainless steel, of cross-sections with small web, flange and lip dimen- sions and are considered to be thick and hence four corner bends are formed with small r/t ratios (<1). These four cor- ners will have an effect on the stress-strain response of the material obtained from a full section test, which could then be compared to that obtained for virgin material from a standard tensile test. Also, most commercially available finite element programs allow for a non-linear analysis and hence the inclu- sion of the actual stress-strain data obtained from tensile testing and from existing theories. The ASCE design specification adopts the modified form of the Ramberg-Osgood formula given by equation (1). It is a three-parameter equation for expressing the relationship between the stress and strain for stresses up to a value slightly greater than the yield strength of the material. � � � � � � � � � � E K E n , (1) where � � unit strain � � unit stress (N/mm2) E � modulus of elasticity (N/mm2) K and n are constants for a given curve, which are evaluated through two secant yield strength values for slopes of 0.7 E and 0.85 E. Equation (1) was modified by Hill [6], and, in- stead of using secant yield strengths, K and n can be evaluated in terms of two yield strength values: (i) �1 at an offset �1; (ii) �2 at an offset �2. Using the most common offset of 0.002 for the yield stress (�2) and assuming that the modulus of elasticity E is equal to the initial value E0, equation (1) becomes: � � � � � � � � � � � � E y n 0 0 002. . (2) The ASCE design code makes use of this modified equa- tion 2, and the three points on the stress-strain curve are defined as: (i) the origin; (ii) the point of 0.2 % proof stress; 92 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 3/2005 Czech Technical University in Prague Finite Element Modelling of Cold Formed Stainless Steel Columns M. Macdonald, J. Rhodes This paper describes the results obtained from a finite element investigation into the load capacity of column members of lipped channel cross-section, cold formed from Type 304 stainless steel, subjected to concentric and eccentric compression loading. The main aims of this investigation were to determine the effects which the non-linearity of the stress-strain behaviour of the material would have on the column behaviour under concentric or eccentric loading. Stress-strain curves derived from tests and design codes are incorporated into non-linear finite element analyses of eccentrically loaded columns and the results obtained are compared with those obtained on the basis of experiments on stainless steel channel columns with the same properties and dimensions. Comparisons of the finite element results and the test results are also made with existing design specifications and conclusions are drawn on the basis of the comparisons. Keywords: finite elements, stainless steel, cold forming. (iii) another offset strength (e.g. 0.01 %). If these points are substituted into equation (2), then ‘n’ can be evaluated. The term ‘n’ is referred to in the ASCE design code as the plasticity factor. The accuracy of the above method is largely based on how well the analytical equation fits the stress-strain relation- ship of the material. The code lists for particular grades of stainless steel, tables of yield stress, tangent modulus and plas- ticity factors. The results obtained for the stress-strain relationship from both virgin material and full section tensile tests will be used for comparison with the results obtained from the above ASCE Ramberg-Osgood approach and by a trial and error ‘best fit’ method using the experimental stress-strain curves. These will then be incorporated into a non-linear finite ele- ment analysis of eccentrically loaded stainless steel columns. 3 Load capacity of stainless steel lipped channel columns subjected to combined bending and axial compression loading Rhodes et. al. [9, 10], investigated both concentric and eccentric loading of cold formed stainless steel lipped chan- nel section columns. The findings showed that the relevant design codes refs [7, 8] provided very accurate predictions of load capacity for the concentric loading case using both virgin and full section material properties when compared to exper- imental results. A finite element analysis also produced a very accurate correlation to both the experimental results and the design code predictions. However, for shorter length eccentrically loaded columns, the design codes were very con- servative in their prediction of load capacity using both virgin and full section material properties. It was concluded that the design codes’ interaction formulae were inadequate in pre- dicting the load capacity of short-to-medium length columns. The Eurocode (8) interaction formula to determine the axial strength Nsd is given by equation (3). N f A N e M y M n M sd sd � � � �1 1 1 � � �� � � � � � �� � � , (3) where all terms are defined in Ref. [9]. The ASCE [7] interaction formula to determine the axial strength Pu is given by equation (4). P P P e M P P u n u n u E � � � � �� � � � � � � � � 1 1 0. . (4) where all terms are defined in Ref. [9]. Both equations produced very conservative estimates of load capacity and an attempt to improve the interaction formulae was proposed by Macdonald Ref. [12]. This modifi- cation to the interaction formulae involved replacing the linear moment capacity Mn with the true moment capacity of the lipped channel cross-section Mexp obtained from bending tests. Hence the ASCE interaction formula was modified as given by equation (5). P P P e M P P u n u u E � � � � �� � � � � � � � � exp . 1 1 0 . (5) The Eurocode interaction formula was modified as given by equation (6). N f A N e M y M M sd sd � � � �1 1 1 � � �� � � � � � �� � � exp . (6) In both equations (5) and (6), Mexp is the cross-section true moment capacity and the 0.2 % proof stress is taken from the full section tensile test results, and all other terms are defined in Ref. [12]. 4 Finite element analysis A finite element analysis was then performed using the ANSYS software package to determine the load capacity of concentrically and eccentrically loaded columns. Two types of buckling analyses are available within the ANSYS package – eigenvalue analysis and non-linear analysis. Both types of analyses were used, however, the eigenvalue analysis was used only to verify that the finite element model boundary con- ditions (i.e. column pin ends and eccentric loading) were accurate, as this type of analysis takes no account of material non-linearity. A full non-linear analysis was conducted using shell elements (ANSYS SHELL181) which are four-noded elements with six degrees of freedom at each node, i.e. trans- lations in x, y and z directions, and rotations about x, y and z axes. Fig. 2 shows a typical displacement and boundary condi- tions plot for an eccentrically loaded column, while Fig. 3 shows a typical nodal stress plot. For the eccentrically loaded columns, the non-linear material properties of the stainless steel were defined in ANSYS using the initial elastic modulus, Poisson’s ratio and stress-strain data obtained from: (i) coupon tensile tests on material cut from section webs; (ii) full-section tensile tests; (iii) ASCE (Ramberg-Osgood) approach; (iv) ‘best fit’ stress- -strain curves. For concentrically loaded columns, the virgin material properties were used. The non-linear solution breaks the load up into a series of load increments that can be applied over several load steps. At the completion of each incremental solution, the program adjusts the stiffness matrix to reflect the non-linear changes in the structural stiffness before proceeding to the next load in- crement. For a non-linear buckling analysis to be accurate us- ing ANSYS, it is necessary to set an initial imperfection in the structure being modelled. This was achieved by modelling a very small mid-span deflection which produced a very large radius of curvature for the lipped channel columns which would approximate any actual imperfections. A parametric model was constructed by defining positions of keypoints to allow for easy alterations to the model, partic- ularly for the variation in column length. A column half-model was modelled using appropriate symmetry com- mands that helped to reduce the considerable computer pro- cessing time. © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 93 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 3/2005 5 Experimental investigations 5.1 Tensile tests Fig. 1 shows a schematic diagram of a typical cold formed stainless steel lipped channel member under investigation. The member is commercially available with specimens cut from supplied lengths and all the specimens were accurately measured at a number of points, with the values averaged to obtain the finished dimensions, and all calculations were based on mid-line dimensions. The details are given in Ta- ble 1. In order to determine the material properties of the sections, tensile tests were set-up where the applied load and gauge specimen elongation were recorded continuously until fracture of the specimen occurred. The measured load and elongation were normalised to give a stress-strain relation- ship. Due to the anisotropy of stainless steel, a full analysis of the material properties would require tensile tests in the longitudinal and transverse directions, as well as compression tests in the same directions. Indeed, provision is made in the ASCE design code to enable use to be made of them in specific applications. However, compression tests were not carried out as there would be difficulty in establishing the true material properties of the material due to likely buckling effects. Also, transverse direction tensile tests could not be carried out because of the limitations in the geometry of the sections. Hence tensile testing was limited to the longitudinal direction. All tensile tests were carried out in accordance with BSEN10002-1 [13]. Standard tensile tests were performed to ascertain the material properties of the stainless steel. Cou- pons were cut from the webs of the columns and tested to obtain the 0.2 % proof stress and the modulus of elasticity. Tensile tests were also performed on full sections to in- clude the effects of the cold formed corners and from these tests, the 0.2% proof stress and the initial modulus of elasticity were determined. For the standard coupons, a total of three specimens were tested and the average results were noted. For the full section tests, two specimens were tested and again, the average results were noted and shown in Table 2. The results obtained for the plasticity factors n from the ASCE design code, i.e. the modified form of the Ramberg- -Osgood equation given by equation (2), are also detailed in ref. [11] and shown in Table 3. Also shown in Table 3 are the plasticity factors obtained from a comparative plots/trial and error (‘best fit’) process using the stress-strain curves obtained from the tensile tests, as reported in Ref. [11]. 5.2 Compression tests In the experimental investigation a series of compression tests to failure were made on stainless steel columns of the lipped channel cross-section as described above. The speci- men parameters investigated were as follows: column lengths varied from 222 mm to 1222 mm in increments of 100 mm (slenderness ratios varying from 42 to 234); twenty-two tests to failure (2 sets of columns, average results noted) were car- ried out with the loading applied concentric to the centroid of the cross-section to give pure axial compression; twenty-two tests to failure (2 sets of columns, average results noted) were carried out with the loading applied 8 mm eccentric to the centroid of the cross-section to give combined bending and axial compression loading. Each length of column tested was cut to the specified length and then milled flat at each end to avoid any possi- ble gripping problems. The end grips were designed such that they would hold the ends of the column and allow the loading to be applied at the required eccentricity through knife edges. The specimens were tested using a Tinius Olsen electro-mechanical testing machine, with the column vertical displacement and mid-span horizontal deflection measured during the tests using displacement transducers. Fig. 1 shows a schematic diagram of the column test configuration. 5.3 Observations Fig. 4 shows the results obtained for the load capacity of concentrically loaded stainless steel lipped channel section columns from tests, design codes (using virgin material prop- erties and full section properties) and from finite element analysis (using virgin material properties). The design code and finite element predictions show excellent correlation with the test results for all but the shortest of the columns. 94 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 3/2005 Czech Technical University in Prague P (kN) P (kN) L e=8mm e=8mm t Fig. 1: Schematic diagram of eccentrically-loaded column test Web b1 (mm) Flange b2 (mm) Lip b3 (mm) Thickness t (mm) Radius r1 (mm) Radius r2 (mm) 28.00 14.88 7.45 2.43 1.10 1.10 Table 1: Average dimensions of lipped channel cross-section Thickness t (mm) Av. Virgin 0.2 % P. S. (N/mm2) Av. Virgin UTS (N/mm2) Av. FS 0.2 % P. S. (N/mm2) Av. FS UTS (N/mm2) 2.43 480 553 520 689 Table 2: Tensile test results: virgin material and full section (FS) mechanical properties Tensile test n (ASCE) n (Best fit) Coupon 3.80 6.22 Full section 5.02 6.65 Table 3: Plasticity factors © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 95 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 3/2005 Column Length = 722 mm Cross-Section Centroid Fig. 2: Finite element plot of displacement and boundary conditions for eccentric column loading Column Length = 222 mm Fig. 3: Finite element plot of nodal stress and boundary conditions for eccentric column loading Fig. 5 shows the results obtained for the load capacity of eccentrically loaded stainless steel lipped channel section columns from tests, design codes (using virgin material prop- erties and full section properties) and from modifications to the design codes as described by equations (5) and (6). All de- sign code predictions show conservatism in prediction of load capacity for the shorter range of columns with improvements gained when full section properties are used and further im- provements are gained in using the modified forms. The best correlation obtained was for columns where the modified de- sign codes predicted accurate load capacities for all column lengths. Fig. 6 shows the results obtained for the load capacity of eccentrically loaded stainless steel lipped channel section col- umns from tests and from finite element analysis using the various stress-strain curves described earlier. The predictions of the finite element analysis show an excellent correlation to the test results and a real improvement on the predictions of 96 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 3/2005 Czech Technical University in Prague 0 10 20 30 40 50 60 70 222 322 422 522 622 722 822 922 1022 1122 1222 Column Length (mm) B u c k li n g L o a d (k N ) ASCE ASCE (FS) Euro.1.4 Euro.1.4 (FS) ANSYS Experimental Fig. 4: Graph of load capacity v. column length: concentric loading (test/design codes/finite element analysis) 0 5 10 15 20 25 30 222 322 422 522 622 722 822 922 1022 1122 1222 Column Length (mm) B u c k li n g L o a d ( k N ) ASCE ASCE (FS) ASCE (FS,Mexp) Euro.1.4 Euro.1.4 (FS) Euro.1.4 (FS,Mexp) Experimental Fig. 5: Graph of load capacity v. column length: eccentric loading (test/design codes) the design codes. However, any differences between finite ele- ment predictions and test results are very slight and occur mainly for very short length columns. 6 Conclusions The load capacity of short-to-medium length cold formed stainless steel lipped channels under pure axial compressive (concentric) loading was shown to be accurately predicted by the Eurocode 1.4 and the ASCE design code. A non-lin- ear finite element analysis using virgin material stress-strain properties also provided accurate load predictions. However, for eccentrically loaded columns with shorter lengths, the design codes predicted very conservative load ca- pacities. Improvements were gained by employing enhanced stress-strain properties and by incorporating the full section moment capacity within the code interaction formulae. Finally, it has been shown in this paper that finite element analysis can be used with a high level of confidence in predict- ing the load capacity of eccentrically loaded cold formed stainless steel, short-to-medium length columns of lipped channel section. This has been shown to be true for various stress-strain curves including that obtained from virgin mate- rial tensile tests. References [1] Johnson, A. L.: The Structural Performance of Austenitic Stainless Steel Members. Department of Structural Engi- neering, Report No 327. Cornell University, New York 1966. [2] Wang, S. T.: Cold-Rolled Austenitic Steel: Material Properties and Structural Performance. Department of Structural En- gineering, Report No. 334. Cornell University, New York 1969. [3] Wang, S. T., Errera, S. J. and Winter, G.: “Behaviour of Cold-Rolled Stainless Steel Members.” Journal of the Structural Division, Proc. ASCE 101 (ST11), 1975, p. 2337–2357. [4] Van Den Berg, G. J., Van Der Merwe, P.: “Prediction of Corner Mechanical Properties for Stainless Steels Due to Cold Forming.” Paper Presented at the 11th Interna- tional Specialty Conference on Cold Formed Steel Struc- tures, St.Louis, Missouri, (USA): October 1992. [5] Ramberg, W., Osgood, W. R.: “Description of Stress- -Strain Curves by Three Parameters.” National Advisory Committee for Aeronautics (NACA), Technical Note No. 902, February 1943. [6] Hill, H. N.: “Determination of Stress-Strain Relations from ‘Offset’ Yield Strength Values.” National Advisory Committee for Aeronautics (NACA), Technical Note No. 927, February 1944. [7] ANSI/ASCE-8-90: Specification for the Design of Cold- -Formed Stainless Steel Structural Members. 1991. [8] ENV 1993-1-3, Eurocode 3: Design of Steel Structures: Part 1.4: General Rules – Supplementary Rules for Stainless Steel. July 1996. [9] Rhodes, J., Macdonald, M., McNiff, W.: “Buckling of Cold Formed Stainless Steel Columns under Concentric and Eccentric Loading.” Proc. 15th Int. Specialty Con- ference on Cold Formed Steel Structures, St Louis, Mis- souri (USA), October 2000. [10] Rhodes, J., Macdonald, M., Kotelko,M., McNiff, W.: “Buckling of Cold Formed Stainless Steel Columns un- der Concentric and Eccentric Loading.” Proc. 3rd Int. © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ 97 Czech Technical University in Prague Acta Polytechnica Vol. 45 No. 3/2005 0 5 10 15 20 25 30 222 322 422 522 622 722 822 922 1022 1122 1222 Column Length (mm) L o a d C a p a c it y (k N ) FEA (Best Fit, n=6.22) FEA (ASCE, n=3.80) FEA (Virgin) FEA (FS) Experimental Fig. 6: Graph of load capacity v. column length: eccentric loading (test/finite element analysis) Conference on Thin Walled Structures, Krakow (Po- land), June 2001. [11] Macdonald, M., Rhodes, J., Taylor, G. T.: “Mechanical Properties of Stainless Steel Lipped Channels.” Proc. 15th Int. Specialty Conference on Cold Formed Steel Structures, St Louis, Missouri (USA), October 2000. [12] Macdonald, M.: “The Effects of Cold Forming on Mate- rial Properties and Post-Yield Behaviour of Structural Sections.” PhD Thesis. Glasgow Caledonian University, Glasgow (Scotland), January 2002. [13] British Standards Institution: Tensile Testing of Metallic Materials. BSEN10002-1. Martin Macdonald Glasgow Caledonian University School of Engineering, Science and Design Cowcaddens Road Glasgow, G4 0BA Scotland, UK Jim Rhodes e-mail: jrhodes@mecheng.strath.ac.uk University of Strathclyde Department of Mechanical Engineering 75 Montrose Street Glasgow, G1 1XJ, Scotland, UK 98 © Czech Technical University Publishing House http://ctn.cvut.cz/ap/ Acta Polytechnica Vol. 45 No. 3/2005 Czech Technical University in Prague