UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES

Size: px
Start display at page:

Download "UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES"

Transcription

1 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 requirements for the degree of MASTER OF SCIENCE August 2008 Major Subject: Mechanical Engineering

2 ii 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 requirements for the degree of MASTER OF SCIENCE Approved by: Chair of Committee Committee Members Head of Department J. N. Reddy Xinlin Gao Zoran Sunik Dennis O Neal August 2008 Major Subject: Mechanical Engineering

3 iii ABSTRACT Unconventional Finite Element Models for Nonlinear Analysis of Beams and Plates. (August 2008) Wooram Kim B.S. Korea Military Academy Chair of Advisory Committee: Dr. J. N. Reddy In this thesis mixed finite element models of beams and plates bending are developed to include other variables (i.e. the membrane forces and shear forces) in addition to the bending moments and vertical deflection and to see the effect of it on the nonlinear analysis. Models were developed based on the weighted residual method. The effect of inclusion of additional variables is compared with other mixed models to show the advantage of the one type of model over other models. For beam problems the Euler-Bernoulli beam theory and the Timoshenko beam theory are used. And for the plate problems the classical plate theory and the first-order shear deformation plate theory are used. Each newly developed model is examined and compared with other models to verify its performance under various boundary conditions. In the linear convergence study solutions are compared with analytical solutions available and solutions of existing models. For non-linear equation solving direct method and Newton-Raphson method are used to find non-liner solutions. Then converged solutions are compared with available solutions of the displacement models. Noticeable improvement in accuracy of force-like variables (i.e. shear resultant membrane resultant and bending moments) at the boundary of elements can be achieved by using present mixed models in both linear and nonlinear analysis. Post processed data of newly developed mixed models show better accuracy than existing displacement based and mixed models in both of vertical displacement and force-like variables. Also present beam and plate finite element models allow use of relatively lower level of interpolation function without causing severe locking problems.

4 iv DEDICATION To my God my mother grandfather grandmother and all of my family members who have always been there for me.

5 v ACKNOWLEDGMENTS I am greatly indebted to my advisor Dr. J.N.Reddy for his valuable comments and advice. Without his help this thesis would not have been completed. Sincere thanks to Dr. Xin-Lin Gao and Dr. Zoran Sunik for their valuable recommendations and consideration as committee members.

6 vi TABLE OF CONTENTS Page ABSTRACT... iii DEDICATION iv ACKNOWLEDGMENTS... v TABLE OF CONTENTS. vi LIST OF FIGURES. viii LIST OF TABLES x CHAPTER I INRODUCTION Review of Euler-Bernoulli Beam Theory Kinematics of EBT Equilibriums of EBT Constitutive Relations and Resultants of EBT Review of Timoshenko Beam Theory (TBT) Kinematics of TBT Equilibriums of TBT Constitutive Relations and Resultants of TBT Review of Classical Plate Theory (CPT) Kinematics of CPT Equilibriums of CPT Constitutive Relations and Resultants of CPT Review of First Order Shear Deformation Theory (FSDT) Kinematics of FSDT Equilibriums of FSDT Constitutive Relations and Resultants of FSDT CHAPTER II DEVELOPMENT OF BEAM BENDING MODELS Model I of Beam Bending Weighted Residual Statements of Model I Finite Element Equations of Model I Model II of Beam Bending Weighted Residual Statements of Model II Finite Element Equations of Model II Model III of Beam Bending Weighted Residual Statements of Model III Finite Element Equations of Model III Model IV of Beam Bending Weighted Residual Statements of the Model IV Finite Element Equations of Model IV Lagrange Type Beam Finite Elements

7 vii Page CHAPTER III DEVELOPMENT OF PLATE BENDING MODELS Model I of Plate Bending Weighted Residual Statements of Model I Finite Element Equations of Model I Model II of Plate Bending Weighted Residual Statements of Model II Finite Element Equations of Model II Model III of Plate Bending Weighted Residual Statements of Model III Finite Element Equations of Model III Model IV of Plate Bending Weighted Residual Statements of Model IV Finite Element Equations of Model IV The Lagrange Type Plate Finite Elements CHAPTER IV NONLINEAR EQUATION SOLVING PROCEDURES Direct Iterative Method Algorithm of Direct Iterative Method Newton-Raphson Iterative Method Algorithm of Newton-Raphson Iterative Method Calculation of Tangent Stiffness Matrices Tangent Stiffness Matrices Load Increment Vector CHAPTER V NUMERICAL RESULTS Numerical Analysis of Nonlinear Beam Bending Description of Problem[1] Numerical Results Numerical Analysis of Nonlinear Plate Bending Description of Problem Non-dimensional Analysis of Linear Solutions Non-linear Analysis CHAPTER VI CONCLUSION REFERENCES VITA. 115

8 viii LIST OF FIGURES Page Fig. 1.1 Undeformed and deformed EBT and TBT beams source from[2] 3 Fig. 1.2 A typical beam element with forces and bending moments... 5 Fig. 1.3 Undeformed and deformed CPT and FSDT plates source from[2] 9 Fig. 1.4 A typical 2-D plate element[3] with forces and moments Fig. 2.1 Fig. 3.1 Node number and local coordinate of the line elements of the Lagrange family Node number and local coordinate of the rectangular elements of the Lagrange family 68 Fig. 4.1 A flow chart[1] of the direct iteration method 70 Fig. 4.2 A flow chart[1] of the Newton iteration method 72 Fig. 5.1 Description of the beam geometry 82 Fig. 5.2 Symmetry boundary conditions of beams. 83 Fig. 5.3 A comparison of the non-linear solutions of beams.. 85 Fig. 5.4 A comparison of the membrane locking in various models.. 89 Fig. 5.5 A comparison of effect of the length-to-thickness ratio on the beam 91 Fig. 5.6 A description of the plate bending problem 93 Fig. 5.7 Fig. 5.8 Fig. 5.9 Fig Symmetry boundary conditions[1 25] of a quadrant of the square plate. 94 Plots of the membrane and normal stress of Model I II and CPT displacement model under SS3 boundary condition Plots of the center deflection normal and membrane stress of Model III with that of the FSDT displacement model under SS1 and SS3 boundary conditions. 106 Post processed quadrant images of the variables in various models SS3 with converged solution at load parameter

9 ix Fig Fig Page Plots of the non-linear membrane stresses of Model III and FSDT displacement model along the x = Plots of the non-linear bending moments of Model III and FSDT displacement model along the x =

10 x LIST OF TABLES Page Table 5.1 Comparison of mixed models and displacement based models.. 84 Table 5.2 Membrane locking in mixed models and the displacement models 87 Table 5.3 Effect of the membrane locking in the mode I and II. 88 Table 5.4 Table 5.5 Effect of the length-to-thickness ratio on the deflections in TBT beam 90 Comparison of the effect of the length-to-thickness ratio in the EBT and the TBT beams 90 Table 5.6 Comparison of Model III with other mixed models.. 92 Table 5.7 Table 5.8 Table 5.9 Table 5.10 Table 5.11 Table 5.12 Comparison of the linear solution of various CPT Models isotropic ( ν 0.3 ) square plate simple supported (SS1) Comparison of the linear solution of various CPT Models isotropic ( ν 0.3 ) square plate clamped (CC) Comparison of the CPT linear solution with that of the displacement model isotropic ( ν 0.25 ) square plate simple supported (SS1). 98 Comparison of the CPT linear solution with that of the displacement Model isotropic ( ν 0.25 ) square plate clamped(cc) Comparison of the current mixed FSDT linear solution with that of the other mixed model (Reddy [7] ) with isotropic(ν 0.25 K 5/6 ) square plate simple supported(ss1). Comparison of the linear solution of the FSDT with isotropic ( ν 0.25 K 5/6 ) square plate simple supported (SS1)... Table 5.13 Effect of reduced integration in Model I and II Table 5.14 Table 5.15 Table Comparison of the center deflection and normal stress of Model I and II with the CPT displacement 102 model Comparison of the convergence of Model I II III and IV under the SS1 and SS3 boundary conditions 104 Comparison of the center deflection and normal stress of Model III with the FSDT displacement model under SS1 and SS3 boundary conditions.. 105

11 1 CHAPTER I INTRODUCTION The objective of this study is to investigate the performance of finite element models based on mixed weighted-residual formulations of beams and plates. In particular the study investigates merits and demerits of the newly developed mixed finite element models of beam and plate bending based on weighted-residual and mixed formulations. The von Karman nonlinear equations[1 2] of beams and plates[1 3] are used to develop alternative finite element models to the conventional displacementbased finite element models[4]. Once the basic models are developed and critically evaluated in comparison to the conventional displacement-based finite element models they can be extended to other beam and plate structures with proper modifications. For example the plate bending models can be extended to the laminated composite structures with proper laminate equations[3 5]. The mixed finite element models of beams and plates were developed more than two decades ago by Putcha and Reddy[6 7] to overcome the drawbacks of the displacement based models. The basic idea of mixed finite element model is to include more than two different types of fields in the finite element model as independent variables. For example the bending moment of the beam element can be included as independent variable in addition to the axial and transverse displacements. The mixed finite element models[7] developed in past only included bending moments as independent variables to reduce the differentiability of the transverse displacement component[7]. This mixed models can provide the same level of accuracy for the bending moment as that for the displacement fields whereas in the displacement based model the bending moment is calculated at points other than nodes in the post processing step[1 4]. Thus the displacement finite element models cannot provide the same level of accuracy for force-like variables as in the mixed finite element models. This thesis follows the style of Computational Method in Applied Mechanics and Engineering.

12 2 In the present study mixed finite element models are developed to include other variables (i.e. the membrane forces and shear forces) in addition to the bending moments and to see the effect of them on the nonlinear analysis. The effect of including other variables will be compared with different mixed models to show the advantage of the one type of model over other models. For the nonlinear beam bending problems[1 8] three different mixed models based on the Euler-Bernoulli beam theory[4 9] and one mixed model based on the Timoshenko beam theory[4 9] are developed. For the nonlinear plate bending problems two different mixed models based on the classical plate theory[1 3] and two mixed models based on the first-order shear deformation plate theory [1 3] are developed. To verify the performance of the newly developed finite element models numerical results of them are compared with those of the existing displacement based finite element models[1 7]. For each beam bending model three types of boundary conditions (i.e. clamped-clamped (CC) hinged-hinged (HH) and pined-pined (PP) boundary conditions.) are examined and the results are compared. For plate bending model three types of boundary conditions (i.e. the simple support I (SS1) the simple support III (SS3) and the clamped (CC) boundary conditions.) are examined and the results are compared with those of the conventional finite element models[1 7]. For each of the beam models three different Lagrange type interpolation functions[4 10] (i.e. linear quadratic and cubic) are used for the approximation of the variables to see the relations between the degree of interpolation functions and the accuracy of the solutions. For each of the plate bending models two different Lagrange type interpolation functions[4 10] (i.e. 4-nodes and 9 nodes) are considered. Then the post-processed data on the stresses and the moments of the equilibrium state in the various models are compared. The finite element Models are implemented using Maple 9.5 [11] for the beam bending models and the Fortran [12] for the plates bending problems. All graphs and data are obtained by using the MS Excel and the Matlab 7.1[13].

13 3 1.1 Review of Euler-Bernoulli Beam Theory Fig. 1.1 Undeformed and deformed EBT and TBT beams source from [2]. To develop new nonlinear mixed finite element models of beam bending the Euler-Bernoulli beam theory and the Timoshenko beam theory are considered. Due to the assumption of moderate rotation[1] of a beam cross section perpendicular to the x- axis a geometric nonlinearity[1] can be considered for the present study. As a consequence the nonlinearity only appears as square of the slope (i.e. / ) in the formulations. Detailed geometry and characteristics of both beam bending theories can be found in the Fig. 1.1.

14 Kinematics of EBT By taking the horizontal axis (i.e. longitudinal direction of the beam) of the beam to be located along the x-axis and the vertical axis (i.e. direction along the height) to be located along the z-axis the displacement field[3 8 14] of the EBT can be given as follow (see Reddy[1]): 0. (1.1) The von Karman strain[1 2] associated with the displacement field of the EBT is given as follow:. (1.2) where and are defined as 1 2. (1.3) Equilibriums of EBT Here the equilibrium equations[11] of the EBT are derived by using the force and the moment equilibrium of the infinitesimal free body diagram[1] given in the Fig The vertical shear force resultant can be defined only in terms of the bending moment and certain portion of the membrane force with the nonlinear assumption.

15 5 Fig. 1.2 A typical infinitesimal beam element with forces and bending moments. By Taylor s expansion[15 16] each of the resultants on the right hand side of the free body diagram (see Fig. 1.2) can be expanded to the left hand side by following equation [15 17] 1 2! 1 1! 1! (1.4) where is an arbitrary resultant on the left side of the element and is the associated value on the right side of the element. Then every term multiplied by from the expansion can be omitted by taking the limit of 0. Then the x- direction and the z-direction force equilibrium can be obtained as follow:

16 (1.5) By taking the positive y-direction to be the direction of going through the board the y-direction moment equilibrium can be written as follow: (1.6) We can obtain a point equilibrium of the forces and the moment by dividing above equations by and taking the limit of 0. Finally following equilibrium equations of EBT[1] can be obtained (1.7) Constitutive Relations and Resultants of EBT In the present study the beam is assumed to obey the linear elastic relation thus the stress of the EBT can be related to the strain by the relation known as Hooke s law as follow:

17 7. (1.8) The stress and moment resultants[1] can be defined as 1 2 (1.9) where the is the cross section area of the beam the [1] is the second moment inertia of the cross section(about the y-axis) and the is the Young s modulus[18 19] or elastic modulus. 1.2 Review of Timoshenko Beam Theory (TBT) Kinematics of TBT By taking the horizontal axis (i.e. longitudinal direction) of the beam to be located along the x-axis and the vertical axis (i.e. direction along the height) along the z-axis the displacement field [1 9 14] of the TBT can be given as follow: 0. (1.10) Note that instead of the slope of the deformed beam axis (i.e. / ) the shear rotation[14] was included to account for the shear rotation of the cross section.

18 8 The von Karman strain[1] associated with the displacement field of TBT can be given as follow: (1.11) where ε ε and ε are defined as (1.12) Equilibriums of TBT By substituting the strains into the virtual work statement[8] the equilibrium equations of the TBT can be obtained. By the principle of the virtual work[1 8] it can be stated if a body is in equilibrium the total virtual work done by actual internal as well as external forces in moving through their respective virtual displacement is zero [1]. It can be expressed by following equation[8] where 0 (1.13)

19 9 Note that are the generalized nodal forces[1] and are the virtual generalized nodal displacements[1]. And the virtual strains and the definitions of the resultant forces can be defined as follow: 2. (1.14) By substituting the force resultants the moment resultants and the virtual strains given in the (1.14) into the (1.13) the following energy equation can be directly obtained (1.15) Then by collecting the coefficients of the variations of the displacement terms in the (1.15) following equilibrium equations[1 8] of the TBT can be obtained(for details see Reddy[1]).

20 (1.16) By comparing the equilibrium equations of the EBT given in the (1.7) and the TBT given in the (1.16) it can be shown that the shear resultant of the EBT can be related to the shear resultant of the TBT by the following equation[1].. (1.17) Essentially the equilibrium equations of the EBT and the TBT are the same but the specific variables involved may have different meanings. In this case is the shear resultant acting on the plane perpendicular to the x axis while is the shear resultant acting on the deformed plane (see Reddy[1]) Constitutive Relations and Resultants of TBT Since there are two non-zero strain components in the TBT we have two stress components from the constitutive relations. By assuming that the beam obeys linear elastic relation the stresses can be related to strains as follow: 2. (1.18) The generalized resultant forces can be calculated by the definition given in the (1.14) as follow:

21 (1.19) where is the cross section area of the beam is the second moment inertia of the cross section 5/6 is the shear correction factor[14] is the shear modulus[19] and is the Young s modulus Review of Classical Plate Theory (CPT) Fig. 1.3 Undeformed and deformed CPT and FSDT plates source from [2].

22 12 The major difference between the CPT and FSDT comes from the displacement field given in the Fig Kinematics of CPT The CPT can be considered as an extended 2-D version of the EBT. Thus the displacement field of the CPT is very similar to that of the EBT. The displacement field of the CPT with Kirchhoff hypoth esis[1 3] can be given by. (1.20) And with the assumption[1] of small strain but moderately large rotation we can simplify the components of the nonlinear strain tensor[20]. Then the components of the strain tensor is given by (1.21) By substituting the displacement field described in the (1.20) into the components of the strain tensor given in the (1.21) the following specific von Karman nonlinear strains[1] of the CPT can be obtained.

23 (1.22) Equilibriums of CPT The CPT can be derived by using vector approach[11] with the infinitesimal free diagram of the Fig Since it is assumed that the plane stress condition is still valid for the in-plane forces the x y and z-direction force equilibriums and x and y-direction moment equilibriums of the infinitesimal plate element[3] can be stated by using the generalized force resultants as described in the (1.23a to e). Fig. 1.4 A typical 2-D plate element[3] with forces and moments.

24 14 With the Taylor s expansion given in the (1.4) we can set the equilibriums of the forces and the moments in the given directions as follow: The x-direction force equilibrium: 0. (1.23a) The y-direction force equilibrium: 0. (1.23b) The z-direction force equilibrium: 0. (1.23c) In the force equilibriums of the CPT it is assumed that the plate is in the plane stress condition[18] for the in-plane forces. So by including some portion of in-plane

25 15 forces only for moment equilibriums following equations of moment equilibrium can be obtained. The y-direction moment equilibrium: 2 0. (1.23d) 2 The x-direction moment equilibrium: 2 0. (1.23e) 2

26 Constitutive Relations and Resultants of CPT In the present plate bending problem it is assumed that every in-plane stress and strain remain plane stress[18] condition. For the orthotropic plane stress condition the constitutive relations of stresses and strains can be given in the matrix form equation as follow: 0 0 (1.24) where the components of the matrix are given by (1.25) Note and are the elastic modulus [19] of the x and y-direction respectively and are the Poisson s ratio and the is the shear modulus. By dividing the (1.23a to e) by and taking the limit of 0 the following equilibrium equations of the CPT can be obtained (1.26)

27 17 By locating the x axis along the mid-plane of the plate element the resultants of the CPT can be defined by the following equations.. (1.27) By using the constitutive relations given in the (1.25) and the definitions of the resultants given in the (1.26) the following equations of the resultants can be obtained (1.28) The matrix and can be defined by / 1. (1.29) / where

28 Review of First Order Shear Deformation Theory (FSDT) The FSDT can be considered as 2-D version of the TBT as the CPT can be said to be the 2-D version of the EBT. In the FSDT the same nonlinearity used in the CPT is assumed. Thus the components of the strain tensor given in the (1.21) can be used to obtain the specific strains of the FSDT Kinematics of FSDT The displacement field[3] of the FSDT with Kirchhoff hypothesis can be given by (see Reddy[1 3] for details). (1.30) By substituting the displacement field of the FSDT given in the (1.30) into components of the strain tensor given in the (1.21) strains[1] of the FSDT can be given by (1.31)

29 Equilibriums of FSDT The equilibrium equations of the FSDT can be derived using the virtual work statement (see Reddy[1] for details) with strains given in the (1.31). The equilibrium equations[3] of the FSDT can be given by (1.32) In the part 1.1 the shear resultant of the EBT was related to the shear resultant of the TBT by the relation given in the (1.17). In similar sense the shear resultants of the CPT can be related to that of the FSDT by the following equations by comparing the equilibrium equations of the CPT and that of the FSDT.. (1.33) Constitutive Relations and Resultants of FSDT The FSDT has two more non-zero strains compared with the strains of the CPT. Additional strains of the FSDT can be related to the corresponding stress components by (1.34)

30 20 where and are the shear modulus and respectively. In addition to the in-plane stresses and strains relations given in the (1.24) we have the following additional relations between the shear stresses and the shear resultants of the FSDT which can be defined by. (1.35) By substituting the constitutive relations given in the (1.24) and (1.34) into the definitions of the resultant forces given in the (1.27) and (1.35) the following equations of the resultants[3] of the FSDT can be obtained (1.36) 2 where 5/6 is the shear correction factor[1] and 1 (see the (1.25) and (1.34) for the specific values of the ). 2

31 21 CHAPTER II DEVELOPMENT OF BEAM BENDING MODELS In this chapter development of various types of the nonlinear mixed finite element models of the beam bending problem is discussed. In current models force like physical variables are included as independent nodal variables with proper weighted residual statements[4]. Four different nonlinear mixed finite element models of beam bending are developed for the numerical analysis. The relation between the participation of a typical variable and the accuracy of the linear and the nonlinear solutions are investigated in the chapter V. To clarify the developing procedure the governing equations of the EBT and the TBT were brought from the chapter I. - Governing equations of the EBT (1.7) 1 2. (1.9) - Governing equations of the TBT (1.16)

32 (1.19) 2.1 Model I of Beam Bending Weighted Residual Statements of Model I The governing equations of the EBT which were derived in the chapter I are used to develop the Model I of beam bending. The displacements (i.e. and ) and generalized forces (i.e. and ) are included as independent variables in the beam bending Model I. By using the equilibrium equation and the resultant equations of the EBT following weighted residual statements can be made (2.1) where and.

33 23 Note that variables with superscripted a (i.e. and ) denote approximated variables 1 5 denotes the i th weight function of the i th weighted residual statement. and are the global coordinates of element region. The boundary terms in the first the second and the third equations can be obtained by conducting the integration by parts[4] of the related terms Finite Element Equations of Model I Next with the weighted residual statements given in the (2.1) the variables can be approximated with the proper interpolation functions. Compared with the EBT displacement based model[1] whose variable (i.e. vertical displacement) should be approximated with the Hermite interpolation functions[1 4] the model I allows the use of the Lagrange interpolation functions for the approximation of all variables of it because weighted residual statements do not include any derivative of variable as primary variable.. (2.2) By observing the boundary terms in the (2.1) which are produced by integration by parts with the chosen weight functions in the (2.2) the primary variables and the secondary variables can be specified as follow:

34 24 <The prim ary variable> <The secondary variable> By substituting the (2.2) into the (2.1) the following nonlinear mixed finite element equations can be obtained (2.3) Note that the boldface letters are used to indicate nonlinear terms. Above mixed finite element equations can be rewritten as algebraic matrix form by collecting coefficients of the unknowns in the form of the coefficient matrix and the rest of the terms as the force vector [4] as follow:

35 25. (2.4) where the of the equation denotes that the coefficient matrix is the function of the unknowns. The sub-matrices and the specific terms of the force vectors are given as follow: 1 2. (2.5) zero. All the sub-matrices and the force vectors which are not specified in the (2.5) are

36 Model II of Beam Bending The governing equations of the EBT which were derived in the chapter I can be also used for the development of the Model II. The Model II includes displacements (i.e. and ) and the generalized resultants (i.e. and ) while the Model I included in addition to those. Thus total number of the independent variables is 4 in the Model II. By eliminating the shear resultant from the governing equations of the EBT following equations can be obtained (2.6) Weighted Residual Statements of Model II The equations given in the (2.6) are mathematically equivalent to the equations given in the (1.7) and (1.9) which were used for the Model I but the effect of the elimination of the can be observed both in the equation solving procedure and in the result of the numerical analysis since it affects both the symmetry of the tangent matrix[1] and the accuracy of the solutions compared with other models. For the Model II the following weighted residual statement can be made (2.7)

37 27 where and. Variables with superscripted a (i.e. and ) denote approximated variables W 1 4 is the i th weight function of the i th weighted residual statement. and are the start and the end global coordinate of the element. The boundary terms and were obtained in the different forms compared with those of the Model I but the physical meaning are the same Finite Element Equations of Model II All of the variables can be approximated by using the Lagrange type interpolations and the weigh functions can be chosen to be as. (2.8) The primary and the secondary variables can be specified as follow:

38 28 <The primary variable> <The secondary variable> By substituting the equation (2.8) into the equation (2.7) the following nonlinear mixed finite element equations of Model (II) can be obtained (2.9) 0. The mixed finite element equations of Model II can be rewritten in algebraic matrix form by collecting the coefficients of the unknowns. Note that the Model II contains 4 variables as unknowns. Thus the size of the of the Model II can be reduced compared with that of the Model I.

39 29. (2.10) The sub-matrices and the specific terms of the force vectors can be given by 1 2. (2.11) The sub-matrices and sub-vectors which are not specified above are zero. 2.3 Model III of Beam Bending In displacement based model of the EBT the slope ( / ) was included as a primary variable with the use of the Hermite type interpolation for the vertical deflection. Because has a physical meaning one can include it as an independent variable. This idea was proposed by Reddy by considering the following

40 30 equation as one of the governing equation of the EBT. It can be seen that the linear part of the coefficient matrix of the Model I given in the (2.4) and (2.5) cannot be symmetry and the tangent matrix of the Model I which will be discussed in the chapter III cannot be symmetry either. In the computational point of view the symmetry of the coefficient matrix of the algebraic equation system is very important because of the computational cost. (2.12) By replacing every in the EBT equations we can obtain the following differential equations which can be modified for the Model III (2.13) Weighted Residual Statements of Model III With the equations of the (2.13) the weighted residual statements of the Model III of beam bending can be written as follow.

41 (2.14) where and Finite Element Equations of Model III In a sense the weighted residual statement of the (2.14) can be fully qualified for the Galerkin method[4] because each of the integral equation represents the work done in virtual work sense with the following choice of the weight functions. The approximations of the variables and the chosen weight functions are as follow:. (2.15)

42 32 The primary and the secondary variables can be specified as follow: <The primary variable> <The secondary variable> By substituting the equation (2.15) into the (2.14) the following nonlinear mixed finite element equations can be given by (2.16)

43 33 The mixed finite element equations of the Model III can be rewritten as the algebraic matrix form by collecting the coefficients of the unknowns. Note that the Model III contains 6 variables as unknowns. Thus the size of the is the largest among three mixed EBT models.. (2.17) The sub-matrices and the specific terms of the force vectors are given as follow: 1 2. (2.18)

44 Model IV of Beam Bending Weighted Residual Statements of the Model IV The governing equations of the TBT are used to develop the Model IV. It can be shown that the mixed Model IV is equivalent to the Model III regarding the numbers and the dimensions of nodal variables. By using the governing equations of the TBT Model IV includes 6 variables as independent variables (i.e. and ). With the equilibrium equation and the resultants equations of the TBT following weighted residual statements can be obtained. where (2.19).

45 Finite Element Equations of Model IV Since no derivative of any variable is involved as a nodal unknown the Lagrange type interpolation function[4] should be used for the approximations of variables of the Model IV. The approximations of the variables and the chosen weight functions are given as follow:. (2.20) The primary and the secondary variables can be specified as follow: <The primary variable> <The secondary variable> Q N dw dx The finite element equations of the Model IV can be obtained by substituting approximations and the chosen weight functions of the (2.20) into the (2.19).

46 (2.21) From the (2.21) the following matrix form of the algebraic equations can be obtained.. (2.22)

47 37 where 1 2. (2.23) All of the sub-matrices and sub-vectors which are not specified above are zero. 2.5 Lagrange Type Beam Finite Elements For present study the mixed formula allows the use of the Lagrange type interpolation functions[4] for the approximations of every variable. Here beam elements that were used for the computer implementation and the numerical analysis are

48 38 mentioned. For the beam problems the Lagrange types of linear quadratic and cubic elements were used. The geometry of the elements and the locations of associated interpolations are given in the Fig (a) linear element (b) a quadratic element (c) a cubic element Fig. 2.1 Node number and local coordinate of the line elements of the Lagrange family. Very well Known interpolation property [2] is known as partition of unity which can be written as So by considering the interpolation properties we can derive interpolation functions which are associated with the nodal points with given set of polynomials. And the specific interpolation functions associated with the nodal points are as follow: - Linear interpolation functions (2.24) 2

49 39 - Quadratic interpolation functions (2.25) 2 - Cubic interpolation functions (2.26)

50 40 CHAPTER III DEVELOPMENT OF PLATE BENDING MODELS The 1-D beam bending problems which were discussed in the chapter II can be extended to the 2-D plate bending problems with simple modifications. Two CPT models and two FSDT models are developed. To clarify the developing procedure the governing equations of the CPT and the FSDT are brought from the chapter I. <The governing equations of the CPT> and 0 0 (1.26) (1.28)

51 41 <The governing equations of the FSDT> (1.33) and (1.36)

52 Model I of Plate Bending Weighted Residual Statements of Model I The governing equations of the CPT which were derived in the chapter I are modified to develop the plate bending Model I. Total eleven variables i.e. u v w N N N V V M M and M are treated as independent variable in the plate Model I. The Green-Gauss theorem[15] can be used to obtain the boundary terms when the integration by parts are conducted. The weighed residual statements of the Model I can be written as follow:

53 (3.1) where Γ is the boundary of the element region Ω and variables with superscripted letter a denote the approximated variables. The denotes the unit normal vector of the i-direction where i = x y Finite Element Equations of Model I With the weighted residual statements given in the (3.1) we can develop the finite element the Model I of the plate bending by approximating variables with known interpolation functions and unknown nodal values. But the choice of the known interpolation functions is not arbitrary. By the same reason that we discussed while developing the EBT models in the chapter II the Lagrange type of the interpolation functions should be used for the approximations of the all variables of the Model I because no derivative of the variable is involved as nodal unknown. To develop the finite element model based on the displacement model[1] one should approximate the vertical displacement with the conforming [1 4] or the nonconforming[1 4] type of the Hermite interpolation functions because displacement based model includes the derivatives of the as nodal values[1]. But for the current Model I only continuity of the variables is required. Thus it allows use of linear interpolation functions as the minimum.

54 44. (3.2) The primary variables and the secondary variable of the Model I can be specified as follow: <The primary variable> <The secondary variable>

55 45 By substituting the equation (3.2) into the (3.1) the finite element Model I of the plate bending can be obtained as follow:

56 (3.3) Above equations given in the (3.3) can be rewritten in the algebraic matrix equation in the following form. (3.4) where

57

58 48.. (3. 5) The rest of sub coefficient matrices and force vectors which are not specified in the (3. 5) are zero. 3.2 Model II of Plate Bending The shear resultant and can be eliminated by substituting the forth and the fifth equilibrium equations into the third equilibrium equation of the CPT. By doing this the symmetry of the linear portion of the coefficient matrix can be achieved and the symmetry of the tangent matrix[1] which will be discussed in the chapter IV can be also archived Weighted Residual Statements of Model II With the equations of the Model II of the plate bending the following weighted residual statements can be obtained.

59 (3.6) Finite Element Equations of Model II All variables can be approximated with the Lagrange type interpolation functions for the same reason which was discussed with the Model I.

60 50. (3.7) By substituting the (3.7) into the (3.6) the finite element model of the Model II can be obtained as follow: 0 0 0

61 (3.8) Above equation (3. 8) can be rewritten as the matrix form of the algebraic equation as follow: K U U K K K K K K K F K F K F. (3.9)

62 52 where

63 53.. (3.10) The rest of the sub coefficient matrices and the force vectors which are not specified in the (3.10) are equal to zero. 3.3 Model III of Plate Bending Weighted Residual Statements of Model III To develop the Model III of the plate bending the governing equations of the FSDT are modified. The Model III will include the shear rotations and to account

64 54 for the shear deformations. The weighed residual statements of this model can be made as follow: (3.11)

65 55 where and denote the element region and the boundary of the element respectively Finite Element Equations of Model III For the Model III the Lagrange type of interpolation functions can be used to approximate the variables. The weight functions and the approximations of the variables can be chosen as follow:. (3.12)

66 56 follow : The primary variables and the secondary variable of the Model III can be specified as <The primary variable> <The secondary variable> The finite element equations can be obtained by substituting the approximations and the weight functions of the (3.12) into the (3.11)

67 (3.13)

68 58 the form of Above equations of the (3.14) can be rewritten as the algebraic matrix equation in K U U K K K K K K F K F K F. (3.14) The specific sub coefficient matrices can be obtained from the (3.13) as follow:

69

70 60 F ψ n N n N ds F ψ n N n N ds F ψ qxdxdy ψ V n V n ds Ω F ψ M n M n ds F ψ M n M n ds (3.15) The rest of sub coefficient matrices and force vectors which are not specified in the (3.15) are equal to zero. 3.4 Model IV of Plate Bending The membrane resultants ( i.e. and ) can be eliminated by substituting the resultant equations of the membrane resultants into the equilibrium equations. By eliminating the in plane force resultants( and ) the size of the coefficient can be reduced while the effect of it will be discussed in the numerical analysis parts of the plate bending Weighted Residual Statements of Model IV The weighted residual statements of the Model IV of the plate bending are given as follow:

71 (3.16)

72 62 The primary variables and the secondary variable of the Model IV can be specified as follow: <The primary variable> <The secondary variable> Finite Element Equations of Model IV Ten weigh functions and the approximations of ten independent variables can be chosen as the Lagrange type interpolation functions as follow:

73 63. (3.17) By substituting the (3.17) into the (3.16) the finite element equations of the Model IV can be obtained as follow:

74

75 (3.18) Above equations of the (3.18) can be rewritten in the algebraic matrix equation by the form of K u U K K K K K K F K F K F. (3.19) where 1 2

76

77 67 (3.20) The rest of the sub coefficient matrices and the force vectors which are not specified above are equal to zero. 3.5 Lagrange Type Plate Finite Elements For present study the mixed formula allows a use of the Lagrange type of interpolation functions[4]. The elements that were used for the computer implementation and the numerical analysis are discussed. For the plate problems 4-node and 9-node Lagrange type of elements were used. The geometry of the elements and the locations of associated interpolations are given in the Fig. 3.1.

78 68 (a) 4-node linear element (b) 9-node quadratic element Fig. 3.1 Node number and local coordinate of the rectangular elements of the Lagrange family. And the associated interpolation functions are as follow: - 4-node linear element (3.21) nod e quadratic element (3.22)

79 69 CHAPTER IV NONLINEAR EQUATION SOLVING PROCEDURES We obtained matrix form of nonlinear finite element model equations in the previous chapters II and III. In this chapter the non-linear equation solving procedures of equation (2.4) (2.10) (2.18) (2.22) (3.4) (3.9) (3.14) and (3.19) were discussed. These equation solving procedures of the nonlinear mixed finite element models can be generally applied for the developed nonlinear finite element models. The Picard Iteration method[1 21] and the Newton-Raphson method[ ] were used for the present numerical analysis. The solutions obtained from the two different methods can be compared to insure the obtained solutions are well converged one because the converging characteristic may vary from one method to the other but the obtained solutions should essentially be the same. 4.1 Direct Iterative Method The direct method is one of the simplest methods available because this method only requires update of the coefficient matrix with obtained solutions from the previous iteration at each iteration step. After updating the coefficient matrix the equation solving procedure which is related to the obtaining of new iterative solutions is just the same as solving linear algebraic equations. The flow chart[1] of the direct iteration method is given in the Fig (for details see Reddy[1]).

80 70 Fig. 4.1 A flow chart[1] of the direct iteration method Algorithm of Direct Iterative Method The element wise matrix form of non-linear finite element equations that we obtained in chapters II and III i.e. (2.4) (2.10) (2.18) (2.22) (3.4) (3.9) (3.14) and (3.19) can be assembled as global equations [1] which can be written as (4. 1) where the denotes the global unknowns of the assembled equations. Because the coefficient matrix is the function of the unknown we need to evaluate the by using initial guess solutions or the previous iterative solutions[1]. To imply this concept the equation (4.1) can be rewritten as

81 71 1 or 1 (4.2) where is the solution obtained from the r th iteration is the updated coefficient matrix using the previous solutions and is the force vector. Now with the equation (4.2) the simple algorithm can be used to solve the nonlinear finite element equations with the direct method which can be stated as 1 1. (4. 3) Above procedure should be repeated until the solutions of r th iteration and the (r-1) th iteration satisfy the following criterion [1]: (4. 4) where is the tolerance[1]. In this study values of tolerance[1] 0.01 ~0.001 were chosen for the most of the problems. In many cases with small values of tolerance (say ) the iterative solutions may not satisfy the criterion regardless of the iteration numbers.

82 Newton-Raphson Iterative Method Usually the Newton-Raphson iterative method[1 22] shows faster convergence compared with the direct method[1 21]. Also in many cases the tangent matrix can be symmetry even though the coefficient matrix is not. And with the Newton method only the tangent matrix is inverted to get the incremental solutions thus only the symmetry solver can be used still the calculations of the tangent matrix and implement of the equation solving procedure are substantial. The flow chart[1] of the Newton iteration method is given in the Fig Fig. 4.2 A flow chart[1] of the Newton iteration method.

83 Algorithm of Newton-Raphson Iterative Method In the Newton-Raphson iterative method the residual[1] or the imbalance force vector[1] of the (4.1) can be written as. (4.5) With Taylor s expansion the residual R can be expanded to the known solution (i.e. the solution of the previous iteration) (4.6) By omitting all the terms after the third term of the right hand side of the (4.6) and by taking the residual to be zero i.e. Ru 0 we can obtain the following relation. where (4.7) 1. Here we define the tangent stiffness matrix[1] as follow:

84 (4.8) By the substitution of (4.8) into the (4.7) and the inversion of the tangent stiffness matrix we can obtain the increment of the solution ( which can be written as (4.9) where the residual can be computed from the previous iterative solution as follow: (4.10) If we can calculate the tangent matrix from the equation (4.8) the solutions can be updated as 1. (4.11) For the check of the convergence criterion it can be computed by using the increment of the solutions vector i.e. as follow: (4.12)

85 Calculation of Tangent Stiffness Matrices In the Newton-Raphson iterative method it is required to compute tangent coefficient matrices to get the incremental solution described in the (4.9). The original form of the equation (4.9) can be rewritten as the matrix form of the equation as follow: 1 1. (4.13) The component form of this tangent coefficient stiffness matrix and the residual vector can be given as (see Reddy [1]) and (4.14) 1 1. (4.15) By substituting the equation (4.15) into the equation (4.14) the following equation can be obtained (see Reddy[1])

86 76. (4.16) Note that the α denotes the equation number which can be matched to the sub matrix of the α th row in the (4.13) the β denotes β th column in the (4.13) denotes the total numbers of unknown variables and is the number of degree of freedom related to the variable. Repeated indices mean summation Tangent Stiffness Matrices The symmetry of the tangent stiffness matrix in Newton iterative method is very important because most of the computational efforts to find the converged solution after obtaining the linear solution(i.e. solutions obtained with zero initial guess solution[1]) are related with the inversion of tangent matrix. By the equation (4.9) and (4.10) the increment of solution can be obtained by inverting the tangent matrix. The inverse of the coefficient matrix is only needed to get the linear solution. Thus the invert of the coefficient matrix does not required after very first step of the iteration. It can be shown that the first linear solution also can be obtained by using the symmetry solver with the choice of a zero initial guess solution. To discuss the symmetry of the tangent matrix the tangent matrices of the newly developed models were calculated by using the equation given in (4.16). For example the of the EBT Model I can be calculated with the coefficient matrices which can be given by

87 If the degree of the freedom for each variable is we can calculate the by using the equation (4.16) as (4.17) where the variables are given by for respectively. Likewise every specific term of the tangent stiffness matrices of newly developed nonlinear beam and plate bending models can be calculated. From the equation (4.15) we can notice that each of the tangent coefficient matrices is consist of the sum of coefficient matrix and the additional terms. So we can express every tangent matrix as the form of + additional terms. The results of the calculations of the tangent matrix of the each model can be given as follow:

88 78 - Model I of the beam bending (4.18) - Model II of the beam bending (4.19) - Model III of the beam bending (4.20) - Model IV of the beam bending (4.21) - Model I of the plate bending

89 (4.22) - Model II of the plate bending (4.23) - Model III of the plate bending (4.24)

90 80 - Model IV of the plate bending (4.25) Rest components of the tangent matrix which are not specified above is the same as the components of the coefficient matrix which can be written as 0. (4.26) As mentioned the symmetry of the tangent matrix can be obtained in every Model except for the beam Model I and the plate Model I. Because of shear terms (i. e. and

91 81 ) included we cannot expect the symmetry of the tangent matrix in the beam Model I and the plate Model I. Then these two models should be solved by asymmetry solver[1] which is based on the Gauss Eliminations[24] to invert a matrix. The numerical results will be discussed in the chapter V. 4.3 Load Increment Vector In the applications of the direct iterative method in nonlinear finite element analysis of the structural problems the load increment is very critical to get converged nonlinear solution under a large applied distributed load ( ). Without proper increment load the solution may not converge with the direct iterative method. But the Newton iterative method can be applied at the more general range of applied distributed load ( ) without load increment to get the converged solutions while more iterative time is substantial. The details of the load increment will be discussed in chapter V with the numerical results of the examples.

92 82 CHAPTER V NUMERICAL RESULTS In this chapter we will discuss the numerical results of the nonlinear finite element models of the beam and plates bending problems. Comparisons of various models are presented with linear analysis and non linear analysis Numerical Analysis of Nonlinear Beam Bending Description of Problem[1] A beam made of steel ( ) whose geometry is given in the Fig. 5.1 was chosen for the study of the 1D nonlinear analysis. Three different boundary conditions i.e. HH PP and CC were considered to see the performance of the beam. Fig. 5.1 Description of the beam geometry. Three types of boundary conditions under the distributed load are considered with 4 nonlinear beam bending models developed in the chapter II. The descriptions of three boundary conditions are given in the Fig Under evenly distributed load and the given boundary conditions we can use the symmetry part of the beam as a computational domain of the finite element analysis. To use the symmetry part of the beam as the computational domain the mathematical boundary conditions at the middle

93 83 point (i.e. x=l/2) of the beam should be specified. By the geometry of the beam bending under the given boundary conditions and evenly distributed load the mathematical boundary conditions of the middle point can be specified in addition to the boundary conditions of the one edge of the beam as shown in the Fig It can be seen that the current EBT mixed nonlinear Model I and II do not includes the slope i.e. / as a primary variable. It is important to specify either primary or secondary variable as a boundary condition. In the same sense we should specify only the moment 0 0 as a primary variable at the edge of the beam because the slope / is not known there. Thus it is clear that if any specified boundary condition exists one should specify either the primary or the secondary variable at the typical nodal point. This can be clarified by using the pairs of the primary and secondary variables that we classified in the chapter III. With the beam Model IV the shear rotation of the beam cross section was included so it can be specified as shown in the Fig A hinged-hinged(h-h) beam /2 0 /2 0(EBT only) /2 0(TBT only) A pined-pined(p-p) beam /2 0 /2 0(EBT only) /2 0(TBT only) (TBT only) 0 0 (EBT only) A clamped- clamped(c-c) beam /2 0 Fig Symmetry boundary conditions of beams. /2 0(EBT only) /2 0(TBT only)

94 Numerical Results First the results of the mixed models and the displacement based models[1] are compared to see the validness of the solutions. The center deflections of the mixed Model I and IV using eight linear elements (8 L) mesh are presented in the Table 5.1 along with the results of the displacement based Models of the EBT and the TBT using eight linear-hermite(8 LH) and eight linear elements(8 L) respectively. Every converged solution was obtained by using the Newton-Raphson iterative method. The graph of the mixed models and the displacement based models which are given in the Fig. 5.3 shows almost the same results for the two different boundary conditions. With the same eight linear elements mesh the difference of the converged solutions is not considerable. And the difference of the solutions between the TBT and the EBT beams is also negligible. Table 5.1 Comparison of mixed models and displacement based models. The The center deflection w 0 (in) load Mixed nonlinear Models Displacement based nonlinear Models[1] The Model (I) - 8 L The Model (IV) - 8 L EBT - 8 LH TBT - 8 L (psi/in) CC PP CC PP CC PP CC PP (3) (5) (3) (5) (3) (5) (3) (5) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4) (3) (4)

95 85 Deflection w0 (in.) EBT(II)(MIXED) PP 8XL EBT(DSPL) PP 8XLH EBT(II)(MIXED) CC 8XL EBT(DSPL) CC 8XLH Load q0 (psi.) (a) Comparison of the Model I with EBT displacement based Model. Deflection w0 (in.) TBT(DSPL) CC 8XL TBT((MIXED) PP 4XL TBT(DSPL) PP 8XL TBT((MIXED) CC 4XL Load q0 (psi.) (b) Comparison of the Model IV with the TBT displacement Model. Fig A comparison of the non-linear solutions of beams.

96 86 But under the hinged-hinged boundary condition current mixed models showed much better results compared with displacement based models. The displacement Model showed the membrane locking[14]. The membrane locking occurs because of the inconsistent presence the polynomial degree in the approximations. To examine it we consider a hinged-hinged boundary condition with the Model I II and IV. For the hinged-hinged boundary condition total applied load should contribute for the bending of the beam element because there is no horizontal constrain to cause membrane strain i.e.. In the finite element models the strain can be expressed as follow: (see Reddy [1]). (5. 1) To satisfy the physics under the given boundary conditions strain should be zero. But because of the use of polynomial approximations there can be inconsistency[1] of the degree of terms in the strain. Especially in the displacement based model was approximated with linear interpolation function and w was approximated by using cubic interpolations functions. For this typical pair of approximations the degree of the each term in the strain can be given by (5. 2) Thus it is not easy for / to make whole strain term to be zero because it is presented as constant. This phenomenon is very well known drawback of the nonlinear EBT and TBT finite element model.

97 87 And for the TBT models another locking can be observed from the shear strain relations[1] which can be given as. (5. 3) To fix these defect the reduced integrations[1 7] use of consistent approximations and use of higher order interpolations can be used. The effects of the locking with full integration in different models are given in the Table 5.2. Table 5.2 Membrane locking in mixed models and misplacement models. The load The center deflection (in) -HH (q ) TBT(Model IV) 4ⅹL TBT(DSPL) 4ⅹL EBT(I)(Model I) 4ⅹL EBT(DSPL) 4ⅹLH The results presented in the Table 5.2 are showing that the membrane locking can be eliminated by using the mixed nonlinear model in both of the TBT and the EBT beams. Usually the locking can be mitigated by using a more refined mesh but the mixed Model I and IV didn t showed any locking even with 2 linear elements mesh as shown in the Table 5.2.

98 88 But among current mixed models the membrane locking appeared in different levels. For example the comparison of the Model I and the Model II shows that the Model I is showing better performance compared with the Model II. But the Model II is still showing better result compared with the displacement based model. The result of the hinged-hinged(hh) boundary condition of the Model I and II with 2 linear elements mesh is given in the Table 5.3. The graph (a) given in the Fig. 5.4 shows that the locking can be eliminated with the mixed Model I and IV. While the graph (b) shows that the Model II still has membrane locking. Even though the Model II has the membrane locking the effect of it is not significant compared to the displacement based model. Table 5. 3 Effect of the membrane locking in the mode I and II. The load (q ) EBT(I) 2ⅹL The center deflection (in)-hh EBT(II) EBT(II) 2ⅹL 4ⅹL EBT(II) 8ⅹL

99 89 Deflection w0 (in.) 6.0 HH-TBT(MIXED)-4XL HH-TBT(DSPL)-4XL HH-EBT(I)(MIXED)-4XL HH-EBT(DSPL)-4XLH Load q0 (psi.) (a) Comparison of the membrane locking of various Models. Deflection w0 (in.) HH- EBT(I)-2XL HH- EBT(II)-2XL HH- EBT(II)-4XL HH- EBT(II)-8XL Load q0 (psi.) (b) Comparison of the membrane locking in the EBT Model I and Model II Fig. 5.4 A comparison of the membrane locking in various models.

100 90 Next the effect of the length-to-thickness ratio on the deflections is presented in the Tables 5.4 and 5.5. The data in the Table 5.4 is showing that as the beam becomes thicker it acts almost linearly while thin beam shows nonlinearity more strongly. Table 5.4 Effect of the length-to-thickness ratio on the deflections in TBT beam. The load (q ) TBT(L/H=100) TBT(L/H=50) TBT(L/H=25) It can be shown that the differences of the solutions between the TBT and the EBT are negligible when the beam is thin but it is not when the beam is thick. Table 5.5 Comparison of the effect of the length-to-thickness ratio in the EBT and the TBT beams. The load (q ) L/H=10 L/H=100 EBT(I) EBT(II) TBT EBT(I) EBT(II) TBT

101 91 Deflection w0 (in.) CC- TBT(L/H=100) CC- TBT(L/H=50) CC- TBT(L/H=25) Load q0 (psi.) (a) Effect of the length-to-thickness ratio on deflections in the TBT beam Deflection w0 (in.) CC- EBT(I)-(L/H=100) CC- EBT(II)-(L/H=100) CC- TBT-(L/H=100) CC- EBT(I)-(L/H=10) CC- EBT(II)-(L/H=10) CC- TBT-(L/H=10) Load q0 (psi.) (b) Comparison of the effect of length-to-thickness ratio on deflections in the EBT and the TBT beams. Fig. 5.5 Comparison of effect of the length-to-thickness ratio on the beam.

102 92 The model III showed poor performance compared with other newly developed models but with cubic element it showed good accuracy and convergence. Some results of the model III presented in the Table Table 5.6 Comparison of Model III with other mixed models. The load (psi/in) The center deflection w 0 (in) The Model (III) - 2 C The Model (III) - 4 C The Model (I) - 2 C The Model (II) - 2 C CC PP CC PP CC PP CC PP The poor performance of Model III can be explained by the (2.12). 0. (5.5) As discussed with the membrane locking and the shear locking this typical relation created other kind of locking because of the inconsistent approximation for the and. Since we included this relation only in the Model III only Model III showed new kind of locking. But this locking was not fixed with reduced integration when lower order interpolation functions (i.e. linear and quadratic) are used. Only higher order interpolation function with reduced integration showed good results.

103 Numerical analysis of Nonlinear Plate Bending Description of Problem[1] Next we consider a non-linear plate bending problems using the newly developed mixed models in the chapter III. A square plate with the following material properties was considered (5.4) The origin of the coordinate was chosen to be located at the center of the plated. The geometry and the coordinate of the plate are described in the Fig Fig A description of the plate bending problem. As it was discussed in the beam bending problem due to the given boundary conditions and the geometry of the plate and the applied load the boundary conditions of

104 94 the rectangular plate with biaxial symmetry were considered. Here three symmetry boundary conditions were considered with common mathematical boundary condition along the symmetry lines of the quadrant of the plate. The specific boundary conditions are given in the Fig Note that for the SS1 at the singular points i.e. point (5 5) both boundary conditions of y = 5 and x = 5 were specified. Fig Symmetry boundary conditions[1 25] of a quadrant of the square plate Non-dimensional Analysis of Linear Solutions To check the accuracy of the newly developed plate bending models solutions of the new models were compared with those of the existing models [ ] and analytic

105 95 solutions. First the linear solutions of the mixed CPT models will be discussed by comparing the data obtained with displacement based model[1 7]. It can be clearly shown that the linear solution of the Model II is the same as that of the mixed model developed by Reddy[7] because both models includes the same variables(i.e. vertical displacement and bending moments) which are related to the bending of the plate while the Model I includes shear resultants also in addition to those. The comparison of the results of the various models under the simple support I (SS1) and clamped (CC) boundary conditions are given in the Tables 5.7 and 5.8. For the simple support (SS1) boundary condition the Model II showed best accuracy for the center vertical deflection while the Model I provided better accuracy for the center bending moment as shown in the Table 5.7. Table 5.7 Comparison of the linear solution of various CPT Models isotropic (. ) square plate simple supported (SS1). Mesh size Liner (4-node) MODEL I Current Models Model II Mixed (Reddy [7]) Mixed (Herrmann[26]) Hybrid ( Allman[27] ) Center deflection (* equivalent quadratic) 10 / ( Exact solution [7] ) Compatible cubic displacement Model[1] ( * - ) ( - ) (0.4154) (0.4154) (0.4067) (0.4067) (0.4063) (0.4063) (0.4063) (0.4063) Liner (4-node) Center bending moment(equivalent quadratic) 10/ ( Exact solution [7]) ( - ) ( - ) (0.4906) (0.4096) (0.4797) (0.4796) (0.4790) (0.4790) (0.4788) (0.4789)

106 96 Table 5.8 Comparison of the linear solution of various CPT Models isotropic (. 3) square plate clamped (CC). Mesh size Liner (4-node) Model I Current Models Model II Mixed (Reddy [7]) Mixed (Herrmann[26]) Center deflection(* equivalent quadratic) Hybrid ( Allman[27] ) 10 / ( Exact [7] ) Compatible cubic displacement Model[1] (* - ) ( - ) (0.1512) (0.1512) (0.1279) (0.1278) (0.1268) (0.1268) (0.1265) (0.1266) Liner (4-node) C enter bending moment(equivalent quadratic) 10/ ( Exact [7]) ( - ) ( - ) (0.2552) (0.2552) (0.2312) (0.2310) (0.2295) (0.2295) (0.2290) (0.2291) For the clamped (CC) boundary condition the Model I showed best accuracy both for the center vertical deflection and the center bending moment as shown in the Table 5.8. The difference of the solution between Model I and Model II was caused by the presence or absence of the shear resultant in the finite element models. Thus by including the shear resultants (i.e. and ) as nodal values in the CPT mixed finite element model more accurate center bending moment and center vertical deflection were obtained. Next current CPT mixed models were compared with the displacement based model. For the CPT displacement based model non-conforming [4] and the conforming[4] elements should be used because of the continuity requirement of the weak formulation[1]. Current mixed models provided better accuracy when the compatible nine-node quadratic element was used. But the four-node liner element also

107 97 provided acceptable accuracy compared with the non-conforming displacement based model. And for the SS1 boundary condition with the Poisson s ratio 0.25 the Model II also showed better accuracy as it did with 0.3. In both cases stresses obtained from the current mixed model showed better accuracy because the stresses can be directly computed by using bending moment or shear resultant obtained at the node not including any derivative. Stresses in the Table 5.9 were obtained by the following equations. See equation (1.25) of the chapter I for specific terms of the matrix Q (i.e. 126 ) while is the vertical shear resultant of the FSDT. And is the component of the invert of matrix [D] given in the (1.24) and (1.29). (5.6) 2 (5.7) 5/6.. (5.8) Not only vertical deflection but also stresses showed better accuracy under simple supported I (SS1) boundary condition when they were compared with those of the displacement based model. In most of the cases results obtained with 9-node quadratic element presented better accuracy. Results of isotropic plate under SS1 boundary condition are given in the Table 5.9.

108 98 Table 5.9 Comparison of the CPT linear solution with that of the displacement model isotropic (. ) square plate simple supported (SS1). / / / / / / Model I Mesh type Linear (4-node) Quadratic (9-node) Exact[1] Model II DSPL. [1] Mesh type Linear(4-node) and Non-conforming (12 - node) Linear (4-node) and Conforming (16 - node) Exact[1] Table 5.10 Comparison of the CPT linear solution with that of the displacement Model isotropic (. ) square plate clamped(cc). / / / / / / Mesh type Linear (4-node) Quadratic (9-node) Exact Model I Model II Mesh type Linear(4-node) and Non-conforming (12 - node) Linear (4-node) and Conforming (16 - node) Exact DSPL. [1]

109 99 Improvement was noticed with clamped boundary condition. The comparison of the results with isotropic plate ( 0.25) under CC boundary condition are given in the Table Next the numerical results of the Model III and IV are compared with the results of the Reddy s mixed model[1]. The mixed model developed by Reddy included bending moments as independent nodal value in the finite element model while current Model III and IV included vertical shear resultants (i.e. and ) as independent nodal value. Note that the difference between Model III and VI comes from the presence or absence of membrane forces (i.e. and ) in the finite element models. Thus the solution of the linear bending of each model is essentially the same as shown in the Table Table 5.11 Comparison of the current mixed FSDT linear solution with that of the other mixed model (Reddy[7]) with isotropic (. / ) square plate simple supported (SS1). Current Models Mixed Current Models Mixed Mesh size Model(III) Model(IV) (Reddy[7]) Model(III) Model(IV) (Reddy[7]) Liner (4-node) Center deflection 10 / (Exact 0.427[8]) Center bending moment 10/ (Exact 0.479[8]) (* - ) ( - ) ( - ) ( - ) (0.4345) (0.4345) (0.4779) (0.4779) (0.4277) (0.4277) (0.4779) (0.4779) (0.4273) (0.4273) (0.4785) (0.4785) - Also the comparison of the center deflection and stresses of current models with those of the displacement based model is presented in the Table In most of cases current models showed better accuracy for both of the center vertical displacement and stresses.

110 100 Table 5.12 Comparison of the linear solution of the FSDT with isotropic (. / ) square plate simple supported (SS1). / / / / / / / / Model type Mesh type Linear (4-node) Quadratic (9-node) Exact[1] Model (III) Model (IV) DSPL.[1] Non-linear Analysis Total 12 load step was used to see the significance of the non-linearity with the following incremental load parameter vector [1] / (5.9) A tolerance 0.01 was used for convergence in the Newton Raphson iteration scheme. Model I and II was compared with the CPT displacement base model to see its non-linear behavior. In non-linear analysis of the CPT center deflection normal stress and membrane stress were compared with the results of the nonconforming and conforming displacement based models.

111 101 First the center defection of the newly developed models are presented in the Table In every load step converged solution was obtained within 4 iterations. To investigate the effect of reduced integration results of full integration and the reduced integration were presented in the Table In both of the model the locking was not severe and the effect of reduced integration was not significant. Table 5.13 Effect of reduced integration in Model I and II. Center deflection w CPT-(SS1) P MODEL I MODEL II 4x4-Linear 2x2-Quadratic 4x4-Linear 2x2-Quadratic FI RI FI RI FI RI FI RI Then using 8 8 linear and 4 4 quadratic elements non-linear normal stresses and center deflection of the Model I were compared with those of the non-conforming and the conforming displacement based models. The results of vertical deflection ans the normal stress is presented in the Table 5.14

112 102 Table 5.14 Comparison of the center deflection and normal stress of Model I and II with the CPT displacement model. P Center deflection w CPT-(SS3) MODEL I MODEL II DSPL DSPL 8x8-L 4x4-Q 8x8-L 4x4-Q 8X8-CF 8x8-UCF P Normal stresses / CPT-(SS3) MODEL I MODEL II DSPL DSPL 8x8-L 4x4-Q 8x8-L 4x4-Q 8X8-CF 8x8-UCF The non linear load versus deflection and load versus stress graphs are given in the Fig.5.8. Under the SS3 boundary condition both of the vertical deflection and stresses of the Model I and II showed very close value when they are compared with the displacement based model. The normal stresses and the membrane stresses were computed at the and 000 respectively. 9-nodel quadratic element showed closer solutions to that of the displacement based FSDT model.

113 103 center vertical deflection w (in) CPT-(SS3) DSPL 8X8-CF CPT-(SS3) DSPL 8x8-UCF CPT-(SS3) MODEL(I) 8x8L CPT-(SS3) MODEL(I) 4x4Q CPT-(SS3) MODEL(II) 8x8L CPT-(SS3) MODEL(II) 4x4Q load parameter P = q 0 x a 4 / (E 22 x h 4 ) (a) Load verses center deflection Stresses σ xx CPT- DSPL-SS3 8X8-CF-Normal CPT- DSPL-SS3 8X8-UCF-Normal CPT- MODEL(I) -SS3 4X4-Q-Normal CPT- MODEL(I) -SS3 4X4-Q-Membrane load parameter P = q 0 x a 4 / (E 22 x h 4 ) (b) Load verses center normal stress Fig Plots of the membrane and normal stress of Model I II and CPT displacement model under SS3 boundary condition. To see the convergence of the various models center deflections of previously developed models with 2 2 quadratic and 4 4 linear meshes under SS1 and SS3

114 104 boundary conditions were compared. Every model showed good convergence with a tolerance 0.01 except for the Model IV. The Model IV showed acceptable convergence with SS3 boundary condition but with SS1 it took more iteration times to converge than other models. The iterative times taken to get converged solutions of the various models are presented in the table Table 5.15 Comparison of the convergence of Model I II III and IV under the SS1 and SS3 boundary conditions. P Center deflection w (*iteration times to converge) SS1 various models Model (III) Model (IV) Model (I) Model (II) 4x4-L 2x2-Q 4x4-L 2x2-Q 2x2-Q 2x2-Q (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (5) (5) (3) (4) (3) (3) (6) (6) (3) (4) (3) (3) (7) (7) (3) (3) (3) (3) (7) (7) (3) (3) (3) (3) (7) (7) (3) (3) (3) (2) (7) (7) (3) (3) (3) (2) (6) (6) (3) (3) (3) (2) (6) (6) (3) (2) (3) (2) (6) (6) (2) (2) (2) (2) (6) (6) (2) (2) (2) (2) (5) (5) (2) (2) Center deflection w (*iteration times to converge) SS3 various models P Model (III) Model (IV) Model (I) Model (II) 4x4-L 2x2-Q 4x4-L 2x2-Q 2x2-Q 2x2-Q (4) (4) (4) (4) (4) (4) (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (3) (4) (4) (3) (4) (3) (3) (4) (4) (3) (4) (3) (3) (4) (4) (3) (4) (2) (2) (4) (4) (3) (4) (2) (2) (3) (3) (2) (3) (2) (2) (3) (3) (2) (3) (2) (2) (3) (3) (2) (3) (2) (2) (3) (3) (2) (3) (2) (2) (3) (3) (2) (3)

115 105 In the nonlinear analysis of the FSDT the non-linear center deflection normal stress and membrane stress of the Model III were compared with the results of displacement based models. The results are presented in the Table A 4x4 quadratic mesh showed the closest result to the displacement FSDT model s result as shown in the Fig Table 5.16 Comparison of the center deflection and normal stress of Model III with the FSDT displacement model under SS1 and SS3 boundary conditions. P Center deflection w FSDT-Model (III) SS1 SS3 DSPL(SS1) DSPL(SS3) 8x8-L 4x4-Q 8x8-L 4x4-Q 4x4-Q 4x4-Q Normal stresses P / FSDT-Model(III) 4x4Q SS1 SS3 D SPL(SS1) D SPL(SS3)

116 center vertical deflection w 0 (in) FSDT-(SS1) DSPL 4x4Q FSDT-(SS3) DSPL 4x4Q FSDT-(SS1) MODEL(III) 8x8L FSDT-(SS1) MODEL(III) 4x4Q FSDT-(SS3) MODEL(III) 8x8L FSDT-(SS3) MODEL(III) 4x4Q load parameter P = q 0 x a 4 / (E 22 x h 4 ) (a) Load verses center deflection Stresses σ xx FSDT- DSPL-SS1 4X4-Q-Normal FSDT- DSPL-SS3 4X4-Q-Normal FSDT- MODEL(III) -SS1 4X4-Q-Normal FSDT- MODEL(III) -SS1 4X4-Q-Membrane FSDT- MODEL(III) -SS3 4X4-Q-Normal FSDT- MODEL(III) -SS3 4X4-Q-Membrane load parameter P = q 0 x a 4 / (E 22 x h 4 ) (b) Load verses center normal and membrane stress Fig. 5.9 Plots of the center deflection normal and membrane stress of Model III with that of the FSDT displacement model under SS1 and SS3 boundary conditions. To see distributions of the variables other than displacements images of the distribution of each variable are presented in the Fig and The data was post processed inside of each element using 10 gauss points ranging from to in both newly developed models (i.e. Model I and III) and FSDT displacement based

117 107 model. Converged solutions of SS3 at the load parameter P = were used for the post processing Displacement based FSDT Model III (FSDT) - 4 ⅹ4 Q Model I (CPT) - 4 ⅹ4 Q Fig Post processed quadrant images of the variables in various models SS3 with converged solution at load parameter.

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

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

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

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

Generic Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials

Generic Strategies to Implement Material Grading in Finite Element Methods for Isotropic and Anisotropic Materials University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Engineering Mechanics Dissertations & Theses Mechanical & Materials Engineering, Department of Winter 12-9-2011 Generic

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

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

Lecture 15 Strain and stress in beams

Lecture 15 Strain and stress in beams Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME

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

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

Theories of Straight Beams

Theories of Straight Beams EVPM3ed02 2016/6/10 7:20 page 71 #25 This is a part of the revised chapter in the new edition of the tetbook Energy Principles and Variational Methods in pplied Mechanics, which will appear in 2017. These

More information

General elastic beam with an elastic foundation

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

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

ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN XFEM APPROACH

ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN XFEM APPROACH VI International Conference on Adaptive Modeling and Simulation ADMOS 213 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) ACCURATE MODELLING OF STRAIN DISCONTINUITIES IN BEAMS USING AN

More information

Presented By: EAS 6939 Aerospace Structural Composites

Presented By: EAS 6939 Aerospace Structural Composites A Beam Theory for Laminated Composites and Application to Torsion Problems Dr. BhavaniV. Sankar Presented By: Sameer Luthra EAS 6939 Aerospace Structural Composites 1 Introduction Composite beams have

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

CHAPTER 5. Beam Theory

CHAPTER 5. Beam Theory CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions

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

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

Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian

Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:

More information

COPYRIGHTED MATERIAL. Index

COPYRIGHTED MATERIAL. Index Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,

More information

Accepted Manuscript. R.C. Batra, J. Xiao S (12) Reference: COST Composite Structures. To appear in:

Accepted Manuscript. R.C. Batra, J. Xiao S (12) Reference: COST Composite Structures. To appear in: Accepted Manuscript Finite deformations of curved laminated St. Venant-Kirchhoff beam using layerwise third order shear and normal deformable beam theory (TSNDT) R.C. Batra, J. Xiao PII: S0263-8223(12)00486-2

More information

Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5

Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5 Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline

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

Chapter 6 2D Elements Plate Elements

Chapter 6 2D Elements Plate Elements Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda

More information

Chapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14

Chapter 2: Rigid Bar Supported by Two Buckled Struts under Axial, Harmonic, Displacement Excitation..14 Table of Contents Chapter 1: Research Objectives and Literature Review..1 1.1 Introduction...1 1.2 Literature Review......3 1.2.1 Describing Vibration......3 1.2.2 Vibration Isolation.....6 1.2.2.1 Overview.

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

The CR Formulation: BE Plane Beam

The CR Formulation: BE Plane Beam 6 The CR Formulation: BE Plane Beam 6 Chapter 6: THE CR FORMUATION: BE PANE BEAM TABE OF CONTENTS Page 6. Introduction..................... 6 4 6.2 CR Beam Kinematics................. 6 4 6.2. Coordinate

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

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

MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM

MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM S. Jasotharan * and I.R.A. Weerasekera University of Moratuwa, Moratuwa, Sri Lanka * E-Mail: jasos91@hotmail.com,

More information

AN ALGORITHM FOR TOPOLOGY OPTIMIZATION

AN ALGORITHM FOR TOPOLOGY OPTIMIZATION AN ALGORITHM FOR TOPOLOGY OPTIMIZATION OF MULTIBODY SYSTEMS TOHEED GHANDRIZ Master s thesis 2014:E52 Faculty of Science Centre for Mathematical Sciences Numerical Analysis CENTRUM SCIENTIARUM MATHEMATICARUM

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

CIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives

CIVL 7/8117 Chapter 4 - Development of Beam Equations - Part 2 1/34. Chapter 4b Development of Beam Equations. Learning Objectives CIV 7/87 Chapter 4 - Development of Beam Equations - Part /4 Chapter 4b Development of Beam Equations earning Objectives To introduce the work-equivalence method for replacing distributed loading by a

More information

Mechanics of Inflatable Fabric Beams

Mechanics of Inflatable Fabric Beams Copyright c 2008 ICCES ICCES, vol.5, no.2, pp.93-98 Mechanics of Inflatable Fabric Beams C. Wielgosz 1,J.C.Thomas 1,A.LeVan 1 Summary In this paper we present a summary of the behaviour of inflatable fabric

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

Free vibration analysis of beams by using a third-order shear deformation theory

Free vibration analysis of beams by using a third-order shear deformation theory Sādhanā Vol. 32, Part 3, June 2007, pp. 167 179. Printed in India Free vibration analysis of beams by using a third-order shear deformation theory MESUT ŞİMŞEK and TURGUT KOCTÜRK Department of Civil Engineering,

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

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

Introduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA

Introduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC

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

CO-ROTATIONAL DYNAMIC FORMULATION FOR 2D BEAMS

CO-ROTATIONAL DYNAMIC FORMULATION FOR 2D BEAMS COMPDYN 011 ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 5-8 May 011 CO-ROTATIONAL

More information

CIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen

CIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall

More information

Module-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates

Module-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates Module-6: Laminated Composites-II Learning Unit-1: M6.1 M 6.1 Structural Mechanics of Laminates Classical Lamination Theory: Laminate Stiffness Matrix To this point in the development of classical lamination

More information

Toward a novel approach for damage identification and health monitoring of bridge structures

Toward a novel approach for damage identification and health monitoring of bridge structures Toward a novel approach for damage identification and health monitoring of bridge structures Paolo Martino Calvi 1, Paolo Venini 1 1 Department of Structural Mechanics, University of Pavia, Italy E-mail:

More information

Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams

Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams Correction of local-linear elasticity for nonlocal residuals: Application to Euler-Bernoulli beams Mohamed Shaat* Engineering and Manufacturing Technologies Department, DACC, New Mexico State University,

More information

Unit 13 Review of Simple Beam Theory

Unit 13 Review of Simple Beam Theory MIT - 16.0 Fall, 00 Unit 13 Review of Simple Beam Theory Readings: Review Unified Engineering notes on Beam Theory BMP 3.8, 3.9, 3.10 T & G 10-15 Paul A. Lagace, Ph.D. Professor of Aeronautics & Astronautics

More information

ON SIMPLE AND ACCURATE FINITE ELEMENT MODELS FOR NONLINEAR BENDING ANALYSIS OF BEAMS AND PLATES. A Dissertation YETZIRAH YKSYA URTHALER LAPEIRA

ON SIMPLE AND ACCURATE FINITE ELEMENT MODELS FOR NONLINEAR BENDING ANALYSIS OF BEAMS AND PLATES. A Dissertation YETZIRAH YKSYA URTHALER LAPEIRA ON SIMPLE AND ACCURATE FINITE ELEMENT MODELS FOR NONLINEAR BENDING ANALYSIS OF BEAMS AND PLATES A Dissertation by YETZIRAH YKSYA URTHALER LAPEIRA Submitted to the Office of Graduate Studies of Texas A&M

More information

Fundamentals of Linear Elasticity

Fundamentals of Linear Elasticity Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy

More information

ENGN2340 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn

ENGN2340 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn ENGN234 Final Project: Implementation of a Euler-Bernuolli Beam Element Michael Monn 12/11/13 Problem Definition and Shape Functions Although there exist many analytical solutions to the Euler-Bernuolli

More information

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

Nonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess

Nonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable

More information

A consistent dynamic finite element formulation for a pipe using Euler parameters

A consistent dynamic finite element formulation for a pipe using Euler parameters 111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,

More information

An Introduction to the Finite Element Method

An Introduction to the Finite Element Method An Introduction to the Finite Element Method Third Edition J. N. REDDY Department 01 Mechanical Engineering Texas A&M University College Station, Texas, USA 77843 11 Boston Burr Ridge, IL Dubuque, IA Madison,

More information

AERSYS KNOWLEDGE UNIT

AERSYS KNOWLEDGE UNIT -7016 1. INTRODUCTION The scope of this document is to provide a clarification and a deeper understanding of the two different ways to move the mid plane of the element out of the nodal plane. Although

More information

EML4507 Finite Element Analysis and Design EXAM 1

EML4507 Finite Element Analysis and Design EXAM 1 2-17-15 Name (underline last name): EML4507 Finite Element Analysis and Design EXAM 1 In this exam you may not use any materials except a pencil or a pen, an 8.5x11 formula sheet, and a calculator. Whenever

More information

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004 Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic

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

Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method

Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method 9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No

More information

Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering

Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering

More information

Geometry-dependent MITC method for a 2-node iso-beam element

Geometry-dependent MITC method for a 2-node iso-beam element Structural Engineering and Mechanics, Vol. 9, No. (8) 3-3 Geometry-dependent MITC method for a -node iso-beam element Phill-Seung Lee Samsung Heavy Industries, Seocho, Seoul 37-857, Korea Hyu-Chun Noh

More information

Basic Energy Principles in Stiffness Analysis

Basic Energy Principles in Stiffness Analysis Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting

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

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

Direct calculation of critical points in parameter sensitive systems

Direct calculation of critical points in parameter sensitive systems Direct calculation of critical points in parameter sensitive systems Behrang Moghaddasie a, Ilinca Stanciulescu b, a Department of Civil Engineering, Ferdowsi University of Mashhad, P.O. Box 91775-1111,

More information

HIGHER-ORDER THEORIES

HIGHER-ORDER THEORIES HIGHER-ORDER THEORIES THIRD-ORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 Third-Order Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +

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

DAMAGE ANALYSIS OF LAMINATED COMPOSITE BEAMS UNDER BENDING LOADS USING THE LAYER-WISE THEORY. A Dissertation WOOK JIN NA

DAMAGE ANALYSIS OF LAMINATED COMPOSITE BEAMS UNDER BENDING LOADS USING THE LAYER-WISE THEORY. A Dissertation WOOK JIN NA DAMAGE ANALYSIS OF LAMINATED COMPOSITE BEAMS UNDER BENDING LOADS USING THE LAYER-WISE THEORY A Dissertation by WOOK JIN NA Submitted to the Office of Graduate Studies of Texas A&M University in partial

More information

PLAT DAN CANGKANG (TKS 4219)

PLAT DAN CANGKANG (TKS 4219) PLAT DAN CANGKANG (TKS 4219) SESI I: PLATES Dr.Eng. Achfas Zacoeb Dept. of Civil Engineering Brawijaya University INTRODUCTION Plates are straight, plane, two-dimensional structural components of which

More information

Beam Models. Wenbin Yu Utah State University, Logan, Utah April 13, 2012

Beam Models. Wenbin Yu Utah State University, Logan, Utah April 13, 2012 Beam Models Wenbin Yu Utah State University, Logan, Utah 843-4130 April 13, 01 1 Introduction If a structure has one of its dimensions much larger than the other two, such as slender wings, rotor blades,

More information

Continuum Mechanics and the Finite Element Method

Continuum Mechanics and the Finite Element Method Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after

More information

Virtual Work and Variational Principles

Virtual Work and Variational Principles Virtual Work and Principles Mathematically, the structural analysis problem is a boundary value problem (BVP). Forces, displacements, stresses, and strains are connected and computed within the framework

More information

NON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil

NON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil Abstract: NON LINEAR BUCKLING OF COLUMNS Dr. Mereen Hassan Fahmi Technical College of Erbil The geometric non-linear total potential energy equation is developed and extended to study the behavior of buckling

More information

Table of Contents. Preface... 13

Table of Contents. Preface... 13 Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...

More information

ME 475 Modal Analysis of a Tapered Beam

ME 475 Modal Analysis of a Tapered Beam ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory

More information

Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING

Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING CHETTINAD COLLEGE OF ENGINEERING AND TECHNOLOGY-KARUR FINITE ELEMENT ANALYSIS 2 MARKS QUESTIONS WITH ANSWER Prepared by M. GUNASHANKAR AP/MECH DEPARTMENT OF MECHANICAL ENGINEERING FINITE ELEMENT ANALYSIS

More information

DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)

DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS) DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types

More information

3D Elasticity Theory

3D Elasticity Theory 3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.

More information

MITOCW MITRES2_002S10nonlinear_lec05_300k-mp4

MITOCW MITRES2_002S10nonlinear_lec05_300k-mp4 MITOCW MITRES2_002S10nonlinear_lec05_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

BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I

BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION

More information

Chapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements

Chapter 12 Plate Bending Elements. Chapter 12 Plate Bending Elements CIVL 7/8117 Chapter 12 - Plate Bending Elements 1/34 Chapter 12 Plate Bending Elements Learning Objectives To introduce basic concepts of plate bending. To derive a common plate bending element stiffness

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

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

Stiffness Matrices, Spring and Bar Elements

Stiffness Matrices, Spring and Bar Elements CHAPTER Stiffness Matrices, Spring and Bar Elements. INTRODUCTION The primary characteristics of a finite element are embodied in the element stiffness matrix. For a structural finite element, the stiffness

More information

Biomechanics. Soft Tissue Biomechanics

Biomechanics. Soft Tissue Biomechanics Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic

More information

Chapter 6: Cross-Sectional Properties of Structural Members

Chapter 6: Cross-Sectional Properties of Structural Members Chapter 6: Cross-Sectional Properties of Structural Members Introduction Beam design requires the knowledge of the following. Material strengths (allowable stresses) Critical shear and moment values Cross

More information

Chapter 2: Deflections of Structures

Chapter 2: Deflections of Structures Chapter 2: Deflections of Structures Fig. 4.1. (Fig. 2.1.) ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 1 (2.1) (4.1) (2.2) Fig.4.2 Fig.2.2 ASTU, Dept. of C Eng., Prepared by: Melkamu E. Page 2

More information

A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE AND MOVING-POLE APPROACHES IN STATIC AND DYNAMIC ANALYSIS OF NON-LINEAR SPATIAL BEAM STRUCTURES

A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE AND MOVING-POLE APPROACHES IN STATIC AND DYNAMIC ANALYSIS OF NON-LINEAR SPATIAL BEAM STRUCTURES European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS 212) J. Eberhardsteiner et.al. (eds.) Vienna, Austria, September 1-14, 212 A NOTE ON RELATIONSHIP BETWEEN FIXED-POLE

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

Unit 18 Other Issues In Buckling/Structural Instability

Unit 18 Other Issues In Buckling/Structural Instability Unit 18 Other Issues In Buckling/Structural Instability Readings: Rivello Timoshenko Jones 14.3, 14.5, 14.6, 14.7 (read these at least, others at your leisure ) Ch. 15, Ch. 16 Theory of Elastic Stability

More information

Lab Exercise #5: Tension and Bending with Strain Gages

Lab Exercise #5: Tension and Bending with Strain Gages Lab Exercise #5: Tension and Bending with Strain Gages Pre-lab assignment: Yes No Goals: 1. To evaluate tension and bending stress models and Hooke s Law. a. σ = Mc/I and σ = P/A 2. To determine material

More information

MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4

MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 MITOCW MITRES2_002S10nonlinear_lec15_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

HIGHER-ORDER THEORIES

HIGHER-ORDER THEORIES HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain

More information

Aircraft Structures Kirchhoff-Love Plates

Aircraft Structures Kirchhoff-Love Plates University of Liège erospace & Mechanical Engineering ircraft Structures Kirchhoff-Love Plates Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/ Chemin

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

Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS

Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms

More information

The Finite Element Method for Mechonics of Solids with ANSYS Applicotions

The Finite Element Method for Mechonics of Solids with ANSYS Applicotions The Finite Element Method for Mechonics of Solids with ANSYS Applicotions ELLIS H. DILL 0~~F~~~~"P Boca Raton London New Vork CRC Press is an imprint 01 the Taylor & Francis Group, an Informa business

More information

Lecture 8. Stress Strain in Multi-dimension

Lecture 8. Stress Strain in Multi-dimension Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]

More information

VIBRATION PROBLEMS IN ENGINEERING

VIBRATION PROBLEMS IN ENGINEERING VIBRATION PROBLEMS IN ENGINEERING FIFTH EDITION W. WEAVER, JR. Professor Emeritus of Structural Engineering The Late S. P. TIMOSHENKO Professor Emeritus of Engineering Mechanics The Late D. H. YOUNG Professor

More information