f\4t3+z

Tt*g

I nrnd,vulLitrn

Cau'

n/nd Afph'co,{r'an4 I.nr (craf,

aq_
4

Wei Cai
Stanford University
Spring 2012-2013

fu^".J, ih apphauhrvv',
1,v,tf"; uv\r.rW,tnl0ilSc,w54
plarsnak;,
atnd {ra.+-nl

Sfu{ ,+ td"c,{'ta;t*uh^il*- it\ *tine*aY yn**,,m-ttrelv*rts

0./^vunds"tsLr.ndy pt"at'v\t rh $-,.l| "l {hs crc.vYfif
"+
TI',i^cMM{ is thg seqwqt' j- ME}40 Ffrr^,h.^>,t\r
we
, in ,J,;c,t,rtho.i tl'* :tt** * Uel;Wd
/rrruVtt{, t^od 5,^ch
sirt6 gnry v;sal .
Ln tl"a CA'r,rM{.,r^Jt_
oJU,r)l,r)
ti\{-. l,o t, lpau{r, f w excee.J
) .yqtd S4,Rv)
"d
an&affflN th{" sfuc,'+wt-qh dt"ar,&Jftf
tpt,mnn-{.^ru{-;
yl*,1.2 s*oin .

'Yfuv+t

7h.L 3u^'f-r'5 [Sh((t h f"il tl.w Slrr,A,

etatt^ttt-Plooht'

60,-^4*l

fL t,

Jkar;- , arno!dnplo,t4n"u^lf f;eldz

N ,owuf "f t("+ sh*,-+tn-n1,,{,&u*trdy.,
U' *
^P'4'!*ai[l' wnduyo ylx*<- *f'-'@,

e{^t*t-

ry*

tfu- ru>,+" yp/rrto-'r.A irt th.q, O{u+ra

fo

04 &,41inJttwrnLffi,

,"irr^,

h]e u-Uo

tY,

Wo{ t" c:lsf?r,n;y\eefu g{ar*(-plm+Z n*Uy

r') t:h&t -@ tU
TlttCa-+rrn
n'ar;4 i-s ttL |illr4
r ' - ' 0{
ol'atfl
u'

TfrL othuv ct
^e_{"&

+hril,f th*
Nt Ehoil N!4a*&
h tM/rJ,
i^

WL^L

pla*z
4,nd

ry*,

-/
Pt a-Ya+L

sdroJ' 'n .

otw>lqZ :;*'a')n ,'-) prvpo-r#*xl

B"l- 4hL plrv"/sctta,i--

tuC

l-d_

Sl,'a,v2 or1

npt yrvVtt*n,nJ h stve+'t. 7A.tflcr>a+cStra,,*t

wn/, frt-L firt"r.\ *!*tl

nT +hL /6)+*'J-'f t,,ni,t6-.

j\E:W

:, CaA'

/n+rcUn,r'cfrr*r-

6lo^I

(frtv\,'-\ i

ilul

"/

t. Tlte 4,*T

,l- ft*d,z;irr:p

t',e,rlr+plursrl.^"a
lgh*til -fl,* yh-ut"",rrnt
na a^d
A"nr ft o(Uw;l,e if m*lJv"'n*4/ CP>e)
( nof in,{,'uerY FEM )
2. Me+/'*J- ,f .-!rhr*t=u-,'r,
A-r'L'W4;c"!
Ienw)--a^#'

vnuQr^d
c,J

r -glr{wsl,,>vlnr y

M*{*

3 . 5ol-,*o" + c(urw;r',i-p{^*++^fy pW{yl,t,^s

A ^y* go& ,f fl,,b c,tutvu-is h oLrra/

f^n
/ |,,|*;livt
,j
, )o +h^f
fr"A'
/w+wVV,nn&xna
" Wlnurt W

vnC,vun'v1rr'A t,^,?b pl at+-,Zrty prv\l*rtt

t, c6/n uW tun wnglitX+C+l"

W.rewly,
Y0
!illw
ttz ^ S;nYl"!'n\wlJtic*L cDd'g h
Sol^*asn
, fr *
I"t"t it
in

"

,''f.r^

wc;y,Wrr(et, &n fty^h dym-a/wt rAJfit^I

cttM) !uLo*rt".W

cowvw"ry

bil

br^rfrhu^.

+lv- go|^Arn / vv1rl,,l) gUnt\4/.

RO(oarct*U,p ^n'fut^ o+htL c-uwt-tt'b

'E{ou*z;%'
' N1t )(0
( w egwtvdtentt^A
frunquuih'te )
/'fuy*+ar+r--tJ- A4^-C,r-,r.r.
' Ce?4U
r-Karrr;u\"
Sgpti$ ?--t9-)
{*rl^

*rfh-

f-o rata+orv

tlt+ r'nzy}neP@
\h*"^'"u-rvq
o'I7t"tftn
na*mut)cJ

Tltr) o{'v-rt^a f ^*

(i1. t+'''
lt'"'y"'t

'f +l^r p;+r^,

-p/'Au*a+r;,

M*th,l

-'q}fud'

I'4E14,L

/n*t &w,r&t.'1

*Wfr,
F*" E7to41i?k4

u.rsLh^,

= K.-t(r*-)
r
rnao{

h': (a-zr) b

Al

Shry4 anncann4va&ic,r
f*^nnz

u,f^!lv s1rp/nB b&,N ,) (trrtd y,?,n)

,rY,\$-!/14

lo

K== J

a , ur/r*

*rt,e"^A>>t:

tfu- >r,u dr$*g"*.^6-rn,

6Aeq-

@*).* v 6v ?
A, Ar*

",1\

a +h!. ry

"/ fAL pl*ttz Bo,l ?

5;L'l- LALL o"r-oqk

ffnnLtnhT
{l'^C Pndiv0
f ^-;!,

\r

,.hd\hrdf
ruit{q^&+do-1-+*.
'uU*a',;lg' srn^[dF.

q,, #fY,^.r
i-s acc-+q,,.,rrf,ed
f"- ?

(d. whatyr's fhL "sV {;hye) "f ile

pl at+e' ;"rrt{, ?

cy 6/v {^rd.i**
d., ,tAho*;t thL
') "firfr I
lt thr,un

ilA Lu*"'ruis ta qqfih.art^I/ q0t0"h'0r14

*&'ry

Th.!-pwyow (

*9

Nirh, 6

fr^ct,ourrurwf'n-{'
fuual{Yt Can'

l't483+2,

J,

lo(un'c ,

isut'.p& mtcJiunq,

hl i-tl+ br;(

l-prttnvh-rrts

u'nL d,lsCtwVW

,
T/u, tAfa ,f rlo&-ou ind-,adp'equ'4'by,'wrrL
ouzJ"^hl,n,frn sfit't's ( FqI
. corvf
/"1--1sfua'; ( f s )
*'b;L-tr, and;lxn f*
' SlMt'Sfrc-n

relo-h'rvt

' ytVll c,md^

l'rrvt
' StTUy)- Slta;a

t'n Ela'ut'z

lui-;\'71L ( { f ;

( f l* .r

Ytbh-tn-

CJtt"^

,rr/-e-) | rt Ptorrt-;t Ry,rrt

( gI )

A Swmrnuy"f 4 1u'r"wtst ,i5 yw.n th gtL.

viLa^)lt' a. 3,nyledl-)plMW +t*y:,
Tl,,e-go& i.s h, prwfu & 'lobauA
'

MWg /^'5c*t,1<iu'nJrrL
M,

Nitl,L aXwm7lzi,

A p"tt*V

t, k*p

l4ieH ,orrthVn- und l:[r-ru r"'/-s- wi(( be

in scabe(ul"vt

le'uhne-t'

( t, '(, n0 haadlrtnJ )
plahfc 114&d,4'tvtw-

+he disCt,,*ltv, tryb.

d^:c-^^ad Laf,n .

rf

a,uwns/

l'lerte-

Vanu'itu'lLr*rle,n;rtj ru/-et w;t( bL

ilE)+r

f*Ao"*'d""9,Equ"'Jtna Cu'

elru+atV., h &L tlry

"J ptaz#c;tn

Siahn t" 4q- tuI

uti'lv sfi%J' s'traun a'nd dnplacune-''*

utg wil,( al:,o be r/4^b/*

iJ +h^hdl- s'|rst5-s+fr-^,ot",l+ntslgis
$o6. fnt na;. l;ft'a.o*tLa(^.1, h;Xn1-/Lfu,J4."t) n tho #*ry
nn,a-Linea'r
"f pleft.
9t

Dis/aa-*e*

no( d
U :
?isplac*"+ta.-e,./.
'l
*

X.- - X
1a,/

l"isft**'r',rn^*
fctJ ! c*l
wo w;[( attT'vvrvwLn +hil cJo.,,us,,
I Al<Z I *'d

A(X)
i1n*- tluu'd;fiu'wrra J'>ezfineonL c X t a , n d
dxyh-r""'^'t {z6 ln owyaat'*ft
uxoi''x'y,Z)
u/

+o*,

ti:,:)(
t Jt 3 t

U? (X,Y'?)

9Z

Ska'i.

F-e&t

ei =!ffi ++.) ==tQt,,)+ui.t)

e,1, xx: U+x
2*/: *-(ur,y +ilr",*)
ThL uW *tY
4 c*1^fut.un*

grat)n , e ,g. fr/ = 1'lx,/-,' ily,* ,

,f

il,u

5f-72uvatpnao-r,

tr'r)''=1.2.3
x,
lt:
dz-: y.
Xe : ?

is not

ME)+2-

Fto"dw^ar*tri

Cao'

ryttn'+'Lrrui

Ji"e,n fiy {d

9l

,<
o;i
-

t -,

u/yt^t Aweq ar1 i,t{,v fzcl,

J-6.rcL
W.
1

<{c'wa'll.Tn

i/l ;-tL
J

G;rrrn,,th+Stya od- yv;,rt ?,

!,

T fqlL

nvL
wrlA+ dJL"s.6,

ry
wllh nd.,;lr*fr*rle;'ffU L iJ

21 ,

f,rlT"&

dpntp-,t
'

a- vel-M

{r'

5e cond o>'c{l-1r tt{tlsrJ

{synme+lsz )

dttt

T/r&t,t*yozt**fra cA"-ye

lt ^L6 tn'w:fwn**sn

Ltl: ai ,F
ej : a^'fA), ?t

*/,.p"J-

coctdn^rt6-*aSy:)+-n a4

Q.-o=(e"'.
[

Qr)

^
/

U,nal- ve*trt
"f nz^l

t7': Q"'rQ't-Tz

wmf

r,^s.trbt*ztVl"r,
"J { . J
"{

I
c.*rtlln&*"

coaA:v'e&SySI+,*>

'iJ]erl

( W frL{w d,t*"Ju,

4-

ME3rt.Ofe,'A^^r AbjitJ t/Jrnte^r2-al3. )

fq,-L.tti*ry- fr"$+.r{3t!,t

jii.t/'l

"

f r 6 ^ ^ + r r o y * ' r a z 6vx
a
aV

- L - ? - / uyl
r'
o' (x.,:z'T
TJ

h(ate;Tl'c.t7w'Abril,rr"taond;+;s+a ra,h; frrd

'
ul

il/r'if

r.idPtttl

A
,
-E u ) ulyt
alt uxy
y

,nsrt+* ^L

PtiL

Y'cr.-F:o
)

C.6?nftu\*,"-f

body $\"a

-..

J- tj:P

t {4 tX.-aSrt' [o:lzr,ktvL
in

f""*

= rry f7i

5
t(r' i5

$r. Jt^vfrn

d at>+t

o-r'' pl a,,td<-.

Fx:t2

'
2
1- Az oul

'y:o

+- a *
>*

-f

F* :o

reymc{@t .J
r

uA#!,r,y th<

31-z
ME=

F,^dwrrw**tQaat;a''l

C*;

t-t' ?,"Yq:!'LV ryd{fr:*Ii{*ie4.-s @n

ei : * (u,,i + uj ,, )

a1:

s4"-c*'-"

is dra- {u+'!

ef + {'

u#"

P:ffi

c-md.^{t.nr-,f +hz ^"/vrd

ho,t t" S^n4U +l* oowy6a't/,;A'hg

e^;

)'J nof

rwy+tz'i"d

--E,.r,jL
UIwq,ufrEt^Wannd^'{rrr,o,' J ; ' l - " + ,or,
- " ' ,J
)
'-lJ"
6
fu,xt.
|
'T1r,,a,
c-crvy6gp1l2,1^.1-?
c"-d;-liuyl- go'.rane.rctz4f +r\^t
Cwnt-

talu+fe,"tin-Fhz'rr

TlLe- *Y*lr$%c.,md^'{+sa
ridrut@
C fu

u;

ttt*-e,- dt +"-fu ,

Vote. lL

*:?*+nL&n-o"trcf*ea
ts

a*,d o6+o;
%

1'

o,r*rrrn*-*tr
g:i

a/\ *C dl

J'1,;x

- c2

( 6 OoF ^r-al
,f

porr',tJ
:

Ut C3oo1s.

Sa4$sfed tf

+ t-lr-,;) .

2-ol) ) .
ME 34O IRtAML,\./o+{p (,\.)tvdt^r

8+?r;. fiuld tl*

eJa,y++c

dou rvul ha,*tat^s^rrfW

*^f *h(^% oorJ,-+.:r',

;J e_,1.'!
;
"
*

I,Ll to +h^APr,\t , W%

r.d tLW SWna-0d rh etar*:A+-,6+ttuvt 14

eqwG*1'0n4,

( Pruv$qJ.,$ : o).
i

I,AEN7
F4,4dwryt&
e,yfy^.
f o Ta,vr;le-S*rel's- {+/a}n Cu.n'r,v

Cori

!':.:,i'M^S&d/',31nArt
+ a d^ sblg ,", "r{-wtoi.

\s+'"^).

t'n & {fi*-Lz

t *.t;"9

{Ut

N WW

Sivw-saro,tn
"yd

\ nn slr.^^;n

( .u-eoL'aas6a )

plannv rncuFvu{

{Orrdr.","it3

:
l
i

7noYg- arrrylerl ( veo"h'sf;c) t

Ng. NiA dn^Su,ttv3

m,od,rLs c1Lylatth'c nov*u'w-L b4q,

Fy. y ie{l s'lruvs
i,s;,
Motx-, 4f*4hL |ie{d
d^!-f,tyw-\,tm)n cmwf,
Vot'ilr reacA-e-d,
berff U

gT furffi

U*t

ara'P'e-a42rtt-'

(SY-tfu.,,t (el/'&vnin t/a Ela*v P$tnw

year/na4'!1v
(
lf ob s-lrt/vt rr.o,uerv
Ttvld arrcL-+rw., t'g tr<6y in W;/aW)
d4-p/Lt/,e S+rvlj-Sltar"t wb4?:l,Lr5 {;nso,rr,*nd

Qp^ur-!;as4 #*k-ph taAil.

6-i :

lft+Y-

irt"f*Aer,fr.

C'iirt Lff

,l."

lo'*+ic S$#Nt''A
4+aa'v'v

zj:

7sotryi
c e,(o'*;r;+
,

irlo,

dt--Q

Clrt - 7 E,1'Ert +)^ ( fir 5, * f,'r{r< )

6iji: ) trr li:j + 2A t).,

tll ,

the^u

'/t4"od.&oa

)) , Part*a

ro*fo

7-#

:
'

ME3+z
I-n

1\' v'a'|ivAA

6r"n d a,'nPrr*oJ

Can

&mFffie/^Jf

6,<x.= ()t zy) frw t )?qv t

6ry =_z,F f^y

) l+'e

A & " t ( l + z / ) L y y+ ) h t

/f

"YY
f++

2"t t

l oxt
-

V (T.
'-YY

Y 6ut

-I6xy+

-:6yo

+t

- + 6 . n / Y 6 y v+ + d z t
,}t{=li*

f*x:

)
E = =z , p \ l + ) )-I*

s+v*)
hld^^t*r"-c Strl{, ( rntaunn0l'rrr^.1,

f6,

t c-+t)
+kn*a6yy

-:6 * + 6 r y + f z z (: I +z/ )

l
t

c1(y

t3o*,ce.6^,,"/,+
z

trz: # 6Yz

6Li,

)z

[ry=F ez ,

t t

:-7)^

G"=

) )v r LlrlA)[o*

tr

Q+

dyv

*d, &Aa

s+rar.^
Gtorn

t&**.

[ = i L ,- i (

''',i-ll
+ &* )

Lnt4y*f<*)

tr

-8",{k

r_-

K ' =t(rx+r7tr)=
4#r:

t1orar,j g-rc,/\^)

3(t->v)

nqirnutt

o = 3 <t
C/*-nrfa,rfo-rt'C S+ratyr te"nZ nf

il-q,u^i6I'ic S+re/vs tt/ruAs-r

Drhw

S,-, =

Jxn: 4," -E ,
\n,/ -

nn/

e.j= [,j

tr6J

fj

<

) vt Vl =

Uyy-"
' t

(^
)ys: D/*,

/,

)'3* ::: 6 r v - 7
?

*'tY

c
1 wg-t q{*+\-rvL bq,trme'n dooi5$rrv{ Stx,t5

for

-'

r/

6rr
ernd

"!.J

( t x x =[ " , - f ,

ole,vinhMl

E3'r'
y y =t r r - { , r a : f . * * }
Yu= ty =

lxl
Stra,vn iJ

"* "8

);ir"irll

r. i = l, -7, 3
)

+;uny
M th eta^A-'2"+U
h,lrLt- Q^-l-oJ:o tlu- 501,1A!t O6vVU\t fhL equzutrbvr,n
,+n ,rr,'/*fo
n;tL

6rr d*"ub'tue Sh-vsl-s-lnvn t'\ W

lr'p"* h

d-etu'a*c-i t.),
p lottm-c- Sf'.a';r^' ( uoLwotv'iJ' "'(t{

bee-

an/fA'aL-

Q*= t+,x

-7

I'qE 3+L

3s

Eq w r',{-iu.l'tj

F*,'.'cla merCol;

Yie\d

fta

C,c'rdt f,-un

fu-lo,,"nl;on)
( ,: n&tt* tf tr:lan*ii
vlz posnJ*e-&*,t tl,r, tlie{ar*,d,;}"s}1
n,,,Ls*":'v\ {q"'t'o{
ce4\ l;te- ,@c"itd rl bJ A ftt"+frT'u

+(/tl'];

:o

At +,tr1) foird, f

( r r x , u v y , ( 1oty: 'r,, 6 2 ^ , 6 r / )o:

J
;
:
"('

'
r) ^ 3er"tt*l?- 1\*"(,f.,tr?L

t'r'=#\)

M o t z t h ) ,tf, n t - o ^ d r i t * . . s( r i , t = * o j ,

d - on o t

onJ-wf

tJ ^q'(1i*r$ed bJ
Ad"
'
^'j^
idz'Jiaohun
Yi/(
Jo ttu i-r
^J''2 l'p' ''nuh-U '
nrylr -^t f"l- sfuzrrtovhwv'; J^*
tl-t t^k^t "l {

-ilr- wia.;te;ri*L iJ isolrrrp

iu, +A-p-',
+hlt
645u-m-eS,,h^-!,J*
e
W
A cnt',,:-rlirLa-t"L+-'z'naf*n"rltrn , l.e .
b'1
'l^*'f
nat
5hlr"1"4

(
1
6r' ) ),,,,, 'n"/os"-t tJ'
;
{
\
}
)
:
t'
f

-fh^\

hJe rnn,,l b" '''/uf'r'6

rr&-attti +hsf

S4tvt: iitv'Nia'att

"f
fuqf,ft-''t'J
Cna-&t7',,r^f-\

6L

A,if A-it fif.

in ++/wU

th^t /'ou nnf fuuW U

franAftunorA\v1'

-n&,,rt\ **nnJ,- s*e-n

f* uXao''ytn,t{'o-

l,(.. 6: t

UXlk'

E' =*

5,';

tJ asfw

invwt;w't

6
..: + 6!"
" t !
)

in kJ?t"tuJ4
bt0+fu:rr s/tt| t'nvaa'^'a/^t4C-a/^
"ttuf'wl
' q ' 6> ' 63
pv inc,ya| sl'v"s vo'h'rL$

tl4'L

m a s7+c''^L
I
Thul a,rt- thz nvan,^'(-s']i+*r vaJN4A
J"""*nuf

rJ,g s\,sknt
u
a

st'rJath^t

'A

rt^alh
ske*'' 5fi-e-+re-4

Nlo*t^r^"+,'*lt ,
Stlt/vl n#*X

v l t

d
a1Lr

c
\
), G ; i \
J )

#it fli,Se,ovatu+;
cI: o.A-e"f +i\"l*
5xl 6xr -l
4$^
[

I*
i"e

&^rg

6r , 6, ,6)

daf

,T' Syv6y*

4x

IL rt - i i i Jl

rrr, - ti

l-3

cir,+
r

tf i
oqwc-*irn
Pd1,1,'ot,in-8-

JrX 1 - h

wQnq."s- ar , 1-.,

u nl

+kL ,lj *. -e4 'at+rlrn

ut;,. 6yy-) 6'r. l

I
I
-i I
I orr, 61 ,ra*

Tl'\L *-S*-etiAAffi"\ DD

6+y dtt J

th.c- t hrr: Jol,q-,tirrrul.+

{
l

Ca/'

Ft,.J*t'.{,t{cL fo ttffia'5

MFYII

(zt

ll

a/rQ- 5;,1"J-re^4
lhuo-nin ^fo1

-e1r'Ln-tiraiart mr,u1.hrvrt

Btc".*rc- 6r ,6>,6-1 one soL\*ikv\a"f +"" dt*

&qlti-o'l(1-o'):o

i - 6rd" 6l : o - --- ( 3)
il - (qt o-"-ttr: ) i - (o'o"+ 6'6: * rtar )

^f*ry

a,,nd
4rl), y3o'f

%rrt

Thrvep s*ren,h uen^an'b ;

(lr:6r-tfr*q-36

Lz=-S, 6r + fa trl r sr dr)

I-3 :6r

Ttg. sfu,t

covqolaoW c'n'';D"t'-rJfn'*- l,t)

irtwuo,,,^bt

Oa>dqvfr

c-omVu^>,1
\tti
'

D}ihcrp^L s{f/o vtrJ-o:r;5

J

tBtt')
lvota: io-no boolc: c]""rv12"F'F"t"'tt J /
of
in tl"t ds'f'r^n*+icrt f'.
\

63

fr

*{ tlu

k)L*^a

a-m n'b;"^'1

L I
(
J'L,,l't.'tStnva:eA dn r,'ot v-rurrrrJ
SL,).4&e-*t
f"l Wt'A"t\

"^r(

7 ( t-l, *n

!tJ'

l'I,

-L,=l6vt 6yrl + lonr 6r- \ + Jo,r5*/ I = Q*6yy+dyrrrr+6.cv)

e>x61+J | 6yx6ry1
_(syi + c*j+ 6yi)
l6v u*l I
6r+
- lfir 6/
T
L3fyo
I q^ 6yy
r<
tyzx

f..,
,tl

f,'t- t*

= :L(t,-. 6ii - t)
"tt)

q^d"-q- k \Ao.*8n4
F.^^dc"rr

ilAe3+L

Co,"

is vxyo',;rxr,,'h!-lnV
ynldco'"&t;a"r
Eeco.^,r*'tt^*
fr"^il t, t ;^g^a"^rtl,vsLv\)-

( i-e /tJd"'sl+^1-"sh+44o , ),

ln +un\AoJtI
O{.ltwJ*-C..,*rfiun
f
kl

5,, St.,I

Lz+ht-

it

iJ

+Lw lnval,i'a.,nfu I
J

/q)tr^tTnz 5-17+ts,!

"a-!"u 't' fuiJ

ffincifo'l-

7hL etfi-Tus'a+tn
a/"t+- the- -Sal^lrrna, o
,l:

i.,c t^1

I,t"*jl

J,x/

'/,

s)Y'*A

Jax

s?

N-_Jrt+f.ltJ;:o
Jr,

wfu-u-

t'

m u-rt cAT vw,v;e/^rf tp

J-,

tt
T3 o"'F-thuon,.-amD

Jx *

: p

-! v*
t -

,^
.!*+-tr

s+re/4 { SJJ
o{-e-r6'tw{L

): o
(1- s,)[tr-s') (tr-L'

5
15- ' J 3= o
S
,
t
S
,
I
3
+
J
r
l
,
)
t
r
P-(Srtt+Sr)t+S,

-"y""Y

%(4)

o^otEt(5), :'rj3et

(6)

(. 3:n*l;:i1-i*,.,+1,'+r'r.tzr,r

s3
IJ, - Jr-ts-t
r /^1, rlr rL\
/+ a*n c.errvt
r* Ji'f
=::=
(srJ,)
s.s-r
t
tr-t
:y:,,1*
iL:-F,
ta,t'J+
h3*- Ho{t
":l
:- sr !. 5; sii,tsi,) ( -f,+st*,j:\t#'H"i'**
J3
i(si*
t
,, J" - s slSr=-J -frJ.(J, -"-,
\
-{
/)
tX(6) , !r,'r.
3er
corfarig LXrH o"Yrd
(ur.^Io-r-irrtkr :o
r^its;+\i*
I o=[t^r+Jyy+jsr;"=
\

fJr

+lJxxi"r+rJry!**r-rrrr

* S"i I
!*o* 5**5^')+( 5f- t Sui
t S.1r
f. --(JxxSyy

tJ
= +(sk+ Ivi + Sti) '' ( tYk{-s& 'o I

1 \
.

,.r, ,1.
\t
rYYJYL\
fn
I
I .J+, Sy 5l* I

b+ ^ vt 'tov \uarrtd
Y\orm
",

l
I

/
I

l/ttwL

obla,,"i'A'd
S+eA rnUnUa,'nf,S6i{) I'rn'rre-

f-tVU- 0/rg thg

--36,

Ir

Jr:c,

tl<t-

lo

(^'

F*d"*'s^rt*-Q. ftvre*tt'n't
,T'') ,
Tr=-i(oi'61*c;,-i

Ij=-duf(E)

fL=|SJJj

J3.:d'e'ttf';)

,,

ca'^ bL cuzi#t,rn.aa
yigtd c,rmdr'-fio-"<-

?.lth'*rr'

{ct,,r,,T3):o
6{

t'

+ ( 6 , f " . , J r) : o
it w,'t M

*f*4,
'4t
Jo tfti- Tleld

''ogeru^!

(a'nd pla*2fut"^^t.*

ltellj-'U**
u'*ldt^

u'*Vu

41""t o"

+ha'f 6 hat ntttryfe

vnl'il^tu)'
6nt2'r'\'et
"'e^\h

tzt'pAsamco'n bz un;lt<'r a't

fcL, rs):o

(l'-7\^^+afu\
con&'*-rw
'1ie(J
u5"4
- thc moSfwida\
I^
tt<r- l.zyead^n
{cf,)--

L*

I,

i5 alto

Jr"\y*l'
-l/o^

L-t-o

M;xt' c'ti}+';o""

-th!, qrc^rteil+s fu*t ttu/A '

) k, ;l
l// al ?u,ru(ii:n^dn*#xn
,"nJ,!,(,,'t
- 63,
- '1, i' t'-fnr7^th*-r,,
TrV6,V6) , 1"h-{'*6t
tl

ot'urtlilnl,

r^ WA

Q
't'k;lfrt;-s'
;'-v*;l[rs'-;;'-vt|]:
, (s,-'ii,)'-

:
= +J-;'-tT Jr'-3tft;'T"+ 16W]' 6+1 o
-) (T,,Js)
+

3? Mohv/sattis

rn'v7t-'-^1'Y
14u'ull't1,,'lr*Irrn-^^h'$'
Sfi+*t n'{tt--L'*"^* J.,'o"I4
,rt)w
h
a-Af
+ f"h"pt{
Mdle d^'ffi
a'ndTreacp'
6*hwp^ vrn /v13u
fuctiu;$un*
L^- pr^c,h'cL ,
?mqT
-(tlrrr; ts"{, utilhll tJl';c^t uyvo*J
Lp*l
is
c,/rfry"itu
-f t\t_ s,w,7l*
(',-F*

l-rr.r^ ,+

vm M\ tu.t c*i{e*'tw'. 7t5l<0arf

'
&x p't'i.ivr'-r'vfo{ "&"-i' rp'}asvl'l

tJ'fii' ftt'iz'\

mui-f

wr'
,Do7''"1

i?-R* l,-e,^r+,"r1M&2 )
C.r:,k;,;s'^n,

ME3+L

V,Jo**

+?.!*'
y{,{'J:e.
tfu pnwmtttttv k t' +W 1'Ylds*az df h +tu"n|c6'tY::::r^f,o* o

fX^ 2 o,

+trrtvl,!- ktvf ,

L"

lt

tl^,^'

F^^d xvw\a'^fu]- 4{L^h]ilt':l

6': * a,
Jxx :;

)_

syy ='' --j- 6^r,

rxxt

&?:-i*'
a

)= i +
Jz= i(si" + ryits*l,k ':yie}pwt

-z
l
--L
6,r, = T ov

fz= * ui : k:'
(V,'n nfitu )

C!-1

K- ;T1i=

6 y : J 3k ,

(d gnaJ^'-*rt'a,kr : 3;
( fu Tvusc.,g'rt:
s'lrttA *7 G!,(6thuL si'lzrt
s^r'
i: :r,.,f'r-ted 6 ,lrn* sh-e.o'"
tlt!Y'
7lx
ry
6- y

th-e/L f z:

4t onn* r{

VNtd,

5^/=T"

I' : TJ =' k'

sM'$ ^rtr"'L ( V,r',Miel)

,.- T. : k- , i.e L riJths- cnlrc*!-.9h-p-ort
A lrtittL ShuarrS*etYt iJ *Yb"eJ

, l<r: T, ols' )
C Fto','Tru.* fl.\{ cn'r$hrsrr
-->

6xY

<-l[][-,-*.

^rg'r^&nA/R'
J

_L_ 6it "1" rTxT

c^:nd^*i6r4
Um\J

V"n Miliw
;

5nl

th{- truug 1.,li.n-}*/ln s&v:l c

a^t
W*.u^- 6v a^d ()i:r:,7,

II

6; + o,li,= k'
Cw-uL *Y*fan*]
l2xtn'o&t=_ &ari tlre-

k=F: ae]ffi+

( t" tl^a- fvga6e, W(d

N^vr-lYt

k, =t''o/

Tivxo'\
5xx

q--i51.
dy=.kf

fxx - gx/

1tl*'r*

cszt$-$7tn 611 +/tt'

*infr' M,'fbL

q-^firy
Fnrla,wt*u17"

MF3+L
Tn ginQ

L;l
f ud

lz

lD,4

br) EXpero,nt^e.,,t'5
On'fe,,,,r-nyr...

t, botl, PMirn f"r*

Co,t*idp'- a +^lx- 9ulje"r"t',eA

tr

a r " df,u I * T
wil( h- s"t\u^redh brtlL
Thr.ulrl"er^tSlwwLa,bw^(L
?tvr-LaQ st rtn

Ex

,und S/rza,+' Sfu'v1 {y7.

a* oo**+,,^zf
T toqlrtlaL,Aht@f*y
+

f^urut";xgF*d

u^ b, ,'hupsa+." cfy-aa," d Qt
a,UnTNS
W

witl,

: :-s>yvft

MeaUnny4,w";hZ{,,rr,-truns
?*
{ 6<x,A* F WttwwS
'l

u i--

a/,%u6 ,,^l lp m^f tts j/htd Swr,(-r-e-o^rl u*Wv t-th flr$-phb!drh .

t**rr wr+l+%PUit^QlfE
b
+he^ftWrai
* Vrn ltth'sp,sortfu'l"n Tp17-,,*"g+.n
Fto;,
o

.5

P.H;U,Tl',nlttodttnc*hUtL

r/-%{ltr,,rn*

law
Misss

law
Trescd

,4

2xN (r;, Prt>7"

/75 0
P.2L

.3
-2
'l

To/n,

o CoPP',
X Aluminium

pld(,Tr^"akj

o MrldSrcel
.l

,z

a^ol Qu;"Y'Y

'3

,4

,7

'8

.9

r.0

Fro. {. Experirnontal recults of Taylor and Qrriruneyfrom combinod toraion snd

tension tsets, each motsl being work.hardonrxl rgetho samo ststn for all 0eet8.
Tho Miees law is ot*3rr : yt, whilo tho I'rcpce law ig orl4tr - Yr, whoro
o : ten-siloatreaa,r -- ahear atresS, .ll' ,: teneile yiold s11,r988.

5oc,
A L3o,
)23 (l?3t)

tq

t
r
A
iwn d o/r'\ti w*tNl fi t t a*t}rrt,S
t

Ca!4,,h-{rl\rq- (Sirer4{-Jnra^4^)&e[ot#in4 in thr

pta,r+Zreyta

TL*^ ;n tb

i. e. ! r : t ,

o,

J,-t:

{Uo;t')):

L rv^"'^

c.anrla4nr-

j ,: o

tlrD "Ine'd naf rnz'a,n tha*

dlau.twe4

Pl@

pl"*tt. -aon^[ cno s1qin havLda,.'v

^ y$^lr\
, h.attwL cl-,nz,,.&.,.u

5" ,^sdb

rj

C""

ME3+L

a// s+ruv ctYr'p'*ta

slalz cni"'"'*

snj"u-' ln tb s+t'+t tface^' ,1^'e[J

I, - kt

7hz s+w ^rf"ffi

d+Lin +A* TVld tW

wU

) + t\in 1*i+l*i )
Rxro,!,(- 7, : 4$h + sy',-rS'-+*
iyy t Szti* t rEzls tr 5+"J*,t.l 5x/ty: sg'sr': o
ir= f^r 5x"t l,1y
-nM,

k a, c.a,t*tral-,t on+Az stre*

mte

Sj

{- w')itxndtryo
NhnL"$w sfre/ SW M +h!-W'td 5v'n{aw,tl,l*s-n^rfu?r
ttf * -,
ll;a^.y:
p{a*frv Strar'o *; , Soth^t
a ' *
r-4 |

,.tiWyt s-1

t"
/e'rna; f of(n+iv^J

g^f

t,tfl E;i+ zrue(

J,) = tU r1

^)

Q-= Haro h

dPlu>rrv;fie- 8,!i'

-ry
5q2"Jtil?5b' f,.d o'^/rth-s-c'uLnL^+
h"*
+iy/\!^
fu
Tt +r\k
a
'
( *lfr, -lnerJ ;nYnq sto"r'6)
colq&

"?L
-J

l
!,

Lt becan'*

J
''f Cannotbg
dLwep'-k,
,]
-?L )r L',.t*mp>
"-l' t"'
piL$rt +n!6 .J is ,hdahrY'
.

UrSthL

fu,twn''ngJ
-

Lnsteaa,

.PL

.ij

UeL!'v\f

ttt
Sl^rai$ S-tatu (
)'

rnwy(be-cl-p'tw,n-.,,^ld,JACCUtrtwl"AnftSnLJltTlcz|-rltpnfi5.
:;

PIod^'so

tj'vt pta'>R'ry
I^ a,Jhm, i? ln'"tfuf' f-^ai'
*f

'
cj\tur\&L
g;bl.!-val^^rN-

iPiL:o,

pla,+ht. Str-alt
"tzufw'1o-n-c-

)'til-o-

i::

- o

-lc,rg- ^-o pl.a":lr''L}k--'u)'

ME3+L

t+

Gu

Eqwahbn'xr

F*^d*^{^t^[

A,ts**yhrrt
it 't'tzaclsd)
co"1d4't56Y\
,'yVt4
Wftlrrr.

2y e;'

ela'+t.'zS+Ylin
NoR. bd\" tlre ds'aatlni<'
slra,)*
?Ad,+lw plt/>'Fv

r"ftffffi\

il!. ch,e,,hw,f
ri,

,rrr.{"-q- }" iJ a. pos;Nt- scoJ*{ {-,i",/

tF bq- rL{",tg:u'ra,i*rd

A Jr'i
J

,,$L 2y nf := -i)' )

*A*+h;t,

i{: Wfi

,a0

: J,'i,.,u
?,1
'
J
)

@t)

( ,51
t

';f*lL
4-

( J,-t')

s+e,et
the- de4^u,'+lYui:c-

+ ls,r
z/ery: ,u(dff" e,f')=: .s,i
T

dl,l*an;nw
,

6l =

inttpduco'

Na?k l-ovti--ax)ocr^fteA^':'
<- Y*tt, ':f
''of

'

S,, rJ'
.

.eI

^/Q' n'a2/ to

X,

d/-2- f-"tt-t

'Dattf

^l';'l-

It
1ot%

shy+ eV

stnzt' elashL %u^3y

^4 t"k''f
p^fi ,f *AnrJ' C^tt;y",ted

- - - I i i f e t i + i'i,\)
,J '
J

(w,r.l,.J'r','!-a''vnaz-tqdwftLwLww c'Wry+ A

5+audoA Ca*z'?'^uw )
ov'ctittL'k
'
'
'-) ,U
.),
:
,
;
t
;
r - a J2**x +)ixy*,
,.i
Sf*ey+
+
d** +
W =Jxx (xx *Syl 4y t Sz+
'

5; (S,,j+),rl)

W:

2/"

-T

kecaft

Jz-2 r

2/N

tc; j - 'C) '^ ,

?

Plw*z r'tr'ne'xa"l
Fr"le^H"U
^'-Sii- zT -'l''-:- zT *
:o

f*

)")

T:ft,;,/
s.;
s,.f'" ' l z{r.{ J ) ,

fi; (
r',-l!.tlo"; t, oklo-*'"n-**-+*.**

J',)--,lt i( - 'y(e",-a!,!)

s*ut
+)ra,h
rr^t\tnstan
f$' ; s''t' >*;rf, rxw
*x:,
o=3k6
a--3K8,
/t+Lr--fu L,,

plt3*7,

Furdo''rraatfxf

t5

G*t

h,,7'
vr,x*ilsrui

clr{rQ ) Jkah, rch i

, t*nd*,,,a. ryec,.t';ud

Sw^n',*ut. Llta 1MAS

i'n t/r-a-plw'u*L teJri,vrtt!

GtW,r,_
'
q

f',J-

a*u"r'nt 9ru4<l

r* d

*-->
J- 9..(-f (

:
7 7\-

i ,
'''J
-.--=-*

-l-..

Oi,' ,

J^
r
|
-f;nd
6;i'J
'

(J':=:=
JK

5 to

Su-rtttran3

p rrien

----__,\

6,j

/) f
/',
( ifn^Aifte- ^ rn-al4-Lf-ai,nntJ

-"' J) : rJ Vj-yf tl
'

l n

, ' t

, /

a'.

"f

Earttn R"ryq =

Att

Ffl"n'"L^'{

t l^ad^m+.rhtpcv/t t

ff:3K
SiiJ : zf

e;i

I cl-il,r.ntuFur,tc?^n/= )
\

)ied

c*trrd^f-s '

V"n Miq'a )
J, : I s,; s) -- kl'- [
oy - Js k-

Ptq*'c Hft*.

ef :

i t

W
; ?t
c-.
eU
Ztr -')

( rlp'elds*,tzf i'n *'nam)

e,\d,/'4'r
sfa,,}^t*fu )
llo{:r,r,{a*mrr( p te*h-c S+rci^,' rc^rtc )

'
: u^(/ e,l- if;"::;5
si
)
G :3K

( d-uriu,.*o,r.ic ylr+44 rdx. )

(l'\^4d*erldpt{L'
S+re'' ta'fu)

; , '=
g
) S)+66!

fuftr"fl fbfrwdarw Cu

^/1EWZ

I,n +Arb k t'|wt, ) t^)o gzW m,tnt dnsu''W-WS Ftr'r +lrt:" t7;eId icoryh'frvw|.

( t/rn l\4|x/ wrtd Tvwc*) irC,ud,ty f 6-Qhgr"yt":'n{ reTrw**un

,
vw to**icr|-,
arnd U pe r'ft<-ot\t+t^l,
f;

4t Uratd
!;4't^*Yy:!*

li:::!{PW

/
If -t/L!-m^teno-Q.rJ riolryic , tlu_ thyeepricryJ r*e;e uafutsa

-tlAL in.fvr,ma.hst'L oefi^t

rtJi

<:antzttM

ntuttltt-

th-L rna*-en;t*L

Aaa rgaJ-e.d tl\e uyield6md;4*rvt {L rwf

waA4,
Tn V+,14trL

th'y

cn1di+",fruiS (a,d,ret trrL ^ zD SrY*,

W1q.td

/\ tfu-,feil.lrf*,

in du-3Dsyaesyto.nnnd
*,V6 q, q
)

+(6, , 6.. ft S :

()

lfuA
rur,++** vnrrTt OYlu.;1<+ntsf,,^ue-sbo,n4tha* lvydnS+^,tl-

no ,,rl*' irt ye ld

6: t@t*'6.* {j) [rlf

6p:6"=d3

-*

-(,^r.,V^
lJ^J^!-

dl

cf;t6i63:o

{q,6,6t )=,o
ai-.

r , 4r? tJSr^^tar&
/

L,

, wittL t/v- a,xet

l>rtgrn$c s ha-1,Y
1ut

-(u-,ts "te-fin J b-t 61=6r:-6;

o-tl^!Nwwtd,A, thL lrt*+'r[*'Ai:vt

( l, ?- tttl d,ie o*")

[

bafwe**

+hL
f Cr,,du, d1): Q n'od"f(l'J s"f"*
Cuii"tI b,Y {rt6v*61 : chrt y:L'v,ruhu*t .\ c*rP I drs Sumuslt"ye-

f'vlls$^fre

*Y

throf tl?-e-Ortd n'r'/rrz

Th,b lrtewa^*!

- 6p6y=67
ur^L

Lrnaa.+n* lr;l *uk

C".nrr," dn +h{

t/4e^ t',,t'L %,4il4Vie ll

^

6ttA,tdr : o

suvJiro;;-h" L4 ti ke-

Q,'-'''a:f' .

I^ ofks,'\-uvr=e{ il,Ll u,e rr*d ra -h /'i /6ft6r+

63= o 7lc',rre'

S1eCI;f,/ft/,t Shaiz+'f tLA c.{.t'ue,fr* yhvtd S,nlzr,-rI"L,t t u

rt-bl;i..u**r"fu.J
hV heTalr-h,ng
+tDi Chtj,f-aI/^1[t1J.

h**fuc-L PA,4ffii

TAN+L

Co;

- i

// i:
,
i

Movfe*rr'{. ;soffif"J reqwil\'- tlrL y'elJ s,^,y'oz* be

!
i

Syrrtmu+*' hJi+h*ff^f
, t D

M
khli

1-

ro{v,*wvt a"rawd +ha-0 |t : j u : d r

"f,

i
I

f)-, S, ,
t E
, s-g, 6 -61 , 61J

"r'fu*tl,*

{aAitu)<JfitJciv4 be inpn-d 6n +h!- yfitd t**

tl,LTetJ wW

) U I L'lA \
l
i

6+ trr, 6z- 4,
urfrilau5 mtvvvt

. - . - . / / ^ i

sy":fb,

iJ w""tfi, nuMtuw-Q,d
t, he c$nvrN, hrlsJn wL

r^rr> 1qd%.
r,^l
i I{ dr'S

9z Vrn Mtre,.t Yietlh'd^'li.,.

J , = + ( s l + s J + s ;=):L: '
Mofs- 6:

ds) ,

|Gtt6,+

: 6t -T ,
.$,r

J'z: 6t- {

, Sg= fS-F

(c,-6)J, : i &"-T ft {6,-6i'-t

_ * G'*oJ+6)'- te,1
=
-

tt6,-t6J+6t*[8,-

+ ( 6v6r)- t ( fr- ci',,,)'

Vn lrl^V-;

ePldw
6l:6-=
$/Le-

- (rr,6Ltc-{r + 6t6)l

q,,J

sr-$^t rJ
TfroVrnl/t'7'u V4Vtd
^ r$"u&tt, witl,t A C4Tu"r!-s,/f

fs

urvl' t0-w1.

Th-L Frr(

X- y)rnn ,h tht- IruxffW.

TWn\*z&W',Cq.,4'
il'ay
Wt^rYag,

NIEH-L
6,4r=(5

[trt1

Iotol

To be mo{Lprcdu,, de sh^ll
417P4,
J in Wnuot' v,PnT6 6t't,

6L

[oor-1

the-

'IlL-?lo'nQ,',i' L'
S,, S,, S,,

|[.ii]

tr+6i,f6r=o

s=!(r -r-t).&Yryl
'g'

| --_-

itlzTJ

LtotJ

ltr r 11

* s .

frziil
5l

N'o&,. Srtjrt-t):o

iz--i(-rL -l)-(?)==FI
\cJ-:)
,
sr: *Cr

lVote +hosf S,,, Sz, 9j ana npl- ,"}rdl?erilr^ti

I
o/re
Sr,
5',
not
S3
;
Ttt&orYrb4
h Qn'l^ "t'hlPr^t*U.,rf.-o"r
i
i
q
a&oL c*x*) ,) 7er7a"d^'c"{ei
in$^dnnc-o
of
S.
+ +{,.{-orYb
q-q
' (, )
r r,

s ,- 5 ( r D - r) " \ f r ) = j '

tligco^ro-e'

Vnn f{LiSg5 yVl'J crrd"'*ibrrt

t{ (, 0 - )t + i r , e , ) l = t ( L 1 l

15 o\ c'Nd* of r adnvtA

t t ( 1o r ) + i ( t2 , ) | = l ( . r 2 )

^?rc: trt

LUe hd'+;e

ln tho- pla,na of 5,, E l-

L\

t ?-

t ) - , /- 1 - - R

tvttS+z

j \

r
f
n
,(.t,
n^
tttmltl,r'co!.
lteyvwr**W_ !0,4'-_
'

'

i'^ rlfunnrfu'vt
Vvn[Arv,sC,,ri'lt"n{n
6,

{t:o

/6
,/ _k
,/

.,---,
/.

-o(

6 , ' - 6 6 r + a , 7 : 6 k - J - - -q

rJ fuAn(a cryt:hdtu
Theinlo,rt"tat{6n

ulrhfr=oflw

rJ a/^el.lyW

Witb

C^d+r^firnin thl,
P|,*"

po{htdd-n-

SpjrYn*frq*N

fletdsv't,*

's
\/-

<5t= 6{

6v:

rf,:

6r-*dY

tfr :dy

6r-=dy

{= k

G= -lt

,f s,5.
{,y dx"frr.nfiiln
, Srt !-,--F$ -_o
5-:: - i I'tI'_)
!.:

.(sl-t IJ+ tr")

--:(tl* Sl+(5,+i")'J

rn (,-!. yla,r*
r."-fa,.u
tytfrJ4, 7ht-t,yrei.i

-R

4-',

5i t '!'r' -t J', S"

,9 6r^pt$*-A.

l_Cei

Trv>co / re(J u^J,r'nvn.

9I

5z

9g=r., I^
l

Urn$
T'U
r

x
l _15
r

g ,*,4 J2

;.1==rE

3 5.
L

z"
l r
- 5u
t ==2fu.
5- :i,
+
+
,
I
l i
J-' { l ' L ko

l
l

t_

S
) r' -

fJ-

TUCt

l's

"tidd

annda?inn-

Wr6."w r7 t;d4-W

*u : i,"

ir1 t l'I-a-yta,rtz tj

S,, J., 9t

2,-

T'\
\

*
F= o'r..7
J3

),

) q

o-6c7 Fv

4f'eld cozrrhhM 15 \

Truco

twr;,,r* irtson kzd in tAL

c-<Drelfwhrf h #L

oircb

Von Yl'W

arnd,.'#iu-,T,

sr
l(re&-a,

/un Mrrs4

irt glo'nZ sjruA

Tyurn ,ne1 cov,dn'L1:w
(dt:o) = 16,-4rl--6'r, It,l -- dy, lryl= 6,
i(
'r

/'

ltpxq

a pol*t A otq- in Sc-'k

cl
J

o!

rh +he-

ul

eflLlt SZ car\nlt Pan'.l"vY t"

a
I
f-'
drw rv{Dezlctr d"v+rra.

1l<L

'

t?epre;utalwt
G,a'
t4e74z
6wlu'col,
\-/ :
'
-.--'--'-'---L
t
I
tL'Js
9+ FLrro
Plarvh^c
l/I^f?il*{, ( no had*nry)
:

,'h"Sh', rzafon/W.

f:;

"-L c..

;r-/L "'J
/
n

d-etlfalwn"C,
OLa/ytiL Strcr;

clg,uMlv,6u

s.h+zvs

an#
^"

t a

:fi
/

Itrce-rerltz-Lag-ula'ffi, darfrL
I'u5Pfn42..

: rtn

i')'*:v

ptatfL sli)a
'fu^lwwhtq

dr"
trrorvn*ala{

dna'aitu;r:,
s*-*tu)

,t;lra.*

,1o2^hz

as;y

L1.

rar!'

de*r*f triL)

,(e6!: 1{
'

inorrnlnfu

Jr,'-(A"A't
J

y'autt- 5.lYu>,
//
t

ii

(;',W;Iw4)

./r

iVwasfu

ol

nd
&)==
def + dti'

'\

,14

' LI t pr"t(el wuw5

dsii

'.vn'a
o
bttu6z
shvtt

%: *ti

k/"* Vtetd

olgj'lYu+o"it'
Sfrai4 fffi

bf,

5, learthr,a
y;a{a(srehy, &J,b p*rlloL U tJ, /A1/ a:5

4t* Wdl cor{thwd /(.tzl}ul, u *ry.Wu"ljila S'

,0,i
t, dSj
dq u nnL*fu"yo,.^Xlz{
==d*f + aeli , dtii,t,'/
ds., dl,i!tl sJ
4orr'

Mr3w

tZ,trvur*atan
6^eWrrL
u
l

Gi

,T*,

L/, nh!, F,[arr+v

tii' 4J /
T*
{nl d.Si'
0t;vwS{

otr*l

dht;

( ry

s+r+'s
clptw'wfn'a)

P^" Ul-[i t" S5

^f pothfP

cnut'wr tt^!,pM+

ryre{t**

graN ou Prt',^tl'JtW

at Vsivt A
pc'rtc)w

/,iu.Jry a^o[Sn^J";rF6fo/''
dmn,PnmtR bvh^}et^
F"',d+tu)rr,{tr}t

= at?} ) ? ?': /s) = 2/'aR

= o,{i-., Q,R
R"4

&,alrl.,,fO$v
5i-ftrr4

y atJ",

-[rnnn--

-l<nwsvv
-L -Cai
tL ::h,LorL

trMTUL

/i

----.-t---..-

--.--

p"**ls,t 4n "tf^'J* il,w

Vls,w^) {^ft?k0-"i-^ r;uytn utrou*yt*
matUrn!. rJ s,fir*eA h b,VL te'*v,^nuo,ndslqn"ttlc^r|,^r.tLCavrLbe aJr;Apd

-,fA^Z
ry^pb
'l

irt ,n tu"{& zvnd.Lru

l<zrul-raattd +aZnn1.

;l.l^^thaWl,u^t h sodvuA-ptat4"^lyF*(Un,

tus.'d^lM (Ar-ln!%:ily_M
fh)

nwl,u*- ptt'*b"ry-,
t'
Pc^ul-Pa8h"/A Lo't*-""e-vtqL!'D

%Nryb ;S Affnc*q& frorr," Dr

o
h Thpss*'ul ^'^d Appkud Pb*+-,2"'%
"To/*dr^^Av?,
<r"fogbhd*-e"'t:qJtutn_
W *ftr"J ?hAD/.D.,Stnrlo,rni-(,h
le^
9l' --['rob

-It

lr-lf

*--ll
il

' rY
[lll

sl^trvrr"P^tr

rw w&k l,*'*V
rt14,{^rrr'9

- ')
,'
v [ ) Y .lJ
t e vuil-e'rrvY't.
yun
Mrw
lvll>{/)
Vnn
W

K* * , e #
o

t nwrrtprwr:bt-r-rnoqfarV{ ( / = 0,5),
l

<-+_-

0\)wf tLrs
Ln,:'owytt>>hhhO
uil,^*i sira,'rteff odn\u^Je-iwl'wml'
(dwra}uzru
sottL^t e,f : ,,r{
)
-l-l+^)),ryhp, o{,,,-s
cuut5,t-ny'.
,5,i,,tttv

trr' =- ,;'

( d.uriov+rh\)

Net
haL
4:e l*+= ht= 1.41
@r:q)

h,lu tlLa,L!-arlwdxrL frry drft1:{

O nA-AD

*"**

W#,fr W"

fd'-ul':t,-

a,Y^\tVthLy'ettpo'ht

O oB- BD FYyl{^"d in shtu 6 4.){q*Wfh{-lie(Jpnvf

tlrun (p^d i"rr lpnt;rn

M["l+Z

k Sheo"{
TqAn^\n

nl0t0 . rrUoha,rre,
^J';,r-rru.Jh*t^?reseilfit^?'U(V -ct-S)
SP l=LlL(l*v)-1i,^

11,egooJ-i-l fr, fi\r{ fihL s:{{ws(o,a" frn,fun + Cr"in)

il^!,14.trJJD

q
9Z

Stro"-^' pa*tV

Shen{

Te'rt*"r't

St\'v'>"Preb @=

!-

e{a,th'c .
pa*h c)A, M,( sfrninsan'Q6 , r x : E . t r x " ( " f o %f r t * o h JT k-:-6r )
6yy:_tr*?,_ e/ = 6y7,:6*N:o

C6ru^rtl4, W

tyy: f**=-DErr(:-t\*r
^ l

AP, d,<Lshvur?wd-t, stj

t'l{,vtwp pffi
0

T
JL-

r \ f.

j,.: o

b'

??L
cti - W

7kt

(d *,j*
rn ,7;e

0i' , ra,fu, *edt?lL zas,r;,'frd

Yv=J)
t

5,j

with,shoye
ery

(In,+lu,o
vXM*, /( ,a ^k, fo yfu
"f
wu,k- d,E71L'.
bgcar^rWvu{*rrw ;@)
[vv:
yy:

!rw*

672-0,,

F:'

c*,

:f

L K
/ ^
unr, Syy=
o^r,
- lSxx=T
56n*, St*=-t*
7

n sni
+ (+i+Jex
)
J,: +(tnir:ri+:#)

= i o^i+ o^i
_:=S

CYv

try

LE

o^/

i
l

l?l"a-rn tI ihtnnr

l!E:t+2

r i

i
I

6x-,

--l___

,,TA
6x7,
-tr:

l- (sir)'
\ k)

i:i
Tl,rLSoh#ion"hfh// equwhvn

2l

t{

frrvturut'(
%l

Tr'-

(r,Y)
W- tu^r
thot

(#

l(

II
I

t-ful,(z;l-

t 1
|
.-v
*'L.
:
6xx
fx,,
+
k"
'
v
+
""
"l
)

I -xL

= (or*xl-6;xtry)
(oat,x )'

%= Jt|{l*
6iLft)

ardanlypll

Gr*

--

(c{il,x)"

^tr
vv''v
1
w|L(341xrtt
\--;r-l

7/*r
k

t\

A
At p;nt *,

6o: i5&- ,Mol.,^r: '

slrun 6r, de-cttc*sts

nlitl,"inuWtrngs|p*, J;t',a'in,#q rtvt*ua!-

Arv*rfffirre,ffi
ktz1rJr"r'\ttu,us^)a>\fr
fi,n01"V, ry
ry th* stra.;n
5\y + c ,,,

6xy -_o'L

( %'nd un&'li'vnltpw+sAua,)i

TrWnrt,g, 9,\'aa'tu

VV:W

Shew.r-Ia"sionn

6 =', S+r."i" ?afrlbtO,

ohl.ufu,

, G;

a4L ca'tt-,^v
ofi, 0,1{s'fra)'^-r
?arth

Mf

(E- firm o b k-)

rry = 2/,' t^/

8x : 6/y:

6]? ,= tTyL:

62, =: o

ex:- 4y : 2r* ==tp - 2tx :Q

h tf"J o7L7utl s,nfv
$k,fl- p^,ahBZ ilL sh'vv),,-e:-oL
I
,^'r- f ai'vL : rc
6
'
^
Jz=' 3
I

i,?L
-) - llt

n
i(xx
P L-zE.

5j,

zh-

w
ily|'ryihr.nnpr#btrq
J

5xt

: N - .J1 n
-lx:
zR-

;*

,.
L x x -:
^

E
L

(t_ I,- \ i
(-w)151

- ^ l

+; - hJi*
T
r

vxx

,r.-tL
}<x

: - - ',:lb
.
)

/t*

'

'r_,x^
;# =

&

.l<-

r-ffi)'

\-

\jrul
frrrt ) _

W*x

)
, 1,
q>
--!L
-r:
*^ &* r t,::
E ' 3k'

w r

ht,(syW)

Ca"
So"*U ga.U"8D,
6xx

: J7ta*"t(.n',,^
Tl

Cx/

cast"
(Gt%l

ornpoud utrcA,
7a*U +D

(* *:9lw.
(*rt*c^-)t, ,

{Az TtU

gr!1,
pt.rb gla'>+icSrl'nr,Dn
*

a^d y
"f &x

an-ere,r*rrted

e^,t fwnc*nnA
"{ &^' )

c^rwifutw a {1ra.^>'y^*'n ,r'!^tJ*f^* 1'% {r*, o t' E

t,73
t

J3

t.JI I

o^v/<

o.641

Et*
JTK4t' p*w b,

- .Er,n
/ - f-xx

o,rl 6^/: k

fxx:o,

l/)rtl" thr'qg.o,/v1Y ft,vt,nta!-s"lrat rt , fu

/r

1/taan-s*,vvl dy

-iltL nv),n^"lrtslrt*s )rtqzavz,

kf,;b

tr f"M
sft"r;,
ry rm,rna,tr
'

Pazpt m l,tult*iny
I

6xx- jtb

oLoav-eraA
-

X J.{+t u,f
{--cn tfu

hlVndP, ar/,Aury

( Wn.,d
:ihw in 1&,t'rirn
)

CV n O .

MEz+L
lAlo

I
Twr<a'm'
Ca.t
Plu,a,Sh'oir,n

no{/) u)v,k thvw,tql'v

a"n lrxryyb.Wk'un th,nl" rc(n;{s4

ilWrnr"q,ri-,r'
Co-[ wutA^I,'t^

t l(u,sh,adB.
lwrNru**r'Ct-L

wtdlwdrl

LaA4

app}ulh ttla*n;tu,
yrcfuu

9|

| roblBrwS+"hYw!^I
6yY=o

{'lo"r.cstn-in cand"-filsrr. C * z : o
I1.z mai-ezrial iJ {,oailed t u Y t
)l S,nly,a

i.t {Ya-uti\h f-tt".

ASsq/wLay6r\ lWs*t: yvld

r(a hmdvn;n1

c.a,ile*{64,

VI o.5;
=
tlu^ thtry V o.7i, l.e. rhcnnyrtr*'bb,

(ovn ba. solvvJ *nl"l,,;4rhr*\'.

Tfu c-*td*(tn Y< o'5 i: "th{d' mo,k-u
ml-biurc{ necz/vtNV

hnd tfu strW-slvatJ^'+'(/*v+

6*,tal

Ty: o

Mt"3(z
Ix

n/

-f^

.t

Vlaru. \trrr,rivt lgiltn,rn'

Cu

ElaNhc P.a4;*,IL
Br4t{u Vrtl

occs*tA,frht-s:htaS-s4n)nreb4;rn ,lt hha.^,L

- luu - t/,r' - t* orr

f,xr,
q^/=-Ioo* t"4 -|lao"
ttt = - t 6*, -

- !e.ar, * o

iaq

6++= v 6x*

Sxx: 6xx- 6 Sz*:

If :,#

or^
gyy:6rr_F: _

\i-orr,

Y urr.

fz*-G - ?tlL cr,

T z : + 6 ; + r r i + I +)i .= g ( v - - v + t )
,f ?.*(d.
At w\,r4t+

j,:

d<x: g
J r'-rt
fn,
--------:>
l-YF'::Jv4.tl

6y
C:

i-

k..': +
/

*-: L@, -ratr*)= )- vt nv z x - f f il-qL,'E

E

fy

rhl rat^

*:;

rvnalnz

uufth^"t tkt- e/a,"tV

{trt1 E

is h^yuuu tlts'vt 61
( Tfu
\ t1,,"c1tnb!^h# t**:o.

u"u+wnzt' ))

,^-"?i

i
-lerw;m-,
7tow. Sfra,ivr
fu

VE3,+L

P(ou+;, kYm-q-

9l

is h d^no1t tlw 'pnvtapo,l-'

ur,k*urr1-5,f nhn
pnbtx,r,t
J hfeqrla,f tt^r*,
, g^/yw uwt ofun u,nkwoo,n
+/<e
an,d f,'^/ swff'cz't'nt-u^le.r "j eTual-w'l a s'ol'ue
fu

0w' W.fV

' pnuuy"!-' *nh,aaY-vt

'y*,u^p,o!,'
llJe Shall ol"&t'L thre<
unk ,.ourM,
AA,r1,t 6fu!/v ,1"^ k

o-xr, E+,

a,4t be e,,Xyrc*xd irt 14r/',u .J t\r^n-

three

wt tl4!,'t ruer{ t^ ,'/t*4r'fy ',f'Arveqan*'vva k so{vgttu

7rc6{on.

R,ec"l,( fyy : o

6 ': ! Q ^ + 6 w l o z i :

Syy: -&'' ,

Sx:.x:6rt- 6,
J:

i (t*^+6zt)

Szz:62+_-;-

Grl\4^.a-K--f;,

is

brlk:nnod^^d^^/

s3=*-+
r*= edtu* _ ft* t+, 41=,ftn4ry,,
t"*-=*sp(,

=-htrr,
;N

tJ':^s++

l\Aaa

t\l

6t,,oo th*f

-!

Ht^rU, )tra^nr Ie4^"'f-/t-

1',/l

tl fofrL, dtf/d!,^f", w x:etd h ,to/re

d*hA%"

pw4l*w n trntJ'l- /iW- lnc,uemtAh

61,'tt

wL fu,n'J ^""l.nl^rtr,nt AMr;.

rt{:;utrrw
(7a(t),6vtt).
o,rytd,,ve n-aJd,t"

f*A

\fc
0t ttat.
iltr Soih,filrn

6y7(tJat ), 62,(trat),

*c .

[,r. tt+2tt ) , {6 k+ ct ) at P@NvLqw"hha

't" 'trwrtu "f tzx,
lecor^l.o-Ne dsl', h etpw), *"Vt/*t
a^d 4,* 1lurt stTa/in amlr*w' reTw{ttott*=.o,

l,t/<tSfrdll o!'so fte^*

knn'rn :

6xrLt),

unk*w:n'

Tfu:-W

dzt-Lt), t77(t'1,

dxxc*1dt),

o*(tt,L)

r<rt),

frrtt-t

crt I

ftttt+at),

(;^+t^elatod[t,trt|

,
]

,ttrt'n -f*^ t to t-tat ut

,f

AO : 6 Gtat t - 6 <t ),
A5^,o= Srrct+at ) - -rrr (+)

^iw"a- 6 Ct+atl: !;(**rttatt t 6-4rfraf,,]

c ct:t : t[ox;-c+r +esct.r]

ASa*= Ss*ct+af) - 5++Ct2

Tlt'U b"fu

b" o\ ,ry

^ c,-g
A
xk

_ 49 *
R+

a^ .6'd* z- =

aG

+,

Slra*:n
,f <>latvhc
alxx
z1t'
ASea

M83,42-

F{wrtp-Slr&;h Tetvnt

1'/& p(a,+hc.stt-ai'/) )ack-.n*a*t

ea.A

t t" ttz$

frt

Ca., b*

btr
aVltwXrn,,ted
nP
A^ C
x *L
^.?L
ciLzZ :

fraf

Sr.rtt) + :in(t+at)

W'
T*t

Jua(t; + $rC11af,

Arbticrth^t il-e-**f

da.uatu,o s*e*5 l*W.

[t, t*ot\
,7U>;oJ
rJ ,'P/W\Aa'fZd +/La-awafiie \,"1'* ""f t on,l tt ot .
ry
Th^z azutpndz t" the T^pr,*od ^*b 'f rtwmut^!- Tua{arvrwu
r7r fu*/

(orryarud ^AYl'va,r\ ^h<J\'^r+1W ,;6yrt7,7rc'Alsvtbatyl

otlt:

epcLr)
, uft : 4_f r'*c.,

+la Trdf++q'd for rr.'d-polf) ruLa iJ nnt

) rrtw tUhU

<lm@ rytu*U

\f

mfiL acwuefte,, 6t

( h-e.ncct4u**t

^lh

/.atty,,-at ).

i4 epfrf+red rrt t+'v,*aot' U,nn(ttat, fn(+*"t,,$+

l,/,nit tl4^k *"fr-t&Y

Nf- ^,re n^V

fulerL

uru?^dofia- f/te tkruz

e7 anl;arul,
I

I E^ct) + ^ * + ^*,
1t*

c + , + o; t; t ; + ^ t I L

trx cfrot.r:. o

(;^pua,+^;)l

tercf+ot):o

(inyot<as*ai-,

(t+ats+ s,++(ttat)f- h-- o

J
l+t tct+ot, * l;' t r
l--

Tlrr4 u a- s<* {
I(

(yvlt
c "dt-trl^)
\y

mn-hxe^rt

,:a'^ 69 6^^trrrttf/fl.

e.a-a*ana frr

coh,e4 /Zf ilz

'fsolre'

3 dunk-rc',tzu.
cm',rturrJ in MNft^rb

bfJo"t (o'h" qr.' Gr^laLulv'tk;.

Tlt& Ma'(l*b f"lu tus. afrld-4.leJ
infune,.h
e(t-S-plaaa - sfrar.V - wnlw6t'ol-plaaf. m
ol"L ) e7u+ltnA
'fsol '
w ; rt ,o/o'tl4- Slrc,ib - urao4Jo!. m
attz svtvd bt ,^/hY
"",tJ^

N\E-W

o'5

Iarten:n
flotr,eS{a,,,rr,

rw.ryru:rziCc-l-SoLln/-Wr,

Fbtn

tl,\l

Ca4

SnL +1noJ|'

ftr",

V b 0.5

\.,

ptnpar .efl^{-

6*t
6xz

A/g

tfarn*-'+nn'ta

d\rrr^d.tl^!-

Plane Strain Tension

4/16/13 11:34 AM

D:\My Documents\Courses\ME342-Inela...\plane_strain_uniaxial.m

% ME342 Theory and Applications of Inelasticity

% Wei Cai, [email protected]
%
% Example: plane strain uniaxial loading
% Material parameters
mu = 100;
nu = 0.3;
sig_Y = 1;

% shear modulus
% Poisson's ratio
% yield stress

K = 2*mu*(1+nu)/3/(1-2*nu); % bulk modulus

E = 2*mu*(1+nu);
% Young's modulus
k = sig_Y/sqrt(3);
% condition at yield
sig_xx = sig_Y/sqrt(nu^2-nu+1); % see HW 1.1(a)
sig_yy = 0;
sig_zz = nu*sig_xx;
% hydrostatic stress
sig_bar
= (sig_xx + sig_yy + sig_zz)/3;
% deviatoric stress
s_xx = sig_xx - sig_bar;
s_yy = sig_yy - sig_bar;
s_zz = sig_zz - sig_bar;
J2 = (s_xx^2 + s_yy^2 + s_zz^2)/2;
% strain
eps_xx =
eps_yy =
eps_zz =

% J2 should equal k^2 at yield

(elastic at onset of yield)

sig_bar/(3*K) + s_xx/(2*mu);
sig_bar/(3*K) + s_yy/(2*mu);
sig_bar/(3*K) + s_zz/(2*mu);

% strain array (going beyond yield point)

eps_xx_data = eps_xx + [0:120]*1e-4;
% initialize arrays to store results
sig_xx_data = zeros(size(eps_xx_data));
sig_zz_data = zeros(size(eps_xx_data));
lambda_dt_over_2mu_data = zeros(size(eps_xx_data));
s_xx_data = zeros(size(eps_xx_data));
s_yy_data = zeros(size(eps_xx_data));
s_zz_data = zeros(size(eps_xx_data));
J2_data = zeros(size(eps_xx_data));
% the first data
sig_xx_data(1) =
sig_zz_data(1) =
s_xx_data(1)
=
s_yy_data(1)
=
s_zz_data(1)
=

point is the yield point considered above

sig_xx;
sig_zz;
s_xx;
s_yy;
s_zz;

7
1 of 3

Plane Strain Tension

4/16/13 11:34 AM
J2_data(1)

D:\My Documents\Courses\ME342-Inela...\plane_strain_uniaxial.m

= J2;

for i = 2:length(eps_xx_data),
param = [mu, nu, k, sig_xx_data(i-1), sig_zz_data(i-1), ...
eps_xx_data(i-1), eps_xx_data(i)];
%options = optimset('Display','iter','TolFun',1e-10);
options = optimset('Display','off','TolFun',1e-10);
trial = [sig_xx_data(i-1), sig_zz_data(i-1), 0];
sol = fsolve('eqns_plane_strain_uniaxial_plast', trial, options, param);
sig_xx = sol(1); sig_zz = sol(2); lambda_dt_over_2mu = sol(3);
% compute hydrostatic stress and deviatoric stress
sig_bar = (sig_xx + sig_zz)/3;
s_xx = sig_xx - sig_bar;
s_yy = sig_yy - sig_bar;
s_zz = sig_zz - sig_bar;
J2 = (s_xx^2 + s_yy^2 + s_zz^2)/2;
% save results to array
sig_xx_data(i) = sig_xx; sig_zz_data(i) = sig_zz;
lambda_dt_over_2mu_data(i) = lambda_dt_over_2mu;
s_xx_data(i) = s_xx; s_yy_data(i) = s_yy; s_zz_data(i) = s_zz;
J2_data(i) = J2;
end
% plot solution
% plot stress
figure(1);
subplot(3,1,1);
plot(eps_xx_data, sig_xx_data, '.-');
xlabel('\epsilon_{xx}');
ylabel('\sigma_{xx}');
subplot(3,1,2);
plot([0 eps_xx_data], [0 sig_zz_data], '.-');
xlabel('\epsilon_{xx}');
ylabel('\sigma_{zz}');
subplot(3,1,3);
plot([0 eps_xx_data], [nu sig_zz_data./sig_xx_data], '.-');
xlabel('\epsilon_{xx}');
ylabel('\sigma_{zz} / \sigma_{xx}');
% plot deviatoric stress
figure(2);
subplot(3,1,1);
plot(eps_xx_data, s_xx_data, '.-');
xlabel('\epsilon_{xx}');
ylabel('s_{xx}');
subplot(3,1,2);
plot(eps_xx_data, s_yy_data, '.-');
xlabel('\epsilon_{xx}');
ylabel('s_{yy}');

8
2 of 3

4/16/13 11:34 AM

Plane Strain Tension

D:\My Documents\Courses\ME342-Inela...\plane_strain_uniaxial.m

3 of 3

subplot(3,1,3);
plot(eps_xx_data, s_zz_data, '.-');
xlabel('\epsilon_{xx}');
ylabel('s_{zz}');
% plot lambda_dt_over_2mu and J2 (should equal k^2)
figure(3);
subplot(2,1,1);
plot(eps_xx_data, [NaN lambda_dt_over_2mu_data(2:end)], '.-');
xlabel('\epsilon_{xx}');
ylabel('\lambda \Delta t / (2\mu)');
subplot(2,1,2);
plot(eps_xx_data, J2_data, '.-');
xlabel('\epsilon_{xx}');
ylabel('J_2');
ylim([0.9 1.1]*k^2);

4/15/13 12:42 AM

Plane Strain Tension

10

D:\My Documents\Courses\...\eqns_plane_strain_uniaxial_plast.m

1 of 2

% ME342 Theory and Applications of Inelasticity

% Wei Cai, [email protected]
%
% equations to be solved by plane_strain_uniaxial.m
function F = eqns_plane_strain_uniaxial_plast(vars, param)
% F(1): eps_xx_tot (specified)
% F(2): eps_zz_tot = 0
% F(3): yield condition
sig_xx = vars(1);
sig_zz = vars(2);
lambda_dt_over_2mu = vars(3);
sig_yy = 0;
mu = param(1);
nu = param(2);
k = param(3);
sig_xx_0 = param(4);
sig_zz_0 = param(5);
sig_yy_0 = 0;
eps_xx_0 = param(6); % total strain at previous time step
eps_xx
= param(7); % total strain at current time step
eps_zz_0 = 0;
eps_zz
= 0;
K = 2*mu*(1+nu)/3/(1-2*nu); % bulk modulus
% change of stress
dsig_xx = sig_xx - sig_xx_0;
dsig_yy = sig_yy - sig_yy_0;
dsig_zz = sig_zz - sig_zz_0;
dsig_bar = (dsig_xx + dsig_yy + dsig_zz)/3;
ds_xx = dsig_xx - dsig_bar;
ds_yy = dsig_yy - dsig_bar;
ds_zz = dsig_zz - dsig_bar;
% increment of elastic strain
deps_elast_xx = dsig_bar/(3*K) + ds_xx/(2*mu);
deps_elast_yy = dsig_bar/(3*K) + ds_yy/(2*mu);
deps_elast_zz = dsig_bar/(3*K) + ds_zz/(2*mu);
% hydrostaic stress
sig_bar
= (sig_xx + sig_yy + sig_zz)/3;
sig_bar_0 = (sig_xx_0 + sig_yy_0 + sig_zz_0)/3;
% deviatoric stress
s_xx = sig_xx - sig_bar;
s_yy = sig_yy - sig_bar;
s_zz = sig_zz - sig_bar;

4/15/13 12:42 AM

Plane Strain Tension

11

D:\My Documents\Courses\...\eqns_plane_strain_uniaxial_plast.m

2 of 2

s_xx_0 = sig_xx_0 - sig_bar_0;

s_yy_0 = sig_yy_0 - sig_bar_0;
s_zz_0 = sig_zz_0 - sig_bar_0;
% average deviatoric stress
ave_s_xx = (s_xx + s_xx_0)/2;
ave_s_yy = (s_yy + s_yy_0)/2;
ave_s_zz = (s_zz + s_zz_0)/2;
% increment of plastic strain
deps_plast_xx = lambda_dt_over_2mu * ave_s_xx;
deps_plast_yy = lambda_dt_over_2mu * ave_s_yy;
deps_plast_zz = lambda_dt_over_2mu * ave_s_zz;
% construct equation
F = [0 0 0]';
% total strain in xx
F(1) = (eps_xx_0 + deps_elast_xx + deps_plast_xx) - eps_xx;
% total strain in zz
F(2) = (eps_zz_0 + deps_elast_zz + deps_plast_zz) - eps_zz;
% yield condition
J2 = (s_xx^2 + s_yy^2 + s_zz^2)/2;
F(3) = J2 - k^2;

1.18

Plane Strain Tension

12

xx

1.16
1.14

xx

10

12

14

16
3

x 10

0.8

zz

0.6
0.4
0.2
0

0.002

0.004

0.006

0.008
xx

0.01

0.012

0.014

0.016

0.002

0.004

0.006

0.008
xx

0.01

0.012

0.014

0.016

zz / xx

0.5
0.45
0.4
0.35
0.3

0.64

Plane Strain Tension

13

sxx

0.62
0.6
0.58
2

xx

10

12

14

16
3

x 10

0.45

syy

0.5
0.55
0.6
0.65

xx

10

12

14

16
3

x 10

szz

0.05
0.1
0.15
0.2

xx

10

12

14

16
3

x 10

Plane Strain Tension

1.75

x 10

14

t / (2)

1.7

1.65

1.6

1.55

xx

10

12

14

16
3

x 10

0.36
0.35

J2

0.34
0.33
0.32
0.31
2

xx

10

12

14

16
3

x 10

Vt)+L

//\

Be,,1d^''ig

L0^

,L

Cn?xi:f
Tl,^Zis w,vfYltWryb in vufnJtrplr.ll/'iearncj( vhtut"hzru-*sy\A
J

in -th{, )1lLoirnerl-

TLutuJ
li book"'fl*Mothrn.h'^{.
'f'trt

Thid *ryqfb ft-ta[v^f"*

'\ P(o*'ufu',
\t,T {,.81

9 |. hobIetn jjfr*Prn?"'t
I

pbn+-Z,y

0,
r ,
r

27422_?

f f i i

('

do**.2 \J

o
" )

-c

'//r?',ZZt

777/77;

/
-O,

For a. ve,@

-tv\.

(/{VY

Se}*l(n/l

bzar" s,,r,f,etJfedh pw* b*nd.fry

"J ft,*WrT^k-rl
frril +fu :tl;cknotvs

a^/, funfunV-atr,am*untz.
k 6A {,^^*;wt ,f

'
1,4.

vt/uea[tr-* fh.L ylot n, \v.,frsy1

+ ^ffre,.* ,f M 2 My,
wf,,pnt,t\Ay h ^ tLre-s'hndvrl^^e-. ,!* 9a2fof 7''old)
Vt/e n(rl,

%yrrt

4,L" /w*i>

vvJ"^' 'M^ax ., of ^e^*

A rT)aV:l'u'A'vrvL
*rc^
e-Yh+rtbs
t, C{t94-t/n,,-o^irW oytry krtW-

arnd 6^,nnot
lh-v bzarru co(b-yst't d-l'-'M = n1',4nnx
4f,
N ,7+*u&t' * *d^ft1 rr*Ynt{*^t

frra NL*,/rAr.
W

ilt-

n^;ta,t"L
th^t-+1rL
a,vttw,w
Stlyt^'r;fph&Shrl,t,^46,;n
( v: (t.5)
is in'.r,''yveatriH,l-

!
I

A\ewz

Bv^cltg

.-__j
,Cu'_
)
|
:lnv
?,1

n*
.9t. e{ra+;c
,f n'ta';"Q.s'*'t^5'y'w
frr," s+Y'e/\fdr

__Zi

'

ll

r r=
-

Qx:

tFnl^r+ro-taxD

) ailoft?ut6fuiwal'tei.;

ry
P

f hr-Sz- exyrc*t ztt": -t' Un^)

ths ^*t01/nfd\T14tha;t plana,r cgTA

YamoJ>)PLa'4a'r'

ry
fnT

e${rW*w 6n&*<aut (lo+t't^llt,?.8J-)

tl,re o,,ttu1re-

| =-o.0"

frd:

=E'+ nt4r
.t'bl
J:-{rx,.{*v

:--(*)5#l : e
ar,fur*:gb).f:,:'c'd^J
M_

E'ji.f2

,K'

-$q. =

Elz )

My
f

r?

^ir
rno,6,no't'odu0 COIJL,IA

I{b't$rwur,"

1j= :b Ot

d-rr==
\'i, ., 4y:4*:c'
-b i Cq-f)'+ (6 *61* (q-d)'-]:E=

ttu

At 6naxt- Vle{d,
d

tfry !r:
ill.=,,

!W'

hu

='*
fr
ki6ah
ry=rt
f(trJ'=
t \

M rt :

a f 'y
L+
-

frAnw

Be^d^ry

a,;?6MQtl- ySeteVU
+

fhe cnVc"J M
.1

Cs,A'

6y

MLv

I
JY

Fa,

tL?

Ec,

/Y :

rJ=+--

Ao M:Y\,:
ul.ot,
h,{v,i

-{Pq"

Y---(

i.e

al c: e.

arr'tlr .
S+ruysfiet,J A th(,tl;*t:.o,rv
!

?3' P(a'vt'tf)^dn'nj

(u> Mv)

tW f
Htil sl,auea
& r , =- { ,

D:'o'5 ,, thltrt'

6/1-9- Sh^il'tflg

Jinrz 6s

*uY/n

ln +!w ,t:,btvtIw
ryftr,

( lu\>tt c ) .

( ttlt, f .B)

,^b

;\ +he-

+ts bewm i: ,h the

H+rnrs

ffi\nu'*.i:
mffiT;*WJ,ff,r

6y
^
t
:
{lo"

,VW r FA, -cJ cvn

)r +hp pla,,>+rL
J I c, o7

I
t f rr-

EY

j'n thv g{,.ul{,. y.yffr' , [-r, c I

J

o'if J= Ct
F-/ cr"+r)n,-r$f.6

L r

L L

1T:

"Y

:(t' -+e/*r\:*
M- !:-"* !,,J{eb)

: F v.i|it(al"r
Qb)

'tdy)-eL)

t\A?3+L

S(rni,hy
C:J;.

ItVVln*^{^u C=o ,

6-- J:4 -..6r'ztb

6Y

-:I
4,/
F -l'2 -

E:L

A'=ryry.
,
fy.zb

{=

)*t

r l

Mmoy

1'"M> Mc,: I, v

. lrrl^^r: q zo?b
L^.L

ay'+

MY

f\'F1r.[".-v/E,.lTJ

Nl^r*= dy' ZaLb

M.-

LU

+N;
l\4

: = 3
2

::

N\E3+2_
'e+

Bwd,rng

i
i

f--;

{ n \
l-t^,l1

F \ ,

r-/

S;,,y\ s,^1yofteJ
be*rn wy&;;t
r wufwrrt LoqJ
',rfuTe*,tlt
u5,roado,bw-u,n be wd
6 S+14 ^ be*m ir, ^l,h-cl,tl4L
D:amt w*ft,,n
, [ !ffi^noL ntmt*t

( lrogw&Ll,ofuo,
67, Pq,+)

PLI

Sht"n 9 Wtnu*

a|.),\tram

V c ) = ? L ,- V c t ) = - ? L

V'(r)=-f)

I
L

l'+1

V r : r )= - f L
M(1.): C

t\A(-L):Q ,

V ' (x )=V o )

->x

M U 1 = : ' fP ( L L - T |)
^U

sLwa{WtW vnetr urj

t

lI ' \ xl l \
"---t

;f e,A*obwrn ii th e)i,a4.,y12
repfu,

lr

p lI

cxy: -'((x) +(or|)'

-KJt/

'91
\t\
r\

beA;*1 ta,tv,e,,,.f,StuU\ ylrr +. dt x .

'f
e,^Trgbt'rn ,.J ,b elanFzy*
,

Q v =- " #

.. Jr= j ur; , tx| :

)r?{dCa-d,ih-rn

ku

rropr<rc,t,cc-.
ltAl <<l6nri

I-.x t s^i

w3 rh !.r,
).{.r'nttTrtldcxrdrrun E {[L l,iorvrt

All'urrl

jh*

exCA,-;drI

nw1nr.rlim
*r,rt^nv+f,J,** (ift:
MX

ad X,=o_)aLnnoy

+hjn ' h ffa pN'a/vl.Z

levh'o-ia ,(.kr f I, y. I I

M{x;- cr ( *- } ) r t, : t r (L,.-x,)
Wwu'x,Fna,{af{*,a-F^
"+ l

aucs
tlqtr*^1> nsr.,rto{

Be,"d,ng

N\E3+L
65.

Ce;

&nl,oa"J;n'7

k6fu2u

af/n

f^f f "n

M /urtwtw /*m tat ba"k h

>5 il-P. rA;a|qaj-

Nhat

r,- F*rL Le.d;n4

hd4 inozap,a-ed
fu^ o +" Pl,, My

We-rwvo ,k*
Z.

o\ lann

sh'?rt5 t1 t/11-len-

N ho* A +;h,l-cuwarfi^4t- tf f/t beot* ,"/h!*-

* aWor,nbo'&d (M:,,)?
Is tt-a-vwJ^*Q- s'trvv tfu,Y e/6 /
(,il^att /M,

bal;ng
o-xr

fW

( M .tl"t"rrueytLlaum,nl )

t4 < Mr

y/

l'f.s:
--I-=,

ff

\
-l[*q"

-.1-\

c/trtrd,/..t lhthL

6n=

qiJ

tL

(t'r ='

Ert

M,) Mr

bey

__I',
o/' _JJ
\1',rro,.
+__-r* Lrq

) T'

'*
6xx= |
{.tt
r

-gs
:
V
Ysc'

iJ

.1= {lvl>c,
lr5

' = E b l f*''b
J

- : : - : : - - =

ft

Fc.

+
EB,l aL-tut,/1oa.a,)

t:=
[

t'=-i

I
l

Mr:X2

t3udinl

Wtlnadrry

-/i

GM'

t :

M( M,

"
cx(

\d-i

q.= r- L*
a,'nbad.rng, 0.8fug *

t+tv4 t1^! beLv,o6 Y .

"*Y
*

Igr^o.r;r\) wrrc/*o,yed

aJI

Stayn

(t^Lood"ng

, t,
oNtn-rqeJ a/re c\cc0rnrnodated oy

"Rf
Cx-r

M:O

,/.{o(y'AutL

<ry

6xI:

!:-:-YA l -L*
*=1 rt Eri

^"1'

t'txrc

It 5* !U
currrniwrr*

r%id^{a.l.tiw^D

ur*i

fr

Mr

F_T_E:

t9=t'

_ _*

--t
F
u ^- - = -! * M ' Y

I!;IZ

f$ifucl

M= l^r

.' M' ct

rz

M'Q,
oo lr=* = - 6Y't- _E

F',[3 Fx.-r&/(6d,t,

T>

6I:

Po

M\fi+L
Rtco$

Wkw

De,nu'rj

C0,,,1'

I'Ay { M, S Mrnax

( M,.o=:i Mv)

Mr: fvtY
*

6,-ly=a= -drr af*

o (

no re,s,'d,^^J,
silz

fl

: *or

Mr ( M ( M..y,

6rr i,_

; _cr

Q.l

)=o.

llv,re- p{nM-",: F{r"

d,

-= o

'6Y

.I" gemm"{", frt

M\.4
I2

(Wyw \^\^;t )

N+\.Lt\ Mr= V\noX

'_6v (

( M " =L * l

CG'W 6rmt)

lr--., : o*l)=o. = --'5r+

"^ly==o

JI

t^fff'eAftw ounbo.&,'\.

sto".l{

1 .*

sho,Id not )

r\L tfpU

,_,*, occq^ t*ry wn(ordnj

^ be-rd^ly mxrv,e^* h ThL

rulPrlW d"h>++n,

o.l+rL t*Uy-

ocu^rv rytu;^ 7

htL",,!- wi(l fYd

ov;[( '/a'vtLt"^,

v't'st*ort" N;
occwLf"*

M
Mnrox

ln"*,t"i\* ofi*r,f {*Vl.Mv)

fAv
l-

I
!
!
!

ToYS;sYL

'-'...-.-,'-,-J

Cai

ill/e,
{ilil{t+"ry +lw t-wfomn'
td c5.rrdt+"'-rnj
lto'y 4lrr:d r6."e
4 prismn.{,c
tr**f lo,, t^Jsw,pe*tiul- 7l^*"Y*
h fil\Lbmdt)\y
Sr,n*L^,r
lnr,nrrvrl.,,r'n t'liU
oto\ q}'urN)
i'{,rthf{.ah tl'ro?ry,,ipl'''.*6

tb

C^{Y5!ed+n.

*ry
"fru,n'Lo.l*i.-tohp{ir"1

itV{.5h#, t{anJvr,u\th^ 1rutk-rwit""}

q^fu,& * JN glovHfifi.^dtUf shrr^a

1q,'&\ils o^d.+l'"qa
ry^*-,
tl,

firublt* 51air,n4r''\#
C/fT"3 te-+l<wt

.T

D! r^'f'e^l'e I + ffiqps1.lctr^ln+,kJ CzI a,/A .

r
u

A p/r:rrtcvt(-W

Flnot+WP\art^zre4uru(Sh^drlir)a,noltwisF pvuu'iF k^fh

a^+"^-ii*Aof

-l-

th-a-ptae,licrgru;rt h $WNL if T >' Iy ,

tri/e e4pe,,;lr(
hsffu.rlTy B o,fhrull^\d 'rqlona- snutr + it \J )
IMA ,^ho Offel,f

o' rnwlrfruurl.\'r;5l'uri' Trnu.I , #

tf,.q-pLorul-t- ne{is

wttt'[J-

trn\rnd,r t" cirfQ/L ti'\L o^r{tt.'Q,c/Nl/\ r{i*il.r\

-fh.a^lr-T:'-1 rnnx a"ttrCannqts\YP*Ltb"rt e.*L\01:xl'

a. T*"h*^f trf?"M-

{-rrabvTy

vE-Wt

l-uf','ort

p e Flathc

Ca,

TorsisrA (-f . Tf: )

'Tr{rr*"
(nt ntn'tj-&s+o-;{-*
, gr.c-ilnf3+o lAl,h4t/UoV WW*- Nott't
. )

-frrd"
5d#rx'4.

.z
}LZt=.}t

,P.

twel Ytwwif *th

s'lra,;nfie (d
7*x: }y : trt :-o

d,aylaw+ic'eu,\tfie{.ci
(

t t

"'-

- EJt : - - f * y

t"1 * Q

{),f,*: p4*

luv:

L x * : + t{ V+f{ * )

t u * =I +tt'Y)

t). : + (pt +f t,,tt

S+"er^
ft:t{
6Y'J-:t-o
Io_*=
)o\-o

+{,xI
i q. - )'',pc-!
I a y " :) , F t * r l , g t

L$d,-^*l-

ftn^dffJ

Gr: f,y

gi"%c

, { C",U,t sur"d^
+ha,f
$r'n c41-tr/1

Qy=- fx

c,s{\d$rsrr 5xA1 f 6Y+,) t G + , r - O

so4{'rilts
thL equnl0ibutrA
D4 c\r;6o'^p\+ *J

t!,o-:K!{Pd .

QCI,J)
(: %lllo"ni e1i,r.a"+rn
)
tro st-\o'^{rue
-->

b"r^rnd+-1 crnd,idi6v\ t^t oi^fe^r t;t^

+:

u1*

f,*

.t b*

c'ir, f

lrv% r+
r,rritho,,rf
,,rrt e-tr^r++*-* : . o
fnh\'L;h,tr,,

" I

,f t\ru-q/^.
Torol,rsu T D n,ln*':o#

- i.l( *P,.- 9P,)Jxd*,

T: If Cu+y-46o,),A:xd,t1
i
:
A
.

r;, * ({,,1)' Sr=f,

f =iW,^)'+
fi

5r=={,St=o,J,=T.'

-[o

VC3+L
fSio-rr!: t Gr.^^hnrC/XT\A 'E"h.trYL

Can

-- ry (x,+1,-,^,)
Qu,,j)
r=ffi
( { r t _- r y f t Y ^ , t l

f =.)uPKt

,r1

)<'

rtg;d;T
K+= twy;rnaL

kt: {or

- .T j:v c^Ye'^Lorr?row se^*4n

r=.1ffi$=1#l:yp, i
(mor,x::,p
pa

zT

.T
: -Eo

l^/^nT=Tv, (na^:p : - - 6,r.

f5.'/

rd
r

' T r- y -

oJF

f:

Ino'

91.2 F(&fuc{ crorcC^,fiona

)- Q(;t,y

&tt
a--,Fb'

7 , =l p K+
,(r= TTdf

a'+b-

-t
(+:"$
)
ractt^ed'fu
{ q,4nlrt7
a=b)
'lf ab
| -l
* T \ f.L,
tD /
r'dA

T
LNN

(rn*Y :

Wlrrrt,T= 1f ,

pr(wr mrrur\,t
af h,evh-a,

'
'ad ll
t
tb
-_/=-tt
_
/f^'1"
-'ly=o
|lf,,.ll
I by I
L,
i ,_" anfb'
'",1*
Y--b

zA\a'b
zT
T z^'b
nF+*=fW-rcabL
-T1r Lnli
Lt__ 6v
_ R==
...l',.=
(n0,,
T eat
E_,

l4Ew

Tofs;trr-

t e.3 Squa,remw

J-.-ii
i
I

gecth-f]rn*

4 fx,D mtxtt*be solwdntJrtf,Lvi*/

4.

Car-

o,

-f :,t^f

Kr

{a >ztf Ar4 _[ : *o*

Kr
N,t-

o .,

fl
r

'frrN/

{-- toq\
u- L
u

J_

',r
x - (0.
t .tu,
a
fnD et(

r( adw&*-

,f
plot c,s*turrv

cV
=lir-

'lu*il

Iil
_ l K 'Yzl
'Yz
i i 'l ll r

8v1

..-/Y
6y

7T
I'

'

{--o

9-!^

&-t7 k-a3

*n :lt{avrsc^Nr)
f CX,:l) f* wu^l.nr' oLtlp-}.,.c

MrjryL

lbrsio-n

9l Ptwn, To'rqrYr

q>Tv)

Ca,i

**'te'""l"( t^ har*-tnfil tlA-A.SLtt*, slrUvSI-.

plo,'*rz

t tr ytf*ty

k-ff cVbrrMt-a)

q/d'cnrd
farwr es&carrinu'L
il hd,LfirLa"

ln if.o-pfartir
Y*

T: h,, i rl I Vftx,Yl: k .

eXcssil#{
pCX,j; cq,nnot
J* o{fo-lrr,
unrd4,+he-rWoJ-fwndr'sn
!:.r ftYfrN

grti\t'4OYuvJ_
c,$Yw*ryL
SDa?L"JrLte-

rt -

x-+1J'

Q t : > c: , y k>( t l - f

c<r<q

Q(r',J):
ry(r=c)+[(a-r) O(rsc
n%r o'/r-I-: cnffiN,tudrl:
"+

li:- fr'"^
-;'\-*
f*;l

P:

UC
/

Tr -: ( ' iz QI r' . z< T

-'T
, .rr d
,rl
).,1

- {:,P(++T}e[c-r:';
= { }.cr +}[(d--c')
c) :-# (+a'-cel
= r c1 , , [
]ot- * )

c : U-'ot+rft,
1y

-lnrx

c:o atfkr,-=!

kat
x"L.fiu|nn-A-

lrrtL
L-J"-'q
U

[EH,

i
I
it / \ I
. I
n

Tbrs'vrt

it-wl

?-L
i
I

Ploot\, elaumu l; <T<1,,'x

t

ry,* Tei,F&

pl'cv
f' t/c,afir

,Y

(-

Tr*

InnnX

(A*a"f etaarryh^/

"

,N)
(fV,*"*4 9oh,',tiln

9l :

t&/s^n*
St,n*e C,x\rA

k+t

Skp'l.J- ULwtkL L Vla,rfz- Yefr|9

a^wwa* Prtil&,no)

I
I
I

I
I

Ty( T (Trn^x
ro,ry,treI)
(lr*nt;g-{ solvhrn

t
l

Tr^*
Ty

- j

_
-

lfqflI

4/26/13 1:42 PM

D:\My Documents\Courses\ME342-Inelas...\torsion_plast_example.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2D Cylinder Example: Torsion with plasticity using triangular elements
%
% Wei Cai [email protected]
% William Kuykendall [email protected]
%
% adapted from FeCalc.m written by Peter Pinsky [email protected]
% universal meshing from Michael Hunsweck
%
% First Adapted 01/03/2013
%
% Last Modified 04/22/2013
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Part 1: Set simulation parameters
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%clear;
global nNodes nElements Coord u IEN K F Nxe Nye EBC ID g penalty_factor
nDoF = 1;
% number of degrees of freedom
nDim = 2;
% number of dimensions
% L: [l w pl pw]
%
l, w:
domain length and width [m]
%
pl, pw: density of triangles along length and width
% Circular cylinder [x0 y0 r]--[center coordinates, radius] [m m m]
%shape = 'circle'; L = [0.8 0.8 10 10]; c = [0 0 0.2];
% Elliptical cylinder [x0 y0 a b]--[(x-x0)/a]^2+[(y-y0)/b]^2=1
%shape = 'ellipse'; L = [1.3 1.3 10 10]; c = [0 0 0.6 0.2];
% Rectangular rod [x0 y0 a b]-- [x_center, y_center, length, width]
%shape = 'rectangle'; L = [0.4 0.4 5 5]; c = [0 0 0.4 0.4];
shape = 'rectangle'; L = [0.6 0.4 5 5]; c = [0 0 0.6 0.4];
E = 200e9;
nu = 0.32;
mu = E/(2*(1+nu));% shear modulus [ N/m^2 ]
tau_max = 4.0e6; % maximum shear stress
%beta
= 3.5e-4;
% twist per unit length [1/m]
beta
= 4e-4;
% twist per unit length [1/m]
% options for solving plasticity problem by iteration
Niter = 2000; dt = 5e6; plotfreq = 100; penalty_factor = 2e-14;
% set axis limits for plotting
if strcmp(shape,'circle')
clim = c(3);
elseif strcmp(shape,'ellipse')
clim = max(c(3:4));
elseif strcmp(shape,'rectangle')
clim = max(c(3:4))/2;
else
fprintf('Unrecognized shape for plotting.\n');
clim = 1;

1 of 8

4/26/13 1:42 PM

D:\My Documents\Courses\ME342-Inelas...\torsion_plast_example.m

end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Part 2: Generate Mesh, Initialize Arrays, Set Up Boundary Conditions
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Step 1: Generate solution mesh
[Coord,IEN,c,shape]
= CreateMesh(L,c,shape); % Generates the mesh
IEN = IEN';
nNodes
= size(Coord,1);
% Number of nodes
nElements
= size(IEN,2);
% Number of elements
nNodesElement = size(IEN,1);
% Number of nodes per element
nEdgesElement = 3;
% Number of edges per element
% Step 2: Set Tolerance for detecting boundary nodes
boundary_tol = 1e-8;
onSurface = logical(zeros(1,nNodes));
% Step 3: Allocate Arrays
C
= spalloc(nElements,nDoF,nElements);
f
= spalloc(nElements,nDoF,nElements);
g
= spalloc(nNodes,nDoF,round(nNodes/4));
EBC = spalloc(nNodes,nDoF,round(nNodes/4));
% Step 4: Set Essential Boundary Condition
% Step 4.1: Get coordinates of all nodes
X = Coord(:,1);
% x-position of all nodes
Y = Coord(:,2);
% y-position of all nodes
% Step 4.2: Define essential boundary condition value(s)
T = [ 0 ]; % Traction
% Step 4.3: Set g, EBC for all nodes on the surface
% Step 4.3a: Find all nodes on external surface
if strcmp(shape,'circle')
% distance squared of node from circle center
D2 = (X-c(1)).^2 + (Y-c(2)).^2;
onSurface = (D2 >= c(3)^2-boundary_tol);
elseif strcmp(shape,'ellipse')
% distance squared of node from ellipse center
D2 = ((X-c(1))/c(3)).^2 + ((Y-c(2))/c(4)).^2;
onSurface = (D2 >= 1-boundary_tol);
elseif strcmp(shape,'rectangle')
for ii = 1:nNodes
onSurface(ii) = abs(X(ii)-c(1))>(c(3)/2-boundary_tol) ...
|| abs(Y(ii)-c(2))>(c(4)/2-boundary_tol) ;
end
else
error('Code should not get here. Shape should have been set.\n');
end
% Adjust the interior mesh coords
[newCoord] = AdjustMesh(Coord, IEN, onSurface, 200, 0e-1, 200);
Coord = newCoord;

2 of 8

4/26/13 1:42 PM

D:\My Documents\Courses\ME342-Inelas...\torsion_plast_example.m

% Step 4.3b: Set g, EBC for node on surface

EBC(onSurface) = 1;
g(onSurface)
= T;
% Step 5: Define the Natural Boundary Condition
% Step 5.1: Allocate h array
h = spalloc(nEdgesElement,nElements,round(nElements/4));
% In this problem there is not any natural boundary condition applied
% Step 6: Define the Load Value (often RHS of equation you are solving)
% Step 1: Set f for all elements
f(:) = 0.1;
f(:) = 10;
f(:) = 1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Part 3: Create Destination Array (ID) and Location Matrix (LM)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Step 1: Create Data Processing Arrays
% Step 1.1: Allocate arrays
ID = zeros(nNodes,nDoF)';
% Destination Array: ID(A) = P if A is not on essential boundary
%
0 if A is on the essential boundary
%
index is global node number (A)
%
result is global equation number (P)
LM = zeros(nNodesElement*nDoF,nElements);
% Location Matrix: P = LM(a,e)
%
row is degree of freedom (a)
%
column is element number (e)
%
result is global equation number (P)
% Step 1.2: Populate ID Array
I = full(EBC' == 0);
% Creates logical indexing array
nEquations = sum(sum(I));
% Count number of DoF's
ID(I) = 1:nEquations;
% Assign Equation numbers
ID
= ID';
% Step 1.3: Populate LM Array
P
= ID(IEN,:)';
LM(:) = P(:);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Part 4: Assembly of K and F Matrices
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Step 1: Allocate K and F
K = spalloc(nEquations,nEquations,20*nDoF*nEquations);
F = spalloc(nEquations,1,nEquations);
% Step 2: Assemble K and F
% Step 2.1: Compute element contributiontions
for e = 1:nElements
% Step 2.2: Get coordinates of element nodes
x = Coord(IEN(:,e),1);
y = Coord(IEN(:,e),2);
% Step 2.3: Define Shorthand

3 of 8

4/26/13 1:42 PM

D:\My Documents\Courses\ME342-Inelas...\torsion_plast_example.m

4 of 8

x12 = x(1) - x(2); x23 = x(2) - x(3); x31 = x(3) - x(1);

y12 = y(1) - y(2); y23 = y(2) - y(3); y31 = y(3) - y(1);
% Step 2.4: Compute element Jacobian
J_e = [ -x31, x23; -y31 y23 ];
% Step 2.5: Compute element size
elementSize = 0.5*det(J_e);
Ae = elementSize;
%==================================================================
% Subsection 1: Computation of k_e
%==================================================================
% For isotropic case.
k_e = 1/(4*2*mu*beta*Ae)* ...
[y23^2+x23^2,
y23*y31+x23*x31,
y23*y12+x23*x12;
y23*y31+x23*x31,
y31^2+x31^2,
y31*y12+x31*x12;
y23*y12+x23*x12,
y31*y12+x31*x12,
y12^2+x12^2];
%==================================================================
% Subsection 2: Computation of f_e
%==================================================================
ff = f(IEN(1,e));
f_e = ff*Ae/3*[1;1;1];
%==================================================================
% Subsection 3: Computation of f_g
%==================================================================
g_e = g(IEN(:,e));
f_g = -k_e*g_e;
%==================================================================
% Subsection 4: Computation of f_h
%==================================================================
% Step 4.1a: Edge load intensity (homogeneous boundary)
f_h = h(:,e);
% Step 4.1b: Edge load intensity (non-homogeneous boundary)
l12 = sqrt(x12^2 + y12^2);
l23 = sqrt(x23^2 + y23^2);
l31 = sqrt(x31^2 + y31^2);
h_e = h(:,e);
f_h = h_e(1)*l12/2*[1;1;0] + h_e(2)*l23/2*[0;1;1] + h_e(3)*l31/2*[1;0;1];
%==================================================================
% End of Subsections
%==================================================================
% Step 2.6: Get Global equation numbers
P = LM(:,e);
% Step 2.7: Eliminate Essential DOFs
I = (P ~= 0);
P = P(I);
% Step 2.8: Insert k_e, f_e, f_g, f_h
K(P,P) = K(P,P) + k_e(I,I);
F(P)
= F(P) + f_e(I) + f_g(I) + f_h(I);
end
% Global matrices for computing spatial gradients and area of elements
Nxe = zeros(nElements,nNodes); Nye = zeros(nElements,nNodes); Ae = zeros(nElements,1);

4/26/13 1:42 PM

D:\My Documents\Courses\ME342-Inelas...\torsion_plast_example.m

for e = 1:nElements
% Get element coordinates
xe = Coord(IEN(:,e),1);
ye = Coord(IEN(:,e),2);
% Define Shorthand
x12 = xe(1) - xe(2); x23 = xe(2) - xe(3); x31 = xe(3) - xe(1);
y12 = ye(1) - ye(2); y23 = ye(2) - ye(3); y31 = ye(3) - ye(1);
% Area of a triangle:
Ae(e) = (xe(1)*(ye(2)-ye(3))+xe(2)*(ye(3)-ye(1))+xe(3)*(ye(1)-ye(2)))/2;
% Compute element jacobian
J_e = [ -x31, x23; -y31 y23 ];
% Calculate natural derivatives of N1,N2,N3
N1es = 0; N2es = 1; N3es = -1;
N1er = 1; N2er = 0; N3er = -1;
d_Nrs = [N1er N1es; N2er N2es; N3er N3es];
d_Nxy = d_Nrs/J_e;
Nxe(e,IEN(:,e)) = d_Nxy(:,1)';
Nye(e,IEN(:,e)) = d_Nxy(:,2)';
end
% compute volume contribution from each node
for ii = 1:nNodes
[row,col] = find(IEN == ii);
V_ave(ii)= sum(Ae(col))/3; % 3 for triangular elements
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Part 5: Compute Finite Element Solution (including plasticity)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Step 1: Solve elasticity problem
d_slash = K\F;
% Step 2: Solve plasticity problem by
% minimize (0.5*d'*K*d - F'*d)_in_elast_region + (penalty)_in_plast_region
% using steepest descent algorithm
% Plast_region is defined for nodes on elements whose slope exceeds tau_max
% Elast_region is defined for remaining nodes
% For nodes in Plast_region, the (elastic) gradient
% from (0.5*d'*K*d - F'*d) term is set to zero, and the gradient from
% penalty function = (tau^2-tau_max^2)*penalty_factor is added.
% initialize d from elasticity solution (scaled)
if ~exist('d')
d = d_slash*0.25;
end
I = (EBC == 0); u = zeros(size(g)); tau_max2 = tau_max^2;
for iter = 1:Niter,
% grad is the residual of Poisson's equation
grad = K*d - F;
u(I) = d(ID(I)); u(~I) = g(~I);
dxe = Nxe*u; dye = Nye*u;

5 of 8

4/26/13 1:42 PM

D:\My Documents\Courses\ME342-Inelas...\torsion_plast_example.m

% total torque
Torque = 2*dot(V_ave, u);
% find elements whose slope exceeds tau_max
tau2 = (dxe).^2+(dye).^2;
element_factor = (tau2 >= tau_max2);
ind_e = find(element_factor);
id_in_elements = reshape(ID(IEN(:,ind_e)),length(ind_e)*3,1);
non_zero_entries = find(id_in_elements>0);
% for these elements, set the residual term from Poisson's equation to zero
grad(id_in_elements(non_zero_entries)) = 0;
% for elements whose slope exceeds tau_max, add penalty function
penalty = sum( (tau2-tau_max2).*element_factor ) * penalty_factor;
dpenalty_du = zeros(size(u)); dpenalty = zeros(size(d));
for e = 1:nElements,
if element_factor(e)
% derivative of slope wrt nodal value
dpenalty_du = dpenalty_du + ...
(2*(Nxe(e,:)*u)*Nxe(e,:)' + 2*(Nye(e,:)*u)*Nye(e,:)');
end
end
dpenalty_du = dpenalty_du * penalty_factor;
dpenalty(ID(I)) = dpenalty_du(I);
grad = grad + dpenalty;
% move d in the steepest descent direction
d = d - grad*dt;
% plot intermediate results during minimization
if (mod(iter,plotfreq)==0)
disp(sprintf('iter = %d/%d Torque = %e sum(d) = %e norm(grad) = %e', ...
iter,Niter,Torque,sum(d),norm(grad)));
figure(1)
x = Coord;
T = trisurf(IEN',x(:,1),x(:,2),u);
xlim([-clim clim]); ylim([-clim clim]);
figure(3)
for ii = 1:nNodes
[row,col] = find(IEN == ii);
du_ave(ii,1)= mean(dxe(col,:),1);
du_ave(ii,2)= mean(dye(col,:),1);
end;
T = trisurf(IEN',x(:,1),x(:,2),sqrt(du_ave(:,2).^2+du_ave(:,1).^2));
xlim([-clim clim]); ylim([-clim clim]);
drawnow
end
end
% Step 3: Arrange results (d) into solution vector (u)

6 of 8

4/26/13 1:42 PM

D:\My Documents\Courses\ME342-Inelas...\torsion_plast_example.m

u = zeros(nNodes,nDoF);
I
= (EBC == 0);
u(I) = d(ID(I));
u(~I) = g(~I);
Torque = 2*dot(V_ave, u);
% Step 4: Create coordinate array
x = Coord;
% Step 5: compute average slope on nodes
for ii = 1:nNodes
[row,col] = find(IEN == ii);
du_ave(ii,1)= mean(dxe(col,:),1);
du_ave(ii,2)= mean(dye(col,:),1);
end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Part 6: Plotting
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Step 1: Plot Prandtl's stress function (u, i.e. phi)
figure(1)
T = trisurf(IEN',x(:,1),x(:,2),u);
title('Torsion: Prandtl Stress Function',...
'Interpreter','Latex','FontSize',12);
xlabel('$x$','Interpreter','Latex','FontSize',12);
ylabel('$y$','Interpreter','Latex','FontSize',12);
zlabel('$\phi$','Interpreter','Latex','FontSize',12);
set(gca,'FontSize',8);
axis([-clim clim -clim clim 0 max(u*1.1)])
co = colorbar;
set(co,'FontSize',8);
colormap jet
shading interp
set(T,'EdgeColor','k');
% Step 2: Plot sigma_xz and sigma_yz (gradient of u)
figure(2)
subplot(2,1,1)
T = trisurf(IEN',x(:,1),x(:,2),du_ave(:,2));
axis equal
axis([-clim clim -clim clim])
view([0 90])
title('FE Solution for $\sigma_{xz}$',...
'Interpreter','Latex','FontSize',12);
xlabel('$x$','Interpreter','Latex','FontSize',12);
ylabel('$y$','Interpreter','Latex','FontSize',12);
zlabel('$\sigma_{xz}$','Interpreter','Latex','FontSize',12);
co = colorbar; set(co,'FontSize',8);
colormap jet; shading interp
set(T,'EdgeColor','k');
subplot(2,1,2)

7 of 8

4/26/13 1:42 PM

D:\My Documents\Courses\ME342-Inelas...\torsion_plast_example.m

T = trisurf(IEN',x(:,1),x(:,2),-du_ave(:,1));
axis equal
axis([-clim clim -clim clim])
view([0 90])
title('FE Solution for $\sigma_{yz}$',...
'Interpreter','Latex','FontSize',12);
xlabel('$x$','Interpreter','Latex','FontSize',12);
ylabel('$y$','Interpreter','Latex','FontSize',12);
zlabel('$\sigma_{yz}$','Interpreter','Latex','FontSize',12);
co = colorbar; set(co,'FontSize',8);
colormap jet; shading interp
set(T,'EdgeColor','k');
% Step 3: Plot maximum shear stress (magnitude of gradient of u)
figure(3)
T = trisurf(IEN',x(:,1),x(:,2),sqrt(du_ave(:,2).^2+du_ave(:,1).^2));
xlim([-clim clim]); ylim([-clim clim]);
title('Torsion: Maximum Shear Stress',...
'Interpreter','Latex','FontSize',12);
xlabel('$x$','Interpreter','Latex','FontSize',12);
ylabel('$y$','Interpreter','Latex','FontSize',12);
zlabel('$\sigma_{xz}$','Interpreter','Latex','FontSize',12);
cbar = colorbar; set(cbar,'FontSize',8);
colormap jet; shading interp
set(T,'EdgeColor','k');
% Step 4: Plot contour of stress function with elast-plast boundary
figure(4);
lxi = linspace(-1,1,300);lyi = linspace(-1,1,300);[xi,yi] = meshgrid(lxi,lyi);
ui = griddata(x(:,1),x(:,2),u,xi,yi,'cubic');
dui = griddata(x(:,1),x(:,2),sqrt(du_ave(:,2).^2+du_ave(:,1).^2),xi,yi,'cubic');
[c h]=contour(xi,yi,ui,20);
hold on
c2=contourc(lxi,lyi,dui,[0.95 0.95]*tau_max,'b--');
plot(c2(1,2:end),c2(2,2:end),'r.');
hold off
title('Torsion: Prandtl Stress Function',...
'Interpreter','Latex','FontSize',12);
axis equal; axis([-clim clim -clim clim])
figure(5);
xi = linspace(-1,1,300);yi = linspace(-1,1,300);[xi,yi] = meshgrid(xi,yi);
contour(xi,yi,dui,20);
title('Torsion: Maximum Shear Stress',...
'Interpreter','Latex','FontSize',12);
axis equal; axis([-clim clim -clim clim])

8 of 8

Plastic solution
5

Torsion: Prandtl Stress Function

x 10

x 10

7
5
6
5
4

4
3
3
2
1
2
0
0.3
0.2

0.3
0.1

0.2
0

0.1
0

0.1
0.2

0.1
0.2

Torsion: Prandtl Stress Function

0.25
0.2
0.15
0.1
0.05
0
0.05
0.1
0.15
0.2
0.25

0.25 0.2 0.15 0.1 0.05

0.05

0.1

0.15

0.2

0.25

Plastic solution
Torsion: Maximum Shear Stress

x 10
4

x 10

3.5

4.5
4

3.5
3

2.5

xz

2.5
2

1.5
1

1.5

0.5
0
0.3

0.2

0.3

0.1

0.2
0

0.5

0.1
0

0.1

0.1

0.2

0.2

FE Solution for xz
0.3

x 10
4
3

0.2

0.1

0
1

0.1

0.2

3
4

0.2 0.1

0.1

0.2

FE Solution for yz
0.3

x 10
4
3

0.2

0.1

0
1

0.1

0.2

3
4

0.2 0.1

0.1

0.2

C^f^n&"*l1T*lt{'

VurW
-fb

C"l

,Xo^ftu. I Ir.&nxi l"p"n^oh2- rna\<q/

a oa/{-iwf LD Plo*-5lvcu},t

prrbtrry arfrrr',t^
zD plaru"-1hrc,;"
tf.a-pt,rtfu^5r$t{, ths'vr{iii" rir,'zrue*L
in *hg ^qD*bJwr n^t{"
tN\\b0dnlcnuAed

gl, Probt
- S+dr'nl,^i-

i'

h15;'*,, o\ +{,;(ktur.(toAW"MI-L

{ffi

rtttv nW,"*L (7atl

+-r^tgQretxhtnwvo,rn(-S

rrwr..T,
p*er4in^&4,rL$*knua + plw>htY*
vrrfL inrntx*lt!. f,
ccr,v,a'f
7lotrtt' flr^) o
rod^ud 0r*tc 6 '34rlo'

D+tewia'-t if

Wu

+he'
tyltn 6 a{y1ir'/L'lf- f, ' f" .,
?la/v4L
WO

b"/nxl.Lf,r b ^ thlP^hHJ v"Lut'' I nntbr '+ y'e| 1

lhfip- rS o.!*s a znoNmww'r^l^v

fffrWw

/M{9 w;!{b''ur'et!

{rnl
, ? P^*,

/r*

, 6'r '#*l'-

tU

}A^en,ll

f,

iI

Ete"rtry'
-$z' -?Lane--stt"'.
Sohtf**
l,lV=0,
*\WdJ,,vf t' Z { , ' fh = o )
ux*
L"rb+u' U x ( X ' 1) ) \r,to'tY
uy(x,y) F*yr,{
)

f'ry

4a
6aa
/

t)

V z -

cz?:-Ew*

tr

6vv+; e,*+
:O

6-*e : r't(orx t 6yf )

-)

LWb &uvtw r,s/rdn{,1rn, 6xx,.t: f 6yr, y i' F- - o

oxy, t,
F

c.;lnTaJthU6

+ 6yy,y 'f r/ :
,

osnd;{i\sn' &*,,,1y+ Eyy,xx .- L txy, y :

H*nka+La$l,
C,p.nu^la.,./

Ex--+crnc.-ryayy
yy= -rry\ 6*-x-r
txy :

l-t+yp d,ur* A.-rJ S+'4^4 6^^'ti;r)sr\,

d ] J L x -h , w

nYl=<l',,xtc

.E:>/( +u)

4 t1 , ll )

, tn*h tf^th#

Q)t--fr.y

Th-a/yL
iv)tht- o,,#+*nw"f ,n"ry*,.{ t l=r:7:,
ho&'4+,:.^ ,> a^fu1^h'u$,

tA rf'!,h)wm
bo
W

rAz

6yy

6!

*,

c-r'r d''r*fk-nr-' .Lec''rrnl-'r

).
so,rl+lf,U,d

( b; h NLm'sn/f C -{t'\xqrn )

V++ = f,** + {,: rvr/ - 7 {'*ry

-- o

= fndne l"rt$ftel L;hu'*rvicr{*rr*,

cvttdunsvo
L^ poloq,
, +(Y,1!,))
^D
M k^y
ll
(;rn

t4E3+o UVi",|W

*,,

,n7: {), - n,
k-,,**',vr-!- Na\pg

FD. ytobb',w: vath M, |*-!- syf;nn"4ltft t

K'r{:

tDd

T*

9.*n,, 7- n .

, " Pol,^4 C':ctil"

7==-o,

t' '*te4

q
)

(ntWt{

* = 4 , r z r A zr L l ^ r + A s / - r + A + 9

+ l;#,,

- *(+-?9,), %- &_
dic>--=
drz

so{,vr:.r,,)

Merlzl-2

W&^LA

Cu"

1l'bL

frr : A , - 2+ A , ( r l r r t l ) + A r . f

A+t

0'-r0

ctu :

1-3) + & ( T)
h t ' 7 + hL'(r-L^t
d-fl' t--ir*u" 'W
Boi^

\i s - .I'r=^
. = -1,

t\

6Yt

t-=L

I
I

O-.re
t:b

l"
f f a , , =A=; : O , A \ : L

l
1

_. \.,

I sru \r:or ==o

t
l

hrJre

pre/twr^trl$

:Q

0'
tLA" ,,

A:=#

,/b

J"oJ (

6rr

il-u,'

V-* \
I

C rg:'

, l;,)

i:)

6ee:+*( r+i"l
'\Qrt.*
6t* -V(6tr
+ Qe )

6j :
-=

g :

/-

* Ae.z.)
T(6M t0rr
,

t<-

+ Grr+dso):ry
ffi

"=,"ffli&
*!:)

J
s,=6>-A=&t-,(ry - ?b't

-\*-a'/
L3

TA\-3

!'':'-1)

ru{'tn
K*

as'hnp ?,= M
b'-a"

.r':f'fry*T j

r*-l

D' I

":,',,,1FFF^'"

, r ,|

s r =6d=f f i 3 , ' + )

5j: 6t-6= !#

$6-

ffi
',?-

-z /',il
Su+x/

- J"'(F
+)
itt
ars*"

]
f

"*

-!"-( r*

{ l

'f+)

b"lk,wt^J,wJ
(t--zv)'

i7 ---

/ ,,

/t r

/^

'i

/ ^

Ut*Lmd'rt{a't' luJg

- MQ3*2'

, ( dl
L_*-_-

+'_l

+na'r{rnr*'n T,

At *nr "{

6rI

tz
t
wJ: t

g+yYz

er)
Lry'J_

c\
'Yv'
F6{ ?o > ?'Y , th-L+-u-btun:f*'ha
fl"v'I-L
^' talt'* yeo'tt't
'rnfig*'t^":;c&-!'*U
+*- 0 b
S;,,'yyos*f , irtwe4
(r:-a)
th* fJ Tt^p" 4* plarfrr- reffw sb4fuJiot +lt!-inrtut r^rnJ,t
l'nd' qrru0dr;rlql,

l2gdrue,(6,ytL*/*ebwth*ry

Lq,t t:f

Tha eh*fi2

^4 lt
c,rytdn*f?Ena

bo"-@

D
17
/
/t-

1s
b
'iz -

'

r
|
j

\\

fr,

r:b

.Q

\{"J

houl l':ltu^(cu{

t'12-fiav'- ! < f <b

Y:r(

f,tu

l r =b

- v

5' ilo sYrnftrid in +hL elw>firz r-n y".m,..tb

Yt
^-.
vr(

g{ant-t

r-L

lf .lt" [ r - L \
- -F:"-lV
.

f,i

Sirn

fi,g:Q
r^t
P^ F't'
rpo= #_r,
(t+ +)
'
b''r

I'l y=f

.fePt

r : -, ^ z t ) ) ' , b * l

b-1"
I , :k r y [ L P ' +
T/"b fo+l\L^T^fuJ,

ns,*J s*uryln, | "r th-Ltler?tc,*t- brrl*.

+ $,trt

{1{,;n{'*tr-L T"ry

YlE::3*Z
ot-a{tne,T"
.+l

!tl-LrrL q4
f ?

iI

'vv b

Co;

f ft-rtr)' '
= &* +
:
F
i
L : Y T
-

rhL *ia&Er*nyitn

?",(r-+,)

lfre= o
I
t f r e o =t " ' ( l +$ l

L+ l-V'

r < {* b .

s'hufo"rb
I vw1
)

t
t

rh
\\ {agua&x
J

- t"*\ a. t
pr+[,t .,Uc,retu!{Y
rYtovqp'u,l
l''
Io, f

*h"\

F.3 /
l,-

b,

I EV+L
t)n
M--uluuuu'cr-[ J"{n^"f,'o""t

Qq

'tlns
Vtli-W daSUtrW;wU yl"*"
iA ttA!-eh,yliL tY.

a4 duut'wd a,6ilr<- ca-,v be o6l-rl*A"

1o/,hwc i,r ftL
tLw rv.m.s.rLi:t^!-

tt4 F

pr,f*%

Th,^

tQ*,\ver-tL,X ,,,,"9!u hN^) +t"a \ome. 5"1"finc

vuJ/tti/tA'ry

Tlrv>+tz

Ef w;!"'briwvvLc''r,d;t'lV"t r\' Folan' atttdn'nr&A

76v'

> S v e+ 6 v t - 6 N
+
I
>o
Y

\-t"T

II r ?qqe+

b6ya
+
A Y

bB

Fr_e

.{. Fo : o

ag:Y

fro :o,

ln ulrv W*flt-.,,
I
)6/y

f r r
-fh!-

6yy -

#:o

fr' *" Fo4,

600

t(eq;0; b rivwtt r^r"lr+*^)

A% a':d;*'",r- (*tt bv u6+rr';,^ud@

-Efb

tl"t''

S4rwovt( ,'n ? o (c,nt cntzd) na;tu )

lQ,=kw
1 t e &+=L u ' * # )

[ 8 , , tu(=i ' # * W - V )
, |n ,>*'<xa*1lr,, l'rohtured,aylau,,^ux
fVH' U{"

Ir,t ,"f**^t

(Lr': *"

I a" = l-ur

\ f"e:o

Wt

o!^+'
0

ort-L r1,on-7er?'o

ilrain ln lslasaaJxE = |/ (?*'* %o+ tzx) =

'

t(#"Y)

- - L) Yr ru ew trr !). t - - , Llal * o +Y W

?t-eq=
?y

ay

- s*)
(e"
t

( c"1^nb|1"7c*n&r*,a)

Mv',?42

ry^h;L^L

Cyt,^tt-rl^'w d

l4vkaL

frr =(Ao*)L r r + ^

-T*6e

L-"/^) i

tw

t-,n + ()+z1,tl [o,g

fgt:

rr*=

^ tw + A E*e

frr - fo<g: zV (W- funr)

PU^tur^l,c eau,ifihrn,notcsdt$Vn
0r

lQ*rru)*
)o U,rcnW*- tN

l)0,*dr,"{

egwd:,:rws{-fv tM AtlO t"vtk*afi*rvuJ ,

*b

( 2+ >J^)

=q

[:''

*,

toy .

osn&rl+ans -.

+ + ( Qt-s'o): Q

. 7W

I
i o d l : - 1"
'y=o

I':ttr

q*|, --, + i f " " /

:12

Marflo'.$UdJ- tA,"t f,n&

-o

I l:1,

0,,,l.f ) tu'1dtee(v) Ihat sun'tff

q{."*'rin
tt^rr r-r"fa,fi.b-l^^6 rO1 r,; A;bnb,n-,,n " J b a-,vd*Ul 6"t^4
is wvttle-y'r

or6,UVJ-,'(rrJr<-*Q-Lart'f.r'")
.z A{ns'"yb^J-+tdeLlasf'n

Ll rna@ u\,w ,+ Mr,Yt*h f s,l,,e f*.**

8
Elastic solution
0.8
num

rr

num

0.6

anl

rr

0.4

anl

/ p0

0.2
0
0.2
0.4
0.6
0

0.5

1.5
r

2.5

2
num

rr

num

1.5

anl

rr

anl

/ p0

0.5

0.5

1
0

0.5

1.5
r

2.5

-Tt
Orf"hd,v,tn'L t)-

vEv*L

Cc^t

Renr
?s Plae+r<,
-b { t),e lrl[o.,tlr'n
[*fuiru. ?" orJf
(f il tl'rst"J*
b**dry ) ,^,|^s,n
lr> lY )
"f ?,luDd't-plna{,r,

0,,'nrnXt {a.Ak uo k

q^rd"h
*nsl sfu;^ lrnld hsr& {h,t plat+ic'ttf'lttl L. f (J
fr'nlt+e ehq,'rrl
!/tl shJt chr"nt a,w wt{4'*lr\t' ^4

(.,rr, too

r-i-__.

+"il

Sry

otevia#nrtuJtF44

S{ra,\

anJ-gso,lrh:h3 PDE'5.
.f i*tnana cwnbg uyfc,'rn,r
d In lur*s ,f fha4 f l"m&.
Att '4'!r(tusm{Thta
S$m;"^
,
Hyero$o+cc

E : j i [,n-t tr"l

dg"rit"tOAC Sta,1rn.

e r r : e r r - -E :
&e:

notq- zz:o

+ tr - | t"s

ru*{ = - ) t r y * i t " o

ti* = t*.- g = - | ? , w - f t $ o
l-lvt^o$t*lt- s+rQrl
I

df^,&+nrlL Sl-'rur].

o - - 3 K i 1 = K (z,v -t ?so)
Srrl

5'{i'0, S

lry t 5o++ hz : e )

'

Svyl u Ct rt loo)

4;z
in fu- plat

: k'
) = sft + tJ, + S,r-!oo
Jz: | ( 5iritSdt sJa

t j-4k,.-:tr)
ie= :,i(--s,,

W]-t'il,&?
Ird,6 ,

Qv= 9rr+ 6' == S,tf K (n t Ue)

trso= SOp+
Ci'
6r, - G := J'rr - 569

G,t,{no
ilr'c^t
a

METM

'

lwb-

I n :

Cw'

ll (v 1j i
--+*

Ws nln/osntr,1 thL ) ?DV'sf"/ U, tgo,Srv

76rr
VY

lpaftbrfuYncqt*{;an,

qr-60 ==o
Y

h/r.o bgOStY\l4

-----o

: o

+rK(41)'*$)+
bhfy rx'ir{r+isn'
CrnTod''

Afue

w=' +(trv-too)

-,-"@

r
P[or,t'cth,J r,^tLwit[ Le"{t -fh,{"ilunr\ e{uro.'hrn
|
2
-(i,iff
- f -rtl
L rzefu\3
=
Af

-P
r,,
i W)0,] #'* ?

?),50,

-l

pl)E Nith,
#

tute-tW {/ tl.,a-4

TVeg >\ dP-veJ )-^ P'to

Gry*

ttw- e(anfiWn;'^nVlLth,eg-i* d'! elu*it fadlne(6+l

Tl,,q-,c^',^f"atf
B,*

b{-f

c-.ndratryv d

tltl

J ^/v'12'

eq4l;$to}r"4 c^m0Yarrr fue,-rt-ca.nr'pf be- b/61ten ,rr Js'rrTra "J

tw , tse o,t^r,re- ( Lnc**w
,r,,/' fLa- /a,Hiz" g**r,v. )

tl**

Jy19+*uA, S*

3" ftse ,i

|'-murvv a'nd- fla

D4 i,1+wd^LLedlre,r*t, (ltz.t d,rz 3'd u..t't

sy1hed

tlw fil,*z
-f

f kN Ywk'

h-0.equa,ltrw' ,y*

rueo{t" Lz

Soliel;n ths 7T syaw I

Y *,,i !'

T[a Sol*v-n r]r the, ela*t Yefr;wtiS

Alrt^['Y lcr*t"tn'
a

lfr
Cyt^^a$c^I-

_W-l*+L

-i C*'

tu=Lsv,
* = *r*
el!: }'t*

------J+

(?h4tk
i{*'"J-)

L'lo=*u'u='
dne^r+rn M 4h!- dturNhruzsMJ
S+Tari.nA rl'rr-{/nl-So,"rrs-

Tlw dgA'o&al^yg.btn,
nel
t
vlt=*

( ;4=4:in
vvt-z;t)
I
'

)?Y,y

; 4 -fi,S*
r
, . i ffi
-.-=)'
[

- & ---)^
) - ,*
c^
co
-el
| .
eE = ;y Jzz

ert -- i^ 57s'=o
',To,{nQ-g.ka,>- rate

(rtstirfr,iu

pwa-:t)

- " ' = eI n t{! =,;!l:lt( Srr+}svy) I

f ,,
",^'nl )
l
i
- L) ) cfst')
^'[(r i*+
:
)
=
*
s6
?"9
A{}
4L,tL
/

^
o[^v'in'r'

)a,

( *r= i#+il::o
:= (,/ e*o - iro) f*
- i'" 5o+1= T,S, sury
erv
)
Qr^
S0o J
,)"t@iBo - eo. \'rr) : S;.yJg4r
"r
;

+W

"

^'l rdfu' P

bouild"^uffis **V
+he vtontit'-Sa!.,u>tiu
a

P
tht- {-t t,t.a4-rsrtb"*
a

?trrd*

n'ffQ- _g*
-_- /\
/r
J

#/*

vSp* n
zJ,(
(:# so$tryr*,)= ryroo-*
w){r

yl'ugin

ry=

*e"

Ne atr44-ok

7 too,

P?E A

e1e:?*- ltvr
o',,',
P,?

a*,
dP
I

IZ

Co.,'

t!t^n&ac^'l* Tt'b{-

ME;3+L

96. N*,w;coLS,6+r"*'w P!yry_!Ef:_

-

tP,,
,g,
,s,?.,
-',1'rf;Y:''-,"'c
"r
I [ - -);',g,
,g)
- ' - ' , ' tfl.
,'l

tlht Sol"krn ,b ^lnv&*

S^f p*

s,9,)
^fi-(r, f ) , (t?,rfr,,
uw,rr,a^d ( rrdr,f+af), (8,t3,s3,S
a,*
f
Hu* W oQ)uiLQ.
**{lr^
u
^tr
,
f"d +M soh^#'he
,
(y,
f+df) , (t*. tN,fn)ri

AY

l2P? f,sgrvzl'5,

fl.L da'swe-,"f
ited
c.n&l+m,
eq*;Abo^r^
'

r co
17-9v
-T+K|,ff+

urrtpo,ttlrl
W

"

(sS*Joo.)l

w)J-o

c^"rrd^*nn

*-&o
Ar

.
flo'a*,t+(N*,14,"t-L

ffi#r

,,nLl(b. 3a " - e , 9 _ r
4P J

: o

. n O
aso- as

$)(q=iif-(;
o('

!,,-sP h*_!iil + sso t$

z - r - a T *

r " O
)Yvt)tr

,-

= i (- S,,n W-t;,:,ii:)
Soo
f\[rr^o; C,& sol^^+fr,., o{+*^*;- by

- Plaat'n,
o@d- +*,\'e

e qn5
- ?los-}:'n.
- 0ylrn,l=+".{lt,

:o

13
Plastic solution
0
0.1

rr / (2k)

0.2

1.0

0.3
0.4

2.0
/a = 1
1.2
1.4
1.6
1.8
2.0

0.5
0.6
0.7
1

1.2

1.4

1.6

1.8

r/a

1
0.9
0.8
0.7
/ (2k)

2.0

/a = 1
1.2
1.4
1.6
1.8
2.0

0.6
0.5
0.4

1.0

0.3
0.2
1

1.2

1.4

1.6

1.8

r/a
0.3

2.0

0.25
0.2

zz / (2k)

0.15
0.1

1.0

0.05
0

/a = 1
1.2
1.4
1.6
1.8
2.0

0.05
0.1
0.15
0.2
1

1.2

1.4

1.6
r/a

1.8

4/29/13 10:04 AM

D:\My Documents\Courses\ME342-Inelastic...\cylind_tube_elast.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2D Example: solve strain field in pressurized cylindrical tube
%
elastic deformation only
%
% Wei Cai [email protected]
% Created
03/24/2013
% Last Modified 04/29/2013
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
mu = 100;
nu = 0.3;
lambda = 2*mu*nu/(1-2*nu);
a = 1; b = 2;
dr = 0.05;
r = [a:dr:b];
p0 = 1;
eps_rr = zeros(size(r));
eps_qq = zeros(size(r));
Niter = 10000;
plotfreq = 1000; % 50;
param = [a b dr mu nu p0];
F0 = eqns_cylind_tube_elast([eps_rr, eps_qq],param);
%options = optimset('Display','iter','TolFun',1e-8);
options = optimset('Display','off','TolFun',1e-8);
sol = fsolve('eqns_cylind_tube_elast', [eps_rr, eps_qq], options, param);
eps_rr = sol(1:end/2); eps_qq = sol(end/2+1:end);
F1 = eqns_cylind_tube_elast([eps_rr, eps_qq],param);
%
%
%
%
%
%

(Alternative method using cgrelax, require gradient evaluation)

[U dU_depsrrqq] = grad_cylind_tube_elast([eps_rr, eps_qq],param);
acc = 1e-6; maxfn = 40000; dfpred = 1e-5; n = length([eps_rr, eps_qq]);
[U,eps_rrqq_rlx] = cgrelax('grad_cylind_tube_elast',n,acc,maxfn,dfpred,...
[eps_rr,eps_qq],param);
eps_rr = eps_rrqq_rlx(1:end/2); eps_qq = eps_rrqq_rlx(end/2+1:end);

% stress
sig_rr =
sig_qq =
sig_zz =

field
(lambda + 2*mu) * eps_rr + lambda * eps_qq;
(lambda + 2*mu) * eps_qq + lambda * eps_rr;
lambda * (eps_rr + eps_qq);

pp=p0*a^2/(b^2-a^2); E = 2*mu*(1+nu);
sig_rr_anl = pp*(1-b^2./r.^2);
sig_qq_anl = pp*(1+b^2./r.^2);
eps_rr_anl = (1-nu^2)/E*sig_rr_anl - nu*(1+nu)/E*sig_qq_anl;
eps_qq_anl = (1-nu^2)/E*sig_qq_anl - nu*(1+nu)/E*sig_rr_anl;

1 of 2

4/29/13 10:04 AM

D:\My Documents\Courses\ME342-Inelastic...\cylind_tube_elast.m

% plot solution
fs = 17;
figure(1);
plot(r, eps_rr /p0*mu, '.', r, eps_qq/p0*mu, 'd', ...
r, eps_rr_anl/p0*mu, r, eps_qq_anl/p0*mu);
xlim([0 max(r)*1.5]);
set(gca,'FontSize',fs);
xlabel('r'); ylabel('\epsilon \mu / p_0');
legend('\epsilon_{rr}^{num}','\epsilon_{\theta\theta}^{num}',...
'\epsilon_{rr}^{anl}','\epsilon_{\theta\theta}^{anl}');
figure(2);
plot(r, sig_rr/p0, '.', r, sig_qq/p0, 'o', r, sig_rr_anl/p0, r, sig_qq_anl/p0);
xlim([0 max(r)*1.5]);
set(gca,'FontSize',fs);
xlabel('r'); ylabel('\sigma / p_0');
legend('\sigma_{rr}^{num}','\sigma_{\theta\theta}^{num}',...
'\sigma_{rr}^{anl}','\sigma_{\theta\theta}^{anl}');

2 of 2

4/29/13 9:54 AM

D:\My Documents\Courses\ME342-Inela...\eqns_cylind_tube_elast.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2D Example: solve strain field in pressurized cylindrical tube
%
elastic deformation only
% Construct the equations to be solved by cylind_tube_elast.m
%
% Wei Cai [email protected]
% Created
03/24/2013
% Last Modified 04/29/2013
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function F = eqns_cylind_tube_elast(vars, param)
% F(1:N-1): equilibrium equation
% F(N:2N-2): compatibility equation
% F(2N-1:2N): boundary condition
a = param(1);
b = param(2);
dr = param(3);
mu = param(4);
nu = param(5);
p0 = param(6);
N = length(vars)/2; F = zeros(1,2*N);
eps_rr = vars(1:N); eps_qq = vars(N+1:end);
lambda = 2*mu*nu/(1-2*nu);
r = [a:dr:b]; r_avg = (r(2:end) + r(1:end-1)) / 2;
% stress
sig_rr =
sig_qq =
sig_zz =

field
(lambda + 2*mu) * eps_rr + lambda * eps_qq;
(lambda + 2*mu) * eps_qq + lambda * eps_rr;
lambda * (eps_rr + eps_qq);

dsigrrdr = (sig_rr(2:end) - sig_rr(1:end-1)) / dr;

sig_rr_avg = (sig_rr(2:end) + sig_rr(1:end-1)) / 2;
sig_qq_avg = (sig_qq(2:end) + sig_qq(1:end-1)) / 2;
% equilibrium equation
F(1:N-1) = dsigrrdr + (sig_rr_avg-sig_qq_avg)./r_avg;
% strain field
depsqqdr = (eps_qq(2:end) - eps_qq(1:end-1)) / dr;
eps_qq_avg = (eps_qq(2:end) + eps_qq(1:end-1)) / 2;
eps_rr_avg = (eps_rr(2:end) + eps_rr(1:end-1)) / 2;
% compatibility condition
F(N:2*N-2) = depsqqdr - (eps_rr_avg-eps_qq_avg)./r_avg;
% boundary condition
F(2*N-1:2*N) = [sig_rr(1) + p0, sig_rr(end)];

1 of 1

5/10/13 8:39 AM

D:\My Documents\Courses\ME342-Inelastici...\cylind_tube_plast.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2D Example: solve strain field in pressurized cylindrical tube
%
pressure large enough to cause plastic deformation
%
% Wei Cai [email protected]
% Created
03/24/2013
% Last Modified 04/29/2013
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
mu = 100;
nu = 0.3;
k = 4;
% yield stress / sqrt(3)
E = 2*mu*(1+nu);
% Young's modulus
K = 2*mu*(1+nu)/3/(1-2*nu); % bulk modulus
a = 1; b = 2;
dr = 0.02; drho = dr;
r = [a:dr:b]; rho = r;
p0 = 1;
S_rr_sol
S_qq_sol
sig_rr_sol
sig_qq_sol
sig_zz_sol
eps_rr_sol
eps_qq_sol

=
=
=
=
=
=
=

zeros(length(r),
zeros(length(r),
zeros(length(r),
zeros(length(r),
zeros(length(r),
zeros(length(r),
zeros(length(r),

length(rho));
length(rho));
length(rho));
length(rho));
length(rho));
length(rho));
length(rho));

% analytic solution in elastic region

for i=1:length(r),
for j=1:length(rho),
% limit to elastic region
if r(i) >= rho(j),
ppp=k*((1-2*nu)^2/3 + b^4/rho(j)^4)^(-1/2);
sig_rr_1 = ppp*(1-b^2./r(i).^2);
sig_qq_1 = ppp*(1+b^2./r(i).^2);
sig_zz_1 = nu*(sig_rr_1 + sig_qq_1);
S_rr_1 = sig_rr_1 - (sig_rr_1+sig_qq_1+sig_zz_1)/3;
S_qq_1 = sig_qq_1 - (sig_rr_1+sig_qq_1+sig_zz_1)/3;
eps_rr_1 = (1-nu^2)/E*sig_rr_1 - nu*(1+nu)/E*sig_qq_1;
eps_qq_1 = (1-nu^2)/E*sig_qq_1 - nu*(1+nu)/E*sig_rr_1;
S_rr_sol(i,j)
= S_rr_1;
S_qq_sol(i,j)
= S_qq_1;
sig_rr_sol(i,j) = sig_rr_1;
sig_qq_sol(i,j) = sig_qq_1;
sig_zz_sol(i,j) = sig_zz_1;
eps_rr_sol(i,j) = eps_rr_1;
eps_qq_sol(i,j) = eps_qq_1;
end
end
end

1 of 3

5/10/13 8:39 AM

D:\My Documents\Courses\ME342-Inelastici...\cylind_tube_plast.m

2 of 3

% find solution in plastic region

disp('find solution in plastic region...');
for j=2:length(rho),
disp(sprintf('j = %d / %d',j, length(rho)));
for i=j-1:-1:1,
% find solution at (r, rho+drho)
eps_rr_1 = eps_rr_sol(i,j-1); eps_qq_1 = eps_qq_sol(i,j-1);
S_rr_1
= S_rr_sol(i,j-1);
eps_rr_2 = eps_rr_sol(i+1,j); eps_qq_2 = eps_qq_sol(i+1,j);
S_rr_2
= S_rr_sol(i+1,j);
param = [r(i), dr, mu, nu, k, eps_rr_1, eps_qq_1, ...
S_rr_1, eps_rr_2, eps_qq_2, S_rr_2];
%options = optimset('Display','iter','TolFun',1e-8);
options = optimset('Display','off','TolFun',1e-8);
%F0 = eqns_cylind_tube_plast([eps_rr_1, eps_qq_1, S_rr_1], param);
%sol = fsolve('eqns_cylind_tube_plast', [eps_rr_1, eps_qq_1, S_rr_1], ...
%
options, param);
[sol,Fval,exitflag] = fsolve('eqns_cylind_tube_plast', ...
[eps_rr_1, eps_qq_1, S_rr_1], options, param);
switch exitflag
case 0, fprintf(['Solution incorrect: ' ...
'Maximum number of function evaluations or iterations reached.\n']);
case -1, fprintf(['Solution incorrect: ' ...
'Algorithm terminated by the output function.\n']);
case -2, fprintf(['Solution incorrect: ', ...
'Algorithm seems to be converging to a point that is not a root.\n']);
case -3, fprintf(['Solution incorrect: ', ...
'Trust region radius became too small.\n']);
case -4, fprintf(['Solution incorrect: ' ...
'Line search cannot sufficiently decrease the residual '...
'along the current search direction.\n']);
end
%F1 = eqns_cylind_tube_plast(sol, param);
eps_rr_sol(i,j) = sol(1); eps_qq_sol(i,j) = sol(2); S_rr_sol(i,j) = sol(3);
S_qq_sol(i,j) = (-S_rr_sol(i,j) + sqrt(4*k^2-3*S_rr_sol(i,j)^2))/2;
sig_rr_sol(i,j) = S_rr_sol(i,j) + K*(eps_rr_sol(i,j)+eps_qq_sol(i,j));
sig_qq_sol(i,j) = S_qq_sol(i,j) + K*(eps_rr_sol(i,j)+eps_qq_sol(i,j));
sig_zz_sol(i,j) = -S_rr_sol(i,j)-S_qq_sol(i,j) ...
+ K*(eps_rr_sol(i,j)+eps_qq_sol(i,j));
end
end
% verify solution
disp('verify solution in plastic region...');
for j=2:length(rho),
for i=j-1:-1:1,
% find solution at (r, rho+drho)
eps_rr_1 = eps_rr_sol(i,j-1); eps_qq_1 = eps_qq_sol(i,j-1);
S_rr_1
= S_rr_sol(i,j-1);
eps_rr_2 = eps_rr_sol(i+1,j); eps_qq_2 = eps_qq_sol(i+1,j);

5/10/13 8:39 AM

D:\My Documents\Courses\ME342-Inelastici...\cylind_tube_plast.m

S_rr_2
= S_rr_sol(i+1,j);
param = [r(i), dr, mu, nu, k, eps_rr_1, eps_qq_1, ...
S_rr_1, eps_rr_2, eps_qq_2, S_rr_2];
F = eqns_cylind_tube_plast([eps_rr_sol(i,j), eps_qq_sol(i,j), ...
S_rr_sol(i,j)], param);
disp(sprintf('i = %d j = %d F = %e %e %e',i,j,F(1),F(2),F(3)));
end
end
% plot solution
fs = 17;
figure(1);
skip = find(r/r(1)==1.2,1,'First')-1;
plot(r/a, sig_rr_sol(:,1:skip:end)/(2*k), '.-');
set(gca,'FontSize',fs);
xlabel('r / a');
ylabel('\sigma_{rr} / (2k)');
legend('\rho/a = 1','1.2','1.4','1.6','1.8','2.0','Location','SouthEast');
figure(2);
plot(r, sig_qq_sol(:,1:skip:end)/(2*k), '.-');
set(gca,'FontSize',fs);
xlabel('r / a');
ylabel('\sigma_{\theta\theta} / (2k)');
legend('\rho/a = 1','1.2','1.4','1.6','1.8','2.0','Location','NorthWest');
figure(3);
plot(r, sig_zz_sol(:,1:skip:end)/(2*k), '.-');
set(gca,'FontSize',fs);
xlabel('r / a');
ylabel('\sigma_{zz} / (2k)');
legend('\rho/a = 1','1.2','1.4','1.6','1.8','2.0','Location','SouthEast');

3 of 3

4/29/13 10:07 AM

D:\My Documents\Courses\ME342-Inel...\eqns_cylind_tube_plast.m

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 2D Example: solve strain field in pressurized cylindrical tube
%
pressure large enough to cause plastic deformation
% Construct the equations to be solved by cylind_tube_plast.m
%
% Wei Cai [email protected]
% Created
03/24/2013
% Last Modified 04/29/2013
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function F = eqns_cylind_tube_plast(vars, param)
% F(1): equilibrium equation
% F(2): compatibility equation
% F(3): plastic flow rule
eps_rr = vars(1);
eps_qq = vars(2);
S_rr
= vars(3);
r = param(1);
dr = param(2); drho = dr;
mu = param(3);
nu = param(4);
k = param(5);
eps_rr_1 = param(6);
eps_qq_1 = param(7);
S_rr_1
= param(8);
eps_rr_2 = param(9);
eps_qq_2 = param(10);
S_rr_2
= param(11);
K = 2*mu*(1+nu)/3/(1-2*nu); % bulk modulus
F = [0 0 0]';
S_qq
= (-S_rr
+ sqrt(4*k^2-3*S_rr^2 ))/2;
S_qq_1 = (-S_rr_1 + sqrt(4*k^2-3*S_rr_1^2))/2;
S_qq_2 = (-S_rr_2 + sqrt(4*k^2-3*S_rr_2^2))/2;
dSrr_dr =
dSqq_dr =
depsrr_dr
depsqq_dr

(S_rr_2-S_rr)/dr;
(S_qq_2-S_qq)/dr;
= (eps_rr_2-eps_rr)/dr;
= (eps_qq_2-eps_qq)/dr;

depsrr_drho
depsqq_drho
dSrr_drho
dSqq_drho

=
=
=
=

(eps_rr
(eps_qq
(S_rr (S_qq -

- eps_rr_1)/drho;
- eps_qq_1)/drho;
S_rr_1)/drho;
S_qq_1)/drho;

1 of 2

4/29/13 10:07 AM

D:\My Documents\Courses\ME342-Inel...\eqns_cylind_tube_plast.m

% equilibrium condition
F(1) = dSrr_dr + K*( depsrr_dr + depsqq_dr ) ...
+ ( (S_rr_2+S_rr)/2 - (S_qq_2+S_qq)/2 ) / (r+dr/2);
% compatibility condition
F(2) = depsqq_dr - ( (eps_rr_2+eps_rr)/2 - (eps_qq_2+eps_qq)/2 ) / (r+dr/2);
% plastic flow rule
F(3) = 2*mu*( (2/3*depsrr_drho - 1/3*depsqq_drho)*(S_qq+S_qq_1)/2
...
-(2/3*depsqq_drho - 1/3*depsrr_drho)*(S_rr+S_rr_1)/2 ) ...
- dSrr_drho*(S_qq+S_qq_1)/2 + dSqq_drho*(S_rr+S_rr_1)/2 ;

2 of 2

i
!

Y'rl

til3w

(LanL

J1-fi:^*tL

plo,n-ts*rai^"

tln'e-x-zsl;.,ol o-tr,-%ta f*

To '\A*ou*
l

C-'

hwtt- h, ^^SwmL t|rg- f {a*c - fgN

p,w(/&{n.s, uJe Sho"LL

nwt-a-{

('veSh^LL
e/a*t+c5'h'6*Y}
fu Jry{;t^Fg. T/',nzr^-e.a'nA
Tfrrrv
Th/> A

4,q.

4"ht+iu'

^tl,!r,
o,rWc S+T?-,h4 r'-s(o'rg.
1ll
Joa.( ryrfloXr'rn,aJ-l'dr,,t-

irv rvreJal {=ur;Y, .

YW*

Yeq"ffi,
o.jb''trtt&

oK'6

'.fd
Np, rnJtt ir*-,"od^r.e{h.Q- sUq" Ly,"t-_yffid_

oonlttwttt Jo&.,tfrtrrn,4

t,

b +Lvec\w,o+ prt${u,J.,'r
?us +b t",ys4 CrrugtrarnfrZ
, {,+g thn^it ou^ SC,upe6 +lv*
'
[ cto//]t^{- ) SWp hh,e-solar.:fi\rrlzJ
"e+cArhlT^
+ exrhry
I
6t

.-

E? 1,1ru{rlTn,,\^ Pta,na Sd-aa,ro PIo*m*- riq,t{ ^oot L

L

--._

flcrruS#crr,rn ,

.'
pta*t -rry;c{

i8z*- o

Ll*-- A,

inPwet'
___i_______)
tfi,

4:o,

t
;l7v'Dya{;b)M
(/
t

/-+

t t .

non-711w):

I tix

| 6**

14y
\ t"j

16yy
I rry
lr* - *( 6r*+6yy
)

I{"

lry

( ttx= Vr,x

t t y y =u y , y

(,,,,n^Q I V : i(WJ* ry,,)

t
'---'
INU)'
t

(phzrc) ttYa";.
ra{Q-

|\-/
a,vht
I^ ffi!,1e{a

6r.x,* t 6y*,y -

0)

c)

(rl

6 * r . , x + 6 y y , Y: ' c

rd) w"d p(a,t.vlrc

{VTEryI
e4
, Sfrrtv),fvld sor,l,t
t'-'
( atWurna &4rw brd-J
f** ;

)'wrQ
Pd

|I

Pt*o Strzr,tr:

\AES+?-

In d'L ylutw'c-rQ71il', Akt Vnd ontC^-l+n,mrrrrtbe Sr,'t$f',Ud

c1
rlrsa.6it,: +Q^rt6,1y)- +(6rt6

k"f

l:

q.

tfl

e - !,0 rr..t):r@r, E):6?_

. ' .S , : 6 l - F - . * G - g )
-lz:

Cf*-f,:,O

53: Cts-F =-L(q-ar)

_I
J L _
,

(sl+),'+J,-)=F (u,*a:\' : ,Q-

zk:"

W;etJcrnchlffn: f , r - 6 s :
U

H+tn"+ t/rt- r,trt /tYlrr)t/ WUll

arnd^'tltYt -Cn\aC-e{tA hn4l'L
"\t-

TrUc^ Wtd cnrd^fir"r

, l^.p-r,
i^ c'>-'yVt>:brh6-

tuy'

6xx,, o,/y,

f''^ kn'r,tA
4

o*!

Bgcan,* th.L da/rtlc ,ffu

d'\r- ?Vd
=?

k-

o)

is r;g,ol, a/'l

t
ilnurhZ slrain ta,Mz,, t/*o['

N{a,^;'^rarb-s

-tu^ f14e {|A-n

N,,tg
^ry

t**
tyr:
t

,;,
orrrdi1i.'irn

I
-

S,,t

try= A S"J

( tina ./+*

h.e/ze-

'',^l&.
ru

b hr)

wewL

o t

K0/nL

1,1

: a

L otl

\tt'il,zn

-i: t@63): *F*rcyy)

: fxr( - e -_
J3g1g
+ (F,* 6W)
Sl, = 6Yv- d - *('6',y->en]
SA =' 6o/
,"Yr,tflvt.tu)
f t^" rr^fu- Ca/,\ be-

tu*

t**-i
;
4Y: X

6i.-64
z

f 6g

4:

tri
Cnnbttod tan-{hilAL }ghJibvr beztrn/ee,tn
.)

orn{. Un', h/e hu^N

\Jy,r-:
i ry
Ul v' I, v :

' '

V,y!W,x- i
z

6w-6*n
L

6 --/x v

Add^*J tt'\a {'YH- tu'ro etl}.o*fl5lpC)T$-/ uug tgF

t)ryxI Vy,f - o
S^btradu thr- ft
th,x - Uy,! :

tottdr{tli.q
r'Ac-nrrn1rezlr'ar/zrJ-v'2

@l

$N-0equulq.rlAd,'^e gtf
t Ctr4..--ory )

ZP,WdAt b1 -tht- tft^vd e7 wt'.|iirn, ilg eAn",ne',te t

Vx,x-Vy,l

6xx-y

ffi:
Cntt't"A?ll- 5

trnYrvavJ16 ,

dr-U,A Afrg* 5

gpuc*-6M.,

_-(s)

,,n;f
tTe:.L,
4I,

noy,

4c,

")

Ct), cL)

f*lilott^,ty

cJ)
&t

7'eH und'^{-"w
lnc*rtfFMt;lthf

6)

fl*

rula

MEV+Z

/+-

&i

Plo,r"q-Sfntl2n

9z lW @-frvt#
Ftr fl@ 5 wrtlerwJna- ,f th* boudarty- andN-,Tn/ ^/&- Jwjf'a'o*,rt
h /otzwr\^^i- tk* -s#e.)d
, fb P^4&"- i4 cat/ed
11,cld
Aof**t*nrtn
,S&'c%,
(utrl ewie'w 6 sol,+*), atno/ fAL

,e%

t'ALirpts fvtds /,rr,tLee,r"6r""h{,

f,VtJsCd4-bo d-fai,uedfu

ft/tprwisB. t/\r- pr"fik-rt il covlk/

o-.,01
st""hi*l.tuy ,\d,e-fezr.'na,f-e

sWA fwtdsad v*J^c"Wd

lntcL:rrurvfbe *bfrd W{t14}
No W;[( {"r-"ra rTt, sn*i ^}f

cL*?:zrn'ncQ p$terr*

/''"tru

thr- s{nu slate

ln tlw pbrtic Wtu,
nWf l;e nrL /re Moltr's o;rcL-t-wrtL

cx,qL+v

tr//

-6

ra&;,,t s k @o,ytwurn
shtn*)

Lel 1f be ttL ry

-l__J - oaz-

rn)U\r!v,kTl++

.t

+ rolv'tfitnfr) rAz

shco*,
rno+Y'tnrum
h' 7

t'r/tult>alt ull,rLil- 51ruy3% nf*J

6xx:

Zf
*:, l.,v-,

6 t t?J^ln>/
k c,$i45

lovY:
("*y:

^/rL

c and f

f..^*w^t cx,,
l)

7/,e. W,rd a'nl'*-rn,1(.D, nl a,*, ^r*rsl$;^/ti1fud.

T/r& efyib'fury

qd"Varv

hrv,+be o^n)lk't* lrt l-?r,nM

* t^'"h)=e
l#rk(*s'n'f
-o
+;+-r*'tf)
k'(,**+
I T.

u
Tft'a'q*rYu{uHr * h^J'}!"
d'9ry0'

7 a/"d #

t
I
I

'

p{W- Skalin

lv\E3+2

hW;dtu,
r

-l

. -4.*. -*

/\

t'c\l'

*_.

^.-

a at'nv;l^wfu6ra ,

,{^.!J* x+Lwg,sv,J
0

f"r.
{",*;,{t"y-+ c,u-},Ju&5 "9p CovolftttWp(anfi:-Ytd'irst^
(Jo, -k') c-arnbe-an*ilu*-J

0,A a A^""IdAnlg,^-.

furLdA'Y1efl.
Wk/r,+.

i
TlrL +""ryet* t ar'-attt'tt W v Ae,,Wf/'a)
]cx'yl
ffi)a,*t,v*h[vL
I
HenrZb6fh {-wrttu a,ndp- auuru a,rz"srheart
daruiirn, f m%ilt*/rn
furC -drtt/r' ar)/^17>'za'v a>tilri^sfuA) +ho
f*rr'
f
tryA (
po.n:h'c*..tr*,[42
bnltru, orzldA*b; (*L

I^

f "k (

o - z k _ )Q: : o

vI
I

l&[a+zkf):.,

IA ot1yJturDT?U )

'kP

{
(a+rLf
( rn,o

o*LL
lpr*

61'sta"i? 4
: 2 o ' . v s t v n , '4t
/npethus *

D{-(;rw
P-j;*u

c/s,r*,*e.rgw cArll\el{ Poe )

t/\,(-sl4^lhla* t%r. N- ^d f- ciutttvt|, a,nd tL-e.

ry txy Rtuv1rT
(A, f) wJ,u,utnY stm'L potlnl-.,rY Ta e6','1t h fud
orl){'"ydry (he/watt,r0A(f , p ) ,tvJ$,!,tor d^L en*t*> tW
TftL cl*lLe,y

i: h

a>rutrt'nt;tt/u-L Urwa in t'ltL plU't;C fY*

n/
.r
\
Ceu'
lLoW Jt1zan
frrWW
93- 6esrr,e:*
c#- lrya.,*;u:f SW- h'K

[,

l,/

(o

)(r

tt
Y/

\ \-

'

S
-C

t \ 4

./A

./

,4

tl4-a,
SM Slatt (6dx.,6yy,6"f
ry .nreF-nrys,1
wY-knro^) +/w hcaL
(X,gl,
af a f,o,*'"F
V?i4v*ivru tf tte sfaT6Y1'04' becatt-

t--_- L

Ae q- /,;n-c-n ursf mnl<u *fu

fhl

x-d,,>,:iA,

-6/ted th-L &- .&'i,,s-S*fiifu

Thr> mz.au+'ui
d

P ,,,,rcL

a tf ' , ,

fhL c-".nh'hEh
_

&-(

s"rt;1f,,{ the C{ndf#fnra,rrol4t-Q- F-At*

AU
F o4:
,col-si,, ,
t
#=btSf
A

Jf

(cnna+e'nt)

k(s+I)

+(/)--t)
?*W

t.

X--&rru ,,tlTl,,,Lclv^Xe
+ oustrylw*n'wt"f-ar" tz irtd=n-1zt'nde'vf
'f

lt*u
|^1aA^,+Wl with- ourrrt^'4tgrv
P
--A;"'r*
thlcl,w:",^&
th!
f
d.a,
+

!1,: 9,

rAfi+L

Cw'
(t.CI'd{"'*) i5.i}'p3td
T+a\ se/rrn*un
"Jtryf,;w-p-"
d,atWW,t^"uns

tl4lnf, 0t j, !., (Txx,6,T,my a^t oov#rv&'Xn6 t?'

wnd a.lj a lu^,ttzE
c,,r,tf
,+l- l,r11th! :w,'a f3 l"y.u *4-e
slrarqJf
and, o.U.Lk ifu -Jwmt lt ^1,t44
-

4L (^vf,

.IraIft!,lAY

arrlo\^hYv) A,rl.-ro,i,tva
t'*n6*S
tW
lir!enrro'evu"rf tl^t'(J-l"W o'f
b"( il\"s-do'stomcs
tm^relsd

,J*y,

1\
{\

&=-l

1,R?

\ a
d04

l 1

a R'<

aJp

Ao
t1t
t\
,'ll

*4 e,yn,ryfutb,; tllwftaa fu/"

l\,

ao,d
Pnpenluyz
Prytuf3,

, t l l \r
t

t l

A
Y'qfrffiThis ''fWV?t-slrwwurzsTonda(t*+-'0b-1"*)
uylu

(,

Carcv'{-prr

arct

J*rfu."

\
'
',

Mt342

P{ant ,11'1ajrn

,
.. *

,
---.*.*-.1-*-*

.a

( A n
\-N-

lttw,*^atu frn,6W, d^!"vu,hahV

fu sahcoJJ*y
{*fu
Vn(t,5),

Ar,rc
fuu^'|,ofrl th"t,*t"i ,L'rwfietdt

vyttr,I)

ar^/, s'ftu'tc trelU haNt be,un o$ta.,>t+t

Vr,x -Vy,J

-ZkS,n$
G^ -Qr =:
==
zku1rl
@

v*,y* vy,,
l,kna

/n'eJ^6

s't'rci)

( (v*,yt vy,r)lw{rfo,yl -r (V,,* vy,y): o

r

/,

t ,

f -

{
I

^' ' (5s

I V,,?Lt Uy,,l : c)
th^> tp*n
; ll

''nsortrTrcuti
h'ffi

4 qto#ffiA r's (,'aq&.t),f \f*bouict1,1,

uAat a ciuri;tr cs

c-ofrr
c.:de u/i fh +l'a- t LP t\rul .

,,tr/,8
f;^**^w vo,
+

u^d

d-!;t z

(V*:

rri- y
N tl,,L Ot^TvTtt
"I
ry

f h*

on, d

rufl'nvehJ

V^caf - UF9rrnf

l y =l* rc"f n vFc'ttf

Shvm StzJual"ry tl^f A't"+sr,c
not" +14.s
ck I
S'h-u,rt,
Pu\nJ-

pt^a lYd""J{'Ft{L t+.ta,5(6)

thrat^ 'atrg "l^i 4 silp hfrr,

So ilr* l.-10-y*"{J i J

D) tiwc C i qww eLa.mcsi-ai-' )

v

fdt,e. ,f

yubd-iVe ,z-Lonqw|irn abry

t kP Lvuoaa'rc Z4/1e

7
ca,nbe "y$lfel

ea

v+J'outrAn
, h u'r |^3q Ct\3oS
lhsAt-a,nl*6*fn*/;XWt cnd^h'wr6
f* fuI

IUP
.v,^'t
v

Tr \tl+
VF

{-rt

To
^

(,d"uo

tt)e
lorw,.ra''ry#' ng(a#xvY:,
-lnwL'- c;wtdi nafu tfw
x'-V'

5o tlv* t{^r- ft'-o"X i) I

-ta +14{ or-f' ruJa,nP+ar6

o
ofrpo,v"C*

.+

ht*rd*
porhf 0,
I,n I*LLnaShbar
- vpfruf'
vo(ciD+'
{r",:
,r, - ud s;rt' + vpcnrf/
[
--Vo< Uy,:Vp
,
&f hf od- Porrrf,Qr, 0':o , 5othnf, 1/i,

q-bw
TW c,rn&ti&\tf +tnnelttyo.litrn"frry
raj'u4k",04
iw{xr

7vx'
AN'

)1'"

T,T, #t^ +, vpcn|,

cd>V:F ,r"fsih{,)
r

)'l'

Q':o ,, @{:1..

At porhf4

?ve'

&
-'.

:=.o

r'P
A ry-: o

ax'

dvq-vpJ:o

Rot-afg b

sn(:o

a , h n jT '

d-ycrrlr".G'r'.i*trgr**,rn Q:Q++1,

dLk- Ypdf -c

C"ry *,- l,;*

ca/^toe

v#+"

Pfun:tS'fo'.ly"n

( - -

lo

LAr

33 SirrrpwSrhenlSlufu-s

y:tuts"
ttnfiflanslwys
omi \-JAW NtL
b,othd- Un|\eL
s'tw;glrt.[,r*M,
b6il^3 -"il \ Atw cmbfu"'il (1",!,)

.r
*""t"^
av.t*
tetrtr
olryaYo,
At[ f#eyi4

'l"t-tG I

Y{trr'fv1
7 : h t 3 " + t l o, . l f = + ( ! " - - 2 r , )
(

SWlr, >lrwt s-Jr{e

(a-('w/'w!)A slraffi

AhB @lV.

nfrhel
rt6a:.
)
f*rAry-( p-In,!- twt-Q.
cvt-L

1: f"

fl/^raed

6 rLf/,wf,'rne(
W. f/,wf'rnek
so

tle,rg- o(..hv-pzan*W,

',1='1"h o.
yvw?
/ )k-rD6
" u- --,qlyl,,
{
!: c4rwf ^\f
St 'nnd," r):ccvtt't

t@,

l:

oph*-

ea-'-h21-!a^n2'f
i+,P "1-ft1
W
a'145"'-],^
p-A-rv

%'

^r-rvr+inuth.zeznti'qt-rYtu

A rryt<-aa/o,wn?r
+ W'f
h alnatA
v J

v w '

' )

SttnAU tlvuD fhh :

@try

-,F-

:4+ i =i Y'
d = zk ( 1,,_
f)

uF t,-f,e- -rh+/r5
f

iYr a. sh,ffi, -f SWb Sfrf/H:J.

' - -

P(aru
Srta,i,'n

NleW

?5 Axts'.,rynynehi,
te,{d
,Hg,lefuviltA ilw- pe*$r,aiteei
&r,rrr.f
,ry(r^H

'11A{
N&,w ( f., f"Y) , *ntr"

t"X h @n/fru* ^ sAy h^Lfreld

.tbh#i5l- ur) +lw plurtZ rVWL

(gsr<p)

n,6 axtrry
n#ff,, rrs: Q

,w

elqlvliz'

tr A

frr w*d dpo rwEPtfinTJ

Lle,vtc-u
S'fi.e/l?/, i'.?

o'p& sJr?v)
pmrr
A, arnl tr ,or'rr-

,flru.'\.tt
)A*-

,r\

'

dl'Q-ufi-6lnd

:'
:
----\ -ll4,smo%irruarn,Sh.ewr -rtr-uY dilutf\l

----_....-

-*l

*x

be ^f +5" M+LLA, a'nJa i

rru^D+

o-

r
;o\

*y'

:r
:

-fi

""

L><

\p

n ^x-L^h-ez i
bdfh
^/nA
fnM
m-l.^/na
thwt Do(
fhw
I l7/U3 7TU
:
?ruu/tf6.1-,ra'W
bZ
At
a"rt,t
j
a
P-/;,,"'A

4so L,\ntLA anl e

lI aJ&AIetlo*

,'!^IW

rrw@ -tht

i
i

JiAmLfrrrr&lt?)s;n1nftLra&,"^!-d{Ye,Mr
-!|
l,
,rl
4
- ^
rt
-/
I
sPir^L
/J ffu Loaarif
Uc-

i
i

Y^,: A eLo

ay e = + J**
arctwnI = Y : 4 5 o

-) b: L
+ - b , + 5 " = !-

12

Logarithmic Spiral
http://en.wikipedia.org/wiki/Logarithmic_spiral

ME3+L

?l*u

t3

Gr'

Slka^""
d - lnn-Ott ,

f t
-

9- 9.

a e

,(-

a g

@rryi/-e!, An a--b'rw,

t/
I

0= *r+.UA

E= 0"-I*X

- (9 -0")

P-!;w.

g: /^ +
a"

n o f eQ
, =Tfg
A(- Ag ^,)wL+rc"e'("
"ry/q{-L't*
cd?u;-t
"ht

kc*N' fr-+_
=

l7

o- :
9rttl-

6ss - 6r, :

o-2n'nL

d- Urw
il':.eLk- AO a,Inng
o

-|":o't 7k't**
(

2-R-

T'd&

L6hdrfiil,"")

6 r vl r = " ^ : * I ,
,%olr--^: --f, +LL
-f ' + k
o l.:^ : *: (cu'+ ro,{v,)r:-o:
,

F
_
U

tl--'

utYr

-1" -f lq-f LW 1'" f

: d-k-_^'T,+zh1^e

6bo =

t r + l ( : * I " + 2 k -J - z k - ! ^ &
&dd'n3frof"U Ploxti-stafu
t t!: pJ*)
-->

0=Wrl-,= -fp ,k"0^*

|

f,l*: Lk!,"*
('wpu-

ttn )t
wi6u numorcgs"t
,

Th^ai,otkp Srlrtiot ht 4w ftt;d- ploxrc fl?tt,it

V+ Yz, / +*

trrvLS+z

',-'fra-'rr
L

2l a*tz

roLulion;:
A lut
Cq,tre
t

SUp{,n-n-,
Jtvt&:
Ql. e{e^Wor*t

bt

| /7r
.- ....-\ '11'

':'.afljlturfuA U3 canL;r,;ng

" {e-"

o(tf"*>qr/ t g,fwe

S"hr*ouZ

*.1,!
%pt<;oX),
\
d
I

twl,/y{l^tt*L
r
I

- -S^il",!t
-l--'

"

Arfrlru ths' ,2ntn4nltn-Rfl-pSM/n u

-

nn U"!- S,frtz:" ponq,,U.& bo Sux$'ra .

Artrn',:,rJ!',yr"t;tlv tnz fa*"';(y ^4 c/t4ce'd6L

hl 0z:5.

p-t;rrzs
)AU-rr:-.
F-

*r-(=i

a'uLslo^;,Vl"t, t"

-*=t" (urwd-)

j - , , -2 L r p + 8 " )
3- /-ogaif

'c.:lir^t

-W,',44,,f-l1,,;
frv,'r SunJaczt

IA

7/a^.r a/?^d/iGL

t\AE
7+z

,,":v-a)n

P ! Arut-

Co.rl

'.
/
,tn,4*al"o{
32. ? I,at - - * " / _( t _;qn'd
- ,t
Ln

r't

It

lflwJr"tw gflPjy

BDt, 6FA,
s-lota,.
Act3

SInq{"s*ut s-fiJ|z, A7c, BCP

Reg;* BpE

lro^J Ll's so(n.,h-w

'\7'

( karJ.^nt", P ztt 1

lv
+x-

t<

"yy-u

&-x:2K

\/

6::trt

Regir" CDFA

6v---zA
\

fiy: o

/
\

Wx:2k

/
\

V.t--

6:-A

A L o gB i l ,f : + . 6 = = - k -

Q=-&

4'!=T
9\'.

r <--f = ,K

Rt 3ior AB C

eklry a-l;oL, G-'.,L+ : cn-^A'tflvt ,-t- 'il-or1 Bc 1i : - k-t zk.(_TI

= - k ( 7 t +t 1
6-: -k-Urt t)

6 a y =e - k t = - A ( l T + z )

5o thl.

Qx:.+lr

- -klt

l,^rr;+{oadis

P:

- ZGqJ

( fi-+-l )
p .=,.Zq R_

( a*sWr',-e L.^br; c"tt d

t
d.rr^.ta+t , tw cl'rzo^/Y

;kvt

rn AB )

lVt73*Z
Ne twrt

d..p*e*w,),.* tl^-{- ve.f,"a\-f (;eU *:1

1$s6ott-f:r*".6'ar={.;7O
c^ndnlnt,,

nu^- Upd,/:

' oluFt va d :o
f
o'nd bw-^O*V

rRe6.-urr
13f f
I

-YE

"J"rJ a-k"r-t

*
"tu"g f3 iui"*

w,d;flrn

D- i",t u\athc-pla*hz bowndtu-*

itrgcanqrr.<)o"t +''>^r.(p{]r4r1 r iqr6.i.
ur.qcA4r1,.?"r[.
J

\,zp=-r:l

'

6rL p C

[3ecs,ur.y';,
dt':

o -

o.i-[ F" h*s

(srra$v1 t

Vp=.,;> rrt e/4"1t &Pg nqit-.

r'J l,\tt-q. DF

So t{t-q *ektftb

cli}'t-rfs"vr

!q,( rzr-eqc{
t f"1 iy1+qtY*1T:l
CD b ztiar\AL*plan+ZfuwnJ-e',,6
Vp =,O

fY1 C D

Ysor^^*,- p );nM an-c-Y[3sgi's1,

Vp--.:) o"n*ha zntiv BCD tejitr^

ic +h{ ur.{o"h d,rrt^*ir'nf otb'rr} a'fc-'l.w ^ffcs

Re3*'+BC

Beca,,n>tr".-bcfA d- 1; n-ptso,,"a p-'!"'utt

*o

tW\
I

Vd

re.rnai"n

Otfrwllavr''k

o1't a - L;y'lz

f g,na,-t

c:.^:uvdv*'f

t*n

A rr"rnff rhd,ic^Ie

1l'

t'<l-01

G'os-L
"d'oU;Fy J.lrt.e'4"n

1/'.0-{"oa't*tM\-e laurvlrr.,q"aura "+ A C onrt CB

lfcrw<,vaqfu-p-ruo'p+4|v'!-u"t'+an ""J'o$ly

mu,,Mle carut*'wYd,^lAlvr(Udh t'lYtl'ry (lv'e.

'
bu:t n^tf ,uc*n"^46- *b fou^+ *Yw{*'t

Al(+4^rr,Z

YefLflfd

Th^n- mearl.f

Vy .==',u?l/r,f (nga'tM)

vP =::V,
J

an AB

., A I

Sfmw: $
P(o.nu"

Mfj{Z

f
]
,
([*,|tt
irnf,rt4ir"*= nl,=, ,p=#

=-i. -7-.

-T:.

''),/ o-'''v7'
.*--\ Jx

z l^
"'

./

.1

'r--04

&;tr-r,"y;t"4riit{
J

b[ork s
Tl+,,:am{o,nl 6jt' t"i'w"pvJnt
J

rrLove
64 rlpic{{llrj'r'%,

B"'t

+{";z Sol,,.fic:rt ;\

nof ytntqv-L .
L.

';offi'or1-,,^l*irh U.a.e*r*ixl'-L'.y
^^ a1-fiemnhve
J-{i{tP*Foio-d
g,jw-tAyan&l
6 lo#nf

= -,E 4
1,,k
',t/p ff vJ

or1 oB

dyr A C>

*,,,,,('t )o[.rrln'orn{-zqda ta tlnilN{

Sl"rv*/-'uUo,

B,ft,,6 F'A.
?il_, u

o ==_K

t:"tc"
9' -

o-= -lq 6\ BY/

6 = -R(/t-+r
) rn Bc"

ra')c"
b', AC'o,- o:-ktn-t\ )
{YY='-k (n+z I
6xx:-kn
kcnda't^ th;- sannL
$ t1;ttf io[nr.A;-o"n
l"tr*t Lo^a p := z^ h-(7t+z I
a,,.d
B*/[- +i^"r-plan'fit-?crLa-is srno.IArrrz,
+h* r,3)d blntl:a c'"tt mov,S 6'f,
{^i./iil +M spmd "+ Pnr^"Jt {'s &t"^ffivl

r-.tatr-Ttxta
NU^td dirf nhmag ",,,{ou.(dsfyzl- t.c C-u:,vg{W
in <i/s-'-'t-7'r'ltg
)rt so f1e
t /, tsm'nqt
,',+v
A , i ( s;y^t*r' s:+nos?

l-tiil

-b;m
WI_'"M,@

a-l ftu*lttS n1^r1hln

**^_".
thott,"fu /rzhrL wtt,{t'll^'llksotufiw^wfil*!

tv\E3W
9g

n/

/-r

frcvf,JL5t1,,u/tS

Ca,t'

Dou,bLe
Crack;
Tu.zrury tlLs f/^* inds,",lt*
i4 /e,r-*ltrn-( wrrL3{"" 6ffr;e* t
arnd tl,'-pt"rcflwcL +h'1,i'nr/t^*-t't

by tluzm')rray'rerLk,taw+
tl,a soa,y(ti-t'dt, Lr)eN
+A!-aLo"ilL ctrark 6"/,

f*.tivYt

nl*rf ion
Md s^o1nirmft$J,,
lhz
it'ifts soin''fi?rt:
i"^ol<",, dot'{o"F-q
{ p{wa mirva:',refe'-*itn }

6yy:7k

6yy:o

6:K

p
Rrpinn
.oc

6- (fitp g

o' * ''*- -tr,

qf, (r+uk

/e\-

ffi:* 7fk
-rh*"
&:*"f Cca*dY

,?(tttz)
L k-

f:

Pri{lrn*iftw i>lltlS>
{,t;fu
Cnn\fu:tLfid
wrc}+hzf.
a,fft"X i.,,vuii*loetct iS

p n* e L " J e
U

Pn-

_ - - 8

f,ftY

r+I-

1 C

B -I
ftt* solutian

va,h/.jtrvsutr1+*e#h1
deeTnotz,te4,)

ME3+L

Plarlu

C*;

Sf-to'^;" lt

OtherSlipLineFieldlSolutions
for Notches
L. M. Ka.chanov,
Fundamentals
of il-'eTheoryof plasticity,Dover 2004

rI
J

II

Fis.107.

' -wr'llllli
*
llllllli|lll]lll
Fis.108.

V\EWz
i=*

,:l

-'

Cc)

?(Mtg- S'tv,t*n lf

l,t)edqe {nd<^-{-artirn

Q<"oa**v,7:a6)

l 7

J, I
Y

FI
r )

yla.,>1",'c
re.!vn
( kacha^orz

lF

-r,'alrtgYL
l' il, ea'l;

tI
I

h,ylx s+rr,o
((qcAaav

Frj t+8 )

,f LlcA= Y;t
pyoof

^fu^
toB&- f -+ zoF$.: 7d- tr
/_AFC:7o'- |

L BPA= tPT'
_ r) ) r 7 f = f c i l o

/DAE_a
='(-f /f

LEAI

AB:AC- I'
oB: lt

= l S o ol d- F - / A F e
+0) - (?f- Y )
=t8oo-'(qoo

=hg

plw.L:ia;l*t"*5

l,/rt;3+z

-_rfi&

' -ZH

6tf

4"so= o
dyr"5a ;- 0

6 = -k

ReVrrtfuAE

f {"nra *''v W&'h'

A
\

(c>*ff)

*-#=5:f"
* r , o n 6 . A Eo, : - k

\
o\\_
-1)r
Y

5o

4 l , a n 6 A DA
, ,*-ktekGS)
=: _k(ite-f"t)

Reg;ott
' S l " ^B[A D"

ei--kb-s+()
7i=6-k-

6trr
\

6loy,:6+k

o/u/2"

\*=

*Lkc>

4 , lN4Jwrt- ozr4g : f

Let . Len6+h,+ll7

('

2kto'rt)

irt&eza J-o=r {"^"

Toto.l,

F
tle Sht( tWd, t, ote*-au.*-rza-& i7"Vfu, /
fft?

aNafoFh) = 0"4e^(AFC1

vah*ryL Awlwuolnv,

a/re,^(c1T;:tt@,t--N y eG
a* t : Z @ e S , G c )

trr\

AF
'r -,._

:,r.(Qno+gt f;n{l/-85
trT

-!->-

ca'e

Atr

sal0l.gl

f'y.rqo"-v)
,{*

ut

{,/tt3+Z
LrL
l'

Aq:

(]

C;*,

Flant- S*raunI"
.! a>$

t"t} {

1s,>r(Y-9)

crYa.(
AFc ) --

!cn9

-{ J,rr(y-8)

!1tnt(t/-o)a>:o

(ey

o2
AF = (- s',uuy

co: y

c6s

l i I , - ,85nG0)

FB: AB_AF -

dny

cF = fE'[*t/-

. { ( p s \ a( rrl ;y, q))s n /

)""'"

- oB: FB-c,ty-:-/ (ql5,"U-S))a

-uL
/
Wiu-.r,h(/-(t))
ADf

<

.{.' -6y
Go).f Jn(tr'o)lct',f

!.' snrl-ol a>c9

z hatr

z c.t>/

<

>R

(a,11

J.>

t.t

5th ( /-0) co)B _

(*,,y- Jtvlu-U)'

s,h/

(zy : 0 r a6i 'b.yLtf-*l

\

f q(Aar',N, p.2y7 /

lrl'l^lrl""h,

2.O

co:&
)- (cor{3.r- 5in$-t l)i Z
f : inL;ns( Sratg-t1.x

.xgfti(l)l7
r t

't,,

,g,),

3 : I o : oo 2 :l . s 1 ;
1 = [ * l ^ r c ( @ c{xut, l t ) , 5 ) i

f (r^d't
Jo'

/ o

00

If

-f"*,
ai-,Ydaadw
It?/h5uva.
P

L
-2&:

0+I

+-1^,LhAba;\g frl^r. n
tr

er:

{(9+l)s,n/

Mu3+L

C*,,

V{o,w Sir,r-n'*

{o

a(v,-,o"Lz* pa*avn,

mifrnut d'rna;-'
_>
OBD
O,EC

ODEd

-)

fr^J

dnrta,rrtt

0+A,

Wsh,ept
AEC

t{rt,fe sh-p*,ro t A

A'DgA

u"jrw

dtlnt',w$r{l

Ve3*z

'

PtarutS*?'*o:f

$'5

iy

tt

'

La4

f l6r)
C Pro3w8-i-tod5t,

I
t
/r\
J.gJ

/al6\
l \

,
t

'
I

al.t fU6"-r"*0, 0, 9. @
\

9
r I

arrc in u,vnfm sJrgry5

\
|

, '

A , n , t . !J J r

lr-nw7ft,08,0O,oP

,-gw

Qen+,,*ac^1'. rehkuYrt,

Y=T-* ft* p**r)

Y'=Tr ;{,
Y+Y'=T

\}*l''^'' /

: . L A o p= l B o C =$
Ac,rutl hnZ ,f JXtUr*'*r&4-, e'9- oB
cnrvlirtuow
s/''/{
D
/ " { 'lrwnt 1'bL
6.yrrt'
u
/,
I
(fyyt ca,n k- /4c.o,v4irww4, / so F -o7,-b"i'e
J
''

v6ttf
lrr.
^l^k
a
4t

-1O

^
io

o,f rardtu ,k .
J

a
-\
_*:

'

o
o(

tV
d,

617;"9

\-)

o*7*,

(
I

[,r

oy"Y,
-r$
=v \/t,l,t

n
t
t
rt-t',
-t-l
4**z k be on MoA,aatoLut,

gO

the Jo+,< hnt ,

:,'
t>
f

i
ol':r'tti+*^rv
I
I

ri''t
-

'

-k
_ u?r,,__,6rtry,r:
-k

R-ry

"x'{
|
i
I

^
' d' drotnRIYL
d
l h, mx.Nl-trl,Lm
A
\
. sh%'/
\
r rt^wffiwt
t
t
i
T/It- A a^qls.
a{\
'T,trrrt+
I
rl
,
,
fi
\^
J }
!

0u , 0'a

*a
re gfact rh^Iff*qu
t'

it' l',fuo,,rPo$h'u
I
t

kerrroe- f.h-e

?(*,Y S*a'r':

lAt3(L

'I"

IL

- -LCl:-i

lr*f 4 6ar""d'iL6!,re(a,fint".

r
I

Jilf

Wr,)n(yrr)+$=x
n

'iI

zy+d=f

Y'/

3\

Y=+-+.

t,-t

(g^,)r[{ -v'lrt .d =T t ,r

o t
iv

1',

->r"t 6+ = u
t
y,,=[+ E
4r

,1frnt

AB

!-t'

^A

U 6 = f f i , . U L = r gC r? , 7

?=tu(T-v't=W,=#

+ {*=t##.=r#,:]#
I Q, o0h4/4

\ ffi''f#

soc
eV+vn@
= -Zk
fl,rt,
\

6Oc^,=O

6{ y, =o

o --*4
-,'\
K22-{a-n tU
"JP
nJe'
't\
I

Mrtl "onctuz6ty) = T+ S
c;!ty,r= ' k- k s,>6
k asE

Altgr.:-

6i,y" --k-ksa":.6
6xo?', :

yof6liurlang.tt -f*r" !' t" x" ot

,Qc",zr6

F = -k-zbs,,nE

Qi-:6+k

--2Psrt"f

6yv:

: /&. i r+J,n
I "t
6x:c: 6-+ lr,. : 2t:s;5
6jy = 7- H = -zk: ( l-9,h5s - -7-

'l;lfJ
?-_.3!-_!4r

7l

tvlE3p

j Grr'

L'..rtif Anoriysr
.J

rnod-uLCal^be trard t^
rrjid -plaa:n*c

Tlu- sot^r.no-rr"
h tks
C,oM+n^,rt
Yf/

nrnd t^ww boi^m{4

f*

{,{04+ir-

rq?.rt?ri0.(,feee-{^;19 a
l:{a,>hc_

PWftc,ttj

5l
-ko*rN.'tq^r,i
TJ ({)

inprra,tp4 mnnofor,j C^\

arvtatiptro,JJ nr,r.l,/-.''pl,^
,
c^.i.W ornlta,,',.d:

A Sfu&. {

hlitf iiQr^*\{al4 be rtq.u,wc\, cnl.-q^ap(orh\ rtrat-

inc.{eaM$
5,1, I

@tMnq

C/ri},.col5+aftt.

(fua3vv-F(odge.$33,P.1t3)

ASlr-r..c.- "f*jr-*

t\L had rru q

trrrh {r}rre bJ

,fugrld^tf p(orf,tif

"^

a/,thor,nt
rhc,,rep,ag

f r.e. }zirreafyivted
pl,r,mtft"*).

Ya,u" af ,i^yetd^,,'g
p(^n+'tf t{,", .
4 stn,aa

Sj :g

alw1h,h^,-e

thrb ca^ be Ptb.Jeduz!;nq Drcukrzr ineq(ol,ifY Pos+,{h&

'

('b, be d.sc^weJ(rra.'v>

6t' &,f 2"

J

i;l**ve'l^'d
Tke- Sjo.trl.wUnt

f*

botfuup",r-Aw"d'r.'n3 a,,rJ ?t4",h\

6i,':o

iS rund-rd b

ps-le fi{

.la+-c chtttt

+/.r0 tkln-ur\a b&Lvrll.

Looe^r Bo, ''d Tt',4!.-arrr

A*j s+*L'c^tly adr..s:rblZ (s4rw

a c.
f.rptJ) -toL^fibr4 nr5,,,r-o,k
i%"t pt"th'"tf,2
Le'.)t/y baw,nd "f filw (,oqd od-.inc,r9
sxbl-s s.|nA fl rtJ r, i s ^t *4,
A -srur"t^Q ^d,r."i
0 eqwi,hbriw" c,cnd^*+.-rneJvvlw{-RlLc .
@ bom'd^n1

cxn&i{.-on,

-"f tu
tgu
/osru,*-n\

A5 hr. :

fl.tL {"thrdnq c,snd^-+,-u.rra

..

gr,j^,,.t

Fj

: o

T).

of B,C.:tf 4
T1o1,--.{.ed)
aN\J S,^ w#,t/LL
t,t; :9

.t th'e(dfneq*o.l^\ tAe"jt"'4'rnr.

(o.^ _ 6yyl+ +oi, _+k

o'n 5r is tpea$'ed uq to a.n o.ra*U

S*qc*- T;Cr.l
u>nztv.rvf , th-u larrqezt
J
t
I
J-'.',
n - ,
. |
^
,,
coY1a."k f-n
^ :+-o*'c*\
ad-rnrp>-;'bLzSf-,*ry> ft ld Cj
c*n" s].;ll
",e^e
i5 ;+r'(( a L&^,Q/r t'"-"J -fo #ul conS+or,* + ,fd
be{--"4
pbuLc ft,u"r

vfi+z
t l

WW
t f

i " -*l
ulrv\r
n

i:

f
"*
(n
h.h,l"v,twu)
. - - - W
a "

* l

boutytd tf.!fi-Lw

ry
^fw

lcircno"h*tb ^C'a,isi,;bl&
fujnqfrtd

AA
I So{p.hort/rrkl,ra

p(affi
b**"td ,f +lu-htado,t thr"'ltlvnt

A !4Ln4At;ceUA/ a**W;uk

*fr S"+nfnt

"ehrt'l;A fyT

'
A, f-U^try os,\U;|fibtuS
O i/to{r"fTvb;Mh
ry"t{"t""e,

,4fr--sr

,,,1,,n ui,l = o

rl t/'t"V*i"f'

l^^* fr"b- s^,bre.lrvna

(^tr^r-4
ry^fuJ r"V/u'n-,,t,"J"*^-,,^yt
,u!""%

@ O^,^-ru

Ctvf" b-a cla)annnv,r"tl.

ar?d{+}lt?
,

Vl: o

pdgtl+ve
N&k\-^&,

T/LL ovz-aU

tuyfill

"",h/^'c.GJ11,Q

Nt g^ tuL'.-e- tr * e
urrcxtJS > o

!.(zr2

[,

t.o"tzft-vf 11

a o61b)nel /^-

TA)

t/,a C{yd^-.+ib.'( t, A'{vt ^ ^,flM bts^d ;

.

N=

r f _

^tJ: i
[,'*'ol

:lrl (^t;;
l; dv+
J

J Tu,&,4s
sr

L'N

\'h"ff*.'d
,*,n
i
L*,1

if :1 ("{.c,)

U^.:c.',+t^^'7
,N}"rrb
Y;fr p*r'J"aAo' P{rt;lrt-Lscona+io ( (-*!^"rr-(t\ra7ftaf ptwr*'"ft,r-)
74,e-^*rrr- equa{it* g|'ez .51,.t* (LQ- uouk d/v'!- bt @ Ofler-*{h^d m,uvf bz s^.fuv^* 6 ,b^!nnue-/<9 rhfut-rr^r' d-[7t;f^th6>L
+AL {^,r.t-,t-y ga.a,ztt)o
+
Mtb

fh^*

So()^hv,t to tla
4 +r^4..e-

iaaVft",'t pldnhztF,

e{a,ttty -/arnz-

noleyf?^l dl-

4d.l N sff?t)t ft'/d rLa* a :n*ut.f/,yy

a' ,)erlno*,y r*td

od-r,n.e*rb@

;r(

nCrn&S;bLi-,

S;rr-+l*"+o,n17

+l'^* A

khum;h'a|L9

ME3+L
9z

.l,.'*;i-Ano..!ur,,>

,^rm/t
of
Prr'"c^
Il{-e- virl*..d

C*,

( hqgw - tio*e",
.[sz p.b? )

drrrurad-r'or",f t{,r!, /^IuW aunl,ryyw bow.,df^ts,oarr4,V**n -tL&

v'Nr]K,^+"'uh.h jtaf,edtLa're( in zD plwru jh.?*^)/
f)yltclni- ,l vV{r,^'n[

Ib

ute-a+e nnt t,,ry

<^.ey^@

ft fu*",e-+lre hwV,y and atrty+rbunl

{Lwt-p.m'1.

6Yr,q&-,- 6, Stri/A$^eH 6xX, 6yl, Ary

^/t4 ^, nLhuh, 6.eU \/x. rl
t/t fe?"v" R {ho* o,tts- ccn*iruwota a,nd l,orrtd,,A; n;u-t^ /-W asa*vova+nJ

Tlp sfrvysfiu ld

a..nd thy

cJ'ne-r :!9e,"fu-*?

"tJa'",W,fUtl ca,t bL

art^n

"f ea*

k4hn4- " tyac,t*m f"-t N, bn*rd.*^n-l- = 6r- o^,

"
Sb"i4 rdfiL c,, : tCv;; + u1,;)
PrinqyLl-"f

vir+n^l' '^tu".,k..

f , * 6 * 2 ' ^ . - rd v l t y v *2 t r r y & yd) {

SrGV.rt Ty"t) dS
q4 tonj

o-)

su,h}f'ea elutAby}r'*

C-.cn*'d*zt, K

MtW

tt"r
F

Ivr.b

\,-o

K1 \ Rp

<",rr1^i-rv.-

tAa+ oL;-r,i,a,4' ,r^rl-',,^ {-t';C'

R,,, Rr,, -"

crn$a1.rxar)n R

, Rn.

4^>\-t*J
,1,$c^-,-tn*1*O
i^'"1* "J staYyl-w) ,

fi^'.)^tAct

Ok2, fVy,

AJL+ C'i>\$'^a.o't^.t a.-vrd Awv<- c..o-trlfa,ryLJvr2

itr aae'L

t*

O-xJ

clfuv+",Att+r+zl

Stt'href'a'vn.

Wo ^nr. ^-I[r"u +o,Vr-.av-L fuub

"tf !^j\t- + {i\,w"4,pu^9";.

b b+ /4ecnn'axary

PnrrqFsof
tr,,
rll:
r*y)lA..I (eqtTtvy)lS+Zl7,ro,(rl_vf>d(tt
*^ Lxt 6yyr^r1tzt
!fR
^
D
,
L
h
k
*> I,cng a-\

-)

S^**f,Vl

eqW'0;bWwu c,{r\dj+tuh rn R

lchRr v\ Shs"r ;h+2.''fro^urnyttdo.crrvra Ltp- fr* na to Rt

Cr+- V/) i.s thz dDcon*irw^? * +"ye^*,'a..luu{*crl1-o,<.r-!piLL,q

r. \' A
L/frht Attwi,y+r*

h/1tr?/n
t:iL F-l-

bJ
f

L^/

NofchedDar in B"l:!:t
Fnd M a*' rnf er<i;agplw>*'c+-6 i

lqux"t bo\" d g5fi'm'f,*

^. s+-* c"-Q-&X.
alr,r,adl*
Sta"*> ffutri

o * ^ =- z A

*rr< a

6ir. = zR

o< x < +

= o

fx'

-Q(xco

lil, 6il'tJ- s't'rUsosryon;,'d:*ur,\

ta'i-oli)'++r"/ -4R'< o)

l v =1 :z ' z k + . * . z v
U?pw

bou.rnde.a+r-rnof-a.

C-o-tr,l+wtcf,a. k rqlmrJr rr.lJy

/ "O^-Xb[

a+y flerd

lsgu/yttL

*:

/
a " m p e[4aq'5rn

)A #r\r{.

vlg id se3nd/\rh rotote

o,u-r-w,n{

a r\^^9e
fabuf

e x t u . n ^ l - u t t e j < ,W : 2 M

t"

vn*;o-te?'fu irl'y^ a:l--'^1e^frr.c

'(rsinq=e

laql:rw2#n

4\

Ll?rQof ,lncx'n+r**ry
. Z LX,.f#+

7)

d&tipai-rTn:
t"+"-0 r'A-Ftrr.n0.I*

u ==
Rffi) ," #* zb=zka'b6#

j47
Mre

DrYit Ans.1a|,.:r,t

-b

flto.

ob+o;,*th-t qh+*f rf
flr- brcynrt.
l,a/ech$s}&dx thoi rrvin\hsqa
ft

t.l8oI

of =i.1555

= l38ol ^.t c**116

lhi,t $
)'h o{

,,;l
f = int* ('x .7surxt.^r'
X o = f m l r p o v . A r( )t , i

f rxot

o<d<lr

_l
il1/- LMrxl:

M - #;.
fkt1fu,

L q,4

ktb ;+*

::L

lz.a'b :

J)h 6\

'

/.Jt karb

t/,A",,"rwutU a* fuytc^d.hfpl^+zfM,

l,oka'b < M

t\AEH-2

Lrruf ftnat,ysis

lat

_5_

frnd F ^f

irrtfr2nJinq
ilaM
J '
/

bwe.r

boun d, ez {;l.nofz
A

{rtfiA{rr^}"t

Sh%. Nd
feca.{l &-r- gr-rneAv

AB

co>ta/..Ji-vv

t- sin
I

TD = i+rnT
A B= 2 ( z a - * f

,l,e

cD: +q

r**"

-stu>b,]
f = t.,L'nc('z*srn(zxx)
{ z e r o( , o . t )
)

tF

rnovQ. art ^.

q *

t;X^J b to et<

L -

r-

\lelpdfiA 1^rr.1.,s,'-.;f';
l^\^-e-

T -

Lh ( t-e;n6)

+ a . L - 3 ( t - s , nk 3q )

F -

5-06 kq

\Lrtnetrboutnd uinbfu.
I t

r-'ontfruct

a, /<ila'rn^'h''r. 0{y. "d/'{t5"11+

a

v"hcf futJ

i
I

atyswmya. ryid bLoe.Kshd"'y

.

f-oJe- of

{;r4tr'ng.{-v$sul<.

W = f . V.F

"ff

tnry*tot '^rh*$j*T , lovrl--v

A-"t I
of Jlir"*="n"^'fu
tPrY{-h

i a' E

, f. {9{ k'tevurrr-c-Q-'
rid}9pnt'rrn,
,

D = "k V . 3 a E
r

\/:

-rA-rfu.

f , r J . " * = k . \ r . . 3 o . J T= D

| = 6 Aq
tLr-f **a-t* rnpaa&;y
frnr*z*N .

f.o6kq < F < 6 A^

I
I

MFgr

I
I

lrnu t Aneii,t1i,;

i Lq^

a- t,TWrrr ryPw pc,rordcqn''bg ofutp"i-nra

|.r

t'

yAnr,q C'f

;hp {;ne Srl^dlrn

-

C spi r *15
L0
3s.'(i-{h^,,
5+t{/6fvf,J,6rr:o

".4 f-G.
6yr=- z k-l,* *

G( fC zo.

coo= zk+ LL!^+

l 1 l v
I

wJ*"+V
f*td.

vF
It

(J"f.lg

on Cen
nscvrutk

|n-* p"?e)

ar"tr+z"4ve/rc"t%v"J+Lrl-r;X"dveg'lstu5
Th-csl?eaflu borr"da'y o6nd"7.6n)
f*
lz'",. fuq
a'rtol

Spi ra,L AB

fAY

ir,'a *V

spiua-tA'B

( \ry^'f,)^l uo{nuYU
csv bL d464h\^^^.o1^4
oLf,d'q- "Td -p/asnv b""rnW)
Luh-r1:

du6a-'l/pdf :o
d u p :* V n d P : o

^1-"?
^ry

4-L"n+
f-C"ro

fh* e^*rv*- /<laot7 PEA W, vy can bz aryulrwvfed

-f-, ,Tn
(
,4 13A/
* e.(+ya>L)
I

Ag a,td fllt9 a/t Awa ,t' dna*',w"y

TIA tLf l|,l2 rrt,rrratfVd ;

k"t't''no,hujbad*W;bl*

b e f f i -o , * r y l r ; l l r y w
fl thaLllvte,-vlebo*<nA^ft
"ilw,t 4=o
0 ( r. /, vr i d J > o
JB

ME3+L

hrrut Anat,yua

Cu

To *"f-Lu t{u d,wpadt>n

ra.ts-if"trt,rzve}na$
{-,utl.,
u$t{,
Ne M!, tl,-oprnot/* vf ,rW'fu,.r,L

w;iln t/te- sLp kl,t!*

Sf,t<,ts
/orc t^hfury

efu;\Lrot* an.l B of,

{ ( f f f / L B = 4 B A / +C B C ,

4i ^rt "f A^ravu^ry 51:AE t A,g+ cE+r/g

ota
+ lo lo'l,ldsJ
cr.)

)\trrrttyvl)
cta,,^<d
------>B'n^J6't'
ltry h f,!^tuL
o',bs-t/+

-ZVI

ru "lS
1
"

ks

tIC

f;et,t),LA\n

lrec*'t* Ai ud;tfnr ?fht&hn"/,^

= 2V ) Tv d5
G

o<f,p-np
5=1-n".
fj a tfu skp k,

DL

= { v 6ay'15
I
PE

= +v
l6oodJ
vB

: f k'v'IrI^I1'^
A

lc"

= ?kw.(z^!"21
= 3 !^Z

'kq

:, S-S;.fka,

(^1fry batnJ )

5.o6ka < tr

( nal,Jfer' LY/W brnd eh'*, l.j-/e

(6 ka )

/\4r}rt

r t

tt

,, CAi

HL.xdyvtr:4 Louut$
J

^-^*',ir*roo{.ela',,+;t--penptt
ptufir. tnndufa ,h uu{*;,,,'r.
h {at" tl.;q-},-a"ixl..
'rnuyhl&;*'fi:i" Tt/rnnrr*u^cA*ny'J, a'g. n q shx;^ ha,ndo,,rg .
\de *t-l

tX"",ro,j,^rU uon^ 6p1, *p

tn d-Al

t"

Straj-.

hw-*d*r,Vy

,f,,}o

Siq- C^'1"/ytitrr*:utl" eq,.t*tnrn,

bt Tlornc Inst*ii{*a h kzn,tnr'

(R,Hit{,p-rr)

. >k:lt- tiorulo"t6'
t
odr4
rn

- 6Pt^r"t'c
Kh-

We, wi[\ srhm.r't

th^f, ^ pryq{$

p(ar"flc rno*en,iv.-l*( no hand'-vuny;

i S u,nSta.ibL in *arnansn.t | . e. Aoflg

*P

o#*

fetdtg.

J-",-fu..t, s*far}^
SVmS

wi tl.. rwrrru-d"-u.i-v]

honfu"rUnrl

t*t;

+" bu s^J{"^e,-"!r
n-c-Ld

{..v st",bt{ft / ,'.t . {c prerrxrnx-vreuk-\ .

MoSt dAc{*L.g- vqu|-e.nnv'{zl to crdgd

m {+znz4rr'1 'v9'^'b''{^L-li^/

ne vk-ig , d/^-a--{rc vio la*irrn of Stb-hh|ft c,-'l-t ,A&

t
C r'.( - rnS*ff-'u e""f hondrr*lf ) .

-trAl

o'/rtA a'-,.d A^t*

* nt"
Larl ft. , )r h +lno- cYfl4 r+-.'A.-art^!ba^r pri'w 1;l &fuurna.tr-o-rt (i.*- w'pt *'nc* sf"/*';.
ylo"*e t L+fnlLVf A, L bS thr- v'nfruq ovt C^ivLt'Z^?t"'l't'\t^w1 )
cs nk'uvel ny'u*''r'
MelV**- +rlarrtt 9fra,^rvr
, becarvr'Q-plw*-c strc''vt't

hlg fr-awq-

A ,I == A, .1.

d G'l)

==A J,i-+ ),dA : o

#-.=-+

ME3+L

&fw,

Hara.e,,r*1
t
|

S+ra./-/L
t . ztlA'inp"pz,fnry
r t
f I
_
. +.rt-s-

n >
e U

St-at'-:

I al-

(\,>
*\/ "

JN

z )

!-L
s - --T
--

dt

et=!*.+ :,A"(t+g)

T )

, . .lf= o , : - t )
, ' ) ct = ) 1 " , f = I /

(^*,
\

L o*"

tu rr\i

f t : - -r c

, l- !.u+t)

ll

t: 'hZ /

S'{"tfbfrd Ye\lwYwlrur &w (a't^lF * F ()n &* c,rursY-,Axw

)nOr.eoant-t

lf

,atfth-

,; ?' - ' df Fr > ' o

d/,,

Pr\fi/e.

d/'A A n,,of 'tr1a2cat>e., ohe. 9,041l.-*caww tf{wl@

b6' S$"ar,h-laro,.l"z+Nvn

//J7e,( d^r- {-r.- WVr)/

'

- ---*q----A-

<-!--

'

@ -

t.-'r

Ffi e-{l!-+
F:--f:--l3ra

I ninz

/nt^XzXa tarc SP-p4;r'^d

"t'
{rht ro"{.

^ uirfua!- ,LJn*^"tm
?\,e. d-L'7<o.

S<^pfose-^ l C
ff.o-,

r,

"(futti- PNf A mdw blL

,/4@ = -{pQ > o

il-sm F,>f,

l/A*,L poft^t 4

Fr<F

^ttL^lil\evt W+v*

k tA, k-J+.

MJlebk-..

A -ootF: d (A'c-) :
F:

: 11(o{r+ *

Adtr'-1'o-dA
nl

_ A (d6-- *

otl Ld

-;

-(I-l )
d F= A ( d o - r # f u ; : A d t ( + fI+t:-

olF- A (d,r- o-dilt):

Adtt(S

dti

T =W ) = q

tr)

i^ luns ol
nX"nur;(5

, , n* * 1
+n^{ st}?,v-

l
I
,

MFI+L

St'*t{,%g6lirv^nrrU
vtTl"'su

'U'$f

S ,",

*8,fa

dl ,o
, o, r-e
dlf
t

EA

t I

A( 80t'viWJt'r6uy
L
0

Tl,D c^tnd,i&un 4^4 ^ grTh ,'cJ

CrA'

Ht-tdr,no;, ,-uMlS

rp.Pr-e'aarta-*ivn

t'r/rr4fdt-c"

-,

WA''W
,"ftg "+
\

./'.//
-4-g

- /
Ite

dtra"^, " {-*;r\

"{

alnt'ui Cd^afudt,1q;

o-* cx'|+ t )
o)
tt,5fuaTArAv.u p"trry ( t'=-r,6-=.

d*r,

Q4J/^{"d UY-{-)

6-

fTw"lkta

c(.

"Jl
SrffifU thL c&dtrhl-ii+

fh4tt'u u'/w{A

ilv\ t\t*

fhr- cd'd^"1itrn

= ue*: o-

At tt^s in{ezk,fitrt Po''lrl|,

,
t tF + +h"-!'|'ed4-stra.h

SJrai^ aarwe" rCt'.) aord

9fu/Y)g,h6"; uNv\/*- ia Lwy*n' tfu'.t

'ftL d",AAd AYq, ahta
gr"f+
tle- glv,t a gh"tl-e

UN,he-i t l6nA//v tht/*

0:-l

lJyr-

c'wv&,,
o4+* ,rlo,rtl.^zd
I{* ,pot\rft,ZJ+abq,

srat hon/i--*,.,q-6r,
o,scoJJed

(rc"te the p,nt

\

{egarrdl,Wut ^,

#;

-!'+Y,s.
tf $F trt*
*:l d"-+'-\d

O= o4I{ L

C r e c o l r (L + - - l n c r + t , l )
U'tnl"eA SMtl

rtt tfu- irtv,lr-oh]r)4"

pfiXl- lon^,te+n
+ht tl&i

[=de*

o ( :f z

#-

," t^^dL, o{ c"r*ert

tl:-

tZ qs'Pvwkt*

tt: -rc

rAESq.L

ra"p$
flandau,,"E

l'fuurtywr,{^r

yetyvttitt\ {.:r. >iaeia,'t";, r^n *aznr,i{n ?

nko,"rdl-r.;na',
4^1#.9Wt+r\*

f
(MPa)

i:-*

(-^'

=-:/!6!

e*Gv

/oo

gla,r+''*r--

f = /ott 1,4fa

E=Zo"6fa

t(a/

S't'min,4 6 '6V

6'r

Et= Y:
Ey =o

0 . x9 t i } . oo f/

^r tl,w" |c*,Lc

a Os * 6, = /<:aplfc\.
Tfu requ;e ho'ad'un,u,1rc'fu,
T h i z m p , a - , v j^ t

1'.(.

tztr iAup-M-L*\ fi"u^tt*

Tha 'urylar;nt ,U
ei7<n;aw

6>/./Cy

f-k1,

"

nal+uJ nr^/U.tr*ft/

,|',Xtn'c*n

t/r'Wr.- y,x4.**0,

{*

(a% plar+1 54la;

d'4-l4'^th s+o'nli$

sh(( be bNX/ Y?tr4tllryt+d

r,',nc/,tJ,
p(o,r4+u

bV +h*

IAr:rz
fZ

Lau;t
P{ondz*t,*
0
ttl,a&a"tr.

Jso {rrttst'c
Had**;rr.1y
-__-----tr

r/z-u 1e'*e""{* ;;{reaa.54r^*."

th-+ plua:c h>dwin'.

I /At'
vvwJ

otlrLA,e i3 n + "s*f/'A'et,'t

pa**

Ng ,n*J- &

J^e-un;k-

h fllf

th* z''-t*p- TvtJ'@

hnil

ana^',|;m
,r/rttv y{a't'l+t ok t'-;rl

"A"'ry".

is.h"P;c. Aandn';nft- -'^/-t-r!' a-,o( +fu- l+'nznt-t-L

TA-.eI

A/YV t^ft'

h,

u"ttd apTnXrw*hu"e

caznn'sJh

h a ta+rt'1'12,r"o/4

VhL c--r,wlta
,/

b<.lvwrva

rea-LrrrJt u^LLl

M
74 DeJrTr,)^whan-tU"-,Y m'ed'e''|, u'/e r"tr'fu' fh!- a7't^eftc<rnd^'{'rrn

*
f (tl
f f r,1';

P'&o)
&"oha"nsv

ia, +ha d-e,utorzfnx-cp*-t 4 ihL st\<,v!

f Si:
T'
s*eot 6 )
, (s,nrt-- ,x..,v(ds),w^td rwt- ,t1u,rrd.cr-y.@s*L

or--af1.,')

Nf^-r/L!- ,y

t >o

-9urzue- mt o'+.wLz-

c-rr'ntn,zc;-t
"r"s^^"*,J, of

,;L h

iJdfnl>tY

h o"t-oLe'^.:a3 ,

fui- w<rre .1

plartt;c ^ryrr^^,+;:cvL"

L: I oy d,jL
, rYr4t-A)-tttt4-, ,\

A*+t^ln

, LUny fi"1"-,rf
plo*4eu lhfte.a.^n

C)tq^*<*ftna

aCc,wrrvv.al-r;tad

L: IttqTf
Wa r*r.

fAft raa.u th.^- rfrn r\4^-:.t,l czvrNfg,rxrsru

t"' fv"a.-r/+sa f

: * s,',-t,.,.
fu^5>
r

(rL*

Fq

t,

( no ^^ o6r,t,,!) cxwAyond.+ ^)
Sy c.,,
dl c*,-a-,t y*-l*- y:tw*z+u6.
tt? 11 = k--

( cnr'wt-a'""t)

h"l,"tUttrry yrroJa-L,
4ro 4+tSl*

4te- c.vrt cA."*X-

S6e'ts-S1va>n otaav.

f (L ) to

feprcJ"ucZ

VE7+u

Ho,nrl-w*r,-ryLaa

t : k-+ h 1,.
w,al f CtL
a'".& (t*'> finl ou-f h",A 1i;hA
h

l*.r {cs,.it:tSjt,j

k o'-t ora.mfat

rtnnoX;o-! det"t;r*

tOvt*+b

Lo.1

? . =J

o^r,

ri

o yl l j -

S'h-eA'-5+rv,rv,h c-ururL

|Ircx) ,7, 6yy- 6yy:a

S^x==qr-G :

r-)

Syy= -

f:

Ut"Y4 kk!

* O^,

)'tt : -t6yy

tt*",

l(^tt = t tyrJ,=+ ( to**ltrn}+to*7)= *d*

'rr^l-(,
flna+rcfkN
',il.Pt
-

*')
^

3
sL
) fi,tik
dr = i; ! a,i*a.r

c.fL
1"- \ 6",.d
c4i
l
)

d,l ==

* tl *j

dtl = h: 6*xd,*

d't

c
A
z ) t trJ

h-anc/'ttu,rg {cu,r/.
J

( , ,
f( I'J')= i t;b : <le)= k"+p L

*ui

- k'+ PL

L ^

i6uJ611

J o^^ =

i;

^ - . + o"j, a*
B
, z/\

F.lZ:

* 6xx4.*

olfux= 3 n

t_ 3 Hana*'-aa1t0iz-,
d
,
6
u
._
= ^|H
*dg}
'
=
? P , +fI
- =& ^ = ( r
t
rffi En g,,
J rTxt
e ft f
*r

No@ th^*

*Ajtn

4A*/at

hand.p'nir7;-i4 1tesa,vf,
wa'zin*' fvWULc,n

cr-*, be dbtttn"-tutJ
{rt
pA'f<^,4{,rfpl**mo+4+,u-"*l
.
X

&, h"d tha A-,/ua6> n,&. A ,,vuroq;>-lLr*^

f)

t't'vtJ

ta*eal'

,L Z - I,F${

f\Aev(ft

L
Lland-t-rrrtv,1.

catuLzlTt-*vt-gLt f
J

K -

Ao fiv,p1-c h a,caar^'ry &,f,'7'

rP/u?/1tr1
l^oo,d^1"4
d -

W(A ftflw,in ra/?/$tlrvtd;tur

a
otao iacne.a&.J

rw Baws..t*y
4fu

93

h'lu,

^Jq

H4"d43\^.,L_ftto"lsL
CTro,^l"nm"!- hn,Ldo.-"ttn

f cti - *i)
oii =

rrn Ka tAon al
)

= R I .L
. t''j. d

L gt^ {_lnt ott- ..ul arn,1{'.e,,nce in,,-ru oxt"r^d]'tmn''\.

Sxr: ? no,
df{t=-

\f =-\ t*", h+- - ! {^n'

} sr lt

a q ' i = o d=l- L d t * '

r#=[J*'=-t#
xipa=c t*,

d y y= - ! c f n l : ,

d*r:

-Ir_ *

tt t'-cri) = i (io-"-, *1)"*(-i o*+;u!)i)'+(- s6rn+i,at";

]
=

tfxx

d
&x
=,'T'*r

o l5d

/3^oPL\a
\,+oC-f .'xx )

na.
rr K

= 'F k + !, t{*
z.
=
Hndanirj
ratr-.P =
J ffi,

lt
\ F

3 r-f
E)

ME34z
NY*nu-

"f

C*'

llo"','rt"^,y L*IIJ

P;as,rnnth'c fLo.?-obli^f nolul-

ra fUJervtt- k,ada'al,

.lW
tBa^s"a."l^?zv

,t

,f
e.suY!t, ,6. c,^t""btn^*r'frn

i>o}ruyiu a^nd h"r,vnlotln! ha,"rrf*nrg-rlroU.rL i

cllr\ b- ok/A\td, u3

f,t;-xi): Hl+p?
x,J= . r;",
Q-,-

WWK o,r.li- t+S oa^/te-f)$^r4-

l:

I o; ,tit

r'^ un"'crrto.l

ft^ru!\n^

h,b''t

/ ,

l1rur{o*'lto1,
L^s

Me%\

ASiouvfrz,{
fU"^t1g{_

9t
-fhe,

'pL

ft-

r,4 d/n

+
t;. : i rii

fu-unttt>"'ift'
rt*[.c-tw-l''a^'fAi6OO,-,nneJ f-t""f

L^''r (a)t+\". +Vrz- u'ur, Ma'rza ,6-t'e\d gf,^thh\r,-

: 9( t )
\r*,, ==lti t)
:

0n t4rL *lfgi, h.and, it t,

nnt a9roc.r"urftilrdith th-r

frtacrr getd annd,,{i-r''/\ lq - Ot l = q rZ) ,

l
6,>.6t >.6t

-To

^ f-^er>gfti)l
'd t^tP\"rtrh m''enn,l
., We 1h"X.l intvod.r't"c-a
(plw>n ftk yh>l)
+hak ttrt- fu* ywLz-6.a'nb'- ovnlk'- ^4
s1^<.h
6t

',

g -- +_

I+

/ t

y th*

\",

h,

^derJ- lL t4 a srala'v .L'|\A/h-',v'

:

i#,

+hzn tlw f h','v rui-a A a44 a c,;a'hv{.

( w;+h thl-

^r,^ r r

1! t,utt r*"fuu

f-w ay uvwyLt-,
f*
l

'

Te

lJ s^n da-{'-}^ )

A fu^V nownJ Io yll s"-l=^* i

4d;cfhw2
' ( - t t'j S)
tlrt vv'4 MrW'l a/'eld cad+t-;urr
r

f ln-J Yw(z tz
:pL , . Ef
: \t ; t ?,
t'Tq

tl,,,L x+t o azuted

t,f

TA^z h ",i+t9*rr,t M

*,?L: *

d.. nl,rt ,) th!- alrtooro$iv(fh^

/'t.
I
1-t r J5 J; s l ; 1h) ': 5J ; ;

rl

Y^b Y

tl *

iou.,,t;t'Vf--

Trpxcol:,?Pld c-.dt-su ?

J-t-d

g 9-

, a
lA'4
"

Lr,*'ts
tlarydww;rlg

hAL2/t4

tvlL-

t {20
stu{ail
Y( ^Cl<-a/(

Peo)
t Ka.-!n6-nov

Cxw$A,e^o, harralz-13 rrrul,,r,rntL,ti*h ^!.|nd

'^it s{pb- ad{.n+tan"l.
strw:
-l-[1-0'
lno[^t3 A

{6:l'

"^d

S$rt,l+

!t,

m ,J^',,4u

tkq^ ygrnrl.{-it,

-l
shr*i e'ru^r,,3'u be c,er.ontgud is<,tharr"nAQ-

^ts ?r ae{<et^ Posh^{

{- In thr- to^Jarg

l)1t.,(p)\, crr\di{l'onA S+*. dciy.doea

P,JSIfi\^e W.-4(

t^nd^vg^& ,:Lo N{,,,g,

F". a csrn?tdr r
1-jr "f
addr\t*n"f S1"Lt1oloc,JFr+ve u'dl'k ,'f plat+tic
cd:-ftu"no.+lsn -l- l-or:, p/".1 ,

(Fa,

a-6rndlt;1c,

orJ\ f-.
o'cPl.
'tJ

p**t

rno\+e-,vft-L,t{..

hrt"/< t.ttl b* *qr}A

'a,\arHc &trzl"'AJ'irh )

(.o-n*fd"a, ho'e;)43 Fofu A + B ' + f

/'unh6d,{Yf^bh C - A

Po*Jott- (?.) ,

?rt,\*-at\

I
ttzrN

4:6i-ur:) dty >o

6lnt^ fkftl- u)-tk d-rr''L- o-n Ca"'tn sfrtu;

v'rr"lilL
tz(l

C,'\ob*:"

61,,Yry th*

ft

U"/JL

LtJ HJI'Q

fi @i- a,) d'f-o

('rk A^L

-c \ tn Prl
: o
6
f V,l j) t ti
| ,-

Br,fr f wr$,

s'l-t^c"-: ot\td- acc^^M ,|Ml"Y t:-> C

(t;

t h,t:l b4',l-n^!- W

- oj")

oorl {i

l{onW rt ,5 a Pr&,f

{t>o

(ktu'!- ryr'"xltnt$lt
fthrfre)

r'ttf;dL +h.v xiYd :

*)a
Strtn4 rEkit"?rtrvrt,

lhfi+Z

lws
hfurlwruy

C1t,t'

Tf * ,l:,'r-rl<-s-2otW,4crthatLL r:Jr^{lrPo rrfr-l, ^}a ha,vz

d*W B-tC:
da^idL,,; ) o
Pruch',,l; y; sl,nle$s t t.,
,-/

ortvnadffiYy
fl\L ,%.fu, &+Cn (;i,

ASSIC; oh ve

Fh.rN

^/\

'Drnc"Wi:
tt ;
ffi1r,fntu

du,t-dtj: >o

c-aa- Atn S1U"g'n

?L

dtu5

rts*,rw, tsled^/---

pottttf B ( 6j )

(^n rgdlz

6/t

+/tl,

afvd*I"-*

A
6;

A/oTt th-r ryVwc,t.na$

G; - oj") , d' to
f*

"fL

o,Jt in$;tr-.a
Z,

o.l,l ponta A, ft',

T l,il
tfu

Vanu-

> ra*vu/

i 5 ^!

polokrl-o rf

d t{t

r'r

t.+- w

'

N/, 4rrt
r^i0r maJ-

d Forttt A .
hb td s'r"fa u- 27,

t"JL gey +W
Dtcuc4+.^r/spos*,(o.,tU
f *..p(,nus a'*4or>r'o-f,r\rOf {"-'J ruk C3t- f )

I
t/O

l,\83+L

lln dsu"iri Lanls

| .-\

t<_

t''1 t

c4r,usx";{.,t:"t \.vA _l*,fu*

:, is u* ctnvgy,
ry flw v1*tdSw.,["rr
A

nb'w NL ca/rL r!^W

(s,;)

|.ofr\tr
A G,i)

Jt^Ll"th,f.

F,i -uJ")dtiL< o

!nt{J

?w4att

fnd o-

ur'o(r't-,{

Dt ^rA<o*' poshJafo

( p\vucnl-\a
. lh,,qv*,u +h^rt-urv/r-tA
\ r 0
u '

qxlr'a ot*A in a c,rmgtltl cyb

LtWC-. A

thovl prod"rc-r5ptaMic &3fovrnarfnn,

is TowSV'f Vw.r ;s 'hltnznc-[
rWb
u.1,r.,1..in t[,rf nno'rftunQ-)

F"-. p v-{.-rt ptw{rzlfit I nt, hoo0t-e,4.\^'T

), Dy vrku"a posf.^la,tt
i5 "\ad$!t fi .
(1"*t

d\6idf : e

^d- [,0.d^1 pr\Q-e2B'-C and fir ^\

1-t *lt

Bd''q)

Oyr*t Pla*+'uq

VlrJ'+l

Clo'

?,\uww
l\ b"{Ok n^dG-A-is ^ poUJu,rX>*,{
{xw..'sfit4 "f

,"ar\j t,qlu- *r1,rfn.|

,.1railts I o-nthr atd-on

"f l- t"l*
rif a.vru,ted O!

th J-

T"/;

\yatn

ry)

bou.r.'doa)et

Sru'3Lau-'S,sf"X gYa,v\^

a,,n
d plwr+t
cJl,'14'/\^tvlrc^1y
rrrnAsotnlD ie

[i,^.t th,L aggregate+ -olJ

gtaxr\^
isolrY-:a
CaJ3
of I ews t0 be-ho^ictt- large*

W4tn

SCa.\-Q-.

we- tUarnf, +,a "d-ptrlvP," ( a. urn&r,c+o,r,nd

) p"!o1s+^-[
G'S+--.{^rrd o ( v ra,4o-7ro^,bt-o
) aplrvo^op1il*r.-%,

If

is t

pt,*

u,rtd-av*wyto *r\e/-4-- aud.r^L

arnd,tfurrv pe"ft/*

+h4- proc-r--s
"f

fl^N+?;t\b.

"tnNvnt-p0^ryt&c^fihn/,.

Tt";a i-S as9v,zrniy th^tr fhJ- 1r"^/\ botn'r'datv"ieA

COel
d-

nrof pb!

clth.Q.rltho^n

4 *f^trh't

orn qul4rQ- {oV

r.l platt*;c- ddF.-A^]r$'4-

pva-"ndiv-,na- ca.,nlfca.jrr,f

g*ri*t

6rL plar$Z

skur').r,

( G z. Txylo.rtq38 )

Tr\^/ 6jlemf1-rorui S c,o',f;r mq{ bJ Wyagr"a.,rta^*or"ln;+-rt C-ond^*i:rn6.

'

l""{L-iz.rryra*-wv,:-GB St"dry ^/r,{ GB d"f*lian'h

bn Cotaapl
crntd oiga i g,"a{.u'"r;f$
plru+Tc dq-l-,vr,ma,h

Aeu*l

O""pd f{r":trc.Q

$A*qL

2_

f o.4

sry hgt{era
5t {nN C.,{,rtut
a,A )
ontr)ntt)-itt7 'ntiu"el,,i ( mdal-l , F-rt",l'crn;Jt-a6

mfrt

hallu cuJvic

,10il

n $rn- c-trteneri c,ulic

l/aa1{uss, uJth- ablr.tn-'

( I:Cc ) | aJf,'cz

o- bo4- c*s,funai' ur-lt'tc CBcrI

( u-H-,;'a

F cc {x*h'ca
f

/\r

[.r,], fltr ftLl.

r.ile p".,t*
(-i tl." S*^t
"n-t
\ft-)^farru
hX- ?rt 6,n FCC u{ BCC Ou$rrt"{ .

I{

cn earl- (.q,+*'',<-r
k&,

BuJ-u.rt- can aizo p"f *xli.('s- atnfl\A aJr eadL /^^hu rtr
t; c+rk d{llllll,',\.r c^yt'l- Eltrv'l'fwo
txrrzrrlrk , t& d^'o,tnnnd c,,"h'L S't-nr.rA,rmC"f #**n1
f...
.Ana 2- J- U-*^ ^,toff\s fn e rt cl' F CC I nltt'ct SiDt

I)

H<,,ttrrrlewitl fcU s-n Sinryeor-'1"i.t \*rd*r,.^-qirufth rv,iqrruo't eo-d^ -{ur.+},

ct tib"
( n^d +hr to,ryn{. ) a-:.itl^tu
W

TAEMz

G,ptn-{-?lu+rc*ty
Nhw o\rlh[,&*fr*{
rnnt
/<,n

e A

b rile,tr,hed,
r'f )

t' t6?

plaru:
"*ptuUq*fh^Y

"I"ry sye;-fc d,Wuh-m'

( tt;r^

l+tra;lz [eU

rn,ruY but fq . ilclrv-c)

Th* S{^'y&irvvhdnA o/r{ th!- Sl+w{uf lo*rnr?ftt

u{A,w5.
i
a/ft.
thu plaure.twirL W
ThL shV pbn-e't
( /a -g4/yf lak/7-f(*
rl,S)a.nu):
in-p(ant- a}^t;tV

-t-.rirr.t
V

I.o B Cc ,rqfcu.[ ,, th.-t drr4Mh/rw+sGV9'6sWtn tA

Ia l= CC "^.-q,t-^l., +h-2Cor'nnarnt I t,.p s.6ai-enn

A

frri\ Pto'r'!t\
(t
d"'w-u+rsrr
o to)

f ttoi 1>lwtt
f<rtr;

f trol(trt) f"v srrx'rf

{ f ft } ( t t o 2s t , 3f , v s i r n t

Q,

( st n {rtz1, ,{i2l} ptwrz hot^,if

^bo Sr- /saued th Tcc xr:,l"L: ")

l-to,,^lvna,v,? sLif s?<krr.

i'o FCc

wnJ

BCc

dmtttttn

t,l

1",r,tiu.s

Me*z
lr

t)

GML ['tnvt'z'\-t,

( r't'

tl

tt

Y \
'- --t

)Cfl,lilh'd tnc,rc^r

(-Driti'd'!* h u*iuX;"l-

t+rrdb {Wt

,f

o\ 9"\'{& untrtt!.
TI

Wni+ vpt"W alny 4+rWUol{tA

.R-

wtt^fvt tM t try Pl*

D :

\ J :

7_

n^Yvrnal

tV d,cyt-"+n"n

wv\4+up-rt'Nt*t
lrrn*;b st'rw)
^"\qU
0

b+t"t'wee'r.
,

-f

a^J

c ' q q : (\ {+ . n )
/

f,

b<,tWetrT ^J
a^ryf-,"oo>f

b*

: (L k_)

L*

Lt tA4- reJo{axd s htwr sfre*t

b.

n"-

7(nn-e-L "I^V b.
tA!- Sftt a"r stra*vu o-"upL^^a * *"t
k"
S .6

Tw *-

( *

c^r>/ uof , SchnsJFr.A*

(/ ^J < { )
( v

K*

V Y :

on tk-p- ve(^h',g xa'e*^.fool<-ttrt"f

,S clef.pnb

u.r t.

lke- k"d">p- Ayt j

t {^{, t

f+e//)
h be cu'e+'ta,lfri,tLfe-nrt- ale-eyc* fy 4 s;r*z
"''"t'fu
'le"t-*lsn
dufud,,r*

( b,a- ,*t

.|",n<oi ort

t/'J- ut'ytw>fx;

*#*.,r!rt;,mtffi#"_*
'lt"zot*U o{,,tci.by
Fcc(,,tcP)

VE:I+Z

)r,V+^!" Dlafi'urr;.
Lon^r
$4' Sc"&r'''.d
Thu-P- arLR-w,-A*-if!e sl"p g4ttua

(A,t'

i^ A unXlt u,,6r+r,!

L I D Ygal,mrlol,p-'tD oxyzttf tk^t

W onyLuurysnk wn&W uo,ior/,"!-slrvj. th-r shy tdslslrv
Wt{h+W -A^zXlatfiS"h,-ld f^"t*u ) acM.fed fnlWnd fUtfrattyt +hl mosf-

ac*.\,,{ -

9&wid

l'fa,rtetw, b, ce*6'/n (q{, yry*tuP

m/?-a-vt"gv,y*e-a-knrt- ,'i-**';1
fo *\

taitl be :m^Wc+tnuy%

wo,,r,fra"myil,
huOr
trl-r\'J {"'-,+-t',,
ruczfirra.te
d.

]'lettcl- tl+t *z"v(t, ^xa cnn b+ c*ugwixd

snftq sA!, {

14, ^ pW
%t^!,
A o^,ttil
t

a,l

y,["li-:6:f

Qr fryk" sA""rth^t h ?e"w*t

gram m)ux auliltxu

6t o"bQ h

La*t

-t6V,w4sUrna

6161nntdm,&- thL Skattn

'
c,urr|fyuu)r,fu

inTo*d by tt1 ru"ghbo+,s

A.. ul!,r+fiw a
-

6,L-f""Xl^./ , J.I,6i nrot"J.l

U
6 z. 3oT C/qa*-t
M^^/ k,*ttnt,
l

tq3&

r'o /V1rt-o-lyl',
o-,rCry.,ll,+e
otvx/ f"t,'to..rlLy' 7-{f-y^krw2!.

"? lr,b.li,

f rr--

I\AeXL
C -

a\,

qo-orl.^'oProig,,+,-ir^
S. S'fefe,o
b--_ts+

*o{. +, *^fl AU,,tle,t?t44

"
(cn ^ H,r'r^*
sph-ua) ctnt" a. f{^nt_
plnrnn- b, f"al'-t'ta- J\c,,ou.asr,a.

$i-g1e^qyo+ty,''cPyuie-rjFrlf^ is
J

fh" *^{2P015
ff{"1t'^eBlr^W

bu# d.^ nof

N-!^ * lpryI1r.
PTLTQn^rt-

1?
[oo11 M (t'lx{.[.Pots)

ott d*eokut/
PYoJ'otf,

W,

gttt) ,t
| -,/

rf tVrtu,,trnvtrn4,h*h.efotoJ

! -)

U/^tf- sn\r-a^-!-.
I

troot-; prq
[oor]S

IJ66f]rlobl

,V
Loql

,"---'f'
o

Second
vt-u\)u,14
pnd, ryoje

tolrl

^tuW M.'fV4
rn lLu ylanaW*tM Cnortl'*Polu
)

* j..
t

*."r /

dram
Thrr4{,ooL
k ++vpl"mu
prfer"fr'furr
(x, v, T)

cX.Y)

0ntfu Fkru ----> m lAL/arufaW-+Y

fr'fu

A^rolzz

luot

av clpts

:
f'-a

i
/

xl{1p:1,
}= ul}].-f

Poro6wf'oje"rm"
flnrt*(X,Y ) - m

ry
\11D
{r{

,#)

g6'+y') --+C(-2.)
-GE=
rr,+yr:
--tT+=jr _: +(L-P)
/\
--1;1
\ \
,L

,{'l/

z
tf q+T = tpr , K . I

RvaremqflP,hj
.
-f -

14-

\r/

' lt4:a4"
l+-**#'
+ ' +

--"=---'= - ----"...t.-

"---'"

Pofifu6\ ',w*fsflteJtL(* y,t) ,

lf/o{fh,r,,^
l-en*'ry1r*,-c,
Z)to,

*=-Z--

r*ff+f

' fo,fitnq a)cj-q- ;s th+ sma,l,LW

wl
pr1'e,uArnplarrw.
J
9pr.rV-utrrJ-,r'-- +l4t-

'EarJ*

ru-/n*w*f a planz {pa*;ry UnwA^A

O)
1*ot- uYctra
thz rvnrrt"!-,f tl,A 6rlarnt-catnbe o6*a*,"ttrby
pln'rta
+*?;yy oynnpndu* ,f- *A rt^n luhnA th'fA,,b

$at1Pl*,*

oio

6$.
i:
Q rul'rlaP
( ttr)ploru,?

i4ElU.
; QUtiq,{flry
I e sternpgT-flf,;i
Trt'**tb
tire-o/rec m':Mt+kL b*dl*
Notics'hhotf

cr*ctL15

d'vrdsd
i^t' 2+ fiituyjuCuir+h
auuued
Srfu)

UW tr'*W hu a. '(oot)nodt,a (httlnod-a,

anda<trr,
WtN 1"4\'fi4fu-clrl9J6 dU, SarnOla-eAon #* wl'rt sphtrc*
J

(il% dnit L^pfu swrra0,/..L*

ayvthL Ne.,fib"rfry)

fu-/ hnnfr,ltis eyw,wlt',tf

,Jt{,!rtbVs,6nry
h eqeJ^
rt'+t"l,"W,tCq^\'c syrr,ryfwy
1
l'lutu, frsAaWt #"v

dyn \o{wtu-,{ fqgi[e-lufs,

NL *.lM read 6 (x'arm,irtQl-urgbu^N^ wlthnv^ bna4ylL

haaq-{"tdlyrnmuy.
6oot1
n
ltot] hnaz-f"tdtyrrrr^yr]A
0
Iltll haa)-J"ldsymnrr{y
n

'nEb

( - , , ' )\
'. )hP -'

rorl

->X

TYT^r'fr,V
ThL S|ordNYd,

Ln gOrw^!-, it' th!- -tr^rsltuarx)t

Li esw A tfu, co,Uu,'

+

tlv Sl"^dirLt, hAw^ifu, flqrn

avJs
Q-, fiird-,Y ihth-t'bn'rnd^n1
m +,hvyzlviru.ble1p-

s,*'ry ,aA,""a

if +ltt '**ulsM4'ti-tJ;P oc-euna

NYth *tlr,

ejrLt tt th-e-

e,p'tuf ttu #ry

!l_i**
-

'

fuyfW
Hr"l oW a +h-a Secnnc.l

r"w*?

S"An,6J
f^,fv

tltu

UwnnuMaf,{^bcodv to a^JW 0,u lL rL"ysy"lo,"d ear.tw

ftr
^/xU rVt tfu Sfa,nd"anel"
+-0vrt6l&
W,

('tcc*sCy-ry>lwr^6f,-)

Io

n
v lrrr

cirrl[ttof

tr T ol
[t o i]
LotTl

lll

,l

c\nt)Y

(tit)fott:

0ir)[rro]

It r o]
[r o T]
[ o r r ]
It I o]

001 < r r i r r o r r l l o l
t)

D<hM,
' o

[rol]

/;X*e,r+ SJr*;t f"^,Y",

t ryfurX tJa'r ,l fuA'/
!,, -f<.c-c-nc/
J I

1
j

8", t/,l,tlrJ"
(r,rr,t"ftz ) t4 arAUofud,
st^pl,X,S+t,rrt
w/,tn 6/* knA;lL
axb tA inrdq tlqY*o^dt+oL1"'W7
&,

/1

n2al-e"'t
& , Ts f/w st"'S>
uf +t*
kcsnd{ti^wF rc*n *-fa*w
+h-o-Sa,rn-uia thy
ahra',X,5
+^^6t
S+a^d,on-A
#,lLz
tl
/ t !

_ - S r= o , 7 . l 1 7
+ -l=4.3]JJ

tl"?
&trerlirn

\ .g'?fr
'lxy

fi
N
tu

l $ l

\l

t \

T':Tfob
+

TW*

O\il/rz

I,rnWnl
{+nn-k

Te,wlp axi: rofotLe.J q ave^t*d-s6JLs|p

tow*rd"frhgd,{ve,-h\w

:
\

A-, Niil +hr- '/(rutillt- aN) lg!,,r,u

.
\

tL,L shndn/, +ttW^VU/

A"

11
Schmid Factor in Single Crystal
highest Schmid factor S1
111

0.7
0.45
0.6

0.5

0.4

112

0.4
213

0.3

0.35

0.2
0.3

0.1
001

101

0.2

0.4
x

0.6

0.8

S1 / S2 1
0.35
111

0.7

0.3

0.6

0.25

0.5

112

0.2

0.4
0.15
213

0.3

0.1
0.2
0.05
0.1
001

101

0.2

0.4
x

0.6

0.8

tA?3+2'

C,y+*L
Fac.,t-nr
{-, P,\

|q_:yhn

4..7

VI

I L

h(dl
5

q/v ,u v4
u-'"

/\

ta;
l/+

eAl

,r:

l^,t1.41 V
L/\-^*
LJI,L{>

C.orrt4d*"r th P r W+ryyto( wltL

l-

(+)
\-

^l
rl z' ln r ) - i r i f r "

::(;y pla o- n*rnrJ- /Pp+s(

f ) r

fl "
-nro,

SUP J,YE"tin"

'*--4
\
(

,!rod,,fta) ovni L*^' ,u* wni? quhy

\

(NfuTjr.nLQSu1p& tht- St^yvp+o,r/v

rfuw /^

\------

fr t!,Pp*d_ra*
l^ := sbfd'tttanrz
vd,u'nt"f o\t^!
Ther'til*(sfual)fiffi"in l,nl.,ral

tUyrftlena ,t

ro,)
,",: +({,r.!,r,o b,o,g
-/- \/

A.

r.{;radic prdw-t/

in Corr1zorutln

albo tp^roJu*t ^ rDfa#*t

Thn shlt
(Jt

(u)'^'=

od'

+lro /eng;(n ^xur aT

-l

"f

,tr

,
llodu,f

ax,\,f r "**,1",

,4*'

poi'rvh,h th-t

/{'o'/

n wnXle.
L4 the rofau,f6
cr ad')

bt*))
I\ r nn,x
/
u
\ 'Unvt

^'4Ju/nutft"t euu4 1x?*41oKyw'one- le-+,d:l*Ta4{or

Sil?;- ;t
I
a-d
T
'
,J"^,L.+^-trrA

,J-ry

r-orrfr<4!"t't tfrz,vrn

o.f tt

L(l'4

.A HA -f
v

2 - ' r u

f*-->tr

I '
I

><

6- iaLadt,f
u

cl
--

v) |
\ " n

i
l
a
4 i

J'fl *fw*r'Jr

lvun

' -

/ - J

c1n
t

\
\ -

vt)

+i &T

MetrJ)tt

-r

+ ilJ
rDN
r \ - v \

nott

ir, frrr"f*il;t-url*"r

/ ^ T ' 1 - -

f\ J l , ' tl ;

- 3 t)

G1s'{-ti F(a,r{',c^'{';

IAES+L
a

ft{

Ya/ndcy^h cLr'rur. ry-,vht, I

J
, ^ ^ +

/)

Hencz tltt

by a nw,"fu,r of
m*r'of b+ 6.CCdzrtt-n."ct.,te(1

F"{ f CC Cryd , t hetu'-ar? I 2

/2-

::: Z
,,
J

^p<t

t,i'

.(

V,=1

-l-ryLsr

)\
d,:l

sAn*d th^* in ?u,u*t

sty rytWnt.

SU,psg>+e.rna

tff'4f!'* njt't?'
)

/o

5 rudh*wt l'n's M-e

(*d s^-fi'c'o-,t
) 6 S"frM rht l,,yd

*u^^T

6".f tke fol^*t{En l;, f/ut

7""J6.

^,

V'ad

,xvA
I

tt

a'titvQ-t\uohxn

d 5h,.twld cL,oat+

lr"J

o-r\ {aofr

nnf un'f

/n 3 5r",t,c,1,,
t/4^*

z
o(

lY*l

z'

72a *wnapd
u

Ul wun *L( --v

T

/'J lL.Tw Ca,/j4/ il!-la

t Lrt jLds,, .
/

fuo""

r) 6)rvirYvi'':ecJ

F IY"l

fdgb/

rhpond f*.o,v.v.:

I*fu a-b;*0",,,6.

D atAo

b a$;{rar;{rg' ntvutfed
0

tl,\L Uup@agtywyilt a,xi.

ygta;tt\re Lo

t^,

/)

LAJ

s+von TL
rt(otfe-atl,'r crr'fr'cr,[ )'E.Jo[,ua
d sh-a-o"a

tiL'p x6svo*t

tuitl'. -fh2 ftt^

sh+ttt ,T{ ,l

Po)r"'Xsh{,

)
l rlitirr t'o,n,Lo"mr7ia-v'n o'i\q-,e,fx'*;'sri

s Ip :

M'T.

Fv Fcc c*pla.d rrith frft](to> Ihilrsyltwt, M * J o5J

6 9*-t laefur,'J)

F". BCc c-.g:t*-{.ai*h frro'1(li>sV ir}u}*,, tLu )l{-j\.I.t

r^ll..tr?
{^, TyW {-"t".. i:: icl*,"A'r^{
J

Oror+"{

Mez+z

nl

r \

-t

8p s,tvi

-fn'{t^i

l^

,
i|

l)ln(4l'rit'tty'
t l L > \ ' t v t l ,v1
'

( n,t
\_wt

t+

b"r h'nt'*
Problp,,nn
*nd tuIa,nnr.e].,Twns. ul*^!,| . ,f oc. Af rve

Cl-"

z3f, t4cm-l{o5 (/f 57 )

To^1[m's
prrbtz"rt ;t ^ slztndcuro(
y ntb-n in L'n&nt haSrw^^ n3 ftp)
h'npt oQ

ThL fA^tla.h fr/'*1n

rJ abb t"

-toht

rh!-

l.U^'V LP pni{*'^
TruirL f'' '

-iucln

t h ^ f ,| A x < b
1

7:

*
I Aq^ [e[
Lb</(ub

vecrtvs\
(*,f,b,lb.ub
rnotln'u-t) /
\ A, ry

v"L^e (t 't"l ) iq Tylm}

Th! ahsoL^nV
noMW a c'couunledf-{t
l-/we-t ut{ ttui((nof u,K

lrt LP

Xu =

xn*,.
fo th^i

= = ( t , ' )(tz.
/"
-f*

Tn^y {
== lnaX

Xx.vo

n
v

lae.y
wlnhi/ c{nah'ru"h( A x < b >

A . =[ ] , -

I_
o-

r-1
L _.Jr.

o 7

'

l,/

/,,,

d,=(,.2,

]
0 i
)

tu'

f-,'Y^L( d = / . , ) , e g

LtL: t o, o.

\{e;"^b

krf7

l
7
r
t
xzv)' Su'h tA ,vl'

X - r Jr

]'{encp-

ynbLorn Ca,,nbe

b,7 DtfnfuuY.'. 14trrvr"tAabb

,

Jo wL ni{( S*

D4* r

f kt^re"
bor,q,11 (wpry hw^c,lI

ub_ [ J.

nrYej)

&ryr+*itPlnth,t+,'y Cu

t\AE3+7

l5

6;-': *{ nlr6;' r nf'bf')

3**

0 l= + G r ' 5- t l )
tr( Sralf-t-1
tlv

c,an':4ra,v"i-

tt : 2, Yaeji"'
J

d,:l

(a7l,untf Str"-*l

Lt,seVo,* nofahltL 6

,frfi

fq/hoC .

A A 6-,7 {= t,,, 0,:0r' 0r=t,

lartJrt

r !
A
:
(
t
'
,
t , t t , L r , , L y , ( g )0,. * : e r ,0 , r : E , ,t r = S "
t'r
('L

A
:
L
1

,'

L--t,
c(.:l

IL

Y
i("tt
do1 C 1

_ Z x4 0i'

x,,taptf,

d:l

be un'(lmr,^tU fnn "rt

Tt",)t Ca,,n

47 ' * : b"f

o.r\'a:!^

n7-

^( z)

(l)

L"(z)
1

\-1

.(t)

.c')

LL

L].

ct L -

r-l

,cD

L6

o c{})

{lt)

1j.",
Lb

L--=-.--------------l'--

-rt')
LI

r . (l l

Lt

-1lt)
LI

.(!

Lz

--)

-f(tt
L I

e(lt1
L Z

- 4. ( l ) - cc I 6L l

I l')

------..v.-.----

l:L Skl?sfte/rn'

'l tl14
"'*91N
Iplrrnr
r> SLip

potycn{ s+"{".t*
5<e- f.. - sUP- )ylie"'^,*
M it ccntgr^S::df*
(:-"'

{Iurt;',n c,txq CI ) in tfur- Sfo"ntt*i d +htLf)il'

(-3 6 nlW
r-.ux}*p&tq )
M i1 lc*1d

en#"

c,rl,.*v\^A'
/a+ ^

( had Traiu ) a^d Smo"Utir(-2" 9) ruoi

ffot J a,rr{[f tt ]
rl'z arv!/-WM b 3'"{7 [T^lbr-f'Ltv).
foof7 (erft qtu,'tt',)

16
Slip Resistance in Polycrystal
mean dissip = 3.0665 (Taylor Factor)
3.6

0.7

3.4

0.5

3.2

0.4

0.6

0.3

2.8

0.2

2.6

0.1
2.4
0
0

0.2

0.4
x

0.6

0.8

mean dissip = 3.0665 (Taylor Factor)

3.5
3
2.5

0.6

0.4

0.2
0.4

0.2
0.6
x

y
0.8

fu,1t_3W

Ple,),\vft;t
Goyyni
i
u * . - a -

Q,u*r,t PnfoJ$n tur Poi,

Elo

Xa, D nt ed 4
[a : (lo

NE {Md
-

toh-{.rw o{

loJaX

u_,,
/'\- : f

"f|'

sf"4
(h;p,rtr in MatrLL)

.{a+p) er
/

d,: |, 2, ' '- / I

bnt;lz auo A

4^:
z [( !,o,* h,",)

i: d:,;fiuo*t /* ea+h I
t^,9,;c.t.
CrrX-uhW4
Th.,a_

t7

Ca,i

,'a fhr ff*',dw'A W

ca'nbe vB^"boeJAA awtws

[_ rv1v,]O,rt

iY\fhL Slamo(aw'LfrYftn
L^ |+nu;m. gru^,": ro1'attesoth"t
ilvnfahrry b,nrnL ilrY"t''l

tAu)

lUt|r 6y ltOol

uttk lls l+rru;G a^tA,

Thw> f tt ol mvttodinn4 m,o'rq"t -^J f^.l" +lrt {vry51Qqx A ,

-fhr-

oyyo5tfr.r-( tru4- ln

[l"n*

oarr'fruttikn

wi+l,ui,\i6aI
r,oUJc,,Ws+"I

rau'dor* ym)n

olrvwlv**i64s

u il /ur<1."'y prcfwzt J f t'dv'n Uue,-"'(v'Jtn ( t e' kxluw

,t*

{ettus

Kp's.kt:
L

pla*hc ,/'zfrton^4rn ( ex*w""

{-'n"*t^l

0v roLLtnfi
)

-t-o^lLor ,1
{',3ou,,+xd het.u uvr\:-crrvvt3k/",r+wtt"k
, U T,
tY)
62. 3u7 ' t73e "florvt;<>":;fra,';n-t
me,+t{.a//lnsr A4ryl-k

18
Crystal Rotation
rotation of tensile axis
0.7

111

0.6
0.5

0.4
0.3
0.2
0.1
0

001

101

0.2

0.4
x

0.6

0.8

6/4/13 12:38 PM

D:\My Documents\Courses\ME342-Inelasti...\fcc_slip_system_map.m

1 of 3

% FCC slip systems for single crystal under uniaxial tension

%
plane
slip vector
slip_system = [ 1 1 1
1 -1 0
1 1 1
1 0 -1
1 1 1
0 1 -1
1 1 -1
1 -1 0
1 1 -1
1 0 1
1 1 -1
0 1 1
1 -1 1
1 1 0
1 -1 1
1 0 -1
1 -1 1
0 -1 -1
-1 1 1
-1 -1 0
-1 1 1
-1 0 -1
-1 1 1
0 1 -1 ];
N_slip_system = length(slip_system(:,1));
%domain_type = 1; % a simple parameterization covering several triangles
domain_type = 2; % parameterization of the standard triangle
% G. Y. Chin and W. L. Mammel,
%
Computer Solutions of the Taylor Analysis for Axisymmetric Flow
%
Trans. Metall. Soc. AIME 239, 1400-1405 (1967).

if domain_type == 1
theta = linspace(1e-4,1-1e-4,50)*acos(1/sqrt(3));
phi
= linspace(-0.2+1e-4,1-1e-4,50)*(pi/2);
elseif domain_type == 2
theta = linspace(1e-4,1-1e-4,50)*(pi/4);
phi
= linspace(1e-4,1-1e-4,50)*(pi/4);
else
disp('unknown domain_type');
end
X = zeros(length(theta), length(phi));
Y = X; x = X; y = X; z = X;
maxS = X; maxSind = X;
secondS = X; secondSind = X;
T_select = [ 0
1
1
0
1
2

0
0
1
1
1
1

1
1
1
1
2
3 ];

for ai = 1:length(theta),
for bi = 1:length(phi),
if domain_type == 1
x(ai,bi) = sin(theta(ai))*cos(phi(bi));
y(ai,bi) = sin(theta(ai))*sin(phi(bi));
z(ai,bi) = cos(theta(ai));

6/4/13 12:38 PM

D:\My Documents\Courses\ME342-Inelasti...\fcc_slip_system_map.m

elseif domain_type == 2
factor = atan(sin(theta(ai)))/(pi/4);
x(ai,bi) = sin(theta(ai))*cos(phi(bi)*factor);
y(ai,bi) =
sin(phi(bi)*factor);
z(ai,bi) = cos(theta(ai))*cos(phi(bi)*factor);
end
X(ai,bi) = 2*x(ai,bi)/(1+z(ai,bi));
Y(ai,bi) = 2*y(ai,bi)/(1+z(ai,bi));
T = [x(ai,bi) y(ai,bi) z(ai,bi)];
S = zeros(N_slip_system,1);
for i = 1:N_slip_system
n = slip_system(i,1:3);
b = slip_system(i,4:6);
S(i) = abs(dot(T,n)/norm(T)/norm(n) * dot(T,b)/norm(T)/norm(b));
end
[sortedS ind] = sort(S,1,'descend');
maxS(ai,bi) = sortedS(1);
maxSind(ai,bi) = ind(1);
secondS(ai,bi) = sortedS(2); secondSind(ai,bi) = ind(2);
end
end
X_select = zeros(length(T_select(:,1)),1);
Y_select = X_select;
maxS_select = X_select;
maxSind_select = X_select;
secondS_select = X_select;
secondSind_select = X_select;
for si=1:length(X_select),
T = T_select(si,:); T = T / norm(T);
X_select(si) = 2*T(1)/(1+T(3));
Y_select(si) = 2*T(2)/(1+T(3));
S = zeros(N_slip_system,1);
for i = 1:N_slip_system,
n = slip_system(i,1:3);
b = slip_system(i,4:6);
S(i) = abs(dot(T,n)/norm(T)/norm(n) * dot(T,b)/norm(T)/norm(b));
end
[sortedS ind] = sort(S,1,'descend');
maxS_select(si) = sortedS(1);
maxSind_select(si) = ind(1);
secondS_select(si) = sortedS(2); secondSind_select(si) = ind(2);
end
%plot results
figure(1);
contourf(X,Y,maxS,20);
hold on
plot3(X_select,Y_select,ones(size(X_select)),'ko','LineWidth',2);
hold off
xlabel('x');

2 of 3

6/4/13 12:38 PM

D:\My Documents\Courses\ME342-Inelasti...\fcc_slip_system_map.m

ylabel('y');
axis equal
title('highest Schmid factor');
figure(2);
mesh(X,Y,maxSind);
hold on
plot3(X_select,Y_select,13*ones(size(X_select)),'ko','LineWidth',2);
hold off
view([0 90]);
xlabel('x');
ylabel('y');
axis equal
title('slip system ID with highest S');
figure(3);
mesh(X,Y,secondSind);
hold on
plot3(X_select,Y_select,13*ones(size(X_select)),'ko','LineWidth',2);
hold off
view([0 90]);
xlabel('x');
ylabel('y');
axis equal
title('slip system ID with second highest S');
figure(4);
contourf(X,Y,maxS./secondS-1,20);
colorbar
hold on
plot3(X_select,Y_select,ones(size(X_select)),'ko','LineWidth',2);
hold off
xlabel('x');
ylabel('y');
axis equal
title('S_2/S_1 - 1');

3 of 3

6/4/13 1:55 PM

D:\My Documents\Courses\ME342-I...\fcc_slip_system_polycrystal.m

% FCC slip systems for poly-crystal under uniaxial tension

% Find the 5 slip systems activated for each tensile axis
% following Taylor 1938
% G. Y. Chin and W. L. Mammel,
%
Computer Solutions of the Taylor Analysis for Axisymmetric Flow
%
Trans. Metall. Soc. AIME 239, 1400-1405 (1967).
%
slip_system = [

plane
slip vector
1 1 1
1 -1 0
1 1 1
1 0 -1
1 1 1
0 1 -1
1 1 -1
1 -1 0
1 1 -1
1 0 1
1 1 -1
0 1 1
1 -1 1
1 1 0
1 -1 1
1 0 -1
1 -1 1
0 -1 -1
-1 1 1
-1 -1 0
-1 1 1
-1 0 -1
-1 1 1
0 1 -1 ];
N_slip_system = length(slip_system(:,1));
theta = [0:0.01:1]*(pi/4);
phi
= [0:0.01:1]*(pi/4);
dtheta = theta(2)-theta(1);
dphi
= phi(2)-phi(1);
Ntrial = 3;
dTfactor = 1e-4;
X = zeros(length(theta), length(phi));
Y = X; x = X; y = X; z = X;
activeS = zeros(length(theta), length(phi), N_slip_system);
dissip
= zeros(length(theta), length(phi)); weight = dissip;
Taxis
= zeros(length(theta), length(phi), 3); netrot = Taxis; dTaxis = Taxis;
Taxisnew = Taxis; Xnew = X; Ynew = Y; dX = X; dY = Y;
strain = zeros(N_slip_system,6);
rotation = zeros(N_slip_system,3);
for i = 1:length(slip_system(:,1)),
n = slip_system(i,1:3); n = n/norm(n);
b = slip_system(i,4:6); b = b/norm(b);
tmp = (n'*b + b'*n)/2;
strain(i,:) = [tmp(1,1) tmp(2,2) tmp(3,3) tmp(2,3) tmp(3,1) tmp(1,2)];
rotation(i,:) = cross(n,b)/2;
end
% compute participation ratio of each slip system
for ai = 1:length(theta),
if mod(ai,10)==0
disp(sprintf('ai = %d / %d ', ai, length(theta)));
end
for bi = 1:length(phi),

1 of 3

6/4/13 1:55 PM

D:\My Documents\Courses\ME342-I...\fcc_slip_system_polycrystal.m

factor = atan(sin(theta(ai)))/(pi/4);
x(ai,bi) = sin(theta(ai))*cos(phi(bi)*factor);
y(ai,bi) =
sin(phi(bi)*factor);
z(ai,bi) = cos(theta(ai))*cos(phi(bi)*factor);
X(ai,bi) = 2*x(ai,bi)/(1+z(ai,bi));
Y(ai,bi) = 2*y(ai,bi)/(1+z(ai,bi));
T = [x(ai,bi) y(ai,bi) z(ai,bi)];
Taxis(ai,bi,:) = T;
tmp = (T'*T - (T*T'/3)*eye(3))*(3/2);
tensile_strain = [tmp(1,1) tmp(2,2) tmp(3,3) tmp(2,3) tmp(3,1) tmp(1,2)];
%options = optimset('Display','iter','TolFun',1e-8);
options = optimset('Display','off','TolFun',1e-8);
% double # of slip systems so that coeff is non-negative
Aeq = [strain', -strain'];
beq = [tensile_strain'];
lb = zeros(length(slip_system(:,1))*2,1);
f= ones(length(slip_system(:,1))*2,1);
coeff = linprog(f, [], [], Aeq, beq, lb, [], [], options);
% half # of coeff to store as positive and negative values
compact_coeff = coeff(1:end/2)-coeff(end/2+1:end);
activeS(ai,bi,:) = compact_coeff;
dissip(ai,bi) = sum(abs(compact_coeff));
netrot(ai,bi,:) = rotation'*compact_coeff;
end
end
% compute change of Taxis
for ai = 1:length(theta),
for bi = 1:length(phi),
dTaxis(ai,bi,:) = cross(netrot(ai,bi,:),Taxis(ai,bi,:));
Taxisnew(ai,bi,:) = Taxis(ai,bi,:) + dTaxis(ai,bi,:)*dTfactor;
Xnew(ai,bi) = 2*Taxisnew(ai,bi,1)/(1+Taxisnew(ai,bi,3));
Ynew(ai,bi) = 2*Taxisnew(ai,bi,2)/(1+Taxisnew(ai,bi,3));
end
end
dX = Xnew - X; dY = Ynew - Y;
% compute area of mesh
area = zeros(size(weight));
for ai = 1:length(theta)-1,
for bi = 1:length(phi)-1,
dr1 = [x(ai+1,bi) y(ai+1,bi) z(ai+1,bi)] - [x(ai,bi) y(ai,bi) z(ai,bi)];
dr2 = [x(ai,bi+1) y(ai,bi+1) z(ai,bi+1)] - [x(ai,bi) y(ai,bi) z(ai,bi)];
dr3 = [x(ai+1,bi) y(ai+1,bi) z(ai+1,bi)] ...
- [x(ai+1,bi+1) y(ai+1,bi+1) z(ai+1,bi+1)];
dr4 = [x(ai,bi+1) y(ai,bi+1) z(ai,bi+1)] ...
- [x(ai+1,bi+1) y(ai+1,bi+1) z(ai+1,bi+1)];
area(ai,bi) = (norm(cross(dr1,dr2))+norm(cross(dr3,dr4)))/2;
end
end

2 of 3

6/4/13 1:55 PM

D:\My Documents\Courses\ME342-I...\fcc_slip_system_polycrystal.m

% compute integration (quadrature) weight

for ai = 1:length(theta),
for bi = 1:length(phi),
weight(ai,bi) = area(ai,bi)/4;
if ai > 1
weight(ai,bi) = weight(ai,bi) + area(ai-1,bi)/4;
end
if bi > 1
weight(ai,bi) = weight(ai,bi) + area(ai,bi-1)/4;
end
if ai > 1 && bi > 1
weight(ai,bi) = weight(ai,bi) + area(ai-1,bi-1)/4;
end
end
end
mean_dissip = sum(sum(dissip.*weight))/(4*pi/48);
% plot results
fs = 17;
figure(1);
contourf(X,Y,dissip,50);
colorbar
set(gca,'FontSize',fs);
xlabel('x'); ylabel('y');
axis equal
title(sprintf('mean dissip = %.4f

(Taylor Factor)',mean_dissip));

figure(2);
mesh(X,Y,dissip);
set(gca,'FontSize',fs);
xlabel('x'); ylabel('y');
axis equal
view([30 80]);
title(sprintf('mean dissip = %.4f

(Taylor Factor)',mean_dissip));

figure(3);
skip = 10;
contour(X,Y,dissip,50); hold on
quiver(X(1:skip:end,1:skip:end), Y(1:skip:end,1:skip:end), ...
dX(1:skip:end,1:skip:end),dY(1:skip:end,1:skip:end));
hold off
set(gca,'FontSize',fs);
xlabel('x'); ylabel('y');
axis equal
title('rotation of tensile axis');
t001 = text(-0.01,0.02,'001');
t101 = text( 0.80,0.02,'101');
t111 = text( 0.70,0.70,'111');

3 of 3