Web Crippling Behaviour of Cold-Formed High Strength Steel Unlipped Channel Beams

: Cold-formed sections (CFS) fabricated using high strength steel have recently been utilised in construction due to their numerous advantages, such as higher load-to-weight ratio, ﬂexibility of shape, and availability in relatively long spans. High strength CFS channel sections can be used as purlins and joists in structural systems; thus, they are vulnerable to different buckling instabilities, including web crippling. Predicting their web crippling capacity using the current design guidelines may be insufﬁcient due to their empirical nature. This study, therefore, aims to investigate the web crippling capacity of high strength unlipped CFS sections under End-Two-Flange (ETF) loading conditions. Numerical simulations were carried out using nonlinear ﬁnite element (FE) analysis. The developed models were ﬁrst validated against available experimental data and then used as a base for conducting an extensive parametric study. The ultimate web crippling capacity obtained from the parametric study was used to assess the accuracy of the available design equations in the standards and those proposed in the relevant studies. The assessment revealed that the existing design equations are not suitable for predicting the ultimate web crippling capacity for high strength CFS channel sections under the ETF loading condition. Thus, a modiﬁed design equation was proposed, following the same technique of current design standards, and a new Direct Strength Method (DSM) approach was developed.


Introduction
Cold-formed steel members have been extensively used in structural frame systems due to their prominent advantages, including light weight, durability, stability, high accuracy in fabrication, efficient erection, high strength-to-weight ratio and high stiffness.Different shapes of CFS sections, such as LiteSteel Beam (LSB) and built-up sections, as well as conventional sections (Z-section, C-section, ∑-section and track sections), have been developed to meet construction needs [1].High strength steel is one of the great additions to the construction industry; it is capable of fulfilling structural design needs such as longer spans, higher strength with reduced construction mass, higher load-carrying capacity and reduced self-weight in high-rise buildings [2].The structural performance of high strength steel sections has been studied by various researchers [2][3][4][5].Web crippling is a vulnerability of CFS members under concentrated reactions and concentrated transverse reactions due to their thin cross sections.Since the 1940s, several experimental investigations of the web crippling behaviour of cold-formed steel sections [6][7][8][9][10][11] have been conducted by researchers, and web crippling design equations and standards have been adopted in design specifications such as Eurocode 3 Part 1-3 [12], AISI S100 [13] and AS/NZS 4600 [14].In addition, the design specifications assign the failure modes of web crippling to four categories-Interior-Two-Flange (ITF), End-Two-Flange (ETF), End-One-Flange (EOF) and Interior-One-Flange (IOF)-according to the locations of the loading, supporting and failure region [13].The web crippling load cases are illustrated in Figure 1.
Buildings 2022, 12, x FOR PEER REVIEW 2 of 29 crippling is a vulnerability of CFS members under concentrated reactions and concentrated transverse reactions due to their thin cross sections.Since the 1940s, several experimental investigations of the web crippling behaviour of cold-formed steel sections [6][7][8][9][10][11] have been conducted by researchers, and web crippling design equations and standards have been adopted in design specifications such as Eurocode 3 Part 1-3 [12], AISI S100 [13] and AS/NZS 4600 [14].In addition, the design specifications assign the failure modes of web crippling to four categories-Interior-Two-Flange (ITF), End-Two-Flange (ETF), End-One-Flange (EOF) and Interior-One-Flange (IOF)-according to the locations of the loading, supporting and failure region [13].The web crippling load cases are illustrated in Figure 1.
The web crippling behaviour of CFS channel sections under all four load cases, considering different supporting conditions (fastened support to the flanges or unfastened support to the flanges), has been investigated, and improvements in the codified design equations have been proposed .Moreover, web crippling studies have been conducted on various sections, such as rectangular hollow flange sections [36][37][38][39] and SupaCee sections [17,40], to check the applicability of the codified design equations available in the standards, whilst some studies have investigated the web crippling behaviour of various sections with web openings [41][42][43][44][45].However, all these findings are limited to a yield strength less than 500 MPa, due to the unavailability of cold-forming technology [5].High strength steel sections are receiving attention from the construction community due to their added benefits in terms of structural and architectural behaviour [5].The researchers Santaputra et al. [3] and Wu et al. [4] studied the web crippling of high strength steel, including brake-pressed hat and built-up sections and hat and deck sections, respectively.The study of Santaputra et al. [3] covered a steel section strength of 401 to 1138 MPa, whilst Wu et al. [4] considered the virgin material yield stresses parallel to the rolling The web crippling behaviour of CFS channel sections under all four load cases, considering different supporting conditions (fastened support to the flanges or unfastened support to the flanges), has been investigated, and improvements in the codified design equations have been proposed .Moreover, web crippling studies have been conducted on various sections, such as rectangular hollow flange sections [36][37][38][39] and SupaCee sections [17,40], to check the applicability of the codified design equations available in the standards, whilst some studies have investigated the web crippling behaviour of various sections with web openings [41][42][43][44][45].However, all these findings are limited to a yield strength less than 500 MPa, due to the unavailability of cold-forming technology [5].High strength steel sections are receiving attention from the construction community due to their added benefits in terms of structural and architectural behaviour [5].The researchers Santaputra et al. [3] and Wu et al. [4] studied the web crippling of high strength steel, including brake-pressed hat and built-up sections and hat and deck sections, respectively.The study of Santaputra et al. [3] covered a steel section strength of 401 to 1138 MPa, whilst Wu et al. [4] considered the virgin material yield stresses parallel to the rolling directions, ranging from 716 to 776 MPa, and perpendicular to the rolling directions, ranging from 825 to 894 MPa.Furthermore, web crippling investigation on high strength tubular sections (square and rectangular) under all four load cases was carried out by Li and Young [2,5].
The researchers considered the yield strengths of 700 MPa and 900 MPa in their study [2,5].Sundarajah et al. [18] investigated the web crippling behaviour of high strength steel lipped channel sections under EOF and IOF loading conditions.The study also limited the investigation to less than 550 MPa.
This paper investigates the web crippling behaviour of unlipped channel sections under the ETF load case, with high strength material (yield strength of 700 MPa, 900 MPa and 1000 MPa).In total, 243 numerical studies, following the validation study, were conducted to investigate the web crippling behaviour and evaluate the applicability of current design equations available in standards.The test specimen length and loading conditions were chosen as denoted in the design specifications [12][13][14].Web crippling strengths for this study were compared to the predictions of codified design equations provided in the design standards [12][13][14] and with previous studies on cold-formed unlipped channel sections [21,22].Finally, this paper proposes a modified design equation for the web crippling capacity of CFS lipped channel sections with high strength under the ETF load case based on the current design standards.A DSM design equation was also proposed for this study.

Finite Element Modelling Description
The web crippling behaviour of channel sections can be predicted using Finite Element (FE) software and using experimental studies available in the literature [17,[19][20][21][22][23][24][25][26][27]36,40].The experimental set-up of channel sections as shown in Figure 2 is simulated using ABAQUS version 6.14 (developed by Dassault Systèmes Simulia Corp., Johnston, RI, USA) [46] to extend the study of the web crippling of unlipped channel sections with high strength material.In this study, full and then half models with axial symmetry assignment along the symmetrical line were analysed during the validation (Figure 3).Similar behaviours in both models were obtained and are discussed in Section 3 in detail.The number of elements and Degree of Freedoms (DoFs) for the half models are less than that of the full model, and thus reduced analysis time was observed for the half model as compared to the full model.In addition, a simplified modelling process was addressed during the half model development.Hence, the half model set-up was adopted during the parametric studies and is described herein.
directions, ranging from 716 to 776 MPa, and perpendicular to the rolling directions, ranging from 825 to 894 MPa.Furthermore, web crippling investigation on high strength tubular sections (square and rectangular) under all four load cases was carried out by Li and Young [2,5].The researchers considered the yield strengths of 700 MPa and 900 MPa in their study [2,5].Sundarajah et al. [18] investigated the web crippling behaviour of high strength steel lipped channel sections under EOF and IOF loading conditions.The study also limited the investigation to less than 550 MPa.
This paper investigates the web crippling behaviour of unlipped channel sections under the ETF load case, with high strength material (yield strength of 700 MPa, 900 MPa and 1000 MPa).In total, 243 numerical studies, following the validation study, were conducted to investigate the web crippling behaviour and evaluate the applicability of current design equations available in standards.The test specimen length and loading conditions were chosen as denoted in the design specifications [12][13][14].Web crippling strengths for this study were compared to the predictions of codified design equations provided in the design standards [12][13][14] and with previous studies on cold-formed unlipped channel sections [21,22].Finally, this paper proposes a modified design equation for the web crippling capacity of CFS lipped channel sections with high strength under the ETF load case based on the current design standards.A DSM design equation was also proposed for this study.

Finite Element Modelling Description
The web crippling behaviour of channel sections can be predicted using Finite Element (FE) software and using experimental studies available in the literature [17,[19][20][21][22][23][24][25][26][27]36,40].The experimental set-up of channel sections as shown in Figure 2 is simulated using ABAQUS version 6.14 (developed by Dassault Systèmes Simulia Corp., Johnston, RI, USA) [46] to extend the study of the web crippling of unlipped channel sections with high strength material.In this study, full and then half models with axial symmetry assignment along the symmetrical line were analysed during the validation (Figure 3).Similar behaviours in both models were obtained and are discussed in Section 3 in detail.The number of elements and Degree of Freedoms (DoFs) for the half models are less than that of the full model, and thus reduced analysis time was observed for the half model as compared to the full model.In addition, a simplified modelling process was addressed during the half model development.Hence, the half model set-up was adopted during the parametric studies and is described herein.

Element Type and Mesh
The FE model was developed with two major parts: the cold-formed channel beams and bearing plate.The 3D deformable shell element type of S4R was assigned to the coldformed steel beam, whilst rigid element R3D4 was used to create the bearing plates.Although various types of elements are available in the ABAQUS, the above-mentioned element types were chosen according to the various previous numerical studies conducted on web crippling [2,16,17,19,[25][26][27].Since the shell offset definition was set to be the middle surface, the centre-line dimensions were simulated in accordance with the ABAQUS.This means the middle surface would act as a reference surface, and an assigned thickness will be added to both sides of the reference surface with half of the thickness [46].Mesh size is rather crucial factor in deciding the accuracy of numerical results.Mesh size of the cold-formed steel beam was adopted in this study as indicated in the previous web crippling studies based on mesh refinement analysis [16][17][18][19].In this study, the mesh size of the flat region and corner regions was 5 mm × 5 mm and 5 mm × 1 mm, respectively.The applied mesh sizes and element types are shown in Figure 4. Plastic deformation at the web-flange juncture and appropriate load transfer to the web from the flange is significantly influenced by the corner region, and thus finer mesh is recommended for the corner radii in the web crippling studies [16][17][18][19].Since the rigid plate does not take part in the finite element analyses, the mesh size of the bearing plate was assigned as 10 mm × 10 mm.

Material Modelling
The two-stage Ramberg-Osgood model is used to describe the stress-strain behaviour of the stainless steels and cold-formed steels.Even though a number of researchers have investigated the behaviour of stainless steel [47][48][49][50][51], only a few researchers have studied the material modelling of cold-formed steel [52,53].Recently, Leroy Gardner and Xiang Yun [54] proposed a two-stage Ramberg-Osgood model (Equations ( 1)-( 6)) to describe the stress-strain curve of the cold-formed steel, which is based upon the material model developed by Mirambell and Real [51] for stainless steel.The researchers [54] gathered many stress-strain models from a literature review by considering various material strengths, shapes and forms of production, such as press-braked or cold-rolled and validated.In addition, the study [54] suggested that the study of Rossi et al. [55] can be used to obtain the material model of the corner regions (Equations ( 7)-( 10)).The material models proposed by Gardner and Yun [54] and Rossi et al. [55] were applied in this study to determine the stress-strain behaviour of the flat and corner regions, respectively.

Element Type and Mesh
The FE model was developed with two major parts: the cold-formed channel beams and bearing plate.The 3D deformable shell element type of S4R was assigned to the cold-formed steel beam, whilst rigid element R3D4 was used to create the bearing plates.Although various types of elements are available in the ABAQUS, the above-mentioned element types were chosen according to the various previous numerical studies conducted on web crippling [2,16,17,19,[25][26][27].Since the shell offset definition was set to be the middle surface, the centre-line dimensions were simulated in accordance with the ABAQUS.This means the middle surface would act as a reference surface, and an assigned thickness will be added to both sides of the reference surface with half of the thickness [46].Mesh size is rather crucial factor in deciding the accuracy of numerical results.Mesh size of the cold-formed steel beam was adopted in this study as indicated in the previous web crippling studies based on mesh refinement analysis [16][17][18][19].In this study, the mesh size of the flat region and corner regions was 5 mm × 5 mm and 5 mm × 1 mm, respectively.The applied mesh sizes and element types are shown in Figure 4. Plastic deformation at the web-flange juncture and appropriate load transfer to the web from the flange is significantly influenced by the corner region, and thus finer mesh is recommended for the corner radii in the web crippling studies [16][17][18][19].Since the rigid plate does not take part in the finite element analyses, the mesh size of the bearing plate was assigned as 10 mm × 10 mm.

Material Modelling
The two-stage Ramberg-Osgood model is used to describe the stress-strain behaviour of the stainless steels and cold-formed steels.Even though a number of researchers have investigated the behaviour of stainless steel [47][48][49][50][51], only a few researchers have studied the material modelling of cold-formed steel [52,53].Recently, Leroy Gardner and Xiang Yun [54] proposed a two-stage Ramberg-Osgood model (Equations ( 1)-( 6)) to describe the stressstrain curve of the cold-formed steel, which is based upon the material model developed by Mirambell and Real [51] for stainless steel.The researchers [54] gathered many stress-strain models from a literature review by considering various material strengths, shapes and forms of production, such as press-braked or cold-rolled and validated.In addition, the study [54] suggested that the study of Rossi et al. [55] can be used to obtain the material model of the corner regions (Equations ( 7)- (10)).The material models proposed by Gardner and Yun [54] and Rossi et al. [55] were applied in this study to determine the stress-strain behaviour of the flat and corner regions, respectively.The applied stress-strain relationship of the flat region in this study is as follows: The suggested values and expressions for the key input parameters in Equation (1) are as follows: n-First strain-hardening exponent: In this study, the first strain-hardening exponent was taken from following table (Table 1) [55].The applied stress-strain relationship of the flat region in this study is as follows: The suggested values and expressions for the key input parameters in Equation ( 1) are as follows: n-First strain-hardening exponent: n = ln(4) ln f y /σ 0.05 (2) In this study, the first strain-hardening exponent was taken from following table (Table 1) [55].
ε u -Strain corresponding to the ultimate tensile strength f u : For an unknown ultimate strength, f u : m-Second strain-hardening exponent: f y,c -Enhanced yield strength of corner regions of cold-rolled or press-braked sections: ε c,av -Average plastic strain induced during the forming of the corner regions of press-braked and cold-rolled sections: p and q-Material parameters calculated directly from the basic properties of the unformed material: Figure 5 shows the applied material stress-strain curve of the flat and corner regions for a yield strength of 700 MPa.The derived engineering stress-strain curve is then converted to a true stress-strain curve, as stated in the ABAQUS [46].Then, a weight density of 7850 kg/m 3 , Young's modulus of 203,000 MPa and Poisson's ratio of 0.3 were applied as the material properties.

Loading and Boundary Conditions
Since full and half models were developed during the validation, different boundary conditions were applied to the two models, as shown in Figure 6.The boundary conditions respective of loading and supporting were applied to the bearing plate via a reference point on the plates, whilst boundary conditions were employed for the symmetry line of the channel section to replicate the symmetrical type used for the cold-formed steel beam in the half models.Translation of the loading plate along the X-axis and Z-axis and rotation about the Y-axis and Z-axis were restrained, whilst a −20 mm displacement was applied in the vertical direction to simulate the experimental set-up of web crippling studies.The supporting plate in the full model was restrained with all movements, except rotation about the X-axis.As half hinges are used in the experimental set-up, loading and supporting plates are allowed to rotate about the X-axis.Vertical displacement at the loading plate was controlled using an amplitude model of smooth step, which ensures the smooth deformation from the initial stage.Various types of boundary conditions are available in the ABAQUS, and the symmetrical (symmetry/antisymmetry/encastre) type of boundary condition was applied along the symmetrical line of the cold-formed beam.The translation along the Y-axis and rotation about the X-axis and Z-axis were restrained along the symmetrical line.Previous studies also used similar boundary conditions for web crippling studies [2,16,17,19,[25][26][27]31,[36][37][38].
smooth deformation from the initial stage.Various types of boundary conditions are available in the ABAQUS, and the symmetrical (symmetry/antisymmetry/encastre) type of boundary condition was applied along the symmetrical line of the cold-formed beam.The translation along the Y-axis and rotation about the X-axis and Z-axis were restrained along the symmetrical line.Previous studies also used similar boundary conditions for web crippling studies [2,16,17,19,[25][26][27]31,[36][37][38].

Contact Formulation and Properties
Interactions between bearing plates and the cold-formed steel beam were assigned by generating the surface-to-surface contacts with contact properties.Even though contact formulation can be described using various methods, such as general contact, contact pair and surface-to-surface contact, the surface-to-surface contact method was employed in this study in accordance with previous research [2,17,19,26,27].Rigid bearing plates were defined as master surfaces, whereas the flange of the beam and corner radii were specified smooth deformation from the initial stage.Various types of boundary conditions are available in the ABAQUS, and the symmetrical (symmetry/antisymmetry/encastre) type of boundary condition was applied along the symmetrical line of the cold-formed beam.The translation along the Y-axis and rotation about the X-axis and Z-axis were restrained along the symmetrical line.Previous studies also used similar boundary conditions for web crippling studies [2,16,17,19,[25][26][27]31,[36][37][38].

Contact Formulation and Properties
Interactions between bearing plates and the cold-formed steel beam were assigned by generating the surface-to-surface contacts with contact properties.Even though contact formulation can be described using various methods, such as general contact, contact pair and surface-to-surface contact, the surface-to-surface contact method was employed in this study in accordance with previous research [2,17,19,26,27].Rigid bearing plates were defined as master surfaces, whereas the flange of the beam and corner radii were specified

Contact Formulation and Properties
Interactions between bearing plates and the cold-formed steel beam were assigned by generating the surface-to-surface contacts with contact properties.Even though contact formulation can be described using various methods, such as general contact, contact pair and surface-to-surface contact, the surface-to-surface contact method was employed in this study in accordance with previous research [2,17,19,26,27].Rigid bearing plates were defined as master surfaces, whereas the flange of the beam and corner radii were specified as slave surfaces during the contact formulation.Normal and tangential behaviours were also assigned in the contact formulation.For the contact properties, a pressure overclosure of "hard" contact in the normal direction and a friction formulation of a penalty with a 0.4 friction coefficient were employed [21,25,36].

Analysis Method
ABAQUS/Implicit and ABAQUS/Explicit are capable of solving the static and quasistatic problems in the ABAQUS [46].However, due to the convergence issues that occurred during the application of the ABAQUS/Implicit analysis method, ABAQUS/Explicit was employed in this study to simulate the web crippling behaviour of the channel section.Previous studies also investigated the application of various analysis methods and revealed that the explicit method requires a huge number of increments; thus, the method can be more efficient and is more suitable for complex nonlinear problems as well as convergence issues generated by the other solver options [19].

Initial Geometric Imperfections
The effect of initial geometric imperfection has been neglected in various research studies, since it does not greatly affect the ultimate web crippling capacity [19,26].Therefore, the effect of geometrical imperfection was not considered in this study.

Residual Stresses
Residual stresses can be categorized into flexural/bending and membrane residual stresses.However, previous studies have stated that the effects of membrane and flexural residual stresses on the web crippling capacity are less than 0.5% [25]; therefore, the residual stress effect was not considered in this study.

Validation of Numerical Models
Web crippling experimental studies on tubular sections with high strength cold-formed steel [2,5] and unlipped channel sections with yield strengths of 500 MPa to 600 MPa [21] under ETF loading conditions were validated in this study to ensure the accuracy of the numerical model.The developed FE models were verified by comparing the (i) web crippling failure loads, (ii) failure modes and (iii) applied load vs. vertical deflection curves with corresponding experimental results [2,5,21].Since half models were used in this study during the parametric studies, both the full and half models generated and validated are presented here.A good agreement was obtained for the models, and less analysis time was required for the half models compared to the full models.Thus, the authors used the half models during the parametric studies.
A comparison based on the web crippling failure load (ultimate load) of the ten experimental studies was carried out for both the full and half models and is presented in Tables 2 and 3 for tubular sections and unlipped channel sections, respectively.The validation was conducted for a material strength of 700 MPa and 900 MPa for tubular sections, whilst for unlipped channel sections, a strength of 500 MPa to 600 MPa was used.The mean and COV values for the ratio between the experimental and FE results were also observed, as shown in Tables 2 and 3.The mean value and COV of tubular sections with high strength material for the full model were 0.94 and 0.05, and for the half model they were 0.95 and 0.05, respectively.Meanwhile the mean and COV of unlipped channel sections for full models were 1.00 and 0.02, and for half models they were 1.00 and 0.03.A summary of the validation results is given in Table 4.The comparison indicated that developed FE models were capable of predicting well the experimental study, with less than 10% deviation under the ETF loading condition in terms of failure load.The numerical models were validated based on the load vs. vertical deformation curve for both tubular and unlipped channel sections [2,21].The load vs. vertical deformation curve for tubular section ETF H120X120X4 with a bearing plate length of 60 mm was plotted against the load-vertical deformation curve of numerical models, whilst three experimental studies of unlipped channel sections, including ETF-UC10012 with 50 mm bearing length, ETF-UC15015 with 100 mm bearing length and ETF-UC20015 with 150 mm bearing length, were compared for unlipped channel sections.The curves were also compared for both the full and half models to ensure the usage of half models in the parametric study.The obtained plots are shown in Figures 7-10; FE models represented the experimental loadvertical deformation curve well.mm bearing length, were compared for unlipped channel sections.The curves were also compared for both the full and half models to ensure the usage of half models in the parametric study.The obtained plots are shown in Figures 7-10; FE models represented the experimental load-vertical deformation curve well.mm bearing length, were compared for unlipped channel sections.The curves were also compared for both the full and half models to ensure the usage of half models in the parametric study.The obtained plots are shown in Figures 7-10; FE models represented the experimental load-vertical deformation curve well.Finally, the validation was performed by comparing the failure modes.The failure mode of the tubular section (ETF H120X120X4 with a bearing plate length of 120 mm) and unlipped channel sections (ETF-UC10012 with 50 mm bearing plate, ETF-UC15015 with 100 mm bearing plate, ETF-UC20015 with 150 mm bearing plate) for both the full and half models was compared.A good agreement was observed and is illustrated in Figures 11-14.The comparison shows that the failure modes of FE models match well with their experimental counterparts.Finally, the validation was performed by comparing the failure modes.The failure mode of the tubular section (ETF H120X120X4 with a bearing plate length of 120 mm) and unlipped channel sections (ETF-UC10012 with 50 mm bearing plate, ETF-UC15015 with 100 mm bearing plate, ETF-UC20015 with 150 mm bearing plate) for both the full and half models was compared.A good agreement was observed and is illustrated in Figures 11-14.The comparison shows that the failure modes of FE models match well with their experimental counterparts.Finally, the validation was performed by comparing the failure modes.The failure mode of the tubular section (ETF H120X120X4 with a bearing plate length of 120 mm) and unlipped channel sections (ETF-UC10012 with 50 mm bearing plate, ETF-UC15015 with 100 mm bearing plate, ETF-UC20015 with 150 mm bearing plate) for both the full and half models was compared.A good agreement was observed and is illustrated in Figures 11-14.The comparison shows that the failure modes of FE models match well with their experimental counterparts.The comparison of FE models with experimental studies showed that the numerical model for both the full and half models is capable of predicting the web crippling capacity of cold-formed beams with high strength material under the ETF loading condition.However, reduced analysis time was observed for the half model in comparison to that of the full model, and thus half models were used to perform the parametric study of the cold-  The comparison of FE models with experimental studies showed that the numerical model for both the full and half models is capable of predicting the web crippling capacity of cold-formed beams with high strength material under the ETF loading condition.However, reduced analysis time was observed for the half model in comparison to that of the full model, and thus half models were used to perform the parametric study of the cold- The comparison of FE models with experimental studies showed that the numerical model for both the full and half models is capable of predicting the web crippling capacity of cold-formed beams with high strength material under the ETF loading condition.However, reduced analysis time was observed for the half model in comparison to that of the full model, and thus half models were used to perform the parametric study of the cold-formed channel section with high strength undergoing web crippling under the ETF loading condition.Since accurate results were obtained in the validation, the parametric study was extended using numerical analyses.

Parametric Study
The validated FE models were incorporated in a comprehensive parametric study to generate a database of high strength steel unlipped channel sections subject to ETF loading to investigate their web crippling behaviour.The geometric and material parameters of section depth (d), width (b f ), thickness (t), corner radius (r i ), bearing length (N) and steel grade (f y ) were varied during the study, covering a range of section slenderness values.Figure 15 shows the notations used for key cross-sectional dimensions of unlipped channel sections.The available web crippling design rules were then assessed, and improved design provisions were proposed to increase the accuracy of web crippling resistance predictions in succeeding sections using the generated database of unlipped channel sections.
Buildings 2022, 12, x FOR PEER REVIEW 14 of 29 formed channel section with high strength undergoing web crippling under the ETF loading condition.Since accurate results were obtained in the validation, the parametric study was extended using numerical analyses.

Parametric Study
The validated FE models were incorporated in a comprehensive parametric study to generate a database of high strength steel unlipped channel sections subject to ETF loading to investigate their web crippling behaviour.The geometric and material parameters of section depth (d), width (bf), thickness (t), corner radius (ri), bearing length (N) and steel grade (fy) were varied during the study, covering a range of section slenderness values.Figure 15 shows the notations used for key cross-sectional dimensions of unlipped channel sections.The available web crippling design rules were then assessed, and improved design provisions were proposed to increase the accuracy of web crippling resistance predictions in succeeding sections using the generated database of unlipped channel sections.Table 5 summarises the key parameters that were considered during the parametric study.Since high strength steel sections have been applied in various areas [56,57], the material strength was chosen by considering those applications and by focusing on the potential applications of high strength channel sections in the construction industry in the near future.Meanwhile, the geometric parameters of the parametric study were selected in accordance with commercially available channel sections.The section heights of unlipped channel sections varied between 150 mm and 250 mm, whilst the section thicknesses varied between 4 mm and 8 mm.The developed FE models had corner radii ranging from 12 mm to 18 mm, whilst bearing lengths varied from 50 mm to 150 mm.Three different yield stress values, from 700 MPa to 1000 MPa, were considered in the study.In total, 243 FE models of unlipped channel sections were generated during the parametric study.Table 5 summarises the key parameters that were considered during the parametric study.Since high strength steel sections have been applied in various areas [56,57], the material strength was chosen by considering those applications and by focusing on the potential applications of high strength channel sections in the construction industry in the near future.Meanwhile, the geometric parameters of the parametric study were selected in accordance with commercially available channel sections.The section heights of unlipped channel sections varied between 150 mm and 250 mm, whilst the section thicknesses varied between 4 mm and 8 mm.The developed FE models had corner radii ranging from 12 mm to 18 mm, whilst bearing lengths varied from 50 mm to 150 mm.Three different yield stress values, from 700 MPa to 1000 MPa, were considered in the study.In total, 243 FE models of unlipped channel sections were generated during the parametric study.The FE results, including web crippling resistances and failure modes of unlipped channel sections under the ETF loading case, are analysed in the coming sections to assess the web crippling behaviour and design provisions, with the aim of providing efficient design rules to increase the range of application of unlipped channel sections in structures.

Parametric Study Results
The ultimate web crippling capacity of lipped channel beams was obtained for 243 numerical studies.Figure 16 shows the typical web crippling failure modes of unlipped channel beams using corresponding pictures of the points indicated in the load vs. vertical deformation curve.The failure modes of the 200 mm section (t = 6 mm, N = 150 mm, r i = 12 mm and f y = 1000 MPa) are illustrated for the (i) initial, (ii) before-peak, (iii) at-peak and (iv) post-peak stages in Figure 16.The web crippling capacities of the 243 parametric studies are given in Table 6.The FE results, including web crippling resistances and failure modes of unlipped channel sections under the ETF loading case, are analysed in the coming sections to assess the web crippling behaviour and design provisions, with the aim of providing efficient design rules to increase the range of application of unlipped channel sections in structures.

Parametric Study Results
The ultimate web crippling capacity of lipped channel beams was obtained for 243 numerical studies.Figure 16 shows the typical web crippling failure modes of unlipped channel beams using corresponding pictures of the points indicated in the load vs. vertical deformation curve.The failure modes of the 200 mm section (t = 6 mm, N = 150 mm, ri = 12 mm and fy = 1000 MPa) are illustrated for the (i) initial, (ii) before-peak, (iii) at-peak and (iv) post-peak stages in Figure 16.The web crippling capacities of the 243 parametric studies are given in Table 6.The effect of various parameters, such as section thickness (t), internal radius (r i ) and bearing length (N) of the channel section, are investigated in the following subsections.

Effect of Thickness (t) on the Web Crippling Capacity
The web crippling behaviour of cold-formed steel sections varying in thickness is illustrated in Figure 17.It can be noticed that the web crippling strength increases with the thickness of the section.From the parametric results, the web crippling strengths of section 150 × 60 (r i = 12 mm, N = 50 mm and f y = 900 MPa) are 31.99kN, 73.20 kN and 126.7 kN when thicknesses are 4 mm, 6 mm and 8 mm, respectively.Therefore, almost a 129% and 296% increment in web crippling strength was observed for the 6 mm and 8 mm thick sections, respectively, compared to the 4 mm thick sections in the corresponding study.The effect of various parameters, such as section thickness (t), internal radius (ri) and bearing length (N) of the channel section, are investigated in the following subsections.

Effect of Thickness (t) on the Web Crippling Capacity
The web crippling behaviour of cold-formed steel sections varying in thickness is illustrated in Figure 17.It can be noticed that the web crippling strength increases with the thickness of the section.From the parametric results, the web crippling strengths of section 150 × 60 (ri = 12 mm, N = 50 mm and fy = 900 MPa) are 31.99kN, 73.20 kN and 126.7 kN when thicknesses are 4 mm, 6 mm and 8 mm, respectively.Therefore, almost a 129% and 296% increment in web crippling strength was observed for the 6 mm and 8 mm thick sections, respectively, compared to the 4 mm thick sections in the corresponding study.

Effect of Internal Radius (ri) on the Web Crippling Capacity
The effect of corner radius on the web crippling strength of cold-formed steel sections is shown in Figure 18. Figure 18 clearly depicts a reduction of the web crippling strength with corner radius.For example, the web crippling strength of section 150 × 60 (t = 4 mm, N = 50 mm and fy = 900 MPa) is decreased from 31.99 kN to 28.37 kN when the corner radius is increased from 12 mm to 18 mm.To the above-mentioned parametric study, web crippling strength is reduced by 5.8% and 11.3% for the sections with 15 mm and 18 mm corner radii, respectively, compared to the section with the 12 mm corner radius.

Effect of Bearing Length (N) on the Web Crippling Capacity
The web crippling capacity of cold-formed steel sections also depends on the bearing length of the loading and supporting plates.Figure 19 shows how the web crippling ca-

Effect of Internal Radius (r i ) on the Web Crippling Capacity
The effect of corner radius on the web crippling strength of cold-formed steel sections is shown in Figure 18. Figure 18 clearly depicts a reduction of the web crippling strength with corner radius.For example, the web crippling strength of section 150 × 60 (t = 4 mm, N = 50 mm and f y = 900 MPa) is decreased from 31.99 kN to 28.37 kN when the corner radius is increased from 12 mm to 18 mm.To the above-mentioned parametric study, web crippling strength is reduced by 5.8% and 11.3% for the sections with 15 mm and 18 mm corner radii, respectively, compared to the section with the 12 mm corner radius.

Effect of Bearing Length (N) on the Web Crippling Capacity
The web crippling capacity of cold-formed steel sections also depends on the bearing length of the loading and supporting plates.Figure 19 shows how the web crippling capacity of the unlipped channel section varies with the bearing length.An upward trend in web crippling capacity is noticed with an increase in the bearing length.The web crippling capacity of section 150 × 60 (t = 4 mm, r i = 12 m and f y = 900 MPa) with 50 mm bearing length is 31.99kN, which is increased to 36.31 kN and 40.23 kN when the bearing length is increased to 100 mm and 150 mm, respectively.This shows an almost 14% and 26% increment in the 100 mm and 150 mm bearing plates, respectively, compared to the 50 mm bearing length in the corresponding study.

AISI S100-16 (NAS 2016) and AS/NZS 4600:2018
The web crippling design equations available in the design standards AISI S100 [13] and AS/NZS 4600 [14] are identical.The AISI S100 [13] and AS/NZS 4600 [14] recommend a unified equation to predict web crippling capacity with different coefficients, which is illustrated by Equation (11).The web crippling capacity of the cold-formed steel section is developed based on various parameters, such as flat portion of web to thickness ratio

AISI S100-16 (NAS 2016) and AS/NZS 4600:2018
The web crippling design equations available in the design standards AISI S100 [13] and AS/NZS 4600 [14] are identical.The AISI S100 [13] and AS/NZS 4600 [14] recommend a unified equation to predict web crippling capacity with different coefficients, which is illustrated by Equation (11).The web crippling capacity of the cold-formed steel section is developed based on various parameters, such as flat portion of web to thickness ratio (d1/t), inside corner radius to thickness ratio (ri/t), bearing length to thickness ratio (N/t), The web crippling design equations available in the design standards AISI S100 [13] and AS/NZS 4600 [14] are identical.The AISI S100 [13] and AS/NZS 4600 [14] recommend a unified equation to predict web crippling capacity with different coefficients, which is illustrated by Equation (11).The web crippling capacity of the cold-formed steel section is developed based on various parameters, such as flat portion of web to thickness ratio (d 1 /t), inside corner radius to thickness ratio (r i /t), bearing length to thickness ratio (N/t), section thickness (t), material strength (f y ) and angle between the plane of the web and the bearing surface (θ).The coefficients vary depending on the (i) profiles, (ii) loading conditions and (iii) supporting conditions.In addition, the following variables are used: C-Coefficient, C R -Coefficient of inside bend radius, C N -Coefficient of bearing length and C h -coefficient of slenderness.The parametric study results were compared with the prediction of the equation available in AISI S100 [13] and AS/NZS 4600 [14] (Figure 20).Mean and COV values were calculated for the fraction of FE to prediction of unified equation as available in AISI S100 [13] and AS/NZS 4600 [14].The mean and COV were 0.62 and 0.11, respectively, and the comparison revealed that the codified equation is not capable of predicting the web crippling capacity of unlipped channel sections with high strength under the ETF loading condition with the unfastened supporting condition.The coefficients drawn from AISI S100 [13] are C = 2.00, C R = 0.11, C N = 0.37 and C h = 0.01 for this study.
Buildings 2022, 12, x FOR PEER REVIEW 20 of 29 section thickness (t), material strength (fy) and angle between the plane of the web and the bearing surface (θ).The coefficients vary depending on the (i) profiles, (ii) loading conditions and (iii) supporting conditions.In addition, the following variables are used: C-Coefficient, CR-Coefficient of inside bend radius, CN-Coefficient of bearing length and Ch-coefficient of slenderness.The parametric study results were compared with the prediction of the equation available in AISI S100 [13] and AS/NZS 4600 [14] (Figure 20).Mean and COV values were calculated for the fraction of FE to prediction of unified equation as available in AISI S100 [13] and AS/NZS 4600 [14].The mean and COV were 0.62 and 0.11, respectively, and the comparison revealed that the codified equation is not capable of predicting the web crippling capacity of unlipped channel sections with high strength under the ETF loading condition with the unfastened supporting condition.The coefficients drawn from AISI S100 [13] are C = 2.00, CR = 0.11, CN = 0.37 and Ch = 0.01 for this study.11)) [13,14].

Eurocode 3 Part 1-3 (CEN 2006)
Eurocode 3 Part 1-3 (CEN 2006) [12] suggests using separate equations to load cases of ETF, ITF, EOF and IOF separately.The equations do not consider the supporting conditions, such as the use of flanges fastened for support.Equations available in the Eurocode 3 Part 1-3 (CEN 2006) were also used to check the applicability to this study (Figure 21).For the comparison, mean and COV values were calculated.The mean and COV were 0.7 and 0.23, respectively.Therefore, the authors concluded that Eurocode 3 Part 1-3 (CEN 2006) cannot be used to predict the web crippling capacity of unlipped channel sections with high strength under the ETF loading condition.
Even though AISI S100 [13] and AS/NZS 4600 [14] recommend a unified equation to predict the web crippling capacity with different coefficients by considering various coldformed steel sections (all four loading conditions and supporting conditions), Eurocode 3   [12] suggests using separate equations to load cases of ETF, ITF, EOF and IOF separately.The equations do not consider the supporting conditions, such as the use of flanges fastened for support.Equations available in the Eurocode 3 Part 1-3 (CEN 2006) were also used to check the applicability to this study (Figure 21).For the comparison, mean and COV values were calculated.The mean and COV were 0.7 and 0.23, respectively.Therefore, the authors concluded that Eurocode 3 Part 1-3 (CEN 2006) cannot be used to predict the web crippling capacity of unlipped channel sections with high strength under the ETF loading condition.
Even though AISI S100 [13] and AS/NZS 4600 [14] recommend a unified equation to predict the web crippling capacity with different coefficients by considering various coldformed steel sections (all four loading conditions and supporting conditions), Eurocode 3 Part 1-3 (ECS, 2006) provides different design equations for each load case.Moreover, Eurocode 3 Part 1-3 (ECS, 2006) does not categorize the supporting conditions.The AISI S100 [13] and AS/NZS 4600 [14] equations are developed based on experimental results of Prabakaran et al. [58].Those experiments used different specimen lengths and test procedures and did not consider the higher yield strength of the material.Meanwhile, the design rules available in the Eurocode 3 Part 1-3 (ECS, 2006) are developed based on the old AISI standard [59].The design equations are not suitable for commercially available high strength sections due to the inconsistencies of past research and the parametric limitations in design standards., 2006) does not categorize the supporting conditions.The AISI S100 [13] and AS/NZS 4600 [14] equations are developed based on experimental results of Prabakaran et al. [58].Those experiments used different specimen lengths and test procedures and did not consider the higher yield strength of the material.Meanwhile, the design rules available in the Eurocode 3 Part 1-3 (ECS, 2006) are developed based on the old AISI standard [59].The design equations are not suitable for commercially available high strength sections due to the inconsistencies of past research and the parametric limitations in design standards.Gunalan and Mahendran [22] studied the web crippling capacity of unlipped channel sections under all four load cases.The study covered a section depth of 75 mm to 300 mm, thickness of 3.8 mm to 6.0 mm and corner radius of 4.0 mm to 8.0 mm and proposed modified coefficients to the codified equation available in AISI S100 [13] and AS/NZS 4600 [14].The proposed coefficients were also compared in this study (Figure 22), and mean and COV of the comparison were 0.78 and 0.08, respectively.This study limited the investigation to 450 MPa, and due to the unsatisfied numerical facts, the study was also not suitable for calculating the web crippling capacity of unlipped channels with high strength material under the ETF load case.

Sundarajah study
Sundararajah [21] also studied the web crippling capacity of unlipped channel sections under all four load cases.The study conducted by Sundararajah [20] accommodated a section depth of 100 mm to 250 mm, thicknesses of 1.0 mm to 2.4 mm and yield strengths of 506 and 581 MPa.The researchers proposed coefficients to the standard equation available in AISI S100 [13] and AS/NZS 4600 [14].The proposed coefficients were also compared in this study, and mean and COV of the comparison for Sundararajah [21] were 0.74 and 0.08, respectively.Due to the parametric limitations of their research and the unsatis-  Gunalan and Mahendran [22] studied the web crippling capacity of unlipped channel sections under all four load cases.The study covered a section depth of 75 mm to 300 mm, thickness of 3.8 mm to 6.0 mm and corner radius of 4.0 mm to 8.0 mm and proposed modified coefficients to the codified equation available in AISI S100 [13] and AS/NZS 4600 [14].The proposed coefficients were also compared in this study (Figure 22), and mean and COV of the comparison were 0.78 and 0.08, respectively.This study limited the investigation to 450 MPa, and due to the unsatisfied numerical facts, the study was also not suitable for calculating the web crippling capacity of unlipped channels with high strength material under the ETF load case.

Sundarajah Study
Sundararajah [21] also studied the web crippling capacity of unlipped channel sections under all four load cases.The study conducted by Sundararajah [20] accommodated a section depth of 100 mm to 250 mm, thicknesses of 1.0 mm to 2.4 mm and yield strengths of 506 and 581 MPa.The researchers proposed coefficients to the standard equation available in AISI S100 [13] and AS/NZS 4600 [14].The proposed coefficients were also compared in this study, and mean and COV of the comparison for Sundararajah [21] were 0.74 and 0.08, respectively.Due to the parametric limitations of their research and the unsatisfactory comparison results to the parametric study, the authors concluded that the Sundarajah study [21] was also not suitable for calculating the web crippling capacity of unlipped channels with high strength material under the ETF load case.The comparison is illustrated in Figure 23.Sundarajah study [21] was also not suitable for calculating the web crippling capacity of unlipped channels with high strength material under the ETF load case.The comparison is illustrated in Figure 23.Based on a comparison with the current design equations, the equations proposed in the literature review were inadequate to predict the web crippling strength of unlipped Sundarajah study [21] was also not suitable for calculating the web crippling capacity of unlipped channels with high strength material under the ETF load case.The comparison is illustrated in Figure 23.

Proposed Equation-1
Based on a comparison with the current design equations, the equations proposed in the literature review were inadequate to predict the web crippling strength of unlipped

Proposed Equations 7.1. Proposed Equation-1
Based on a comparison with the current design equations, the equations proposed in the literature review were inadequate to predict the web crippling strength of unlipped channel sections with high strength under ETF loading conditions with unfastened supporting conditions.Therefore, the authors proposed coefficients to the design equation available in the design standards [13,14].The proposed Equation-1 is illustrated by Equation (12).The mean, COV and capacity reduction factors (φw) were calculated for the comparison of the FE values to the proposal and were 1.00, 0.07 and 0.90, respectively.Therefore, the proposed equation (Equation ( 12)) can be used with a reduction factor (φw) of 0.9.The comparison of FE values to the prediction of the proposed equation is shown in Figure 24.
Buildings 2022, 12, x FOR PEER REVIEW 23 of 29 channel sections with high strength under ETF loading conditions with unfastened supporting conditions.Therefore, the authors proposed coefficients to the design equation available in the design standards [13,14].The proposed Equation-1 is illustrated by Equation (12).The mean, COV and capacity reduction factors (ϕw) were calculated for the comparison of the FE values to the proposal and were 1.00, 0.07 and 0.90, respectively.Therefore, the proposed equation (Equation ( 12)) can be used with a reduction factor (ϕw) of 0.9.The comparison of FE values to the prediction of the proposed equation is shown in Figure 24.12)).

Proposed Equation-2
Previous studies on web crippling strength accommodate a strength factor (Cf) addition to the equation available in the design standards [19].Therefore, strength factor was incorporated in this study to evaluate the web crippling strength of unlipped channel beams with high strength material.In addition, the incorporation of strength factor predicted the web crippling strength more accurately than previous equations (Proposed Equation-1 (Equation ( 12))).Proposed Equation-2 is illustrated by Equation ( 13).Mean, COV and capacity reduction factors (ϕw) were calculated to compare the parametric study to the prediction of proposed Equation-2.The mean and COV were 1.00 and 0.05, with a capacity reduction factor (ϕw) of 0.91.The comparison is illustrated in Figure 25.Therefore, both the proposed equations are capable of determining the web crippling capacity of unlipped channel beams with given reduction factors (ϕw) for high strength material under ETF loading conditions and unfastened support.12)).

Proposed Equation-2
Previous studies on web crippling strength accommodate a strength factor (Cf) addition to the equation available in the design standards [19].Therefore, strength factor was incorporated in this study to evaluate the web crippling strength of unlipped channel beams with high strength material.In addition, the incorporation of strength factor predicted the web crippling strength more accurately than previous equations (Proposed Equation-1 (Equation ( 12))).Proposed Equation-2 is illustrated by Equation ( 13).Mean, COV and capacity reduction factors (φw) were calculated to compare the parametric study to the prediction of proposed Equation-2.The mean and COV were 1.00 and 0.05, with a capacity reduction factor (φw) of 0.91.The comparison is illustrated in Figure 25.Therefore, both the proposed equations are capable of determining the web crippling capacity of unlipped channel beams with given reduction factors (φw) for high strength material under ETF loading conditions and unfastened support.13)).

Direct Strength Method (DSM)
The design equations available in the current design rules were proposed and modified (Proposed Equation-1 and Proposed Equation-2) for this study to predict the web crippling capacity of unlipped channel sections with high strength material under the ETF load case with unfastened flanges to supporting plates.Few researchers [21,22,39] have proposed web crippling-based DSM equations.Keerthan and Mahendran [39] conducted crippling studies for LSB, whereas Gunalan and Mahendran [22] and Sundarajah [21] studied the web crippling behaviour of unlipped channel sections.The researchers obtained the critical buckling load (Rb,cr) and yield load (Rb,y) to develop their DSM equation.The yield load was derived based on the 45 ° load distribution, which was determined from the end of the bearing plate to middle of the web.A similar method was incorporated in this study to propose a DSM-based equation to predict the web crippling behaviour of unlipped channel sections under ETF loading conditions with high strength material.The DSM equation (Equations ( 14) and ( 15)) was proposed by Sundarajah [21], which was compared with this study to ensure the applicability of the equation.The mean and COV were observed as 1.28 and 0.13, respectively.The comparison revealed that Equations ( 14) and (15) were not suitable for predicting the web crippling strength of unlipped channel sections with high strength material under ETF loading conditions.The comparison is shown in Figure 26.

Direct Strength Method (DSM)
The design equations available in the current design rules were proposed and modified (Proposed Equation-1 and Proposed Equation-2) for this study to predict the web crippling capacity of unlipped channel sections with high strength material under the ETF load case with unfastened flanges to supporting plates.Few researchers [21,22,39] have proposed web crippling-based DSM equations.Keerthan and Mahendran [39] conducted crippling studies for LSB, whereas Gunalan and Mahendran [22] and Sundarajah [21] studied the web crippling behaviour of unlipped channel sections.The researchers obtained the critical buckling load (R b,cr ) and yield load (R b,y ) to develop their DSM equation.The yield load was derived based on the 45 • load distribution, which was determined from the end of the bearing plate to middle of the web.A similar method was incorporated in this study to propose a DSM-based equation to predict the web crippling behaviour of unlipped channel sections under ETF loading conditions with high strength material.The DSM equation (Equations ( 14) and ( 15)) was proposed by Sundarajah [21], which was compared with this study to ensure the applicability of the equation.The mean and COV were observed as 1.28 and 0.13, respectively.The comparison revealed that Equations ( 14) and (15) were not suitable for predicting the web crippling strength of unlipped channel sections with high strength material under ETF loading conditions.The comparison is shown in Figure 26.
where, λ = R b,y R b,cr .A new DSM equation was proposed following the comparison of available equations.In this study, the critical buckling load (Rb,cr) was calculated based on Equation ( 16).The buckling coefficient (k) and yield load (Rb,y) were determined from Equations ( 17) and ( 18), respectively.The coefficients in the derivation of k value were determined from the Sundarajah study on unlipped channel sections [21].Furthermore, the supporting equations to calculate the yield load (Rb,y) are given by Equations ( 19)-( 21).Based on the calculated critical buckling load (Rb,cr) and yield load (Rb,y), a DSM equation was developed and is given by Equations ( 22) and (23).Figure 27 compares the proposed DSM and FE results of unlipped channel sections with high strength under ETF loading conditions.The mean, COV and reduction factor (ϕw) of the comparison were 1.00, 0.06 and 0.91, respectively.Therefore, the proposed DSM can be used with the reduction factor (ϕw) to predict the web crippling capacities of unlipped channel sections with high strength material under the ETF loading case.A new DSM equation was proposed following the comparison of available equations.In this study, the critical buckling load (R b,cr ) was calculated based on Equation ( 16).The buckling coefficient (k) and yield load (R b,y ) were determined from Equations ( 17) and ( 18), respectively.The coefficients in the derivation of k value were determined from the Sundarajah study on unlipped channel sections [21].Furthermore, the supporting equations to calculate the yield load (R b,y ) are given by Equations ( 19)-( 21).
Based on the calculated critical buckling load (R b,cr ) and yield load (R b,y ), a DSM equation was developed and is given by Equations ( 22) and (23).Figure 27 compares the proposed DSM and FE results of unlipped channel sections with high strength under ETF loading conditions.The mean, COV and reduction factor (φw) of the comparison were 1.00, 0.06 and 0.91, respectively.Therefore, the proposed DSM can be used with the reduction factor (φw) to predict the web crippling capacities of unlipped channel sections with high strength material under the ETF loading case.

Conclusions
This paper presents a detailed investigation of the web crippling behaviour of high strength unlipped channel sections under the ETF load case.Earlier, a developed numerical model was validated with experiment results.Since the validation results showed great agreement with the experimental results, a parametric plan was developed to analyse the web crippling behaviour of unlipped channel sections with high strength under the ETF load case using ABAQUS.In total, 243 numerical models were developed to obtain the ultimate web crippling capacity of the sections by considering various parameters.The web crippling results were compared with the existing design code equations as well as equations proposed by various researchers in this respective area.Since the comparison discloses that available design equations are too conservative to predict the ultimate web crippling capacity under the ETF load case, a web crippling design equation according to the current design standards was proposed.Moreover, new design equations were proposed based on the DSM.This paper concludes with a modified equation with respect to current design standards and new DSM design equations, which accurately predicts the ultimate web crippling capacity of unlipped channel sections with high strength under the ETF load case.

Figure 3 .
Figure 3. Finite Element model for full and half models.

Figure 3 .
Figure 3. Finite Element model for full and half models.

Figure 4 .
Figure 4. Applied element types and mesh sizes.

Figure 4 .
Figure 4. Applied element types and mesh sizes.

Figure 6 .
Figure 6.Boundary conditions of full and half models.

Figure 6 .
Figure 6.Boundary conditions of full and half models.

Figure 6 .
Figure 6.Boundary conditions of full and half models.

Figure 11 .
Figure 11.Failure mode comparison for full and half models of ETF H120X120X4 with a bearing plate length of 120 mm [2].

Figure 12 .
Figure 12.Failure mode comparison for full and half models of ETF-UC10012 with 50 mm bearing plate [21].

Figure 11 . 29 Figure 11 .
Figure 11.Failure mode comparison for full and half models of ETF H120X120X4 with a bearing plate length of 120 mm [2].

Figure 12 .
Figure 12.Failure mode comparison for full and half models of ETF-UC10012 with 50 mm bearing plate [21].Figure 12. Failure mode comparison for full and half models of ETF-UC10012 with 50 mm bearing plate [21].

Figure 12 . 29 Figure 13 .
Figure 12.Failure mode comparison for full and half models of ETF-UC10012 with 50 mm bearing plate [21].Figure 12. Failure mode comparison for full and half models of ETF-UC10012 with 50 mm bearing plate [21].

Figure 14 .
Figure 14.Failure mode comparison for full and half models of ETF-UC20015 with 150 mm bearing plate [21].

Figure 14 .
Figure 14.Failure mode comparison for full and half models of ETF-UC20015 with 150 mm bearing plate [21].

Figure 14 .
Figure 14.Failure mode comparison for full and half models of ETF-UC20015 with 150 mm bearing plate [21].

Figure 15 .
Figure 15.Notations used for key cross-sectional dimensions.

Figure 15 .
Figure 15.Notations used for key cross-sectional dimensions.

s 2022 ,
12, x FOR PEER REVIEW 15 of 29

Figure 16 .
Figure 16.Typical web crippling failure modes of unlipped channel beam under progressive loading.(a) Initial stage.(b) Before peak stage.(c) At peak stage.(d) Post peak stage.

Figure 17 .
Figure 17.Effect of thickness (t) on the web crippling strength.

Figure 17 .
Figure 17.Effect of thickness (t) on the web crippling strength.

Buildings 2022 ,
12, x FOR PEER REVIEW 19 of 29in web crippling capacity is noticed with an increase in the bearing length.The web crippling capacity of section 150 × 60 (t = 4 mm, ri = 12 m and fy = 900 MPa) with 50 mm bearing length is 31.99kN, which is increased to 36.31 kN and 40.23 kN when the bearing length is increased to 100 mm and 150 mm, respectively.This shows an almost 14% and 26% increment in the 100 mm and 150 mm bearing plates, respectively, compared to the 50 mm bearing length in the corresponding study.

Figure 18 .
Figure 18.Effect of internal radius (ri) on the web crippling strength.

Figure 19 .
Figure 19.Effect of bearing length (N) on the web crippling strength.

Figure 18 .
Figure 18.Effect of internal radius (r i ) on the web crippling strength.

Figure 18 .
Figure 18.Effect of internal radius (ri) on the web crippling strength.

Figure 19 .
Figure 19.Effect of bearing length (N) on the web crippling strength.

Figure 22 .
Figure 22.Comparison of parametric results and equation developed by [22].

Figure 22 .
Figure 22.Comparison of parametric results and equation developed by [22].

Figure 22 .
Figure 22.Comparison of parametric results and equation developed by [22].

Figure 23 .
Figure 23.Comparison of parametric results and equation developed by [21].

Figure 23 .
Figure 23.Comparison of parametric results and equation developed by [21].
Cb-General coefficient = 0.59, Cb,r-Coefficient of inside bent radius to thickness ratio = 0.01, Cb,w-Coefficient of web slenderness ratio = 0.05, Cb,l-Coefficient of bearing length to thickness ratio = 0.4 and Cb,b-Coefficient of flange width to thickness ratio = 0.01.

Figure 27 .
Figure 27.Comparison of parametric results and proposed DSM equation.
Coefficient of inside bent radius to thickness ratio C b,w Coefficient of web slenderness ratio C b,l Coefficient of bearing length to thickness ratio C b,b Coefficient of flange width to thickness ratio N m Yield mechanism length r m Inside bent radius measured along the middle of the section r ext External bent radius b f Flange width θ Angle between the plane of the web and the bearing surface FE Finite Element COV Coefficient of variation
Note: H and V denotes the experimental study with a material strength of 700 MPa and 900 MPa respectively.

Table 3 .
[21]arison of web crippling test strengths with FE results for unlipped channel sections under the ETF loading condition[21].

Table 4 .
Summary of the statistical analysis of validation.

Table 5 .
Summary of key parameters considered.

Table 5 .
Summary of key parameters considered.
Strain at the ultimate strength fu (ultimate strain)ε c,avThrough thickness averaged plastic strain induced during the forming of the corner regions of press-braked and cold-rolled sections f y