Finite element lower bound “yield line” analysis of isotropic slabs using rotation-free elements

• Published in: Engineering structures 2013-08, Vol.53, p.38-51
• Main Authors: Salam Al-Sabah, A., Falter, Holger
• Format: Article
• Language: English
• Subjects: Applied sciences; Building structure; Buildings; Buildings. Public works; Computing time; Construction (buildings and works); Convergence; Curvature; Exact sciences and technology; External envelopes; Failure; Finite element; Finite element method; Floor. Ceiling; Isotropic slab; Lower bounds; Mathematical analysis; Reinforced concrete; Reinforced concrete structure; Rotation-free element; Slabs; Stresses. Safety; Structural analysis. Stresses; Yield line analysis
• ISSN: 0141-0296; 1873-7323
• DOI: 10.1016/j.engstruct.2013.03.005

•Non-linear analysis of isotropic slabs was made using rotation-free elements.•Yield lines can pass through the element when yield condition is satisfied.•Element stiffness was modifies to account for internal yield lines.•New lower bound method is simpler than full non-linear finite element solution.•New method was used to solve four slab examples with reasonable accuracy. A new lower bound finite element method for slab analysis is presented as a practical substitute to full, non-linear, finite element methods that require expert knowledge and long running times. The method provides a general, safe and efficient lower bound solution for the analysis of reinforced concrete slabs up to failure. As it is finite element based, the method is more general than the yield line and strip methods currently in use. Furthermore, its lower bound nature makes it safer than the yield line method. The method uses a rotation-free, plate finite element modified to allow plastic “yield lines” to pass through at any direction. Yield lines are generated at the principal moment directions when the plastic moment capacity is attained. The material is assumed to be elastic perfectly-plastic. Following the general spirit of yield line analysis, the effects of a yield line are projected to the sides of the triangular element and then used to calculate the bending curvatures. The method’s efficiency is achieved by using rotation-free plate elements with a single degree of freedom per node and by the incremental solution that does not require iterations. The method’s accuracy and convergence are assessed by comparing standard cases with known results. In all cases, results were close to the theoretical values with difference of less than 1%. It is also used to solve a practical sized flat slab problem in order to demonstrate the method’s efficiency, convergence, and speed.