Interactive Analysis of General Beam Configurations Using Finite Element Methods and Javascript by CHRISTOPHER HERNANDEZ Presen

Total Page:16

File Type:pdf, Size:1020Kb

Interactive Analysis of General Beam Configurations Using Finite Element Methods and Javascript by CHRISTOPHER HERNANDEZ Presen Interactive Analysis of General Beam Configurations using Finite Element Methods and JavaScript By CHRISTOPHER HERNANDEZ Presented to the Faculty of the Graduate School of The University of Texas at Arlington in Partial Fulfillment of the Requirements for the Degree of MASTER OF SCIENCE IN AEROSPACE ENGINEERING THE UNIVERSITY OF TEXAS AT ARLINGTON December 2016 i Copyright © by Christopher Hernandez 2016 All Rights Reserved ii Acknowledgements I would like to thank my family for always being there for me and supporting me through all the years I’ve been in school. Without them I would not have been able to achieve all that I have up to this point. I would also like to express my gratitude to Dr. Lawrence, who gave me the opportunity to not only learn new skills but was also there to guide and help me make my work the best that it could be, and was instrumental in shaping this work into what it is today. Finally, I would like to thank my friends who I have been with me since my undergraduate studies. November 21, 2016 iii Abstract DEVELOPMENT OF FINITE ELEMENT BEAM ANALYSIS APPLICATION USING JAVASCRIPT Christopher Hernandez, MS The University of Texas at Arlington, 2016 Supervising Professor: Kent Lawrence Advancements in computer technology have contributed to the widespread practice of modelling and solving engineering problems through the use of specialized software. The wide use of engineering software comes with the disadvantage to the user of costs from the required purchase of software licenses. The creation of accurate, trusted, and freely available applications capable of conducting meaningful analysis of engineering problems is a way to mitigate to the costs associated with every-day engineering computations. Writing applications in the JavaScript programming language allows the applications to run within any computer browser, without the need to install specialized software, since all internet browsers are equipped with virtual machines (VM) that allow the browsers to execute JavaScript code. iv The objective of this work is the development of an application that performs the analysis of a completely general beam through use of the finite element method. The app is written in JavaScript and embedded in a web page so it can be downloaded and executed by a user with an internet connection. This application allows the user to analyze any uniform or non-uniform beam, with any combination of applied forces, moments, distributed loads, and boundary conditions. Outputs for this application include lists the beam deformations and slopes, as well as lateral and slope deformation graphs, bending stress distributions, and shear and a moment diagrams. To validate the methodology of the GBeam finite element app, its results are verified using the results from obtained from two other established finite element solvers for fifteen separate test cases. v Table of Contents Acknowledgements ...................................................................................................................................... iii Abstract ........................................................................................................................................................ iv List of Illustrations ..................................................................................................................................... viii List of Tables ................................................................................................................................................ x 1. Introduction ............................................................................................................................................... 1 2. Installed Software Licensing ..................................................................................................................... 1 2.1 SOFTWARE LICENSING BEFORE THE INTERNET .............................................................................. 1 2.2 Types of Licenses ............................................................................................................................... 2 3. JavaScript Programming Language .......................................................................................................... 5 3.1 JavaScript: A Brief Overview ............................................................................................................. 5 3.2 Developing Computational Applications with JavaScript .................................................................. 6 4. The Finite Element Method .................................................................................................................... 10 4.1 Overview of FEM: History, Advantages, and Concepts ................................................................... 10 4.2 Beam Models: Euler-Bernoulli Beam ............................................................................................... 14 4.3 Formulation of Beam Element Stiffness Matrix: Euler-Bernoulli .................................................... 16 4.4 Formulation of Beam Global Stiffness Matrix: Euler-Bernoulli ...................................................... 24 5. JavaScript Finite Element Beam Analysis .............................................................................................. 32 5.1 An Application in Three Parts .......................................................................................................... 32 5.1.1 Input Spreadsheet: The Main Menu ........................................................................................... 33 5.1.2 Solution Output Window ........................................................................................................... 41 5.1.3 Plot Output Window .................................................................................................................. 42 5.2 Sorting Data Into Arrays ................................................................................................................... 43 5.3 Creation of Element & Global Stiffness Matrices ............................................................................ 47 5.3.1 Stiffness Matrix Creation ........................................................................................................... 47 5.3.2 Effects of Boundary Conditions On Global Stiffness Matrix .................................................... 50 5.4 Spring Support .................................................................................................................................. 63 5.5 Applied Forces and Moments ........................................................................................................... 64 5.5.1 Application of Distributed Forces .............................................................................................. 65 5.6 Beam Refinement .............................................................................................................................. 70 6. Solution and Outputs ............................................................................................................................... 71 6.1 Displacements ................................................................................................................................... 72 vi 6.2 Displacement Output Window .......................................................................................................... 75 6.3 Plot Output Window ......................................................................................................................... 77 6.3.1 Lateral Displacement Diagram .................................................................................................. 78 6.3.2 Slope Diagram ........................................................................................................................... 79 6.3.3 Shear and Moment Diagram ...................................................................................................... 81 6.3.4 Bending Stress Diagram ............................................................................................................. 88 7 Output Comparison .................................................................................................................................. 89 8 Conclusions ............................................................................................................................................ 105 9 Future Work ........................................................................................................................................... 107 Appendix A- Additional Application Images ........................................................................................... 111 Appendix B- Beam Problem Outputs and Comparisons........................................................................... 113 vii List of Illustrations Figure 1 Circle Area Application Example 1 ................................................................................................ 6 Figure 2 Circle Area Application Example 1 Script ..................................................................................... 7 Figure 3 Circle Area Application Example 3 ................................................................................................ 7 Figure 4 Circle Area Application Example 2 Script ..................................................................................... 8 Figure 5 Client-Side Application [7] .........................................................................................................
Recommended publications
  • Structural Behavior of Strand Reinforced
    STRUCTURAL BEHAVIOR OF STRAND REINFORCED CONCRETE MEMBERS WITHOUT PRESTRESSING by Alaa Maher Al Hawarneh A Thesis Presented to the Faculty of the American University of Sharjah College of Engineering in Partial Fulfillment of the Requirements for the Degree of Master of Science in Civil Engineering Sharjah, United Arab Emirates December 2016 © 2016 Alaa Maher Al Hawarneh. All rights reserved. Approval Signatures We, the undersigned, approve the Master’s Thesis of Alaa Maher Al Hawarneh. Thesis Title: Structural Behavior of Strand Reinforced Concrete Members without Prestressing Signature Date of Signature (dd/mm/yyyy) Dr. Sami W. Tabsh Professor, Department of Civil Engineering Thesis Advisor Dr. Farid H. Abed Associate Professor, Department of Civil Engineering Thesis Committee Member Dr. Basil Darras Associate Professor, Department of Mechanical Engineering Thesis Committee Member Dr. Robert John Houghtalen Head, Department of Civil Engineering Dr. Mohamed El Tarhuni Associate Dean, College of Engineering Dr. Richard T. Schoephoerster Dean, College of Engineering Dr. Khaled Assaleh Interim Vice Provost for Research and Graduate Studies Acknowledgement In the name of Allah, I would like to start my thesis by appreciating everybody who was with me during the research program. The first and foremost thanks go to Allah for His help and blessing during the whole period. I want then to thank everybody who supported me in order to accomplish this work. It has been a great period of intensive learning on the scientific as well as the personal level. I would like first to express my sincere gratitude to my supervisor who gave me the full educational and technical support in order to conduct my research and finish my thesis.
    [Show full text]
  • Stiffness Matrix Method Beam Example
    Stiffness Matrix Method Beam Example Yves soogee weak-kneedly as monostrophic Christie slots her kibitka limp east. Alister outsport creakypreferentially Theobald if Uruguayan removes inimitablyNorm acidify and or explants solder. harmfully.Berkeley is Ethiopic and eunuchized paniculately as The method is much higher modes of strength and beam example Consider structural stiffness methods, and beam example is to familiarise yourself with beams which could not be obtained soluti one can be followed to. With it is that requires that might sound obvious, solving complex number. This method similar way load displacement based models beams which is stiffness? The direct integration of these equations will result to exact general solution to the problems as the case may be. Which of the following facts are true for resilience? Since this is a methodical and, or hazard events, large compared typical variable. The beam and also reduction of a methodical and some troubles and consume a variable. Gthe shear forces. Can you more elaborate than the fourth bullet? This method can find many analytical program automatically generates element stiffness matrix method beam example. In matrix methods, stiffness method used. Proof than is defined as that stress however will produce more permanent extension of plenty much percentage in rail gauge the of the standard test specimen. An example beam stiffness method is equal to thank you can be also satisfied in? The inverse of the coefficient matrix is only needed to get thcoefficient matrix does not requthat the first linear solution center can be obtained by using the symmetry solver with the lution. Appendices c to stiffness matrix for example.
    [Show full text]
  • Unit 17 Approximate Methods of Analysis
    UNIT 17 APPROXIMATE METHODS OF ANALYSIS Structure 17.1 Introduction Objectives 17.2 Substitute Frame Method 17.3 Horizontal or Lateral Loading 17.3.1 Portal Method 17.3.2 Cantilever Method 17.4 Mixed Approximate Methods Based on Existing Exact and Approximate Methods 17.5 Summary 17.6 Keywords 17.7 Answers to SAQs INTRODUCTION In statically determinate structures, the analysis of a structure is made without the knowledge of the cross-sectional dimensions of the members, whereas cross section and members' elastic properties are prerequisite for the analysis of indeterminate structures. Therefore, it becomes necessary to perform approximate analysis for deciding preliminary selection of member sizes. Approximate methods are also useful for quick checking the results of exact analysis. It is also a powerful tool for carrying out spontaneous scrutiny of design, which involves large amount of analysis. In approximateanalysis, the statically indeterminate structure is simprified to a statically determinate structure by making suitable assumptions based on the experience of the analyst. Then, the analysis is carried out by using principles of statics. The validity of the results is based upon assumptions made in the analysis. In this unit, substitute frame method is described in Section 17.2. The approximate methods such as portal and cantileyer method are explained in Section 17.3.1 and 17.3.2, respectively. Objectives After studying this unit, you should be able to conceptualise and appreciate the use of approximate methods of analysis, find out internal forces, i.e. axial force, shear force and bending moment in any member of a building frame by method of substitute frame method due to vertical loading, compute internal forces in beams and columns of a plane frame subjected to lateral loading, by using portal and cantilever methods, and find out shear and moment due to transverse loading by generalising and approximating the facts based on existing exact methods.
    [Show full text]
  • A New Stiffness Matrix for Analysis of Flexural Line Continuum 1, O
    The International Journal Of Engineering And Science (Ijes) ||Volume|| 2 ||Issue|| 2 ||Pages|| 57-62 ||2013|| Issn: 2319 – 1813 Isbn: 2319 – 1805 A New Stiffness Matrix for Analysis Of Flexural Line Continuum 1, 2, 3, O. M. Ibearugbulem, L. O. Ettu, J. C. Ezeh 1,2,3,Civil Engineering Department, Federal University of Technology, Owerri, Nigeria ---------------------------------------------------------Abstract------------------------------------------------------ The analysis of line continuum by equilibrium approach to obtain exact results requires direct integration of the governing equations and satisfying the boundary conditions. This procedure is cumbersome and 4x4 stiffness matrix system with 4x1 load vector has been used as practical simplification for a long time. The problem with this system is that it is inefficient and does not give exact results, especially for classical studies that analyse only one element as a whole. This paper presents 5x5 stiffness matrix system as a new stiffness system to generate exact results. The stiffness matrix was developed using energy variational principle. The approach differs from the 4x4 stiffness matrix approach by considering a deformable node at the centre of the continuum which brings the number of deformable nodes to five. By carrying out direct integration of the governing differential equation of the line continuum, a five-term shape function was obtained. This was substituted into the strain energy equation and the resulting functional was minimized, resulting in a 5x5 stiffness matrix used herein. The five term shape function was also substituted into load-work equation and minimized to obtain load vector for uniformly distributed load on the continuum. Classical studies of beams of four different boundary conditions under uniformly distributed load and point load were performed using the new 5x5 stiffness matrix with 5x1 load vector for one element analysis.
    [Show full text]
  • Euler-Bernoulli Beam Theory
    Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 24 February 2021 doi:10.20944/preprints202102.0559.v1 EULER−BERNOULLI BEAM THEORY USING THE FINITE DIFFERENCE METHOD Article EULER−BERNOULLI BEAM THEORY: FIRST-ORDER ANALYSIS, SECOND-ORDER ANALYSIS, STABILITY, AND VIBRATION ANALYSIS USING THE FINITE DIFFERENCE METHOD Valentin Fogang Abstract: This paper presents an approach to the Euler−Bernoulli beam theory (EBBT) using the finite difference method (FDM). The EBBT covers the case of small deflections, and shear deformations are not considered. The FDM is an approximate method for solving problems described with differential equations (or partial differential equations). The FDM does not involve solving differential equations; equations are formulated with values at selected points of the structure. The model developed in this paper consists of formulating partial differential equations with finite differences and introducing new points (additional points or imaginary points) at boundaries and positions of discontinuity (concentrated loads or moments, supports, hinges, springs, brutal change of stiffness, etc.). The introduction of additional points permits us to satisfy boundary conditions and continuity conditions. First-order analysis, second-order analysis, and vibration analysis of structures were conducted with this model. Efforts, displacements, stiffness matrices, buckling loads, and vibration frequencies were determined. Tapered beams were analyzed (e.g., element stiffness matrix, second-order analysis). Finally, a direct
    [Show full text]
  • Structural Analysis Program Using Direct Stiffness Method for the Design of Reinforced Concrete Column
    Journal of Electrical Engineering and Computer Sciences ISSN : 2528 - 0260 Vol. 2, No.1, June 2017 STRUCTURAL ANALYSIS PROGRAM USING DIRECT STIFFNESS METHOD FOR THE DESIGN OF REINFORCED CONCRETE COLUMN AHMAD FAZA AZMI Engineering Faculty, Bhayangkara University Of Surabaya Jalan Ahmad Yani 114, Surabaya 60231, Indonesia E-mail: [email protected] ABSTRACT Most of the analyzes and calculations in civil engineering have been using computer. But many commercial structural analysis programs are relatively expensive. Hence one of the alternative way is using a self developed structural analysis program. The computer program is constructed using Visual Basic 6 programming language. Using direct stiffness method and theory of ACI 318-02, a numerical procedure along with a computer program was developed for the structural analysis and design of reinforced concrete member. From the direct stiffness method,in terms of structural analysis, the output that obtained from the program are internal element forces and support reactions node displacements. From ACI 318-02 codes, in terms of reinforced concrete member analysis, the program generate the amount of longitudinal and lateral reinforcement of concrete member based on the output in the previous structure analysis program. The program also reveals the biaxial bending-axial interaction diagrams in each direction of rectangular and circular column. The interaction diagram of the columns is a visual representation to determine the maximum axial load and moment exceeded the capacity of the column. Keywords: Structural, Analysis, Reinforced, Concrete, Column 1. INTRODUCTION Structural analysis and design is known to human beings since early civilizations. Structural analysis and design of a reinforced concrete column with biaxial bending and axial force is complicated without the assistance of a computer program.
    [Show full text]
  • DIRECT STIFFNESS METHOD for ANALYSIS of SKELETAL STRUCTURES Dr. Suresh Bhalla
    DIRECT STIFFNESS METHOD FOR ANALYSIS OF SKELETAL STRUCTURES http://web.iitd.ac.in/~sbhalla/cvl756.html Dr. Suresh Bhalla Professor Department of Civil Engineering Indian Institute of Technology Delhi 1 The audio of this lecture is based on classroom recording from CVL 756 class for batch 2019-20. Use earphone/ amplifier for best audio experience. 2 SLOPE DEFLECTION METHOD Applied moment M θ Carry over moment 4EI 2EI L 6EI 6EI L 2 L2 L 12EI 3 6EI L 2 L 6EI 12EI L2 L3 3 Generalized slope deflection equations 2EI 3( ) θ 2 1 M1 2 M1 21 2 2 L L θ1 2EI 3( 2 1) 1 M2 M 2 1 2 2 L L F F1 2 M M F F 1 2 1 2 L 2EI 6( ) 2EI 6( 2 1) 2 1 F 3 3 F2 31 3 2 1 2 1 2 2 L L L L 4 PUT IN MATRIX FORM 12EI 6EI 12EI 6EI 3 2 3 2 F L L L L 1 6EI 4EI 6EI 2EI 1 M1 2 2 1 • Symmetric L L L L F 12EI 6EI 12EI 6EI 2 2 • K +ive 3 2 3 2 ii M L L L L 2 6EI 2EI 6EI 4EI 2 L2 L L2 L Force Stiffness Displacement vector matrix vector (symmetric) 5 F1 K11 K12 . K1n 1 1 2 3 F K K . K 2 21 22 2n 2 . . . Fn Kn1 Kn1 . Knn n n = Degrees of freedom (DKI) th Fi = Force along i degree of freedom th i = Displacement along i degree of freedom th j col.
    [Show full text]
  • Structural Analysis of Continuous Beam by Slope Deflection
    Shear and Moment Diagrams for a Continuous Beam The slope-deflection method is used to determine the shear and moment diagram for the beam shown below. A comparison between results obtained from the hand solution and spSlab/spBeam software is provided to illustrate the features and capabilities of the spBeam and spSlab software programs. Figure 1 – Continuous Beam Scope 1. Determine the Fixed-End Moments (FEM) ............................................................................................................. 1 2. Slope-Deflection Equations ...................................................................................................................................... 2 2.1. Span AB (End Span with Far End Fixed) ......................................................................................................... 2 2.2. Span BC (End Span with Far End Pinned)........................................................................................................ 2 3. Equilibrium Equations.............................................................................................................................................. 3 4. Shear and Moment Diagrams ................................................................................................................................... 3 5. spSlab/spBeam Software Program Model Solution ................................................................................................. 4 6. Summary and Comparison of Results .....................................................................................................................
    [Show full text]
  • Behavior of Two-Span Continous Reinforced Concrete Beams
    Behavior Of Two-Span Continous Reinforced Concrete Beams A thesis presented to the faculty of the Russ College of Engineering and Technology of Ohio University In partial fulfillment of the requirements for the degree Master of Science Colin Michael McCarty August 2008 2 This thesis titled Behavior of Two-Span Continuous Reinforced Concrete Beams by COLIN MICHAEL MCCARTY has been approved for the Department of Civil Engineering and the Russ College of Engineering and Technology by Eric P. Steinberg Associate Professor of Civil Engineering Dennis Irwin Dean, Russ College of Engineering and Technology 3 ABSTRACT MCCARTY, COLIN MICHAEL, M.S., August 2008, Civil Engineering Behavior of Two-Span Continuous Reinforced Concrete Beams (127 pp.) Director of Thesis: Eric P. Steinberg Nine two-span continuous reinforced concrete beams with point loads 4.5’ on each side of the center reaction were tested. This tested the effects of both longitudinal and transverse reinforcement on the shear strength of the beams. Results of testing the beams showed the effects of maximum shear and moment occurring at the center reaction. Three sets of three beams were constructed to examine the affect varying the amount of longitudinal and transverse reinforcement had on the flexural and shear strength of a beam. Within a set of beams, the longitudinal reinforcement consisted of either No. 3, 8, or 10 reinforcing bars. Transverse reinforcement was consistent for each set of beams. The first set of beams did not have any shear reinforcement and were used to study the effects of longitudinal reinforcement on the shear strength of a beam.
    [Show full text]
  • Direct Stiffness Method: Beams
    Module 4 Analysis of Statically Indeterminate Structures by the Direct Stiffness Method Version 2 CE IIT, Kharagpur Lesson 27 The Direct Stiffness Method: Beams Version 2 CE IIT, Kharagpur Instructional Objectives After reading this chapter the student will be able to 1. Derive member stiffness matrix of a beam element. 2. Assemble member stiffness matrices to obtain the global stiffness matrix for a beam. 3. Write down global load vector for the beam problem. 4. Write the global load-displacement relation for the beam. 27.1 Introduction. In chapter 23, a few problems were solved using stiffness method from fundamentals. The procedure adopted therein is not suitable for computer implementation. In fact the load displacement relation for the entire structure was derived from fundamentals. This procedure runs into trouble when the structure is large and complex. However this can be much simplified provided we follow the procedure adopted for trusses. In the case of truss, the stiffness matrix of the entire truss was obtained by assembling the member stiffness matrices of individual members. In a similar way, one could obtain the global stiffness matrix of a continuous beam from assembling member stiffness matrix of individual beam elements. Towards this end, we break the given beam into a number of beam elements. The stiffness matrix of each individual beam element can be written very easily. For example, consider a continuous beam ABCD as shown in Fig. 27.1a. The given continuous beam is divided into three beam elements as shown in Fig. 27.1b. It is noticed that, in this case, nodes are located at the supports.
    [Show full text]
  • Use of Mathcad in Computing Beam Deflection by Conjugate Beam Method
    Session 2649 Use of Mathcad in Computing Beam Deflection by Conjugate Beam Method Nirmal K. Das Georgia Southern University Abstract The four-year, ABET-accredited Civil Engineering Technology curriculum at Georgia Southern University includes a required, junior-level course in Structural Analysis. One of the topics covered is the conjugate beam method for computing slope and deflection at various points in a beam. The conjugate beam method is a geometric method and it relies only on the principles of statics. The usefulness of this method lies in its simplicity. The students can utilize their already acquired knowledge of shearing force and bending moment to determine a beam’s slope and deflection. An approach to teaching this important method of structural analysis that complements the traditional lecturing through inclusion of a powerful, versatile and user-friendly computational tool, is discussed in this paper. Students will learn how to utilize Mathcad to perform a variety of calculations in a sequence and to verify the accuracy of their manual solutions. A Mathcad program is developed for this purpose and examples to illustrate the computer program are also included in this paper. The integration of Mathcad will enhance students’ problem-solving skills, as it will allow them to focus on analysis while the software performs routine calculations. Thus it will promote learning by discovery, instead of leaving the student in the role of a passive observer. Introduction With the objective of enhanced student learning, adoption of various instructional technology and inclusion of computer-aided problem-solving modules into the curriculum has been a trend for civil engineering and civil engineering technology programs.
    [Show full text]
  • Beam Element
    Finite Element Method: Beam Element Beam Element ¨ Beam: A long thin structure that is subject to the vertical loads. A beam shows more evident bending deformation than the torsion and/or axial deformation. ¨ Bending strain is measured by the lateral deflection and the rotation -> Lateral deflection and rotation determine the number of DoF (Degree of Freedom). Mechanics and Design SNU School of Mechanical and Aerospace Engineering Finite Element Method: Beam Element ¨ Sign convention 1. Positive bending moment: Anti-clock wise rotation. 2. Positive load: y$ -direction. 3. Positive displacement: y$ -direction. Mechanics and Design SNU School of Mechanical and Aerospace Engineering Finite Element Method: Beam Element ¨ The governing equation dVdV -w-wdxdx$ +ˆ -dV dV= 0 =0or or w w= =- - dxdx$ ˆ dMdM VVdxdx$ ˆ+-dM dM= =0 0or or VV == dxdx$ˆ Mechanics and Design SNU School of Mechanical and Aerospace Engineering Finite Element Method: Beam Element ¨ Relation between beam curvature (k ) and bending moment 1 M d 2v$ M k = = or = r EI dx$ 2 EI d 2v$ Curvature for small slope (q = dv$ / dx$ ): k = ¨ dx$ 2 r : Radius of the deflection curve, v$: Lateral displacement function along the y$ -axis direction E : Stiffness coefficient, I : Moment of inertia along the z$-axis Mechanics and Design SNU School of Mechanical and Aerospace Engineering Finite Element Method: Beam Element Solving the equation with M , d 2 F d 2v$I EI = - wax$f dx$ 2 HG dx$ 2 KJ When EI is constant, and force and moment are only applied at nodes, d 4v$ EI = 0 dx$ 4 Mechanics and Design SNU School of Mechanical and Aerospace Engineering Finite Element Method: Beam Element Step 1: To select beam element type Step 2: To select displacement function Assumption of lateral displacement 3 2 vˆ() xˆ = a1 xˆ + a 2 xˆ + a 3 xˆ + a 4 - Complete 3-order displacement function is suitable, because it has four degree of freedom (one lateral displacement and one small rotation at each node) - The function is proper, because it satisfies the fundamental differential equation of a beam.
    [Show full text]