Microsoft Word - 01-3231_s_ETASR_V10_N1_pp5099-5103 Engineering, Technology & Applied Science Research Vol. 10, No. 1, 2020, 5099-5103 5099 www.etasr.com Tran & Nguyen: Reliability Assessment of Steel Plane Frame's Buckling Strength Considering … Reliability Assessment of Steel Plane Frame's Buckling Strength Considering Semi-rigid Connections Ngoc Long Tran Department of Civil Engineering Vinh University, Vinh, Vietnam longtn@vinhuni.edu.vn Trong Ha Nguyen Department of Civil Engineering Vinh University, Vinh, Vietnam trongha@vinhuni.edu.vn Abstract—The buckling strength of a structure depends hugely on material properties, geometry dimensions, and boundary conditions that are potentially random. This paper presents an assessment of the safety of buckling strength of steel plane frame semi-rigid beam-column connections. The Newton-Raphson method is used to solve the nonlinear equation for the buckling limit of the column. The reliability of the structure is evaluated by the Monte Carlo simulation method. The effects of input parameters are also investigated. Keywords-plane frames steel; semi-rigid connection; reliability; Monte Carlo simulation I. INTRODUCTION The analysis of steel frame with semi-rigid connections is always an interesting research topic. Straight-line common beam method was used to classify semi-rigid connection in [1]. Authors in [2] considered the stiffness of the connection in the research on the moment distribution method. Author in [3] investigated the influence of semi-rigid connections on the frame subjecting to shake. Research and application of the stiffness matrix method to analyze the plane frame with semi- rigid connections, in which the stiffness matrix and nodal force vector of each element depend on the linear stiffness of the connection, was conducted in [4]. A novel method to analyze the spatial frame considering the material nonlinearity and large deflection effect was presented in [5]. Linear and nonlinear semi-rigid connection models have been presented in [6]. There are many studies about buckling analysis of steel structures considering semi-rigid connection, like the in-plane buckling behavior of pitched-roof steel frames with semi-rigid connections and elastic buckling [7], second-order behavior of pitched-roof steel frames [8], and the buckling analysis of non- prismatic columns using slope-deflection method [9]. The buckling and post-buckling behavior of gabled frames was studied in [10]. Stability analysis of semi-rigid composite frames has been studied in [11]. Stiffness analysis on semi- rigid joints in gabled frames was studied in [12], and the elastic buckling of end-loaded, tapered, cantilevered beams with initial curvature in [13]. In this research, buckling analysis of steel gabled frames with tapered members and flexible connections was used. Concerning the reliability assessment of the structure, authors in [14] assessed the reliability of steel frame under dynamic loads generated from the El Centro earthquake in 1940 using the β probability index method. The same method was used in [15] to study the steel frame according to the plastic limit state. The reliability of the steel structure under corrosion was studied in [16]. The Monte Carlo simulation method to investigate steel frame with flexible beam-column joint was utilized in [6]. Rotational stiffness and moment- resisting of the connection were referred in [17]. This model of flexible joint is not correct [18]. Instead a novel model based on empirical tests and Eurocode-3 was proposed. This paper performs reliability assessment of the buckling strength for steel plane frame considering semi-rigid connections. An algorithm using Monte Carlo simulation is proposed to be used in analysis and assessment. In addition, the effects of input parameters, which are safety factor and coefficient of variation, on the reliability of structures have been examined. II. THEORETICAL FRAMEWORK A. Stiffness of the Beam-column Joint The rotational stiffness of flexible beam-column joint of an I-section steel portal frame is determined by the experimental formula proposed in [18]: ini ϕ ϕ η = (1) where ϕ is the elastic rotational stiffness of the joint (kNm/rad), η is the hardness adjustment coefficient depending on the connection type and presented in Eurocode-3, ini ϕ is the initial stiffness (kNm/rad) determined as [18]: 0.44 1.2 0.35 0.005 1 2ini c b p K h h t d Kϕ = − (2) where c h (mm) is the height of the column section (HEB), b h (mm) is the height of the beam section (IPE), p t (mm) is the thickness of the end plate, and d (mm) is the bolt diameter, Corresponding author: Trong Ha Nguyen Engineering, Technology & Applied Science Research Vol. 10, No. 1, 2020, 5099-5103 5100 www.etasr.com Tran & Nguyen: Reliability Assessment of Steel Plane Frame's Buckling Strength Considering … and K1=1.5 and K2=19211 are empirical coefficients from the experimental results in [18]. Because of their empirical origin, these coefficients possess potential randomness. B. Monte Carlo Simulation Method Monte Carlo simulation method is based on the use of pseudo-random numbers and the law of large number to assess the reliability of any system. If the safe domain is defined by the condition ( ) 0,f >X where X is a random vector containing all the input random variables, the unsafe probability of the system is determined by: ( ) ( ) ( )XX 0 X 0f f fP I f x dx E I< < = =  ∫ (3) where ( )X 0fI < is the indicator function defined by: ( ) ( ) ( )X 0 1 if X 0 0 if X 0 f f I f < < =  ≥ (4) According to the theory of statistics, if we have N realizations of the random vector ,X by propagating the randomness, we obtain a sample of N realizations of the indicator function. The expected value of the indicator function can be approximatively determined by taking the mean of the samples: ( ) ( )X 0 X 0 1 1 N i f f f i P E I I N < < =  = =  ∑ � (5) A 95% confidence interval of the estimation is defined by [20]: 1 1 1 200 1 200 f f f f f f f P P P P P NP NP    − −    − ≤ ≤ +         � � � � � � (6) C. Displacement Method This method is commonly used when analyzing stability structures by analytical method. This method is determined in [19] and the stability equations are shown by: 11 1 1 ... det 0 ... j i ij r r r r     =      � � � (7) where ij r is the reaction from a generalized relocation of ties. III. METHODOLOGY A. Reliability of Buckling Strength of the Steel Plane Frame Considering Semi-rigid Connections 1) Safe Condition Although the buckling strength of the steel plane frame considering semi-rigid connections beam-column is determined by (13), in reality we often introduce a safety factor n and the safe condition of the steel plane frame considering semi-rigid connections beam-column will be written as: cr P P n ≤ (8) where P is the external load which can fluctuate randomly. 2) Deterministic Model and Uncertainty Model The deterministic model deals with the above buckling analysis problem, in which the input parameters are those of geometry ( )1 2 1 2, , , , , , , , , , , , , ,c fc fc wc b fb fb wb pL L H H h b t t h b t t t d material (E) load (P) and flexible joint (K1, K2). This model can be written as: ( ) cr P = ℑ X (9) with 1 2 1 2 1 2, , , , , , , , , , , , , , , , ,c fc fc wc b fb fb wb pL L H H h b t t h b t t t d E P K K =  X . In this paper, we consider the steel plane frame considering semi- rigid connections with the deterministic input parameters shown in Table I. The uncertainty model is constructed based on the deterministic model by taking into account the randomness of some input paramters. In this paper, we distinct two vectors of input parameters: the first one of the parameters assumed to be deterministic 1 1 2 1 2, , , , , , , , , , , , ,c fc fc wc b fb fb wb pL L H H h b t t h b t t t d =  X and the second one of the parameters assumed to be ramdom [ ]2 1 2( ) ( ), ( ), ( ), ( )P K K Eω ω ω ω ω=X with ω representing the randomness of the parameters. This model can be written as: ( )1 2( ) , ( )crP ω ω= ℑ X X (10) TABLE I. DETERMINISTIC INPUT PARAMETERS Variable Value Variable Value H1 (cm) 400 L1 (cm) 400 H2 (cm) 800 L2 (cm) 600 hc (mm) 400 hb (mm) 400 bfc (mm) 200 bfb (mm) 200 tfc (mm) 15 tfb (mm) 15 twc (mm) 10 twb (mm) 10 tp (mm) 20 d (mm) 22 TABLE II. RANDOM INPUT VARIABLES AND THEIR REPRESENTATIVE PARAMETERS Random variable ( )P kN Law of probability Norm Representative parameters P µ P CV 1.5088×10 4 0.1 Random variable 2 (kN/cm )E Law of probability Norm Representative parameters E µ ECV 2×10 4 0.1 K1 K2 Uniform Uniform Reference Interval Reference Interval 1.5 [0.95 1.05]÷ 19211 [0.95 1.05]÷ Engineering, Technology & Applied Science Research Vol. 10, No. 1, 2020, 5099-5103 5101 www.etasr.com Tran & Nguyen: Reliability Assessment of Steel Plane Frame's Buckling Strength Considering … B. Column Reliability Assessment by Monte Carlo Simulation By introducing the uncertainty model in the Monte Carlo simulation method, we obtain the scheme of the reliability assessment of the steel plane frame considering semi-rigid connections shown in Figure 1. Fig. 1. Scheme of the reliability assessment of the steel frame considering semi-rigid connection using Monte Carlo simulation method IV. NUMERICAL RESULTS A. Validation of the Steel Frame Buckling Strength Considering Semi-rigid Connections In order to validate the Matlab code, the steel frame considering semi-rigid connection shown in Figure 2 was considered with the input parameters shown in Table III. The validation of the program was performed in articulation joints (stiffness of the joint is zero) and rigid joints (stiffness of the joint is infinitely high). The result of the computer program and its comparison with the result of [19] is shown in Table IV. TABLE III. INPUT PARAMETERS Geometry and section details Beam L (cm) hwb (cm) bfb (cm) tfb (cm) twb (cm) 600 40.0 20.0 2.0 1.0 Column H (cm) hwc (cm) bfc (cm) tfc (cm) twc (cm) 600 40.0 20.0 2.0 1.0 Fig. 2. Plane frame semi-rigid connections for validation TABLE IV. VALIDATION BY COMPARISON WITH [19] Articulation joints Rigid joints Matlab code [2] Error (%) Matlab code [2] Error (%) Pcr (kN) 545.47 545.16 0.31 759.86 761.00 1.14 We can observe that the error is very small in the case of button joint (<0.31%) and slightly increases in the other case (<1.14%). This result can be explained by the fact that in our program the rigid joint is modeled by setting a very high value to the stiffness of the joints. This result confirms the reliability of our program. B. Convergence of the Monte Carlo Simulation Consider the steel plane frame semi-rigid beam-column connections shown in Figure 3 with the deterministic input parameters and the random input parameters presented in Table I and Table II respectively. The compression load P (kN) and Young’ modulus are assumed to be normal variables with mean , P E µ µ and coefficients of variation , , P E CV CV whereas the empirical coefficients K1, K2 are assumed to be uniform variables in the interval [0.95, 1.05] of the reference value [18]. Figure 4 shows the convergence of the safe probability of the frame in the Monte Carlo simulation to the value of 0.9312 or 93.12% after about 1051 samples in 30 minutes. The used convergence criterion of 2.5% justifies the confidence of the estimated reliability. This result also shows that although the safety factor was taken as 1.1 in the analysis, because of the randomness of some input parameters the reliability of the structure is only 93.12%. Thus the assessment of the reliability of the structure is necessary. Fig. 3. Semi-rigid connection plane frame analysis Fig. 4. Convergence of the safe probability in the Monte Carlo simulation C. Coefficient of Variation and Safety Factor Effects We know that the variation of input random variables and the safety factor influence directly but inversely the safe probability of the structure. Thus in order to clear the effect of L H P I I c d d d d d 1 2 2 2 1 1 2 P 4P 4P P L L L L L H H I I 2I I I I I 2I 2I c c c Engineering, Technology & Applied Science Research Vol. 10, No. 1, 2020, 5099-5103 5102 www.etasr.com Tran & Nguyen: Reliability Assessment of Steel Plane Frame's Buckling Strength Considering … these parameters, one reconsiders the above column with different coefficients of variation of the compression load CV=0.05; 0.1; 0.15; 0.2; 0.25 and different safety factors n=1.1; 1.15; 1.2; 1.25; 1.3. The randomness of the empirical coefficients K1, K2 is assumed to be unchanged in the interval [0.95, 1.05] of the reference value. The results of the safe probability are listed in Table V and presented in Figure 5. TABLE V. EFFECTS OF THE COEFFICIENT OF VARIATION AND OF THE SAFETY FACTOR ON THE SAFE PROBABILITY OF THE FRAMES n CV 0.05 0.10 0.15 0.20 0.25 1.10 0.0276 0.1514 0.2566 0.2833 0.3500 1.15 0.0000 0.0676 0.1302 0.2222 0.2553 1.20 0.0000 0.0225 0.1166 0.1313 0.2166 1.25 0.0000 0.0067 0.0572 0.0964 0.1722 1.30 0.0000 0.0000 0.0278 0.0609 0.1031 (a) (b) Fig. 5. (a) Effect of the coefficient of variation and (b) effect of the safety factor on the safe probability of the column We can easily observe the inverse effect of the coefficient of variation and of the safety factor in Figure 5. The safe probability of the column decreases when the coefficient of variation increases, whereas it increases when the safety factors increase. This result seems to be obvious, but it has a very important significance. It shows that if there are many input random parameters or furthermore, with high randomness in the structural design or in the optimization problem, the use of local coefficients such as the overload coefficient seems to be not sufficient. The structure can be in a dangerous state. In this case, it is necessary to determine a global safety factor, as is done in this study, to assure the absolute safety of the structure. For example, in this test, if the coefficient of variation is 0.05, the global safety factor needs to be only 1.15 to obtain a safe probability of 100%. If the coefficient of variation is 0.1, the global safety factor needs to be 1.3. V. CONCLUSIONS This paper studies the reliability assessment of a steel plane frame with semi-rigid beam-column joint connections. The authors have successfully established the deterministic model using the Newton-Raphson method to solve the transcendental equation of the column’s equilibrium. This model was then combined with the Monte Carlo simulation method to construct the column’s reliability assessment program. Parametric tests were also performed to study the input parameter effects on the reliability of the structure. The obtained results show the significance of this research. ACKNOWLEDGMENT This research was funded by Vinh University. REFERENCE [1] C. Batho, H. C. Rowan, “Investigations on beam and stanchion connections”, in: Second report of the steel structures research committee, pp. 61-137, Department of Scientific and Inductrial Research, His Majesty’s Stationery Office, London, 1934 [2] J. C. Rathburn, “Elastic properties of riveted connections”, Transactions of the American Society of civil Engineers, Vol. 101, pp. 523-596, 1936 [3] B. Sourochnikoff, Wind stresses in semi-rigid connections of steel framework, Transactions of the American Society of Civil Engineers, Vol. 115, No. 1, pp. 382–402. 1950 [4] G. R. Monforton, T. S. Wu, “Matrix analysis of semi-rigidly connected frames”, Journal of Structural Division, Vol. 89, No. 6, pp. 13-42, 1963 [5] S. E. Kim, S. H. Choi, “Practical advanced analysis for semi-rigid space frames”, International Journal of Solids and Structures, Vol. 38, pp. 9111-9131, 2001 [6] M. A. Hadianfard, R. Razani, “Effects of semi-rigid behavior of connections in the reliability of steel frames”, Structural Safety, Vol. 25, No. 2, pp. 123-138, 2003 [7] N. Silvestre, A. Mesquita, D. Camotim, L. Silva, “In-plane buckling behavior of pitched-roof steel frames with semi-rigid connections”, in: Frames with partially restrained connections, SSRC 1998, Theme Conference Workshop Volume, pp. 21-34, 2000 [8] N. Silvestre, D. Camotim, “Elastic buckling and second-order behavior of pitched-roof steel frames”, Journal of Constructional Steel Research, Vol. 63, No. 6, pp. 804-818, 2007 [9] H. Tajmir Riahi, A. Shojaei Barjoui, S. Bazazzadeh, S. M. A. Etezady, “Buckling analysis of non-prismatic columns using slope-deflection method”, 15th World Conference on Earthquake Engineering, Lisbon, Portugal, September 24-28, 2012 [10] A. S. Vlahinos, A. Cervantes, “Buckling and postbuckling behaviour of gabled frames”, Mathematical and Computer Modelling, Vol. 14, pp. 873-876, 1990 [11] J. F. Wang, G. Q. Li, “Stability analysis of semi-rigid composite frames”, Steel and Composite Structures, Vol. 7, No. 2, pp. 119-133, 2007 [12] Y. Wang, Y. J. Liu, Y. F. Xu, “Stiffness analysis on semi-rigid joints in gabled frames”, International Conferemce on Civil Engineering, Architecture and Building Materials, Haikou, China, June 18-20, 2011 [13] J. F. Wilson, D. J Strong, “Elastic buckling of end-loaded, tapered, cantilevered beams with initial curvature”, Structural Engineering and Mechanics, Vol. 5, No. 3, pp. 257-268, 1997 Engineering, Technology & Applied Science Research Vol. 10, No. 1, 2020, 5099-5103 5103 www.etasr.com Tran & Nguyen: Reliability Assessment of Steel Plane Frame's Buckling Strength Considering … [14] A. D. Kiureghian,“Structural reliability methods for seismic safety assessment: A review”, Engineering Structures, Vol. 18, No. 6, pp. 412- 424, 1990 [15] H. P. Hon, S. Wang, Reliability of steel frame systems with semi-rigid connections, ICLR Research Paper Series No. 35, University of Western Ontario, 2003 [16] R. E. Melchers, “The effect of corrosion on the structural reliability of steel offshore structures”, Corrosion Science, Vol. 47, No. 10, pp. 2391– 2410, 2005 [17] N. Kishi, W. F Chen, Y. Goto, K. G Matsuoka, “Design aid of semi- rigid connections for frame analysis”, Engineering Journal, Vol. 30, No. 3, pp. 90-107, 1993 [18] A. Kozlowski, R. Kowalczyk,M. Gizejowsk, “Estimation of the initial stiffness and moment resistance of steel and composite joints”, CTBUH 8th World Congress, Dubai, March 3-5, 2008 [19] V. Q Anh, Advance analysis and design frame semi-rigid connections, Publishing of Contruction, 2013 (in Vietnamese) [20] M. Lemaire, Fiabilite des Structures: Couplage mescano-fiabiliste, European Journal of Computational Mechanics, Vol. 15, No. 7-8, pp. 989-992, 2005