Determination of the reinforced concrete slabs ultimate load using (cid:12)nite element method and mathematical programming

In the present paper, the ultimate load of the reinforced concrete slabs [16] is determined using the ﬁnite element method and mathematical programming. The acting eﬀorts and displacements in the slab are obtained by a perfect elasto-plastic analysis developed by ﬁnite element method. In the perfect elasto-plastic analysis the Newton-Raphson method [20] is used to solve the equilibrium equations at the global level of the structure. The relations of the plasticity theory [18] are resolved at local level. The return mapping problem in the perfect elasto-plastic analysis is formulated as a problem of mathematical programming [12]. The Feasible Arch Interior Points Algorithm proposed by Herskovits [8] is used as a return mapping algorithm in the perfect elasto-plastic analysis. The proposed algorithm uses Newton’s method for solving nonlinear equations obtained from the Karush-Kuhn-Tucker conditions [11] of the mathematical programming problem. At the end of this paper, it is analyzed six reinforced concrete slabs and the results are compared with available ones in literature.


INTRODUCTION
Reinforced concrete slabs [16] are among the most common structural elements.Despite of the large number of slabs designed and built, the details of their elastic and plastic behavior are not fully appreciated or properly taken into account.Although almost all the technical standards approach to slab design is basically one of using elastic moment distributions, it is also possible to design slabs using plastic analyses (perfect elasto-plastic analyses) to provide the required moments.In this paper, the ultimate load of the reinforced concrete slabs and the ultimate moment distributions are determined using the finite element method [20] and mathematical programming [12].
The resistance of a point within a structural element of reinforced concrete subjected to a multiaxial stress state depends on the interaction between the stresses which is subject.In the Trial stress tensor; slab, the strength criterion is defined for each midsurface point in terms of bending moments and torsion moment.With the strength criterion at each point and the elastic properties of the material is possible to perform a model analysis to determine the acting efforts from the applied load considering the elasto-plastic behavior of the material.The strength criterion proposed by Johansen [10] is used in this paper.
In the present paper, the acting efforts and displacements in the slab are obtained by a perfect elasto-plastic analysis developed by finite element method.The perfect elasto-plastic analysis of the slabs, described by their midsurface and discretized by the finite element method, is performed under the hypothesis of small displacements with consistent formulation in displacements.In the perfect elasto-plastic analysis the Newton-Raphson method [20] is used to solve the equilibrium equations at the global level of the structure.The relations of the plasticity theory [18] are resolved at local level, that is, for each Gauss point of the discretized structure.The return mapping problem in the perfect elasto-plastic analysis is formulated as a problem of mathematical programming [12].The Feasible Arch Interior Points Algorithm proposed by Herskovits [8] is used as a return mapping algorithm in the perfect elasto-plastic analysis.
The Feasible Arc Interior Point Algorithm [8] proposes to solve mathematical program-ming problems [12] with nonlinear objective function and nonlinear constraints quickly and efficiently.
The Feasible Arc Interior Point Algorithm is a new technique for nonlinear inequality and equality constrained optimization and was first developed by Herskovits [8].
This algorithm requires an initial point at the interior of the inequality constraints, and generates a sequence of interior points.When the problem has only inequality constraints, the objective function is reduced at each iteration.An auxiliary potential function is employed when there are also equality constraints.
The fact of giving interior points, even when the constraints are nonlinear, makes of proposed algorithm an efficient tool for engineering design optimization, where functions evaluation is in general very expensive.
Since any intermediate design can be employed, the iterations can be stopped when the objective reduction per iteration becomes small enough.
At each point, the Feasible Arc Interior Point Algorithm defines a "feasible descent arc".Then, it finds on the arc a new interior point with a lower objective.
The proposed algorithm uses Newton's method for solving nonlinear equations obtained from the Karush-Kuhn-Tucker conditions [11] of the mathematical programming problem.
The proposed algorithm requires at each iteration a constrained line search looking for a step-length corresponding to a feasible point with a lower objective.Herein we have implemented the Armijo's line search technique [8].
The implementation of the Feasible Arc Interior Point Algorithm was developed using the programming language C++ [19] that uses the technique of object-oriented programming.This technique allows quickly and located implementation of the proposed methods and also facilitates the code expansion.
In this paper will be presented: the strength criterion proposed by Johansen, the elastoplastic analysis of plates using finite element method and mathematical programming, the Feasible Arc Interior Point Algorithm and six examples of reinforced concrete slabs whose results are compared with results available in literature.

STRENGTH CRITERION
The strength criteria are characterized by a yield surface, defined as the geometric locus of the independent combinations of the stress tensor components or of the stress resultants that provoke the material plastification.Mathematically the yield surface can be defined by the Equation (1) presented as follows: The plasticity postulates define the yield surface as a continuous, convex region which could be regular or not.The yield surface implemented in this work was proposed by Johansen [2,10,14,16] and is of specific application for reinforced concrete slabs.

Strength criterion of Johansen
According to Johansen, the yield condition is based on the following physical criterion proposed by Massonet [14]: "The yield happens when the applied moment of flexion in a cross-section of inclination θ in relation to the x axis reaches a certain value that just depends on the angle θ and on the resistant moments in the reinforcement directions."The basic parameters of the Johansen criterion are presented in the Figure 1.The yield surface proposed by Johansen is frequently used to determine the ultimate resistance in the design of reinforced concrete slabs [10,14,16].The mathematical equations that define this surface are presented as follows.
In Equations ( 2) and (3), M + ux , M + uy , M − ux and M − uy are respectively the positive and negative ultimate moments of resistance.These ultimate moments of resistance are moments per unit of length in the x and y directions.
The Equation ( 2) is associated to the positive yield line [16] and Equation ( 3) is associated to the negative yield line [16].The Equations ( 2) and (3) represent two conical surfaces that combined define the yield surface of Johansen.The surface of Johansen is presented in the Figure 2.

ELASTO-PLASTIC ANALYSIS USING THE FINITE ELEMENT METHOD AND MATH-EMATICAL PROGRAMMING
The equations presented in this item are valid for materials with perfect elasto-plastic behavior.
In the determination of the efforts in a structure through a perfect elasto-plastic analysis is necessary to consider the plastic behavior of the material depending on the applied loading history.In this work the constitutive model applied to the reinforced concrete is perfect elastoplastic with associative flow rule.The condition that limits the stresses space [3] is presented as follows: Latin American Journal of Solids and Structures 9(2012) 69 -93 Figure 2 Yield surface of Johansen.
The Equation ( 4) represents a convex surface in the generalized stresses space.The interior region of this surface is formed by points belonging to the elastic regime.The plasticity theory allows only the existence of points in the interior (f < 0 ) or in the frontier of the yield surface (f = 0 ).Points placed out of this surface are inadmissible.
The perfect elasto-plasticity equations that govern the plastic behavior of the material after the yield are presented as follows: dλ ≥ 0 (8) The Equation (5) assumes that the strain increment tensor dε ε ε can be decomposed into an elastic part and a plastic part, indicated for dε ε ε e and dε ε ε p , respectively.The Equation (6) represents the incremental relation between stress and elastic strain where the stress increment tensor d σ is related with the elastic strain increment tensor dε ε ε e .The Equation (7) represents the dependence of the plastic strain increment tensor dε ε ε p with the associative flow rule.The Equation (8) represents that the consistency parameter is non-negative.The Equation (9) represents the complementarity condition that indicates that either the consistency parameter (d λ) or the strength criterion equation must be zero, so that the product of the two is vanished.The Equation (10) is the consistency condition which means that if dε ε ε p is different from zero, the stress state must "persist" in the surface, that is, d f (σ(t)) = 0.This expression is already represented by the Equation ( 9) for the perfect plasticity, once the yield surface does not change.
Analyzing the Equations ( 9) and (10) we verified the existence of three possible loading conditions that are indicated in the Table 1.
Table 1 Loading conditions in the perfect plasticity.
In this work the elasto-plastic consistent tangent stiffness matrix is used [17,18] for the update of the elasto-plastic stiffness matrix with the goal of improving the convergence of the Newton-Raphson method [20].

Used algorithms in nonlinear analyses with discretization by the Finite Element Method
The return mapping algorithms used in nonlinear physical analyses generally intend to solve the following problem: given a field of displacements, strains and stresses (u 0 , ε ε ε 0 and σ 0 ) of a structure, determine the new field of displacements, strains and stresses (u, ε and σ) due to a load increment.In this solution the compatibility and equilibrium relations defined in the solid mechanics and the constitutive relations defined in the plasticity theory are used.
The formulation generally adopted is to solve the equilibrium equations, deduced starting from the principle of virtual work and discretized by the finite element method [1,4], at global level of the structure.The relations of the plasticity theory are solved at local level, that is, for each Gauss point of the discretized structure.
The global equilibrium is obtained through algorithms that solve nonlinear equations systems and that determine the load-displacement curves [20].Among other algorithms, there are the Newton-Raphson algorithm, the displacement control and the arch length.In this work the algorithm of Newton-Raphson is used.
The return mapping algorithm intends to solve the following problem: given an initial state of allowable stress σ 0 , determine the new stress state σ due to a strain increment d ε, respecting the equations of the plasticity theory.The integration of the plasticity equations can be viewed as the solution of an initial value problem with constraints.This problem can be transformed in a mathematical programming problem with constraints through the use of an integration method of differential equations like implicit backward Euler [18].
The trial of a stress state σ e considering the behavior of the material in elastic regime, obtained starting from an allowable stress state σ 0 after a strain increment dε ε ε, is determined through the equation presented as follows: Applying the strength criterion: If f (σ e ) = 0, the assumption of linear-elastic regime is correct.There was not plastic strain and the total strain increment is elastic, dε ε ε = dε ε ε e .

Latin American Journal of Solids and Structures 9(2012) 69 -93
The trial elastic stress σ e is the solution stress state σ.If f (σ e ) > 0, the regime is perfect elasto-plastic.The solution stress state is on the yield surface, that is, f (σ) = 0 and it can be determined solving the mathematical programming problem presented as follows: Where C is the elastic stiffness modules matrix and f (σ) is the adopted strength criterion.The function s(σ) represents an ellipsoid in the stresses space.The problem consists of determining the smallest ellipsoid with center in σ e that touches f in σ.The graphic representation of this problem is presented in the Figure 3.The Lagrangian function [21] of the mathematical programming problem presented by the Equation ( 12) is presented as follows: The necessary 1st order conditions for the existence of a local minimum or Karush-Kuhn-Tucker conditions [11] are determined starting from the Equation (13) and are described as follows: Where f (σ) and ∂f ∂σ are appraised in the solution stress state σ.This point minimizes the mathematical programming problem indicated by the Equation (12).
Substituting the Equations ( 6) and (11) in the Equation ( 14) and after some algebraic manipulations, it is obtained the equation presented as follows: Latin American Journal of Solids and Structures 9(2012) 69 -93 An important observation to make is that the Karush-Kuhn-Tucker [11] equations of this mathematical programming problem correspond exactly to the perfect elasto-plasticity equations with associative flow rule developed previously.The Equation ( 16) corresponds to the non-negativity condition of the consistency parameter represented by the Equation ( 8).The Equation ( 18) corresponds to the associative flow rule represented by the Equation (7).The Equations ( 15) and ( 17) correspond respectively to the strength criterion condition represented by the Equation ( 4) and to the complementarity condition represented by the Equation (9).
The solution of this mathematical programming problem [12] is the stress state σ and the Lagrange multipliers d λ.Therefore, starting from a viable stress state σ 0 , it is obtained the stress solution state σ and the Lagrange multipliers d λ after a strain increment dε ε ε.The solution of this problem satisfies the Karush-Kuhn-Tucker conditions and consequently the perfect elasto-plasticity equations.The plastic strain is also determined since that the Lagrange multipliers correspond to the consistency parameters.
In this work the strength criterion of Johansen is used.Using this criterion, the mathematical programming problem results in a nonlinear programming problem with constraints.Both the objective function and the constraints of this problem are nonlinear.For the solution of the mathematical programming problem represented by the Equation ( 12), the Feasible Arch Interior Points Algorithm is used [8].This algorithm uses the Karush-Kuhn-Tucker equations of the mathematical programming problem indicated by the Equation (12).The advantage of this algorithm in relation to the others is its efficiency [8] for solving directly the Karush-Kuhn-Tucker equations, thereby solving a system of nonlinear equations.In addition, the unbounded number of constraints can be used in this problem without the need of significant change in the computational code, facilitating the treatment of yield multi-surfaces.

FEASIBLE ARCH INTERIOR POINTS ALGORITHM
Feasible Arch Interior Points Algorithm [8] is an iterative algorithm to solve the nonlinear programming problem [12] Where The following assumptions on the problem are required: 1.The functions f(x), g(x) and h(x) are continuous in Ω, as well as their first derivatives.
At each point the proposed algorithm defines a "feasible descent arc".A search is then performed along this arc to get a new interior point with a lower potential function.
Latin American Journal of Solids and Structures 9(2012) 69 -93 We denote ∇g(x) ∈ R nxm and ∇h(x) ∈ R nxp the matrix of derivatives of g and h respectively and call λ ∈ R m and µ ∈ R p the corresponding vectors of Lagrange multipliers.G(x) denotes a diagonal matrix such that G ii (x) = g i (x).The Lagrangian is presented as follows: The Hessian of the Lagrangian is presented as follows: Let us consider Karush-Kuhn-Tucker, (KKT), first order optimality conditions: A point x * is a stationary point if there exists λ * and µ * such that the Equations ( 22), ( 23) and ( 24) are true and is a KKT Point if KKT conditions (Equations ( 22), ( 23), ( 24), ( 25) and (26)) hold.
KKT conditions constitute a nonlinear system of equations and inequations on the unknowns (x,λ,µ).It can be solved by computing the set of solutions of the nonlinear system of Equations ( 22), ( 23) and (24) and then, looking for those solutions such that Equations (25) and (26) are true.However, this procedure is useless in practice.
The proposed algorithm makes Newton-like iterations to solve the nonlinear Equations ( 22), (23) and (24) in the primal and the dual variables.With the object of ensuring convergence to KKT points, the system is solved in such a way as to have the inequalities Equations ( 25) and (26) satisfied at each iteration.
Iterative methods for nonlinear problems in general include a local search procedure to force global convergence to a solution of the problem.This is the case of line search and trust region algorithms for nonlinear optimization [12,15].The present method includes a line search procedure, in the space of the primal variables x only, that enforces the new iterate to be closer from the solution.
Let d 0 ∈ R n such that d 0 = x 0 -x.From Equation (27), we have Which is independent of the current value of µ.Then Equation (28) gives a direction in the space of primal variables x and new estimates of the Lagrange multipliers.
Let the potential function be Where, at the iteration k, Where sg(.) = (.)/|(.)|.d k 0 is a descent direction of ϕ(c k ,x) [5,7].However, d 0 is not useful as a search direction since it is not necessarily feasible.This is due to the fact that as any constraint goes to zero, d 0 goes to a direction tangent to the feasible set [7].
To obtain a feasible direction, a negative vector -ρλ is added in the right hand side of the Equation (28).A perturbed linear system in d and λ is then obtained: Where ρ ∈ R is positive.The new direction is d and λ and μ are the new estimate of the Lagrange multipliers.We have now that d is a feasible direction, since ∇g t (x) d = −ρ < 0 for the active constraints.
The addition of a negative number in the right hand side of the Equation (28) produces a deflection on d 0 , proportional to ρ, in the sense of the interior of the feasible region.To ensure that d is also a descent direction, we establish an upper bound on ρ in order to have With α ∈ (0,1), that implies d t ∇ϕ (c, x) < 0. Thus, d is a descent direction of the potential function.
In general, the rate of descent of ϕ along d will be smaller than along d 0 .This is a price that we pay to get a feasible descent direction.
To obtain the upper bound on ρ, we solve the auxiliary linear system in (d 1 ,λ 1 ,µ 1 ).
Latin American Journal of Solids and Structures 9(2012) 69 -93 It follows from Equations ( 28), ( 31) and (33) that d = d 0 + ρd 1 .Then, we have that Equation (32) is true for any ρ > 0, if d t 1 ∇ϕ (c, x) < 0. Otherwise, we take And Equation (32) holds.However, when there are highly nonlinear constraints, the length of the feasible segment supported by the feasible descent search direction d can be not enough to accept a step equal to one.
The basic idea to avoid this problem consists on making the line search along a second order arc, tangent to the feasible descent direction d and with a curvature "close" to the curvature of the feasible set boundary.
The arc at x is defined by the following expression: Where d is obtained by solving: This linear system is similar to the system presented in Equation (33) being wI i and wE i computed as follows: The arc employed in the proposed algorithm is represented in Figure 4 for the case when the constraint g i (x) = 0 is active at the iterate x k .Since d k 0 and d k are descent directions of the potential function ϕ(x) at x k , their angle with -∇ϕ(x k ) is acute.

ARMIJO'S LINE SEARCH
The Feasible Arc Interior Point Algorithm requires at each iteration a constrained line search looking for a step-length corresponding to a feasible point with a lower potential.The first idea consists on solving the following optimization problem on t: Latin American Journal of Solids and Structures 9(2012) 69 -93 Figure 4 The feasible arc.
Instead of making an exact minimization on t, it is more efficient to employ inexact line search techniques [6,9,12,15].In this work the Armijo's line search technique was implemented.
Armijo's line search defines a procedure to find a step-length ensuring a reasonable decrease of the potential function.In our case, we add the condition of feasibility of the inequality constrains.Armijo's line search is stated as follows: Define the step-length t as the first number of the sequence {1,ν, ν 2 , ν 3 ,. . .}, satisfying Where η 1 ∈(0,1) and ν∈(0,1).

STATEMENT OF THE FEASIBLE ARC INTERIOR POINT ALGORITHM
The pseudocode of the Feasible Arc Interior Point Algorithm is presented on the next page.The present algorithm is very general in the sense that it converges to a Karush-Kuhn-Tucker point of the problem for any initial interior point.We work with the following updating rule for λ.
Data.Initial values for x ∈ R n such that g(x) < 0, λ ∈ R m , λ > 0, S ∈ R nxm symmetric and positive definite and c ∈ R p , c ≥ 0.
Step 1. Computation of a feasible descent direction.(i) Solve the linear systems: Let the potential function be Otherwise, set (iv) Compute the feasible descent direction: d=d 0 +ρd 1 Step 2. Computation of a "restoring direction" d Compute: Solve: Step 3. Arc search -Employ a line search procedure to get a step-length t based on the potential function ϕ c (x + td + t 2 d).In this work the Armijo's line search technique was implemented.
Step 4. Updates (i) Set the new point: The parameters ϵ, ḡ and λ I are taken positive.

EXAMPLE -DETERMINATION OF THE REINFORCED CONCRETE SLABS ULTI-MATE LOAD
In the determination of the slabs ultimate load using the finite element method and mathematical programming we used the following values.
In this paper, it was used the strength criterion of Johansen and also the finite element software FEMOOP [13].The Feasible Arc Interior Point Algorithm was implemented in FEMOOP using the programming language C++ [19] that uses the technique of object-oriented programming.This technique allows quickly and located implementation of the proposed methods and also facilitates the code expansion.

Square simply supported slab
In this example the ultimate load of the square simply supported slab on all edges with uniformly distributed load is determined.The slab is solid concrete.The span of the slab is l=5m.Figure 5 presents the reinforced concrete slab.The slab is isotropically reinforced with ultimate positive moments of resistance per unit width presented as follows: According to the yield line theory the ultimate load [16] of the slab without considering the corner effects is: Latin American Journal of Solids and Structures 9(2012) 69 -93 Considering the corner effects and assuming that the yield line pattern is composed of corners levers in the form of circular fans the ultimate load [16] is: Considering the corner effects and assuming that the yield line pattern is composed of corners levers in the form of hyperbolic fans the ultimate load [16] is: In the perfect elasto-plastic analysis [18], it was used isoparametric elements with eight nodes, Q8 [4].The used mesh was a bilinear-quadrilateral mesh with 20 elements in x direction and 20 elements in y direction.The integration order used was 2x2.
Using the perfect elasto-plastic analysis, it is obtained the ultimate load presented as follows: q u = 24.0KN / m 2 (perf ect elasto − plastic analysis) (54) Figure 6 presents the distribution of the principal moments (M 2 ) in the ultimate configuration.

Rectangular simply supported slab
In this example the ultimate load of the rectangular simply supported slab on all edges with uniformly distributed load is determined.The slab is solid concrete.The spans of the slab in the direction x is l x =7m and in the direction y is l y =5m. Figure 7 presents the reinforced concrete slab.
The slab is isotropically reinforced with ultimate positive moments of resistance per unit width presented as follows: Latin American Journal of Solids and Structures 9(2012) 69 -93 Figure 7 Reinforced concrete slab.
According to the yield line theory the ultimate load [16] of the slab is: In the perfect elasto-plastic analysis [18], it was used isoparametric elements with eight nodes, Q8 [4].The used mesh was a bilinear-quadrilateral mesh with 28 elements in x direction and 20 elements in y direction.The integration order used was 2×2.
Using the perfect elasto-plastic analysis, it is obtained the ultimate load presented as follows: Figure 8 presents the distribution of the principal moments (M 2 ) in the ultimate configuration.

Hexagonal slab fixed around the edges
In this example the ultimate load of the hexagonal slab fixed around the edges with uniformly distributed load is determined.The slab is solid concrete.The length of each side is l=5m.The inclined sides have an inclination of the 45 ○ with respect to the axis x. Figure 9 presents the reinforced concrete slab.
The slab is isotropically reinforced in the top and in the bottom with ultimate positive and negative moments of resistance per unit width presented as follows: According to the yield line theory the ultimate load [16] of the slab is:  In the perfect elasto-plastic analysis [18], it was used DKT elements [4].The used mesh was a triangular mesh with 554 elements.The integration order used was 2×2.
Using the perfect elasto-plastic analysis, it is obtained the ultimate load presented as follows: Figure 10 presents the distribution of the principals moments (M 2 ) in the ultimate configuration.

Rectangular slab bridge
In this example the ultimate concentrated load of the rectangular slab bridge is determined.The concentrated load is acting alone anywhere on the transverse centerline at midspan.The self-weight of the slab is neglecting.The slab is simply supported at two opposite edges and is free at the remaining two edges.The slab is solid concrete.The spans of the slab in the  The slab is isotropically reinforced in the top and in the bottom with ultimate positive and negative moments of resistance per unit width presented as follows: There are a number of possible yield line patterns, the critical pattern depending on the aspect ratio of the slab and the position of the load on the transverse centerline.The coefficient λ defines the position of the load on the transverse centerline.In this example is adopted λ=0.25. Figure 12 presents all possible yield line patterns.
According to the yield line theory the ultimate load [16] of the slab for mode 1 is: According to the yield line theory the ultimate load [16] of the slab for mode 2a is: and 28 elements in y direction.The integration order used was 2×2.
Using the perfect elasto-plastic analysis, it is obtained the ultimate load presented as follows: P u = 140.0KN(perf ect elasto − plastic analysis) (67) Figure 13 presents the distribution of the principals moments (M 2 ) in the ultimate configuration.

Rectangular slab with three edges supported and one edge free
In this example the ultimate load of the rectangular slab with three edges supported and one edge free with uniformly distributed load is determined.The slab is solid concrete.The spans of the slab in the direction x is l x =5m and in the direction y is l y =8m. Figure 14 presents the reinforced concrete slab.Where: In the perfect elasto-plastic analysis [18], it was used isoparametric elements with eight nodes, Q8 [4].The used mesh was a bilinear-quadrilateral mesh with 20 elements in x direction and 32 elements in y direction.The integration order used was 2×2.
Using the perfect elasto-plastic analysis, it is obtained the ultimate load presented as follows: q u = 21.9KN / m 2 (perf ect elasto − plastic analysis) (73) Figure 16 presents the distribution of the principals moments (M 2 ) in the ultimate configuration.

Square slab with openings
In this example the ultimate load of the square uniformly loaded slab with a central square opening is determined.The slab is fixed around the outside edges.The span of the slab is l=5m.The size of the opening is defined by the value of k.In this example we use k=0.2.The slab is solid concrete.Figure 17 presents the reinforced concrete slab.
The slab is isotropically reinforced with positive and negative ultimate moments of resistance per unit width presented as follows: Latin American Journal of Solids and Structures 9(2012) 69 -93

CONCLUSIONS
The values of ultimate loads of the examples presented above are presented in Table 2.This table shows the values using the perfect elasto-plastic analysis [18] and the yield line theory [16].The percentage error between these two values is also presented.In all the examples presented in this paper the stress distribution in the ultimate configuration determined using the perfect elasto-plastic analysis is according to the collapse mechanism predicted by the yield line theory.In example 7.1, the corner of the slab was held down and sufficient top steel was provided to avoid the appearance of the corner effects.The ultimate load found in the perfect elastoplastic analysis and the ultimate load predicted by the yield line theory both have the same value.
In examples 7.2, 7.3 and 7.4, the ultimate load found in the perfect elasto-plastic analysis and the ultimate load predicted by the yield line theory both have the same value.In example 7.5, the percentage error between the ultimate load found in the perfect elastoplastic analysis and the ultimate load predicted by the yield line theory is 4.4.The value provided by perfect elasto-plastic analysis is in favor of safety.In example 7.6, the percentage error between the ultimate load found in the perfect elastoplastic analysis and the ultimate load predicted by the yield line theory is 4.592.The value provided by perfect elasto-plastic analysis is in favor of safety.
Taking into account the previous results, we can conclude that the values using the perfect elasto-plastic analysis are very close to the values predicted by the yield line theory.Due to the use of the Feasible Arc Interior Point Algorithm [8] the computational cost of the analyses of the reinforced concrete slabs presented previously was not high.The perfect elasto-plastic analysis [18] allows the determination of the stresses and displacements at each gauss point in all loading stages.The yield line theory [16] does not allow to obtain these values.Therefore, we can assert that the tool developed is efficient and robust to determine the ultimate load of reinforced concrete slabs.

Figure 1
Figure 1 Basic parameters of the Johansen criterion.

Figure 3
Figure 3 Graphic representation of the return mapping algorithm.
ii) Define new values for λ>0 and S symmetric and positive definite.(iii) Go back to Step 1.

Figure 6
Figure 6 Distribution of the principal moments.

Figure 8
Figure 8 Distribution of the principal moments.

Figure 10
Figure 10 Distribution of the principal moments.

Figure 13
Figure 13 Distribution of the principal moments.

Figure 16
Figure 16 Distribution of the principal moments.