Services on Demand
Article
Indicators
 Cited by SciELO
 Access statistics
Related links
 Similars in SciELO
Share
Latin American Journal of Solids and Structures
Online version ISSN 16797825
Lat. Am. j. solids struct. vol.9 no.1 Rio de Janeiro 2012
http://dx.doi.org/10.1590/S167978252012000100004
ARTICLES
Determination of the reinforced concrete slabs ultimate load using finite element method and mathematical programming
A.M. Mont'Alverne^{I,} ^{*}; E.P. Deus^{II}; S.C. Oliveira Junior^{III}; A.S. Moura^{IV}
^{I}Center of Technology, Federal University of Ceará, Pici Campus, 60455900, Fortaleza – Brazil
^{II}Deptt. of Metallurgical and Materials Engg., Center of Technology, Federal University of Ceará, Fortaleza – Brazil
^{III}PostDoctoral Researcher of ARMTEC (CNPq/PNPD), Fortaleza, Ceará – Brazil
^{IV}Deptt. of Hydraulic and Environmental Engg., Center of Technology, Federal University of Ceará – Brazil
ABSTRACT
In the present paper, the ultimate load of the reinforced concrete slabs [16] is determined using the finite element method and mathematical programming. The acting efforts and displacements in the slab are obtained by a perfect elastoplastic analysis developed by finite element method. In the perfect elastoplastic analysis the NewtonRaphson 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 elastoplastic 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 elastoplastic analysis. The proposed algorithm uses Newton's method for solving nonlinear equations obtained from the KarushKuhnTucker 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.
Keywords: optimization, Finite Element Method, plates theory, reinforced concrete.
NOTATIONS
C  Elastic stiffness modules matrix; 
dε  Strain increment tensor; 
dλ  Consistency parameter or Lagrange multiplier; 
dσ  Stress increment tensor; 
E  Young modulus; 
f  Yield surface or strength criterion; 
f  Objective function; 
g  Equality constraint; 
h  Inequality constraint or slab height; 
L  Lagrangian function; 
l  Span; 
Positive and negative ultimate moment of resistance;  
P_{u}  Ultimate concentrated load; 
q_{u}  Ultimate Uniformly distributed load; 
m  Acting moment; 
u  Displacement tensor; 
ε  Strain tensor; 
λ  Lagrange multiplier; 
µ  Lagrange multiplier; 
ν  Poisson ratio; 
σ  Stress tensor; 
σ_{e}  Trial stress tensor; 
1 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 elastoplastic 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 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 elastoplastic 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 elastoplastic analysis developed by finite element method. The perfect elastoplastic 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 elastoplastic analysis the NewtonRaphson 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 elastoplastic 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 elastoplastic analysis.
The Feasible Arc Interior Point Algorithm [12] proposes to solve mathematical programming problems [] 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 KarushKuhnTucker conditions [11] of the mathematical programming problem.
The proposed algorithm requires at each iteration a constrained line search looking for a steplength 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 objectoriented 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.
2 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.
2.1 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 crosssection 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), and 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.
3 ELASTOPLASTIC ANALYSIS USING THE FINITE ELEMENT METHOD AND MATHEMATICAL PROGRAMMING
The equations presented in this item are valid for materials with perfect elastoplastic behavior. In the determination of the efforts in a structure through a perfect elastoplastic 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:
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 elastoplasticity equations that govern the plastic behavior of the material after the yield are presented as follows:
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 nonnegative. 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, df(σ(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.
In this work the elastoplastic consistent tangent stiffness matrix is used [17, 18] for the update of the elastoplastic stiffness matrix with the goal of improving the convergence of the NewtonRaphson method [20].
3.1 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 loaddisplacement curves [20]. Among other algorithms, there are the NewtonRaphson algorithm, the displacement control and the arch length. In this work the algorithm of NewtonRaphson 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 linearelastic regime is correct. There was not plastic strain and the total strain increment is elastic, dε = dε^{e}. The trial elastic stress σ_{e} is the solution stress state σ. If f(σ_{e}) > 0, the regime is perfect elastoplastic. 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 KarushKuhnTucker conditions [11] are determined starting from the Equation (13) and are described as follows:
Where f(σ) and 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:
An important observation to make is that the KarushKuhnTucker [11] equations of this mathematical programming problem correspond exactly to the perfect elastoplasticity equations with associative flow rule developed previously. The Equation (16) corresponds to the nonnegativity 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 KarushKuhnTucker conditions and consequently the perfect elastoplasticity 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 KarushKuhnTucker 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 KarushKuhnTucker 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 multisurfaces.
4 FEASIBLE ARCH INTERIOR POINTS ALGORITHM
Feasible Arch Interior Points Algorithm [8] is an iterative algorithm to solve the nonlinear programming problem [12]
Where x ∈ ℜ^{n} and f ∈ ℜ, g ∈ ℜ^{m}, h ∈ ℜ^{p}. Ω = {x ∈ ℜ^{n}/g(x) = 0}. 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.
2. For all x ∈ Ω the vectors ∇g_{i}(x), for i = 1,2,...,m such that g_{i}(x)=0 and ∇h_{i}(x), for i = 1,2,...,p are linearly independent.
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.
We denote ∇g(x) ∈ ℜ^{nxm} and ∇h(x) ∈ ℜ^{nxp} the matrix of derivatives of g and h respectively and call λ ∈ ℜ^{m} and µ ∈ ℜ^{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 KarushKuhnTucker, (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 Newtonlike 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.
Let S = Lh(x, λ, µ). A Newton iteration for the solution of the Equations (22), (23) and (24) is defined by the following linear system:
Where (x, λ, µ) is the current point and (x_{0}, λ_{0}, µ_{0}) is a new estimate. We call Λ = diag(λ). We can also take S = B, a quasiNewton approximation of Lh(x, λ, µ), or S = I (identity) [7].
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} ∈ ℜ^{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, is such that
Where sg(.) = (.)/ (.) . 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 ρ ∈ ℜ 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}).
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 ∇ϕ(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 is obtained by solving:
This linear system is similar to the system presented in Equation (33) being and 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 and d^{k} are descent directions of the potential function ϕ(x) at x^{k}, their angle with ∇ϕ(x^{k}) is acute.
5 ARMIJO'S LINE SEARCH
The Feasible Arc Interior Point Algorithm requires at each iteration a constrained line search looking for a steplength corresponding to a feasible point with a lower potential. The first idea consists on solving the following optimization problem on t:
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 steplength 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 steplength t as the first number of the sequence {1,ν, ν^{2}, ν^{3},... }, satisfying
And
Where η_{1} ∈ (0,1) and ν ∈ (0,1).
6 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 λ.
Set, for i=1,...,m.
If g_{i} (x) – and λ_{i} < λ^{I} , set λ_{i} = λ^{I}.
Parameters. α ∈ (0,1) and φ > 0.
Data. Initial values for x ∈ ℜ^{n} such that g(x) < 0, λ ∈ ℜ^{m}, λ > 0, S ∈ ℜ^{nxm} symmetric and positive definite and c ∈ ℜ^{p}, c > 0.
Step 1. Computation of a feasible descent direction.
(i) Solve the linear systems:
And
Let the potential function be
(ii) If c_{i} < 1.2µ_{0}(i), then set c_{i}=1.2µ_{0}(i); i=1,...,p
(iii) If ∇ϕ_{c} (x) > 0, set
Otherwise, set
(iv) Compute the feasible descent direction: d=d_{0 }+ ρd_{1}
Step 2. Computation of a "restoring direction"
Compute:
Solve:
Step 3. Arc search  Employ a line search procedure to get a steplength t based on the potential function ϕ_{c} (x + td + t^{2}). In this work the Armijo's line search technique was implemented.
Step 4. Updates
(i) Set the new point:
(ii) Define new values for λ > 0 and S symmetric and positive definite.
(iii) Go back to Step 1.
The parameters є, and λ^{I} are taken positive.
7 EXAMPLE  DETERMINATION OF THE REINFORCED CONCRETE SLABS ULTIMATE LOAD
In the determination of the slabs ultimate load using the finite element method and mathematical programming we used the following values.

Slab height: h=0.1 m;

Young modulus: E=23800 MPa;

Poisson ratio: ν=0.2.
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 objectoriented programming. This technique allows quickly and located implementation of the proposed methods and also facilitates the code expansion.
7.1 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:
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 elastoplastic analysis [18], it was used isoparametric elements with eight nodes, Q8 [4]. The used mesh was a bilinearquadrilateral mesh with 20 elements in x direction and 20 elements in y direction. The integration order used was 2x2.
Using the perfect elastoplastic analysis, it is obtained the ultimate load presented as follows:
Figure 6 presents the distribution of the principal moments (M_{2}) in the ultimate configuration.
7.2 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:
According to the yield line theory the ultimate load [16] of the slab is:
In the perfect elastoplastic analysis [18], it was used isoparametric elements with eight nodes, Q8 [4]. The used mesh was a bilinearquadrilateral mesh with 28 elements in x direction and 20 elements in y direction. The integration order used was 2×2.
Using the perfect elastoplastic 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.
7.3 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 elastoplastic 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 elastoplastic 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.
7.4 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 selfweight 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 direction x is l_{x}=5m and in the direction y is l_{y}=7m. Figure 11 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:
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:
According to the yield line theory the ultimate load [16] of the slab for mode 2b is:
According to the yield line theory the ultimate load [16] of the slab for mode 3a is:
According to the yield line theory the ultimate load [16] of the slab for mode 3b is:
In the perfect elastoplastic analysis [18], it was used isoparametric elements with eight nodes, Q8 [4]. The used mesh was a bilinearquadrilateral mesh with 20 elements in x direction and 28 elements in y direction. The integration order used was 2×2.
Using the perfect elastoplastic analysis, it is obtained the ultimate load presented as follows:
Figure 13 presents the distribution of the principals moments (M_{2}) in the ultimate configuration.
7.5 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.
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 two possible yield line patterns. The governing alternative collapse mode is the one giving the lowest ultimate load. Figure 15 presents the two possible yield line patterns.
According to the yield line theory the ultimate load [16] of the slab for mode 1 is:
Where:
According to the yield line theory the ultimate load [16] of the slab for mode 2 is:
Where:
In the perfect elastoplastic analysis [18], it was used isoparametric elements with eight nodes, Q8 [4]. The used mesh was a bilinearquadrilateral mesh with 20 elements in x direction and 32 elements in y direction. The integration order used was 2×2.
Using the perfect elastoplastic analysis, it is obtained the ultimate load presented as follows:
Figure 16 presents the distribution of the principals moments (M_{2}) in the ultimate configuration.
7.6 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:
According to the yield line theory the ultimate load [16] of the slab is:
In the perfect elastoplastic analysis [18], it was used isoparametric elements with eight nodes, Q8 [4]. The used mesh was a bilinearquadrilateral mesh with 384 elements. The integration order used was 2×2.
Using the perfect elastoplastic analysis, it is obtained the ultimate load presented as follows:
Figure 18 presents the distribution of the principals moments (M_{2}) in the ultimate configuration.
8 CONCLUSIONS
The values of ultimate loads of the examples presented above are presented in Table 2. This table shows the values using the perfect elastoplastic 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 elastoplastic 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 elastoplastic 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 elastoplastic 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 elastoplastic analysis is in favor of safety.
Taking into account the previous results, we can conclude that the values using the perfect elastoplastic 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 elastoplastic 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.
References
[1] K.J. Bathe. Finite Elements Procedures. KlausJurgen Bathe, 2007. [ Links ]
[2] M.W. Braestrup and M.P. Nielsen. Plastic methods of analysis and design. In F. K. Kong et al., editors, Handbook of Structural Concrete. Pitman Publishing, 1983. [ Links ]
[3] W.F. Chen. Plasticity in Reinforced Concrete. J. Ross Publishing, 2007. [ Links ]
[4] R.D. Cook, D.S. Malkus, M.E. Plesha, and R.J. Witt. Concepts and Applications of Finite Element Analysis. John Wiley & Sons, 4th edition, 2001. [ Links ]
[5] J. Herskovits. A twostage feasible directions algorithm for nonlinear constrained optimization, volume 36. Math. Program., 1986. [ Links ]
[6] J. Herskovits. A view on nonlinear optimization. In J. Herskovits, editor, Advances in Structural Optimization, pages 71117, Dordrecht, Holland, 1995. Kluwer Academic Publishers. [ Links ]
[7] J. Herskovits. A feasible directions interior point technique for nonlinear optimization. JOTA, J. Optimiz. Theory Appl., 99(1):121146, 1998. [ Links ]
[8] J. Herskovits, P. Mappa, E. Goulart, and C.M. Mota Soares. Mathematical programming models and algorithms for engineering design optimization. Computer Methods in Applied Mechanics and Engineering, 194(3033):32443268, 2005. [ Links ]
[9] J.B. HiriartUrruty and C. Lemaréchal. Convex Analysis and Minimization Algorithms I and II. SpringerVerlag, Berlin, Germany, 2010. [ Links ]
[10] K.W. Johansen. Yield Line Theory. Cement and Concrete Association, 1962. [ Links ]
[11] H.W. Kuhn and A.W. Tucker. Nonlinear programming. Proc. 2º Berkeley Symp. Math. Statistics and Probability, 481, 1950. [ Links ]
[12] D.G. Luenberger. Linear and Nonlinear Programming. AddisonWesley Publishing Company, 2nd edition, 1984. [ Links ]
[13] L.F. Martha, I.F.M Menezes, E.N. Lages, E. Parente Jr., and R.L.S. Pitangueira. An OOP class organization for materially nonlinear nite element analysis. In Proceedings of the XVII CILAMCE, pages 229232, Venice, Italy, 1996. [ Links ]
[14] C.E. Massonet and M.A. Save. Plastic Analysis and Design of Plates Shells and Disks. NorthHolland Publishing Company, 1972. [ Links ]
[15] J. Nocedal and S.J. Wright. Numerical Optimization. Springer Series in Operations Research. Springer, New York, 2nd edition, 2006. [ Links ]
[16] R. Park and W.L. Gamble. Reinforced Concrete Slabs. John Wiley & Sons, 2nd edition, 2000. [ Links ]
[17] J.C. Simo and T.J.R. Hughes. A return mapping algorithm for plane stress elastoplasticity. International Journal for Numerical Methods in Engineering, 5:649670, 1986. [ Links ]
[18] J.C. Simo and T.J.R. Hughes. Computational Inelasticity, volume 7 of Interdisciplinary Applied Mathematics. SpringerVerlag, 1998. [ Links ]
[19] B. Stroustrup. The C++ Programming Language: Special Edition. AddisonWesley Publishing Company, 3rd edition, 2000. [ Links ]
[20] O.C. Zienkiewicz R.L. Taylor. The Finite Element Method for Solid and Structural Mechanics. ButterworthHeinemann, 6th edition, 2005. [ Links ]
[21] G.N. Vanderplaats. Numerical Optimization Techniques for Engineering Design. Vanderplaats Research and Development, 3rd edition, 1999. [ Links ]
Received 26 May 2011;
In revised form 15 Sep 2011
* Author email: montalverne@ufc.br