A Simple Elasto-Plastic Iterative Method by Fem for the Analysis of Plane Articulated Truss: Case of a 3-Bar Truss ()
1. Introduction
Structures consisting of truss elements are widely used in many engineering fields [1]. These trusses structures are fast to build, lightweight and allow large spans. The design of trusses structures has been the subject of numerous studies [2] [3] [4]. Many methods have been developed for the dimensioning of trusses. For small structures, classical analytical methods such as the node method (Culmann method) and the Ritter method are often used [5] [6]. For large structures, the graphical method was one of the first adopted techniques.
As a result of the increasing performance of computers, finite element software is now widely used for the calculation of trusses. Using numerical simulation tools, some authors have worked on methods for structural optimization of trusses [7] [8] [9]. McKeown [7] presented a displacement optimization approach for the design of linear elastic trusses. The mechanical modeling of truss structures is simple under ordinary conditions to which the classical assumption of linear elasticity applies. However, modern applications of truss structures often require plasticity material analysis, which can accurately predict the structural responses. This observation is justified by structural failures, such as the collapse of the I-35W Bridge over the Mississippi River in Minneapolis [10] [11]. Fan and et al. [12] studied the static and dynamic elasto-plastic stability of cross-linked hulls, while Ferrari et al. [13] numerically studied the non-linear elasto-plastic response of the historical Pardermo truss bridge. All these approaches are strongly non-linear.
In the present work, we propose a numerical simulation approach by FEM for the plane articulated truss. Our modeling will be based on a linearized material law of isotropic plastic strain-hardening (linearization of the elasto-plastic behavioral law: bilinear law) and on a projection method using the modified Newton-Raphson approach.
2. Simplified Behaviour Model
Selecting a Template
In this study, we are interested in the elasto-plastic behaviour of articulated truss structures. A simplified linear isotropic strain-hardening model of the Johson Cook type [14] is implemented. Thus, the simplified model is bilinear and characterized by a Young’s modulus E and a tangent modulus ET, the yield strength
, yield stress projection
(see Figure 1). The material constituting our truss is made of a metal alloy, therefore the actual behaviour law is shown in Figure 1(a). The corresponding simplified model is shown in Figure 1(b). In this study, we do not consider the phenomena of damage and failure [15].
This approximation of the real material behaviour law (Figure 1(a)) by a bilinear behaviour law of elasto-plastic type with linear strain hardening (Figure 1(b)) permits to write the following incremental plastic flow law (Equation (1)):
(1)
With
,
: elastic deformation increment;
: plastic deformation increment.
We deduce the evolution law of linear strain-hardening (Equation (2)):
(2)
With
, H is called the strain-hardening modulus,
: Current
Figure 1. Tensile stress-strain curves. (a) Actual experimental curve, (b) Linear approximation of strain-hardening (simplified model).
yield strength.
Equation (2) can also be expressed as a
function as follows (Equation (3)):
(3)
The linearized elasto-plastic evolution allows a simplified writing of the discretized elementary equilibrium, obtained analytically by the FEM an articulated truss (Equation (4)):
(4)
: local stiffness matrix,
: strain-displacement relationship matrix,
: Nodal displacement field.
3. Elasto-Plastic Numerical Resolution
The numerical solving technique used in our elasto-plastic problem is based on the classical incremental, iterative approach and on the minimization of a convergence criterion. The modified Newton Raphson-type projection method, which uses the elastic stiffness matrix at each iteration. Figure 2 presents the algorithm followed in MATLAB program [2].
Figure 2. Summary of numerical resolution algorithm.
4. Projection Method
The projection method on the solution sought is the modified Newton-Raphson method. The linear approximation of the plasticity law allows an easy calculation of the hardening, thus saving the CPU processing time. Figure 3 shows the principle of projection on the boundary of the current elasticity domain for an increase in dε deformation. We have represented the most general case of an element being plasticized. We consider a history of initial stress equal to σ and the state of current strain hardening is denoted σh.
The linear approximation of the strain hardening allows to obtain by a simple calculation, the plastic evolution law. From Figure 3(b), are geometrically deduced the relations presented in equations 5 to 8. First, we define a quantity R, as follows (Equation (5)):
(5)
We deduce from the Figure 3(b) (Equation (6)):
(6)
The law of evolution (Equation (7)):
(7)
The increments of elastic and plastic deformations are deduced from the stress increment as follows (Equation (8)):
Figure 3. Representation of the projection method [2].
(8)
The numerical implementation is done through the projection algorithm of Figure 4.
5. Case of Application
The application case structure in this article is a truss consisting of three (03) bars of identical cross-section S (see Figure 5). The inclined bar is of length
and the other of length L. The material behaviour law described is the one presented in Section 2. The numerical values considered are: L = 500 mm, Young’s modulus E= 70.000 MPa, tangent modulus ET = 510.8 MPa, yield strength σ0 = 281.559 MPa, bar section S = 50 mm2. Table 1 shows the table of connectivities of the 3-bar truss structure.
The resolution by FEM involves either a matrix or a discrete formulation of the problem. In the displacement approach, the main unknowns of the problem are the nodal displacements of the structure. In this article, the stiffness matrix is used to construct the FEM model. The elementary stiffness matrix (for each bar element) is written. In the case of the articulated trusses we consider 1 nodal degree of freedom in the local reference and 2 nodal degrees of freedom in the global reference.
The global stiffness matrix is obtained after assembling the elementary stiffness matrices. Thus, the main unknown, which is the nodal displacement vector is calculated by solving the global matrix equilibrium equation [K]{Q} = {F}. In this equation, [K] is the global stiffness matrix, {F} the nodal force vector and {Q} the nodal displacement vector. When the bar element is in the plastic behaviour phase, the local matrix equilibrium equation is modified as shown in Equation (4). Once the nodal displacement vector has been determined, the internal forces in the bars can be deduced. In the application case of 3-bar truss in this article, the global stiffness matrix is presented in Equation (9). A computational code is developed in MATLAB for the application case.
Figure 4. Projection algorithm (modified Newton-Raphson) [2].
(9)
6. Results and Validation
The results presented in this section are obtained from the finite element computation code developed under MATLAB and, having the architecture of the numerical resolution algorithm presented in Figure 2 and Figure 4. The solving method uses an incremental and iterative method based on minimizing a convergence criterion. Different values of the approximation error of the convergence criterion are used to study its impact on the algorithm used. The MATLAB program is illustrated in the case of the articulated truss shown in Figure 5. A load F = 35 kN is applied to the structure, the load increment dF is chosen to be equal to 100N.
6.1. Simulation and Validation
In this paragraph, examples of numerical simulation are presented in order to show the reliability of our algorithm to correctly simulate an elasto-plastic behavior of an articulated truss. Then, the validation is made by comparing with the results obtained from the ABAQUS/CAE computation code.
Tables 2-4 show the main results (N1, N2, N3: internal forces in each bar element and V1: vertical displacement of node 1, see Figure 5) for some load increments obtained from the calculation code developed in MATLAB. The results are presented in each table according to the selected value of the approximation error ξ of the convergence test. The ξ values shown in the tables are: 1, 0.01, and 0.00001.
The loads chosen as examples in these different tables correspond to the 3 phases of behaviour of the structure: elastic structure (elastic bars), partial plastic structure (bar 2 in the plastic domain), plastic structure (all bars in the plastic domain).
Table 2. Internal forces (N1, N2 and N3) in the bars and vertical displacement V1 (of node 1) for different loads (approximation error ξ = 1).
Table 3. Internal forces (N1, N2 and N3) in the bars and vertical displacement V1 (of node 1) for different loads (approximation error ξ = 0.01).
Table 4. Internal forces (N1, N2 and N3) in the bars and vertical displacement V1 (of node 1) for different loads (approximation error ξ = 0.00001).
The results presented in Tables 2-4 are compared with the results obtained using the ABAQUS/CAE calculation code. For example, we choose two load cases of 20 kN and 30 kN corresponding to partial and total plasticization phases of the structure. Figure 6 shows the stress states in the truss for the two load cases obtained by ABAQUS. The internal forces in the bars are respectively the couple (N1 = N3 = 5858 N, N2 = 11715.5 N) and (N1 = N3 = 11216 N, N2 = 14138.5 N) for loads F of 20 kN and 30 kN respectively. Whatever the load considered, the results are quite close to ours (see the corresponding row in Tables 2-4).
Figure 7 shows the evolution of the internal forces as a function of the load. It can be seen that each curve obtained has 3 areas of linearity corresponding to the 3 phases of behaviour of the following truss structure: elastic phase (bars in the elastic), partial plastic phase (bar 2 in the plastic) and complete plasticization phase (bars 1, 2 and 3 in the plastic).
Figure 8 shows the evolution of the vertical displacement V1 of node 1 as a function of the load F. As in Figure 7, a tri-linear curve of the displacement V1 is observed. An almost exponential increase of the displacement can be seen when all the bars are plasticized. This is explained by the low value of the tangent modulus (ET ≈ 510 MPa) of the bilinear material behaviour law used (see Figure 1).
Figure 6. Stress mapping of the 3-bar truss from ABAQUS, for loads (a) F = 20 kN and (b) F = 30 kN.
Figure 7. Evolution of the internal forces in the each bar as a function of the load for ξ = 0.00001.
Figure 8. Displacement of the load application point as a function of the load F.
Figure 9. Influence of the approximation error ξ of the convergence test on the internal forces in the bars: (a) bar 1 (b) bar 2.
Figure 10. Influences of the approximation error ξ of the convergence test on the vertical displacement V1 of node 1.
6.2. Influence of the Value of Approximation Error x
In this section, we are interested in the influence of the ξ value on the final result obtained. Figure 9 and Figure 10 illustrate respectively the influence of the ξ value on the solution of the internal forces in the bars and the vertical displacement V1 of node 1. The illustrative load cases are: 5000, 24,400 and 34,600 N corresponding to the 3 phases of behaviour of the structure. Whatever the load level, it can be seen that the value of ξ has little influence on the final result. Nevertheless, for high loads (F = 34,600 N), a non-negligible influence of the ξ value on the displacement V1 is observed (Figure 10). This weak influence is explained by the bilinear nature of the implemented behaviour law in our study.
7. Conclusion
This article presented the development of a FEM calculation code for the simple calculation of an articulated truss structure. A finite element program was developed in the MATLAB environment for a three-bar truss structure. The special feature of the numerical resolution of our computational code is the use of a modified Newton-Raphson projection method, by using a linearized elastoplastic material law. Simulations on a 3-bar articulated truss structure and a comparison with a similar model under ABAQUS have proven the reliability of our approach.