Failure Load of Plane Steel Frames Using the Yield Surface Method

Similar documents
Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

Elementary Differential Equations and Boundary Value Problems

MEM 355 Performance Enhancement of Dynamical Systems A First Control Problem - Cruise Control

Ratio-Product Type Exponential Estimator For Estimating Finite Population Mean Using Information On Auxiliary Attribute

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

CSE 245: Computer Aided Circuit Simulation and Verification

The Model and Numerical Analysis of Hardening Phenomena for Hot-work Tool Steel

Double Slits in Space and Time

I) Title: Rational Expectations and Adaptive Learning. II) Contents: Introduction to Adaptive Learning

Ministry of Education and Science of Ukraine National Technical University Ukraine "Igor Sikorsky Kiev Polytechnic Institute"

The transition:transversion rate ratio vs. the T-ratio.

Modelling of three dimensional liquid steel flow in continuous casting process

Decline Curves. Exponential decline (constant fractional decline) Harmonic decline, and Hyperbolic decline.

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

COMPUTATIONAL VISCOELASTICITY OF AGING MATERIALS

Instability Analysis of Laminated Composite Beams Subjected to Parametric Axial Load

On the Speed of Heat Wave. Mihály Makai

Let s look again at the first order linear differential equation we are attempting to solve, in its standard form:

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

Chapter 4 Longitudinal static stability and control Effect of acceleration (Lecture 15)

A HAMILTON-JACOBI TREATMENT OF DISSIPATIVE SYSTEMS

ELASTO-PLASTIC ANALYSIS OF STRUCTURES USING HEXAHEDRICAL ELEMENTS WITH EIGHT NODES AND ONE-POINT QUADRATURE

2.1. Differential Equations and Solutions #3, 4, 17, 20, 24, 35

4.3 Design of Sections for Flexure (Part II)

The model proposed by Vasicek in 1977 is a yield-based one-factor equilibrium model given by the dynamic

Impulsive Differential Equations. by using the Euler Method

Spring 2006 Process Dynamics, Operations, and Control Lesson 2: Mathematics Review

Equation For non-self Energizing Gasket

Wave Equation (2 Week)

Phys463.nb Conductivity. Another equivalent definition of the Fermi velocity is

H is equal to the surface current J S

Transfer function and the Laplace transformation

STRUCTURAL ANALYSIS OF MONO-SYMMETRIC PLATE GIRDERS IN COMPOSITE BRIDGES

7.4 QUANTUM MECHANICAL TREATMENT OF FLUCTUATIONS *

= x. I (x,y ) Example: Translation. Operations depend on pixel s Coordinates. Context free. Independent of pixel values. I(x,y) Forward mapping:

An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

Charging of capacitor through inductor and resistor

a dt a dt a dt dt If 1, then the poles in the transfer function are complex conjugates. Let s look at f t H t f s / s. So, for a 2 nd order system:

Applied Statistics and Probability for Engineers, 6 th edition October 17, 2016

A Condition for Stability in an SIR Age Structured Disease Model with Decreasing Survival Rate

1. Inverse Matrix 4[(3 7) (02)] 1[(0 7) (3 2)] Recall that the inverse of A is equal to:

General Article Application of differential equation in L-R and C-R circuit analysis by classical method. Abstract

CHAPTER CHAPTER14. Expectations: The Basic Tools. Prepared by: Fernando Quijano and Yvonn Quijano

fiziks Institute for NET/JRF, GATE, IIT JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES MATEMATICAL PHYSICS SOLUTIONS are

Midterm exam 2, April 7, 2009 (solutions)

The finite element models of thin-walled branched structures in heat transfer problems

International Conference on Energy and Environmental Protection (ICEEP 2016)

Control System Engineering (EE301T) Assignment: 2

Lecture 4: Laplace Transforms

whereby we can express the phase by any one of the formulas cos ( 3 whereby we can express the phase by any one of the formulas

A THREE COMPARTMENT MATHEMATICAL MODEL OF LIVER

UNSTEADY FLOW OF A FLUID PARTICLE SUSPENSION BETWEEN TWO PARALLEL PLATES SUDDENLY SET IN MOTION WITH SAME SPEED

CHAPTER. Linear Systems of Differential Equations. 6.1 Theory of Linear DE Systems. ! Nullcline Sketching. Equilibrium (unstable) at (0, 0)

CHAPTER 9 Compressible Flow

ON THE USP CALCULATION UNDER SOLVENCY II AND ITS APPROXIMATION WITH A CLOSED FORM FORMULA

Micromechanical Modeling of Concrete at Early Age

Estimation of Mean Time between Failures in Two Unit Parallel Repairable System

Single Electron Devices for Logic Applications

Lecture 2: Current in RC circuit D.K.Pandey

( ) 2! l p. Nonlinear Dynamics for Gear Fault Level. ( ) f ( x) ( ),! = sgn % " p. Open Access. Su Xunwen *,1, Liu Jinhao 1 and Wang Shaoping 2. !

Institute of Actuaries of India

A MATHEMATICAL MODEL FOR NATURAL COOLING OF A CUP OF TEA

B) 25y e. 5. Find the second partial f. 6. Find the second partials (including the mixed partials) of

Chapter 3: Fourier Representation of Signals and LTI Systems. Chih-Wei Liu

A STUDY OF MINDLIN PLATE FINITE ELEMENTS

Laplace Transforms recap for ccts

PFC Predictive Functional Control

G. =, etc.

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

LaPlace Transform in Circuit Analysis

EXERCISE - 01 CHECK YOUR GRASP

Lagrangian for RLC circuits using analogy with the classical mechanics concepts

CPSC 211 Data Structures & Implementations (c) Texas A&M University [ 259] B-Trees

A Numerical Simulation and New Traveling Wave Solutions of Convection-Diffusion Equation with Reaction

Lecture 1: Growth and decay of current in RL circuit. Growth of current in LR Circuit. D.K.Pandey

A Simple Procedure to Calculate the Control Limit of Z Chart

3.4 Repeated Roots; Reduction of Order

On General Solutions of First-Order Nonlinear Matrix and Scalar Ordinary Differential Equations

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

A UNIFIED APPROACH FOR FIRE RESISTANCE PREDICTION OF STEEL COLUMNS AND FRAMES

Integrity Control in Nested Certificates

Experimental and Computer Aided Study of Anisotropic Behavior of Material to Reduce the Metal Forming Defects

Combined Bending with Induced or Applied Torsion of FRP I-Section Beams

CONTROL SYSTEM DESIGN OF GUIDED MISSILE

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS I: Introduction and Linear Systems

Voltage v(z) ~ E(z)D. We can actually get to this wave behavior by using circuit theory, w/o going into details of the EM fields!

Chapter 12 Introduction To The Laplace Transform

Logistic equation of Human population growth (generalization to the case of reactive environment).

FWM in One-dimensional Nonlinear Photonic Crystal and Theoretical Investigation of Parametric Down Conversion Efficiency (Steady State Analysis)

3(8 ) (8 x x ) 3x x (8 )

Feedback Control and Synchronization of Chaos for the Coupled Dynamos Dynamical System *

Almost power law : Tempered power-law models (T-FADE)

14.02 Principles of Macroeconomics Problem Set 5 Fall 2005

ERROR ANALYSIS A.J. Pintar and D. Caspary Department of Chemical Engineering Michigan Technological University Houghton, MI September, 2012

The Markovian portfolio selection model with GARCH volatility dynamics

Title. Author(s)ANG, K. K.; THI, T. M.; HAI, L. V. Issue Date Doc URL. Type. Note. File Information.

S.Y. B.Sc. (IT) : Sem. III. Applied Mathematics. Q.1 Attempt the following (any THREE) [15]

GENERALIZATION OF NON-ITERATIVE NUMERICAL METHODS FOR DAMAGE-PLASTIC BEHAVIOUR MODELING

Transcription:

ISBN 978-93-84422-22-6 Procdings of 2015Inrnaional Confrnc on Innovaions in Civil and Srucural Enginring (ICICSE'15) Isanbul (Turky), Jun 3-4, 2015. 206-212 Failur Load of Plan Sl Frams Using h Yild Surfac Mhod Smail Boukloua 1, Mohamd Laid Samai 2, Abdlhadi Tkkouk 2 1 Civil nginring darmn, Univrsiy of Bordj Bou Arréridj, Algria 2 Darmn of Civil Enginring, Univrsiy of Consanin 1, Algria Absrac. In h rsn work, a nonlinar analysis mhod of lan sl frams using yild surfac mhod is roosd. Th yild funcion, considring h inracion of bnding momn and axial forc can b usd o drmin h laso-lasic siffnss marix of bam lmn usd for such srucural. Th collas load and h collas mchanism of lan sl frams ar drmind by a numrical rogram in Malab using h incrmnal dirc mhod (s by s) and h fini lmn for which h yild surfacs of Eurocod3 is adod. Analysis rsuls show ha h roosd mhod is saisfacory. Kywords: nonlinar analysis, sl fram, yilds surfac, collas loads, marial nonlinar. 1. Inroducion Whn h srucur is subjcd o loads ha xcd h roorional limi of h marial, h marial sars o yild, hus h abov mniond assumions bcom inadqua and canno rrsn h ral bhaviour of h srucur. In his cas lasic analysis is rquird. Plasic analysis mhods can b classifid in wo grous: disribud lasiciy mhods ha accoun for srading of lasic zons wihin h whol volum of h srucur (Plasic zon mhods) and lumd lasiciy mhods ha assum lasic zons o b formd wihin small ara sa h nds of fram mmbrs calld lasic hings, whil fram mmbrs xhibi lasic bhaviour bwn lasic hings (Plasic hing mhods). Marial yilding, h ffcs of gomrical non linariis, rsidual srss and h yild funcion includs h ffc of h srss comonns acing in h sysm o rdic h yilding of h marial ar major aramrs ha conrol h load-carrying caaciy of h srucur, and hav bcom a ar of many naional Sandards and Cods (Eurocod3, AISC, Briish Sandards, c). Morovr, fas-sd rsonal comurs dvlod in h las 20 yars mad h us of nonlinar analysis rocdurs mor availabl for racical uross. Howvr, in his cas, h discussion will b limid o h marial nonlinar, h ffcs of rsidual srss and h lasiciy is suosd o b concnrad only in h cross scion of h nds of h bams (lasic hing mhod) and i is in lasic sa by h combinaion of srss ha saisfis h yilding condiion, inracion of h bnding momn wih axial forc. Digial rogram in Malab has bn dvlod for h load facor calculaion by adoing h aroach of yild surfac. This rogram uss h fini lmn mhod o succssiv linar analyzs and basd on h s by s mhod. Corrsonding auhor, Docora sudn, E-mail: smail_ing097@yahoo.fr h://dx.doi.org/10.17758/ur.u0615319 206

2. Yild Surfac Th variaion of h bnding momn wih axial forc in a cross scion can b lod in rms of h dimnsionlss quaniis N/N and M/M. Th rsuling curv is calld h yild surfac bcaus any oin on h yild surfac rrsns a sa of h fully yildd cross scion. Th I or H-shad scions ar ofn usd in sl frams, for which h yild surfacs of Eurocod3 [1] is adod in rsn work. Th quaions ar rsnd bllows: - Yild surfac of Eurocod3 [1] Figur 1 : Wih: M 0.9M 1 M 1 = 0.9M 1-0.5a a =min [A w /A, 0.5] A: ara of scion. A w : ara of wb. for N N a/2 N N (2) 1- for a/2 0.8N N (1) M/M P 1.0 Fig. 1: Yild surfacs of sl I-H scions of Eurocod3 In which, n=n/0.8n is raio of h axial forc ovr h squash load, m=m/0.9m is h raios of h majoraxis momns o h corrsonding lasic momns, and h numbrs 0.8 and 0.9 in h dnominaor accoun for rsidual srsss. 2.1. Normaliy Rul Th yild surfacs dscribd for various cross-scional shas can b rsnd using a yild funcion Ф such ha for a scion in a fully yildd sa undr forc inracion Ф =0. Whn h ffcs of bnding momn and axial forc ar akn ino accoun on h yild surfac, h associad gnralizd srains ar h roaion and h axial dislacmn of scion. Th normaliy rul was originally roosd by Von Miss in 1928, i may b alid for his cas as follows: N/N P h://dx.doi.org/10.17758/ur.u0615319 207

U { } N M Or, symbolically: f { dp } (4) Whr d } rrsn h vcor of h lasic dformaion incrmns, λ is h lasic dformaion magniud, { P f is a gradin vcor a a oin of h yild surfac Ф. Whn h lasic loading occurring, h oin forc is on h yild surfac (or subsqunc yild surfac) Ф=0. In aking h drivaiv of his rlaionshi, w obain: (3) N M 0 N M Whr h arial drivaivs mus b akn a h original sa of srss rsulan. Equaion (5) can b wrin in vcorial form as { P} 0 f (6) Whr {P} rrsn h vcor for h incrmns of srss rsulans. Th orhogonal condiion can b alid o h rlaionshi bwn h incrmns of srss rsulans and lasic dformaion as imlid by Pagr s [2] samn ha for lasic-rfcly lasic marial, h srss incrmn dos no work on h incrmn of lasic srain. Whn alying o fram mmbrs, his samn mans ha { d}{ P} 0 (7) For marials in h lasic sa, h lasic flow always occurs in associaion wih a dissiaion of mchanical nrgy. Thus, for an incrmn of lasic dformaion{ d}, h dissiaiv nrgy ΔW is always osiiv and is givn by: W {P} { d} {P} f 0 (8) 2.2. Elasolasic Siffnss Marix For a scion in lasic sa, h incrmnal dformaion vcor, { d}, consiss of boh lasic and lasic dislacmns, dnding on which forc comonns ar aciv in h yild funcion. Hnc: { d} { d} { d} (9) Whr h incrmnal lasic dislacmn vcor K { d} d f K (5) { d} is rlad o h incrmnal forc vcor by: { P} (10) Whr K is h lasic siffnss marix. Using Equaion (10) in Equaion (6), h lasic mulilir λ can b found o b: h://dx.doi.org/10.17758/ur.u0615319 208

f K d f K f Subsiuing Equaion (13) ino Equaion (12), h lasolasic siffnss marix, K K { d} (11), can b found: { P} (12) Whr h lasolasic siffnss marix is: K K K f f K f K f (13) Equaion (12) is a gnral xrssion for a yildd bam lmn. Sinc a bam lmn may b subjcd o diffrn combinaions of yilding sas a is nds, h form of K varis according o h sa of yilding and h yild funcion adod for lasic analysis. 3. Numrical Vrificaion 3.1. Vogl Poral Fram Th oral fram shown in Figur2 was analysd numrically in 1985 by Vogl [3], and his fram has bn usd by svral rsarchrs (Chn 1993[4], Chn and Kim 1997 [5], Kim and Lu 1992[6], Toma and Chn 1992[7]) as a bnchmark soluion for including marial non-linariis including rsidual srsss, gradual yilding and full lasiciy. Th fram siz, marial roris and load informaion ar illusrad in Figur 2, and h fram mmbr sizs ar lisd in Tabl 1. Th horizonal dislacmn of righ ur cornr (nod A) vrsus load facor curv by h roosd aroach is comard wih h lasic zon mhod and wih h lasic hing mhod of Vogl (1985) wih h lasic hing mhod in Figur3. Th ulima load facor obaind by h mhod roosd is 0.9192 whras ha by Vogl s lasic hing analysis is λ=1.017 and Vogl s lasic zon analysis is λ=1.02. V=2800kN V=2800kN H=35kN HEA340 A Δ HEB300 HEB300 5m 4m Fig. 2: Gomric configuraions and loading arn of Vogl s fram h://dx.doi.org/10.17758/ur.u0615319 209

1.2 1 0.9192 1.017 1.02 0.8 Load facor 0.6 0.4 0.2 Vogl s lasic zon mhod Vogl s lasic hing mhod Prsn analysis (yild surfac Eurocod3) 0 0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 Laral dislacmns(m) Fig. 3: Load-dislacmn curv a h o of Vogl s oral fram TABLE I: Mmbr sizs and scional roris of h Vogl oral fram Scion A(mm 2 ) I(x10 6 mm4) h(mm) b(mm) roris W f W l (x10 3 mm 3 HEA340 330 300 9.5 16.5 13300 276.9 1850 HEB300 300 300 11.0 19.0 14900 251.7 1869 3.2. Six Sory Fram Th fram siz and load informaion ar illusrad in Figur 4 and h fram mmbr sizs ar lisd in Tabl 2. Th marial lasic modulus E of sl is 206kN/mm2 and h yild srngh f y is 235N/mm2. Th horizonal dislacmn of righ-ur cornr (Nod A) vrsus load facor curv by h yild surfac modl roosd in his ar is rsnd in Figur 3. Th ulima load facor λ obaind by h mhod roosd is 1.2888. 1,4 1.2888 1,2 1 Load facor 0,8 0,6 0,4 0,2 0 0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 Laral dislacmns(m) a nod A Fig. 3: Load-dislacmn curv a h o of six sory fram h://dx.doi.org/10.17758/ur.u0615319 210

TABLE II: Mmbr sizs and scional roris of h six-sory fram Scion A(mm 2 ) I(x10 6 mm4) h(mm) b(mm) roris W (mm) f (mm) W l (x10 3 mm 3 HEA340 330 300 9.5 16.5 13300 276.9 1850 HEB160 160 160 8.0 13.0 5430 24.92 354 HEB200 200 200 9.0 15.0 7810 56.96 643 220 220 9.5 16 9100 80.91 827 HEB240 240 240 10.0 17 10600 112.6 1053 HEB260 260 260 10.0 17.5 11800 149.2 1283 HEB300 300 300 11.0 19.0 14900 251.7 1869 IPE240 240 120 6.2 9.8 3910 38.92 367 IPE300 300 150 7.1 10.7 5380 83.56 628 IPE330 330 160 7.5 11.5 6260 117.7 804 IPE360 360 170 8.0 12.7 7270 162.7 1019 IPE400 400 180 8.6 13.5 8450 231.3 1307 P 2 P 2=63.4kN 2P 2 P 2 P 2 IPE240 HEB200 HEB160 A 19 IPE300 IPE300 HEB200 HEB240 HEB160 IPE330 IPE360 HEB240 HEB260 IPE400 H 2 =10.23kN HEB260 6m 6m =20.44kN Fig. 4: Six-sory sl fram. h://dx.doi.org/10.17758/ur.u0615319 211

4. Conclusion An aroach for nonlinar analysis of sl frams using yild surfac mhod is roosd in his ar. This aroach us h yild funcion includs h ffc of h axial forc and bnding momn acing in h sysm o rdic h yilding of h marial. This aroach also considrs h influncs of h marial nonlinar including h rsidual srss. Th numrical rsuls show ha h roosd aroach is saisfacorily, and is suiabl for h nonlinar analysis of sl frams. 5. Rfrncs [1] Eurocod3, Dsign of Sl Srucurs (2003), Euroan Commi for Sandardisaion. [2] W. Pragr,An inroducion o lasiciy, London. Addison-Wsly Pub. Co, Inc (1959). [3] U.Vogl, Calibraing frams. Dr Sahlbau 1985; 10: 296-301. [4] S. Toma, WF. Chn, Euroan calibraion frams for scond-ordr inlasic analysis.enginring Srucurs 1992; 14(1): 7-14. h://dx.doi.org/10.1016/0141-0296(92)90003-9 [5] ZH.Zhou,SL. Chan, Elasolasic and larg dflcion analysis of sl frams by on lmnr mmbr. I: On hing along mmbr, Journal of Srucural Enginring, ASCE 2004;130(4): 538-544. h://dx.doi.org/10.1061/(asce)0733-9445(2004)130:4(538) [6] W. F. Chn, Advancd Analysis in Sl Frams, (1993) CRC Prss, Boca Raon, FL. [7] W. F. Chn ands-e. Kim, LRFD Sl Dsign Using Advancd Analysis,(1997) CRC Prss, Boca Raon, FL. h://dx.doi.org/10.17758/ur.u0615319 212