23.03.2013 Views

CONTINUUM MECHANICS for ENGINEERS

CONTINUUM MECHANICS for ENGINEERS

CONTINUUM MECHANICS for ENGINEERS

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Solving Eq 4.10-9 <strong>for</strong> and using Eq 4.10-11, we may write<br />

˙n i<br />

( ) −<br />

˙ n= v n− nΛ˙ / Λ = D + W n D nnn<br />

i i,j j i ij ij j qk q k i<br />

(p) If now ni is chosen along a principal direction of D so that Dn = ij j D(p) ni<br />

(p = 1, 2, 3) where D (p) represents a principal value of D, then<br />

˙ ( p)<br />

( p)<br />

i ( p) i ij j<br />

n = D n + W n − D n n n = W n<br />

( p)<br />

( p)<br />

( p)<br />

( p)<br />

( p) q q i ij j<br />

(4.10-20)<br />

(p) (p)<br />

since n = 1. Because a unit vector can change only in direction, Eq 4.10-<br />

q nq<br />

20 indicates that Wij gives the rate of change in direction of the principal axes<br />

of D. Hence the names vorticity or spin given to W. Additionally, we associate<br />

with W the vector<br />

1<br />

1<br />

wi = εijkvk,j or w = curl v (4.10-21)<br />

2<br />

2<br />

called the vorticity vector, by the following calculation,<br />

( )<br />

1<br />

1<br />

pqi i 2 pqi ijk k,j 2 pj qk pk qj k,j<br />

ε w = ε ε v = δ δ −δ<br />

δ v<br />

1<br />

= v v W<br />

(4.10-22)<br />

2 ( q,p − p,q )= qp<br />

Thus if D ≡ 0 so that L ij = W ij, it follows that dv i = L ijdx j = W ijdx j = ε jikw kdx j<br />

and since ε jik = –ε ijk = ε ikj,<br />

dv i = ε ijkw j dx k or dv = w × dx (4.10-23)<br />

according to which the relative velocity in the vicinity of p corresponds to a<br />

rigid body rotation about an axis through p. The vector w indicates the<br />

angular velocity, the direction, and the sense of this rotation.<br />

To summarize the physical interpretation of the velocity gradient L, we<br />

note that it effects a separation of the local instantaneous motion into two parts:<br />

1. The so-called logarithmic rates of stretching, D (p), (p = 1,2,3), that is, the<br />

eigenvalues of D along the mutually orthogonal principal axes of D, and<br />

Λ˙ / Λ<br />

d lnΛ<br />

dt<br />

= ( ) ( p) ( p) ( p)<br />

( p)<br />

= = =<br />

Dnn nD n D<br />

ij i j i ( p)<br />

i<br />

2. A rigid body rotation of the principal axes of D with angular velocity w.<br />

( p)<br />

(p)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!