University of Illinois at Urbana-Champaign College of Engineering

Size: px
Start display at page:

Download "University of Illinois at Urbana-Champaign College of Engineering"

Transcription

1 University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the pages must be stapled at all times. PS: You are allowed to use pen, pencil, eraser, and a calculator. Mobile phones are NOT allowed during the examination. Please turn your cell phones off. PS3.: This is a closed book, closed notes, open minds exam. Show A work on the exam sheets SCORE: Problem : 0 / 0 points Problem : 0 / 0 points TOTA : 00 Problem 3: 30 / 30 points Problem 4: 30 / 30 points Page of

2 Problem (0 points total) In each question below, define or describe the term given and explain its significance in relation to finite element analysis, theory, or practice. Use a few sentences, several phrases, or an outline form. You will be graded for content and logical presentation of the topics. Topic (): Boundary versus domain methods (e.g. BEM x FEM) Numerical techniques, such as boundary and domain methods, are approximate methods to solve the governing differential equations of a problem. For a linear problem, the boundary element method (BEM) requires discretization only on the boundary of the domain, e.g. 3D problems are discretized on the D surface (using plane elements), and D problems are discretized along the D contour (using line elements). The most common form of the BEM is based on the collocation method and BEM equations can be derived either using Green s theorem or the method of weighted residuals (MWR). The BEM is based on the so called fundamental solution of the problem, which in the case of elasticity, is a solution of a point load in an infinite medium (Kelvin solution). On the other hand, the finite element method (FEM) is a domain method that requires discretization both on the boundary and interior of the domain. The FEM relies on variational methods, e.g. MWR, from which the most popular is the Galerkin Method. The standard BEM leads to matrices that are full and unsymmetric, while the FEM leads to matrices that are banded and symmetric. For the same problem and similar level of accuracy, BEM matrices tend to be smaller than FEM matrices. The advantages of the BEM in comparison with FEM is that the BEM is especially suited to model infinite domain problems (without the need for infinite elements) and BEM can treat incompressible materials without any difficulty associated with ill-conditioning. However, the BEM is not as general as the FEM in the sense that BEM is not applicable to problems where the fundamental solution is unknown. Moreover, the mathematics of BEM is comparatively more difficult than the FEM. Other domain methods include the finite volume method (FVM) and the finite difference method (FDM). Finite volume refers to the small volume surrounding each node point. Because the flux entering a given volume is identical to that leaving the adjacent volume, these methods are conservative. The FVM is popular in modeling computational fluid dynamics (CFD) problems. The FDM consists of approximating the solutions to differential equations using finite difference equations to approximate derivatives. The domain methods are alternatives to boundary methods. Moreover, it is not necessary to do an analysis based entirely on a single method. There are ways to couple two methods in a single analysis so that each is used to represent portions of the model to which it is best suited. Thus the best properties of each method can be explored in the solution of a problem. Page of

3 Problem (continued) Topic (): Self-Adjoint Operator The method of weighted residuals (MWR) relies on the inner product <R,W>=0 where R is the residual and W is the weight. The residual is associated to a linear operator (u). If we integrate the inner product <(u),w>=0 by parts until all derivatives in u have been eliminated, this will lead to the transposed form on the inner product. As a result, the integration by parts produces domain <*(w), u> and boundary terms S and G which are differential operators. The term S*(w) contains the w terms resulting from the initial phase of the integration, and S(u) contains the corresponding u terms. The operator * is the adjoint of. If *=, then the operator is selfadjoint. In this case, S=S* and G=G*. In addition to determining if the operator is self-adjoint, the integration by parts also generates two types of boundary conditions (BCs): the set S(u) is called essential BCs and the set G(u) is called non-essential or natural BCs. One can specify either type of BC on the boundary of the domain. The essential BCs must be enforced for the solution to be unique (e.g. FEM). We notice that self-adjointness of an operator (in the continuum sense) is similar to symmetry of a matrix (in the discrete sense). In the context of the FEM, the most popular MWR is the Galerkin method, which leads to a symmetric matrix, when the associated operator is self-adjoint. Page 3 of

4 Problem (0 points) Consider the differential equation: u, xx + 4u = in the range 0<x<, with essential boundary conditions u=3 at x=0 and u= at x=. There are no natural boundary conditions. An admissible one-parameter approximating polynomial is u* = 3 - x + a(x -x). (a) Determine the single parameter a by the east Squares Method. (b) Determine the Residual and the Error at x=0.6 and compare the quantities. The exact solution is given by u (exact ) = sin (x) [Acknowledgement: this question is based on Homework #4, Problem 5.-3 of Cook et al. 00] Solution) Part (a) Exact solution at: at x 0.6 u 3.5sin( 0.6) 0.48 Residual methods: exact u 3 xa( x x) ; u, x a(x ) ; u, xx a R 4 ax (4a 8) x a R a 4x 4x east Squares Method: 0 R R dx 0 a [4 ax (4a 8) x a](4x 4x ) dx 0 0 yields: 0 a and 0 34 u x x Thus R x x Part (b) Rx ( 0.6).885 ux ( 0.6).43 Error: u u exact % Error = u u exact exact u % 0.48 Page 4 of

5 Problem 3 (30 points total) - Part I (0 points) Consider a uniformly distributed transverse loading q acting on the beam element shown in the figure below. Determine the consistent nodal load vector using the equivalent nodal load concept of FEM. q The beam shape-functions are (0 < x < ): N = 3x / + x 3 / 3 ; N = x x / + x 3 / N 3 = 3x / x 3 / 3 ; N 4 = x / + x 3 / [ Hint: Draw a FBD to verify that your answer in indeed correct! ] [ Acknowledgment: This question is similar to the one in the Practice Quiz #.] Solution) Equivalent oad oads: N T T T N Wext ue Peq ue N dx u u u3 u4 qdx N 0 3 N 4 N N N4 Peqq dx N x x Ndx x x x x N dx x x N3 dx x x N dx Page 5 of

6 Peq q q q q q Fy q q q q q M () Page 6 of

7 Problem 3 (30 points total) - Part II (0 points) continued Draw the quadratic strain triangle (QST) and determine its corner node shape functions using area coordinates (,, 3 ). Please do a clear drawing. [ Hint: Use the concept of SHAPE FUNCTION MAGIC! Note: This problem was discussed in class and was suggested as an optional HW. ] Solution) Recall the triangle: x x y xy y x xyxy y 3 3 N c(3)(3 ) At Node :, N c c N (3)(3) Ni i(3i )(3i ) ( i,,3) N4 c(3 ) At Node 4:,, N4 c c N4 (3 ) N0 c 3 At Node 0:,, 3 N N5 (3 ) N6 3(3 ) N7 3(33 ) N8 3(33) N 3(3), N0 c c Page 7 of

8 Problem 4 (30 points total) [Patran/ABAQUS] [Acknowledgement: This question is related to your project, and is intended to assess your knowledge of FEM to solve practical engineering problems. It is also a continuation of the last problem of Quiz #]. We want to conduct finite element analysis of the -D plate with a hole, shown in the first figure below, using linear triangular elements, i.e. T3 (you programmed this element in HW#5). The bar has a unit thickness, and is subjected to an axial stress applied over the top and bottom (complete model). A plane-stress idealization is adopted for the -D model. A quarter-symmetric FE model is shown in the second figure below along with the corresponding finite FE mesh. y 3 R= (0,0) x C A B Figure. Geometry of a plate with a hole Figure. Finite element mesh Table : FEM Stress distribution along the A-B cross section. Nodal position Averaged stress values X Y σ x σ y τ xy Page 8 of

9 (a) Do a clear sketch of the quarter-symmetric model and indicate the displacement boundary conditions to be applied to the nodes along the lines Y=0 and X=0. y (0, x (b) Describe the load boundary conditions (*COAD) that should be applied to the nodes along the line Y=. The six element edges at Y= have the same length *COAD,, 0.5,, 0.5 3,, 0.5 4,, 0.5 5,, 0.5 6,, 0.5 7,, 0.5 (c) The averaged nodal stresses across the section A-B are listed in Table. Using those values, check equilibrium along the Y-direction (i.e. quantify how well equilibrium is satisfied). Use the trapezoidal rule to obtain the net force along the Y=0. ( )( ) / ( )( ) / ( )(.83.5) ( )(.5.70) / (.5.5)(.70.) / (.375.5)(. 0.) (.5.375)(0. 0.7) /.5.5 Page of

10 (d) In general, how can the FEM solution for stresses be improved? h -, p -, r -, s - refinements or a combination (e) Estimate quantitatively the quality of the solution at point B. Note that point B is located at a traction-free boundary. (Hint: Recall Cauchy relation from CEE47). n n x x xy x n n 0.43 y y xy xy The solution quality at point B is not very good since it does not satisfy the traction-free boundary (f) From a mechanics point of view (CEE47 is a pre-requisite for this class), what are the stresses σ Y and τ XY at point C? What do you expect the FEM solution at point C to be and why? σ Y and τ XY should be zero because of the traction-free boundary condition at point C. However FEM solution may show non-zero results since results of stress are obtained numerically. (g) Based on your knowledge of CEE47 and using the stress values provided in Table, what are the principal stresses at (X,Y) = (.375,0)? - Mohr s circle, , Stress tensor Eigenvalues of the stress tensor are,.068, Page 0 of

11 (h) The stress concentration factor K for flat bars with circular holes is shown in the below figure. In the figure below, circle the value of the stress concentration factor (K) of the bar with a hole (see Figure ). max P K, nom, t thickness nom ct 3.8 P b c/ d P.6 c/ K d/b Figure Stress concentration factor K for flat bars with circular holes. (Ref. Roark, R.J., and Young, W. C. (8). Formulars for stress and strain, 6th Ed., Mcgraw- Hill Book co., Inc., Newyork.) (i) Calculate the stress concentration factor (K) from the actual finite element analysis results. nom P 3, 4.5 max 4.5, K=.5 ct.5 Page of

University of Illinois at Urbana-Champaign College of Engineering

University of Illinois at Urbana-Champaign College of Engineering University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #2 April 14, 2014 Name: SOLUTION ID#: PS1.:

More information

Contents as of 12/8/2017. Preface. 1. Overview...1

Contents as of 12/8/2017. Preface. 1. Overview...1 Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite

More information

FINAL EXAMINATION. (CE130-2 Mechanics of Materials)

FINAL EXAMINATION. (CE130-2 Mechanics of Materials) UNIVERSITY OF CLIFORNI, ERKELEY FLL SEMESTER 001 FINL EXMINTION (CE130- Mechanics of Materials) Problem 1: (15 points) pinned -bar structure is shown in Figure 1. There is an external force, W = 5000N,

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

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

Lecture 15: Revisiting bars and beams.

Lecture 15: Revisiting bars and beams. 3.10 Potential Energy Approach to Derive Bar Element Equations. We place assumptions on the stresses developed inside the bar. The spatial stress and strain matrices become very sparse. We add (ad-hoc)

More information

COMPUTATIONAL ELASTICITY

COMPUTATIONAL ELASTICITY COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY

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

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST

More information

ME 1401 FINITE ELEMENT ANALYSIS UNIT I PART -A. 2. Why polynomial type of interpolation functions is mostly used in FEM?

ME 1401 FINITE ELEMENT ANALYSIS UNIT I PART -A. 2. Why polynomial type of interpolation functions is mostly used in FEM? SHRI ANGALAMMAN COLLEGE OF ENGINEERING AND TECHNOLOGY (An ISO 9001:2008 Certified Institution) SIRUGANOOR, TIRUCHIRAPPALLI 621 105 Department of Mechanical Engineering ME 1401 FINITE ELEMENT ANALYSIS 1.

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

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

Discontinuous Distributions in Mechanics of Materials

Discontinuous Distributions in Mechanics of Materials Discontinuous Distributions in Mechanics of Materials J.E. Akin, Rice University 1. Introduction The study of the mechanics of materials continues to change slowly. The student needs to learn about software

More information

ME FINITE ELEMENT ANALYSIS FORMULAS

ME FINITE ELEMENT ANALYSIS FORMULAS ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness

More information

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric

More information

CRITERIA FOR SELECTION OF FEM MODELS.

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

More information

Back Matter Index The McGraw Hill Companies, 2004

Back Matter Index The McGraw Hill Companies, 2004 INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See

More information

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

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

More information

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

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

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

More information

Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems

Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.

More information

GATE SOLUTIONS E N G I N E E R I N G

GATE SOLUTIONS E N G I N E E R I N G GATE SOLUTIONS C I V I L E N G I N E E R I N G From (1987-018) Office : F-16, (Lower Basement), Katwaria Sarai, New Delhi-110016 Phone : 011-65064 Mobile : 81309090, 9711853908 E-mail: info@iesmasterpublications.com,

More information

ME C85/CE C30 Fall, Introduction to Solid Mechanics ME C85/CE C30. Final Exam. Fall, 2013

ME C85/CE C30 Fall, Introduction to Solid Mechanics ME C85/CE C30. Final Exam. Fall, 2013 Introduction to Solid Mechanics ME C85/CE C30 Fall, 2013 1. Leave an empty seat between you and the person (people) next to you. Unfortunately, there have been reports of cheating on the midterms, so we

More information

Nonconservative Loading: Overview

Nonconservative Loading: Overview 35 Nonconservative Loading: Overview 35 Chapter 35: NONCONSERVATIVE LOADING: OVERVIEW TABLE OF CONTENTS Page 35. Introduction..................... 35 3 35.2 Sources...................... 35 3 35.3 Three

More information

Chapter 3: Stress and Equilibrium of Deformable Bodies

Chapter 3: Stress and Equilibrium of Deformable Bodies Ch3-Stress-Equilibrium Page 1 Chapter 3: Stress and Equilibrium of Deformable Bodies When structures / deformable bodies are acted upon by loads, they build up internal forces (stresses) within them to

More information

Code No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)

More information

Math 1132 Practice Exam 1 Spring 2016

Math 1132 Practice Exam 1 Spring 2016 University of Connecticut Department of Mathematics Math 32 Practice Exam Spring 206 Name: Instructor Name: TA Name: Section: Discussion Section: Read This First! Please read each question carefully. Show

More information

Introduction to Finite Element Method

Introduction to Finite Element Method Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,

More information

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is

More information

IV B.Tech. I Semester Supplementary Examinations, February/March FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80

IV B.Tech. I Semester Supplementary Examinations, February/March FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 www..com www..com Code No: M0322/R07 Set No. 1 IV B.Tech. I Semester Supplementary Examinations, February/March - 2011 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer

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

A short review of continuum mechanics

A short review of continuum mechanics A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material

More information

Name (Print) ME Mechanics of Materials Exam # 1 Date: October 5, 2016 Time: 8:00 10:00 PM

Name (Print) ME Mechanics of Materials Exam # 1 Date: October 5, 2016 Time: 8:00 10:00 PM Name (Print) (Last) (First) Instructions: ME 323 - Mechanics of Materials Exam # 1 Date: October 5, 2016 Time: 8:00 10:00 PM Circle your lecturer s name and your class meeting time. Gonzalez Krousgrill

More information

Computation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads

Computation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 215 Paper o. 31 Computation Time Assessment of a Galerkin Finite Volume

More information

Failure from static loading

Failure from static loading Failure from static loading Topics Quiz /1/07 Failures from static loading Reading Chapter 5 Homework HW 3 due /1 HW 4 due /8 What is Failure? Failure any change in a machine part which makes it unable

More information

Due Tuesday, September 21 st, 12:00 midnight

Due Tuesday, September 21 st, 12:00 midnight Due Tuesday, September 21 st, 12:00 midnight The first problem discusses a plane truss with inclined supports. You will need to modify the MatLab software from homework 1. The next 4 problems consider

More information

Interpolation Functions for General Element Formulation

Interpolation Functions for General Element Formulation CHPTER 6 Interpolation Functions 6.1 INTRODUCTION The structural elements introduced in the previous chapters were formulated on the basis of known principles from elementary strength of materials theory.

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

Exam paper: Biomechanics

Exam paper: Biomechanics Exam paper: Biomechanics Tuesday August 10th 2010; 9.00-12.00 AM Code: 8W020 Biomedical Engineering Department, Eindhoven University of Technology The exam comprises 10 problems. Every problem has a maximum

More information

ELASTICITY AND FRACTURE MECHANICS. Vijay G. Ukadgaonker

ELASTICITY AND FRACTURE MECHANICS. Vijay G. Ukadgaonker THEORY OF ELASTICITY AND FRACTURE MECHANICS y x Vijay G. Ukadgaonker Theory of Elasticity and Fracture Mechanics VIJAY G. UKADGAONKER Former Professor Indian Institute of Technology Bombay Delhi-110092

More information

Curriculum Catalog

Curriculum Catalog 2017-2018 Curriculum Catalog 2017 Glynlyon, Inc. Table of Contents INTEGRATED MATH I COURSE OVERVIEW... 1 UNIT 1: FOUNDATIONS OF ALGEBRA... 1 UNIT 2: THE LANGUAGE OF ALGEBRA... 2 UNIT 3: GEOMETRY... 2

More information

Q1 Q2 Q3 Q4 Tot Letr Xtra

Q1 Q2 Q3 Q4 Tot Letr Xtra Mathematics 54.1 Final Exam, 12 May 2011 180 minutes, 90 points NAME: ID: GSI: INSTRUCTIONS: You must justify your answers, except when told otherwise. All the work for a question should be on the respective

More information

CE 715: Advanced Strength of Materials

CE 715: Advanced Strength of Materials CE 715: Advanced Strength of Materials Lecture 1 CE 715 Course Information Instructor: Tasnim Hassan Office: Mann Hall 419 Office Hours: TTh 2:00-4:00 pm Phone: 515-8123 Email: thassan@eos.ncsu.edu 1 Course

More information

Department of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Spring Semester 2019

Department of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Spring Semester 2019 Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Spring Semester 2019 1. Please write your 1- or 2-digit exam number

More information

ENGR-1100 Introduction to Engineering Analysis. Lecture 17

ENGR-1100 Introduction to Engineering Analysis. Lecture 17 ENGR-1100 Introduction to Engineering Analysis Lecture 17 CENTROID OF COMPOSITE AREAS Today s Objective : Students will: a) Understand the concept of centroid. b) Be able to determine the location of the

More information

CIVL4332 L1 Introduction to Finite Element Method

CIVL4332 L1 Introduction to Finite Element Method CIVL L Introduction to Finite Element Method CIVL L Introduction to Finite Element Method by Joe Gattas, Faris Albermani Introduction The FEM is a numerical technique for solving physical problems such

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

Lecture No. 5. For all weighted residual methods. For all (Bubnov) Galerkin methods. Summary of Conventional Galerkin Method

Lecture No. 5. For all weighted residual methods. For all (Bubnov) Galerkin methods. Summary of Conventional Galerkin Method Lecture No. 5 LL(uu) pp(xx) = 0 in ΩΩ SS EE (uu) = gg EE on ΓΓ EE SS NN (uu) = gg NN on ΓΓ NN For all weighted residual methods NN uu aaaaaa = uu BB + αα ii φφ ii For all (Bubnov) Galerkin methods ii=1

More information

If the solution does not follow a logical thought process, it will be assumed in error.

If the solution does not follow a logical thought process, it will be assumed in error. Please indicate your group number (If applicable) Circle Your Instructor s Name and Section: MWF 8:30-9:20 AM Prof. Kai Ming Li MWF 2:30-3:20 PM Prof. Fabio Semperlotti MWF 9:30-10:20 AM Prof. Jim Jones

More information

Practice Final Examination. Please initial the statement below to show that you have read it

Practice Final Examination. Please initial the statement below to show that you have read it EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use

More information

L e c t u r e. D r. S a s s a n M o h a s s e b

L e c t u r e. D r. S a s s a n M o h a s s e b The Scaled smteam@gmx.ch Boundary Finite www.erdbebenschutz.ch Element Method Lecture A L e c t u r e A1 D r. S a s s a n M o h a s s e b V i s i t i n g P r o f e s s o r M. I. T. C a m b r i d g e December

More information

Physical Science and Engineering. Course Information. Course Number: ME 100

Physical Science and Engineering. Course Information. Course Number: ME 100 Physical Science and Engineering Course Number: ME 100 Course Title: Course Information Basic Principles of Mechanics Academic Semester: Fall Academic Year: 2016-2017 Semester Start Date: 8/21/2016 Semester

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

EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE OF A RECTANGULAR ELASTIC BODY MADE OF FGM

EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE OF A RECTANGULAR ELASTIC BODY MADE OF FGM Proceedings of the International Conference on Mechanical Engineering 2007 (ICME2007) 29-31 December 2007, Dhaka, Bangladesh ICME2007-AM-76 EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE

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

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

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

More information

Stress, Strain, Mohr s Circle

Stress, Strain, Mohr s Circle Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected

More information

Module 3 : Equilibrium of rods and plates Lecture 15 : Torsion of rods. The Lecture Contains: Torsion of Rods. Torsional Energy

Module 3 : Equilibrium of rods and plates Lecture 15 : Torsion of rods. The Lecture Contains: Torsion of Rods. Torsional Energy The Lecture Contains: Torsion of Rods Torsional Energy This lecture is adopted from the following book 1. Theory of Elasticity, 3 rd edition by Landau and Lifshitz. Course of Theoretical Physics, vol-7

More information

Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE

Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE 1 Chapter 3 Load and Stress Analysis 2 Chapter Outline Equilibrium & Free-Body Diagrams Shear Force and Bending Moments in Beams Singularity Functions Stress Cartesian Stress Components Mohr s Circle for

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

Lecture 12: Finite Elements

Lecture 12: Finite Elements Materials Science & Metallurgy Part III Course M6 Computation of Phase Diagrams H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of continuous quantities such as temperature

More information

Version: A. Earth s gravitational field g = 9.81 N/kg Mass of a Proton m p = kg

Version: A. Earth s gravitational field g = 9.81 N/kg Mass of a Proton m p = kg PHYS 2212 G & J Quiz and Exam Formulæ & Constants Fall 2017 Fundamental Charge e = 1.602 10 19 C Mass of an Electron m e = 9.109 10 31 kg Earth s gravitational field g = 9.81 N/kg Mass of a Proton m p

More information

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

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

More information

ME Final Exam. PROBLEM NO. 4 Part A (2 points max.) M (x) y. z (neutral axis) beam cross-sec+on. 20 kip ft. 0.2 ft. 10 ft. 0.1 ft.

ME Final Exam. PROBLEM NO. 4 Part A (2 points max.) M (x) y. z (neutral axis) beam cross-sec+on. 20 kip ft. 0.2 ft. 10 ft. 0.1 ft. ME 323 - Final Exam Name December 15, 2015 Instructor (circle) PROEM NO. 4 Part A (2 points max.) Krousgrill 11:30AM-12:20PM Ghosh 2:30-3:20PM Gonzalez 12:30-1:20PM Zhao 4:30-5:20PM M (x) y 20 kip ft 0.2

More information

Preprocessor Geometry Properties )Nodes, Elements(, Material Properties Boundary Conditions(displacements, Forces )

Preprocessor Geometry Properties )Nodes, Elements(, Material Properties Boundary Conditions(displacements, Forces ) در برنامه يك تدوين براي بعدي دو يك سازه محيط MATLAB Preprocessor Geometry Properties )Nodes, Elements(, Material Properties Boundary Conditions(displacements, Forces ) Definition of Stiffness Matrices

More information

ME325 EXAM I (Sample)

ME325 EXAM I (Sample) ME35 EXAM I (Sample) NAME: NOTE: COSED BOOK, COSED NOTES. ONY A SINGE 8.5x" ORMUA SHEET IS AOWED. ADDITIONA INORMATION IS AVAIABE ON THE AST PAGE O THIS EXAM. DO YOUR WORK ON THE EXAM ONY (NO SCRATCH PAPER

More information

Mechanics of Earthquakes and Faulting

Mechanics of Earthquakes and Faulting Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body

More information

THE UPPER BOUND PROPERTY FOR SOLID MECHANICS OF THE LINEARLY CONFORMING RADIAL POINT INTERPOLATION METHOD (LC-RPIM)

THE UPPER BOUND PROPERTY FOR SOLID MECHANICS OF THE LINEARLY CONFORMING RADIAL POINT INTERPOLATION METHOD (LC-RPIM) International Journal of Computational Methods Vol. 4, No. 3 (2007) 521 541 c World Scientific Publishing Company THE UPPER BOUND PROPERTY FOR SOLID MECHANICS OF THE LINEARLY CONFORMING RADIAL POINT INTERPOLATION

More information

Lecture 27: Structural Dynamics - Beams.

Lecture 27: Structural Dynamics - Beams. Chapter #16: Structural Dynamics and Time Dependent Heat Transfer. Lectures #1-6 have discussed only steady systems. There has been no time dependence in any problems. We will investigate beam dynamics

More information

Final Exam Ship Structures Page 1 MEMORIAL UNIVERSITY OF NEWFOUNDLAND. Engineering Ship Structures

Final Exam Ship Structures Page 1 MEMORIAL UNIVERSITY OF NEWFOUNDLAND. Engineering Ship Structures Final Exam - 53 - Ship Structures - 16 Page 1 MEMORIA UNIVERSITY OF NEWFOUNDAND Faculty of Engineering and Applied Science Engineering 53 - Ship Structures FINA EXAMINATION SONS Date: Wednesday April 13,

More information

MEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy

MEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy MEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy e-mail: jnreddy@tamu.edu Texas A&M University College Station, TX 77843-3123 USA

More information

Contents. Prologue Introduction. Classical Approximation... 19

Contents. Prologue Introduction. Classical Approximation... 19 Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................

More information

I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. NAME: ME 270 Fall 2012 Examination No. 3 - Makeup Please review the following statement: Group No.: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

More information

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Circuits & Electronics Spring 2007.

Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Circuits & Electronics Spring 2007. Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.002 Circuits & Electronics Spring 2007 Quiz #2 25 April 2007 Name: There are 20 pages in this quiz, including

More information

Mohr's Circle for 2-D Stress Analysis

Mohr's Circle for 2-D Stress Analysis Mohr's Circle for 2-D Stress Analysis If you want to know the principal stresses and maximum shear stresses, you can simply make it through 2-D or 3-D Mohr's cirlcles! You can know about the theory of

More information

Practice Questions Exam 1/page1. PES Physics 2 Practice Exam 1 Questions. Name: Score: /.

Practice Questions Exam 1/page1. PES Physics 2 Practice Exam 1 Questions. Name: Score: /. Practice Questions Exam 1/page1 PES 110 - Physics Practice Exam 1 Questions Name: Score: /. Instructions Time allowed for this is exam is 1 hour 15 minutes 5 multiple choice (5 points) 3 to 5 written problems

More information

Form #221 Page 1 of 7

Form #221 Page 1 of 7 Version Quiz #2 Form #221 Name: A Physics 2212 GH Spring 2016 Recitation Section: Print your name, quiz form number (3 digits at the top of this form), and student number (9 digit Georgia Tech ID number)

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

Mathematics Qualifying Exam Study Material

Mathematics Qualifying Exam Study Material Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors

More information

Math 164-1: Optimization Instructor: Alpár R. Mészáros

Math 164-1: Optimization Instructor: Alpár R. Mészáros Math 164-1: Optimization Instructor: Alpár R. Mészáros First Midterm, April 20, 2016 Name (use a pen): Student ID (use a pen): Signature (use a pen): Rules: Duration of the exam: 50 minutes. By writing

More information

Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis

Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods

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

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

Introduction to Heat and Mass Transfer. Week 7

Introduction to Heat and Mass Transfer. Week 7 Introduction to Heat and Mass Transfer Week 7 Example Solution Technique Using either finite difference method or finite volume method, we end up with a set of simultaneous algebraic equations in terms

More information

McGill University April 20, Advanced Calculus for Engineers

McGill University April 20, Advanced Calculus for Engineers McGill University April 0, 016 Faculty of Science Final examination Advanced Calculus for Engineers Math 64 April 0, 016 Time: PM-5PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer Student

More information

ME185 Introduction to Continuum Mechanics

ME185 Introduction to Continuum Mechanics Fall, 0 ME85 Introduction to Continuum Mechanics The attached pages contain four previous midterm exams for this course. Each midterm consists of two pages. As you may notice, many of the problems are

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

IDE 110 Mechanics of Materials Spring 2006 Final Examination FOR GRADING ONLY

IDE 110 Mechanics of Materials Spring 2006 Final Examination FOR GRADING ONLY Spring 2006 Final Examination STUDENT S NAME (please print) STUDENT S SIGNATURE STUDENT NUMBER IDE 110 CLASS SECTION INSTRUCTOR S NAME Do not turn this page until instructed to start. Write your name on

More information

46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference April 2005 Austin, Texas

46th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference April 2005 Austin, Texas th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference - April, Austin, Texas AIAA - AIAA - Bi-stable Cylindrical Space Frames H Ye and S Pellegrino University of Cambridge, Cambridge,

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

Mechanics PhD Preliminary Spring 2017

Mechanics PhD Preliminary Spring 2017 Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n

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

CONSTITUTIVE RELATIONS FOR LINEAR ELASTIC SOLIDS

CONSTITUTIVE RELATIONS FOR LINEAR ELASTIC SOLIDS Chapter 9 CONSTITUTIV RLATIONS FOR LINAR LASTIC SOLIDS Figure 9.1: Hooke memorial window, St. Helen s, Bishopsgate, City of London 211 212 CHAPTR 9. CONSTITUTIV RLATIONS FOR LINAR LASTIC SOLIDS 9.1 Mechanical

More information

Introduction to Engineering Analysis - ENGR1100 Course Description and Syllabus Tuesday / Friday Sections. Spring '13.

Introduction to Engineering Analysis - ENGR1100 Course Description and Syllabus Tuesday / Friday Sections. Spring '13. Introduction to Engineering Analysis - ENGR1100 Course Description and Syllabus Tuesday / Friday Sections Spring 2013 Back exams, HW solutions, and other useful links can be found at the following website:

More information

Solution Methods. Steady State Diffusion Equation. Lecture 04

Solution Methods. Steady State Diffusion Equation. Lecture 04 Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.

More information

BOOK OF COURSE WORKS ON STRENGTH OF MATERIALS FOR THE 2 ND YEAR STUDENTS OF THE UACEG

BOOK OF COURSE WORKS ON STRENGTH OF MATERIALS FOR THE 2 ND YEAR STUDENTS OF THE UACEG BOOK OF COURSE WORKS ON STRENGTH OF MATERIALS FOR THE ND YEAR STUDENTS OF THE UACEG Assoc.Prof. Dr. Svetlana Lilkova-Markova, Chief. Assist. Prof. Dimitar Lolov Sofia, 011 STRENGTH OF MATERIALS GENERAL

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

Abstract. keywords: finite elements, beams, Hermitian interpolation, shape functions

Abstract. keywords: finite elements, beams, Hermitian interpolation, shape functions Generation of shape functions for straight beam elements Charles E. Augarde Originally published in Computers and Structures, 68 (1998) 555-560 Corrections highlighted. Abstract Straight beam finite elements

More information

Multi-Point Constraints

Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Single point constraint examples Multi-Point constraint examples linear, homogeneous linear, non-homogeneous linear, homogeneous

More information

Advanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian

Advanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz

More information