nonlinear introduction

Upload: san

Post on 05-Feb-2018

215 views

Category:

Documents


0 download

TRANSCRIPT

  • 7/21/2019 Nonlinear Introduction

    1/16

    Basic Introduction to

    Nonlinear AnalysisRonald D. ZiemianBucknell University

    The function of a structural engineer isto design not to analyze

    Norris and Wilbur1960

    Analysis is a means to an endrather than the end itself.

    Role of the analysis:

    forces, moments and deflections design equations

    insight into the behavior of a structure better the understanding, better the design

    Limit States Design:Prior to limit of resistance, significant nonlinearresponse, including- geometrical effects (P-, P-)- material effects (yielding, cracking, crushing)

    - combined effects

    AISC Ch. C: P-, P- (App. 7)App. 1: Inelastic Design

    Seismic: Pushover AnalysisOther: Progressive Collapse

    Impetus:

    Nonlinear Analysis

    Limit States Design

    Available Software

    Research

    Education

    Nonlinear AnalysisHand methods

    Second-order effects i.e. Moment Amplification Factors (B1 and B2 factors)

    Material nonlinear effects i.e. plastic analysis (upper and lower bound theories)

    Computer Methods (focus of todays lecture) Lots of variations

    all use same basic concepts (most important to today)

    one approach will be presented (basis for MASTAN2)

    Please keep in mind All methods are approximate

    Not a substitute, but a complement to good engineering

    Lecture Overview

    Brief Introduction (done!)Computer Structural Analysis (Review?)

    Basis for Material Nonlinear ModelsIncorporating Geometric Nonlinear Behavior

    Critical Load Analysis

    Overview of MASTAN2 software

    Summary and Concluding Remarks

  • 7/21/2019 Nonlinear Introduction

    2/16

    How does the computer getthese results?

    State-of-the-Art Crystal Ball? Not quite.

    By applying 2 requirements and 1 translator Two Requirements:

    Equilibrium (equations in terms of Fs and Ms, 1 per d.o.f.)

    Compatibility (equations in terms of s and s, 1 per d.o.f.)

    Translator apples to oranges Constitutive Relationship (i.e. Hookes Law, = E )

    Generalized to Force-to-Displacement (i.e. F=k)

    Re-write equilibrium eqs. in terms of unknown displacements

    # of Equil. Eqs. = # of Unknown Displacements

    Equilibrium Equations

    A B C

    D 40 kips

    CDxD

    BDxD

    ADxD

    xD

    fff

    Fux-

    ++=

    =

    40

    0:d.o.f.

    40 kips

    BDxDf

    BDyDf

    Membe

    rBD

    BDxDf

    BDyDf

    MemberCD

    CDxDf

    CDyDf

    CDyDf

    CDxDf

    ADyDf

    Membe

    rAD

    ADxDf

    ADyDf

    ADxDf

    D

    Free BodyDiagram

    CDyD

    BDyD

    ADyD

    yD

    fff

    Fvy-

    ++=

    =

    0

    0:d.o.f.

    Translator: Forces Displacements

    ui

    vi

    vj

    uj

    i

    j

    (Deformed)

    i

    j(Undeformed)

    fxi

    fxj

    fyj

    fyi

    Big Question:Where do these known stiffness coefficients kscome from? Little Answer:

    Function of members material and geometricproperties, including its orientation.

    jjiiyj

    jjiixj

    jjiiyi

    jjiixi

    vkukvkukf

    vkukvkukf

    vkukvkukf

    vkukvkukf

    44434241

    34333231

    24232221

    14131211

    +++=

    +++=

    +++=

    +++=

    F for all members

    BDD

    BDBDD

    BDBDB

    BDBDB

    BDBDyD

    BD

    D

    BDBD

    D

    BDBD

    B

    BDBD

    B

    BDBD

    xD

    BDD

    BDBDD

    BDBDB

    BDBDB

    BDBDyB

    BDD

    BDBDD

    BDBDB

    BDBDB

    BDBDxB

    vkukvkukfvkukvkukf

    vkukvkukf

    vkukvkukf

    44434241

    34333231

    24232221

    14131211

    +++=+++=

    +++=

    +++=

    Member BD:

    ADD

    ADADD

    ADADA

    ADADA

    ADADyD

    ADD

    ADADD

    ADADA

    ADADA

    ADADxD

    ADD

    ADADD

    ADADA

    ADADA

    ADADyA

    ADD

    ADADD

    ADADA

    ADADA

    ADADxA

    vkukvkukf

    vkukvkukf

    vkukvkukf

    vkukvkukf

    44434241

    34333231

    24232221

    14131211

    +++=

    +++=

    +++=

    +++=

    Member AD:

    A B C

    DuD

    vD

    uA

    vAuB

    vBuC

    vC

    Member CD:

    CDD

    CDCDD

    CDCDC

    CDCDC

    CDCDyD

    CD

    D

    CDCD

    D

    CDCD

    C

    CDCD

    C

    CDCD

    xD

    CDD

    CDCDD

    CDCDC

    CDCDC

    CDCDyC

    CDD

    CDCDD

    CDCDC

    CDCDC

    CDCDxC

    vkukvkukfvkukvkukf

    vkukvkukf

    vkukvkukf

    44434241

    34333231

    24232221

    14131211

    +++=+++=

    +++=

    +++=

    Substituting into Equil. Eqs.

    A B C

    D40 kN

    CD

    xD

    BD

    xD

    AD

    xD

    xD

    fff

    Fux-

    ++==

    40

    0:d.o.f.

    ADD

    ADADD

    ADADA

    ADADA

    ADADxD vkukvkukf 34333231 +++=

    Member AD:

    BDD

    BDBDD

    BDBDB

    BDBDB

    BDBDxD vkukvkukf 34333231 +++=

    Member BD:

    Member CD:CDD

    CDCDD

    CDCDC

    CDCDC

    CDCDxD vkukvkukf 34333231 +++=

    ( )( )( )CD

    D

    CDCD

    D

    CDCD

    C

    CDCD

    C

    CD

    BD

    D

    BDBD

    D

    BDBD

    B

    BDBD

    B

    BD

    AD

    D

    ADAD

    D

    ADAD

    A

    ADAD

    A

    AD

    D

    vkukvkuk

    vkukvkuk

    vkukvkuk

    u

    34333231

    34333231

    3433323140

    :d.o.f.

    +++

    ++++

    ++++=

    Substituting into Equil. Eqs. (cont.)

    A B C

    D40 kN

    CD

    yD

    BD

    yD

    AD

    yD

    yD

    fff

    Fvy-

    ++==

    0

    0:d.o.f.

    ( )( )( )CDD

    CDCD

    D

    CDCD

    C

    CDCD

    C

    CD

    BD

    D

    BDBD

    D

    BDBD

    B

    BDBD

    B

    BD

    AD

    D

    ADAD

    D

    ADAD

    A

    ADAD

    A

    AD

    D

    vkukvkuk

    vkukvkuk

    vkukvkuk

    v

    44434241

    44434241

    44434241

    0

    :d.o.f.

    +++

    ++++

    ++++=

    Member AD:

    Member BD:

    Member CD:

    ADD

    ADADD

    ADADA

    ADADA

    ADADyD vkukvkukf 44434241 +++=

    BDD

    BDBDD

    BDBDB

    BDBDB

    BDBDyD vkukvkukf 44434241 +++=

    CDD

    CDCDD

    CDCDC

    CDCDC

    CDCDyD vkukvkukf 44434241 +++=

  • 7/21/2019 Nonlinear Introduction

    3/16

    So, where are we at?

    We have two equilibrium equations (1 perd.o.f.) in terms of a lot of displacements:

    ( )( )( )

    ( )( )( )CD

    D

    CDCD

    D

    CDCD

    C

    CDCD

    C

    CD

    BD

    D

    BDBD

    D

    BDBD

    B

    BDBD

    B

    BD

    AD

    D

    ADAD

    D

    ADAD

    A

    ADAD

    A

    AD

    D

    CD

    D

    CDCD

    D

    CDCD

    C

    CDCD

    C

    CD

    BD

    D

    BDBD

    D

    BDBD

    B

    BDBD

    B

    BD

    AD

    D

    ADAD

    D

    ADAD

    A

    ADAD

    A

    AD

    D

    vkukvkuk

    vkukvkuk

    vkukvkukv

    vkukvkuk

    vkukvkuk

    vkukvkuku

    44434241

    44434241

    44434241

    34333231

    34333231

    34333231

    0:

    40:

    +++

    ++++

    ++++=

    +++

    ++++

    ++++=

    What card havent we played yet?

    Compatibility Eqs. (consistent deflections)

    CD

    CC 0vv ==

    CD

    CC 0uu ==AD

    AA 0vv ==

    AD

    AA 0uu ==

    BD

    BB 0vv ==

    BD

    BB 0uu ==

    Member-to-Support

    A B C

    D

    D

    MemberCDM

    embe

    rAD

    Membe

    rAD

    Membe

    rBD

    Membe

    rBD

    CD

    D

    BD

    D

    AD

    DD vvvv ===

    CD

    D

    BD

    D

    AD

    DD uuuu ===

    Member-to-MemberDv

    AD

    Dv

    ADDu

    BD

    Dv

    BD

    Du

    CD

    Dv

    CD

    Du

    Du

    Time for some serious simplifying

    Applying Compatibility to Equil. Eqs.:

    ( )( )( )

    ( )( )( )CD

    D

    CDCD

    D

    CDCD

    C

    CDCD

    C

    CD

    BD

    D

    BDBD

    D

    BDBD

    B

    BDBD

    B

    BD

    AD

    D

    ADAD

    D

    ADAD

    A

    ADAD

    A

    AD

    D

    CD

    D

    CDCD

    D

    CDCD

    C

    CDCD

    C

    CD

    BD

    D

    BDBD

    D

    BDBD

    B

    BDBD

    B

    BD

    AD

    D

    ADAD

    D

    ADAD

    A

    ADAD

    A

    AD

    D

    vkukvkuk

    vkukvkuk

    vkukvkukv

    vkukvkuk

    vkukvkuk

    vkukvkuku

    44434241

    44434241

    44434241

    34333231

    34333231

    34333231

    0:

    40:

    +++

    ++++

    ++++=

    +++

    ++++

    ++++=

    Which simplifies to

    All = 0

    All = uD

    All = vD

    After simplifying

    ( ) ( ) DCDBDADDCDBDADDD

    CDBDAD

    D

    CDBDAD

    D

    vkkkukkkv

    vkkkukkku

    0:

    40:

    444444434343

    343434333333

    +++++=

    +++++=

    A B C

    D40 kN

    Solve for Unknown Displacements

    Since ks are known, we have2 Equations and 2 Unknowns

    ? ?

    ? ?

    uD = # and vD = ##

    ##

    #

    0

    0

    0

    0

    0

    0

    ====

    ====

    ==

    ==

    ==

    ==

    ==

    ==

    CDD

    BDD

    ADDD

    CD

    D

    BD

    D

    AD

    DD

    CD

    CC

    CD

    CC

    BD

    BB

    BD

    BB

    AD

    AA

    AD

    AA

    vvvv

    uuuu

    vv

    uu

    vv

    uu

    vv

    uu

    With all displacements, solve formember forces

    ADD

    ADADD

    ADADA

    ADADA

    ADADyD

    ADD

    ADADD

    ADADA

    ADADA

    ADADxD

    ADD

    ADADD

    ADADA

    ADADA

    ADADyA

    ADD

    ADADD

    ADADA

    ADADA

    ADADxA

    vkukvkukf

    vkukvkukf

    vkukvkukf

    vkukvkukf

    44434241

    34333231

    24232221

    14131211

    +++=

    +++=

    +++=

    +++=

    Member AD:

    Member BD:

    BDD

    BDBDD

    BDBDB

    BDBDB

    BDBDyD

    BDD

    BDBDD

    BDBDB

    BDBDB

    BDBDxD

    BDD

    BDBDD

    BDBDB

    BDBDB

    BDBDyB

    BDD

    BDBDD

    BDBDB

    BDBDB

    BDBDxB

    vkukvkukf

    vkukvkukf

    vkukvkukf

    vkukvkukf

    44434241

    34333231

    24232221

    14131211

    +++=

    +++=

    +++=

    +++=

    Member CD:

    CDD

    CDCDD

    CDCDC

    CDCDC

    CDCDyD

    CDD

    CDCDD

    CDCDC

    CDCDC

    CDCDxD

    CDD

    CDCDD

    CDCDC

    CDCDC

    CDCDyC

    CDD

    CDCDD

    CDCDC

    CDCDC

    CDCDxC

    vkukvkukf

    vkukvkukf

    vkukvkukf

    vkukvkukf

    44434241

    34333231

    24232221

    14131211

    +++=

    +++=

    +++=

    +++=

    Summary of Computer Approach For each d.o.f., write an equilibrium equation:

    Fexternal = fmember Re-write (translate) each member force in terms of

    its end displacements (Stiffness Eqs.)

    fmember = k member member end Substitute Stiffness Eqs. into above Equil. Eqs.

    Simplify Equil. Eqs. by applying member-to-memberand member-to-support compatibility conditions

    Solve n Equil. Eqs. for the n unknown displacements

    Use Stiffness Eqs. to calculate member forces

    Apply Equil. Eqs. to solve for reactions

  • 7/21/2019 Nonlinear Introduction

    4/16

    Lots of Questions

    So, this is how most commercial programs such

    as SAP2000, RISA, STAAD, etc. get theanswer? Yes! Known as Direct Stiffness Method

    So, all such programs will give the same answer? Yes, as long as it is a static 1st-order elastic analysis.

    Wait a minuteIs this the basic analysisprocedure for the finite element method? Yes! Bit more tricky to get ks, s, and s

    Two Big Questions

    Where do those stiffness coefficients

    come from? You mean the ones that relate member endforces to member end displacements?

    Yeah, those ks !

    What happens when we go static nonlinearor even dynamic? Same basic procedure, but apply loads in

    increments and perform a series of analyses.Then, sum incremental results.

    < Much more to come on this! >

    Important PointsThe only opportunity for most computer

    analysis software to model the actualbehavior of the structure is through themember stiffness terms.

    So, to include first-order effects

    second-order effects material nonlinear behavior

    Must modify member stiffness!!!

    Lets review member stiffness

    Stiffness Coefficients, ks

    Lets start with high school physics Extension Spring Experiment

    After:

    F

    F

    k

    1

    F = k

    Before:Force

    Displacement

    After:

    M

    k

    1

    M = k

    M

    Stiffness Coefficients, ks (cont.)

    More advanced high school physics lab Rotational Spring Experiment

    Before:Moment

    Rotation

    How about real structural members?Axial force member

    F = k(A,L,E)

    Stiffness k function of: Geometry: Area and Length (A,k & L, k)

    Material: Elastic Modulus (E,k)

    After:F

    F

    k1

    Before:

  • 7/21/2019 Nonlinear Introduction

    5/16

    How about real members? (cont.)Flexural members

    Before:

    Stiffness k function of: Geometry: Moment of Inertia & Length (I,k& L, k) Material: Elastic Modulus (E,k)

    Before:

    After:

    F

    After:M

    F = k(I,L,E)M = k(I,L,E)

    Other factor impacting stiffnessOrientation of member

    consider axial force member:

    kv=EA/Lkh=0

    VerticalMember kh=EA/L

    kv=0

    HorizontalMember

    Orientation of axial force member

    Important Point: Less vertical a member,the less stiffness to resist vertical loads.

    kv=(sin2)EA/L

    kh=(cos2)EA/L

    Summary: Three Perspectives

    Reality: What you see

    F1

    Three Perspectives (cont.)What you see on your computer screen:

    DL

    WL

    Collection of elementsconnected by sharingcommon nodes

    Three PerspectivesWhat your computer actually sees:

    Assemblage ofequivalent springs {F} = [K]{}

    k

    DLWL

  • 7/21/2019 Nonlinear Introduction

    6/16

    Analysis Review: Key Points

    Reviewed the Direct Stiffness Method Equilibrium Translator F() Compatibility

    Response of structure controlled bystiffness of members (a.k.a. springs)

    First-order elastic stiffness of memberfunction of: Material Property (E)

    Geometric Properties (A, I, L, and orientation)

    Time to go nonlinearlets begin with material nonlinear

    Material Nonlinear (Inelastic)

    Best place to start is with a tensile test

    = P/A

    = /L

    E

    1

    yield Perfectly Plastic (E = 0)

    Ela

    stic

    Normal Stress: Structural MembersFor typical structural steel members

    (L/d>10), elastic/inelastic behaviorcontrolled by normal stresses s actingalong the length axis of the member.

    Normal stress produced by: Axial force (P/A)

    Major and/or minor axis flexure (Mc/I)

    Combination of above effects (i.e. P/A + Mc/I) Warping (not today!)

    We will assume elastic-perfectly-plasticmaterial (often done for steel)

    E=0 k=0

    Post-Yield:

    = y

    Inelastic Behavior: Axial ForceOriginally:

    P=P/A=0

    Yield:

    Py

    Py=Ay

    = yPlastic Hingeat P = Py or

    when P/Py = 1.0

    k=EA/L1

    Elastic:

    P

    < y

    Inelastic Behavior:Flexure

    Mp=Zy

    M

    My=Sy

    MA

    A

    k=4EI/L

    1EIelastic

    -y +y

    M < My

    Section A-A

    Inelastic Behavior:Flexure (cont.)

    MA

    A

    Mp=Zy

    k=4EI/L

    1

    M

    My=Sy

    Section A-A

    EIelastic

    -y +y

    M=Sy=My

  • 7/21/2019 Nonlinear Introduction

    7/16

    Inelastic Behavior:Flexure (cont.)

    Mp=Zy

    k=4EI/L

    1

    M

    My=Sy

    Section A-A -y +y

    E=0

    EI

  • 7/21/2019 Nonlinear Introduction

    8/16

    Inelastic Behavior:Combination P & M

    Mp=Zy

    k=4EI/L

    1

    M

    MA

    AP

    -y

    +y

    M < Mp

    = Mc/I + P/A

    +y

    -y

    Plastic

    Elas

    tic

    Inelastic Behavior:Combination P & Mfor Plastic Hinge

    MA

    AP

    +y

    -y

    +y

    -y

    +y

    -y= +

    M < Mp

    M/Mp < 1

    P < Py

    P/Py < 1

    Fully yielded

    section when:

    Elastic

    YieldSurface

    Plastic Hinge Criterion:P/Py

    MA

    AP

    M/Mp

    1.0

    1.0

    -y P=PyM=0

    +y

    -y P

  • 7/21/2019 Nonlinear Introduction

    9/16

    NoAxial Force

    Second-Order Effects

    A.K.A. Geometric Nonlinear BehaviorEquilibrium Equations

    Reality: Should be formulated on deformedshape

    Difficulty: Deformed shape (deformations) isa function of the member forces, which are inturn a function of the deformations(Chicken n Egg)

    Remedy: Perform a series of analyses withloads applied in small increments and updategeometry after each load increment.

    Truss is susceptible to2nd-Order effects,luckily is oftenquite small.

    Equilibrium Equations Formulated on

    Undeformed Shape Formulated on

    Deformed Shape

    HP P

    H

    Different reactions and member forces.

    k=3EI/L3

    1

    H

    Equilibrium Equations

    H

    P

    M=HL

    k=3EI/L3

    1

    H

    keff < k1

    Effectivelateralstiffnessis reduced!

    Formulated onUndeformed Shape

    H

    P

    M=HL+P

    Formulated onDeformed Shape

    Focus on Lateral Stiffness Formulated on Undeformed Shape: Linear Response

    Before:

    kspringP

    H

    After:

    klateral = kspring1

    H

    Lateral Stiffnessis slope of H-response curve

  • 7/21/2019 Nonlinear Introduction

    10/16

    Focus on Lateral Stiffness (cont.) Formulated on Deformed Shape: Nonlinear Response

    Before:

    kspring P

    H

    After:

    klateral < kspring1

    H kspring1

    Effective lateralstiffness is reduced

    Focus on Lateral Stiffness (cont.) Equilibrium Formulated on Deformed Shape

    PH kspring

    R=kspring

    Mo = 0 RL = HL + P

    R = H + P/L

    H = klateral with klateral = kspring P/L

    Lateral Stiffness (slope of response curve)

    Pt. o

    L L

    Lets start by assuming L L,

    kspring = H + P/L

    H = kspring P/L

    H = (kspring P/L)

    Some thoughts here

    This simple analysis becomes less accurate as/L becomes large (i.e. /L >> 1/5) Remedy: Perform an incremental analysis and update

    geometry after each load incrementhence, limit/L in each step to some small amount

    Keep in mind serviceability limits are oftensomething like /L < 1/400

    Most importantly, klateral = kspring P/L takes onthe form:

    k2nd-Order El. = k1st-Order El. + kgGeometric Stiffness

    Geometric Stiffness

    Effective lateral stiffness of a member: decreases as a member is compressed

    kg is negative for compressive P

    backpacker example

    increases when subjected to tension kg is positive for tensile P guitar string example

    Employing geometric stiffness approach Other methods exist (i.e. stability functions)

    How about real members? (recall)Flexural members subjected to axial force

    F

    Stiffness k function of: Geometry: Moment of Inertia & Length (I,k& L, k)

    Material: Elastic Modulus (E ,k)

    Axial Force: Compressive (P, k)

    M

    P

    P

    F = k(I,L,E,P) M = k(I,L,E,P)

    Closer look at stiffness terms

    Flexural members subjected to axial force

    M = k(I,L,E,P) withk = 4EI/L 2PL/15

    M

    P

    F = k(I,L,E,P) withk = 12EI/L3 6P/5L

    F

    P

    Again, basic form:

    k2nd-Order El. = k1st-Order El. + kg

  • 7/21/2019 Nonlinear Introduction

    11/16

    Geometric Nonlinear AnalysisEmploy Direct Stiffness Method applying loads

    in increments: Solve Equil. Eqs. {dF} = [K]{d}

    At start of increment, modify member stiffnessto account for presence of member forces (suchas axial force): k = kelastic + kg with kg = geometric stiffness

    At end of increment, update model of structuralgeometry to include displacements

    Continue to accumulate results of load increments(i = i-1 + d and fi = fi-1 + df) until all of load isapplied or elastic instability is detected.

    Comparison: 1st- and 2nd-Order Analysis Results

    Moments increase by ~10%

    2nd-Order Inelastic AnalysisEmploy Direct Stiffness Method applying loads

    in increments: Solve Equil. Eqs. {dF} = [K]{d}At start of increment, modify member stiffness

    to account for presence of member forces andany yielding:

    k = kelastic + kgeometric + kplasticAt end of increment, update model of structural

    geometry to include displacementsContinue to accumulate results of load increments

    (i = i-1 + d and fi = fi-1 + df) until all of load isapplied or inelastic instability is detected.

    Critical Load Analysis (Basics)Definition: Critical or buckling load is the

    load at which equilibrium may be satisfiedby more than one deformed shape.

    Big Q: How does computer software calculate this?

    P Solution

    #1 P

    Solution

    #2

    P

  • 7/21/2019 Nonlinear Introduction

    12/16

    Critical Load Analysis (Background)

    Elastic stiffness of a member k = kel + kg kel is f(A or I, L, and E)

    kg is f(P,L), also note directly proportional to PElastic stiffness of structure [K] = k

    [K] = [Kel] + [Kg] [Kg] directly proportional to applied force

    i.e. Double applied forces, hence, double internal forcedistribution and double [Kg]

    To the computer, buckling will occur whenour equilibrium equations {F} = [K]{} permitnon-unique solutions, e.g. det[K] = 0.

    Example

    Demonstrate computational

    method for calculating theelastic critical load (bucklingload) for the structuralsystem shown.

    A, IL, E

    PRigid Beam

    Example: Key Stiffness Terms

    vert

    A, IL, E

    HP

    Lateral Stiffness:

    klateral

    H = klaterallatklateral = 12EI/L3 6P/5L

    kvertical

    Vertical Stiffness:

    P = kverticalvert

    latPRigid BeamH

    Example: Solution

    1. Apply reference load, and use 1st-order elastic analysis to obtaininternal force distribution.

    A, IL, E

    PRigid Beam

    klateral = 12EI/L3 6P/5L

    klateral = 0 when P = 10EI/L2

    (Ptheory=9.87EI/L2)Pcr = P = 10EI/L2

    Pcr

    2. Determine load factor at whichsystem stiffness degrades topermit buckling.

    P = 400 kipsW1265

    100 kips

    =9.4

    =1.7

    Thoughts on Critical Load AnalysisComputer analysis for a large system:

    First, apply reference and perform analysis Solve equilibrium eqs. {Fref} = [K]{}

    With displacements solve for member forces

    Second, assemble [Kel

    ] and [Kg

    ] based on {Fref

    }

    Finally, determine load factor causing instability;computationally this means find load factor at which[K]=[Kel]+[Kg] becomes singular Determine at which det([Kel]+[Kg) = 0

    Eigenvalue problem: Eigenvalues = Critical Load Factors, sEigenvectors = Buckling modes

    Accuracy increases with more elements percom ression members (2 often ade uate)

  • 7/21/2019 Nonlinear Introduction

    13/16

    Basic Introduction Complete

    Acquire nonlinear analysis software Commercial programs Educational software (i.e. MASTAN2)

    Where do I go from here? (Learning to drive)

    Review the slides (Read the drivers manual) Acquire nonlinear software (Borrow a friends car)

    Work lots of examples (Go for a drive, scary at first) Apply nonlinear analysis in design (NASCAR? not quite)

    MASTAN2:

    -Educational software-GUI commercial programs

    -Limited # of pre- and post-processing optionsto reduce learning curve

    -Suite of linear and nonlinear 2D and 3Danalysis routines

    -Available with textbookor online at no cost

    www.mastan2.com orwww.aisc.org [Steel Tools]

    MASTAN2

    1st-Order Elastic: [Ke]{}={F}

    Levels of Analysis:

    2nd-Order Elastic: [Ke + Kg]{d}={dF}

    1st-Order Inelastic: [Ke + Kp]{d}={dF}

    2nd-Order Inelastic: [Ke + Kg + Kp]{d}={dF}

    Critical Load: [Ke + Kg]{d}={0}

    Yield Surface:

    Function of P, Mmajor, and Mminor

    1

    st

    -Order

    Elastic

    Lateral displacement,

    2

    nd

    -Order

    Elastic

    Elastic Stability Limit

    Hes

    Lateral

    load,

    H

    Inelastic Limit Load

    st

    -Order

    Inelastic

    Hp

    H= P

    H= P

    P

    P

    Elastic Critical Load

    Inelastic Critical Load

    He

    Hi

    2

    nd

    -Order

    Inelastic

    Inelastic Stability

    Limit

    His

    Actual

    MASTAN2

    1

    st

    -Order

    Elastic

    u

    in.)

    Load

    Factor

    First

    hinge

    Second hinge

    8783 1

    st

    -Order)

    8783 2

    nd

    -Order)

    14

    1831

    1665

    1720

    Mechanism Moments at Limit in. kips)

    Planar Frame:

    2

    nd

    -Order

    Elastic

    2

    nd

    -Order

    Inelastic

    1

    st

    -Order

    Inelastic

    E= 29,000 ksi

    y = 36 ksi

    Summary and Conclusions

    Provided an introduction to nonlinear analysis Review of direct stiffness method

    Material nonlinear analysis (Inelastic hinge)

    Geometric nonlinear analysis (2nd-Order) 2nd-Order inelastic analysis (combine above)

    Critical load analysis (eigenvalue analysis)

    Nonlinearthink modifying member stiffness!Overview and availability of MASTAN2Now, its your turn to take it for a spin

  • 7/21/2019 Nonlinear Introduction

    14/16

    Appendix

    Several examples to try out

    Solutions by MASTAN2Need a reference text with many

    examples? see Matrix Structural Analysis,2nd Ed., by McGuire, Gallagher, and Ziemian(Wiley, 2000)

    See tutorial that comes with MASTAN2OK, time to jump in and start driving

    Multi-story

    Frame:

    Demonstrate:

    2

    nd

    -Order Inelastic Analysis

    Non-proportional loading

    Gravity Load

    a

    t

    e

    r

    a

    l

    o

    a

    d

    E= 29,000 ksi

    y = 36 ksi

    P

    0.5P

    P

    P

    P

    P

    P

    P

    P

    50 kips 25 kips25 kips

    MASTAN2

    Model:

    Gravity Load

    2nd-Order

    Inelastic

    Analysis:

    P = 30 kips

    2nd-Order

    Inelastic

    Analysis:

    P = 37.7 kips

    2nd-Order

    Inelastic

    Analysis:

    Limit State:

    P

    lim

    = 40.3 kips

  • 7/21/2019 Nonlinear Introduction

    15/16

    2nd-Order

    Inelastic

    Analysis:

    Post-limit State:

    P = 40 kips

    1st

    Plastic Hinge (P=30 kips)

    2nd-Order

    Inelastic

    Analysis:

    Strength Limit State (P=40.3 kips)

    All members:

    A = 9.34810-2 in2

    I = 6.954 10-4 in4

    E = 29,000 ksi

    P = 100 lbs

    20

    3 @ 20

    3P

    3P

    Truss Hoff et al.):

    Demonstrate:

    Elastic Critical Load

    2

    nd

    -Order Elastic

    Experimental P

    limit

    =220 lbs)

    MASTAN2

    Model:

    Elastic Critical Load: P

    cr

    = 210.7 lbs 2

    nd

    -Order Elastic: P

    lim

    = 210 lbs

  • 7/21/2019 Nonlinear Introduction

    16/16

    Response

    Curves:

    Elastic Critical Load

    2

    nd

    -Order Elastic

    1

    st

    -Order

    Elastic

    P = 210 lbs

    Beam-Column:

    Demonstrate:

    Elastic Critical Load Analysis

    1. Flexural Buckling

    =0.0)

    2. Torsional Flexural Buckling

    =0.04)

    L = 24

    PM =

    PL

    M = PL

    W2476E = 29,000 ksi

    MASTAN2

    Model:

    Elastic Critical Load

    = 0.0)

    Isometric View Elevation View Plan View

    Elastic Critical Load Analysis = 0.04)

    Isometric View Elevation View Plan View

    v

    mi

    in.)

    Load

    Factor

    Suspension System:

    2

    nd

    -Order

    Inelastic

    1

    st

    -Order

    Inelastic

    A = 5.40 in2

    y = 150 ksi

    A = 50 in2

    I = 20,000 in4

    Z = 1,000 in3

    y = 50 ksi

    1.5

    1

    0.5

    2

    0

    0 50 100 150 25000

    1

    st

    -Order Inelastic

    2

    nd

    -Order Inelastic

    Hinge Formation