Thanks to visit codestin.com
Credit goes to www.scribd.com

Scribd
Upload
60 views4 pages

Elements of Fracture Mechanics: Birla Institute of Technology & Science Pilani (Rajasthan)

The document describes deriving an expression for the approximate shape of the plastic zone around the tip of a Mode II crack in an infinite plate under plane stress conditions, using the Mises yield criterion. It shows the derivation of the expression in several steps. Plots of the plastic zone shape generated using Matlab and Excel are also included.

Uploaded by

Gopinath Chakraborty
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
60 views4 pages

Elements of Fracture Mechanics: Birla Institute of Technology & Science Pilani (Rajasthan)

The document describes deriving an expression for the approximate shape of the plastic zone around the tip of a Mode II crack in an infinite plate under plane stress conditions, using the Mises yield criterion. It shows the derivation of the expression in several steps. Plots of the plastic zone shape generated using Matlab and Excel are also included.

Uploaded by

Gopinath Chakraborty
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

ELEMENTS OF FRACTURE MECHANICS

Derive the expression for approximate plastic zone shape around the
crack tip of a Mode II crack in an infinite plate for plane stress situation
using Mises yield criterion. Plot the shape of the plastic zone.
.

By
GROUP 5
PRAKHAR VERMA(2018h1410087G)
VIPUL RITWIK(2018H1410183G)
GOPINATH CHAKRABORTY(2017H1410060G)

(Second semester, Design engineering, ME)

BIRLA INSTITUTE OF TECHNOLOGY & SCIENCE


PILANI (RAJASTHAN)
Solution :
   3 
  sin 2 .(2  cos 2 .cos 2 ) 
 x   
  K II    3  ....................(1)

 y  sin .cos .cos
2 r  2 2 2 
   
   cos  .(1  sin  .sin 3 ) 
xy

 2 2 2 

x  y    x  y 
2

 1,2        xy ............................(2)
2

 2   2 

x  y  1    3   3 
    2sin  sin .cos .cos  sin .cos .cos  ............(3)
 2  2 2 2 2 2 2 2 2

x  y  
    sin  C...............(4)
 2  2

  x  y  1    3 
    2sin  2sin .cos .cos  ...............(5)
 2  2 2 2 2 2

  x  y      3 
    sin  sin .cos .cos  ...........................(6)
 2   2 2 2 2

   y 
2

R  x   xy ........................................................(7)
2

 2 

   3    3  2 2 2 3   3
R  sin 2  sin 2 .cos 2 .cos 2  2  sin 2 .cos .cos   cos  sin cos sin  2sin cos 2 sin ....(8)
2 2 2 2  2 2 2  2 2 2 2 2 2 2

    3  3 
R  1  sin 2 .cos 2  sin 2  sin .cos  cos .sin ................(9)
2 2  2 2 2 2 
sin 2   
R  1  sin 2 (1  2sin 2   cos  )  cos 2 1  2sin 2   cos  
4 2 2

K II
 1  C  R  ....................(10)
2 r
K
 2  C  R  II .....................(11)
2 r
 3  0.......(12)
ByMises yield criterion
( 1   2 ) 2  ( 2   3 ) 2  ( 3   1 ) 2  2 ys ....................(13)
(C  R  C  R) 2   C  R   (C  R) 2  2 ys ....................(14)
2

4 R 2  2 R 2  2C 2   ys ....................(15)
3R 2  C 2   ys ....................(16)
𝐾𝐼𝐼 2 𝑠𝑖𝑛2 𝜃 𝜃 𝜃 𝜃
x [3 x [1 + + 𝑠𝑖𝑛2 (1 − 2sin2 θ + 𝑐𝑜𝑠θ) + 𝑐𝑜𝑠 2 (1 − 2sin2 θ + 𝑐𝑜𝑠θ)] + 𝑠𝑖𝑛2 ] = 𝜎𝑦𝑠 2 11
2𝛱𝑟 4 2 2 2
....................(17)
𝐾𝐼𝐼 2 𝜃 𝜃 𝜃
x [3 x [4 + 𝑠𝑖𝑛2 𝜃 + 4 − 8𝑠𝑖𝑛2 𝜃 + 4𝑐𝑜𝑠θ(𝑠𝑖𝑛2 − 𝑐𝑜𝑠 2 )] + 4𝑠𝑖𝑛2 ] = 𝜎𝑦𝑠 2 ....................(18)
8𝛱𝑟 2 2 2

𝐾𝐼𝐼 2 𝜃
x [3 x [4 − 3𝑠𝑖𝑛2 𝜃] + 4𝑠𝑖𝑛2 ] = 𝜎𝑦𝑠 2 ....................(19)
8𝛱𝑟 2

𝐾𝐼𝐼 2 𝜃
x [ [12 − 9𝑠𝑖𝑛2 𝜃] + 4𝑠𝑖𝑛2 ] = 𝜎𝑦𝑠 2 ....................(20)
8𝛱𝑟 2

𝐾𝐼𝐼 2 1−𝑐𝑜𝑠θ 𝜃
x [ [12 − 9𝑠𝑖𝑛2 𝜃] + 4( )𝑠𝑖𝑛2 ] = 𝜎𝑦𝑠 2 ....................(21)
8𝛱𝑟 2 2

𝐾𝐼𝐼 2
x [ [12 − 9𝑠𝑖𝑛2 𝜃 − 2𝑐𝑜𝑠θ]] = 𝜎𝑦𝑠 2 ....................(22)
8𝛱𝑟

𝐾𝐼𝐼 2
rpz = x [ [12 − 9𝑠𝑖𝑛2 𝜃 − 2𝑐𝑜𝑠θ]] ....................(23)
8𝛱𝜎𝑦𝑠

SHAPE OF PLASTIC ZONE

Figure1: Shape of Plastic Zone in Mode2(Matlab)


Matlab code:
theta = 0:0.01:2*pi;
rho = 14-9*sin(theta).^2-2*cos(theta);
polarplot(theta,rho)

Figure 2: Shape of Plastic Zone in Mode2(Excel)

You might also like