Fabric Evolution and Its Effect on Strain Localization in Sand

Size: px
Start display at page:

Download "Fabric Evolution and Its Effect on Strain Localization in Sand"

Transcription

1 Fabric Evolution and Its Effect on Strain Localization in Sand Ziwei Gao and Jidong Zao Abstract Fabric anisotropy affects importantly te overall beaviour of sand including its strengt and deformation caracteristics. Wile bot experimental and numerical evidence indicates tat soil fabric evolves steadily wit te applied stress/strain, ow evolving fabric influences te initiation and development of sear band in sand remains an intriguing question to be fully addressed. In tis paper, we present a numerical study on strain localization in sand, igligting te special role played by soil fabric and its evolution. In particular, a critical state sand plasticity model accounting for te effect of fabric and its evolution is used in te finite element analysis of plane strain compression tests. It is found tat te initiation of sear band is controlled by te initial fabric, wile te development of sear band is governed by two competing pysical mecanisms, namely, te structural constraint and te evolution of fabric. Te evolution of fabric generally makes te sand response more coaxial wit te applied load, wile te structural constraint induced by te sample ends leads to more inomogeneous deformation witin te sand sample wen te initial fabric is non-coaxial wit te applied stress. In te case of smoot boundary condition, structural constraint dominates over te fabric evolution and leads to te formation of a single sear band. Wen te boundary condition is roug, te structural constraint may play a comparable role wit fabric evolution, wic leads to symmetric cross-sape sear bands. If te fabric is proibited from evolving in te latter case, a cross-sape sear band pattern is found wit te one initiated first by te structural constraint dominant over te second one. Z. Gao ( ) Scool of Engineering, University of Glasgow, Rankine Building, Oakfield Avenue, Glasgow G12 8LT, UK Ziwei.Gao@glasgow.ac.uk J. Zao Department of Civil and Environmental Engineering, Hong Kong University of Science and Tecnology, Clearwater Bay, Kowloon, Hong Kong Springer International Publising Switzerland 2015 K.-T. Cau and J. Zao (eds.), Bifurcation and Degradation of Geomaterials in te New Millennium, Springer Series in Geomecanics and Geoengineering, DOI / _4 21

2 22 Z. Gao and J. Zao 1 Introduction Strain localization is frequently observed in sand and is considered an important precursor of te failure of soil and relevant geostructures. Numerous investigations ave been carried out on strain localization in sand but relatively less attention as been paid to te correlation of strain localization wit te presence of an evolving fabric (e.g., Borja et al. 2013). Fabric anisotropy as been widely regarded to affect te key beaviour of sand including dilatancy, liquefaction and critical state. Based on plain strain compression tests, Tatsuoka et al. (1990) found tat te sear band development in sand was dependent on te initial bedding plane orientation, or te fabric, of te sample. Meanwile, micromecanical studies indicate tat fabric evolves steadily wit deformation and te evolving fabric exibits unique caracteristics wen a sand sample reaces liquefaction, critical state and oter states (Guo and Zao 2013; Zao and Guo 2013a, b). As an important type of instabilities, strain localization as been commonly observed in sand. Te correlation between fabric and fabric evolution is an interesting topic but remains less explored. Inerent fabric anisotropy as recently been considered in a ypoplastic model by Tejcman et al. (2007) in te simulation of sear band development in sand but te interplay between fabric evolution and te development of sear band as not been properly considered. In tis study, we employ an anisotropic sand model developed recently by te autors (Gao et al. 2014) and finite element metod to investigate strain localization in sand under plane strain compression werein te special role of fabric and its evolution is igligted. 2 Finite Element Analysis of Strain Localization in Plane Stain Compression Te constitutive model used in tis study was developed by Gao et al. (2014). It as been implemented in te finite element package ABAQUS troug te user-material interface. Te test data presented by Tatsuoka et al. (1990) will be employed to bencmark te model simulations. Te model parameters ave been calibrated based on te plane strain test results on Toyoura sand. Te initial degree of anisotropy F 0 is set to be Te model simulations for single element tests can be found in Gao et al. (2014) and Gao and Zao (2013). In te finite element analysis, te sample setup is te same as tat in Tatsuoka et al. (1990) (Fig. 1). Uniform 4-noded plane strain elements of 2.5 mm 2.5 mm in size are used and constant confining pressure is applied in te orizontal direction of te sample. Vertical displacement is applied to te top end of te sample by increment to ensure quasi-static loading. Bot smoot (te top and bottom ends are free to move orizontally) and roug (te orizontal displacements of te top and bottom ends are restricted) boundary conditions wit and witout fabric evolution are considered and te initial void ratio distribution is assumed to be uniform trougout te entire sample.

3 Fabric Evolution and Its Effect on Strain 23 =10.5cm Bedding plane y z w=4cm x α Fig. 1 Sample setup, element size, orientation of bedding plane and te reference coordinate system (a) (b) Type-b sear band Fig. 2 Predicted strain localization for te cases wit smoot boundary condition and a evolving fabric; b constant fabric at vertical strain ε = / = 12 % Figure 2 sows te predicted sear band pattern in smoot boundary condition cases wit initial void ratio e 0 = 0.7, confining pressure σ c = 400 kpa and bedding plane orientation α = 45. It can be seen tat te simulated sear band and bedding plane lie on te same side (left side in te upper part of te sample for te present study) of te major principal stress direction (Fig. 2a) and Tatsuoka et al. (1990) terms similar sear band pattern tey observed in laboratory tests as Type-b sear band. Indeed, Type-b sear band is te most commonly observed pattern in laboratory tests wit smoot boundary (Tatsuoka et al. 1990). Wen te fabric is assumed constant (te fabric is proibited from evolving wit deformation), te

4 24 Z. Gao and J. Zao (a) Type-a sear band (b) Type-b sear band Fig. 3 Predicted strain localization for te cases wit roug boundary condition and a evolving fabric; b constant fabric at vertical strain ε = / = 15.5 % strain concentration in te sear band (Fig. 2b) will be more intense tan te fabric evolution case (Fig. 2a). In addition, it is observed tat te predicted vertical strain level corresponding to sear band initiation is independent of initial bedding plane orientation if te fabric is assumed constant, wic appears to be inconsistent wit te observations in Tatsuoka et al. (1990). Tis indicates tat fabric evolution sould be properly accounted for in te strain localization analysis. Figure 3 sows te predicted sear band pattern in roug boundary condition cases (e 0 = 0.7, σ c = 400 kpa and α = 45 ). Two symmetric and asymmetric sear bands are observed for te cases wit and witout fabric evolution, respectively. Te sear strain level in Type-a sear bands [according to te definition by Tatsuoka et al. (1990)] is iger tan tat in te Type-b sear bands in bot cases. 3 Mecanisms Governing te Sear Band Patterns Our study sows tat te development of a sear band is essentially governed by two competing mecanisms te fabric evolution wic reduces te non-coaxial strain increment wen te fabric and applied stress are initially non-coaxial and alleviates te strain localization and te structural constraint imposed by te boundary conditions wic promotes te sear band development. Te structural constraint can be better described by te reaction forces at te top and bottom sample ends as sown in Fig. 4. In te smoot boundary condition cases, te initial bedding plane constitutes a natural weakened plane along wic te sample can develop sear strain concentration and te vertical reaction forces on te two ends drive te upper and lower

5 Fabric Evolution and Its Effect on Strain 25 (a) (b) (c) (d) R Positive orizontal reaction force b a Bedding plane α R Total reaction force R wen R < 0 f Total reaction force R wen R = 0 f Total reaction force R wen R > 0 f Fig. 4 Definition of a te positive direction of te orizontal reaction force R, c d tree cases of te total reaction force imposed on te sample by te boundary Fig. 5 Evolution of te orizontal reaction force at te top end of te sand sample wit α = 45 and roug boundary condition (e 0 = 0.7, σ 3 = 400 kpa): a considering fabric evolution; b considering constant fabric parts of te sample to te rigt and left respectively. Tis causes te occurrence of a single Type-b sear band. Since te strain localization initiates at relatively low strain level, te effect of fabric evolution is not large enoug to prevent te sear band development in tis case. In te roug boundary condition cases, negative R develops to prevent te orizontal displacement of te top and bottom ends due to te non-coaxial strain increment at te initial loading stage (Fig. 5), wic leads to a Type-a sear band initially. As te sample is weakened along direction-a due to stain localization (Fig. 5a), te magnitude of R decreases and becomes positive subsequently. Tis causes te initiation of Type-b sear band. Meanwile, te evolution of fabric can eventually lead to rater symmetric geometry of te two bands as well as sample sape (Fig. 4a). If te fabric is fixed witout evolution, te final sample sape and sear bands are asymmetric as R will become negative again later on (Fig. 5b).

6 26 Z. Gao and J. Zao 4 Conclusions Te effect of fabric evolution on strain localization in plane strain compression as been studied numerically. It is found tat te evolution of fabric generally tends to make te sand response more coaxial wit te applied stress and alleviate te strain localization, wile te structural constraint induced by te sample ends leads to more inomogeneous deformation witin te sand sample wen te initial fabric is non-coaxial wit te applied stress. In smoot boundary condition cases, structural constraint dominates over te fabric evolution and leads to te formation of a single sear band. Wen te boundary condition is roug, te structural constraint may play a comparable role wit fabric evolution, wic leads to symmetric cross-sape sear bands. If te soil fabric is not allowed to evolve, asymmetric cross-sape sear bands develop. Acknowledgements Te study was financially supported by RGC/GRF and DAG08/09.EG04. References Borja RI, Song XY, Recenmacer AL, Abedi S, Wu W (2013) Sear band in sand wit spatially varying density. J Mec Pys Solids 61: Gao ZW, Zao JD (2013) Strain localization and fabric evolution in sand. Int J Solid Struct 50: Gao ZW, Zao JD, Li XS, Dafalias YF (2014) A critical state sand plasticity model accounting for fabric evolution. Int J Numer Anal Met Geomec 38: Guo N, Zao JD (2013) Te signature of sear induced anisotropy in granular media. Comput Geotec 47:1 15 Tatsuoka F, Nakamura S, Huang CC, Tani K (1990) Strengt anisotropy and sear band direction in plane strain tests of sand. Soils Found 30(1):35 54 Tejcman J, Bauer E, Wu W (2007) Effect of fabric anisotropy on sear localization in sand during plane strain compression. Acta Mec 189:23 51 Zao JD, Guo N (2013a) Unique critical state caracteristics in granular media considering fabric anisotropy. Géotecnique 63(8): Zao JD, Guo N (2013b) A new definition on critical state of granular media accounting for fabric anisotropy. In: Powders and grains 2013: AIP conference proceedings, vol 1542, pp doi: /

Path to static failure of machine components

Path to static failure of machine components Pat to static failure of macine components Load Stress Discussed last week (w) Ductile material Yield Strain Brittle material Fracture Fracture Dr. P. Buyung Kosasi,Spring 008 Name some of ductile and

More information

Constitutive modelling of fabric anisotropy in sand

Constitutive modelling of fabric anisotropy in sand Geomechanics from Micro to Macro Soga et al. (Eds) 2015 Taylor & Francis Group, London, ISBN 978-1-138-02707-7 Constitutive modelling of fabric anisotropy in sand Z.W. Gao School of Engineering, University

More information

Discrete Element Modelling of Methane Hydrate Soil Sediments using. Elongated Soil Particles

Discrete Element Modelling of Methane Hydrate Soil Sediments using. Elongated Soil Particles Discrete Element Modelling of Metane Hydrate Soil Sediments using Elongated Soil Particles Yanxin Yu a, 1 (yanxin.yu.@alumni.ucl.ac.uk) Yi Pik Ceng a (yi.ceng@ucl.ac.uk) Xiaomin Xu b (xuxiaominzju@gmail.com)

More information

CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY

CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY CALCULATION OF COLLAPSE PRESSURE IN SHALE GAS FORMATION AND THE INFLUENCE OF FORMATION ANISOTROPY L.Hu, J.Deng, F.Deng, H.Lin, C.Yan, Y.Li, H.Liu, W.Cao (Cina University of Petroleum) Sale gas formations

More information

6. Non-uniform bending

6. Non-uniform bending . Non-uniform bending Introduction Definition A non-uniform bending is te case were te cross-section is not only bent but also seared. It is known from te statics tat in suc a case, te bending moment in

More information

3. Using your answers to the two previous questions, evaluate the Mratio

3. Using your answers to the two previous questions, evaluate the Mratio MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0219 2.002 MECHANICS AND MATERIALS II HOMEWORK NO. 4 Distributed: Friday, April 2, 2004 Due: Friday,

More information

Journal of Engineering Science and Technology Review 7 (4) (2014) 40-45

Journal of Engineering Science and Technology Review 7 (4) (2014) 40-45 Jestr Journal of Engineering Science and Tecnology Review 7 (4) (14) -45 JOURNAL OF Engineering Science and Tecnology Review www.jestr.org Mecanics Evolution Caracteristics Analysis of in Fully-mecanized

More information

Effect of surface tension on swell-induced surface instability of substrate-confined hydrogel layers

Effect of surface tension on swell-induced surface instability of substrate-confined hydrogel layers PAPER CREATED USING THE RSC ARTICLE TEMPLATE (VER. 3.) - SEE WWW.RSC.ORG/ELECTRONICFILES FOR DETAILS www.rsc.org/xxxxxx XXXXXXXX 5 5 Effect of surface tension on swell-induced surface instability of substrate-confined

More information

Comment on Experimental observations of saltwater up-coning

Comment on Experimental observations of saltwater up-coning 1 Comment on Experimental observations of saltwater up-coning H. Zang 1,, D.A. Barry 2 and G.C. Hocking 3 1 Griffit Scool of Engineering, Griffit University, Gold Coast Campus, QLD 4222, Australia. Tel.:

More information

AN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES

AN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES Ice in te Environment: Proceedings of te 1t IAHR International Symposium on Ice Dunedin, New Zealand, nd t December International Association of Hydraulic Engineering and Researc AN ANALYSIS OF AMPLITUDE

More information

Distribution of reynolds shear stress in steady and unsteady flows

Distribution of reynolds shear stress in steady and unsteady flows University of Wollongong Researc Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 13 Distribution of reynolds sear stress in steady

More information

Finite Element Analysis of J-Integral for Surface Cracks in Round Bars under Combined Mode I Loading

Finite Element Analysis of J-Integral for Surface Cracks in Round Bars under Combined Mode I Loading nternational Journal of ntegrated Engineering, Vol. 9 No. 2 (207) p. -8 Finite Element Analysis of J-ntegral for Surface Cracks in Round Bars under Combined Mode Loading A.E smail, A.K Ariffin 2, S. Abdulla

More information

NCCI: Simple methods for second order effects in portal frames

NCCI: Simple methods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames Tis NCC presents information

More information

Upper And Lower Bound Solution For Dynamic Active Earth Pressure on Cantilever Walls

Upper And Lower Bound Solution For Dynamic Active Earth Pressure on Cantilever Walls Upper And Lower Bound Solution For Dynamic Active Eart Pressure on Cantilever Walls A. Scotto di Santolo, A. Evangelista Università di Napoli Federico, Italy S. Aversa Università di Napoli Partenope, Italy

More information

Evaluation and Accurate Estimation from Petrophysical Parameters of a Reservoir

Evaluation and Accurate Estimation from Petrophysical Parameters of a Reservoir American Journal of Environmental Engineering and Science 2016; 3(2): 68-74 ttp://www.aascit.org/journal/ajees ISSN: 2381-1153 (Print); ISSN: 2381-1161 (Online) Evaluation and Accurate Estimation from

More information

FROM THEORY TO PRACTICE IN GEOTECHNICAL ENGINEERING: THE MISSING LINKS. Yannis F. Dafalias

FROM THEORY TO PRACTICE IN GEOTECHNICAL ENGINEERING: THE MISSING LINKS. Yannis F. Dafalias FROM THEORY TO PRACTICE IN GEOTECHNICAL ENGINEERING: THE MISSING LINKS Yannis F. Dafalias Department of Civil and Environmental Engineering, University of California at Davis, & Department of Mechanics,

More information

Velocity distribution in non-uniform/unsteady flows and the validity of log law

Velocity distribution in non-uniform/unsteady flows and the validity of log law University of Wollongong Researc Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 3 Velocity distribution in non-uniform/unsteady

More information

A Multiaxial Variable Amplitude Fatigue Life Prediction Method Based on a Plane Per Plane Damage Assessment

A Multiaxial Variable Amplitude Fatigue Life Prediction Method Based on a Plane Per Plane Damage Assessment American Journal of Mecanical and Industrial Engineering 28; 3(4): 47-54 ttp://www.sciencepublisinggroup.com/j/ajmie doi:.648/j.ajmie.2834.2 ISSN: 2575-679 (Print); ISSN: 2575-66 (Online) A Multiaxial

More information

PECULIARITIES OF THE WAVE FIELD LOCALIZATION IN THE FUNCTIONALLY GRADED LAYER

PECULIARITIES OF THE WAVE FIELD LOCALIZATION IN THE FUNCTIONALLY GRADED LAYER Materials Pysics and Mecanics (5) 5- Received: Marc 7, 5 PECULIARITIES OF THE WAVE FIELD LOCALIZATION IN THE FUNCTIONALLY GRADED LAYER Т.I. Belyankova *, V.V. Kalincuk Soutern Scientific Center of Russian

More information

Deviation from Linear Elastic Fracture in Near-Surface Hydraulic Fracturing Experiments with Rock Makhnenko, R.Y.

Deviation from Linear Elastic Fracture in Near-Surface Hydraulic Fracturing Experiments with Rock Makhnenko, R.Y. ARMA 10-237 Deviation from Linear Elastic Fracture in Near-Surface Hydraulic Fracturing Experiments wit Rock Maknenko, R.Y. University of Minnesota, Minneapolis, MN, USA Bunger, A.P. CSIRO Eart Science

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

Friction Coefficient s Numerical Determination for Hot Flat Steel Rolling at Low Strain Rate

Friction Coefficient s Numerical Determination for Hot Flat Steel Rolling at Low Strain Rate American ournal of Applied Matematics 05; 3(5: -8 Publised online September 6, 05 (ttp://www.sciencepublisinggroup.com/j/ajam doi: 0.648/j.ajam.050305.3 ISS: 330-0043 (Print; ISS: 330-006X (Online riction

More information

Bending analysis of a functionally graded piezoelectric cantilever beam

Bending analysis of a functionally graded piezoelectric cantilever beam Science in Cina Series G: Pysics Mecanics & Astronomy 7 Science in Cina Press Springer-Verlag Bending analysis of a functionally graded pieoelectric cantilever beam YU Tao & ZHONG Zeng Scool of Aerospace

More information

Physics 121, April 1, Equilibrium. Physics 121. April 1, Physics 121. April 1, Course Information. Discussion of Exam # 2

Physics 121, April 1, Equilibrium. Physics 121. April 1, Physics 121. April 1, Course Information. Discussion of Exam # 2 Pysics 121, April 1, 2008. Pysics 121. April 1, 2008. Course Information Discussion of Exam # 2 Topics to be discussed today: Requirements for Equilibrium Gravitational Equilibrium Sample problems Pysics

More information

Nonlinear correction to the bending stiffness of a damaged composite beam

Nonlinear correction to the bending stiffness of a damaged composite beam Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. Nonlinear correction to te bending stiffness of a damaged composite beam W.

More information

INTERNAL RESISTANCE OPTIMIZATION OF A HELMHOLTZ RESONATOR IN NOISE CONTROL OF SMALL ENCLOSURES. Ganghua Yu, Deyu Li and Li Cheng 1 I.

INTERNAL RESISTANCE OPTIMIZATION OF A HELMHOLTZ RESONATOR IN NOISE CONTROL OF SMALL ENCLOSURES. Ganghua Yu, Deyu Li and Li Cheng 1 I. ICV4 Cairns Australia 9- July, 7 ITAL ITAC OPTIMIZATIO OF A HLMHOLTZ OATO I OI COTOL OF MALL CLOU Gangua Yu, Deyu Li and Li Ceng Department of Mecanical ngineering, Te Hong Kong Polytecnic University Hung

More information

Static Response Analysis of a FGM Timoshenko s Beam Subjected to Uniformly Distributed Loading Condition

Static Response Analysis of a FGM Timoshenko s Beam Subjected to Uniformly Distributed Loading Condition Static Response Analysis of a FGM Timoseno s Beam Subjected to Uniformly Distributed Loading Condition 8 Aas Roy Department of Mecanical Engineering National Institute of Tecnology Durgapur Durgapur, Maatma

More information

Chapters 19 & 20 Heat and the First Law of Thermodynamics

Chapters 19 & 20 Heat and the First Law of Thermodynamics Capters 19 & 20 Heat and te First Law of Termodynamics Te Zerot Law of Termodynamics Te First Law of Termodynamics Termal Processes Te Second Law of Termodynamics Heat Engines and te Carnot Cycle Refrigerators,

More information

Differential Settlement of Foundations on Loess

Differential Settlement of Foundations on Loess Missouri Uniersity of Science and Tecnology Scolars' Mine International Conference on Case Histories in Geotecnical Engineering (013) - Seent International Conference on Case Histories in Geotecnical Engineering

More information

Seepage Analysis through Earth Dam Based on Finite Difference Method

Seepage Analysis through Earth Dam Based on Finite Difference Method J. Basic. Appl. Sci. Res., (11)111-1, 1 1, TetRoad Publication ISSN -44 Journal of Basic and Applied Scientific Researc www.tetroad.com Seepage Analysis troug Eart Dam Based on Finite Difference Metod

More information

FULL SHEAR DEFORMATION FOR ANALYSIS OF THICK PLATE

FULL SHEAR DEFORMATION FOR ANALYSIS OF THICK PLATE FULL SHAR DFORMATION FOR ANALYSIS OF THICK PLAT 1 IBARUGBULM, OWUS M., 2 NJOKU, KLCHI O., 3 ZIFULA, UCHCHI G. 1,2,3 Civil ngineering Department, Federal University of Tecnology, Owerri, Nigeria -mail:

More information

Non-linear Analysis Method of Ground Response Using Equivalent Single-degree-of-freedom Model

Non-linear Analysis Method of Ground Response Using Equivalent Single-degree-of-freedom Model Proceedings of te Tent Pacific Conference on Eartquake Engineering Building an Eartquake-Resilient Pacific 6-8 November 25, Sydney, Australia Non-linear Analysis Metod of Ground Response Using Equivalent

More information

Empirical models for estimating liquefaction-induced lateral spread displacement

Empirical models for estimating liquefaction-induced lateral spread displacement Empirical models for estimating liquefaction-induced lateral spread displacement J.J. Zang and J.X. Zao Institute of Geological & Nuclear Sciences Ltd, Lower Hutt, New Zealand. 2004 NZSEE Conference ABSTRACT:

More information

How to Find the Derivative of a Function: Calculus 1

How to Find the Derivative of a Function: Calculus 1 Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te

More information

The entransy dissipation minimization principle under given heat duty and heat transfer area conditions

The entransy dissipation minimization principle under given heat duty and heat transfer area conditions Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area

More information

232 Calculus and Structures

232 Calculus and Structures 3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE

More information

MVT and Rolle s Theorem

MVT and Rolle s Theorem AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state

More information

Click here to see an animation of the derivative

Click here to see an animation of the derivative Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,

More information

Introduction to Derivatives

Introduction to Derivatives Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))

More information

1 Power is transferred through a machine as shown. power input P I machine. power output P O. power loss P L. What is the efficiency of the machine?

1 Power is transferred through a machine as shown. power input P I machine. power output P O. power loss P L. What is the efficiency of the machine? 1 1 Power is transferred troug a macine as sown. power input P I macine power output P O power loss P L Wat is te efficiency of te macine? P I P L P P P O + P L I O P L P O P I 2 ir in a bicycle pump is

More information

1. State whether the function is an exponential growth or exponential decay, and describe its end behaviour using limits.

1. State whether the function is an exponential growth or exponential decay, and describe its end behaviour using limits. Questions 1. State weter te function is an exponential growt or exponential decay, and describe its end beaviour using its. (a) f(x) = 3 2x (b) f(x) = 0.5 x (c) f(x) = e (d) f(x) = ( ) x 1 4 2. Matc te

More information

Analysis of Static and Dynamic Load on Hydrostatic Bearing with Variable Viscosity and Pressure

Analysis of Static and Dynamic Load on Hydrostatic Bearing with Variable Viscosity and Pressure Indian Journal of Science and Tecnology Supplementary Article Analysis of Static and Dynamic Load on Hydrostatic Bearing wit Variable Viscosity and Pressure V. Srinivasan* Professor, Scool of Mecanical

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

lecture 26: Richardson extrapolation

lecture 26: Richardson extrapolation 43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)

More information

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx. Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions

More information

Available online at ScienceDirect. Procedia Engineering 125 (2015 )

Available online at   ScienceDirect. Procedia Engineering 125 (2015 ) Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 15 (015 ) 1065 1069 Te 5t International Conference of Euro Asia Civil Engineering Forum (EACEF-5) Identification of aerodynamic

More information

CHAPTER 2 MODELING OF THREE-TANK SYSTEM

CHAPTER 2 MODELING OF THREE-TANK SYSTEM 7 CHAPTER MODELING OF THREE-TANK SYSTEM. INTRODUCTION Te interacting tree-tank system is a typical example of a nonlinear MIMO system. Heiming and Lunze (999) ave regarded treetank system as a bencmark

More information

Shrinkage anisotropy characteristics from soil structure and initial sample/layer size. V.Y. Chertkov*

Shrinkage anisotropy characteristics from soil structure and initial sample/layer size. V.Y. Chertkov* Srinkage anisotropy caracteristics from soil structure and initial sample/layer size V.Y. Certkov* Division of Environmental, Water, and Agricultural Engineering, Faculty of Civil and Environmental Engineering,

More information

Taylor Series and the Mean Value Theorem of Derivatives

Taylor Series and the Mean Value Theorem of Derivatives 1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential

More information

Reynolds number effects in stratified turbulent wakes

Reynolds number effects in stratified turbulent wakes Abstract Reynolds number effects in stratified turbulent wakes Qi Zou 1 and Peter Diamessis 2 1 Department of Applied Matematics and Teoretical Pysics, University of Cambridge 2 Scool of Civil and Environmental

More information

A Modified Distributed Lagrange Multiplier/Fictitious Domain Method for Particulate Flows with Collisions

A Modified Distributed Lagrange Multiplier/Fictitious Domain Method for Particulate Flows with Collisions A Modified Distributed Lagrange Multiplier/Fictitious Domain Metod for Particulate Flows wit Collisions P. Sing Department of Mecanical Engineering New Jersey Institute of Tecnology University Heigts Newark,

More information

FINITE ELEMENT STOCHASTIC ANALYSIS

FINITE ELEMENT STOCHASTIC ANALYSIS FINITE ELEMENT STOCHASTIC ANALYSIS Murray Fredlund, P.D., P.Eng., SoilVision Systems Ltd., Saskatoon, SK ABSTRACT Numerical models can be valuable tools in te prediction of seepage. Te results can often

More information

f a h f a h h lim lim

f a h f a h h lim lim Te Derivative Te derivative of a function f at a (denoted f a) is f a if tis it exists. An alternative way of defining f a is f a x a fa fa fx fa x a Note tat te tangent line to te grap of f at te point

More information

Exercise 19 - OLD EXAM, FDTD

Exercise 19 - OLD EXAM, FDTD Exercise 19 - OLD EXAM, FDTD A 1D wave propagation may be considered by te coupled differential equations u x + a v t v x + b u t a) 2 points: Derive te decoupled differential equation and give c in terms

More information

An optimum design of robotic food handling by using Burger model

An optimum design of robotic food handling by using Burger model DOI 10.1007/s11370-008-0032-5 ORIGINAL RESEARCH PAPER An optimum design of robotic food andling by using Burger model Naoki Sakamoto Mitsuru Higasimori Tosio Tsuji Makoto Kaneko Received: 28 February 2008

More information

Sample Problems for Exam II

Sample Problems for Exam II Sample Problems for Exam 1. Te saft below as lengt L, Torsional stiffness GJ and torque T is applied at point C, wic is at a distance of 0.6L from te left (point ). Use Castigliano teorem to Calculate

More information

UNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics

UNIVERSITY OF SASKATCHEWAN Department of Physics and Engineering Physics UNIVERSITY O SASKATCHEWAN Department of Pysics and Engineering Pysics Pysics 117.3 MIDTERM EXAM Regular Sitting NAME: (Last) Please Print (Given) Time: 90 minutes STUDENT NO.: LECTURE SECTION (please ceck):

More information

COMPARISON OF FUZZY LOGIC CONTROLLERS FOR A MULTIVARIABLE PROCESS

COMPARISON OF FUZZY LOGIC CONTROLLERS FOR A MULTIVARIABLE PROCESS COMPARISON OF FUZZY LOGIC CONTROLLERS FOR A MULTIVARIABLE PROCESS KARTHICK S, LAKSHMI P, DEEPA T 3 PG Student, DEEE, College of Engineering, Guindy, Anna University, Cennai Associate Professor, DEEE, College

More information

The causes of the observed head-in-pillow (HnP) soldering defects in BGA

The causes of the observed head-in-pillow (HnP) soldering defects in BGA ORIGINAL ARTICLE Solder joints in surface mounted IC assemblies: Relief in stress and warpage owing to te application of elevated stand-off eigts E. Suir E. Suir. Solder joints in surface mounted IC assemblies:

More information

Effects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates with Ramped Wall Temperature

Effects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates with Ramped Wall Temperature Volume 39 No. February 01 Effects of Radiation on Unsteady Couette Flow between Two Vertical Parallel Plates wit Ramped Wall Temperature S. Das Department of Matematics University of Gour Banga Malda 73

More information

Optimization of the thin-walled rod with an open profile

Optimization of the thin-walled rod with an open profile (1) DOI: 1.151/ matecconf/181 IPICSE-1 Optimization of te tin-walled rod wit an open profile Vladimir Andreev 1,* Elena Barmenkova 1, 1 Moccow State University of Civil Engineering, Yaroslavskoye s., Moscow

More information

A general articulation angle stability model for non-slewing articulated mobile cranes on slopes *

A general articulation angle stability model for non-slewing articulated mobile cranes on slopes * tecnical note 3 general articulation angle stability model for non-slewing articulated mobile cranes on slopes * J Wu, L uzzomi and M Hodkiewicz Scool of Mecanical and Cemical Engineering, University of

More information

MAT Calculus for Engineers I EXAM #1

MAT Calculus for Engineers I EXAM #1 MAT 65 - Calculus for Engineers I EXAM # Instructor: Liu, Hao Honor Statement By signing below you conrm tat you ave neiter given nor received any unautorized assistance on tis eam. Tis includes any use

More information

1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).

1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t). . Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use

More information

The development of contact and noncontact technique to study the heat dissipation in metals under loading

The development of contact and noncontact technique to study the heat dissipation in metals under loading Te development of contact and noncontact tecnique to study te eat dissipation in metals under loading More info about tis article: ttp://www.ndt.net/?id=73 Abstract * ICMM UB RAS, Ac. Koroleva Str., 63

More information

Study of Convective Heat Transfer through Micro Channels with Different Configurations

Study of Convective Heat Transfer through Micro Channels with Different Configurations International Journal of Current Engineering and Tecnology E-ISSN 2277 4106, P-ISSN 2347 5161 2016 INPRESSCO, All Rigts Reserved Available at ttp://inpressco.com/category/ijcet Researc Article Study of

More information

Research Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic Spline Quasi Interpolation

Research Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic Spline Quasi Interpolation Advances in Numerical Analysis Volume 204, Article ID 35394, 8 pages ttp://dx.doi.org/0.55/204/35394 Researc Article Error Analysis for a Noisy Lacunary Cubic Spline Interpolation and a Simple Noisy Cubic

More information

Large eddy simulation of turbulent flow downstream of a backward-facing step

Large eddy simulation of turbulent flow downstream of a backward-facing step Available online at www.sciencedirect.com Procedia Engineering 31 (01) 16 International Conference on Advances in Computational Modeling and Simulation Large eddy simulation of turbulent flow downstream

More information

Mass Lumping for Constant Density Acoustics

Mass Lumping for Constant Density Acoustics Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial

More information

Research Article Solution to Two-Dimensional Steady Inverse Heat Transfer Problems with Interior Heat Source Based on the Conjugate Gradient Method

Research Article Solution to Two-Dimensional Steady Inverse Heat Transfer Problems with Interior Heat Source Based on the Conjugate Gradient Method Hindawi Matematical Problems in Engineering Volume 27, Article ID 286342, 9 pages ttps://doiorg/55/27/286342 Researc Article Solution to Two-Dimensional Steady Inverse Heat Transfer Problems wit Interior

More information

Microstrip Antennas- Rectangular Patch

Microstrip Antennas- Rectangular Patch April 4, 7 rect_patc_tl.doc Page of 6 Microstrip Antennas- Rectangular Patc (Capter 4 in Antenna Teory, Analysis and Design (nd Edition) by Balanis) Sown in Figures 4. - 4.3 Easy to analyze using transmission

More information

Optimal parameters for a hierarchical grid data structure for contact detection in arbitrarily polydisperse particle systems

Optimal parameters for a hierarchical grid data structure for contact detection in arbitrarily polydisperse particle systems Comp. Part. Mec. 04) :357 37 DOI 0.007/s4057-04-000-9 Optimal parameters for a ierarcical grid data structure for contact detection in arbitrarily polydisperse particle systems Dinant Krijgsman Vitaliy

More information

CONTINUUM TOPOLOGY OPTIMIZATION OF BUCKLING-SENSITIVE STRUCTURES

CONTINUUM TOPOLOGY OPTIMIZATION OF BUCKLING-SENSITIVE STRUCTURES I-2002-575 CONINUUM OPOLOGY OPIMIZION OF BUCLING-SENSIIVE SRUCURES Salam Ramatalla, Graduate Researc ssistant Colby C. Swan, ssociate Professor, member I Civil and Environmental Engineering, Center for

More information

Quasiperiodic phenomena in the Van der Pol - Mathieu equation

Quasiperiodic phenomena in the Van der Pol - Mathieu equation Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van

More information

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.

Lecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator. Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative

More information

HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS

HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3

More information

The Dynamic Range of Bursting in a Model Respiratory Pacemaker Network

The Dynamic Range of Bursting in a Model Respiratory Pacemaker Network SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 4, No. 4, pp. 117 1139 c 25 Society for Industrial and Applied Matematics Te Dynamic Range of Bursting in a Model Respiratory Pacemaker Network Janet Best, Alla Borisyuk,

More information

The structure of the atoms

The structure of the atoms Te structure of te atoms Atomos = indivisible University of Pécs, Medical Scool, Dept. Biopysics All tat exists are atoms and empty space; everyting else is merely tougt to exist. Democritus, 415 B.C.

More information

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a?

(a) At what number x = a does f have a removable discontinuity? What value f(a) should be assigned to f at x = a in order to make f continuous at a? Solutions to Test 1 Fall 016 1pt 1. Te grap of a function f(x) is sown at rigt below. Part I. State te value of eac limit. If a limit is infinite, state weter it is or. If a limit does not exist (but is

More information

Stability of Smart Beams with Varying Properties Based on the First Order Shear Deformation Theory Located on a Continuous Elastic Foundation

Stability of Smart Beams with Varying Properties Based on the First Order Shear Deformation Theory Located on a Continuous Elastic Foundation Australian Journal of Basic and Applied Sciences, 5(7): 743-747, ISSN 99-878 Stability of Smart Beams wit Varying Properties Based on te First Order Sear Deformation Teory ocated on a Continuous Elastic

More information

The total error in numerical differentiation

The total error in numerical differentiation AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and

More information

Stress analysis of laminated glass with different interlayer materials

Stress analysis of laminated glass with different interlayer materials Alexandria Engineering Journal (01) 51, 61 67 Alexandria University Alexandria Engineering Journal www.elsevier.com/locate/aej www.sciencedirect.com Stress analysis of laminated glass wit different interlayer

More information

REVIEW LAB ANSWER KEY

REVIEW LAB ANSWER KEY REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g

More information

Optimal Shape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow

Optimal Shape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow THE 5 TH ASIAN COMPUTAITIONAL FLUID DYNAMICS BUSAN, KOREA, OCTOBER 7 ~ OCTOBER 30, 003 Optimal Sape Design of a Two-dimensional Asymmetric Diffsuer in Turbulent Flow Seokyun Lim and Haeceon Coi. Center

More information

DEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT: BEHAVIOUR AND MODELLING

DEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT: BEHAVIOUR AND MODELLING Asia-Pacific Conference on FRP in Structures (APFIS 2007) S.T. Smit (ed) 2007 International Institute for FRP in Construction DEBONDING FAILURES OF RC BEAMS STRENGTHENED WITH EXTERNALLY BONDED FRP REINFORCEMENT:

More information

These errors are made from replacing an infinite process by finite one.

These errors are made from replacing an infinite process by finite one. Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat

More information

Journal of Applied Science and Agriculture. The Effects Of Corrugated Geometry On Flow And Heat Transfer In Corrugated Channel Using Nanofluid

Journal of Applied Science and Agriculture. The Effects Of Corrugated Geometry On Flow And Heat Transfer In Corrugated Channel Using Nanofluid Journal o Applied Science and Agriculture, 9() February 04, Pages: 408-47 AENSI Journals Journal o Applied Science and Agriculture ISSN 86-9 Journal ome page: www.aensiweb.com/jasa/index.tml Te Eects O

More information

Discriminate Modelling of Peak and Off-Peak Motorway Capacity

Discriminate Modelling of Peak and Off-Peak Motorway Capacity International Journal of Integrated Engineering - Special Issue on ICONCEES Vol. 4 No. 3 (2012) p. 53-58 Discriminate Modelling of Peak and Off-Peak Motorway Capacity Hasim Moammed Alassan 1,*, Sundara

More information

Polynomial Interpolation

Polynomial Interpolation Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc

More information

Finding and Using Derivative The shortcuts

Finding and Using Derivative The shortcuts Calculus 1 Lia Vas Finding and Using Derivative Te sortcuts We ave seen tat te formula f f(x+) f(x) (x) = lim 0 is manageable for relatively simple functions like a linear or quadratic. For more complex

More information

EFFECTS OF LINE AND PASSIVATION GEOMETRY ON CURVATURE EVOLUTION DURING PROCESSING AND THERMAL CYCLING IN COPPER INTERCONNECT LINES

EFFECTS OF LINE AND PASSIVATION GEOMETRY ON CURVATURE EVOLUTION DURING PROCESSING AND THERMAL CYCLING IN COPPER INTERCONNECT LINES Acta mater. 48 (000) 3169±3175 www.elsevier.com/locate/actamat EFFECTS OF LINE AND PASSIVATION GEOMETRY ON CURVATURE EVOLUTION DURING PROCESSING AND THERMAL CYCLING IN COPPER INTERCONNECT LINES T.-S. PARK

More information

APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction

APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS. Introduction Acta Mat. Univ. Comenianae Vol. LXVII, 1(1998), pp. 57 68 57 APPROXIMATION OF CRYSTALLINE DENDRITE GROWTH IN TWO SPACE DIMENSIONS A. SCHMIDT Abstract. Te pase transition between solid and liquid in an

More information

Chapter 2 Ising Model for Ferromagnetism

Chapter 2 Ising Model for Ferromagnetism Capter Ising Model for Ferromagnetism Abstract Tis capter presents te Ising model for ferromagnetism, wic is a standard simple model of a pase transition. Using te approximation of mean-field teory, te

More information

Mechanical Properties of Cement Mortar: Development of Structure-Property Relationships

Mechanical Properties of Cement Mortar: Development of Structure-Property Relationships International Journal of Concrete Structures and Materials Vol.5, No.1, pp.3~10, June 011 DOI 10.4334/IJCSM.011.5.1.003 Mecanical Properties of Cement Mortar: Development of Structure-Property Relationsips

More information

Chapter 5 FINITE DIFFERENCE METHOD (FDM)

Chapter 5 FINITE DIFFERENCE METHOD (FDM) MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential

More information

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator

Simulation and verification of a plate heat exchanger with a built-in tap water accumulator Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation

More information

The Verlet Algorithm for Molecular Dynamics Simulations

The Verlet Algorithm for Molecular Dynamics Simulations Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical

More information

Reflection of electromagnetic waves from magnetic having the ferromagnetic spiral

Reflection of electromagnetic waves from magnetic having the ferromagnetic spiral Reflection of electromagnetic waves from magnetic aving te ferromagnetic spiral Igor V. Bycov 1a Dmitry A. Kuzmin 1b and Vladimir G. Savrov 3 1 Celyabins State University 51 Celyabins Br. Kasiriny Street

More information

Continuous formulation for bottom friction in free surface flows modelling

Continuous formulation for bottom friction in free surface flows modelling River Basin Management V 81 Continuous formulation for bottom friction in free surface flows modelling O. Maciels 1, S. Erpicum 1, B. J. Dewals 1, 2, P. Arcambeau 1 & M. Pirotton 1 1 HACH Unit, Department

More information

Continuity and Differentiability

Continuity and Differentiability Continuity and Dierentiability Tis capter requires a good understanding o its. Te concepts o continuity and dierentiability are more or less obvious etensions o te concept o its. Section - INTRODUCTION

More information