S.P. Timoshenko Institute of Mechanics, National Academy of Sciences, Department of Fracture Mechanics, Kyiv, Ukraine

Similar documents
Mechanics of Earthquakes and Faulting

Mechanics of Earthquakes and Faulting

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

INFLUENCE OF TEMPERATURE ON BEHAVIOR OF THE INTERFACIAL CRACK BETWEEN THE TWO LAYERS

Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens.

Homework Problems. ( σ 11 + σ 22 ) 2. cos (θ /2), ( σ θθ σ rr ) 2. ( σ 22 σ 11 ) 2

Mechanical Properties of Materials

17th European Conference on Fracture 2-5 September,2008, Brno, Czech Republic. Thermal Fracture of a FGM/Homogeneous Bimaterial with Defects

INTRODUCTION TO STRAIN

Mechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002

NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

STANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

Module-4. Mechanical Properties of Metals

Chapter 7. Highlights:

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

Figure 2-1: Stresses under axisymmetric circular loading

Critical applied stresses for a crack initiation from a sharp V-notch

Stress intensity factor analysis for an interface crack between dissimilar isotropic materials

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

Finite Element Method in Geotechnical Engineering

3D Elasticity Theory

Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method

Chapter Two: Mechanical Properties of materials

Chapter 12. Static Equilibrium and Elasticity

Two semi-infinite interfacial cracks between two bonded dissimilar elastic strips

Module 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression

UNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich

Mode III Stress Singularity Analysis of Isotropic and Orthotropic Bi-material near the Interface End

ME 243. Mechanics of Solids

Members Subjected to Torsional Loads

Mechanics PhD Preliminary Spring 2017

Crack Tip Plastic Zone under Mode I Loading and the Non-singular T zz -stress

Linear Elasticity ( ) Objectives. Equipment. Introduction. ε is then

Chapter 3. Load and Stress Analysis

Mechanical Design in Optical Engineering

SOIL MODELS: SAFETY FACTORS AND SETTLEMENTS

Static and Time Dependent Failure of Fibre Reinforced Elastomeric Components. Salim Mirza Element Materials Technology Hitchin, UK

Interaction of newly defined stress intensity factors for angular corners in two diamondshaped

Introduction to Engineering Materials ENGR2000. Dr. Coates

VERIFICATION OF BRITTLE FRACTURE CRITERIA FOR BIMATERIAL STRUCTURES

Stress-Strain Behavior

Module 9 : Foundation on rocks

Mechanics of Earthquakes and Faulting

Geology 229 Engineering Geology. Lecture 5. Engineering Properties of Rocks (West, Ch. 6)

Lecture 15 Strain and stress in beams

Torsion of Shafts Learning objectives

20. Rheology & Linear Elasticity

Integral equations for crack systems in a slightly heterogeneous elastic medium

Mechanics of Earthquakes and Faulting

Chapter 4-b Axially Loaded Members

Structural Analysis I Chapter 4 - Torsion TORSION

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES

Chapter 1 Fracture Physics

Stresses in Curved Beam

VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS

Fundamentals of Linear Elasticity

Lecture 4 Honeycombs Notes, 3.054

Laminated Beams of Isotropic or Orthotropic Materials Subjected to Temperature Change

The science of elasticity

Rocking behaviour of a rigid foundation with an arbitrary embedment

MECHANICS OF 2D MATERIALS

PURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.

Predicting crack initiation in composite material systems due to a thermal expansion mismatch

MECHANICS OF MATERIALS

Using the Timoshenko Beam Bond Model: Example Problem

2/28/2006 Statics ( F.Robilliard) 1

ELASTICITY (MDM 10203)

Basic Energy Principles in Stiffness Analysis

Introduction to fracture mechanics

Module III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina

Stress intensity factors for an inclined and/or eccentric crack in a finite orthotropic lamina

Chapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING )

Outline. Tensile-Test Specimen and Machine. Stress-Strain Curve. Review of Mechanical Properties. Mechanical Behaviour

A parametric study on the elastic-plastic deformation of a centrally heated two-layered composite cylinder with free ends

2.1 Strain energy functions for incompressible materials

Failure surface according to maximum principal stress theory

CHAPTER 6 MECHANICAL PROPERTIES OF METALS PROBLEM SOLUTIONS

STATICALLY INDETERMINATE STRUCTURES

CRACK INITIATION CRITERIA FOR SINGULAR STRESS CONCENTRATIONS Part I: A Universal Assessment of Singular Stress Concentrations

Initiation de fissure dans les milieux fragiles - Prise en compte des contraintes résiduelles

RHEOLOGY & LINEAR ELASTICITY. B Importance of fluids and fractures in deformation C Linear elasticity for homogeneous isotropic materials

A non-linear elastic/perfectly plastic analysis for plane strain undrained expansion tests

Samantha Ramirez, MSE. Stress. The intensity of the internal force acting on a specific plane (area) passing through a point. F 2

EE C245 ME C218 Introduction to MEMS Design Fall 2007

Fracture mechanics fundamentals. Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design

M. Vable Mechanics of Materials: Chapter 5. Torsion of Shafts

Mechanics of Materials Primer

MECHANICS OF MATERIALS

Analysis of the stress field in a wedge using the fast expansions with pointwise determined coefficients

Available online at ScienceDirect. 20th European Conference on Fracture (ECF20) Yu.G. Matvienko*

4. BEAMS: CURVED, COMPOSITE, UNSYMMETRICAL

Lecture 7, Foams, 3.054

CRACK-TIP DIFFRACTION IN A TRANSVERSELY ISOTROPIC SOLID. A.N. Norris and J.D. Achenbach

Brittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm

Elasticity and Plasticity. 1.Basic principles of Elasticity and plasticity. 2.Stress and Deformation of Bars in Axial load 1 / 59

Computational Analysis for Composites

BREAKING OF A FALLING SPAGHETTI Hamid Ghaednia a, Hossein Azizinaghsh b

Elastic Properties of Solid Materials. Notes based on those by James Irvine at

Transcription:

CALCULATION OF THE PREFRACTURE ZONE AT THE CRACK TIP ON THE INTERFACE OF MEDIA A. Kaminsky a L. Kipnis b M. Dudik b G. Khazin b and A. Bykovtscev c a S.P. Timoshenko Institute of Mechanics National Academy of Sciences Department of Fracture Mechanics Kyiv Ukraine b Uman Pedagogical University Physical and Mathematical Department Uman Ukraine c Regional Academy of Natural Sciences Santa-Barbara USA g_khazin@yahoo.com Abstract The plane symmetrical problem on calculation of the prefracture zone at the tip of a crack reaching the interface of isotropic elastic media is considered. The prefracture zone is modelled by lines of rupture of normal displacement located on the interface. An exact solution of the corresponding static problem of the theory of elasticity for piece-homogeneous plane with the mediaseparating boundary in the form of the sides of angle which contains a semi-infinite crack and two straight lines of rupture emerging from the corner point is constructed by the Wiener-Hopf method. The length of the prefracture zone is determined on the base of this solution. The stress near the corner point is investigated. It is shown that in definite intervals of parameter variation the corner point is the singular point of the above problem of the theory of elasticity. It represents the stress concentrator itself. The exponent of singularity of stresses depends on the angle Young s modulus ratio and on the Poisson s ratios. Dependence of the exponent of singularity of stresses on the angle for various intervals of Young s modulus ratio variation has been investigated. Parameter variation intervals where the corner point is not the stress concentrator have been determined. Introduction There is a lot of specific literature aimed at plane problems of prefracture zone development near the crack tip in a piece-homogeneous body located on the interface of two different media [-5]. Prefracture zones are modelled by displacement rupture lines emerging from the crack tips. Considerable interest in connection with its possible usage when the problem of braking of the crack in composite materials are solutions of such problems in an other case of mutual position of the interface of media and the crack in that case when on that interface there is only the crack tip (the crack reaching the interface of media). For isotropic elastoplastic body plane symmetric problems on the calculation of the prefracture zone at the crack tip reaching the rough interface of media in the model frame with a displacement rupture lines have been solved in [6]. Below is given the solution of an analogical problem for the elastic body under the condition that prefracture zone is modelled by a normal displacement rupture lines located on the interface of media. Formulation of the problem Let a piece-homogeneous isotropic elastic body being under the conditions of plane strain be composed from different homogeneous parts connected in between themselves by a thin connecting layer the material of which is more brittle than the materials of the given parts. Let s consider that the straight crack reaches the rough interface of two media and its tip coincides with the corner point of that interface. The region under consideration is considered to be symmetrical in relation with the straight line on which the crack is positioned. Prefracture zone appears and develops with the increase of the outer loading near the crack tip. We will study only the initial stage of its development considering the outer loading being small enough. Then the size of the prefracture zone will be considerably smaller than the crack length and all the other body sizes. Taking into account the small portion of the prefracture zone and the gral being to make calculations we come to the plane static symmetrical problem of the theory of elasticity for the piecehomogeneous isotropic plane with the interface of media in the form of angle sides having a semiinfinite crack emerging from the vertex when the prefracture zone is available. Asymptotic is being realized at infinity. It is solution of the analogical problem without the prefracture zone born the smallest at the interval ]-;[ root of its characteristic equation. That problem (problem K) was VOLUME ISSUE 3 79

ISSN 9-7478 scrutinized in [6]. Mentioned solution contains arbitrary constant that is considered to be given. It characterizes the intensity of the external field and must be determined from the solution of the external problem. In accordance with the hypothesis of localization the initial prefracture zones near the crack tip and other corner points concentrators of stresses are in themselves thin layers of material narrow strips emerging from concentrators. Following the hypothesis of localization and taking into consideration properties of the binding material we ll consider that the initial prefracture zone represents a pair of narrow strips emerging from the crack tip and located on the interface of media. It s supposed that in the problem of the theory of elasticity for the finite body which corresponds to the stage of deformation process when the prefracture zone has not appeared yet (problem I) on the interface of media near the crack tip the normal stress is stretching (condition T). Limitation on the problem parameter enabling the solving of the defined condition is given below. Due to the fact that the binding material is elastic main deformations in the prefracture zone are developing on the mechanism of normal rupture. Consequently the strip-zone will be modelled by the normal displacement rupture line on which the normal stress is equal to the given constant of the binding material. Thus the boundary conditions of the problem of the theory of elasticity which model the process under consideration are as follows (Figure ): E ν r θ l α E ν Fig. Normal displacement ruptures lines at the tip of a crack reaching the interface of media θπ-α σ θ τ rθ ; θ-α τ rθ u θ ; θ <σ θ > <τ rθ > < u r >; () θ r<l σ θ σ; θ r>l< u θ >; () θ r σ θ Сgr o(/r) (3) Where α θ π α; a is the jump of value a; σ C > are given constants; is the smallest on the interval ]-;[ root of equation ( x ) ( b ( b ( e b ( e ( (sin zα z sin α ){( χ) 4[ χ sin z( z sin α ( ( χ)( χ)sinzπ 4( χ sin zα zsin α )[sin z( z (4) b ]} b (sinzα zsinα ){( χ )( χ ) 4[ χ sin z( z sin α]} sin α ] ν E b ( 4( χ sin z z sin α )[sin z( z sin α ] χ 3 4ν e e e ν E E E are Young s modulus ν ν are Poisson s ratios; g ( α e ν ν ) is the known function. The values of (root of the characteristic equation (4)) for some values of α and e under ν ν 3 are shown in the Table. As shown in the calculation results function g ( α e ν ν ) under fixed α ν ν is positive if < e < e ( α ν ν ) ( e > is the point where the function is zero). Considering that e satisfies the (5). Then in the problem I mentioned above σ θ ( r o) if r and the condition T pointed above will be fulfilled. Condition T is fulfilled under all α ν ν if e is less than the smallest value of the function e ( α ν ν ). In particular under ν ν 3 for the fulfillment of the condition T it is sufficient that e does not exceed three. (5) 8 ADVANCEDSCIENCE.ORG

α e Table. Some values of the root of the characteristic equation 3 5 3 5 5-4398 -4679-4793 -4868-583 -5338-56 -637 3-434 -4483-46 -474-5383 -5688-66 -69 45-45 -435-4476 -46-5588 -6-664 -744 6-489 -488-438 -456-574 -644-688 -7665 75-386 -397-473 -449-5795 -635-6957 -775 9-3467 -3757-466 -437-5737 -6-679 -7546 5-34 -389-45 -445-5594 -5969-6475 -79-479 -46-445 -4636-548 -567-645 -6647 35-46 -4679-4748 -484-57 -538-563 -599 5-489 -496-494 -4945-587 -553-554 -546 65-4987 -4988-499 -4993-53 -55-547 -589 Near the tip of the side rupture line taking into account general things about the stress behaviour in the vicinity of corner points of elastic body the asymptotic is realized. It represents the solution of a homogeneous problem of the theory of elasticity for piece-homogeneous plane which contains on the straightlined interface of media a semi-infinite normal displacement rupture line born the root / of its characteristic equation. In particular χ θ r l σ θ ~ e χe ki (6) ( χ e) π ( r l) Where k I is the stress intensity factor at the tip of the normal displacement rupture line which is to be defined. The task is to define the length l of rupture lines which model the prefracture zone and to investigate the stress near the tip O of the crack. The length l of the rupture line is defined from the condition of the stress limitation near its tip. Solution of Winer-Hopf equation and defining the stress intensity factor The solution of the formulated problem of the theory of elasticity with boundary conditions ()- (3) (Figure ) represents the sum of solutions of the following two problems. The first (problem I) differs from it in the following way: instead of the first condition () we have θ r< l σ θ σ Cgr (7) and at infinity the stresses decrease as o ( / r)( σθ o(/ r) in (3)). The second problem is problem K. Using Mellin s integral transform with a complex parameter p to the equilibrium equations compatibility equation Hooke s law conditions () and considering the second condition () and condition (7) we come to the following Wiener-Hopf functional equation of problem : σ σ Ô p Atgpπ G( Ô ( (8) p p ( χ)[ χ ( χ) e] G p A G p G ( [ χ ( χχ ) e χe ][ a( a ( e] cos pπ [ χ ( χχ) e χe ] G( G ( [ χ ( χ) e][ b ( b ( e b ( e ] sin pπ a ( ( χ)(sin pα psin α )[sin p( psin α ] a( ( χ)(cos pα cosα )[sin p( π α ) p sin α ] p E u σ σ ( ρ ) ρ ρ θ p Cgl Ô p θ l d Ô p 4( ρ dρ ν ) r r θ ρ l Where - ε < Re p <ε ε are sufficiently small positive numbers. Simslar equations are solved in [5-7]. Solution of equation (8) is as follows: VOLUME ISSUE 3 8

ISSN 9-7478 pg ( σ K ( K ( ) σ { [ ] K p K Ô p [ ]} (Re p< ) ( ) K p p pg p G p pg ( ( ) G ( ) K ( G ( σk ( ) σ [ K Ô p ] (Re p> ) A ( p ) G ( ) ( p )( ) G ( ) i ln G( G ( Re p< ± Γ(m exp[ dz] K ( πi z p i G ( Re p> Γ( m ( Γ( is gamma function). Using (6)(9) we can find the stresses and stress intensity factor. Stress intensity factor at the tip of the rupture line is represented by the formula ( χ Γ e) g ( ) k I l Cl σ () χ χ 3 e Γ( ) G ( ) π G ( ) Results Equalling k I to zero (see ()) we get the following formula aimed at defining the length of the prefracture zone: C π gγ( ) I ln G( it) ln G( it) l L L [ ] I σ Γ( 3 exp[ ) I dt I exp[ dt]. () π t π t ( ) / Đn case when ( C / σ ) weakly changes with the changing α and e if e < e and e < e 3 the length of the prefracture zone will be the most if α α m. Đn the above discussed case the length of the prefracture zone increases with a decreasing α if e e e ( e 7 e 45). Đt is known that part of the prefracture zone located at the crack tip represents the destruction of material region where the level of stresses is extremely high. The given region differs at its highest maximum level of deformations resulting in availability of pores and microcracks. Đn that connectionparticular interest lies in the analysis of stress behaviour near the corner point O i.e.in the study of influence of changing the problem parameters on the exponent of stress singularity at the corner point. With aim of realization the above mentioned investigation using (9) and Mellin s formula we define the stresses in the problem. Đn particular we get F( p θ ) M ( p M( τ rθ π r dp D( a ( a( e M ( i ( p )( p ) D( Γ( γ M( p D( a ( a( e M ( M( p[ m( p ) m ( p ) ] Γ( G ( l Γ( gγ( ) m Cl m σ () Γ( 3 ) G ( ) π G ( ) where - α < θ < π α; F( p θ ) is the know entire function of p having zero of first order in the point p ; γ is the parallel to the imaginary axis line that lies in the strip - ε < Re p<. Function D ( has zeroes of firsty order at the points p and p where is the only on the interval ]-;[ root of equation D ( x ). Function D ( has not other zeroes in the strip < Re p <. Thus in the strip < Re p < the integrand in () has three singularities simple poles at the points p ; ;. Using the data about the singular points of the integrand in () applying to the integral () the residue theorem and adding the solutions of problems and K we find the principal terms of the expansion of stress τ r θ in asymptotic series if r in the problem of the theory of elasticity under consideration with boundary conditions () (3) (Figure ). The below given formula takes its place: τ r f θ α e ν ν ) C f ( θ α e ν ν ) σ f ( r θ α e ν ν l σ C ) (3) rθ ( ( ) ( ) ϕ α e νν σl ϕ α e νν Cl C (9) 8 ADVANCEDSCIENCE.ORG

( f f ϕ ϕ are the known functions; f if r ; < Re p <. Under the values of the problem parameters where the equation D ( x ) on the interval ]-;[ has not root in (3) addendum containing r is not available. The result analysis drives us to the following conclusions. Đn definite intervals of changing of the parameters the corner point O represents in itself the singular point of the considered problem of the theory of elasticity. Đt itself represents the stress concentrator. The exponent of singularity of stresses ] ;[ depends on the angle α Young s modulus ratio e E / E and on the Poısson s ratıos ν ν. Conclusions Thus the plane symmetrical problem on calculation of the prefracture zone at the tip of a crack reaching the interface of isotropic elastic media is investigated. The prefracture zone is modelled by lines of rupture of normal displacement. An exact solution of the corresponding problem of the theory of elasticity is constructed. It is shown that in definite intervals of parameter variation the corner point is the singular point of this problem. It represents the stress concentrator itself. The tend of the stresses to infinity corresponding to that part of the prefracture zone located near the crack tip is the region of material destruction where the stress level is extremely high. References:. Comninou M. The Interfase Crack Journal of Applied Mechanics Vol.44 4 p.63-636 977.. Comninou M. Schmueser D. The Interface Crack in a Combined Tension-Compression and Shear Field Journal of Applied Mechanics Vol.46 p.345-348 979. 3. Gautesen A.K. Dunders J. The Interface Crack in a Tension Field Journal of Applied Mechanics Vol.54 p.93-98 987. 4. Goldstein R.V. Perelmuter M. Modeling of Bonding at an Interface Crack International Jornal of Fracture Vol.99 p.53-79 999. 5. Kamınsky A.A. Dudik M.V. Kıpnıs L.A. On the Direction of Development of a Thin Fracture Process Zone at the Tip of an Interfacial Crack Between Dissimilar Media International Applied Mechanics Vol.4 p.36-44 6. 6. Kamınsky A.A. Kıpnıs L.A. Dudik M.V. Initial Development of the Prefracture Zone Near the Tip of a Crack Reaching the Interface Between Dissimilar Media International Applied Mechanics Vol.4 p.76-8 4. 7. Kamınsky A.A. Kıpnıs L.A. Khazin G.A. On Criterion of the Start of Two Shear Cracks in an Elastic Body under the Plane Strain International Applied Mechanics Vol.4 4 p.8-87 6. VOLUME ISSUE 3 83