Lecture 9. Stability of Elastic Structures. Lecture 10. Advanced Topic in Column Buckling

Similar documents
Lecture 6: Moderately Large Deflection Theory of Beams


Module 22: Simple Harmonic Oscillation and Torque

3.10 Implications of Redundancy

VTU-NPTEL-NMEICT Project

Work and energy method. Exercise 1 : Beam with a couple. Exercise 1 : Non-linear loaddisplacement. Exercise 2 : Horizontally loaded frame

Nonlinear Analysis of Spatial Trusses

THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS

1D Heat Propagation Problems

CHAPTER 9. Columns and Struts

Separation of Variables and a Spherical Shell with Surface Charge

Torsion and shear stresses due to shear centre eccentricity in SCIA Engineer Delft University of Technology. Marijn Drillenburg

Finite element method for structural dynamic and stability analyses

221B Lecture Notes Notes on Spherical Bessel Functions

Unit 48: Structural Behaviour and Detailing for Construction. Deflection of Beams

Lecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String

Technical Data for Profiles. Groove position, external dimensions and modular dimensions

UI FORMULATION FOR CABLE STATE OF EXISTING CABLE-STAYED BRIDGE

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur

CS229 Lecture notes. Andrew Ng

SEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l

Instructional Objectives:

Gauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law

Application of the Finite Fourier Sine Transform Method for the Flexural-Torsional Buckling Analysis of Thin-Walled Columns

The EM Algorithm applied to determining new limit points of Mahler measures

David Eigen. MA112 Final Paper. May 10, 2002

Math 124B January 31, 2012

LECTURE NOTES 8 THE TRACELESS SYMMETRIC TENSOR EXPANSION AND STANDARD SPHERICAL HARMONICS

APPENDIX C FLEXING OF LENGTH BARS

4 1-D Boundary Value Problems Heat Equation

Lecture Notes for Math 251: ODE and PDE. Lecture 32: 10.2 Fourier Series

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

Malaysian Journal of Civil Engineering 30(2): (2018)

FOURIER SERIES ON ANY INTERVAL

СРАВНИТЕЛЕН АНАЛИЗ НА МОДЕЛИ НА ГРЕДИ НА ЕЛАСТИЧНА ОСНОВА COMPARATIVE ANALYSIS OF ELASTIC FOUNDATION MODELS FOR BEAMS

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form

Lecture Notes 4: Fourier Series and PDE s

1 Equivalent SDOF Approach. Sri Tudjono 1,*, and Patria Kusumaningrum 2

The Bending of Rectangular Deep Beams with Fixed at Both Ends under Uniform Load

STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION

More Scattering: the Partial Wave Expansion

AST 418/518 Instrumentation and Statistics

arxiv: v1 [hep-th] 10 Dec 2018

Lecture 8 February 18, 2010

Laboratory Exercise 1: Pendulum Acceleration Measurement and Prediction Laboratory Handout AME 20213: Fundamentals of Measurements and Data Analysis

Lobontiu: System Dynamics for Engineering Students Website Chapter 3 1. z b z

Math 124B January 17, 2012

STRUCTURAL ANALYSIS - I UNIT-I DEFLECTION OF DETERMINATE STRUCTURES

Strain Energy in Linear Elastic Solids

MECHANICAL ENGINEERING

Math 220B - Summer 2003 Homework 1 Solutions

Add Math (4044/02) (+) x (+) 2. Find the coordinates of the points of intersection of the curve xy 2 the line 2y 1 x 0. [5]

22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: The Moment Equations

HYDROGEN ATOM SELECTION RULES TRANSITION RATES

4 Separation of Variables

Bending Analysis of Continuous Castellated Beams

> 2 CHAPTER 3 SLAB 3.1 INTRODUCTION 3.2 TYPES OF SLAB

THINKING IN PYRAMIDS

Related Topics Maxwell s equations, electrical eddy field, magnetic field of coils, coil, magnetic flux, induced voltage

A PROCEDURE ON STUDYING A CRITICAL FORCE ON FEM AS PLASTIC DEFORMATIONS ARISE WITHIN STRAIGHT TRUSSES

1 Equations of Motion 3: Equivalent System Method

Problem Set 6: Solutions

C1B Stress Analysis Lecture 1: Minimum Energy Principles in Mechanics Prof Alexander M. Korsunsky Hilary Term (January 08)

Chemical Kinetics Part 2

Buckling of elastic structures under tensile loads

KINGS COLLEGE OF ENGINEERING DEPARTMENT OF CIVIL ENGINEERING. Question Bank. Sub. Code/Name: CE1303 Structural Analysis-I

An approximate method for solving the inverse scattering problem with fixed-energy data

Week 6 Lectures, Math 6451, Tanveer

6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7

Version 2.2 NE03 - Faraday's Law of Induction

$, (2.1) n="# #. (2.2)

Nonperturbative Shell Correction to the Bethe Bloch Formula for the Energy Losses of Fast Charged Particles

Structural Analysis III Revised Semester 2 Exam Information. Semester /9

SVM: Terminology 1(6) SVM: Terminology 2(6)

Physics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions

Automobile Prices in Market Equilibrium. Berry, Pakes and Levinsohn

Introduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.

Solution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima

In-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017

Quantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18

Delhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web: Ph:

C. Fourier Sine Series Overview

7. CREST-TO-TROUGH WAVE HEIGHT DISTRIBUTION

First-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries

Strauss PDEs 2e: Section Exercise 1 Page 1 of 7

A Brief Introduction to Markov Chains and Hidden Markov Models

Legendre Polynomials - Lecture 8

2.1. Cantilever The Hooke's law

Lecture 17 - The Secrets we have Swept Under the Rug

Vibrations of beams with a variable cross-section fixed on rotational rigid disks

Chemical Kinetics Part 2. Chapter 16

Supporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers

Course 2BA1, Section 11: Periodic Functions and Fourier Series

V.B The Cluster Expansion

Modal analysis of a multi-blade system undergoing rotational motion

Model Solutions (week 4)

SIMULATION OF TEXTILE COMPOSITE REINFORCEMENT USING ROTATION FREE SHELL FINITE ELEMENT

Methods for Ordinary Differential Equations. Jacob White

DYNAMIC RESPONSE OF CIRCULAR FOOTINGS ON SATURATED POROELASTIC HALFSPACE

Cantilever Beam Static and Dynamic Response Comparison with Mid-Point Bending for Thin MDF Composite Panels

Transcription:

Lecture 9 Stabiity of Eastic Structures Lecture 1 Advanced Topic in Coumn Bucking robem 9-1: A camped-free coumn is oaded at its tip by a oad. The issue here is to find the itica bucking oad. a) Suggest a simpe form of the bucked of the coumn, satisfying kinematic boundary conditions. b) Use the Rayeigh-Ritz quotient to find the approximate vaue of the bucking oad. c) Come up with another bucking shape which woud give you a ower vaue for the bucking oad. d) Find the exact soution of the probem and show the convergence of the approximate soution to the exact soution. Foow the exampe of a pin-pin coumn, which is presented in the notes of Lecture 9. 1

robem 9-1 Soution: a) Kinematic boundary condition, in term of shape function I( x ), for a camped-free coumn is () '() Choose a bucking shape ( x) x '( x ) x ''( x ) b) Use Rayeigh-Ritz Quotient, the itica bucking oad is N c I'' I'' dx II ' 'dx u dx x u xdx N c L

c) Choose a bucking shape simiar to a cantiever beam I( x ) x Lx I '( x ) x 6Lx I ''( x ) 6x 6L N c I'' I'' dx II ' 'dx ³ 6x 6L dx x 6Lx dx ³ 1L 4 L 5 5 Compare to the resut in b), N c.5 L N c, this bucking shape givers a ower vaue L d) Choose bucking shape I( x) S 1 cos x L

S S I '(x) sin x L L S I ''( x ) cos x L L N c I'' I''dx II ' 'dx ª S S º «cos x» dx «L L ¼» S sin x dx L L ³ S 4L N c.47 L Check for oca equiibrium of the soution ª S 4 S S w IV N c w '' A «cos x cos x «L L 4L L L»¼ This is the exact soution to the camped-free bucking 4

robem 9-: Consider a camped-free coumn oaded by a compressive force at the free end. a) Determine the itica senderness ratio E it distinguishing between the eastic and pastic bucking response. What is the bucking stress and strain? b) Cacuate the itica pastic bucking oad for E.5E it and the corresponding stress and strain. c) Cacuate the itica eastic bucking oad for E E it and the corresponding stress and strain. d) Compare a three resuts. robem 9- Soution: a) First, find the bending oad For Camed-Free coumn S L S 4L Second, find the bucking stress and strain Reca that S A 4AL V bucking 5

I I r r A A Then Reca that Er E V bucking 4L 4L r L r Er E V 4L 4 V S E 4E H Third, find when V bucking V yied E E it when V y, which is S E V y 4E it S E 4V y E it E V y E it b) E.5E it, the coumn yieds + hits pastic bucking 6

S E S t E t p (From Lecture Note, equation 9.7) 1 E it 4 Ei t ¹ p S S n n 1 E it 4 Eit c) E it, the coumn wi bucke easticay E E E V 4 4E 16E it it V S E H 16E it d) Compare the three resuts Yied y y Eastic Bucking E.67.7 y.67 it it astic Bucking 19.74 E 4 it.5 7.9 1.98 y 4 it it E To simpy our comparison, assume n., E t.5e (*) and reca (*)In order to compare pastic bucking to eastic+yied, we need to make future assumption about the materia properties. V y E it 7

robem 9-: Consider the pin-pin coumn. a) Suggest a poynomia bucking shape function ( x) to improve the approximate soution derived in ecture note. Note that the one used in cass was the paraboic shape. b) Determine the accuracy reative to the exact soution. E I L x robem 9- Soution: S x x a) The exact soution is w sin, use the none-dimensioned vaue F, the Tayor L L series expansion is So we know the shape function must be SF sin SF SF 6 IF ( ) C 1 F C F For x L, the boundary conditions are The first boundary condition gives this doesn t hep. The second boundary condition gives I 1 I' F ¹ I C 1 C 8

d x 1 where d F dx dx C C ' 1 1 ' C1 C 4 C 4 C 1 So we have 4 ( ) C1 C1 4 ( ) C1 We can us the Rayeigh-Ritz Quotient N '' dx ' dx '( ) 1 4 x C1 d dx 4 '( ) C 1 8 16 d dx 1 ''( x) 8C d dx ''( ) 64C d dx 1 4 1 Since we have considered the shape function for d x d L, we must adjust the imits on the integra N ' '' 64C1 d dx dx C 1 8 16 d dx dx 1 L dx dx 4 1 4 64 x 1 dx 18 16 1 after engthy agebra 4 x x dx 4 9

N 1 L b) The resut are compared with the poynomia used in cass and the exact soution Exact Soution araboic Cubic C sin C C 1 C C 1 Coefficient 9.87 1 1 Error N/A 1.5% 1.% Notice how we significanty reduce the error by incuding a higher order term. 1

robem 9-4: resent a step-by-step derivation of the bucking soution of the pin-camped coumn from the oca equiibrium equation. E I L x robem 9-4 Soution: Boundary condition for this probem Start with 4 th order ODE ww wl w' w''l w IV w'' IV w w'' We have an eigenvaue probem 4, i 1 4 Define K 4 ik wc C x C sin Kx C cos Kx 1 4 11

Use the boundary conditions to sove for constants C, 1 C, C and C 4 w w ' w L C 1 C L w L Substitute C C 1, into the above expression w'' L w C1C4 C C 1 4 w' x C KC cos KxKC sin Kx 4 w' C KC C KC C sin KL C 4 cos KL C KLsin KL C 1cos KL 4 w'' x K C sin KxK C cos Kx w'' L K C sin KLK C cos KL 4 4 det> KL sin KL 1cos KLC K sin KL K cos KL C4 K KL KL KL K KL KL cos sin sin 1cos KL cos KL sin KL sin KL KL tan KL cos KL So the equation to sove in order to find is tan KL KL The smaest roots are KL and KL 4.49, we choose KL 4.49 1

L 4.49.16 L.7L.7L 1

robem 9-5: a) Derive the soution for an imperfect camped-free coumn (ike that considered in probem 9-1, foowing a simiar derivation given in the notes for a pin-pin coumn in the notes. b) Find the ratio of current defection ampitude to the ampitude of the initia imperfection such that the resuting oad is 8% of the theoretica bucking oad of a perfect coumn. robem 9-5 Soution: a) wx : shape of initia imperfection wx : actua bucked shape w o x : ampitude of initia imperfection w o : end ampitude of actua imperfection Moment equiibrium of imperfect coumn w w'' w w o erfect coumn w x Assume that the initia imperfection is in the same shape as the bucking shaper From boundary condition o 1cos w x w Ox w x w o 1 cos Ox wl 14

w o O cos OL From moment equiibrium of imperfect coumn erfect coumn L n 1 w wo cosx wo 1 1cosx wo w wo w o S E I O 4L Imperfect coumn w o w o w o w o 1 4L wo wo 1 w o b) When.8 w w o o 1. w w o o 5 15

robem 9-6: The pin-pin eastic coumn of ength L (shown beow) is an I section can bucke in either pane. a) Determine the bucking oad in terms of L, b 1,b, t and E. Assume that t<<b. b) What shoud the ratio of b 1 /b be in order for the probabiity of bucking in either of the bucking panes to be the same? Bonus: What coud happen for very arge width to thickness ratio? b b 1 L x robem 9-6 Soution: a) The moment s of inertia for an I shape oss-section is I I zz 1 b b t 1 1 ¹ 1 tb b 1 6b 1 yy tb 1 1 tb 1 1 tb 1 tb 6 1 1 16

If I yy I zz, the coumn wi bucke in x-z pane S yy S E tb b 6b 1 1 If I yy I zz, the coumn wi bucke in x-y pane S S zz E tb 6 1 b) For the probabiity of bucking in either of the panes to be the same, we want I yy I zz The ony physica soution is 1 tb b 6b 1 1 tb 1 1 6 b 1 b1 1 b b b 1 1.81 b c) If b 1!! t, b!! t, then oca pate bucking my deveop. 17

MIT OpenCourseWare http://ocw.mit.edu.8j / 1.57J Structura Mechanics Fa 1 For information about citing these materias or our Terms of Use, visit: http://ocw.mit.edu/terms.