Nonlinear stability analysis of the frame structures

Abstract

In this paper the phenomenon of instability of frames in elasto-plastic domain was investigated. Numerical analysis was performed by the finite element method. Stiffness matrices were derived using the trigonometric shape functions related to exact solution of the differential equation of bending according to the second order theory. When the buckling of structure occurs in plastic domain, it is necessary to replace the constant modulus of elasticity E with the tangent modulus Et. Tangent modulus is stress dependent function and takes into account the changes of the member stiffness in the inelastic range. For the purposes of numerical investigation in this analysis, part of the computer program ALIN was created in a way that this program now can be used for elastic and elasto-plastic stability analysis of frame structures. This program is developed in the C++ programming language. Using this program, it is possible to calculate the critical load of frames in the elastic and inelastic domain. In this analysis, the algorithm for the calculation of buckling lengths of compressed columns of the frames was also established. The algorithm is based on the calculation of the global stability analysis of frame structures. Results obtained using this algorithm were compared with the approximate solutions from the European (EC3) and national (JUS) standards for the steel structures. By the given procedure in this paper it is possible to follow the behavior of the plane frames in plastic domain and to calculate the real critical load in that domain.

U ovom radu istraživan je fenomen gubitka stabilnosti okvirnih nosača u elasto-plastičnoj oblasti. Numerička analiza je sprovedena primenom metode konačnih elemenata. Matrice krutosti su izvedene korišćenjem trigonometrijskih interpolacionih funkcija koje se odnose na tačno rešenje diferencijalne jednačine savijanja štapa prema teoriji drugog reda. U slučaju kada se izvijanje konstrukcije dešava u plastičnoj oblasti, konstantan modul elastičnosti E u matrici krutosti zamenjen je tangentnim modulom Et koji prati promenu krutosti štapa u neelastičnoj oblasti i funkcija je nivoa opterećenja u štapu. Za potrebe ove analize formiran je deo računarskog programa ALIN koji može da se koristi za elastičnu i elasto-plastičnu analizu stabilnosti okvirnih konstrukcija. Program je napisan u C++ programskom jeziku. Primenom ovog programa omogućeno je i određivanje kritičnog opterećenja okvirnih nosača u elastičnoj i neelastičnoj oblasti. U ovom istraživanju formiran je i algoritam za proračun dužina izvijanja pritisnutih štapova stubova okvirnih nosača, a koji se bazira na proračunu globalne analize stabilnosti okvirne konstrukcije. Rezultati dobijeni primenom ovog algoritma upoređeni su s rešenjima koja se dobijaju korišćenjem evropskih EC3 i domaćih JUS standarda za okvirne čelične konstrukcije, a koja su približnog karaktera. Na osnovu postupka koji je dat u ovom radu moguće je praćenje fenomena gubitka stabilnosti okvirnog nosača u plastičnoj oblasti i direktno određivanje njegove kritične sile u toj oblasti.