Bending-Induced Cross-Sectional Deformation of Cold-Formed Steel Channel-Section Beams

In practice, the cold-formed steel structures are widely used, and their geometrical nonlinearity is of great significance in structural design. However, the current analytical solutions are only applied to linear problems, and the nonlinear problems are mostly considered by commercial finite element software.)is paper presents an analytical model to describe the bending-induced crosssectional deformation in cold-formed steel channel-section beams when subjected to transverse bending loads. )e nonlinear effect of cross-sectional deformation on the relationship between the overall bending curvature and applied moment is examined. A parametric study is also carried out to examine how the flange width, lip length, and thickness of the section affect the crosssectional deformation and corresponding nonlinear bending behavior of the beams.


Introduction
in-walled beams frequently exhibit flattened cross-sectional deformation when they are subjected to bending loads. A typical example is the ovalization of the cross-section of circular cylindrical shells of infinite length under pure bending, which was reported by Brazier in 1927 [1]. e cross-sectional flattening leads to a reduction of the overall bending stiffness of the shells, resulting in a nonlinear relationship between the curvature and applied moment. For the circular cylindrical shells of finite length, the ovalization of the cross-section is also affected by the length of the shells and the restraint conditions at the two ends of the shells [2]. e similar phenomenon of the flattened cross-sectional deformation has been also found in angle-section beams [3]; channel-section beams bent about the axis parallel to web line [4] and rectangular box-section beams [5], in which the bending-induced cross-sectional deformation leads to a reduction in the bending rigidity of the beams and thus results in a nonlinear relationship between the bending deflection and applied transverse load. e flattened cross-sectional deformation observed in hat-section beams when subjected to bending is referred to as flange curling, which, in principle, is similar to the ovalization of the cross-section of circular cylindrical shells. When the flange curling is considered, the crosssection properties of the beam need to be calculated based on the deformed cross-section with curved geometry, representing the second-order effect. e phenomenon of flange curling was first reported by Winter in 1940 [6]. In Eurocode 3-Part 1.3 a formula is suggested to estimate the deformation caused by the flange curling [7]. Experimental investigation has been carried out on the flange curling in profiled steel decks with different slender geometries [8].
e results of the experiments were compared with those obtained using nonlinear finite element analysis method. Analytical models were also developed to study the flange curling phenomenon in wide flange trapezoidal panels [9] and wide single-flange steel panels [10].
e analytical calculation methods were demonstrated using experimentally obtained data.
Compared to hot-rolled steel beams, cold-formed steel beams may be even easier to exhibit cross-sectional deformation because of its thin thickness and open-section nature. To predict the effect of the cross-sectional deformation on the overall bending behavior of the beams, one has to know how individual segments deform within the cross-section. In this paper, an analytical model is presented to describe the bending-induced flange deformation in coldformed steel channel-section beams when subjected to transverse bending loads. e nonlinear effect of the local flange deformation within the cross-section on the relationship between the overall bending curvature and applied moment is examined. A parametric study is also carried out to examine how the flange width, lip length, and thickness of the section affect the cross-section deformation and corresponding nonlinear bending behavior of the beam.

Analytical Model Describing Flange Curling in Channel-Section Beams
When a beam is subjected to a pure bending, the bending stresses of the beam on the top and bottom layers are in opposite directions (see the stresses presented by dash lines in Figure 1). Beam curvature causes a vertical imbalance of the longitudinal tension and compression bending stresses.
As it is shown in Figure 1, the combination of the two compressive stresses on the top layer and the two tensile stresses on the bottom layer leads to a vertical compression of the cross-section, which causes an in-plane flattening of the cross-section. For solid-section beams, the deformation generated by this vertical compression is negligible and thus there is no need to consider the flattening of the crosssection. However, for thin-walled open-or close-section beams, the deformation generated by the vertical compression may be significant and thus there is a need to consider the cross-sectional shape changing and/or the flattening deformation of the cross-section. Herein, we consider a thin-walled, cold-formed steel channel-section beam with web depth h, flange width b, lip length c, and thickness t. Assume that the beam is subjected to a bending load and results in deflection v (x) and bending compressive stress σ (x) in its upper flange, where x is the longitudinal coordinate of the beam. e equivalent compressive forces in the upper flange and lip generated by the bending compressive stress thus are σbt and σct (1 − c/h), respectively. For a beam element of length dx, the corresponding flattening force is σbtv"dx, which can be assumed to distribute uniformly on the upper flange over a surface area bdx, and σct (1 − c/h) v"dx, which acts on the lip over a line length dx. Hence, the former provides a surface load of σtv" on the upper flange and the latter gives a line load of σct (1− c/h) v" on the joint line between the upper flange and lip. e similar flattening forces can be also derived for the lower flange and lip, which are shown in Figure 2(a). e deformation of the cross-section of the beam under the action of the flattening forces can be calculated using a frame model, as shown in Figure 2(b), in which the web line is subjected a pure bending, and the upper and lower flanges are subjected to a uniformly distributed load plus a concentrated load acting at the lip. Based on the theory of plates and shells [11], the local deflection of the flange thus can be expressed as follows: where D � Et3121-]2 is the bending rigidity and z is the coordinate along with flange line, E is the Young's modulus, and ] is the Poisson's ratio. e first term in equation (1) is the deflection due to the web rotation; the 2 nd and 3 rd terms in equation (1) are caused by the uniform load on the flange and concentrated load on the lip, respectively [12]. Based on the theory of plates and shells [11], the deflection and rotation of the lip are expressed as follows: (2) e second moment of area of the deformed channelsection beam about its neutral axis thus can be calculated as follows:

Advances in Civil Engineering
For small local deformation, the second-order terms of y (x, z) would be small and thus can be neglected. In this case, the reduction of the second moment of area caused by crosssectional flattening can be approximated as follows: Substituting equations (1) and (2) into (5), it yields Note that the bending stress in the flanges can be expressed as Substituting equation (6) into equation (5), it yields where ∆I o is a section dimension-dependent parameter and is expressed as follows: Bending stresses

Flattening stresses
Flattening stresses Bending stresses e bending equation of the channel-section beam established on the deformed framework can be expressed as follows: where I o is the second moment of area of the undeformed channel-section beam about its major axis and M (x) is the internal moment of the beam. Since ∆I is the quadratic function of curvature, equation (9) is a nonlinear differential equation about the curvature. To solve these nonlinear differential equations, one has to use an approximate method or a numerical method. Assume the solution of equation (9) can be expressed as the sum of two parts, that is, where v o (x) is the deflection function without considering flange curling effect and v 1 (x) is the additional deflection function when the flange curling is taken into account. It is obvious that v o (x) can be determined from the classical bending theory of beams as follows: Substituting equation (10) into equation (9) and using equation (11) and noticing the fact that v 1 In equation (12), ∆I has been evaluated by using v o (x) only. Equation (12) indicates that the additional deflection function v 1 (x) is generated due to an equivalent moment expressed by the right-hand-side term of equation (12), which is proportional to (d 2 v o /dx 2 ) 3 . Hence, to find the solution of v(x), we first find v o (x) from equation (11), then we approximate ∆I (x) using equation (7) and then the additional moment of Eh 2 (d 2 v o /dx 2 ) 3 ΔI o , and finally, we calculate v 1 (x) from equation (12).

Numerical Examples
Two examples are given herein to demonstrate the model developed above. One is the beam subjected to a pure bending, and the other is the beam subjected to a uniformly distributed load at shear centre. In both cases, the beam is assumed to be simply supported at its two ends and have the dimensions of beam length l � 7000 mm, web depth e deflection function obtained by solving equation (12) can be expressed as follows: Substituting equation (13) into equation (7), it yields Substituting equation (13) into equation (12) and solving the differential equation, it yields Hence, the nonlinear deflection function of the beam when the flange curling is taken into account is expressed as Figure 3 shows the variation of the deflection of the channel-section beam at mid-span with the applied moment, calculated using equations (13) and (16), in which the three straight lines and three curves represent the beams with three different flange widths. It can be seen from the figure that the three straight lines, obtained from the linear solution equation (13), coincide together for the three beams of different flange widths. is is because of the use of the dimensionless moment, M o /M y , in the plot. However, the three curves, obtained from the nonlinear solution equation (16), for the three beams of different flange widths, are different. For the same moment ratio, the wider the flange, the larger the deflection, indicating that the flange curling is more significant in the beams with wider flanges. Figure 4 shows the similar moment-deflection curves of the channelsection beams with three different section thicknesses. It is observed from the figure that the three linear solutions coincide together, whereas the three nonlinear solutions are slightly different. e thicker the section thickness is, the closer the corresponding nonlinear solution is to the linear solution, indicating the flange curling is less important in the beams with thicker section thickness. Figure 5 shows the moment-deflection curves of the channel-section beams with three different lip lengths. It can be seen from the figure that the lip length has a considerable influence on flange curling of channel-section beam. e effects of nonlinearity are more significant as the lip length increases. Similar to those shown in Figures 3 and 4, all three linear solutions come together, but three nonlinear solutions are dispersed. e longer the lip length is, the more the nonlinear solution is different from the linear solution, indicating that the flange curling is more n the beams with longer lip lengths. e deflection function v o (x) can be obtained by solving equation (11) and is expressed as follows: Substituting equation (17) into equation (7), it yields Substituting equation (18) into equation (12) and solving the differential equation, it yields

Advances in Civil Engineering
Hence, the nonlinear deflection function of the beam when the flange curling is taken into account is expressed as e maximum deflection at the mid of the beam, where x � l/2, can be expressed as follows: Comparing equation (21) with equation (16), it is obvious that the flange curling in the beam with a uniformly distributed load is very similar to that in the beam under pure bending. Only difference is the precoefficient of the nonlinear term, which is 279/350 instead of 1.0 in the beam under the pure bending, indicating that the flange curling in the former would be slightly larger than that in the latter.
Furthermore, the nonlinear finite element analysis method [13][14][15] is used to demonstrate the present analytical model. e finite element model of the channel-section beam under pure bending is established using ANSYS software. Considering the beam is symmetric about its midsection, only half of the beam length is modelled in the finite element analysis. Figure 6 shows the finite element mesh of the half beam and corresponding boundary conditions employed in the analysis. e analysis is performed using 4noded thin shell elements with a maximum element size not exceeding 10 mm, the choice of which is based on our previous study on the similar analysis [16][17][18]. e boundary conditions are assumed as that all nodes on the symmetric section (x � l/2) have zero axial displacement and zero rotations about lateral and transverse axes, and all nodes on the end section (x � 0) have zero lateral and transverse displacements and zero rotation about longitudinal axis. To apply the pure bending, a fictitious rigid plate is used to cover the end section, on which a concentrated moment is applied. To avoid the lateral-torsional buckling of the beam, the lateral displacement of the upper and lower flanges is restrained in the analysis. Linear elastic material model is used in the analysis with Young's modulus E � 200 GPa and Poisson's ration ] � 0.3. In other words, it is only the geometric nonlinearity that is considered in the analysis. Figure 7 shows the comparison of the moment-deflection curves obtained from the nonlinear finite element analysis and presents analytical model for two channelsection beams. It can be seen from the figure that the moment-deflection curves provided by the present model are very close to those predicted by the nonlinear finite element analysis, indicating that the present approach for modelling flange curling is adequate, although it is very simple.

Conclusions
Cold-formed steel beams have been increasingly used in buildings and other load-carrying structures. Due to thin thickness cold-formed steel beams may exhibit cross-sectional deformation when they are subjected to bending loadings, the interaction between the transverse bending of the beam and the in-plane flattening of the beam crosssection leads the beams to have a nonlinear bending behavior. In this paper, an analytical model has been developed to describe the interactive effect between the bending and flattening for lipped channel-section beams. e present model has been demonstrated by using nonlinear finite element analysis results. From the present study, the following conclusions can be drawn: (i) e cross-sectional deformation in channel-section beams when subjected to transverse bending loadings is mainly caused by the flange curling, which results in a reduction of the overall flexural rigidity of the beams.
(ii) e remarkable flange curling mainly occurs in the channel-section beams, which have wide flanges, long lips, and thin thickness, and the greater flange curling occurs with higher bending stresses. (iii) e magnitude of flange curling occurred in lipped channel-section beams increases with the increase of flange width or lip length, but decreases with increased section thickness. (iv) e flange curling results in an increase of the bending deflection of the beam and the nonlinear deflection may be significantly different from that calculated using the linear bending theory of beams.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
No potential conflicts of interest were reported by the authors. Advances in Civil Engineering 7