Progressive failure of masonry shear walls a distinct element approach *

Similar documents
Temperature Memory Effect in Amorphous Shape Memory Polymers. Kai Yu 1, H. Jerry Qi 1, *

Linear Approximating to Integer Addition

Research on structural optimization design for shield beam of hydraulic support. based on response surface method

European Journal of Mathematics and Computer Science Vol. 5 No. 2, 2018 ISSN

On the energy of complement of regular line graphs

On a Truncated Erlang Queuing System. with Bulk Arrivals, Balking and Reneging

European Journal of Mathematics and Computer Science Vol. 5 No. 2, 2018 ISSN

Simple Linear Regression Analysis

Trignometric Inequations and Fuzzy Information Theory

Solving Constrained Flow-Shop Scheduling. Problems with Three Machines

International Journal of Pure and Applied Sciences and Technology

Manipulator Dynamics. Amirkabir University of Technology Computer Engineering & Information Technology Department

Reaction Time VS. Drug Percentage Subject Amount of Drug Times % Reaction Time in Seconds 1 Mary John Carl Sara William 5 4

Dynamic Analysis of Axially Beam on Visco - Elastic Foundation with Elastic Supports under Moving Load

T-DOF PID Controller Design using Characteristic Ratio Assignment Method for Quadruple Tank Process

Nargozy T. Danayev*, Darkhan Zh. Akhmed-Zaki* THE USAGE OF MATHEMATICAL MLT MODEL FOR THE CALCULATION OF THERMAL FILTRATION

Multiple Choice Test. Chapter Adequacy of Models for Regression

A Method for Damping Estimation Based On Least Square Fit

Functions of Random Variables

CS473-Algorithms I. Lecture 12b. Dynamic Tables. CS 473 Lecture X 1

EFFECT OF MODAL TRUNCATION IN MULTIPLE SUPPORT RESPONSE SPECTRUM ANALYSIS OF BRIDGES.

Chapter 2 - Free Vibration of Multi-Degree-of-Freedom Systems - II

Beam Warming Second-Order Upwind Method

REVIEW OF SIMPLE LINEAR REGRESSION SIMPLE LINEAR REGRESSION

Simple Linear Regression

SURFACE STRESS EFFECT IN THIN FILMS WITH NANOSCALE ROUGHNESS

2C09 Design for seismic and climate changes

Estimation of Stress- Strength Reliability model using finite mixture of exponential distributions

Fourth Order Four-Stage Diagonally Implicit Runge-Kutta Method for Linear Ordinary Differential Equations ABSTRACT INTRODUCTION

ON THE MOTION OF PLANAR BARS SYSTEMS WITH CLEARANCES IN JOINTS

A Result of Convergence about Weighted Sum for Exchangeable Random Variable Sequence in the Errors-in-Variables Model

End of Finite Volume Methods Cartesian grids. Solution of the Navier-Stokes Equations. REVIEW Lecture 17: Higher order (interpolation) schemes

C-1: Aerodynamics of Airfoils 1 C-2: Aerodynamics of Airfoils 2 C-3: Panel Methods C-4: Thin Airfoil Theory

KR20 & Coefficient Alpha Their equivalence for binary scored items

Analysis of Lagrange Interpolation Formula

Some distances and sequences in a weighted graph

PGE 310: Formulation and Solution in Geosystems Engineering. Dr. Balhoff. Interpolation

System Reliability-Based Design Optimization Using the MPP-Based Dimension Reduction Method

The Mathematical Appendix

Computational Modeling for Acoustic Wave Propagation in a Layered Media

2006 Jamie Trahan, Autar Kaw, Kevin Martin University of South Florida United States of America

8 The independence problem

ROOT-LOCUS ANALYSIS. Lecture 11: Root Locus Plot. Consider a general feedback control system with a variable gain K. Y ( s ) ( ) K

ANOVA with Summary Statistics: A STATA Macro

Analysis of von Kármán plates using a BEM formulation

DATE: 21 September, 1999 TO: Jim Russell FROM: Peter Tkacik RE: Analysis of wide ply tube winding as compared to Konva Kore CC: Larry McMillan

Part 4b Asymptotic Results for MRR2 using PRESS. Recall that the PRESS statistic is a special type of cross validation procedure (see Allen (1971))

2.28 The Wall Street Journal is probably referring to the average number of cubes used per glass measured for some population that they have chosen.

Comparison of Analytical and Numerical Results in Modal Analysis of Multispan Continuous Beams with LS-DYNA

A Remark on the Uniform Convergence of Some Sequences of Functions

A note on testing the covariance matrix for large dimension

FRP. ***Download is permitted to HPC-Lab members and for educational purposes only*** (combined finite element /discrete element FRP

Centroids & Moments of Inertia of Beam Sections

Cubic Nonpolynomial Spline Approach to the Solution of a Second Order Two-Point Boundary Value Problem

The numerical simulation of compressible flow in a Shubin nozzle using schemes of Bean-Warming and flux vector splitting

Quiz 1- Linear Regression Analysis (Based on Lectures 1-14)

Acoustics Field and Active Structural Acoustic Control Modeling in ANSYS

A New Family of Transformations for Lifetime Data

Wu-Hausman Test: But if X and ε are independent, βˆ. ECON 324 Page 1

best estimate (mean) for X uncertainty or error in the measurement (systematic, random or statistical) best

Block-Based Compact Thermal Modeling of Semiconductor Integrated Circuits

A Conventional Approach for the Solution of the Fifth Order Boundary Value Problems Using Sixth Degree Spline Functions

Aitken delta-squared generalized Juncgk-type iterative procedure

Comparison of Dual to Ratio-Cum-Product Estimators of Population Mean

Newton s Power Flow algorithm

Further Results on Pair Sum Labeling of Trees

Module 7: Probability and Statistics

DEVELOPMENT OF FLUID-SEDIMENT-SEABED INTERACTION MODEL AND ITS APPLICATION. Tomoaki Nakamura 1 and Norimi Mizutani 1

Introduction to local (nonparametric) density estimation. methods

Numerical simulation of sand mould manufacture for Lost Foam Casting process

On the Interval Zoro Symmetric Single Step. Procedure IZSS1-5D for the Simultaneous. Bounding of Real Polynomial Zeros

ANALYSIS ON THE NATURE OF THE BASIC EQUATIONS IN SYNERGETIC INTER-REPRESENTATION NETWORK

10.2 Series. , we get. which is called an infinite series ( or just a series) and is denoted, for short, by the symbol. i i n

u 1 Figure 1 3D Solid Finite Elements

PROJECTION PROBLEM FOR REGULAR POLYGONS

HOOKE'S LAW. THE RATE OR SPRING CONSTANT k.

p y A modified p-y curve method considering rotation of soil resistance

Physics 114 Exam 2 Fall Name:

Ordinary Least Squares Regression. Simple Regression. Algebra and Assumptions.

Unimodality Tests for Global Optimization of Single Variable Functions Using Statistical Methods

Collapsing to Sample and Remainder Means. Ed Stanek. In order to collapse the expanded random variables to weighted sample and remainder

4. Standard Regression Model and Spatial Dependence Tests

ECE 595, Section 10 Numerical Simulations Lecture 19: FEM for Electronic Transport. Prof. Peter Bermel February 22, 2013

1. Linear second-order circuits

means the first term, a2 means the term, etc. Infinite Sequences: follow the same pattern forever.

ECE606: Solid State Devices Lecture 13 Solutions of the Continuity Eqs. Analytical & Numerical

Application of Laplace Adomian Padé approximant to solve exponential stretching sheet problem in fluid mechanics

Centers of Gravity - Centroids

The internal structure of natural numbers, one method for the definition of large prime numbers, and a factorization test

Problem Set 3: Model Solutions

Lecture 07: Poles and Zeros

Third handout: On the Gini Index

Arithmetic Mean and Geometric Mean

Application of Calibration Approach for Regression Coefficient Estimation under Two-stage Sampling Design

Research Article A New Iterative Method for Common Fixed Points of a Finite Family of Nonexpansive Mappings

Reflection of Plane Waves from a Free Surface of an Initially Stressed Transversely Isotropic Dissipative Medium

Ahmed Elgamal. MDOF Systems & Modal Analysis

Reliability and Cost Analysis of a Series System Model Using Fuzzy Parametric Geometric Programming

G S Power Flow Solution

A Note on Ratio Estimators in two Stage Sampling

Transcription:

Joural of Appled Mathematc ad Phyc, 2016, *, *-* http://www.crp.org/joural/jamp ISSN Ole: 2327-4379 ISSN Prt: 2327-4352 Progreve falure of maory hear wall a dtct elemet approach * (Afflato): School of Cvl Egeerg ad Surveyg, Uverty of Souther Queelad, Brbae, Autrala Emal: ya.zhuge@uq.edu.au How to cte th paper: Author 1, Author 2 ad Author 3 (2016) Paper Ttle. ******, *, *- *. http://dx.do.org/10.4236/***.2016.***** Receved: **** **, *** Accepted: **** **, *** Publhed: **** **, *** Copyrght 2016 by author() ad Scetfc Reearch Publhg Ic. Th work lceed uder the Creatve Commo Attrbuto Iteratoal Lcee (CC BY 4.0). http://creatvecommo.org/lcee/by/4.0/ Ope Acce Abtract I th paper a umercal model ha bee developed to tudy the behavor of maory wall uder plae hear. Maory ot a mple materal, the fluece of mortar jot a a plae of weake a gfcat feature ad th make the umercal modellg of maory very dffcult epecally whe dyamc (emc) aaly volved. Although fte elemet method a very powerful umercal method for the aaly of tructure ad t ha bee ued to model maory, mcro modellg ature ha made the aaly complex. It wa foud that a dtct elemet method (DEM) could be ued to mulate the behavor of dcotute. I DEM, a old repreeted a a aembly of dcrete block. Jot are modelled a terface betwee dtct bode. It a dyamc proce ad could be appled to mulate the dyamc repoe of maory. I the paper, the umercal oluto obtaed from DEM are valdated by comparg the reult wth thoe obtaed from extg expermet ad fte elemet modellg. Keyword Dtct elemet method (DEM), Ureforced maory (URM), Shear wall, Progreve falure, Numercal modellg 1. Itroducto Fte elemet method a very powerful umercal method for the aaly of tructure ad t ha bee ued to model the behavour of ureforced maory (URM) for more tha 20 year. For tatc aaly, ether a oe-phae (macro modellg of compote) materal model or a two-phae (mcro modellg of brck ad mortar) materal model ha bee ued. Maory ot a mple materal, t compoed of two materal a geometrc array - a aemblage of brck et a mortar matrx. The fluece of mortar jot ad bod a a plae of weake a gfcat feature whch ot preet cocrete ad th make the umercal modellg of maory very dffcult epecally whe emc aaly volved. To mplfy the problem, a homogeeou oe-phae fte elemet model whch ** Specal decrpto of the ttle. (dpeable) DOI: 10.4236/***.2016.***** **** **, 2016

corporate the olear materal behavour ha bee adopted by may vetgator, epecally for dyamc aaly [1], [2]. I the model, ome vetgator adopted a homogeato procedure to corporate the effect of mortar jot [3]. However, the weak jot behavour, whch domate the local falure of maory may cae, could ot be mulated well by a cotuum model. Due to the ature of the problem ad mproved computatoal trategy, a more complcated two-phae fte elemet model ha bee ued by may reearcher recet year [4], [5]. I the model, the brck are modelled a cotuum elemet ad the jot are modelled by le terface elemet. However, uch model udoubtedly made aaly more complex ad t early mpoble to apply thee model to the dyamc/emc aaly of maory. There are other umercal model avalable whch are maly ued for Rock Mechac, amely the Boudary Elemet method (BEM) ad the Dtct Elemet Method (DEM), although they are ot commoly ued for tructural aaly. DEM wa frt troduced by [6] ad ha bee progrevely developed over the pat two decade. The method wa mplemeted the computer program [7] ad t major applcato ha bee the feld of rock mechac ad mg egeerg. However, t teretg to otce the mlarte the geeral materal behavour of rock ad maory, a both could be clafed a o-homogeeou, dcotuou materal. Therefore, t poble to apply the ewly developed aalytcal techque for rock egeerg to mcro modellg of maory ad DEM foud t uefule here. The frt applcato of the DEM to model maory, foud the lterature, wa by [8], where DEM wa appled to vetgate the tablty of freetadg clacal colum uder earthquake exctato. A mplfed 2D model wa adopted ad the deformablty of the block wa eglected. Ther tude proved that the DEM very effcet mulatg the progreve collape of blocky type tructure. More recetly, [9] appled DEM to tudy materal parameter detfcato for maory cottutve model. The tructural repoe of low bod tregth maory wall pael wa mulated by DEM. Materal parameter for DEM model were tued ug a optmzato proce order to acheve mlar repoe to thoe obtaed the laboratory. I th paper, a umercal model baed o the DEM ha bee appled to mulate the repoe of a ureforced maory hear wall pael wth ad wthout a opeg, where the expermetal tetg reult are avalable. The DEM formato troduced frt, followed by the materal model for brck ad jot a well a the falure crtera. The aaly performed ad the reult betwee the dtct elemet model ad expermet ad fte elemet model are compared ad dcued. 2. 2D dtct elemet modellg of URM The Dtct Elemet Method ha bee progrevely developed over the pat three decade. It wa frt troduced by [6] to mulate progreve movemet blocky rock ytem ad the model ha bee mplemeted to a computer program UDEC ce the. DEM mulate the repoe of dcotuou meda ubjected to ether tatc or dyamc loadg. 2

Author, Author I the DEM method, a old repreeted a a aembly of dcrete block. Jot are modelled a terface betwee dtct bode. The cotact force ad dplacemet at the terface of a treed aembly of block are foud through a ere of calculato, whch trace the movemet of the block [7]. At all the cotact, ether rgd or deformable block are coected by prg lke jot wth ormal ad hear tffe k ad k, repectvely (Fg 1), that repreet the force-dplacemet relatohp of the jot. Smlar to the Fte Elemet Method (FEM), the ukow the DEM are alo the odal dplacemet ad rotato of the block. However, ulke FEM, whch baed o a cotuum mechac formulato, DEM provde the capacty to repreet the moto of multple, terectg dcotute explctly. DEM ue a explct tme-marchg cheme to olve the equato of moto drectly. The calculato alterate betwee applcato of a force-dplacemet law at all jot ad Newto ecod law at all block. The cottutve law of jot ued to fd cotact force/tree from kow dplacemet. Newto ecod law the gve the moto of the block reultg from the kow force/tree actg o them. Fgure 1. Cotuum ad dcotuum elemet DEM. 2.1. Cottutve law ad falure crtero of jot I geeral, all tele ad hear damage have bee aumed to take place the relatvely weak jot whch ha bee proved to be true for mot maory hear wall pael. I the model, the mortar jot are repreeted umercally a a cotact urface formed betwee two block edge. The cottutve law appled to the cotact are: k u (1) k u (2) where k ad k are the ormal ad hear tffe of the cotact, ad are the effectve ormal ad hear tre cremet, ad u ad u are the ormal ad hear dplacemet cremet. Stree calculated at grd pot located alog cotact are ubmtted to the elected falure crtero. For the propoed model, there a lmtg tele tregth f t for the jot. If the tele tregth exceeded, the = 0. For hear, the Coulomb frcto formulated: C ta max (3) 3

where C the coheo ad the frcto agle. For the jot falure model, the dlato retrcted a follow: f max, The = 0 (4) f max ad c u u, the = 0 (5) where u c crtcal hear dplacemet. 2.2. Cottutve law ad falure crtero of block The ormal ad hear tree calculated Equato (1) ad (2) are the ued a the boudary tree to ether rgd or deformable block. If the block are deformable, they wll be terally dcreted to fte dfferece tragular elemet frt. The ue of tragular elemet elmate the problem of hourgla deformato that may occur wth cotat-tra fte dfferece quadrlateral. The vertce of the tragular elemet are grd pot (Fg 1). Newto ecod law appled at each grd pot a follow: ( t) t j jd F ( ) u g (6) t m where the urface eclog the ma m lumped at the grd pot, j the ut ormal (t) to, F the reultat of all exteral force appled to the grd pot ad g the gravtatoal accelerato. A et odal force vector, F (t) calculated at each grd pot. Th vector clude (g ) cotrbuto from appled load ad from body force due to gravty. Gravty force, F are computed from: ( g) F g m (7) where m g the lumped gravtatoal ma at the grd pot, defed a the um of oe-thrd of the mae of tragle coected to the grd pot. If the body at equlbrum, F (t) o the ode wll be zero; otherwe, the ode wll be accelerated accordg to the fte dfferece form of Newto ecod law of moto ad followg the cetral dfferece tegrato cheme: ( t) ( tt/ 2) ( tt/ 2) F u u t (8) m Where u the velocty, m the ma, ad t the tme. Durg each tme tep, tra ad rotato are related to odal dplacemet the uual form: 1 j ( u, j u j, ) 2 (9) 1 j ( u, j u j, ) 2 The the elected cottutve law for the block ued a cremetal form to determe tree at each grd pot. A olear model ca be ealy corporated a the oluto cheme ued by DEM the explct tme marchg cheme. To mply the problem, the Mohr-Coulomb falure crtero wth teo cut-off adopted the model. The ew poto of the block duce ew codto at block boudare ad thu ew cotact force. 3. Comparo of umercal ad expermetal reult g 4

Author, Author The umercal example preeted th paper were performed wth the umercal code UDEC [7]. A ere of URM hear wall were teted by [4] whch cluded reult for mcromodellg. Thee tetg reult are ued to valdate the umercal model. The hear wall pael ha a wdth/heght rato of 990 mm x1000mm, bult up wth 18 coure wth ad wthout a cetral opeg. The wall wa made of wre-cut old clay brck wth dmeo 210 mm x 52 mm x 100 mm ad 10 mm thck mortar. The wall wa ubjected to a vertcal compreve tre m = 0.3 MPa. Durg tetg, a horzotal dplacemet wa mootocally creaed utl falure. The expermetal crack patter ad FEM reult of the fal falure patter of the wall are how Fgure 2(a) ad 2(b) repectvely [4]. A the load creaed, a tepped dagoal crack through head ad bed jot wa formed, the crack progree the drecto of the upport ad fally a collape mecham formed (Fgure 2). The tele ad hear crack developmet ad the progreve falure of the wall modelled by DEM are dcated Fgure 3. It ca be ee that the behavour of the wall well captured by the propoed model. (a) (b) Fgure 2. (a). Expermetal crack at a dplacemet of 2 mm; (b). FEM model of the fal falure patter [4]. JOB TITLE : Maory w all; UDEC (Vero 3.10) LEGEND 1.000 2-Apr-02 17:10 cycle 1459135 block plot 0.800 0.600 0.400 0.200 0.000 (a) (b) 0.000 0.200 0.400 0.600 0.800 1.000 Fgure 3. DEM Model: (a). Crack at a dplacemet of 2 mm; (b). Fal falure patter. The tetg reult of the wall wth the ame geometry ad loadg codto, but wth a cetral opeg how Fgure 4(a). The cetral opeg of the wall produce two mall weak per. The collape mecham wth falure of the compreed toe, located at the bottom ad top of the wall ad at the bottom ad top of the mall per. The wall behave mlarly to four rgd coected by the hge how. The DEM reult how Fgure 5

4(b) clearly dcated the two compreve trut have bee formed ad the locato of the maxmum compreve tree agreed well wth tetg reult. (a) (b) Fgure 4. (a). Collape mecham wth four hged rgd block [4]; (b). DEM reult of mmum prcpal tree at peak load. 4. Cocluo Maory ot a mple materal, the fluece of mortar jot a a plae of weake a gfcat feature. Due to the ature of the problem, two-phae mcro-modellg of brck ad mortar requred. However, th make the umercal modellg of maory very complex. Th paper dcued a alteratve ad mple way of modellg maory, whch ug the Dtct Elemet Method. The prelmary reult obtaed the curret reearch have how the great potetal of the method, epecally for dyamc aaly. Referece [1] Vrataou, V. (1991) Determato of the behavour factor for brck maory pael ubjected to earthquake acto. Proc. of the Iteratoal Cof. o Sol Dyamc ad Earthquake Egeerg, Germay: 23-26. [2] Zhuge, Y., Thambratam, D. P. ad Corderoy, J. (1998) Nolear dyamc aaly of ureforced maory. Joural of Structural Egeerg, ASCE, 124(3), 270-277. [3] Gambarotta, L. ad Lagomaro, S. (1997) Damage model for the emc repoe of brck maory hear wall. Part II: the cotuum model ad t applcato. Earthquake Egeerg & Structural Dyamc, 26, 441-462. [4] Loureco, P. B. (1996) Computatoal tratege for maory tructure, PhD The, Delft Uverty of Techology, The Netherlad. [5] Wag, C., Forth, J.P., Nkta, N. ad Sarho, V. (2016) Retrofttg of maory wall by ug a mortar jot techque: expermet ad umercal valdato. Egeerg Structure, 117, 58-70. [6] Cudall, P. A. (1971) A Computer model for mulatg progreve large cale movemet blocky rock ytem. Proceedg of the Sym. of the Iteratoal Socety for Rock mechac, Nacy, Frac, Vol. 1, II-8: 11-18. [7] ITASCA Coultg Group (2000) Uveral Dtct Elemet Code, ITASCA coultg Group, Ic., Meapol, Meota, USA. [8] Pychar, I., Papatamatou, D. ad Alexadr, A. (2000) Parametrc vetgato of the tablty of clacal colum uder harmoc ad earthquake exctato. Eart. Eg. Ad Struct Dyamc, 29, 1093-1109. [9] Sarho, V. ad Sheg, Y. (2014). Idetfcato of materal parameter for low bod tregth maory. Egeerg Structure, 60, 100-110. 6