A MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES

Size: px
Start display at page:

Download "A MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES"

Transcription

1 COMPUTATIONAL MECHANICS New Trends and Applications S. Idelsohn, E. Oñate and E. Dvorkin (Eds.) c CIMNE, Barcelona, Spain 1998 A MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES Paola Caracciolo and Emilio Turco Dipartimento di Strutture, Università della Calabria Rende (Cs), Italy emilio@virgo.labmec.unical.it Key words: Mixed finite elements, Elasto-plastic problems, Kirchhoff plates. Abstract. A triangular mixed finite element model for the elasto-plastic analysis of Kirchhoff plates is presented. The proposed discretization is based on rigid in bending elements connected by means of rotational springs between the inter-element sides and constant bending moments on the influence area of each node. The continuity of the normal bending moment through the sides allows a standard Hellinger-Reissner formulation. The advantages are related to a reasonable accuracy and a very simple element algebra. Some numerical results, obtained by using extremal path theory and the arc-length strategy, allow an analysis of the element s performance in the evaluation of the collapse load. 1

2 1 INTRODUCTION Finite elements based on high order interpolation, in spite of their algebraic complexity and the number of variables bound to the element construction, are suitable for linear elastic problems characterized by solutions with a high degree of continuity. If the analysis is extended in plasticity the discontinuities in the constitutive law require a totally different background. In particular, the localization of the plastic strains, when the load is near to the collapse load, leads to the formation of isolated sliding surfaces. This consideration suggests the use of simple finite elements for which the only way to obtain an accurate solution remains related to a mesh with a high number of elements. Furthermore, since the nonlinear behaviour of the material is strongly influenced by the stress level, compatible elements do not appear to be particularly convenient 1. This reason leads to a preference for mixed finite elements, since when using them, generally, a more accurate reconstruction of the stress field should be expected. This paper proposes a simple triangular mixed finite element developed to analyse elasto-plastic Kirchhoff plates. This element was suggested in Kawai s paper 2 and Argyris and his collaborators paper 3 and characterized by a linear interpolation of the displacement field between the nodes of the triangle making it rigid in bending. Therefore, the flexural deformation is concentrated on the sides of the triangle and measured by the relative rotation between the elements adjacent to the side considered. The further hypothesis of constant bending moments on the influence area of each node and the continuity of the normal bending moment through the sides allows a standard Hellinger-Reissner formulation. The above hypotheses produce an element which has a simple algebra and therefore make it suitable for discretization with many elements. In order to reconstruct the equilibrium path of the structure the load evolution law is discretized into a sequence of finite steps. Starting with a known equilibrium point the extremal path theory developed by Ponter and Martin 4 is used. This theory, in the case of perfect elasto-plasticity, leads to the Haar-Karman holonomic solution 5. Along with these assumptions the arc-length strategy is used to evaluate the next equilibrium point and therefore the algorithm keeps its accuracy even near the limit point which is one of the most important aims of the analysis. In the next sections the triangular mixed elements will be presented and discussed along with the incremental approach particularized for the underlying problem. In order to assess the performance of the proposed mixed element some numerical tests, referring to the elastic perfectly plastic behaviour, are performed and their results are compared with some classic results for plates on regular domains. 2 A SIMPLE TRIANGULAR MIXED FINITE ELEMENT The features of the finite element presented below derive essentially from the Kawai s paper 1 which suggests the use of rigid finite elements connected by means of rotational springs over the contact area of two neighbouring elements. A particular type of Kawai s 2

3 element was proposed by Argyris and his collaborators 3 in the simulation of car crash phenomena. This element is precious in its field of application but the elasto-plastic problem of Kirchhoff plates requires a less simplified element. Following Argyris paper 3 Figure 1: Geometric, kinematic and static variables for the triangular element the plate is subdivided into a suitable number of rigid in bending triangular elements connected with rotational springs on the sides of elements. Owing to the indeformability of elements, the bending deformation is measured by the change of angle between two elements which have a common side. Compatibility equations link the displacements of the element nodes to the change of angle along each side of the triangle. By referring to Figure 1 for the notation, for each element, in the hypothesis of small displacements, the side rotation for the i-th side is given by ϕ i = 1 2A e (w i l i w j l j cos α k w k l k cos α j ) (1) where A e represents the element area and w i, l i and α i are, respectively, the out of plane displacement, the length and the angle relative to the i-node. Two similar equations are obtained for the side j and k, by opportunely exchanging the indices. Therefore collecting side rotations ϕ i, ϕ j, ϕ k in vector ϕ e and nodal displacements w i, w j, w k in vector w e, the following relation is obtained ϕ e = D e w e (2) where D e represents the element compatibility matrix. Now, if vector ϕ collects the change of angle on each side of the mesh and w the out of plane displacements, the compatibility is ensured if ϕ = Dw (3) 3

4 where D is the global compatibility matrix having n l rows and n n columns, n l and n n are, respectively, the number of sides and the number of nodes of the mesh. This matrix is obtained assembling the contributions deriving from the element compatibility matrix D e. It is important to note that the side rotation contributions derive only from the elements adjacent to the side considered. When the choice of the kinematic variables is made, the static variables can be derived in such a way that they respect the duality of the virtual work theorem. For this reason, a vector f which collects the external nodal force and a bending moment vector m which collects the normal bending moment on the element sides are considered. The link between the static variables can be still derived from the virtual work theorem which gives equilibrium equations D T m = f (4) At this point it remains to specify the constitutive law of this discrete model. This can be done, for example, by defining the complementary elastic energy. Indeed, this energy referred to the surface S of the continuum problem can be written as Φ S = 1 m T E 1 mda (5) 2 S where m collects m x, m y and m xy in any points of the part S and E = Eh 3 12(1 ν 2 ) 1 ν 0 ν (1 ν) (6) where, following the usual notation, E represents the Young modulus, ν the Poisson ratio and h the plate thickness. Now, two possible ways can be followed introducing some hypotheses on the bending moment distribution. The first way assumes the bending moments m x, m y and m xy constant on the e-th element. This hypothesis define the link between the bending moments on the element sides in a unequivocal way. Indeed, looking at Figure 1 results m i = m x cos 2 β i + m y sin 2 β i +2m xy sin β i (7) where β i is the angle between x-axis and the outward normal to the i-side. Similar relations can be obtained for m j and m k. Collecting m i, m j and m k in vector m l rotation matrix R is automatically defined by means of the relation m l = R m (8) By using the hypothesis on constant bending moments on the element, the complementary elastic energy for the e-th element can be written as Φ e = 1 2 A e m T E 1 m (9) 4

5 or using the normal bending moments on sides of the e-th element Φ e = 1 2 A em T l (R 1 ) T E 1 R T m l (10) This approach uses as variables the vector of the nodal displacements w and the bending moments on the sides of the mesh collected in vector m. The second way assumes bending moments m x, m y and m xy as constants on the influence area of each node of the mesh, being the influence area of the i-th node defined as A i = 1 3 A e i (11) and the sum has to be intended as extended on all elements around the node i, A e i is the area of each of these elements. Furthermore, the normal bending moment on side l between the i-th and the j-th node of the mesh is defined as the average value m l = 1 2 RT l ( m i + m j ) (12) where R l is the column of the matrix R related to side l. Now, using the global equivalence relation between the deformation work written using thesidevariables,m l and ϕ l, and node variables m n and ϕ n which collect, respectively, m x, m y, m xy and ϕ x, ϕ y, ϕ xy on the n-node n l n n m l ϕ l = m T n ϕ n (13) l=1 n=1 and substituting in it the Eq.(12) n l l=1 ( ) 1 2 RT l ( m i + m j ) ϕ l = n n n=1 m T n ϕ n (14) the following relation is obtained once contributions related to each node are assembled n n ( ( )) 1 m T n n=1 2 RT l ϕ n n l = m T n ϕ n (15) n=1 The relation immediately above furnishes the definition of ϕ n and completes this second approach where the nodal displacements and the bending moments on nodes are assumed as independent variables. The two ways delineated above define the basic ingredient to generate a mixed finite element discretization. Indeed, by using the quantities defined above, the standard form of the Hellinger-Reissner functional can be written as π[w, m] = 1 2 mt Fm + m T Dw f T w (16) 5

6 where F is the flexibility matrix. It has to be pointed out that in the immediately above functional vectors and matrices have a different meaning for the two finite elements above-mentioned. Indeed, by referring in particular to the hypothesis of constant bending moments on each node, vector w collects the nodal displacements and vector m the constant bending moments on each node. Consequently, flexibility matrix F can be obtained by assembling the contributions of each node using Eq.(5) referred to the node influence area. Vector Dw represents the dual quantity to m and it can be obtained using Eq.(15) and Eq.(3). Once the stationarity condition is imposed compatibility equations and equilibrium equations are derived Fm + Dw = 0 (17) D T m + f = 0 (18) In order to solve the immediately above algebraic problem, the bending moment variables contained in vector m can be condensed by generating the system of equations D T FDw = f (19) which solved furnish the nodal displacement vector w and successively the bending moment vector m. The two elements delineated before, corresponding respectively to the hypothesis of moment constants on elements and on the influence area of each node, appear somewhat similar. However, the second way appears to be richer, since the bending moment on each side is obtained as an average between the nodal values and therefore represents a intermediate level between a constant and a linear function. Furthermore, since each node is independent of the other nodes and therefore the flexibility matrix F presents a diagonal block structure which, as will be made clear in the next Section, is useful in the solution of the elasto-plastic step of the incremental analysis. Finally, it has to be pointed out that the proposed mixed element respects the necessary condition for the stiffness matrix to have a sufficient rank. Indeed, this condition requires that the number of independent stress parameters be greater than or equal to the number of nodal displacement diminished by the number of rigid-body degrees of freedom. This condition is usually respected for both the finite elements presented above. 3 ELASTO-PLASTIC INCREMENTAL ANALYSIS The finite element analysis of elasto-plastic bodies subjected to monotonous loading can be performed following two ways. The first one derives essentially from the extension of the limit analysis which point to an accurate evaluation of the collapse load by completely ignoring the behaviour of the body before the collapse. The second approach follows an incremental iterative procedure which numerically simulates the elasto-plastic response to successive load increments. This approach appears 6

7 more attractive since it reconstructs the complete equilibrium path giving information not only on the limit point but also on the displacements and deformations reached. Furthermore, it is easily feasible for the most common computer. The incremental iterative approach requires that the loading path has to be discretized into a sequence of finite load increments. In each step, defined by the initial state and the load increment, a time interpolation law must be assumed, explicitely or implicitly, for the state variables: stresses and/or strains. This interpolation for the state variables is required since the elasto-plastic behaviour is path dependent. Naturally, various choices of the interpolation law produce different algorithms which give more or less accurate results. No interpolation law can be considered absolutely the best. Therefore, it could be useful to choose the interpolation law on the basis of computational convenience. A suggestive way, which is also convenient from the computational point of view, to define the loading path was suggested by Ponter and Martin 4, which proposed to follow the so-called extremal paths, which fulfil both the maximum strain work and the minimum complementary work in the loading step. For material which follows the elasto-perfectly plastic behaviour the solution on the extremal path is given by the Haar-Karman principle 5 which, for the plate bending problem, can be expressed as Φ[m m E ]= 1 (m m E ) T F(m m E )ds = min (20) 2 S where m represents the bending moment field in equilibrium with the external loads and which satisfy the yield condition, m E is the elastic solution which corresponds to the load increment and to the assigned initial state. By solving Equation (20) the elasto-plastic solution at the end of the incremental step can be obtained. Generally, the solution of this problem is very complicated, essentially because it is difficult to satisfy a priori all the equilibrium equations. However, this is not a complicated problem when finite elements are used. Indeed, in this case the Haar-Karman problem has to be solved in the particular case of nodal displacements assigned. This makes the calculation of m E straightforward by using the compatibility equations and the elastic constitutive law. Furthermore, being the displacement field assigned, no equilibrium equation has to be satisfied and the elasto-plastic bending moment m is bound to remain in the yield surface. These reasons lead, remembering the assumptions of moments constants on the influence area of each node, to solving the Haar-Karman problem for each node separately and making its solution simple. In other words, for each node the geometric problem of the contact between the Von Mises cylinder and the elastic complementary energy function has to be solved. Having defined the constitutive law in the incremental step it remains to satisfy the nonlinear equilibrium equations which can be written in general as p[λ] s[u[λ]] = 0 (21) 7

8 where p[λ] =p 0 + λˆp is the load vector described in a linear way by the initial load p 0, the λ parameter and a linear part ˆp. The term s[u[λ]] represents the structural reaction depending, clearly, on the nodal displacements vector and so on the λ parameter. The immediately above nonlinear equation, written for the generic step, can be tackled by means of a residual iterative scheme as r j = p s[u j ] (22) u j+1 = u j + K 1 r j (23) where K is a symmetric and positive definite matrix chosen in such a way that the iterative process converges. In order to precise the features of the iterative operator K the matrix K j is defined as s[u j+1 ] s[u j ]=K j (u j+1 u j ) (24) and substituting Eq.(24) in Eq.(23) the following equation is obtained r j+1 =(I K j K 1 )r j (25) So that the iterative process converges, i.e. r j 0, the spectral radius of the matrix I K K 1 j has to be less than unity. Further algebraic manipulation leads to the equivalent convergence condition 6 0 < K j < 2 K (26) Indicating with K E the elastic stiffness matrix and choosing the iterative matrix as K = K E 1 ω<2 (27) ω inequalities (26) are satisfied apart from the collapse zone where K j 0 and therefore a slow convergence has to be expected. This gap can be efficiently eliminated using the arc-length strategy proposed by Riks 7 for nonlinear elastic structures. Indeed, the loss of convergence of the numerical scheme derives from the parameter chosen to describe the equilibrium path, i.e. λ, sincethecurve is not analytical in this parameter. The main idea of the arc-length strategy consists in the choice of another parameter, the arc-length η, to describe the equilibrium curve: u = u[η] (28) λ = λ[η] (29) From this choice, the equilibrium equation (21) has to be completed with a constraint equation in order to be able to solve the problem. Different arc-length schemes are derived from different constraint equations. A simple and efficient way consists of choosing 8

9 the orthogonality between the extrapolation and the correction as constraint equation 6, deriving in this manner the following iteration scheme λ j+1 = λ j ω ut K 1 E r j (30) λ + ω u T û u j+1 = u j + ω ( K 1 E r j +(λ j+1 λj)û ) (31) where r j = p 0 + λ j ˆp s[u j ] (32) is the equilibrium error obtained by using λ j and u j and û = K 1 E ˆp. The quantities λ and u are the increments of the load parameter and displacements in the last step. The behaviour of the above described iterative process can be improved by choosing the parameter ω and the arc-length on the base of the previous steps as exposed in the Casciaro and Cascini paper 6. The mixed element exposed in Section 2 can be easily coded in the arc-length iteration scheme. Indeed, at the j-step the displacements vectors w j is assigned and the use of the compatibility equations provides the side rotations ϕ j = Dw j (33) and by means of the elastic constitutive law the elastic solution can be evaluated m Ej = m 0 + F 1 (ϕ j ϕ 0 ) (34) Successively, by using the Haar-Karman principle, the elasto-plastic solution m j can be found for each node. Finally, the structural reaction s can be derived from the equilibrium equations s[w j ]=D T m j (35) and the stiffness matrix of the mixed model can be used as iteration matrix K E = D T F 1 D (36) 4 NUMERICAL RESULTS In order to analyse the performance of the triangular mixed element presented in Section 2 a series of numerical tests, relative to simple form plates, was executed. Although, plates defined on complex domains could be analysed, these simple schemes are chosen since they allow a large comparison with results reported in technical literature. In Table 1 the normalized collapse load for a square plate varying the boundary conditions and the applied load are reported. In particular, the results for a simply supported and a clamped plate are reported and as load condition a distributed load q and a concentrated load P in the middle of the plate are considered. This Table reports the normalized 9

10 Figure 2: Mesh sequence used in the test problems collapse load (m u is the ultimate bending moment) obtained using meshes generated by an algorithm which subdivide each triangle in four triangles of the same form following the scheme reported in Figure 2. In particular, meshes having 41, 145 and 545 nodes are used. In the same Table are reported the results published by other authors using analytical calculation, lower and upper bounds, and numerical methods 8,9,10,11,12,13. The results reported in this Table suggest that the proposed finite element give quite accurate results. Figure 3 reports the normalized bending moments m x /m u and m y /m u on the middle line of a simply supported square plate in the collapse point. As can be seen from this Figure the curves obtained with the proposed model are very close to the curves reported in Casciaro and Di Carlo paper 12 that can be considered an accurate model as shown in Table 1. Finally, in Figure 4 are reported the equilibrium path and the evolution of the Von Mises function for a simply supported square plate loaded uniformly. The equilibrium curve is described by the load multiplier λ and the normalized middle point displacement w/w E where the index E refers to the elastic limit point. 5 CONCLUDING REMARKS The triangular mixed element proposed in this paper appears interesting when it is used in elasto-plastic analysis of plate bending problems. Indeed, elasto-plastic problems alway require discretization with many elements and therefore finite elements based on very simple algebra are advantageous if they are compared to more complicate finite elements. This advantage is raised when the analysis process requires a large number of operations on each element such as in the incremental iterative approach to the elastoplastic problems. The incremental approach to the elasto-plastic analysis becomes still more convenient when it is implemented by using the extremal path theory and the arclength strategy to reconstruct the loading path as far as the collapse load. 10

11 Clamped Simply supported Nodes ql 2 /m u P/m u ql 2 /m u Hodge and Belytschko Ranaweera and Leckie Lubliner Ang and Lopez Casciaro and Di Carlo Owen and Hinton exact Table 1: Collapse load for square plates The present work belongs to an initial stage of the investigation of the proposed model. Further extensions regard the study of plates defined on complex domains and the development of shell elements both in static and dynamic field. Acknowledgments The authors wish to thank Prof. Raffaele Casciaro for his suggestions and comments. 11

12 Figure 3: Bending moments for a square clamped plate on the A-B line 12

13 Figure 4: Equilibrium path and evolution of the Von Mises function 13

14 REFERENCES 1. R. Casciaro and L. Cascini, A Mixed Formulation and Mixed Finite Elements for Limit Analysis, Int. j. numer. methods eng., 18, , T. Kawai, New Discrete Structural Models and Generalization of the Method of Limit Analysis, Finite Elements in Nonlinear Mechanics, P. G. Bergan et al. eds, Tapir publishers, , J. Argyris, H. Balmer and I. S. Doltsinis, Some Thoughts on Shell Modelling for Crash Analysis, Computer methods in applied mechanics and engineering, 71, , A. R. S.Ponter and G. B. Martin, Some extremal properties and energy theorems for inelastic materials and their relationship to the deformation theory of plasticity, Int. j. mech. phys. solids, 20/5, , A. Haar and T. von Karman, Zur Theorie der Spannungzustande in Plastischen und Sandartingen Medien, Gottinger Nachrichten, math-phys. K R. Casciaro and L. Cascini, Limit analysis by incremental-iterative procedure, IU- TAM Conf. on Deformation and Failure of Granular Materials, Delft, E. Riks, An Incremental Approach to the Solution of Snapping and Buckling Problems, Int. j. solids struct., 15, , P. G. Jr. Hodge and T. Belyschko, Numerical Methods for the Limit Analysis of Plates, J. Applied Mechanics, , Dec M. P. Ranaweera and F. A. Leckie, Bound Methods in Limit Analysis, in Finite Element Techniques in Structural Mechanics, H. Tottenham and C. A. Brebbia eds. Proceedings of a seminar at the University of Southampton, Apr J. Lubliner, Plasticity Theory, Macmillan, New York A. H. Ang and L. A. Lopez, Discrete Model Analysis of Elasto-Plastic Plates, Proc. ASCE, Eng. Mech. Div., 94, , R. Casciaro e A. Di Carlo, Formulazione dell Analisi Limite delle Piastre come Problema di Minimax, Giornale del Genio Civile, 5, May D. R. J. Owen and E. Hinton, Finite Element in Plasticity: Theory and Practice, Pineridge Press Limited,

MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING

MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING 144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes

More information

CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES

CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES CHAPTER 14 BUCKLING ANALYSIS OF 1D AND 2D STRUCTURES 14.1 GENERAL REMARKS In structures where dominant loading is usually static, the most common cause of the collapse is a buckling failure. Buckling may

More information

Nonlinear bending analysis of laminated composite stiffened plates

Nonlinear bending analysis of laminated composite stiffened plates Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in

More information

The Finite Element Method for Solid and Structural Mechanics

The Finite Element Method for Solid and Structural Mechanics The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in

More information

Finite Element Method in Geotechnical Engineering

Finite Element Method in Geotechnical Engineering Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps

More information

International Journal of Advanced Engineering Technology E-ISSN

International Journal of Advanced Engineering Technology E-ISSN Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,

More information

Large Thermal Deflections of a Simple Supported Beam with Temperature-Dependent Physical Properties

Large Thermal Deflections of a Simple Supported Beam with Temperature-Dependent Physical Properties Large Thermal Deflections of a Simple Supported Beam with Temperature-Dependent Physical Properties DR. ŞEREF DOĞUŞCAN AKBAŞ Civil Engineer, Şehit Muhtar Mah. Öğüt Sok. No:2/37, 34435 Beyoğlu- Istanbul,

More information

Esben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer

Esben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics

More information

CRITERIA FOR SELECTION OF FEM MODELS.

CRITERIA FOR SELECTION OF FEM MODELS. CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.

More information

PILE SOIL INTERACTION MOMENT AREA METHOD

PILE SOIL INTERACTION MOMENT AREA METHOD Pile IGC Soil 2009, Interaction Moment Guntur, INDIA Area Method PILE SOIL INTERACTION MOMENT AREA METHOD D.M. Dewaikar Professor, Department of Civil Engineering, IIT Bombay, Mumbai 400 076, India. E-mail:

More information

Analytical formulation of Modified Upper Bound theorem

Analytical formulation of Modified Upper Bound theorem CHAPTER 3 Analytical formulation of Modified Upper Bound theorem 3.1 Introduction In the mathematical theory of elasticity, the principles of minimum potential energy and minimum complimentary energy are

More information

BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS

BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department

More information

GEOMETRIC NONLINEAR ANALYSIS

GEOMETRIC NONLINEAR ANALYSIS GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.

More information

VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS

VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.

More information

A *69>H>N6 #DJGC6A DG C<>C::G>C<,8>:C8:H /DA 'D 2:6G - ( - ) +"' ( + -"( (' (& -+" % '('%"' +"-2 ( -!"',- % )% -.C>K:GH>IN D; AF69>HH>6,-+

A *69>H>N6 #DJGC6A DG C<>C::G>C<,8>:C8:H /DA 'D 2:6G - ( - ) +' ( + -( (' (& -+ % '('%' +-2 ( -!',- % )% -.C>K:GH>IN D; AF69>HH>6,-+ The primary objective is to determine whether the structural efficiency of plates can be improved with variable thickness The large displacement analysis of steel plate with variable thickness at direction

More information

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009

More information

Discrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method

Discrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method 131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using

More information

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES

UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the

More information

Department of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil

Department of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE

More information

INFLUENCE OF FLANGE STIFFNESS ON DUCTILITY BEHAVIOUR OF PLATE GIRDER

INFLUENCE OF FLANGE STIFFNESS ON DUCTILITY BEHAVIOUR OF PLATE GIRDER International Journal of Civil Structural 6 Environmental And Infrastructure Engineering Research Vol.1, Issue.1 (2011) 1-15 TJPRC Pvt. Ltd.,. INFLUENCE OF FLANGE STIFFNESS ON DUCTILITY BEHAVIOUR OF PLATE

More information

Using MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,

Using MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup, Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal

More information

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi

More information

Transactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X

Transactions on Modelling and Simulation vol 9, 1995 WIT Press,  ISSN X An alternative boundary element formulation for plate bending analysis J.B. Paiva, L.O. Neto Structures Department, Sao Carlos Civil Engineering School, f, BrazzY Abstract This work presents an alternative

More information

Elastic plastic bending of stepped annular plates

Elastic plastic bending of stepped annular plates Applications of Mathematics and Computer Engineering Elastic plastic bending of stepped annular plates JAAN LELLEP University of Tartu Institute of Mathematics J. Liivi Street, 549 Tartu ESTONIA jaan.lellep@ut.ee

More information

INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS

INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS SHEN R. WU and LEI GU WILEY A JOHN WILEY & SONS, INC., PUBLICATION ! PREFACE xv PARTI FUNDAMENTALS 1 1 INTRODUCTION 3

More information

202 Index. failure, 26 field equation, 122 force, 1

202 Index. failure, 26 field equation, 122 force, 1 Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic

More information

Chapter 2 Finite Element Formulations

Chapter 2 Finite Element Formulations Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are

More information

Analysis of a portal steel frame subject to fire by use of a truss model

Analysis of a portal steel frame subject to fire by use of a truss model Analysis of a portal steel frame subject to fire by use of a truss model P. G. Papadopoulos & A. Mathiopoulou Department of Civil Engineering, Aristotle University of Thessaloniki, Greece Abstract A plane

More information

LINEAR AND NONLINEAR SHELL THEORY. Contents

LINEAR AND NONLINEAR SHELL THEORY. Contents LINEAR AND NONLINEAR SHELL THEORY Contents Strain-displacement relations for nonlinear shell theory Approximate strain-displacement relations: Linear theory Small strain theory Small strains & moderate

More information

Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method

Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method 9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached

More information

COUPLED FINITE-INFINITE ELEMENTS MODELING OF BUILDING FRAME-SOIL INTERACTION SYSTEM

COUPLED FINITE-INFINITE ELEMENTS MODELING OF BUILDING FRAME-SOIL INTERACTION SYSTEM VOL. 4, NO. 10, DECEMBER 009 SSN 1819-6608 006-009 Asian Research Publishing Network (ARPN. All rights reserved. COUPLED FNTE-NFNTE ELEMENTS MODELNG OF BULDNG FRAME-SOL NTERACTON SYSTEM Ramakant Agrawal

More information

Theoretical Manual Theoretical background to the Strand7 finite element analysis system

Theoretical Manual Theoretical background to the Strand7 finite element analysis system Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter

More information

Torsional Buckling of Double-Walled Carbon Nanotubes

Torsional Buckling of Double-Walled Carbon Nanotubes Torsional Buckling of Double-Walled Carbon Nanotubes S. H. Soong Engineering Science Programme, National University of Singapore Kent Ridge, Singapore 119260 ABSTRACT This paper is concerned with the torsional

More information

MODELLING MIXED-MODE RATE-DEPENDENT DELAMINATION IN LAYERED STRUCTURES USING GEOMETRICALLY NONLINEAR BEAM FINITE ELEMENTS

MODELLING MIXED-MODE RATE-DEPENDENT DELAMINATION IN LAYERED STRUCTURES USING GEOMETRICALLY NONLINEAR BEAM FINITE ELEMENTS PROCEEDINGS Proceedings of the 25 th UKACM Conference on Computational Mechanics 12-13 April 217, University of Birmingham Birmingham, United Kingdom MODELLING MIXED-MODE RATE-DEPENDENT DELAMINATION IN

More information

Back Calculation of Rock Mass Modulus using Finite Element Code (COMSOL)

Back Calculation of Rock Mass Modulus using Finite Element Code (COMSOL) Back Calculation of Rock Mass Modulus using Finite Element Code (COMSOL) Amirreza Ghasemi 1. Introduction Deformability is recognized as one of the most important parameters governing the behavior of rock

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how

More information

1 Static Plastic Behaviour of Beams

1 Static Plastic Behaviour of Beams 1 Static Plastic Behaviour of Beams 1.1 Introduction Many ductile materials which are used in engineering practice have a considerable reserve capacity beyond the initial yield condition. The uniaxial

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 11 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module - 01 Lecture - 11 Last class, what we did is, we looked at a method called superposition

More information

Plates and Shells: Theory and Computation. Dr. Mostafa Ranjbar

Plates and Shells: Theory and Computation. Dr. Mostafa Ranjbar Plates and Shells: Theory and Computation Dr. Mostafa Ranjbar Outline -1-! This part of the module consists of seven lectures and will focus on finite elements for beams, plates and shells. More specifically,

More information

POST-BUCKLING BEHAVIOUR OF IMPERFECT SLENDER WEB

POST-BUCKLING BEHAVIOUR OF IMPERFECT SLENDER WEB Engineering MECHANICS, Vol. 14, 007, No. 6, p. 43 49 43 POST-BUCKLING BEHAVIOUR OF IMPERFECT SLENDER WEB Martin Psotný, Ján Ravinger* The stability analysis of slender web loaded in compression is presented.

More information

ABHELSINKI UNIVERSITY OF TECHNOLOGY

ABHELSINKI UNIVERSITY OF TECHNOLOGY ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural

More information

First-Order Solutions for the Buckling Loads of Euler-Bernoulli Weakened Columns

First-Order Solutions for the Buckling Loads of Euler-Bernoulli Weakened Columns First-Order Solutions for the Buckling Loads of Euler-Bernoulli Weakened Columns J. A. Loya ; G. Vadillo 2 ; and J. Fernández-Sáez 3 Abstract: In this work, closed-form expressions for the buckling loads

More information

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS

A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,

More information

NETWORK FORMULATION OF STRUCTURAL ANALYSIS

NETWORK FORMULATION OF STRUCTURAL ANALYSIS Chapter 4 NETWORK FORMULATION OF STRUCTURAL ANALYSIS 4.1 INTRODUCTION Graph theoretical concepts have been widely employed for the analysis of networks in the field of electrical engineering. Kirchhoff

More information

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element Fourth International Conference on FRP Composites in Civil Engineering (CICE8) 22-24July 8, Zurich, Switzerland Dynamic and buckling analysis of FRP portal frames using a locking-free finite element F.

More information

INFLUENCE OF WEB THICKNESS REDUCTION IN THE SHEAR RESISTANCE OF NON-PRISMATIC TAPERED PLATE GIRDERS

INFLUENCE OF WEB THICKNESS REDUCTION IN THE SHEAR RESISTANCE OF NON-PRISMATIC TAPERED PLATE GIRDERS INFLUENCE OF WEB THICKNESS REDUCTION IN THE SHEAR RESISTANCE OF NON-PRISMATIC TAPERED PLATE GIRDERS Paulo J. S. Cruz 1, Lúcio Lourenço 1, Hélder Quintela 2 and Manuel F. Santos 2 1 Department of Civil

More information

Effect of Mass Matrix Formulation Schemes on Dynamics of Structures

Effect of Mass Matrix Formulation Schemes on Dynamics of Structures Effect of Mass Matrix Formulation Schemes on Dynamics of Structures Swapan Kumar Nandi Tata Consultancy Services GEDC, 185 LR, Chennai 600086, India Sudeep Bosu Tata Consultancy Services GEDC, 185 LR,

More information

1 Nonlinear deformation

1 Nonlinear deformation NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected

More information

Variational principles in mechanics

Variational principles in mechanics CHAPTER Variational principles in mechanics.1 Linear Elasticity n D Figure.1: A domain and its boundary = D [. Consider a domain Ω R 3 with its boundary = D [ of normal n (see Figure.1). The problem of

More information

IMPROVING LATERAL STIFFNESS ESTIMATION IN THE DIAGONAL STRUT MODEL OF INFILLED FRAMES

IMPROVING LATERAL STIFFNESS ESTIMATION IN THE DIAGONAL STRUT MODEL OF INFILLED FRAMES IMPROVING LATERAL STIFFNESS ESTIMATION IN THE DIAGONAL STRUT MODEL OF INFILLED FRAMES I.N. Doudoumis 1 1 Professor, Dept. of Civil Engineering, Aristotle University of Thessaloniki, Thessaloniki, Greece

More information

Software Verification

Software Verification EXAMPLE 1-026 FRAME MOMENT AND SHEAR HINGES EXAMPLE DESCRIPTION This example uses a horizontal cantilever beam to test the moment and shear hinges in a static nonlinear analysis. The cantilever beam has

More information

Chapter 3 Variational Formulation & the Galerkin Method

Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation

More information

Numerical Solution of Non-Linear Biharmonic. Equation for Elasto-Plastic Bending Plate

Numerical Solution of Non-Linear Biharmonic. Equation for Elasto-Plastic Bending Plate Applied Mathematical Sciences, Vol. 9, 5, no. 6, 769-78 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ams.5.575 Numerical Solution of Non-Linear Biharmonic Equation for Elasto-Plastic Bending Plate

More information

Bending of Simply Supported Isotropic and Composite Laminate Plates

Bending of Simply Supported Isotropic and Composite Laminate Plates Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,

More information

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation

More information

Stress analysis of a stepped bar

Stress analysis of a stepped bar Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.

More information

FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES

FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES th International Conference on Earthquake Geotechnical Engineering June 5-8, 7 Paper No. 11 FINITE GRID SOLUTION FOR NON-RECTANGULAR PLATES A.Halim KARAŞĐN 1, Polat GÜLKAN ABSTRACT Plates on elastic foundations

More information

LINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL

LINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL LINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL SREELATHA P.R * M.Tech. Student, Computer Aided Structural Engineering, M A College of Engineering, Kothamangalam 686 666,

More information

DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS

DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the

More information

Sborník vědeckých prací Vysoké školy báňské - Technické univerzity Ostrava číslo 1, rok 2014, ročník XIV, řada stavební článek č. 18.

Sborník vědeckých prací Vysoké školy báňské - Technické univerzity Ostrava číslo 1, rok 2014, ročník XIV, řada stavební článek č. 18. Sborník vědeckých prací Vysoké školy báňské - Technické univerzity Ostrava číslo, rok, ročník XIV, řada stavební článek č. 8 Martin PSOTNÝ NONLINER NLYSIS OF BUCKLING & POSTBUCKLING bstract The stability

More information

Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling

Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling F. Portioli, L. Cascini, R. Landolfo University of Naples Federico II, Italy P. Foraboschi IUAV University,

More information

Chapter 2. The constitutive relation error method. for linear problems

Chapter 2. The constitutive relation error method. for linear problems 2. The constitutive relation error method for linear problems Chapter 2 The constitutive relation error method for linear problems 2.1 Introduction The reference problem considered here is the linear problem

More information

PLAXIS. Scientific Manual

PLAXIS. Scientific Manual PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3

More information

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems

The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 1-21 September, 2017 Institute of Structural Engineering

More information

Finite Element Method

Finite Element Method Finite Element Method Finite Element Method (ENGC 6321) Syllabus Objectives Understand the basic theory of the FEM Know the behaviour and usage of each type of elements covered in this course one dimensional

More information

Structural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed

Structural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed Structural Analysis of Truss Structures using Stiffness Matrix Dr. Nasrellah Hassan Ahmed FUNDAMENTAL RELATIONSHIPS FOR STRUCTURAL ANALYSIS In general, there are three types of relationships: Equilibrium

More information

On Nonlinear Buckling and Collapse Analysis using Riks Method

On Nonlinear Buckling and Collapse Analysis using Riks Method Visit the SIMULIA Resource Center for more customer examples. On Nonlinear Buckling and Collapse Analysis using Riks Method Mingxin Zhao, Ph.D. UOP, A Honeywell Company, 50 East Algonquin Road, Des Plaines,

More information

Methods of Analysis. Force or Flexibility Method

Methods of Analysis. Force or Flexibility Method INTRODUCTION: The structural analysis is a mathematical process by which the response of a structure to specified loads is determined. This response is measured by determining the internal forces or stresses

More information

JEPPIAAR ENGINEERING COLLEGE

JEPPIAAR ENGINEERING COLLEGE JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III

More information

Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure

Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure ISSN : 48-96, Vol. 6, Issue 8, ( Part -4 August 06, pp.3-38 RESEARCH ARTICLE Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure Dr.Ragbe.M.Abdusslam Eng. Khaled.S.Bagar ABSTRACT

More information

Steps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen

Steps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,

More information

Ultimate shear strength of FPSO stiffened panels after supply vessel collision

Ultimate shear strength of FPSO stiffened panels after supply vessel collision Ultimate shear strength of FPSO stiffened panels after supply vessel collision Nicolau Antonio dos Santos Rizzo PETROBRAS Rio de Janeiro Brazil Marcelo Caire SINTEF do Brasil Rio de Janeiro Brazil Carlos

More information

NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS

NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.

More information

NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS

NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS IGC 009, Guntur, INDIA NUMERICAL ANALYSIS OF A PILE SUBJECTED TO LATERAL LOADS Mohammed Younus Ahmed Graduate Student, Earthquake Engineering Research Center, IIIT Hyderabad, Gachibowli, Hyderabad 3, India.

More information

Institute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I

Institute of Structural Engineering Page 1. Method of Finite Elements I. Chapter 2. The Direct Stiffness Method. Method of Finite Elements I Institute of Structural Engineering Page 1 Chapter 2 The Direct Stiffness Method Institute of Structural Engineering Page 2 Direct Stiffness Method (DSM) Computational method for structural analysis Matrix

More information

Schur decomposition in the scaled boundary finite element method in elastostatics

Schur decomposition in the scaled boundary finite element method in elastostatics IOP Conference Series: Materials Science and Engineering Schur decomposition in the scaled boundary finite element method in elastostatics o cite this article: M Li et al 010 IOP Conf. Ser.: Mater. Sci.

More information

EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A SPHERE AGAINST A RIGID FLAT A FINITE ELEMENT STUDY

EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A SPHERE AGAINST A RIGID FLAT A FINITE ELEMENT STUDY Proceedings of the International Conference on Mechanical Engineering 2009 (ICME2009) 26-28 December 2009, Dhaka, Bangladesh ICME09- EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A

More information

Constitutive models: Incremental plasticity Drücker s postulate

Constitutive models: Incremental plasticity Drücker s postulate Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager

More information

MITOCW MITRES2_002S10linear_lec07_300k-mp4

MITOCW MITRES2_002S10linear_lec07_300k-mp4 MITOCW MITRES2_002S10linear_lec07_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources

More information

Plate analysis using classical or Reissner- Mindlin theories

Plate analysis using classical or Reissner- Mindlin theories Plate analysis using classical or Reissner- Mindlin theories L. Palermo Jr. Faculty of Civil Engineering, State Universiv at Campinas, Brazil Abstract Plates can be solved with classical or Reissner-Mindlin

More information

PREDICTION OF COLLAPSED LOAD OF STEEL COLUMNS USING FINITE STRIP METHOD

PREDICTION OF COLLAPSED LOAD OF STEEL COLUMNS USING FINITE STRIP METHOD International Journal of Civil Engineering and Technology (IJCIET) Volume 9, Issue 8, August 2018, pp. 347 357, Article ID: IJCIET_09_08_035 Available online at http://www.iaeme.com/ijciet/issues.asp?jtype=ijciet&vtype=9&itype=8

More information

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These

More information

Determining Slope Stability by Rigid Finite Element Method and Limit Analysis Methods

Determining Slope Stability by Rigid Finite Element Method and Limit Analysis Methods Journal of Novel Applied Sciences Available online at www.jnasci.org 2015 JNAS Journal-2015-4-5/579-589 ISSN 2322-5149 2015 JNAS Determining Slope Stability by Rigid Finite Element Method and Limit Analysis

More information

On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity

On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity 63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain

More information

NONLOCAL PLASTICITY APPLIED TO BENDING OF BEAMS

NONLOCAL PLASTICITY APPLIED TO BENDING OF BEAMS IX International Conference on Computational Plasticity COMPLAS IX E. Oñate and D. R. J. Owen (Eds) CIMNE, Barcelona, 2007 NONLOCAL PLASTICITY APPLIED TO BENDING OF BEAMS L. Strömberg Div of Structural

More information

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e

More information

FETI domain decomposition method to solution of contact problems with large displacements

FETI domain decomposition method to solution of contact problems with large displacements FETI domain decomposition method to solution of contact problems with large displacements Vít Vondrák 1, Zdeněk Dostál 1, Jiří Dobiáš 2, and Svatopluk Pták 2 1 Dept. of Appl. Math., Technical University

More information

DUCTILITY BEHAVIOR OF A STEEL PLATE SHEAR WALL BY EXPLICIT DYNAMIC ANALYZING

DUCTILITY BEHAVIOR OF A STEEL PLATE SHEAR WALL BY EXPLICIT DYNAMIC ANALYZING The 4 th World Conference on arthquake ngineering October -7, 008, Beijing, China ABSTRACT : DCTILITY BHAVIOR OF A STL PLAT SHAR WALL BY XPLICIT DYNAMIC ANALYZING P. Memarzadeh Faculty of Civil ngineering,

More information

THE USE OF DYNAMIC RELAXATION TO SOLVE THE DIFFERENTIAL EQUATION DESCRIBING THE SHAPE OF THE TALLEST POSSIBLE BUILDING

THE USE OF DYNAMIC RELAXATION TO SOLVE THE DIFFERENTIAL EQUATION DESCRIBING THE SHAPE OF THE TALLEST POSSIBLE BUILDING VII International Conference on Textile Composites and Inflatable Structures STRUCTURAL MEMBRANES 2015 E. Oñate, K.-U.Bletzinger and B. Kröplin (Eds) THE USE OF DYNAMIC RELAXATION TO SOLVE THE DIFFERENTIAL

More information

PLANE STRAIN SHAKEDOWN ANALYSIS

PLANE STRAIN SHAKEDOWN ANALYSIS PLANE STRAIN SHAKEDOWN ANALYSIS BAHER L. ABOUSTAIT Head, Construction Engineering Misr University For Science And Technology Six of October City, Egypt baboustait@must.edu.eg, baher_aboustait@hotmail.com

More information

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation

More information

Large deflection analysis of planar solids based on the Finite Particle Method

Large deflection analysis of planar solids based on the Finite Particle Method yuying@uiuc.edu 10 th US National Congress on Computational Mechanics Large deflection analysis of planar solids based on the Finite Particle Method 1, 2 Presenter: Ying Yu Advisors: Prof. Glaucio H. Paulino

More information

Ph.D. Preliminary Examination Analysis

Ph.D. Preliminary Examination Analysis UNIVERSITY OF CALIFORNIA, BERKELEY Spring Semester 2014 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials Name:......................................... Ph.D.

More information

M.S Comprehensive Examination Analysis

M.S Comprehensive Examination Analysis UNIVERSITY OF CALIFORNIA, BERKELEY Spring Semester 2014 Dept. of Civil and Environmental Engineering Structural Engineering, Mechanics and Materials Name:......................................... M.S Comprehensive

More information

Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams.

Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Outline of Continuous Systems. Introduction to Continuous Systems. Continuous Systems. Strings, Torsional Rods and Beams. Vibrations of Flexible Strings. Torsional Vibration of Rods. Bernoulli-Euler Beams.

More information

Kirchhoff Plates: Field Equations

Kirchhoff Plates: Field Equations 20 Kirchhoff Plates: Field Equations AFEM Ch 20 Slide 1 Plate Structures A plate is a three dimensional bod characterized b Thinness: one of the plate dimensions, the thickness, is much smaller than the

More information

Structural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).

Structural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free

More information

. D CR Nomenclature D 1

. D CR Nomenclature D 1 . D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the

More information

Chapter 5 Structural Elements: The truss & beam elements

Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations

More information