UNIVERSITY OF CALIFORNIA BERKELEY Structural Engineering,Department of Civil Engineering Mechanics and MaterialsFall 2010 Professor: S. GovindjeeThe Symmetric Identity: Isym→12(δikδjl+ δilδjk)The derivative of a second order tensor with itself is the fourth order iden-tity tensor I. In component form this reads, for a tensor A with componentsAij, asIijkl=∂Aij∂Akl= δikδjl; (1)i.e. if i = k and j = l then the result is unity and if one of the equalitiesis violated, then the result is zero. So, for example, ∂A12/∂A12= 1 and∂A12/∂A21= 0 etc.When the tensor A ∈ S (the set of symmetric tensors), then we cometo the result that the identity is no longer given by (1). Rather it is givencomponent-wise by the expression in the header. This then returns the some-what non-intuitive result that ∂A12/∂A12= Isym1212= 1/2 (not unity!) eventhough ∂A11/∂A11= Isym1111= 1. How can this be?Consider a function f : S → R. Its derivative1is a linear mapping Df :S → R. As such we can represent it by a tensor:B : C = Df[C] . (2)The tensor B characterizing the derivative of f(·) in the (symmetric) direc-tion C is non-unique since one can add an arbitrary skew-symmetric tensorto B without changing the result. The symmetric part of B, viz.12(B +BT),is however unique. This is the tensor which we choose to use to representthe derivative Df. We do this for two reasons: (1) it is unique, thus there isno ambiguity, and (2) it will always produce the correct result for the rate ofchange of f(·) in a symmetric direction. Note it only makes sense to discussrates of change of f(·) in symmetric directions, since f(·) is only defined overthe space of symmetric tensors.Let us now apply this result to the component function f(A) = (ˆei⊗ˆej) : A = Aij. To compute the derivative of f(·) we apply the directional1Assuming it exists1derivative formuladdαα=0(ˆei⊗ˆej) : (A + αH) = Df[H] , (3)where H ∈ S. Taking the derivative with respect to α, setting α to zero,leads to:(ˆei⊗ˆej) : H = Df[H] , (4)which allows us to identify an expression for B asˆei⊗ˆej. The unique partof B which provides the derivative information is the symmetric part; viz.12(ˆei⊗ˆej+ˆej⊗ˆei). This expression gives us the derivative ∂Aij/∂A. If wenow compute its components we get∂Aij∂Akl=∂Aij∂A: (ˆek⊗ˆel) =12(δikδjl+ δilδjk) . (5)Thus we arrive at the desired expression – an expression for the fourth orderidentity tensor over the space of symmetric tensors. Observe that this ex-pression yields the results ∂A11/∂A11= Isym1111= 1, ∂A12/∂A12= Isym1212= 1/2,as well as ∂A12/∂A21= Isym1221=