SME 3033 FINITE ELEMENT METHOD. Bending of Prismatic Beams (Initial notes designed by Dr. Nazri Kamsah)

Similar documents
Trusses. Introduction. Plane Trusses Local and global coordinate systems. Plane Trusses Local and global coordinate systems

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

AS 5850 Finite Element Analysis

Keywords: Sectional forces, bending moment, shear force, finite element, Bernoulli-Euler beam, concentrated moving loads

VSMN30 FINITA ELEMENTMETODEN - DUGGA

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

Introduction to Condensed Matter Physics

Elements of Statistical Thermodynamics

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

Dynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA *

Nonlinear Bending of Strait Beams

Finite element discretization of Laplace and Poisson equations

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

u 3 = u 3 (x 1, x 2, x 3 )

FINITE BEAM ELEMENT WITH PIEZOELECTRIC LAYERS AND FUNCTIONALLY GRADED MATERIAL OF CORE

Section 11.6: Directional Derivatives and the Gradient Vector

Active Displacement Feedback Control of a Smart Beam: Analytic and Numerical Solutions

Addition of angular momentum

682 CHAPTER 11 Columns. Columns with Other Support Conditions

CHAPTER 1. Introductory Concepts Elements of Vector Analysis Newton s Laws Units The basis of Newtonian Mechanics D Alembert s Principle

Lecture Outline. Skin Depth Power Flow 8/7/2018. EE 4347 Applied Electromagnetics. Topic 3e

Division of Mechanics Lund University MULTIBODY DYNAMICS. Examination Name (write in block letters):.

PHYSICS 489/1489 LECTURE 7: QUANTUM ELECTRODYNAMICS

Hydrogen Atom and One Electron Ions

The sum-product algorithm on simple graphs

Addition of angular momentum

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers:

4.2 Design of Sections for Flexure

the output is Thus, the output lags in phase by θ( ωo) radians Rewriting the above equation we get

Calculus concepts derivatives

Intro to QM due: February 8, 2019 Problem Set 12

Why is a E&M nature of light not sufficient to explain experiments?

On the Hamiltonian of a Multi-Electron Atom

Constants and Conversions:

Sundials and Linear Algebra

Zero Point Energy: Thermodynamic Equilibrium and Planck Radiation Law

Integration by Parts

1 Isoparametric Concept

Direct Approach for Discrete Systems One-Dimensional Elements

Chapter 5. Introduction. Introduction. Introduction. Finite Element Modelling. Finite Element Modelling

Method for including restrained warping in traditional frame analyses

Abstract Interpretation: concrete and abstract semantics

Einstein Equations for Tetrad Fields

Combinatorial Networks Week 1, March 11-12

Derivation of Electron-Electron Interaction Terms in the Multi-Electron Hamiltonian

Potential energy of a structure. Vforce. joints j

Instantaneous Cutting Force Model in High-Speed Milling Process with Gyroscopic Effect

Dynamic analysis of a Timoshenko beam subjected to moving concentrated forces using the finite element method

Lecture 28 Title: Diatomic Molecule : Vibrational and Rotational spectra

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

Numerical Analysis of Transient Responses for Elastic Structures Connected to a Viscoelastic Shock Absorber Using FEM with a Nonlinear Complex Spring

That is, we start with a general matrix: And end with a simpler matrix:

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Massachusetts Institute of Technology Department of Mechanical Engineering

Engineering 323 Beautiful HW #13 Page 1 of 6 Brown Problem 5-12

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

MATH 319, WEEK 15: The Fundamental Matrix, Non-Homogeneous Systems of Differential Equations

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

The Finite Element Method

2.3 Matrix Formulation

The Transmission Line Wave Equation

The Matrix Exponential

dt d Chapter 30: 1-Faraday s Law of induction (induced EMF) Chapter 30: 1-Faraday s Law of induction (induced Electromotive Force)

PHYS-333: Problem set #2 Solutions

Chemistry 342 Spring, The Hydrogen Atom.

The Matrix Exponential

Differentiation of Exponential Functions

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.

Analysis of Structural Vibration using the Finite Element Method

Math 34A. Final Review

The graph of y = x (or y = ) consists of two branches, As x 0, y + ; as x 0, y +. x = 0 is the

MCE503: Modeling and Simulation of Mechatronic Systems Discussion on Bond Graph Sign Conventions for Electrical Systems

MEEN 617 Handout #12 The FEM in Vibrations A brief introduction to the finite element method for modeling of mechanical structures

Keywords- Active vibration control, cantilever composite beam, Newmark-β method

Background: We have discussed the PIB, HO, and the energy of the RR model. In this chapter, the H-atom, and atomic orbitals.

5. B To determine all the holes and asymptotes of the equation: y = bdc dced f gbd

Dynamic response of a finite length euler-bernoulli beam on linear and nonlinear viscoelastic foundations to a concentrated moving force

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

AerE 344: Undergraduate Aerodynamics and Propulsion Laboratory. Lab Instructions

Voltage, Current, Power, Series Resistance, Parallel Resistance, and Diodes

Basic Polyhedral theory

Chemical Physics II. More Stat. Thermo Kinetics Protein Folding...

Finite Strain Elastic-Viscoplastic Model

San José State University Aerospace Engineering AE 138 Vector-Based Dynamics for Aerospace Applications, Fall 2016

ELECTRON-MUON SCATTERING

Application of Vague Soft Sets in students evaluation

are given in the table below. t (hours)

VII. Quantum Entanglement

Properties of Quarks ( ) Isospin. π = 1, 1

Classical Magnetic Dipole

Exercise 1. Sketch the graph of the following function. (x 2

Instructions for Section 1

Intro to Nuclear and Particle Physics (5110)

5.80 Small-Molecule Spectroscopy and Dynamics

16. Electromagnetics and vector elements (draft, under construction)

INC 693, 481 Dynamics System and Modelling: Linear Graph Modeling II Dr.-Ing. Sudchai Boonto Assistant Professor

The Fourier Transform Solution for the Green s Function of Monoenergetic Neutron Transport Theory

GEOMETRICAL PHENOMENA IN THE PHYSICS OF SUBATOMIC PARTICLES. Eduard N. Klenov* Rostov-on-Don, Russia

Transcription:

Bnding of Prismatic Bams (Initia nots dsignd by Dr. Nazri Kamsah)

St I-bams usd in a roof construction.

5- Gnra Loading Conditions For our anaysis, w wi considr thr typs of oading, as iustratd bow. Not: Positiv oads ar shown. Ths ar: a) Distributd oad, p, ovr a ngth ; b) Concntratd oad, P m, at point m; c) Concntratd momnt, M k, at point k. Goba coordinat systm.

5- Dformation of th Nutra Axis Du to th appid oads, th bam xprincs an in-pan (x-y pan) dformation. Th dformd shap of th bam s nutra axis is shown. Not: Positiv dformation is shown. Dformation at any point is charactrizd by two paramtrs: a) Vrtica dfction, v b) Rotation or Sop, v

5-5 Emntary Bam Euations For sma dformations, mntary bam thory stats, strss, M I y strain, E My EI curvatur, d dx v M EI Not: Th bam s cross-sction is considrd symmtric with rspct to th pan of oading (x-y pan).

5- Finit Emnt Formuation Th bam is subdividd into four sctions or mnts, connctd by nods. Each nod has dgrs of frdom. At nod i, thy ar dispacmnt componnts: Q i- (vrtica dfction) and Q i (sop or rotation). Positiv dformations x Goba coordinat systm Loca coordinat systm ' ' v v v v T

5-7 Hrmit Shap Function W dfin Hrmit shap functions, so as to obtain th dformation at any point aong a sing bam mnt. It is givn in a gnra form, H i a i b c i i d i i,,, Not: a, b, c, and d ar cofficints.

Th shap functions ar xprssd in trms of, as H H H H H H H H = - = or or or or H H H H Vaus of th H i at oca nods ( = -) and ( = +) ar,

5-8 Dispacmnt Function Th dformation at any point on a bam mnt can b xprssd in trms of th Hrmit shap functions as, dv dv v Hv H Hv H d d W nd to rpac th brackt trms. Using chain ru of diffrntiation, Rca, dv d dv dx x x x x dx d dv dx x x x x x dx dξ dv Not: Th trm at oca nod is v and at nod is v. dx

Substituting th brackt trm, v() can now b xprssd as, v H v H v H v H v In a condnsd matrix form, ' ' v H whr H H H H H v ' v v ' v

5-9 Th Potntia Enrgy Approach W now dvop finit mnt formuation for a sing bam mnt using th potntia nrgy approach. Th tota potntia nrgy, p p of th ntir bam is givn by p p L EI d v dx dx L pv dx m P m v m k M k v ' k Intrna strain nrgy Potntia nrgy du to distributd oad, concntratd forc and concntratd momnt Not: p is th distributd oad; P m is point oad at point m; M k is th momnt appid at point k; v m is dfction at point m; v k is sop at point k.

5- Emnt Stiffnss Matrix Th stiffnss matrix [k] for a sing bam mnt is now drivd using th potntia nrgy approach. Th intrna strain nrgy of a sing bam mnt, U L d v EI dx dx with dx d..(i) W nd to xprss (d v/dx ) in trms of H, and. Rca, dv dv d v d v dx d dx d Substitut v = [H]{} and simpifying, yids T T d v d H d H dx d d Assignmnt: Driv th trm U =f() from n (i)..(ii)

Aso, d H d..(iii) Substitut Es.(ii), and (iii) into E.(i) and substituting dx = ( /)d, w obtain, 9 9 8 8 9 T 8EI 8 U d 9 8 Symmtric

Intgrating ach trm in th matrix and noting that, d ; d ; d ; th intrna strain nrgy can now b writtn in condnsd matrix form as, U T k whr [k] rprsnts th mnt stiffnss matrix, givn by k EI Connctivity with oca noda DOFs

Exrcis 5- An ovrhang bam carris transvrs oads as shown. Cross-sction: b = h = 9. mm; A = 8.85 mm Mod th bam using two mnts. Us E = GPa, and I = x mm. a) Writ th stiffnss matrix, [k] for ach mnt; b) Assmb th goba stiffnss matrix, [K] for th ntir bam.

5- Emnt Load Vctor Th distributd oad, p acts aong th ngth of an mnt. It has to b transformd into uivant concntratd forcs, acting at th nods. Th potntia nrgy du to th distributd oad p is L p pv dx Substitut for v H p H d and dx p d w obtain, Substituting for [H] and rarranging yids, p p H H H H d

Expanding H i, intgrating vry trms with rspct to w gt, p p p p p which can b writtn in a condnsd matrix form as, p f T whr {f} T is th uivant noda forc vctor du to th distributd oad p, givn by f p, p, p, p T

Physica Intrprtation Actua oading condition p Euivant noda forcs p p p p Vctor of th uivant noda forcs f p p p p T

Exrcis 5- Rconsidr th ovrhang bam in Exrcis 5-. Cross-sction: b = h = 9. mm; A = 8.85 mm Us E = GPa, I = x mm. For th bam undr th oading shown, a) Writ th vctor of uivant noda forcs for sction AB; b) Assmb th goba oad vctor for th ntir bam.

For a sing bam mnt with distributd oad p, th systm of inar uation can b xprssd as, 5- Systm of Linar Euations Expanding, w hav f k p p p p ' ' v v v v EI Not: Connctivity with oca noda dgrs of frdom is shown.

Exrcis 5- Rconsidr th ovrhang bam in Exrcis 5-. Cross-sction: b = h = 9. mm; A = 8.85 mm Us E = GPa, I = x mm. a) Assmb th goba systm of inar uations (SLEs); b) Impos th boundary conditions, and writ th rducd SLEs; c) Sov for th unknown noda dgrs of frdom; d) Find th vrtica dfction, v at th midd of sction AB.

Bnding Momnt 5- Shar Forc and Bnding Momnt Bnding momnt M at any point aong th bam mnt is, It was shown that, with dx v d EI M 8 8 dx v d d v d dx v d (i) (ii) v H

EI M Substituting.(ii) into.(i) and simpifying yids, Not: Rca that, ' ' v v v v

Shar Forc Shar forc V at any point aong th bam mnt, V dm dx EI It can b shown that, with d v dx 8 d dx d v d v v H (iii) d v dx 8 8 8 (iv) Substituting.(iv) into.(iii) and simpifying yids, V EI

Exrcis 5- Rconsidr th ovrhang bam in Exrcis -. Cross-sction: b = h = 9. mm; A = 8.85 mm Us E = GPa, I = x mm. a) Estimat th bnding momnt and shar forc at th midpoint of sction AB. b) Comput th raction forcs at support B.

Examp 5- For a bam with th oading shown bow, dtrmin: (a) th sops at nods and, and (b) th vrtica dfction at midpoint of sction -. Euivant noda forcs for sction -.

Soution. Stiffnss matrix for ach mnt, k 8 Q Q Q Q Q Q Q 5 Q k 8

8 8 5 K Nm N, f. Assmb goba stiffnss matrix for ntir bam,. Euivant noda forc vctor du to distributd oad p, on mnt Q Q Q Q Q 5 Q m N

F. Assmb goba oad vctor, 8 8 5 5 Q Q Q Q Q Q 5. Writ th goba systm of inar uations,

8 8 5 Q Q..79 Q Q. Impos th boundary conditions. W hav, 7. Soving for th unknown dgrs of frdom, w obtain, 5 Q Q Q Q Using th imination mthod, th goba SLEs rducd to, which ar th sops at nods and rspctivy. radians 8. Th vrtica dfction at any point on th mnt is givn by, H H H H v

At th midpoint of th mnt, =. So that, H H H H m 5 8 9 79. v.. H H H H v Aso, for mnt, =, = Q = -.79-, =, = Q =.-. Substituting ths into v(), w gt

- Pan Fram Structur Pan structurs is mad up of rigidy connctd mmbrs. Each mmbr can b tratd as a bam. Axia oads xist, thus axia dformations occur in fram structurs. Emnts of th fram structur can b in diffrnt orintation.

Considr a sing fram mnt. Not: Each oca nod has thr dgrs of frdom, i.. two dispacmnt componnts and on rotationa componnt. Noda dispacmnt vctor in goba coordinat systm is givn by: T 5 Noda dispacmnt vctor in oca coordinat systm is givn by: ' ' ' ' ' ' T ' 5

Not that th rotationa componnts, = and =, Th oca DOFs transform into th goba DOFs according to a ration whr [L] is th transformation matrix givn by, L ' m m m m L Not: and m ar th dirction cosins for th fram mnt.

Not:,, 5, and ar th bam dgrs of frdom, whi and ar simiar to th dispacmnts of a -D mnt. Combining th stiffnss matrics for a bam and -D mnts givs stiffnss matrix for a fram mnt, i.. k ' 5 EA EA EI EI EI EI EI EI EI EI EA EA EI EI EI EI EI EI EI EI

As shown in th drivation of a truss mnt, th mnt strain nrgy is givn by: T U ' k ' ' T T U L k ' L Thrfor th mnt stiffnss matrix for a sing fram mnt, in goba coordinats is: k L T k ' L

Emnt Forc Vctor Suppos a uniformy distributd oad p acts on a fram mnt. W nd to trans-form p into uiva-nt noda forcs. It can b shown that: ' T f ' T L T f ' whr f ' p p p p T (in oca coordinat) Rca, f T L f '

Examp - For a porta fram structur shown, dtrmin th dispacmnts and rotations at joints and.

Soution Stp : Emnt connctivity Emnt No st Nod nd Nod Stp : Emnt stiffnss matrics Emnt k k ' Q Q Q Q Q 5 Q. 7. 7. 78 5.. 78 5. 5. 57 5. 78. 7. 7. 78. 78 5. 5. 5. 78 5. 57

Emnt & (Loca coordinat) For mnt k ' k ' For mnt Q 7 Q 8 Q 9 Q Q Q.5.5.5 7.5 7 7 85 7.5.5.5 7.5 7 7 7 85 Q Q Q Q Q 5 Q Th transformation matrix [L] for mnt and, whr = and m =, is as foows: T L k L k ' L

Emnt & (Goba coordinat) For mnt Q 7 Q 8 Q 9 Q Q Q.5 7.5 7.5.5 7 85 7.5 7.5 7.5.5 7 7 85 For mnt Q Q Q Q Q 5 Q k k Stp : Assmb th goba stiffnss matrix K. 7.7. 5..78 5. 7 5. 5 5. 78.7. 7.78 5.. 5. 5. 78 7 5. 5

Stp : Estabish th oad vctor F 7 7 Stp 5: Writ th systm of inar uations K Q F. 7.7. 5..78 5. 7 5. 5 5. 78.7. 7.78 5.. 5. Q 5. Q 78 Q 7 7 Q 5.Q 5 5 Q 7

Stp : Soving th systm of inar uations givs us th unknown dispacmnts and rotations, i.. Q.9 in. in.9 rad.9 in.8 in 5.88 rad