Development of P-Y Criterion for Anisotropic Rock And

Total Page:16

File Type:pdf, Size:1020Kb

Development of P-Y Criterion for Anisotropic Rock And DEVELOPMENT OF P-Y CRITERION FOR ANISOTROPIC ROCK AND COHESIVE INTERMEDIATE GEOMATERIALS A Dissertation Presented to The Graduate Faculty of the University of Akron In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy Ehab Salem Shatnawi August, 2008 DEVELOPMENT OF P-Y CRITERION FOR ANISOTROPIC ROCK AND COHESIVE INTERMEDIATE GEOMATERIALS Ehab Salem Shatnawi Dissertation Approved: Accepted: ______________________________ ______________________________ Advisor Department Chair Dr. Robert Liang Dr. Wieslaw Binienda ______________________________ ______________________________ Committee Member Dean of the College Dr. Craig Menzemer Dr. George K. Haritos ______________________________ ______________________________ Committee Member Dean of the Graduate School Dr. Daren Zywicki Dr. George R. Newkome ______________________________ ______________________________ Committee Member Date Dr. Xiaosheng Gao ______________________________ Committee Member Dr. Kevin Kreider ii ABSTRACT Rock-socketed drilled shaft foundations are commonly used to resist large axial and lateral loads applied to structures or as a means to stabilize an unstable slope with either marginal factor of safety or experiencing continuing slope movements. One of the widely used approaches for analyzing the response of drilled shafts under lateral loads is the p-y approach. Although there are past and ongoing research efforts to develop pertinent p-y criterion for the laterally loaded rock-socketed drilled shafts, most of these p-y curves were derived from basic assumptions that the rock mass behaves as an isotropic continuum. The assumption of isotropy may not be applicable to the rock mass with intrinsic anisotropy or the rock formation with distinguishing joints and bedding planes. Therefore, there is a need to develop a p-y curve criterion that can take into account the effects of rock anisotropy on the p-y curve of laterally loaded drilled shafts. A hyperbolic non-linear p-y criterion for rock mass that exhibit distinguished transverse isotropy is developed in this study based on both theoretical derivations and numerical (finite element) parametric analysis results. Evaluations based on parametric study on full-scale lateral load test on fully instrumented drilled shaft have shown the insights on the influences of rock anisotropy on the predicted response of the rock socketed drilled shaft under the lateral load. Both, the orientation of the plane of transversely isotropy, and the degree of anisotropy (E/E’) has influences on the main two iii parameters required to characterize the p-y curve, the subgrade modulus (Ki) and the ultimate lateral resistance (pu). In addition to the development of a hyperbolic p-y criterion of transversely isotropic rock, another p-y criterion for cohesive intermediate geomaterials (IGM) using hyperbolic mathematical formulation is developed herein by employing the results of a series of finite element (FE) simulations and the results of two full scale lateral load test for drilled shaft socketed into IGM. iv DEDICATION To my mother, my father, my wife, my sisters and brothers, my friends, and to anyone who would read this dissertation. Also, I would like to dedicate this work to my expected baby. v ACKNOWLEDGEMENTS All praise is due to Allah (S.W.T) without whom nothing good of this work can be ever accomplished. Though only my name appears on the cover of this dissertation, a great many people have contributed to its production. I owe my gratitude to all those people who have made this dissertation possible and because of whom my graduate experience has been one that I will cherish forever. My deepest gratitude is to my advisor, Dr. Robert Liang. I have been amazingly fortunate to have an advisor who gave me the freedom to explore on my own and at the same time the guidance to recover when my steps faltered. I hope that one day I would become as good an advisor to my students as Prof. Robert Liang has been to me. Also of a great importance are the help and constructive comments and contributions I received from my committee members, Dr. Craig Menzemer, Dr. Daren Zywicki, Dr. Yueh- Jaw Lin, Dr. Xiaosheng Gao and Dr. Kevin Kreider. I would also like to acknowledge the support I received from personnel in our department, particularly Mrs. Kimberly Stone and Mrs. Christina. I must acknowledge as well the many friends, and colleagues, who assisted, advised, and supported me during these three and a vi half year. Especially, I need to express my gratitude and deep appreciation to my former roommate, Dr. Mohammad Yamin, whose friendship, knowledge, and wisdom have supported, enlightened, and entertained me over the ten years of our friendship. I must also acknowledge my former advisor, Dr. Abdulla Malkawi, I have learned so much from his keen insight, his research and problem solving abilities, and his amazing energy. Special thanks also to Dr. Jamal Nusairat, Dr. Diya Azzam, Dr. Inmar Badwan, Dr. Firas Hasan, Dr. Sami Khorbatli, Dr. Samer Rababa’a, Dr. Khalid Alakhras, Dr. Wael Khasawneh, Dr. Mohammad Khasawneh, Dr. Qais Khasawneh, Saleh Khasawneh, Abdallah Sharo, Wassel Bdour, Madhar Tamneh, Mohannad Aljarrah, Khalid Elhindi, Jamal Tahat, Khalid Mustafa, and my friends in Jordan, Eng. Kifah Ewisat, Eng. Mamoun Shatnawi, Eng. Mohammad Al-sakran, and Eng. Sameer Mousa. They have consistently helped me keep perspective on what is important in life and shown me how to deal with reality. Most importantly, none of this would have been possible without the love and patience of my family. My parents (Salem and Salmeh), my darling and lovely wife SAHAR, my sisters (Heba, Ruba, Waed, and the little and lovely one Lina) and brothers (Dr. Mohammad, Abdulla, Hussain, Dr. Murad, and Moad) to whom this dissertation is dedicated to, has been a constant source of love, concern, support and strength all these years. I would like to express my heart-felt gratitude to them. vii TABLE OF CONTENTS Page LIST OF TABLES.......................................................................................................... xiii LIST OF FIGURES .........................................................................................................xiv CHAPTER I. INTRODUCTION.........................................................................................................1 1.1 Statement of Problem ............................................................................................1 1.2 Objectives ..............................................................................................................8 1.3 Work Plan ..............................................................................................................9 1.4 Dissertation Outlines ...........................................................................................13 II. LITERATURE REVIEW...........................................................................................15 2.1 Analysis Methods of Laterally Loaded Rock-Socketed Drilled Shafts...............15 2.1.1 Elastic continuum methods........................................................................16 2.1.2 Winkler Method (Subgrade reaction approach) and P-Y method .............17 2.2 Analysis Methods for Estimating Ultimate Lateral Rock Reaction ....................22 2.2.1 Carter and Kulhawy (1992).......................................................................24 2.2.2 Zhang et al. (2000) ....................................................................................25 2.2.3 To, Ernst, and Einstein (2003)...................................................................26 2.2.4 Yang (2006)...............................................................................................27 2.3 Initial Modulus of Subgrade Reaction.................................................................30 2.4 Ultimate Side Shear Resistance...........................................................................32 viii 2.5 Discussion of Analytical Models for Laterally Loaded Sockets .........................35 2.6 Bedrocks ..............................................................................................................38 2.6.1 Rocks in Ohio............................................................................................38 2.6.2 Rock Characterization ...............................................................................39 2.6.3 Rock Categories.........................................................................................40 2.7 Laboratory and In-situ testing for Transversely Isotropic Rock..........................44 III. TRANSVERSE ISOTROPY EFFECTS ON THE INITIAL TANGENT TO P-Y CURVE.............................................................................................................................49 3.1 Abstract................................................................................................................49 3.2 Introduction..........................................................................................................50 3.3 Sensitivity Analysis .............................................................................................51 3.4 FE Analysis..........................................................................................................54 3.4.1 Description of FE Model...........................................................................54 3.4.2 Constitutive Models...................................................................................55 3.4.3 FE Analysis
Recommended publications
  • Transversely Isotropic Plasticity with Application to Fiber-Reinforced Plastics
    6. LS-DYNA Anwenderforum, Frankenthal 2007 Material II Transversely Isotropic Plasticity with Application to Fiber-Reinforced Plastics M. Vogler1, S. Kolling2, R. Rolfes1 1 Leibniz University Hannover, Institute for Structural Analysis, 30167 Hannover 2 DaimlerChrysler AG, EP/SPB, HPC X271, 71059 Sindelfingen Abstract: In this article a constitutive formulation for transversely isotropic materials is presented taking large plastic deformation at small elastic strains into account. A scalar damage model is used for the approximation of the unloading behavior. Furthermore, a failure surface is assumed taking the influence of triaxiality on fracture into account. Regular- ization is considered by the introduction of an internal length for the computation of the fracture energy. The present formulation is applied to the simulation of tensile tests for a molded glass fibre reinforced polyurethane material. Keywords: Anisotropy, Transversely Isotropy, Structural Tensors, Plasticity, Anisotropic Damage, Failure, Fiber Reinforced Plastics © 2007 Copyright by DYNAmore GmbH D - II - 55 Material II 6. LS-DYNA Anwenderforum, Frankenthal 2007 1 Introduction Most of the materials that are used in the automotive industry are anisotropic to some degree. This material behavior can be observed in metals as well as in non-metallic ma- terials like molded components of fibre reinforced plastics among others. In metals, the anisotropy is induced during the manufacturing process, e.g. during sheet metal forming and due to direction of rolling. In fibre-reinforced plastics, the anisotropy is determined by the direction of the fibres. The state-of-the-art in the numerical simulation of structural parts made from fibre re- inforced plastics represents Glaser’s model [1].
    [Show full text]
  • Constitutive Relations: Transverse Isotropy and Isotropy
    Objectives_template Module 3: 3D Constitutive Equations Lecture 11: Constitutive Relations: Transverse Isotropy and Isotropy The Lecture Contains: Transverse Isotropy Isotropic Bodies Homework References file:///D|/Web%20Course%20(Ganesh%20Rana)/Dr.%20Mohite/CompositeMaterials/lecture11/11_1.htm[8/18/2014 12:10:09 PM] Objectives_template Module 3: 3D Constitutive Equations Lecture 11: Constitutive relations: Transverse isotropy and isotropy Transverse Isotropy: Introduction: In this lecture, we are going to see some more simplifications of constitutive equation and develop the relation for isotropic materials. First we will see the development of transverse isotropy and then we will reduce from it to isotropy. First Approach: Invariance Approach This is obtained from an orthotropic material. Here, we develop the constitutive relation for a material with transverse isotropy in x2-x3 plane (this is used in lamina/laminae/laminate modeling). This is obtained with the following form of the change of axes. (3.30) Now, we have Figure 3.6: State of stress (a) in x1, x2, x3 system (b) with x1-x2 and x1-x3 planes of symmetry From this, the strains in transformed coordinate system are given as: file:///D|/Web%20Course%20(Ganesh%20Rana)/Dr.%20Mohite/CompositeMaterials/lecture11/11_2.htm[8/18/2014 12:10:09 PM] Objectives_template (3.31) Here, it is to be noted that the shear strains are the tensorial shear strain terms. For any angle α, (3.32) and therefore, W must reduce to the form (3.33) Then, for W to be invariant we must have Now, let us write the left hand side of the above equation using the matrix as given in Equation (3.26) and engineering shear strains.
    [Show full text]
  • Eulerian Formulation and Level Set Models for Incompressible Fluid-Structure Interaction
    ESAIM: M2AN 42 (2008) 471–492 ESAIM: Mathematical Modelling and Numerical Analysis DOI: 10.1051/m2an:2008013 www.esaim-m2an.org EULERIAN FORMULATION AND LEVEL SET MODELS FOR INCOMPRESSIBLE FLUID-STRUCTURE INTERACTION Georges-Henri Cottet1, Emmanuel Maitre1 and Thomas Milcent1 Abstract. This paper is devoted to Eulerian models for incompressible fluid-structure systems. These models are primarily derived for computational purposes as they allow to simulate in a rather straight- forward way complex 3D systems. We first analyze the level set model of immersed membranes proposed in [Cottet and Maitre, Math. Models Methods Appl. Sci. 16 (2006) 415–438]. We in particular show that this model can be interpreted as a generalization of so-called Korteweg fluids. We then extend this model to more generic fluid-structure systems. In this framework, assuming anisotropy, the membrane model appears as a formal limit system when the elastic body width vanishes. We finally provide some numerical experiments which illustrate this claim. Mathematics Subject Classification. 76D05, 74B20, 74F10. Received July 25, 2007. Revised December 20, 2007. Published online April 1st, 2008. 1. Introduction Fluid-structure interaction problems remain among the most challenging problems in computational mechan- ics. The difficulties in the simulation of these problems stem form the fact that they rely on the coupling of models of different nature: Lagrangian in the solid and Eulerian in the fluid. The so-called ALE (for Arbitrary Lagrangian Eulerian) methods cope with this difficulty by adapting the fluid solver along the deformations of the solid medium in the direction normal to the interface. These methods allow to accurately account for continuity of stresses and velocities at the interface but they are difficult to implement and time-consuming in particular for 3D systems undergoing large deformations.
    [Show full text]
  • On Transversely Isotropic Media with Ellipsoidal Slowness Surfaces
    On Transversely Isotropic Elastic Media with Ellipsoidal Slowness Surfaces a; ;1 b;2 Anna L. Mazzucato ∗ , Lizabeth V. Rachele aDepartment of Mathematics, Pennsylvania State University, University Park, PA 16802, USA bInverse Problems Center, Rensselaer Polytechnic Institute, Troy, New York 12180, USA Abstract We consider inhomogeneous elastic media with ellipsoidal slowness surfaces. We describe all classes of transversely isotropic media for which the sheets associated to each wave mode are ellipsoids. These media have the property that elastic waves in each mode propagate along geodesic segments of certain Riemannian metrics. In particular, we study the intersection of the sheets of the slowness surface. In view of applications to the analysis of propagation of singularities along rays, we give pointwise conditions that guarantee that the sheet of the slowness surface corresponding to a given wave mode is disjoint from the others. We also investigate the smoothness of the associated polarization vectors as functions of position and direction. We employ coordinate and frame-independent methods, suitable to the study of the dynamic inverse boundary problem in elasticity. Key words: Elastodynamics, transverse isotropy, anisotropy, ellipsoidal slowness surface 2000 MSC: 74Bxx, 35Lxx 1 Introduction and Main Results In this paper we consider inhomogeneous, transversely isotropic elastic media for which the associated slowness surface has ellipsoidal sheets. Transverse isotropy is characterized by the property that at each point the elastic response of the medium is isotropic in the plane orthogonal to the so-called fiber direction. We consider media in which the fiber direction may vary smoothly from point to point. Examples of transversely isotropic elastic media include hexagonal crystals and biological tissue, such as muscle.
    [Show full text]
  • Invariant Formulation of Hyperelastic Transverse Isotropy Based on Polyconvex Free Energy Functions Joorg€ Schrooder€ A,*, Patrizio Neff B
    International Journal of Solids and Structures 40 (2003) 401–445 www.elsevier.com/locate/ijsolstr Invariant formulation of hyperelastic transverse isotropy based on polyconvex free energy functions Joorg€ Schrooder€ a,*, Patrizio Neff b a Institut f€ur Mechanik, Fachbereich 10, Universit€at Essen, 45117 Essen, Universit€atsstr. 15, Germany b Fachbereich Mathematik, Technische Universita€t Darmstadt, 64289 Darmstadt, Schloßgartenstr. 7, Germany Received 7November 2001; received in revised form 2 August 2002 Abstract In this paper we propose a formulation of polyconvex anisotropic hyperelasticity at finite strains. The main goal is the representation of the governing constitutive equations within the framework of the invariant theory which auto- matically fulfill the polyconvexity condition in the sense of Ball in order to guarantee the existence of minimizers. Based on the introduction of additional argument tensors, the so-called structural tensors, the free energies and the anisotropic stress response functions are represented by scalar-valued and tensor-valued isotropic tensor functions, respectively. In order to obtain various free energies to model specific problems which permit the matching of data stemming from experiments, we assume an additive structure. A variety of isotropic and anisotropic functions for transversely isotropic material behaviour are derived, where each individual term fulfills a priori the polyconvexity condition. The tensor generators for the stresses and moduli are evaluated in detail and some representative numerical examples are pre- sented. Furthermore, we propose an extension to orthotropic symmetry. Ó 2002 Elsevier Science Ltd. All rights reserved. Keywords: Non-linear elasticity; Anisotropy; Polyconvexity; Existence of minimizers 1. Introduction Anisotropic materials have a wide range of applications, e.g.
    [Show full text]
  • ON the ELASTIC MODULI and COMPLIANCES of TRANSVERSELY ISOTROPIC and ORTHOTROPIC MATERIALS Vlado A
    Journal of Mechanics of Materials and Structures ON THE ELASTIC MODULI AND COMPLIANCES OF TRANSVERSELY ISOTROPIC AND ORTHOTROPIC MATERIALS Vlado A. Lubarda and Michelle C. Chen Volume 3, Nº 1 January 2008 mathematical sciences publishers JOURNAL OF MECHANICS OF MATERIALS AND STRUCTURES Vol. 3, No. 1, 2008 ON THE ELASTIC MODULI AND COMPLIANCES OF TRANSVERSELY ISOTROPIC AND ORTHOTROPIC MATERIALS VLADO A. LUBARDA AND MICHELLE C. CHEN The relationships between the elastic moduli and compliances of transversely isotropic and orthotropic materials, which correspond to different appealing sets of linearly independent fourth-order base tensors used to cast the elastic moduli and compliances tensors, are derived by performing explicit inversions of the involved fourth-order tensors. The deduced sets of elastic constants are related to each other and to common engineering constants expressed in the Voigt notation with respect to the coordinate axes aligned along the directions orthogonal to the planes of material symmetry. The results are applied to a transversely isotropic monocrystalline zinc and an orthotropic human femural bone. 1. Introduction There has been a significant amount of research devoted to tensorial representation of elastic constants of anisotropic materials, which are the components of the fourth-order tensors in three-dimensional geometrical space, or the second-order tensors in six-dimensional stress or strain space. This was par- ticularly important for the spectral decompositions of the fourth-order stiffness and compliance tensors and the determination of the corresponding eigenvalues and eigentensors for materials with different types of elastic anisotropy. The spectral decomposition, for example, allows an additive decomposition of the elastic strain energy into a sum of six or fewer uncoupled energy modes, which are given by trace products of the corresponding pairs of stress and strain eigentensors.
    [Show full text]
  • Behaviour of Anisotropic Rocks Authors: M
    Behaviour of anisotropic rocks Authors: M. Sc. Mohamed Ismael, M. Sc. Lifu Chang & Prof. Dr.-Ing. habil. Heinz Konietzky (Geotechnical Institute, TU Bergakademie Freiberg) 1 Introduction ............................................................................................................. 2 2 Features causing rock anisotropy ........................................................................... 2 2.1 Primary structures ............................................................................................ 2 2.2 Secondary structures ........................................................................................ 4 2.3 Discontinuity frequency .................................................................................... 4 3 Observation and measurement of rock anisotropy ................................................. 5 3.1 Strength anisotropy .......................................................................................... 5 3.2 Stiffness anisotropy ........................................................................................ 11 3.3 Permeability anisotropy .................................................................................. 11 3.4 Seismic anisotropy ......................................................................................... 12 4 Classification of rock anisotropy ........................................................................... 13 4.1 Point load index .............................................................................................. 13 4.2
    [Show full text]
  • Structure and Frictional Properties of the Leg Joint of the Beetle Pachnoda Marginata (Scarabaeidae, Cetoniinae) As an Inspiration for Technical Joints
    biomimetics Article Structure and Frictional Properties of the Leg Joint of the Beetle Pachnoda marginata (Scarabaeidae, Cetoniinae) as an Inspiration for Technical Joints Steffen Vagts 1,* , Josef Schlattmann 1, Alexander Kovalev 2 and Stanislav N. Gorb 2 1 Department of System Technologies and Engineering Design Methodology, Hamburg University of Technology, Denickestr. 22, D-21079 Hamburg, Germany; [email protected] 2 Department of Functional Morphology and Biomechanics, Kiel University, Am Botanischen Garten 9, D-24118 Kiel, Germany; [email protected] (A.K.); [email protected] (S.N.G.) * Correspondence: steff[email protected]; Tel.:+49-40-428-784-422 Received: 21 February 2020; Accepted: 15 April 2020; Published: 20 April 2020 Abstract: The efficient locomotion of insects is not only inspiring for control algorithms but also promises innovations for the reduction of friction in joints. After previous analysis of the leg kinematics and the topological characterization of the contacting joint surfaces in the beetle Pachnoda marginata, in the present paper, we report on the measurement of the coefficient of friction within the leg joints exhibiting an anisotropic frictional behavior in different sliding directions. In addition, the simulation of the mechanical behavior of a single microstructural element helped us to understand the interactions between the contact parts of this tribological system. These findings were partly transferred to a technical contact pair which is typical for such an application as joint connectors in the automotive field. This innovation helped to reduce the coefficient of friction under dry sliding conditions up to 17%. Keywords: locomotion; walking; leg; joints; insects; Arthropoda; friction; coefficient of friction; biotribology; biomimetics 1.
    [Show full text]
  • Universal Relations for Transversely Isotropic Elastic Materials
    UniversalRelations for 'ftansverselyIsotropic Elastic Materials R. C. BATRA Department of Engineering Science and Mechanics, M/C 0219, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061, USA (Received 8 March 2002; Final version 13 March 2002) Dedicated to Millard E Beatty, respected colleague and friend. Abstract: We adoptHayes and Knops'sapproach and deriveuniversal relations for finite deformationsof a transverselyisotropic elastic material. Explicit universalrelations are obtained for homogeneousdeformations correspondingto triaxial stretches,simple shear,and simultaneousshear and extension.Universal relations are also derivedfor five families of nonhomogeneousdeformations. 1. INTRODUcnON Truesdell and Noll [1, section 54] and Beatty [2] have lucidly stated the importance of universal relations for finite deformations of isotropic elastic materials; the same remarks apply to transverselyisotropic elastic materials. Beatty [2, 3] has derived a class of universal relations for constrainedand unconstrainedisotropic elastic materials. Pucci and Saccomandi [4], have given a general approach for finding universal relations in continuum mechanics. Saccomandi and Vianello [5] have given a variational characterization of the universal relations for hemitropic materials. They have proved that the classof transverselyhemitropic hyperelasticbodies is characterizedby the condition that the inner product betweenSC - CS and W N vanishes for all deformations of the body; S is the second Piola-Kirchhoff stress tensor, C the right
    [Show full text]
  • PDF (Chapter 15. Anisotropy)
    Theory of the Earth Don L. Anderson Chapter 15. Anisotropy Boston: Blackwell Scientific Publications, c1989 Copyright transferred to the author September 2, 1998. You are granted permission for individual, educational, research and noncommercial reproduction, distribution, display and performance of this work in any format. Recommended citation: Anderson, Don L. Theory of the Earth. Boston: Blackwell Scientific Publications, 1989. http://resolver.caltech.edu/CaltechBOOK:1989.001 A scanned image of the entire book may be found at the following persistent URL: http://resolver.caltech.edu/CaltechBook:1989.001 Abstract: Anisotropy is responsible for the largest variations in seismic velocities; changes in the preferred orientation of mantle minerals, or in the direction of seismic waves, cause larger changes in velocity than can be accounted for by changes in temperature, composition or mineralogy. Therefore, discussions of velocity gradients, both radial and lateral, and chemistry and mineralogy of the mantle must allow for the presence of anisotropy. Anisotropy is not a second-order effect. Seismic data that are interpreted in terms of isotropic theory can lead to models that are not even approximately correct. On the other hand a wealth of new information regarding mantle mineralogy and flow will be available as the anisotropy of the mantle becomes better understood. Anisotropy And perpendicular now and now transverse, Pierce the dark soil and as they pierce and pass Make bare the secrets of the Earth's deep heart. nisotropy is responsible for the largest variations in tropic solid with these same elastic constants can be mod- A seismic velocities; changes in the preferred orienta- eled exactly as a stack of isotropic layers.
    [Show full text]
  • The Elastic and Electric Fields for Three-Dimensional Contact For
    International Journal of Solids and Structures 37 (2000) 3201±3229 www.elsevier.com/locate/ijsolstr The elastic and electric ®elds for three-dimensional contact for transversely isotropic piezoelectric materials Hao-Jiang Ding*, Peng-Fei Hou, Feng-Lin Guo Department of Civil Engineering, Zhejiang University, Hangzhou, 310027, People's Republic of China Received 30 May 1998; in revised form 14 January 1999 Abstract Firstly, the extended Boussinesq and Cerruti solutions for point forces and point charge acting on the surface of a transversely isotropic piezoelectric half-space are derived. Secondly, aiming at a series of common three-dimensional contact including spherical contact, a conical indentor and an upright circular ¯at punch on a transversely isotropic piezoelectric half-space, we solve for their elastic and electric ®elds in smooth and frictional cases by ®rst evaluating the displacement functions and then dierentiating. The displacement functions can be obtained by integrating the extended Boussinesq or Cerruti solutions in the contact region. Then, when only normal pressure is loaded, the stresses in the half-space of PZT-4 piezoelectric ceramic are compared in the ®gures with those of the transversely isotropic material which are assumed to have the same elastic constants as those of PZT-4. Meanwhile, the electric components in the half-space of PZT-4 piezoelectric ceramic are also shown in the same ®gures. # 2000 Elsevier Science Ltd. All rights reserved. 1. Introduction Since Hertz (1882) published his classic article `On the contact of elastic solids', the research on the contact of elastic materials has been conducted for more than one hundred years, and a lot of scientists including mathematicians, physicists and civil engineers all contributed to this area.
    [Show full text]
  • Contact Mechanics
    International Journal of Solids and Structures 37 (2000) 29±43 www.elsevier.com/locate/ijsolstr Contact mechanics J.R. Barber a,*, M. Ciavarella b, 1 aDepartment of Mechanical Engineering and Applied Mechanics, University of Michigan, Ann Arbor, MI 48109-2125, USA bDepartment of Engineering Science, University of Oxford, Parks Road, Oxford, OX1 3PJ, UK Abstract Contact problems are central to Solid Mechanics, because contact is the principal method of applying loads to a deformable body and the resulting stress concentration is often the most critical point in the body. Contact is characterized by unilateral inequalities, describing the physical impossibility of tensile contact tractions (except under special circumstances) and of material interpenetration. Additional inequalities and/or non-linearities are introduced when friction laws are taken into account. These complex boundary conditions can lead to problems with existence and uniqueness of quasi-static solution and to lack of convergence of numerical algorithms. In frictional problems, there can also be lack of stability, leading to stick±slip motion and frictional vibrations. If the material is non-linear, the solution of contact problems is greatly complicated, but recent work has shown that indentation of a power-law material by a power law punch is self-similar, even in the presence of friction, so that the complete history of loading in such cases can be described by the (usually numerical) solution of a single problem. Real contacting surfaces are rough, leading to the concentration of contact in a cluster of microscopic actual contact areas. This aects the conduction of heat and electricity across the interface as well as the mechanical contact process.
    [Show full text]